This Article Abstract Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Google Scholar Articles by Busch, J. W. Articles by Schoen, D. J. PubMed PubMed Citation Articles by Busch, J. W. Articles by Schoen, D. J.

Molecular Characterization of Lal2, an SRK-Like Gene Linked to the S-Locus in the Wild Mustard Leavenworthia alabamica

Department of Biology, McGill University, Montreal, Quebec H3A 1B1, Canada

¹ Corresponding author: Department of Biology, McGill University, Montreal, Quebec H3A 1B1, Canada.
E-mail: jeremiah.busch{at}mcgill.ca

	ABSTRACT

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS LITERATURE CITED

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

Single-locus sporophytic SI has been studied in a number of genera in the mustard family (Brassicaceae) (BATEMAN 1954, 1955; THOMPSON 1957; LLOYD 1967; SAMPSON 1967; SCHIERUP et al. 2001). Molecular genetic characterization of SI in Brassica and Arabidopsis species has shown that the S-locus consists of a tightly linked cluster of genes encoding proteins that function together as a receptor-ligand system (SCHOPFER et al. 1999; TAKASAKI et al. 2000; for a review see FOBIS-LOISY et al. 2004). One of these genes, the S-locus receptor kinase (SRK), codes for a membrane-bound protein expressed in the stigma that binds specifically to the S-locus cysteine-rich ligand (SCR/SP11), which is expressed in the tapetum of anthers (KACHROO et al. 2001; TAKAYAMA et al. 2001). The successful binding of SCR by SRK initiates a signaling cascade that prevents pollen tubes from penetrating the stigmatic surface, although the roles of other genes in this pathway are still under investigation (MURASE et al. 2004; LIU et al. 2007). It is thought that the linked system of co-evolved SRK and SCR genes has been maintained by natural selection for at least 25–40 million years (UYENOYAMA 1995). As would be expected in a system in which selection is strongly negatively frequency dependent, nucleotide polymorphism is pronounced among S-alleles in Brassica campestris and B. oleracea (HINATA et al. 1995), Raphanus sativus (OKAMOTO et al. 2004), Arabidopsis lyrata and A. halleri (SCHIERUP et al. 2001; CASTRIC and VEKEMANS 2004; PRIGODA et al. 2005), and Capsella grandiflora (PAETSCH et al. 2006). Many of these investigations have also revealed dominance interactions among S-alleles in the pollen and/or the stigma, a phenomenon consistent with a sporophytic mechanism of SI.

Sequence-level variation among S-alleles sampled within and between species provides insight into the historical evolution of a locus under strong balancing selection (RICHMAN et al. 1996; IGIC and KOHN 2001). Since S-alleles are expected to be favored when rare by negative frequency dependence, and thus be maintained in populations for long periods, coalescence times at these loci might often predate speciation events and could provide insight into the long-term effective sizes of species (RICHMAN and KOHN 2000). In fact, alleles of SRK sampled from A. lyrata and C. grandiflora show a pronounced pattern of trans-specific and transgeneric polymorphism and are exceedingly diverse at synonymous and nonsynonymous sites (SCHIERUP et al. 2001; PAETSCH et al. 2006). SRK alleles sampled from species of Brassica and Raphanus also show a pattern of transgeneric polymorphism, although the amount of sequence variation found within each species is much more modest (HINATA et al. 1995; OKAMOTO et al. 2004). A historical bottleneck in the lineage leading to Brassica and Raphanus may have caused the loss of sequence-level variation within species. There is also a great deal of divergence between SRK sequences sampled from species of Arabidopsis and Capsella compared to those sampled from species of Brassica and Raphanus (FOBIS-LOISY et al. 2004). This marked divergence between lineages, coupled with the fact that sequences flanking the S-locus in Arabidopsis and Brassica are derived from alternative chromosomal regions, suggests that a translocation of the entire S-locus region may have occurred (KUSABA et al. 2001).

