Hitch-hiking to a locus under balancing selection: high sequence diversity and low population subdivision at the S-locus genomic region in Arabidopsis 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 plant Arabidopsis 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 with SRK, 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.     

Evolution of the S-locus region in Arabidopsis relatives

The S locus, a single polymorphic locus, is responsible for self-incompatibility (SI) in the Brassicaceae family and many related plant families. Despite its importance, our knowledge of S-locus evolution is largely restricted to the causal genes encoding the S-locus receptor kinase (SRK) receptor and S-locus cysteine-rich protein (SCR) ligand of the SI system. Here, we present high-quality sequences of the genomic region of six S-locus haplotypes: Arabidopsis (Arabidopsis thaliana; one haplotype), Arabidopsis lyrata (four haplotypes), and Capsella rubella (one haplotype). We compared these with reference S-locus haplotypes of the self-compatible Arabidopsis and its SI congener A. lyrata. We subsequently reconstructed the likely genomic organization of the S locus in the most recent common ancestor of Arabidopsis and Capsella. As previously reported, the two SI-determining genes, SCR and SRK, showed a pattern of coevolution. In addition, consistent with previous studies, we found that duplication, gene conversion, and positive selection have been important factors in the evolution of these two genes and appear to contribute to the generation of new recognition specificities. Intriguingly, the inactive pseudo-S-locus haplotype in the self-compatible species C. rubella is likely to be an old S-locus haplotype that only very recently became fixed when C. rubella split off from its SI ancestor, Capsella grandiflora.     

Construction of a binary bacterial artificial chromosome library of Petunia inflata and the isolation of large genomic fragments linked to the self-incompatibility (S-) locus.

The Solanaceae family of flowering plants possesses a type of self-incompatibility mechanism that enables the pistil to reject self pollen but accept non-self pollen for fertilization. The pistil function in this system has been shown to be controlled by a polymorphic gene at the S-locus, termed the S-RNase gene. The pollen function is believed to be controlled by another as yet unidentified polymorphic gene at the S-locus, termed the pollen S-gene. As a first step in using a functional genomic approach to identify the pollen S-gene, a genomic BAC (bacterial artificial chromosome) library of the S2S2 genotype of Petunia inflata, a self-incompatible solanaceous species, was constructed using a Ti-plasmid based BAC vector, BIBAC2. The average insert size was 136.4 kb and the entire library represented a 7.5-fold genome coverage. Screening of the library using cDNAs for the S2-RNase gene and 13 pollen-expressed genes that are linked to the S-locus yielded 51 positive clones, with at least one positive clone for each gene. Collectively, at least 2 Mb of the chromosomal region was spanned by these clones. Together, three clones that contained the S2-RNase gene spanned approximately 263 kb. How this BAC library and the clones identified could be used to identify the pollen S-gene and to study other aspects of self-incompatibility is discussed.     

Molecular characterization of DNA sequences from the Primula vulgaris S-locus

Primula species provide possibly the best known examples of heteromorphic flower development and this breeding system has attracted considerable attention, including that of Charles Darwin. However, despite considerable recent advances in molecular genetics, nothing is known about the molecular basis of floral heteromorphy. The first molecular marker for the Primula S-locus is reported here. This DNA sequence was identified by random amplification of polymorphic DNA (RAPD)-PCR, further defined as a sequence characterized amplified region (SCAR) marker, and subsequently shown to correspond to a restriction fragment length polymorphism (RFLP) that is linked to the thrum allele of the Primula S-locus. The sequence of 8.8 kb of genomic DNA encompassing this thrum-specific RFLP is presented. Analysis of this DNA reveals a highly repetitive sequence structure similar to that found at the S-locus in other species; it also contains sequences similar to elements of a Gypsy-like retrotransposon. The identification of a specific DNA sequence associated with the thrum allele of the Primula S-locus provides the first molecular probe with which to investigate the molecular basis of heteromorphic flower development in Primula.     

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.     

Molecular bases and evolutionary dynamics of self-incompatibility in the Pyrinae (Rosaceae)

The molecular bases of the gametophytic self-incompatibility (GSI) system of species of the subtribe Pyrinae (Rosaceae), such as apple and pear, have been widely studied in the last two decades. The characterization of S-locus genes and of the mechanisms underlying pollen acceptance or rejection have been topics of major interest. Besides the single pistil-side S determinant, the S-RNase, multiple related S-locus F-box genes seem to be involved in the determination of pollen S specificity. Here, we collect and review the state of the art of GSI in the Pyrinae. We emphasize recent genomic data that have contributed to unveiling the S-locus structure of the Pyrinae, and discuss their consistency with the models of self-recognition that have been proposed for Prunus and the Solanaceae. Experimental data suggest that the mechanism controlling pollen-pistil recognition specificity of the Pyrinae might fit well with the collaborative 'non-self' recognition system proposed for Petunia (Solanaceae), whereas it presents relevant differences with the mechanism exhibited by the species of the closely related genus Prunus, which uses a single evolutionarily divergent F-box gene as the pollen S determinant. The possible involvement of multiple pollen S genes in the GSI system of Pyrinae, still awaiting experimental confirmation, opens up new perspectives to our understanding of the evolution of S haplotypes, and of the evolution of S-RNase-based GSI within the Rosaceae family. Whereas S-locus genes encode the players determining self-recognition, pollen rejection in the Pyrinae seems to involve a complex cascade of downstream cellular events with significant similarities to programmed cell death.     