This article describes sequence-level variation in SRK-like genes in a species of Leavenworthia, a small genus from a portion of the Brassicaceae that has served as a model for evolutionary studies on the loss of SI (ROLLINS 1963; LLOYD 1965). In terms of molecular characterization of SRK and SRK-like genes, the genus Leavenworthia is unstudied. We expect that a locus orthologous to SRK in L. alabamica should satisfy the following criteria on the basis of patterns observed in both Arabidopsis and Brassica (CHARLESWORTH et al. 2000): (1) sequences of putative S-alleles should exhibit high sequence identity to known alleles of SRK; (2) sequences should cosegregate with S-alleles as inferred from diallele crosses between siblings; (3) the putative SRK ortholog must exhibit high levels of synonymous and especially nonsynonymous diversity and a pattern of selective constraint similar to SRK; and (4) sequences of the putative SRK ortholog should be expressed in the stigmas of maturing buds (although we note that some SRK alleles in A. lyrata are also expressed in leaves as well as in buds (PRIGODA et al. 2005). We present evidence for each of these expectations and conclude that one of the studied loci (Lal2) either is orthologous to SRK or is closely linked to SRK and therefore may be a useful tool for typing diversity at the S-locus in natural populations of this species (CHARLESWORTH et al. 2006; HAGENBLAD et al. 2006).

	MATERIALS AND METHODS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS LITERATURE CITED

 
Study system and rationale:
Leavenworthia is thought to have diverged from the common ancestor of Arabidopsis 13–19 MYA and from the common ancestor of Brassica 16–21 MYA (MITCHELL-OLDS 2001; BEILSTEIN et al. 2006). There are eight extant species of Leavenworthia, and within the group there have been at least three independent derivations of self-fertilization (ROLLINS 1963; BECK et al. 2006). Leavenworthia is endemic to the southeastern United States and is found under a restricted range of ecological conditions (ROLLINS 1963). A long history of work in this group points to a role of reproductive assurance in the evolution of self-compatibility (LLOYD 1965; BUSCH 2005a) and the loss of genetic variation in self-compatible populations (LIU et al. 1998, 1999). Leavenworthia alabamica and L. crassa possess considerable intraspecific variation in the presence or absence of SI and are therefore interesting models for the study of S-locus variation.

Plant growth, diallele crosses, and pollen tube visualization:
Seeds were collected from two natural SI populations of L. alabamica during the spring of 2005; these are referred to as the "Isbell" and "Waco" populations (see BUSCH 2005b for geographic locations). Plants were raised from seed originally collected from 20 distinct plants in these populations. Seeds were germinated according to established protocols (BUSCH 2005b), and the plants were grown to flowering in the McGill University Phytotron greenhouse under a combination of natural and supplemental light. Plants collected from the same population were paired haphazardly to generate 10 full-sib families; the parents of families 1–5 were of "Waco" origin, whereas the parents of families 6–10 were collected from "Isbell."

Sixteen or more offspring from each of the 10 full-sib families were germinated in March of 2006 and raised to flowering for the diallele crosses. Within each diallel, replicates of three or more reciprocal crosses were conducted for each parental combination. To carry out a cross, flowers were emasculated, anthers were removed with forceps, and pollen was placed directly on the stigmas of newly opened flowers. Forceps were cleaned with a 70% ethanol solution between individual crosses. Seeds were collected 1 month later as they ripened. Individual plants were considered compatible if the majority of crosses between them yielded large fruits (with four or more seeds). Single-seeded fruits were considered to reflect an incompatible cross with some degree of leakiness (SCHIERUP et al. 2001). Two- and three-seeded fruits were also sometimes produced; these were enumerated separately. Individuals within a diallele cross were placed in the same phenotypic group when they exhibited identical patterns of compatibility.

For visualization of the SI reaction at the level of pollen tube growth, compatible or incompatible pollen was applied to stigmas. Stigmas were harvested after 24 hr and fixed in acetic acid–ethanol for 1 week, then rinsed, cleared in NaOH, and stained using a 0.1% aniline blue solution (KHO and BAER 1968). Fluorescence microscopy with a DAPI filter and a digital camera was used to document the SI reaction.

Cloning and sequencing of SRK gene family members:
Leaves were collected from parents and offspring and flash frozen in liquid nitrogen. DNA was extracted from leaves using DNeasy plant mini kits (QIAGEN, Mississauga, ON, Canada). Degenerate primers that amplify SRK in Arabidopsis and Brassica were used to amplify the S-domain of SRK-like sequences (CHARLESWORTH et al. 2000; SCHIERUP et al. 2001; MABLE et al. 2003; supplemental Table 1). Primers were also designed to amplify conserved regions of the kinase domain, which are less variable than the extracellular S-domain (CHARLESWORTH et al. 2003a; kinase 1F: RCTTCARCAYATHAAYCTTG; kinase 1R: TCTTGGTGAAGRTAYARAAG; kinase 2F: ACGGGTGTGTGTATTTGGACTGGA; kinase 2R: AGAATATY CCRTCCATYGC). Kinase 1F and kinase 1R are found at the 3'-end of exon 4 and the middle of exon 5, respectively; kinase 2F and kinase 2R were also designed from the middle of exon 5 and the 5'-end of exon 6. PCRs were conducted at low stringency under the following conditions: 94° for 5 min with 32 cycles of a 30-sec denaturation step at 94°, annealing at 40° for 1 min, and extension at 72° for 1 min, followed by a final 10-min extension step at 72°. Amplicons were cloned using TOPO TA cloning kits (Invitrogen, San Diego) and grown overnight on LB–agarose plates. At least eight clones were randomly selected from each transformed sample. Plasmid DNA was isolated with standard mini-prep procedures. Insert sizes of clones were quantified by digestion with EcoRI and gel electrophoresis. All positive clones near the expected size were sequenced in both the forward and reverse direction using universal M13 primers. DNA templates were purified and cycle sequenced using the v. 3.1 BigDye terminator ready reaction mix and analyzed using an ABI 3730 automated sequencer (Applied Biosystems, Foster City, CA). The GenBank accession numbers for sequence data are EU394446–EU394520. Sequences were aligned using Codon Code Aligner (Codoncode). Nucleotide sequences and their translated protein sequences were screened against known SRK sequences found in closely related species by using the BLAST algorithm (ALTSCHUL et al. 1990).

Phylogenetic analyses, sequence diversity estimates, and detecting selection:
Within each family of DNA sequences, alignments were manually edited on the basis of published amino acid alignments of S-domain genes. Phylogenetic relationships among sequences were studied using the neighbor-joining algorithm with Jukes–Cantor correction across all sites. Synonymous and nonsynonymous nucleotide diversities were estimated using DNAsp v. 4.0 (ROZAS et al. 2003). Nucleotide diversity among Lal1 sequences was estimated ignoring pseudogenes. The seven longest Lal2 DNA fragments were translated into the predicted amino acid sequences and aligned, and a sliding window of nucleotide variability was calculated using DNAsp. This sliding window provided an aid for nucleotide sequence alignments, with a window size of 50 nucleotides and a step size of one nucleotide. The nucleotide diversities of SRK alleles sampled within A. lyrata, C. grandiflora, R. sativus, and Brassica species were also calculated for purposes of comparison.

To analyze patterns of selection at the codon level in the Lal2 sequences, we used the phylogenetically based, codon-level analysis of nonsynonymous and synonymous substitution (NIELSEN and YANG 1998; YANG et al. 2005). This approach calculates posterior probabilities of positive selection (as determined from the ratio of nonsynonymous substitution per nonsynonymous site to synonymous substitution per synonymous site) for individual codons, using an empirical Bayesian approach. The maximum-likelihood tree topology for the Lal2 sequences was first obtained using the branch-and-bound algorithm as implemented in PAUP for the Kimura two-parameter model (KIMURA 1980; SWOFFORD 2002), and the selection analysis was conducted using the CODEML program of YANG (1997), as decribed by SAINUDIIN et al. (2005). The number of Lal2 sequences available to us for analysis (10 sequences) was much smaller than those analyzed by SAINUDIIN et al. (2005), who used the same method to detect specific classes of positive selection (e.g., polar- or volume-changing amino acid substitutions) in SRK sequences of Brassica and Arabidopsis. Hence, our focus was on evidence of any type of positively selected amino acid substitution, regardless of the particular type of property-changing characteristic. The Lal2 codons available for analysis by this method correspond to amino acid residues 154–410 of B. oleracea.

Cosegregation of SRK-like sequences with the SI reaction:
SRK-like sequences were tested for cosegregation with incompatibility alleles determined by diallele crosses; these were identified by analyzing seed set resulting from crosses between full siblings. All incompatibility groups were classified independently and without prior knowledge of sequence data. Sequence-specific primers were designed and used to amplify SRK-like DNA sequences in the offspring of a full-sib family (supplemental Table 2). SRK-like sequences were inferred to have cosegregated with S-alleles if the sequence was found in individuals from the dialleles that were incapable of mating with one another. PCR products were visualized on 1x TBE gels containing 1% agarose and verified by DNA sequencing.

Gene expression of the Lal2-3 allele:
L. alabamica individuals of known genotype were selected for gene expression analysis of the Lal2-3 allele. Tissue was harvested from the leaves, styles, and stigmas. Style tissue was isolated from stigmas and was collected from open flowers. Stigma tissue was harvested from buds 3 days (day –3) and 1 day prior (day –1) to flower opening as well as from open flowers (day 0). Tissue samples from 70 flowers harvested from the same individual were pooled and flash frozen in liquid nitrogen. Frozen tissue was ground to a fine powder in liquid nitrogen using a Kontes pestle (VWR, Mississauga, ON, Canada). Total RNA was extracted using the RNeasy plant mini kit (QIAGEN). Extracted RNA was treated with RNAse-free DNAse (Promega, Madison, WI) to eliminate all possible DNA contamination prior to the reverse transcriptase (RT)–PCR) experiments and the enzyme was removed using the cleanup protocol specified for the QIAGEN RNeasy kit. For each sample, 1 µg of RNA was reverse transcribed using oligo(dT) primer (Invitrogen Canada, Burlington, ON, Canada) and Superscript II RT (Invitrogen).

Primers for the Lal2-3 allele (W14F: TCAGCTTGAATACCGTACTGACTTA; W14R: GTCCCTGCTTGTGTTGATGGA) and an actin control (actinF: TATGCACTTCCACATGCTAT; actinR: CTTTGCGATCCACATCTGCTG) were designed. Both of the actin primers were intron spanning to ensure that only cDNAs would be amplified by PCR. PCR products were visualized on a 1% agarose gel by ethidium bromide staining and product identity was also confirmed by sequencing (Genome Quebec, Montreal). No reverse transcriptase controls were used for each sample to confirm the absence of DNA contamination.

	RESULTS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS LITERATURE CITED

 
Cross-incompatibility and interactions among inferred S-alleles:
Tables 1–5 summarize the results of crosses conducted between full siblings in each of the different diallele crosses. In families 4, 5, and 10, four phenotypic groups were identified on the basis of the patterns of cross-incompatibility. Tables 1 and 2 show that family 4 and family 5 exhibited identical patterns of cross-incompatibility and therefore have identical types of interactions among S-alleles. In particular, in both families, phenotypic group II is incompatible with group I when acting as a maternal parent, while the reciprocal cross successfully produces fruit. This supports the existence of alleles (S_i in family 4 and S_m in family 5) that act codominantly in the pistil yet are recessive to some alleles in pollen. Moreover, matings between groups I and III were compatible in both families, suggesting that there are two alleles (S_k in family 4 and S_o in family 5) that are recessive in both the stigma and pollen. Table 3 shows that, in family 10, most alleles were codominant with one another, but that groups I and III were reciprocally compatible. This pattern suggests the existence of an allele (S_s) with recessive gene action in stigmas and pollen. The overall level of incompatibility between plants sharing S-alleles was lower in family 10, indicating some degree of leakiness in the self-incompatibility mechanism.

View larger version (50K): [in this window] [in a new window] [Download PPT slide]   FIGURE 1.— Pollen tube germination and growth after compatible and incompatible crosses. A compatible cross (a) with pollen tube growth into the style and an incompatible cross (b) with pollen tube growth inhibited.

View larger version (17K): [in this window] [in a new window] [Download PPT slide]   FIGURE 2.— Unrooted neighbor-joining tree with Jukes–Cantor correction of SRK-like sequences amplified in L. alabamica and several relatives. The tree was generated using a pairwise deletion algorithm to handle missing sites. Lal2 and SRK sequences marked with asterisks were shown in this study to be linked to the SI reaction in L. alabamica. Dots denote nodes with 95% bootstrap support. The scale denotes the number of substitutions.

Cosegregation between SRK-like sequences and inferred S-alleles:
Sequence-specific PCR allowed us to amplify six SRK-like sequences in the L. alabamica plants used for diallele analysis of SI (LalSRK, Lal2-1, Lal2-3, Lal2-4, Lal2-5, Lal2-6). Sequences could not be amplified from every plant used in the dialleles, likely because of lack of homology with the PCR primers employed (based mostly on Brassica and Arabidopsis sequence information) and the possibility that Lal2 may not be found in all Leavenworthia S-haplotypes (see below). These six sequences cosegregate with S-alleles inferred on the basis of compatibility patterns in all 113 cases where data were available from crossing studies (Table 6). All Lal2 sequences tested for cosegregation with inferred S-locus alleles were found to segregate in offspring according to Mendelian expectations (² = 4.75, d.f. = 8; P > 0.10; Table 6), as would be expected if these alleles were inherited at a single locus. Given our sample size of n = 97 for tests of cosegregation of Lal2 with inferred S-alleles, we can state with 95% confidence that Lal2 lies at most 3.1 cM from the S-locus in L. alabamica (STEVENS 1942).