Characterization of expressed genes in the SLL2 region of self-compatible Arabidopsis thaliana.

Self-incompatibility in Brassica species is regulated by a set of S-locus genes: SLG, SRK, and SP11/SCR. In the vicinity of the S-locus genes, several expressed genes, SLL2 and SP2/ClpP, etc., were identified in B. campestris. Arabidopsis thaliana is a self-compatible Brassica relative, and its complete genome has been sequenced. From comparison of the genomic sequences between B. campestris and A. thaliana, microsynteny between gene clusters of Arabidopsis and Brassica SLL2 regions was observed, though the S-locus genes, SLG, SRK, and SP11/SCR were not found in the region of Arabidopsis. Almost all genes predicted in this region of Arabidopsis were expressed in both vegetative and reproductive organs, suggesting that the genes in the SLL2 region might not be related to self-incompatibility. Considering the recent speculation that the S-locus genes were translocated as a single unit between Arabidopsis and Brassica, the translocation might have occurred in the region between the SLL2 and SP7 genes.     

Structure of the male determinant factor for Brassica self-incompatibility

Many flowering plants possess a self-incompatibility system to prevent inbreeding. In Brassica rapa, self/non-self recognition in mating is established through S-haplotype-specific interactions between stigma receptors and S-locus protein 11 (SP11, also called S-locus cysteine-rich protein) that is encoded at the highly polymorphic S-locus. Here we describe the solution structure of the SP11 protein of the S8-haplotype (S8-SP11), which specifically binds to the stigma factor of the same haplotype. It folds into an alpha/beta sandwich structure that resembles those of plant defensins. Residues important for structural integrity are highly conserved among the allelic SP11s, suggesting the existence of a common folding pattern. Structure-based sequence alignment and homology modeling of allelic SP11 identified a hyper-variable (HV) region, which is thought to form a loop that bulges out from the body of the protein that is amenable to solvent exposure. We suggest that the HV region could serve as a specific binding site for the stigma receptor.     

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

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

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.     

Genomic organization of the Papaver rhoeas self-incompatibility S(1) locus

The self-incompatibility (SI) response in Papaver rhoeas depends upon the cognate interaction between a pollen-expressed receptor and a stigmatically expressed ligand. The genes encoding these components are situated within the S-locus. In order for SI to be maintained, the genes encoded by the S-locus must be co-inherited with no recombination between them. Several hypotheses, including sequence heterogeneity and chromosomal position, have been put forward to explain the maintenance of the S-locus in the SI systems of the Brassicaceae and the Solanaceae. A region of the Papaver rhoeas genome encompassing part of the self-incompatibility S(1) locus has been cloned and sequenced. The clone contains the gene encoding the stigmatic component of the response, but does not contain a putative pollen S-gene. The sequence surrounding the S(1) gene contains several diverse repetitive DNA elements. As such, the P. rhoeas S-locus bears similarities to the S-loci of other SI systems. An attempt to localize the P. rhoeas S-locus using fluorescence in situ hybridization (FISH) has also been made. The potential relevance of the findings to mechanisms of recombination suppression is discussed.     

Positional cloning of the s haplotype determining the floral and incompatibility phenotype of the long-styled morph of distylous Turnera subulata

Heterostyly is a plant breeding system occurring in approximately 28 plant families and it has often been used as a model system in plant genetics and evolution. Although heterostyly has been studied for over a century beginning with Charles Darwin, the genes determining floral architecture and incompatibility are still unknown. To identify the genes residing at the S-locus of distylous Turnera subulata, we used a positional cloning strategy and assembled three BAC contigs across the S-locus region. In total, 31 overlapping BAC clones were assembled into contigs 1, 2 and SL. We developed and mapped numerous co-dominant markers from the ends of BAC clones across the S-locus region and assayed X-ray deletion mutants to delimit the region of the contig containing the S-locus. Deletion mapping revealed that a single BAC clone (L22s) within contig-SL contains the s haplotype, while two additional BAC clones (I1 and K15) may contain parts of the dominant S haplotype. Furthermore, we exploited the contigs assembled and investigated the rates of recombination at the S-locus as well as in two regions on either side of the S-locus. We found that recombination rates (estimated in kb/cM) are 2-5 times lower at the S-locus relative to flanking regions, although they are not statistically significant. The present study represents a landmark in the molecular characterization of the S-locus of a heterostylous species. We are now on the verge of identifying the genes that have remained elusive since Darwin's comprehensive study of heterostylous systems more than 125 years ago.     