We amplified a Lal2 sequence (Lal2-3) together with LalSRK from the same individual parent used to generate family 3 (Table 6). Twenty-four offspring of family 3 were genotyped for both of these sequences. Thirteen offspring possessed LalSRK while the other 11 plants possessed Lal2-3; this pattern is consistent with the 1:1 inheritance of these sequences, as would be expected if they were alleles segregating at a single locus or were alleles segregating in repulsion phase at two tightly linked loci.

In addition to these Lal2 sequences for which there is direct evidence of cosegregation with inferred S-locus alleles, there are three other Lal2 sequences for which crossing data were insufficient to test for cosegregation with inferred SI alleles (Lal2-7, Lal2-8, and Lal2-9), as well as two sequences (Lal2-2a and Lal2-2b) that are 99% identical at the nucleotide level. The 1-bp difference between Lal2-2a and Lal2-2b may be the result of sequencing error, as multiple clones were not sequenced. Lal2-2a and Lal2-2b were amplified from family 1 and family 3, in which cosegregation of incompatibility patterns with genomic sequences cannot be investigated because of recessive gene action. Additional evidence, apart from their nearly complete sequence identity, suggests that Lal2-2a and Lal2-2b either are copies of the same S-allele or are linked in coupling phase to the same S-allele. Lal2-2a and Lal2-2b were found to be reciprocally cross-incompatible and had to be paired with other plants to generate families 1 and 3. Sequence-specific PCR shows that Lal2-2a and Lal2-2b were segregating in both families 1 and 3. Lal2-2a was found in 4 of 16 offspring in family 1 and Lal2-2b in 10 of 16 offspring in family 3. Therefore, both of the originally cross-incompatible parents used to generate families 1 and 3 were heterozygous for these Lal2 alleles. Lal2-2 was the only Lal2 sequence amplified from multiple parental individuals, a finding that is consistent with the hypothesis that recessive S-alleles should equilibrate at higher frequencies in populations (SCHIERUP et al. 1997; UYENOYAMA 2000).

Tests of cosegregation with inferred S-locus alleles were also conducted for sequences of Lal3 in family 1. A sequence-specific PCR conducted with the offspring of family 1 successfully amplified the Lal3 variant in all offspring, an inheritance pattern that indicates that Lal3 does not cosegregate with S-locus alleles. Tests for cosegregation with SI were not conducted for Lal1 variants as this locus had relatively low total nucleotide diversity and is therefore not a viable candidate for a S-linked locus, as has been observed for Lal2.

Patterns of nucleotide variability and patterns of selection at the codon level:
Figure 3 summarizes diversities among alleles at synonymous and nonsynonymous sites of Lal2 and of the other members of this gene family. The Lal2 diversity levels are significantly higher than those observed at all other gene family members for similar lengths of sequence. The average synonymous diversity of Lal2 in two populations of L. alabamica (_s = 0.216) was at least 10 times higher than the average diversity observed at loci in populations of outcrossing plants (GLÉMIN et al. 2006). More modest levels of synonymous diversity were observed at Lal1 (_s = 0.069) and Lal3 (_s = 0.076), although the levels of variability were also high compared to what is expected for loci not experiencing balancing selection. Perhaps most striking were the extremely high estimates of nonsynonymous diversity at the Lal2 locus (_a = 0.160). Lal2 sequences possess synonymous and nonsynonymous diversity levels similar overall to those observed at the SRK locus in Raphanus and Brassica, although they were lower than the levels reported in Arabidopsis and Capsella (Table 7).

View larger version (10K): [in this window] [in a new window] [Download PPT slide]   FIGURE 3.— Estimates of synonymous (solid bars) and nonsynonymous (open bars) nucleotide diversity () at five loci studied in two natural populations of L. alabamica.

View larger version (23K): [in this window] [in a new window] [Download PPT slide]   FIGURE 4.— Alignment of the predicted amino acid sequences of Lal2 alleles. Dots denote amino acid identity with the topmost reference allele. Triangles point to the 12 conserved cysteine residues typically found in S-domain loci. The underlined region of 10 amino acids indicates a strongly conserved motif characteristic of receptor-like kinases. The Lal2-1 allele does not possess four of the last five conserved cysteine residues.

View larger version (14K): [in this window] [in a new window] [Download PPT slide]   FIGURE 5.— Sliding-window analysis of nucleotide diversities in the S-domain of SRK and Lal2. The nucleotide diversities () of Arabidopsis (solid line) and Brassica (dashed line) SRK were estimated using 8 and 11 alleles, respectively. Leavenworthia Lal2 diversity (stippled line) was estimated using seven sequences. The window size is 50 nucleotides with a step size of 1 nucleotide.

Expression of Lal2 in floral and vegetative tissues:
The Lal2-3 allele is expressed at low levels in leaves, at an intermediate level in the styles of mature flowers, and at higher levels in stigmatic tissue (Figure 6). The expression of Lal2-3 at low levels in leaves does not necessarily rule out this locus as being orthologous to SRK, as work in A. lyrata has shown that a class of SRK alleles is transcribed in vegetative tissue (PRIGODA et al. 2005). Lal2 transcripts are most prevalent in stigmatic tissue immediately prior to anthesis but decline following anthesis, a pattern that is consistent with the temporal expression of SRK in Brassica (STEIN et al. 1996; KUSABA et al. 2001). The decline in Lal2 expression in open flowers is most likely caused by these products being translated into functional proteins during or immediately prior to anthesis. These results suggest that the gene is expressed predominantly in stigmatic tissue, where SRK and SLG have been shown to interact with SCR to block the growth of self-pollen grains (KACHROO et al. 2001; KUSABA et al. 2001; TAKAYAMA et al. 2001).

View larger version (61K): [in this window] [in a new window] [Download PPT slide]   FIGURE 6.— Expression of Lal2-3 in vegetative and floral tissues. Samples were harvested 1 (–1) and 3 (–3) days before flowering and on the day the flowers opened (0). PCR amplified actin and Lal2-3 bands are from cDNAs generated by RT–PCR; gDNA stands for a positive control of genomic DNA with the Lal2-3 sequence. RT–PCR controls demonstrate that samples were not contaminated with genomic DNA.

	DISCUSSION

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS LITERATURE CITED

Our finding of one allele of SRK and many sequences of Lal2 that are less similar to SRK implies either that (1) successful amplification of more than one L. alabamica SRK was not achieved because of marked divergence between SRK-specific primers (based on sequences from Arabidopsis and Brassica) and their target sites in L. alabamica or (2) Lal2 sequences are, in fact, divergent albeit functional S-alleles and are therefore orthologous to SRK. Experiments involving gain of function may help to establish whether Lal2 is necessary and sufficient to cause SI in L. alabamica (TAKASAKI et al. 2000; TAKAYAMA et al. 2001). For example, experiments whereby self-compatible plants are transformed with an intact Lal2-SCR haplotype can be used to determine whether Lal2 codes for SI in L. alabamica (NASRALLAH et al. 2002, 2004).

If Lal2 sequences are functional S-domain receptors affecting SI in L. alabamica, their divergence from SRK alleles would demand explanation. It is possible that a proto-Lal2 sequence may have been captured from another receptor kinase locus, generating divergence between the proto-Lal2 sequence and the SRK sequences in the population. If this gene conversion event caused a bottleneck in S-allele number, the diversification of Lal2 would then simultaneously generate novel S-alleles and explain the similar diversity characteristics of Lal2 and SRK if the same amino acids influence specificity. The relatively high K_a/K_s ratio among Lal2 sequences is compatible with a recent and rapid divergence among Lal2 sequences in their specificities following a bottleneck in S-allele number caused by gene conversion. A scenario such as this has been invoked to explain the ancient divergence between Brassica and Arabidopsis SRK sequences (KUSABA et al. 2001) and the evolution of a divergent clade of closely related SRK alleles in A. lyrata (e.g., SRK3, SRK6, SRK8, and SRK14) that are expressed in leaves (PRIGODA et al. 2005). This group of SRK alleles exhibits relatively shallow genetic divergence in both the S and kinase domains, consistent with a recent radiation (CHARLESWORTH et al. 2003a).

Another possibility is that Lal2 is a recently duplicated locus linked to SRK. The finding of a highly polymorphic SRK-like gene linked to SRK would suggest the existence of a locus analogous to SLG of Brassica (STEIN et al. 1991). SLG is a S-domain protein lacking a kinase domain that is tightly linked to SRK (BOYES and NASRALLAH 1993), facilitates the incompatibility reaction (SUZUKI et al. 2000; TAKASAKI et al. 2000), and is found in most S-haplotypes investigated to date in Brassica (SUZUKI et al. 2000). SLG exhibits signatures of balancing selection and patterns of constraint similar to SRK even though it is not required for SI (HINATA et al. 1995; KUSABA et al. 1997), most likely because of gene conversion between SRK and SLG of the same S-haplotype (WATANABE et al. 1994; FUJIMOTO et al. 2006). There are also other S-domain loci that play no known role in SI that are embedded in the gene-dense S-locus regions of mustard species (SUZUKI et al. 1999). In particular, there are at least three S-linked receptor kinases that are not required for the SI response in Brassica (STEIN et al. 1991; SUZUKI et al. 1997), and a receptor kinase (Aly8) with high homology to SRK is linked to the S-locus in A. lyrata (CHARLESWORTH et al. 2003b). If Lal2 is the product of recent gene duplication at the S-locus, the elevated diversity of Lal2 and its similar diversity characteristics to SRK may involve infrequent gene conversion with SRK or parallel evolutionary constraints on the S-domain portion of the molecule.

Although it is not yet known with certainty whether Lal2 sequences are orthologous to SRK, the work reported here broadens the diversity of mustard species in which the diversity of the S-locus has been studied at the sequence level. Moreover, that Lal2 is S-linked means that it may be employed to infer S-locus genotypes and to investigate selection at the S-locus in L. alabamica. There are several loci closely linked to SRK in A. lyrata that can, in theory, be used to type S-alleles within populations, given the restricted recombination in this region (CHARLESWORTH et al. 2006; KAMAU et al. 2007). However, it should also be mentioned that attempts to exploit the close linkage of Aly8 to SRK in A. lyrata to infer S-haplotypes have been complicated by the presence of a tightly linked paralog of Aly8, which causes S-haplotypes to harbor more than one Aly8-like sequence (HAGENBLAD et al. 2006). If Lal2 is to be applied to infer dynamics at the S-locus in L. alabamica, it will be important to test each Lal2 sequence for linkage with an S-haplotype inferred through controlled dialleles and to continue tests of Mendelian single-locus inheritance.

Identification of S-linked loci in SI species and their close relatives that are self-compatible should provide insights into the evolution of self-fertilization from the outcrossing condition (BECHSGAARD et al. 2006; LIU et al. 2007; SHERMAN-BROYLES et al. 2007), which is a pervasive evolutionary trend in angiosperms (DARWIN 1876; BAKER 1955; STEBBINS 1974; IGIC et al. 2006). Studies of S-linked Lal2 variation within this and other species of Leavenworthia may be especially valuable for the insights that they may provide into the population–genetic mechanisms underlying the maintenance and loss of SI in natural populations of flowering plants (BUSCH and SCHOEN 2008).

	ACKNOWLEDGEMENTS

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS LITERATURE CITED

 
We thank V. Koelling for generously providing seed collected in wild populations of L. alabamica, H. Jensen for laboratory assistance, and A-M. L'Heureux and K. Vitale for conducting crosses in the greenhouse. We acknowledge Genome Quebec for conducting sequencing reactions and D. Charlesworth, B. Mable, and J. Nasrallah for helpful discussions. Comments by S. Hoebee, E. Newbigin, and two anonymous reviewers improved the article. D.J.S. acknowledges continuing support from the Natural Sciences and Engineering Research Council of Canada (NSERC). J.S. was supported by NSERC. J.W.B. was supported with funds from McGill University's Tomlinson Fellowship program and by NSERC.

	FOOTNOTES

 