Contrasted patterns of molecular evolution in dominant and recessive self-incompatibility haplotypes in Arabidopsis

Self-incompatibility has been considered by geneticists a model system for reproductive biology and balancing selection, but our understanding of the genetic basis and evolution of this molecular lock-and-key system has remained limited by the extreme level of sequence divergence among haplotypes, resulting in a lack of appropriate genomic sequences. In this study, we report and analyze the full sequence of eleven distinct haplotypes of the self-incompatibility locus (S-locus) in two closely related Arabidopsis species, obtained from individual BAC libraries. We use this extensive dataset to highlight sharply contrasted patterns of molecular evolution of each of the two genes controlling self-incompatibility themselves, as well as of the genomic region surrounding them. We find strong collinearity of the flanking regions among haplotypes on each side of the S-locus together with high levels of sequence similarity. In contrast, the S-locus region itself shows spectacularly deep gene genealogies, high variability in size and gene organization, as well as complete absence of sequence similarity in intergenic sequences and striking accumulation of transposable elements. Of particular interest, we demonstrate that dominant and recessive S-haplotypes experience sharply contrasted patterns of molecular evolution. Indeed, dominant haplotypes exhibit larger size and a much higher density of transposable elements, being matched only by that in the centromere. Overall, these properties highlight that the S-locus presents many striking similarities with other regions involved in the determination of mating-types, such as sex chromosomes in animals or in plants, or the mating-type locus in fungi and green algae.     

Molecular mechanisms of self-incompatibility in Brassica

In Brassica species, self-incompatibility has been mapped genetically to a single chromosomal location. In this region several closely linked genes have been identified. One of them, S-locus receptor kinase (SRK), determines S haplotype specificity of the stigma and it's the key protein for SI reaction. The role of the S locus glycoprotein (SLG) gene remains unclear. In the last decade approximately 15 additional genes linked to S-locus have been found. Recently, a gene has been identified (SCR) that encodes a small cysteine-rich protein which is a candidate for the pollen ligand. In addition to S locus linked genes there are unlinked SLRgenes (S-locus related genes). In this review, we discuss the role of these genes and the current view on the self-incompatibility mechanism in Brassica.     

Identification of self-incompatibility (S-) locus linked pollen cDNA markers in Petunia inflata.

Solanaceous type self-incompatibility (SI) is controlled by a single polymorphic locus, termed the S-locus. The only gene at the S-locus that has been characterized thus far is the S-RNase gene, which controls pistil function, but not pollen function, in SI interactions between pistil and pollen. One approach to identifying additional genes (including the pollen S-gene, which controls pollen function in SI) at the S-locus and to study the structural organization of the S-locus is chromosome walking from the S-RNase gene. However, the presence of highly repetitive sequences in its flanking regions has made this approach difficult so far. Here, we used RNA differential display to identify pollen cDNAs of Petunia inflata, a self-incompatible solanaceous species, which exhibited restriction fragment length polymorphism (RFLP) for at least one of the three S-haplotypes (S1, S2, and S3) examined. We found that the genes corresponding to 10 groups of pollen cDNAs are genetically tightly linked to the S-RNase gene. These cDNA markers will expedite the mapping and cloning of the chromosomal region of the Solanaceae S-locus by providing multiple starting points.     

Construction of a first genetic map of distylous Turnera and a fine-scale map of the S-locus region

As a prelude to discovery of genes involved in floral dimorphism and incompatibility, a genetic map of distylous Turnera was constructed along with a fine-scale map of the S-locus region. The genetic map consists of 79 PCR-based molecular markers (48 AFLP, 18 RAPD, 9 ISSR, 4 RAMP), 5 isozyme loci, one additional gene, and the S-locus, spanning a total distance of 683.3 cM. The 86 markers are distributed in 5 linkage groups, corresponding to the haploid chromosome number. Molecular markers tightly linked or co-segregating with the S-locus in an initial mapping population of 94 individuals were used to assay an additional 642 progeny to construct a map of the S-locus region. The fine-scale map consists of 2 markers (IS864a and RP45E9) flanking the S-locus at distances of 0.41 and 0.54 cM, respectively, and 3 additional markers (OPK14c, RP45G18, and RP81E18) co-segregating with the S-locus in the total mapping population of 736 individuals. The genetic map constructed will serve as a framework for localization of genes outside the S-locus affecting distyly, while molecular markers of the fine-scale map will be used to initiate chromosome walking to find the genes residing at the S-locus.     