Sequence data from this article have been deposited with the EMBL/GenBank Data Libraries under accession nos. EU394446–EU394520.

	LITERATURE CITED

TOP ABSTRACT MATERIALS AND METHODS RESULTS DISCUSSION ACKNOWLEDGEMENTS LITERATURE CITED

 

ALTSCHUL, S. F., W. GISH, W. MILLER, E. W. MYERS and D. J. LIPMAN, 1990 Basic local alignment search tool. J. Mol. Biol. 215: 403–410.[CrossRef][Medline]

BAKER, H. G., 1955 Self-compatibility and establishment after long distance dispersal. Evolution 9: 347–349.[CrossRef]

BATEMAN, A. J., 1954 Self-incompatibility systems in angiosperms. II. Iberis amara. Heredity 8: 305–322.[Medline]

BATEMAN, A. J., 1955 Self-incompatibility systems in angiosperms. III. Cruciferae. Heredity 9: 53–68.

BECHSGAARD, J. S., V. CASTRIC, X. VEKEMANS, D. CHARLESWORTH and M. H. SCHIERUP, 2006 The transition to self-compatibility in Arabidopsis thaliana and evolution within S-haplotypes over 10 million years. Mol. Biol. Evol. 23: 1741–1750.[Abstract/Free Full Text]

BECK, J. B., I. A. AL-SHEHBAZ and B. A. SCHAAL, 2006 Leavenworthia (Brassicaceae) revisited: testing classic systematic and mating system hypotheses. Syst. Bot. 31: 151–159.[CrossRef]

BEILSTEIN, M. A., I. A. AL-SHEHBAZ and E. A. KELLOGG, 2006 Brassicaceae phylogeny and trichome evolution. Am. J. Bot. 93: 607–619.[Abstract/Free Full Text]

BOYES, D. C., and J. B. NASRALLAH, 1993 Physical linkage of the SLG and SRK genes at the self-incompatibility locus of Brassica oleracea. Mol. Gen. Genet. 236: 369–373.[CrossRef][Medline]

BUSCH, J. W., 2005a The evolution of self-compatibility in geographically peripheral populations of Leavenworthia alabamica (Brassicaceae). Am. J. Bot. 92: 1503–1512.[Abstract/Free Full Text]

BUSCH, J. W., 2005b Inbreeding depression in self-incompatible and self-compatible populations of Leavenworthia alabamica. Heredity 94: 159–165.[CrossRef][Medline]

BUSCH, J. W., and D. J. SCHOEN, 2008 The evolution of self-incompatibility when mates are limiting. Trends Plant Sci. 13: 128–136.[CrossRef][Medline]

CASTRIC, V., and X. VEKEMANS, 2004 Plant self-incompatibility in natural populations: a critical assessment of recent theoretical and empirical advances. Mol. Ecol. 13: 2873–2889.[CrossRef][Medline]

CHARLESWORTH, D., P. AWADALLA, B. K. MABLE and M. H. SCHIERUP, 2000 Population-level studies of multiallelic self-incompatibility, with particular reference to Brassicaceae. Ann. Bot. 85: 227–239.[Abstract/Free Full Text]

CHARLESWORTH, D., C. BARTOLOMÉ, M. H. SCHIERUP and B. K. MABLE, 2003a Haplotype structure of the stigmatic self-incompatibility gene in natural populations of Arabidopsis lyrata. Mol. Biol. Evol. 20: 1741–1753.[Abstract/Free Full Text]

CHARLESWORTH, D., B. K. MABLE, M. H. SCHIERUP, C. BARTOLOMÉ and P. AWADALLA, 2003b Diversity and linkage of genes in the self-incompatibility gene family in Arabidopsis lyrata. Genetics 164: 1519–1535.[Abstract/Free Full Text]

CHARLESWORTH, D., E. KAMAU, J. HAGENBLAD and C. L. TANG, 2006 Trans-specificity at loci near the self-incompatibility loci in Arabidopsis. Genetics 172: 2699–2704.[Abstract/Free Full Text]

DARLINGTON, C. D., and K. MATHER, 1949 The Elements of Genetics. Allen & Unwin, London.

DARWIN, C., 1876 The Effects of Cross and Self-fertilisation in the Vegetable Kingdom. John Murray, London.

FOBIS-LOISY, I., C. MEIGE and T. GAUDE, 2004 Molecular evolution of the S locus controlling mating in the Brassicaceae. Plant Biol. 6: 109–118.[CrossRef][Medline]

FUJIMOTO, R., T. SUGIMURA and T. NISHIO, 2006 Gene conversion from SLG to SRK resulting in self-compatibility in Brassica rapa. FEBS Lett. 580: 425–430.[CrossRef][Medline]

GLÉMIN, S., E. BAZIN and D. CHARLESWORTH, 2006 Impact of mating systems on patterns of sequence polymorphism in flowering plants. Proc. Biol. Sci. 273: 3011–3019.[CrossRef][Medline]

HAGENBLAD, J., J. BECHSGAARD and D. CHARLESWORTH, 2006 Linkage disequilibrium between incompatibility locus region genes in the plant Arabidopsis lyrata. Genetics 173: 1057–1073.[Abstract/Free Full Text]

HINATA, K., M. WATANABE, S. YAMAKAWA, Y. SATTA and A. ISOGAI, 1995 Evolutionary aspects of the S-related genes of the Brassica self-incompatibility system: synonymous and nonsynonymous base substitutions. Genetics 140: 1099–1104.[Abstract]

HISCOCK, S. J., and D. A. TABAH, 2003 The different mechanisms of sporophytic self-incompatibility. Philos. Trans. R. Soc. Lond. B Biol. Sci. 358: 1037–1045.[CrossRef][Medline]

IGIC, B., and J. R. KOHN, 2001 Evolutionary relationships among self-incompatibility RNases. Proc. Natl. Acad. Sci. USA 98: 13167–13171.[Abstract/Free Full Text]

IGIC, B., and J. R. KOHN, 2006 The distribution of plant mating systems: study bias against obligately outcrossing species. Evolution 60: 1098–1103.[Medline]

IGIC, B., L. BOHS and J. R. KOHN, 2006 Ancient polymorphism reveals unidirectional breeding system shifts. Proc. Natl. Acad. Sci. USA 103: 1359–1363.[Abstract/Free Full Text]

KACHROO, A., C. R. SCHOPFER, M. E. NASRALLAH and J. B. NASRALLAH, 2001 Allele-specific receptor-ligand interactions in Brassica self-incompatibility. Science 293: 1824–1826.[Abstract/Free Full Text]

KAMAU, E., B. CHARLESWORTH and D. CHARLESWORTH, 2007 Linkage disequilibrium and recombination rate estimates in the self-incompatibility region of Arabidopsis lyrata. Genetics 176: 2357–2369.[Abstract/Free Full Text]

KHO, Y. O., and J. BAER, 1968 Observing pollen tubes by means of fluorescence. Euphytica 17: 298–302.

KIMURA, M., 1980 A simple method for estimating evolutionary rate of base substitution through comparative studies of nucleotide sequences. J. Mol. Evol. 16: 111–120.[CrossRef][Medline]

KUMAR, V., and M. TRICK, 1993 Sequence complexity of the S receptor kinase gene family in Brassica. Mol. Gen. Genet. 241: 440–445.[CrossRef][Medline]

KUSABA, M., T. NISHIO, Y. SATTA, K. HINATA and D. OCKENDON, 1997 Striking sequence similarity in inter- and intra-specific comparisons of class I SLG alleles from Brassica oleracea and Brassica campestris: implications for the evolution and recognition mechanism. Proc. Natl. Acad. Sci. USA 94: 7673–7678.[Abstract/Free Full Text]

KUSABA, M., K. DWYER, J. HENDERSHOT, J. VREBALOV, J. B. NASRALLAH et al., 2001 Self-incompatibility in the genus Arabidopsis: characterization of the S locus in the outcrossing A. lyrata and its autogamous relative A. thaliana. Plant Cell 13: 627–643.[Abstract/Free Full Text]

LAWRENCE, M. J., 2000 Population genetics of the homomorphic self-incompatibility polymorphisms in flowering plants. Ann. Bot. 85: 221–226.[Abstract/Free Full Text]

LIU, F., L. ZHANG and D. CHARLESWORTH, 1998 Genetic diversity in Leavenworthia populations with different inbreeding levels. Proc. Biol. Sci. 265: 293–301.[CrossRef][Medline]

LIU, F., D. CHARLESWORTH and M. KREITMAN, 1999 The effect of mating system differences on nucleotide diversity at the phosphoglucose isomerase locus in the plant genus Leavenworthia. Genetics 151: 343–357.[Abstract/Free Full Text]

LIU, P., S. SHERMAN-BROYLES, M. E. NASRALLAH and J. B. NASRALLAH, 2007 A cryptic modifier causing transient self-incompatibility in Arabidopsis thaliana. Curr. Biol. 17: 734–740.[CrossRef][Medline]

LLOYD, D. G., 1965 Evolution of self-compatibility and racial differentiation in Leavenworthia (Cruciferae). Contr. Gray Herb. Harv. Univ. 195: 3–134.

LLOYD, D. G., 1967 The genetics of self-incompatibility in Leavenworthia crassa Rollins (Cruciferae). Genetica 38: 227–242.[CrossRef]

LUU, D-T., S. HUGUES, E. PASSELÈGUE and P. HEIZMANN, 2001 Evidence for orthologous S-locus related I genes in several genera of Brassicaceae. Mol. Gen. Genet. 264: 735–745.[CrossRef][Medline]

MABLE, B. K., M. H. SCHIERUP and D. CHARLESWORTH, 2003 Estimating the number, frequency, and dominance of S-alleles in a natural population of Arabidopsis lyrata. Heredity 90: 422–431.[CrossRef][Medline]

MITCHELL-OLDS, T., 2001 Arabidopsis thaliana and its wild relatives: a model system for ecology and evolution. Trends Ecol. Evol. 16: 693–700.[CrossRef]

MURASE, K. H., H. SHIBA, M. IWANO, F. S. CHE, M. WATANABE et al., 2004 A membrane-anchored protein kinase involved in Brassica self-incompatibility signalling. Science 303: 1516–1519.[Abstract/Free Full Text]

NASRALLAH, M. E., P. LIU and J. B. NASRALLAH, 2002 Generation of self-incompatible Arabidopsis thaliana by transfer of two S locus genes from A. lyrata. Science 297: 247–249.[Abstract/Free Full Text]

NASRALLAH, M. E., P. LIU, S. SHERMAN-BROYLES, N. A. BOGGS and J. B. NASRALLAH, 2004 Natural variation in expression of self-incompatibility in Arabidopsis thaliana: implications for the evolution of selfing. Proc. Natl. Acad. Sci. USA 101: 16070–16074.[Abstract/Free Full Text]

NIELSEN, R., and Z. H. YANG, 1998 Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics 148: 929–936.[Abstract/Free Full Text]

OKAMOTA, S., Y. SATO, K. SAKAMOTO and T. NISHIO, 2004 Distribution of similar self-incompatibility (S) haplotypes in different genera, Raphanus and Brassica. Sex Plant Reprod. 17: 33–39.[CrossRef]

PAETSCH, M., S. MAYLAND-QUELLHORST and B. NUEFFER, 2006 Evolution of the self-incompatibility system in the Brassicaceae: identification of S-locus receptor kinase (SRK) in self-incompatible Capsella grandiflora. Heredity 97: 283–290.[CrossRef][Medline]