Evaluation of candidate F-box genes for the pollen S of gametophytic self-incompatibility in the Pyrinae (Rosaceae) on the basis of their phylogenomic context

The recent analysis of the S-locus region of apple and Japanese pear, two species of Pyrinae (Rosaceae), suggested multiple and different F-box genes (called SFBBs) as candidates for the male determinant (pollen S) of RNase-based gametophytic self-incompatibility in these two species. Here, we followed a phylogenetic approach to take advantage of the pattern of molecular evolution of the S-locus of Pyrinae in characterizing SFBB homologs belonging to S-haplotypes of apple and three species of Pyrus (European, Japanese, and Chinese pears). Our results suggested that the S-locus region of Pyrinae contains no less than six SFBB members and that its structure seems to be rather conserved between apple and pear species. In accordance with the prevailing theory on S-haplotype evolution, the pollen S is expected to have coevolved with the S-RNase and to show some common features derived from the long-term evolution under frequency-dependent balancing selection, i.e., high sequence diversity, evidence of positive selection, and shared ancestral polymorphisms. Using this conceptual framework, we present evidence that some SFBB genes may be better candidates for pollen S in Pyrinae than others. Overall, the SFBB genes analyzed exhibited much lower sequence diversity than their associated S-RNases; likewise, they showed little or no evidence of positive selection. However, evidence of coevolution with the S-RNase clearly emerged for two of them. Altogether our results suggested different evolutionary histories for different SFBBs putatively derived from their distinct involvement in self-incompatibility.     

S-LOCUS EARLY FLOWERING 3 is exclusively present in the genomes of short-styled buckwheat plants that exhibit heteromorphic self-incompatibility

The different forms of flowers in a species have attracted the attention of many evolutionary biologists, including Charles Darwin. In Fagopyrum esculentum (common buckwheat), the occurrence of dimorphic flowers, namely short-styled and long-styled flowers, is associated with a type of self-incompatibility (SI) called heteromorphic SI. The floral morphology and intra-morph incompatibility are both determined by a single genetic locus named the S-locus. Plants with short-styled flowers are heterozygous (S/s) and plants with long-styled flowers are homozygous recessive (s/s) at the S-locus. Despite recent progress in our understanding of the molecular basis of flower development and plant SI systems, the molecular mechanisms underlying heteromorphic SI remain unresolved. By examining differentially expressed genes from the styles of the two floral morphs, we identified a gene that is expressed only in short-styled plants. The novel gene identified was completely linked to the S-locus in a linkage analysis of 1,373 plants and had homology to EARLY FLOWERING 3. We named this gene S-LOCUS EARLY FLOWERING 3 (S-ELF3). In an ion-beam-induced mutant that harbored a deletion in the genomic region spanning S-ELF3, a phenotype shift from short-styled flowers to long-styled flowers was observed. Furthermore, S-ELF3 was present in the genome of short-styled plants and absent from that of long-styled plants both in world-wide landraces of buckwheat and in two distantly related Fagopyrum species that exhibit heteromorphic SI. Moreover, independent disruptions of S-ELF3 were detected in a recently emerged self-compatible Fagopyrum species and a self-compatible line of buckwheat. The nonessential role of S-ELF3 in the survival of individuals and the prolonged evolutionary presence only in the genomes of short-styled plants exhibiting heteromorphic SI suggests that S-ELF3 is a suitable candidate gene for the control of the short-styled phenotype of buckwheat plants.     

Expression of a self-incompatibility gene in a self-compatible line of Brassica oleracea.

In cruciferous plants, self-pollination is prevented by the action of genes situated at the self-incompatibility locus or S-locus. The self-incompatibility reaction is associated with expression of stigma glycoproteins encoded by the S-locus glycoprotein (SLG) gene. Only a few cases of self-compatible plants derived from self-incompatible lines in the crucifer Brassica have been reported. In these cases, self-compatibility was generally ascribed to the action of single genes unlinked to the S-locus. In contrast, we report here a line of Brassica oleracea var acephala with a self-compatible phenotype linked to the S-locus. By means of both biochemical and immunochemical analyses, we showed that this self-compatible (Sc) line nonetheless possesses stigmatic SLGs (SLG-Sc) that are expressed with a similar spatial and temporal pattern to that described for the SLGs of self-incompatible Brassica plants. Moreover, the SLG-Sc products segregate with the self-compatibility phenotype in F2 progeny, suggesting that changes at the S-locus may be responsible for the occurrence of the self-compatibility character. A cDNA clone encoding the SLG-Sc product was isolated, and the deduced amino acid sequence showed this glycoprotein to be highly homologous to the pollen recessive S2 allele glycoprotein. Hence, self-compatibility in this Brassica Sc line correlates with the expression of a pollen recessive-like S allele in the stigma.