PRIGODA, N. L., A. NASSUTH and B. K. MABLE, 2005 Phenotypic and genotypic expression of self-incompatibility haplotypes in Arabidopsis lyrata suggests unique origin of alleles in different dominance classes. Mol. Biol. Evol. 22: 1609–1620.[Abstract/Free Full Text]

RICHARDS, A. J., 1986 Plant Breeding Systems. Allen & Unwin, London.

RICHMAN, A. D., and J. R. KOHN, 2000 Evolutionary genetics of self-incompatibility in the Solanaceae. Plant Mol. Biol. 42: 169–179.[CrossRef][Medline]

RICHMAN, A. D., M. K. UYENOYAMA and J. R. KOHN, 1996 Allelic diversity and gene genealogy at the self-incompatibility locus in the Solanaceae. Science 273: 1212–1216.[Abstract]

ROLLINS, R. C., 1963 The evolution and systematics of Leavenworthia (Cruciferae). Contr. Gray Herb. Harv. Univ. 192: 1–98.

ROZAS, J., J. C. SANCHEZ-DELBARRIO, X. MESSEGUER and R. ROZAS, 2003 DnaSP, DNA polymorphism analyses by coalescent and other methods. Bioinformatics 19: 2496–2497.[Abstract/Free Full Text]

SAINUDIIN, R., W. S. W. WONG, K. YOGEESWARAN, J. B. NASRALLAH, Z. H. YANG et al., 2005 Detecting site-specific physiochemical selective pressures: applications to the class IHLA of the human major histocompatibility complex and the SRK of the plant sporophytic self-incompatibility system. J. Mol. Evol. 60: 315–326.[CrossRef][Medline]

SAMPSON, D. R., 1967 Frequency and distribution of self-incompatibility alleles in Raphanus raphanistrum. Genetics 56: 241–251.[Free Full Text]

SCHIERUP, M. H., X. VEKEMANS and F. B. CHRISTIANSEN, 1997 Evolutionary dynamics of sporophytic self-incompatibility alleles in plants. Genetics 147: 835–846.[Abstract]

SCHIERUP, M. H., B. K. MABLE, P. AWADALLA and D. CHARLESWORTH, 2001 Identification and characterization of a polymorphic receptor kinase gene linked to the self-incompatibility locus of Arabidopsis lyrata. Genetics 158: 387–399.[Abstract/Free Full Text]

SCHOPFER, C. R., M. E. NASRALLAH and J. B. NASRALLAH, 1999 The male determinant of self-incompatibility in Brassica. Science 286: 1697–1700.[Abstract/Free Full Text]

SHERMAN-BROYLES, S., N. BOGGS, A. FARKAS, P. LIU, J. VREBALOV et al., 2007 S locus genes and the evolution of self-fertility in Arabidopsis thaliana. Plant Cell 19: 94–106.[Abstract/Free Full Text]

STEBBINS, G. L., 1974 Flowering Plants: Evolution Above the Species Level. Harvard University Press, Cambridge, MA.

STEIN, J. C., B. HOWLETT, D. C. BOYES, M. E. NASRALLAH and J. B. NASRALLAH, 1991 Molecular cloning of a putative receptor kinase gene encoded at the self-incompatibility locus of Brassica oleracea. Proc. Natl. Acad. Sci. USA 88: 8816–8820.[Abstract/Free Full Text]

STEIN, J. C., R. DIXIT, M. E. NASRALLAH and J. B. NASRALLAH, 1996 SRK, the stigma-specific S locus receptor kinase of Brassica, is targeted to the plasma membrane in transgenic tobacco. Plant Cell 8: 429–445.[Abstract]

STEVENS, W. L., 1942 Accuracy of mutation rates. J. Genet. 43: 301–307.

SUZUKI, G., M. WATANABE, N. KAI, N. MATSUDA, K. TORIYAMA et al., 1997 Three members of the S multigene family are linked to the S locus of Brassica. Mol. Gen. Genet. 256: 257–264.[CrossRef][Medline]

SUZUKI, G., N. KAI, T. HIROSE, K. FUKUI, T. NISHIO et al., 1999 Genomic organization of the S locus: identification and characterization of genes in SLG/SRK region of S9 haplotype of Brassica campestris (syn. rapa). Genetics 153: 391–400.[Abstract/Free Full Text]

SUZUKI, T., M. KUSABA, M. MATSUTHITA, K. OKAZAKI and T. NISHIO, 2000 Characterization of Brassica S-haplotypes lacking S-locus glycoprotein. FEBS Lett. 482: 102–108.[CrossRef][Medline]

SWOFFORD, D. L., 2002 PAUP*. Phylogenetic Analysis Using Parsimony, Version 4. Sinauer Associates, Sunderland, MA.

TAKASAKI, T., K. HATAKEYAMA, G. SUZUKI, M. WATANABE, A. ISOGAI et al., 2000 The S receptor kinase determines self-incompatibility in Brassica stigmas. Nature 403: 913–916.[CrossRef][Medline]

TAKAYAMA, S., H. SHIMOSATO, H. SHIBA, M. FUNATO, F. S. CHE et al., 2001 Direct ligand-receptor complex interaction controls Brassica self- incompatibility. Nature 413: 534–538.[CrossRef][Medline]

THOMPSON, K. F., 1957 Self-incompatibility in marrow-stem kale, Brassica oleracea var. acephala. I. Demonstration of a sporophytic system. J. Genet. 55: 45–60.

UYENOYAMA, M. K., 1995 A generalized least-squares estimate for the origin of sporophytic self-incompatibility. Genetics 139: 975–992.[Abstract]

UYENOYAMA, M. K., 2000 Evolutionary dynamics of self-incompatibility alleles in Brassica. Genetics 156: 351–359.[Abstract/Free Full Text]

WALKER, J. C., 1994 Structure and function of the receptor-