Evidence that intragenic recombination contributes to allelic diversity of the S-RNase gene at the self-incompatibility (S) locus in Petunia inflata

Analyses of Cullin1 homologs reveal functional redundancy in S-RNase-based self-incompatibility and evolutionary relationships in eudicots

In Petunia (Solanaceae family), self-incompatibility (SI) is regulated by the polymorphic S-locus, which contains the pistil-specific S-RNase and multiple pollen-specific S-Locus F-box (SLF) genes. SLFs assemble into E3 ubiquitin ligase complexes known as Skp1-Cullin1-F-box complexes (SCFSLF). In pollen tubes, these complexes collectively mediate ubiquitination and degradation of all nonself S-RNases, but not self S-RNase, resulting in cross-compatible, but self-incompatible, pollination. Using Petunia inflata, we show that two pollen-expressed Cullin1 (CUL1) proteins, PiCUL1-P and PiCUL1-B, function redundantly in SI. This redundancy is lost in Petunia hybrida, not because of the inability of PhCUL1-B to interact with SSK1, but due to a reduction in the PhCUL1-B transcript level. This is possibly caused by the presence of a DNA transposon in the PhCUL1-B promoter region, which was inherited from Petunia axillaris, one of the parental species of Pe. hybrida. Phylogenetic and syntenic analyses of Cullin genes in various eudicots show that three Solanaceae-specific CUL1 genes share a common origin, with CUL1-P dedicated to S-RNase-related reproductive processes. However, CUL1-B is a dispersed duplicate of CUL1-P present only in Petunia, and not in the other species of the Solanaceae family examined. We suggest that the CUL1s involved (or potentially involved) in the SI response in eudicots share a common origin.

Wolfberry genomes and the evolution of Lycium (Solanaceae)

Wolfberry Lycium, an economically important genus of the Solanaceae family, contains approximately 80 species and shows a fragmented distribution pattern among the Northern and Southern Hemispheres. Although several herbaceous species of Solanaceae have been subjected to genome sequencing, thus far, no genome sequences of woody representatives have been available. Here, we sequenced the genomes of 13 perennial woody species of Lycium, with a focus on Lycium barbarum. Integration with other genomes provides clear evidence supporting a whole-genome triplication (WGT) event shared by all hitherto sequenced solanaceous plants, which occurred shortly after the divergence of Solanaceae and Convolvulaceae. We identified new gene families and gene family expansions and contractions that first appeared in Solanaceae. Based on the identification of self-incompatibility related-gene families, we inferred that hybridization hotspots are enriched for genes that might be functioning in gametophytic self-incompatibility pathways in wolfberry. Extremely low expression of LOCULE NUBER (LC) and COLORLESS NON-RIPENING (CNR) orthologous genes during Lycium fruit development and ripening processes suggests functional diversification of these two genes between Lycium and tomato. The existence of additional flowering locus C-like MADS-box genes might correlate with the perennial flowering cycle of Lycium. Differential gene expression involved in the lignin biosynthetic pathway between Lycium and tomato likely illustrates woody and herbaceous differentiation. We also provide evidence that Lycium migrated from Africa into Asia, and subsequently from Asia into North America. Our results provide functional insights into Solanaceae origins, evolution and diversification.

Sequence analysis of the Petunia inflata S-locus region containing 17 S-Locus F-Box genes and the S-RNase gene involved in self-incompatibility

Self-incompatibility in Petunia is controlled by the polymorphic S-locus, which contains S-RNase encoding the pistil determinant and 16-20 S-locus F-box (SLF) genes collectively encoding the pollen determinant. Here we sequenced and assembled approximately 3.1 Mb of the S₂ -haplotype of the S-locus in Petunia inflata using bacterial artificial chromosome clones collectively containing all 17 SLF genes, SLFLike1, and S-RNase. Two SLF pseudogenes and 28 potential protein-coding genes were identified, 20 of which were also found at the S-loci of both the S_6a -haplotype of P. inflata and the S_N -haplotype of self-compatible Petunia axillaris, but not in the S-locus remnants of self-compatible potato (Solanum tuberosum) and tomato (Solanum lycopersicum). Comparative analyses of S-locus sequences of these three S-haplotypes revealed potential genetic exchange in the flanking regions of SLF genes, resulting in highly similar flanking regions between different types of SLF and between alleles of the same type of SLF of different S-haplotypes. The high degree of sequence similarity in the flanking regions could often be explained by the presence of similar long terminal repeat retroelements, which were enriched at the S-loci of all three S-haplotypes and in the flanking regions of all S-locus genes examined. We also found evidence of the association of transposable elements with SLF pseudogenes. Based on the hypothesis that SLF genes were derived by retrotransposition, we identified 10 F-box genes as putative SLF parent genes. Our results shed light on the importance of non-coding sequences in the evolution of the S-locus, and on possible evolutionary mechanisms of generation, proliferation, and deletion of SLF genes.

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.

Distinguishing CPT gene family members and vetting the sequence structure of a putative rubber synthesizing variant in Hevea brasiliensis

cis-Prenyltransferases (cis-PTs) constitute a large family of enzymes conserved during evolution and present in all domains of life. cis-PTs catalyze the cis-1,4-polymerization of isoprene units to generate isoprenoids with carbon skeletons varying from C10 (neryl pyrophosphate) to C > 10,000 (natural rubber). Though the previously reported CPTs in Hevea are designated based on sequence variations, their classification was done mostly by phylogenetic analysis using a mixture of partial as well as full length sequences often excluding the UTRs. In this context an attempt was made to reclassify the CPTs strictly based on their sequence similarity and distinguish the members putatively associated with rubber biosynthesis from the others. Extensive computational analysis was carried out on CPT sequences obtained from public resources and whole genome assemblies of Hevea. Based on the results from BLAST analysis, multiple sequence alignments of protein, nucleotide and untranslated regions, open reading frame analysis, gene prediction analysis and sequence length variations, we conclude that there exists mainly three CPTs namely RubCPT1, RubCPT2 and RubCPT3 putatively associated with rubber biosynthesis in Hevea brasiliensis. The rest were categorised as variants of dehydrodolichyl diphosphate synthase (DHDDS) involved in the synthesis of dolichols having short chain isoprenoids. Analysis of the sequence structure of the most highly expressed RubCPT1 in latex revealed the allele richness and diversity of this important variant prevailing in the popular rubber clones. Haplotypes consisting of SNPs with high degree of heterozygosity were also identified. Segregation and linkage disequilibrium analysis confirmed that recombination is the major contributor towards the generation of allelic diversity rather than point mutations. Alternatively, gene expression analysis indicated the possibility of association between specific haplotypes and RubCPT1 expression in Hevea clones which may have downstream impact up to the level of rubber production. The conclusions from this study may pave way for the identification and better understanding of CPTs directly involved with natural rubber biosynthesis in Hevea and the SNP data generated may aid in the development of molecular markers putatively associated with yield in rubber.

S-Locus F-Box Proteins Are Solely Responsible for S-RNase-Based Self-Incompatibility of Petunia Pollen

Self-incompatibility (SI) in Petunia is regulated by a polymorphic S-locus. For each S-haplotype, the S-locus contains a pistil-specific S-RNase gene and multiple pollen-specific S-locus F-box (SLF) genes. Both gain-of-function and loss-of-function experiments have shown that S-RNase alone regulates pistil specificity in SI. Gain-of-function experiments on SLF genes suggest that the entire suite of encoded proteins constitute the pollen specificity determinant. However, clear-cut loss-of-function experiments must be performed to determine if SLF proteins are essential for SI of pollen. Here, we used CRISPR/Cas9 to generate two frame-shift indel alleles of S₂ -SLF1 (SLF1 of S₂ -haplotype) in S₂S₃ plants of P. inflata and examined the effect on the SI behavior of S₂ pollen. In the absence of a functional S₂-SLF1, S₂ pollen was either rejected by or remained compatible with pistils carrying one of eight normally compatible S-haplotypes. All results are consistent with interaction relationships between the 17 SLF proteins of S₂ -haplotype and these eight S-RNases that had been determined by gain-of-function experiments performed previously or in this work. Our loss-of-function results provide definitive evidence that SLF proteins are solely responsible for SI of pollen, and they reveal their diverse and complex interaction relationships with S-RNases to maintain SI while ensuring cross-compatibility.

CRISPR/Cas9-mediated knockout of PiSSK1 reveals essential role of S-locus F-box protein-containing SCF complexes in recognition of non-self S-RNases during cross-compatible pollination in self-incompatible Petunia inflata

Function of Petunia PiSSK1. Self-incompatibility (SI), an inbreeding-preventing mechanism, is regulated in Petunia inflata by the polymorphic S-locus, which houses multiple pollen-specific S-locus F-box (SLF) genes and a single pistil-specific S-RNase gene. S ₂-haplotype and S ₃-haplotype possess the same 17 polymorphic SLF genes (named SLF1 to SLF17), and each SLF protein produced in pollen is assembled into an SCF (Skp1-Cullin1-F-box) E3 ubiquitin ligase complex. A complete suite of SLF proteins is thought to collectively interact with all non-self S-RNases to mediate their ubiquitination and degradation by the 26S proteasome, allowing cross-compatible pollination. For each SCF^SLF complex, the Cullin1 subunit (named PiCUL1-P) and Skp1 subunit (named PiSSK1), like the F-box protein subunits (SLFs), are pollen-specific, raising the possibility that they also evolved specifically to function in SI. Here we used CRISPR/Cas9-meditated genome editing to generate frame-shift indel mutations in PiSSK1 and examined the SI behavior of a T ₀ plant (S ₂ S ₃) with biallelic mutations in the pollen genome and two progeny plants (S ₂ S ₂) each homozygous for one of the indel alleles and not carrying the Cas9-containing T-DNA. Their pollen was completely incompatible with pistils of seven otherwise-compatible S-genotypes, but fully compatible with pistils of an S ₃ S ₃ transgenic plant in which production of S₃-RNase was completely suppressed by an antisense S ₃-RNase gene, and with pistils of immature flower buds, which produce little S-RNase. These results suggest that PiSSK1 specifically functions in SI and support the hypothesis that SLF-containing SCF complexes are essential for compatible pollination.

Use of Domain-Swapping to Identify Candidate Amino Acids Involved in Differential Interactions between Two Allelic Variants of Type-1 S-Locus F-Box Protein and S3-RNase in Petunia inflata

Petunia inflata possesses a self-incompatibility (SI) mechanism, which involves S-RNase and multiple S-locus F-box (SLF) genes at the polymorphic S-locus. For a given S-haplotype, each SLF is thought to interact with some of its non-self S-RNases, but not with its self S-RNase. In this work, we studied an allelic pair of SLF1, S2-SLF1 and S3-SLF1, which differ in 44 amino acids and show differential interactions with S3-RNase. We first used an in vivo transgenic assay to determine whether four chimeric proteins of S2-SLF1 and S3-SLF1, each with one of the three functional domains swapped, interact with S3-RNase. The results narrowed the candidate amino acids for specific interaction of S2-SLF1 with S3-RNase to the 16 in domain FD3. We then examined seven additional chimeric proteins by dividing FD3 into two subdomains and four mini-domains (A, B, C and D). The results further narrowed the candidate amino acids to four in mini-domain A and four in mini-domain D. Molecular modeling of interactions between S3-RNase and S2-SLF1 revealed that three of these eight are at the interaction surface, and all three are conserved in S1-SLF1 and S6a-SLF1, both of which interact with S3-RNase based on the in vivo transgenic assay. Three of the chimeric proteins were used for the in vivo transgenic assay to determine whether FD3 alone contains the amino acids required for S2-SLF1 to interact with S7-RNase and S13-RNase. The results revealed the diversity and complexity of interactions between SLF proteins and S-RNases.

Overcoming self-incompatibility in grasses: a pathway to hybrid breeding

Allogamous grasses exhibit an effective two-locus gametophytic self-incompatibility (SI) system, limiting the range of breeding techniques applicable for cultivar development. Current breeding methods based on populations are characterized by comparably low genetic gains for important traits such as biomass yield. To implement more efficient breeding schemes, the overall understanding of the SI system is crucial as are the mechanisms involved in the breakdown of SI. Self-fertile variants in outcrossing grasses have been studied, and the current level of knowledge includes approximate gene locations, linked molecular markers and first hypotheses on their mode of action. Environmental conditions increasing seed set upon self-pollination have also been described. Even though some strategies were proposed to take advantage of self-fertility, there have, so far, not been changes in the methods applied in cultivar development for allogamous grasses. In this review, we describe the current knowledge about self-fertility in allogamous grasses and outline strategies to incorporate this trait for implementation in synthetic and hybrid breeding schemes.

Molecular evolution and functional characterisation of haplotypes of an important rubber biosynthesis gene in Hevea brasiliensis

Hydroxy-methylglutaryl coenzyme-A synthase (HMGS) is a rate-limiting enzyme in the cytoplasmic isoprenoid biosynthesis pathway leading to natural rubber production in Hevea brasiliensis (rubber). Analysis of the structural variants of this gene is imperative to understand their functional significance in rubber biosynthesis so that they can be properly utilised for ongoing crop improvement programmes in Hevea. We report here allele richness and diversity of the HMGS gene in selected popular rubber clones. Haplotypes consisting of single nucleotide polymorphisms (SNPs) from the coding and non-coding regions with a high degree of heterozygosity were identified. Segregation and linkage disequilibrium analysis confirmed that recombination is the major contributor to the generation of allelic diversity, rather than point mutations. The evolutionarily conserved nature of some SNPs was identified by comparative DNA sequence analysis of HMGS orthologues from diverse taxa, demonstrating the molecular evolution of rubber biosynthesis genes in general. In silico three-dimensional structural studies highlighting the structural positioning of non-synonymous SNPs from different HMGS haplotypes revealed that the ligand-binding site on the enzyme remains impervious to the reported sequence variations. In contrast, gene expression results indicated the possibility of association between specific haplotypes and HMGS expression in Hevea clones, which may have a downstream impact up to the level of rubber production. Moreover, haplotype diversity of the HMGS gene and its putative association with gene expression can be the basis for further genetic association studies in rubber. Furthermore, the data also show the role of SNPs in the evolution of candidate genes coding for functional traits in plants.

Molecular evolution of candidate male reproductive genes in the brown algal model Ectocarpus

BACKGROUND

Evolutionary studies of genes that mediate recognition between sperm and egg contribute to our understanding of reproductive isolation and speciation. Surface receptors involved in fertilization are targets of sexual selection, reinforcement, and other evolutionary forces including positive selection. This observation was made across different lineages of the eukaryotic tree from land plants to mammals, and is particularly evident in free-spawning animals. Here we use the brown algal model species Ectocarpus (Phaeophyceae) to investigate the evolution of candidate gamete recognition proteins in a distant major phylogenetic group of eukaryotes.

RESULTS

Male gamete specific genes were identified by comparing transcriptome data covering different stages of the Ectocarpus life cycle and screened for characteristics expected from gamete recognition receptors. Selected genes were sequenced in a representative number of strains from distant geographical locations and varying stages of reproductive isolation, to search for signatures of adaptive evolution. One of the genes (Esi0130_0068) showed evidence of selective pressure. Interestingly, that gene displayed domain similarities to the receptor for egg jelly (REJ) protein involved in sperm-egg recognition in sea urchins.

CONCLUSIONS

We have identified a male gamete specific gene with similarity to known gamete recognition receptors and signatures of adaptation. Altogether, this gene could contribute to gamete interaction during reproduction as well as reproductive isolation in Ectocarpus and is therefore a good candidate for further functional evaluation.

Gene duplication and genetic exchange drive the evolution of S-RNase-based self-incompatibility in Petunia

Self-incompatibility (SI) systems in flowering plants distinguish self- and non-self pollen to prevent inbreeding. While other SI systems rely on the self-recognition between specific male- and female-determinants, the Solanaceae family has a non-self recognition system resulting in the detoxification of female-determinants of S-ribonucleases (S-RNases), expressed in pistils, by multiple male-determinants of S-locus F-box proteins (SLFs), expressed in pollen. It is not known how many SLF components of this non-self recognition system there are in Solanaceae species, or how they evolved. We identified 16-20 SLFs in each S-haplotype in SI Petunia, from a total of 168 SLF sequences using large-scale next-generation sequencing and genomic polymerase chain reaction (PCR) techniques. We predicted the target S-RNases of SLFs by assuming that a particular S-allele must not have a conserved SLF that recognizes its own S-RNase, and validated these predictions by transformation experiments. A simple mathematical model confirmed that 16-20 SLF sequences would be adequate to recognize the vast majority of target S-RNases. We found evidence of gene conversion events, which we suggest are essential to the constitution of a non-self recognition system and also contribute to self-compatible mutations.

Supergenes and their role in evolution

Adaptation is commonly a multidimensional problem, with changes in multiple traits required to match a complex environment. This is epitomized by balanced polymorphisms in which multiple phenotypes co-exist and are maintained in a population by a balance of selective forces. Consideration of such polymorphisms led to the concept of the supergene, where alternative phenotypes in a balanced polymorphism segregate as if controlled by a single genetic locus, resulting from tight genetic linkage between multiple functional loci. Recently, the molecular basis for several supergenes has been resolved. Thus, major chromosomal inversions have been shown to be associated with polymorphisms in butterflies, ants and birds, offering a mechanism for localised reduction in recombination. In several examples of plant self-incompatibility, the functional role of multiple elements within the supergene architecture has been demonstrated, conclusively showing that balanced polymorphism can be maintained at multiple coadapted and tightly linked elements. Despite recent criticism, we argue that the supergene concept remains relevant and is more testable than ever with modern molecular methods.

Self-incompatibility in Petunia inflata: the relationship between a self-incompatibility locus F-box protein and its non-self S-RNases

The highly polymorphic S (for self-incompatibility) locus regulates self-incompatibility in Petunia inflata; the S-RNase regulates pistil specificity, and multiple S-locus F-box (SLF) genes regulate pollen specificity. The collaborative non-self recognition model predicts that, for any S-haplotype, an unknown number of SLFs collectively recognize all non-self S-RNases to mediate their ubiquitination and degradation. Using a gain-of-function assay, we examined the relationships between S2-SLF1 (for S2-allelic product of Type-1 SLF) and four S-RNases. The results suggest that S2-SLF1 interacts with S7- and S13-RNases, and the previously identified S1- and S3-RNases, but not with S5- or S11-RNase. An artificial microRNA expressed by the S2-SLF1 promoter, but not by the vegetative cell-specific promoter, Late Anther Tomato 52, suppressed expression of S2-SLF1 in S2 pollen, suggesting that SLF1 is specific to the generative cell. The S2 pollen with S2-SLF1 suppressed was compatible with S3-, S5-, S7-, S11-, and S13-carrying pistils, confirming that other SLF proteins are responsible for detoxifying S5- and S11-RNases and suggesting that S2-SLF1 is not the only SLF in S2 pollen that interacts with S3-, S7-, and S13-RNases. Petunia may have evolved at least two types of SLF proteins to detoxify any non-self S-RNase to minimize the deleterious effects of mutation in any SLF.

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.

Origin and diversification dynamics of self-incompatibility haplotypes

Self-incompatibility (SI) is a genetic system found in some hermaphrodite plants. Recognition of pollen by pistils expressing cognate specificities at two linked genes leads to rejection of self pollen and pollen from close relatives, i.e., to avoidance of self-fertilization and inbred matings, and thus increased outcrossing. These genes generally have many alleles, yet the conditions allowing the evolution of new alleles remain mysterious. Evolutionary changes are clearly necessary in both genes, since any mutation affecting only one of them would result in a nonfunctional self-compatible haplotype. Here, we study diversification at the S-locus (i.e., a stable increase in the total number of SI haplotypes in the population, through the incorporation of new SI haplotypes), both deterministically (by investigating analytically the fate of mutations in an infinite population) and by simulations of finite populations. We show that the conditions allowing diversification are far less stringent in finite populations with recurrent mutations of the pollen and pistil genes, suggesting that diversification is possible in a panmictic population. We find that new SI haplotypes emerge fastest in populations with few SI haplotypes, and we discuss some implications for empirical data on S-alleles. However, allele numbers in our simulations never reach values as high as observed in plants whose SI systems have been studied, and we suggest extensions of our models that may reconcile the theory and data.

Functional gametophytic self-incompatibility in a peripheral population of Solanum peruvianum (Solanaceae)

The transition from self-incompatibility to self-compatibility is a common transition in angiosperms often reported in populations at the edge of species range limits. Geographically distinct populations of wild tomato species (Solanum section Lycopersicon (Solanaceae)) have been described as polymorphic for mating system with both self-incompatible and self-compatible populations. Using controlled pollinations and sequencing of the S-RNase mating system gene, we test the compatibility status of a population of S. peruvianum located near its southern range limit. Pollinations among plants of known genotypes revealed strong self-incompatibility; fruit set following compatible pollinations was significantly higher than following incompatible pollinations for all tested individuals. Sequencing of the S-RNase gene in parents and progeny arrays was also as predicted under self-incompatibility. Molecular variation at the S-RNase locus revealed a diverse set of alleles, and heterozygosity in over 500 genotyped individuals. We used controlled crosses to test the specificity of sequences recovered in this study; in all cases, results were consistent with a unique allelic specificity for each tested sequence, including two alleles sharing 92% amino-acid similarity. Site-specific patterns of selection at the S-RNase gene indicate positive selection in regions of the gene associated with allelic specificity determination and purifying selection in previously characterized conserved regions. Further, there is broad convergence between the present and previous studies in specific amino-acid positions inferred to be evolving under positive selection.

Evolution and diversity of a fungal self/nonself recognition locus

BACKGROUND

Self/nonself discrimination is an essential feature for pathogen recognition and graft rejection and is a ubiquitous phenomenon in many organisms. Filamentous fungi, such as Neurospora crassa, provide a model for analyses of population genetics/evolution of self/nonself recognition loci due to their haploid nature, small genomes and excellent genetic/genomic resources. In N. crassa, nonself discrimination during vegetative growth is determined by 11 heterokaryon incompatibility (het) loci. Cell fusion between strains that differ in allelic specificity at any of these het loci triggers a rapid programmed cell death response.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we evaluated the evolution, population genetics and selective mechanisms operating at a nonself recognition complex consisting of two closely linked loci, het-c (NCU03493) and pin-c (NCU03494). The genomic position of pin-c next to het-c is unique to Neurospora/Sordaria species, and originated by gene duplication after divergence from other species within the Sordariaceae. The het-c pin-c alleles in N. crassa are in severe linkage disequilibrium and consist of three haplotypes, het-c1/pin-c1, het-c2/pin-c2 and het-c3/pin-c3, which are equally frequent in population samples and exhibit trans-species polymorphisms. The absence of recombinant haplotypes is correlated with divergence of the het-c/pin-c intergenic sequence. Tests for positive and balancing selection at het-c and pin-c support the conclusion that both of these loci are under non-neutral balancing selection; other regions of both genes appear to be under positive selection. Our data show that the het-c2/pin-c2 haplotype emerged by a recombination event between the het-c1/pin-c1 and het-c3/pin-c3 approximately 3-12 million years ago.

CONCLUSIONS/SIGNIFICANCE

These results support models by which loci that confer nonself discrimination form by the association of polymorphic genes with genes containing HET domains. Distinct allele classes can emerge by recombination and positive selection and are subsequently maintained by balancing selection and divergence of intergenic sequence resulting in recombination blocks between haplotypes.

Two neutral variants segregating at the gametophytic self-incompatibility locus of European pear (Pyrus communis L.) (Rosaceae, Pyrinae)

Extensive survey of the S-locus diversity of plant species with RNase-based gametophytic self-incompatibility has failed to identify neutral variation segregating within S-allele specificities. Although this is the expected result according to population genetics theory, it conflicts with recent models of S-allele evolution, which suggest that new specificities might arise by a continuous process of subtle changes that individually do not alter the specificity of the S-genes, but whose cumulative effects result in new S-allele functions. Genomic analysis of S-RNase sequences associated with the S(104) (=S(4), =S(b)) allele of European pear (Pyrus communis L.) cultivars yielded two distinct variants (named herein S(104-1) and S(104-2)) that differed at five nucleotide positions within the open reading frame, two of which resulted in changes in the predicted protein sequence. Test-cross experiments indicated that the S-alleles associated with the S(104-1) and S(104-2)RNases exhibit the same pollen and pistil functions, suggesting that they are two neutral variants segregating within the S(104) haplotype of European pear. These allelic forms might represent transitional states in the process of generating new specificities in the species, in accordance with models that predict S-function transition through neutral intermediates. This possibility was further evaluated through the pattern of molecular evolution of functionally distinct European pear S-RNases, which indicated that most recent S-allele diversification in this species proceeded in the absence of adaptive selective pressure.

Sexual recombinants make a significant contribution to epidemics caused by the wheat pathogen Phaeosphaeria nodorum

We conducted a 2-year mark-release-recapture field experiment to quantify the relative contributions of immigration and sexual and asexual reproduction to epidemics of Stagonospora nodorum blotch caused by Phaeosphaeria nodorum. The epidemic was initiated using nine genetically distinct P. nodorum isolates. Infected plants were sampled four times across two growing seasons. In total, 1,286 isolates were recovered and assayed with 10 microsatellite markers and 1 minisatellite marker. The proportion of isolates having multilocus haplotypes (MLHTs) identical to the inoculated isolates decreased steadily from 86% in the first collection to 25% in the fourth collection. The novel isolates that had different MLHTs compared with the marked inoculants originated through immigration and sexual recombination. By the end of the experiment, nearly three-quarters of the novel isolates originated from sexual recombination. Our results indicate that recombinant offspring and airborne immigrant ascospores can make significant contributions to epidemics of Stagonospora nodorum blotch during a growing season.

Molecular and genetic characterization of novel S-RNases from a natural population of Nicotiana alata

Self-incompatibility in the Solanaceae is mediated by S-RNase alleles expressed in the style, which confer specificity for pollen recognition. Nicotiana alata has been successfully used as an experimental model to elucidate cellular and molecular aspects of S-RNase-based self-incompatibility in Solanaceae. However, S-RNase alleles of this species have not been surveyed from natural populations and consequently the S-haplotype diversity is poorly known. Here the molecular and functional characterization of seven S-RNase candidate sequences, identified from a natural population of N. alata, are reported. Six of these candidates, S ( 5 ), S ( 27 ), S ( 70 ), S ( 75 ), S ( 107 ), and S ( 210 ), showed plant-specific amplification in the natural population and style-specific expression, which increased gradually during bud maturation, consistent with the reported S-RNase expression. In contrast, the S ( 63 ) ribonuclease was present in all plants examined and was ubiquitously expressed in different organs and bud developmental stages. Genetic segregation analysis demonstrated that S ( 27 ), S ( 70 ), S ( 75 ), S ( 107 ), and S ( 210 ) alleles were fully functional novel S-RNases, while S ( 5 ) and S ( 63 ) resulted to be non-S-RNases, although with a clearly distinct pattern of expression. These results reveal the importance of performing functional analysis in studies of S-RNase allelic diversity. Comparative phylogenetic analysis of six species of Solanaceae showed that N. alata S-RNases were included in eight transgeneric S-lineages. Phylogenetic pattern obtained from the inclusion of the novel S-RNase alleles confirms that N. alata represents a broad sample of the allelic variation at the S-locus of the Solanaceae.

Evolutionary patterns at the RNase based gametophytic self - incompatibility system in two divergent Rosaceae groups (Maloideae and Prunus)

Background

Within Rosaceae, the RNase based gametophytic self-incompatibility (GSI) system has been studied at the molecular level in Maloideae and Prunus species that have been diverging for, at least, 32 million years. In order to understand RNase based GSI evolution within this family, comparative studies must be performed, using similar methodologies.

Result

It is here shown that many features are shared between the two species groups such as levels of recombination at the S-RNase (the S-pistil component) gene, and the rate at which new specificities arise. Nevertheless, important differences are found regarding the number of ancestral lineages and the degree of specificity sharing between closely related species. In Maloideae, about 17% of the amino acid positions at the S-RNase protein are found to be positively selected, and they occupy about 30% of the exposed protein surface. Positively selected amino acid sites are shown to be located on either side of the active site cleft, an observation that is compatible with current models of specificity determination. At positively selected amino acid sites, non-conservative changes are almost as frequent as conservative changes. There is no evidence that at these sites the most drastic amino acid changes may be more strongly selected.

Conclusions

Many similarities are found between the GSI system of Prunus and Maloideae that are compatible with the single origin hypothesis for RNase based GSI. The presence of common features such as the location of positively selected amino acid sites and lysine residues that may be important for ubiquitylation, raise a number of issues that, in principle, can be experimentally addressed in Maloideae. Nevertheless, there are also many important differences between the two Rosaceae GSI systems. How such features changed during evolution remains a puzzling issue.

'A life or death decision' for pollen tubes in S-RNase-based self-incompatibility

Mate choice is an essential process during sexual plant reproduction, in which self-incompatibility (SI) is widely adopted as an intraspecific reproductive barrier to inhibit self-fertilization by many flowering plants. Genetic studies show that a single polymorphic S-locus, encoding at least two components from both the pollen and pistil sides, controls the discrimination of self and non-self pollen. In the Solanaceae, Plantaginaceae, and Rosaceae, an S-RNase-based SI mechanism is involved in such a discrimination process. Recent studies have provided some important clues to how a decision is made to accept cross pollen or specifically to reject self pollen. In this review, the molecular features of the pistil and pollen S-specificity factors are briefly summarized and then our current knowledge of the molecular control of cross-pollen compatibility (CPC) and self-pollen incompatibility (SPI) responses, respectively, is presented. The possible biochemical mechanisms of the specificity determinant between the pistil and pollen S factors are discussed and a hypothetical S-RNase endosome sorting model is proposed to illustrate the distinct destinies of pollen tubes following compatible and incompatible pollination.

RNase-based gametophytic self-incompatibility evolution: Questioning the hypothesis of multiple independent recruitments of the S-pollen gene

Multiple independent recruitments of the S-pollen component (always an F-box gene) during RNase-based gametophytic self-incompatibility evolution have recently been suggested. Therefore, different mechanisms could be used to achieve the rejection of incompatible pollen in different plant families. This hypothesis is, however, mainly based on the interpretation of phylogenetic analyses, using a small number of divergent nucleotide sequences. In this work we show, based on a large collection of F-box S-like sequences, that the inferred relationship of F-box S-pollen and F-box S-like sequences is dependent on the sequence alignment software and phylogenetic method used. Thus, at present, it is not possible to address the phylogenetic relationship of F-box S-pollen and S-like sequences from different plant families. In Petunia and Malus/Pyrus the putative S-pollen gene(s) show(s) variability patterns different than expected for an S-pollen gene, raising the question of false identification. Here we show that in Petunia, the unexpected features of the putative S-pollen gene are not incompatible with this gene's being the S-pollen gene. On the other hand, it is very unlikely that the Pyrus SFBB-gamma gene is involved in specificity determination.

Identification of major lysine residues of S(3)-RNase of Petunia inflata involved in ubiquitin-26S proteasome-mediated degradation in vitro

S-RNase-based self-incompatibility has been identified in three flowering plant families, including the Solanaceae, and this self/non-self recognition mechanism between pollen and pistil is controlled by two polymorphic genes at the S-locus, S-RNase and S-locus F-box (SLF). S-RNase is produced in the pistil and taken up by pollen tubes in a non-S-haplotype-specific manner. How an allelic product of SLF interacts with self and non-self S-RNases to result in growth inhibition of self pollen tubes is not completely understood. One model predicts that SLF targets non-self S-RNases for ubiquitin/26S proteasome-mediated degradation, thereby only allowing self S-RNase to exert cytotoxic activity inside a pollen tube. To test this model, we studied whether any of the 20 lysine residues in S(3)-RNase of Petunia inflata might be targets for ubiquitination. We identified six lysines near the C-terminus for which mutation to arginine significantly reduced ubiquitination and degradation of the mutant S(3)-RNase, GST:S(3)-RNase (K141-164R) in pollen tube extracts. We further showed that GST:S(3)-RNase (K141-164R) and GST:S(3)-RNase had similar RNase activity, suggesting that their degradation was probably not caused by an ER-associated protein degradation pathway that removes mis-folded proteins. Finally, we showed that PiSBP1 (P. inflata S-RNase binding protein 1), a potential RING-HC subunit of the PiSLF (P. inflata SLF)-containing E3-like complex, could target S-RNase for ubiquitination in vitro. All these results suggest that ubiquitin/26S proteasome-dependent degradation of S-RNase may be an integral part of the S-RNase-based self-incompatibility mechanism.

A tale of two continents: Baker's rule and the maintenance of self-incompatibility in Lycium (Solanaceae)

Over 50 years ago, Baker (1955, 1967) suggested that self-compatible species were more likely than self-incompatible species to establish new populations on oceanic islands. His logic was straightforward and rested on the assumption that colonization was infrequent; thus, mate limitation favored the establishment of self-fertilizing individuals. In support of Baker's rule, many authors have documented high frequencies of self-compatibility on islands, and recent work has solidified the generality of Baker's ideas. The genus Lycium (Solanaceae) has ca. 80 species distributed worldwide, and phylogenetic studies suggest that Lycium originated in South America and dispersed to the Old World a single time. Previous analyses of the S-RNase gene, which controls the stylar component of self-incompatibility, have shown that gametophytically controlled self-incompatibility is ancestral within the genus, making Lycium a good model for investigating Baker's assertions concerning reproductive assurance following oceanic dispersal. Lycium is also useful for investigations of reproductive evolution, given that species vary both in sexual expression and the presence of self-incompatibility. A model for the evolution of gender dimorphism suggests that polyploidy breaks down self-incompatibility, leading to the evolution of gender dimorphism, which arises as an alternative outcrossing mechanism. There is a perfect association of dimorphic gender expression, polyploidy, and self-compatibility (vs. cosexuality, diploidy, and self-incompatibility) among North American Lycium. Although the association between ploidy level and gender expression also holds for African Lycium, to date no studies of mating systems have been initiated in Old World species. Here, using controlled pollinations, we document strong self-incompatibility in two cosexual, diploid species of African Lycium. Further, we sequence the S-RNase gene in 15 individuals from five cosexual, diploid species of African Lycium and recover 24 putative alleles. Genealogical analyses indicate reduced trans-generic diversity of S-RNases in the Old World compared to the New World. We suggest that genetic diversity at this locus was reduced as a result of a founder event, but, despite the bottleneck, self-incompatibility was maintained in the Old World. Maximum-likelihood analyses of codon substitution patterns indicate that positive Darwinian selection has been relatively strong in the Old World, suggesting the rediversification of S-RNases following a bottleneck. The present data thus provide a dramatic exception to Baker's rule, in addition to supporting a key assumption of the Miller and Venable (2000) model, namely that self-incompatibility is associated with diploidy and cosexuality.

Extraordinary intraspecific diversity in oyster sperm bindin

In free-spawning invertebrates sperm-egg incompatibility is a barrier to mating between species, and divergence of gamete recognition proteins (GRPs) can result in reproductive isolation. Of interest are processes that create reproductive protein diversity within species, because intraspecific variants are potentially involved in mate choice and early speciation. Sperm acrosomes of the Pacific oyster Crassostrea gigas contain the protein bindin that bonds sperm to egg during fertilization. Oyster bindin is a single-copy gene encoding a diversity of protein variants. Oyster bindins have a conserved N-terminal region followed by one to five tandem fucose-binding lectin (F-lectin) domains. These repeats have diversified by positive selection at eight sites clustered on the F-lectin's fucose binding face. Additional bindin variants result from recombination in an intron in each F-lectin repeat. Males also express alternatively spliced bindin cDNAs with one to five repeats, but typically translate only one or two isoforms into protein. Thus, positive selection, alternative splicing, and recombination can create thousands of bindin variants within C. gigas. Models of sexual conflict predict high male diversity when females are diverse and sexual conflict is strong. The amount of intraspecific polymorphism in male GRPs may be a consequence of the relative efficiency of local (molecular recognition) and global (electrical, cortical, and physical) polyspermy blocks that operate during fertilization.

Different positively selected sites at the gametophytic self-incompatibility pistil S-RNase gene in the Solanaceae and Rosaceae (Prunus, Pyrus, and Malus)

In this work we perform a comparative study on the location of positively selected sites (those likely responsible for defining specificity differences) at the S-RNase gene, the pistil component of the gametophytic self-incompatibility system. For Plantaginaceae and Rosaceae (Prunus and Pyrus/Malus) this is the first study of this kind. A clear sign of positive selection was observed for 13, 17, and 27 amino acid sites in Solanaceae, Prunus, and Pyrus/Malus, respectively, using two different methodologies. In Plantaginaceae no clear positively selected sites were identified. Possible reasons for this result are discussed. Indirect experimental evidence suggests that the identified positively selected amino acid sites play a role in specificity determination. The percentage of positively selected sites is similar in Solanaceae and Rosaceae but the location of those sites is different.

Studies of self-incompatibility in wild tomatoes: I. S-allele diversity in Solanum chilense (Dun.) Reiche [corrected] (Solanaceae)

We characterized the molecular allelic variation of RNases at the self-incompatibility (SI) locus of Solanum chilense Dun. We recovered 30 S-RNase allele sequences from 34 plants representing a broad geographic sample. This yielded a species-wide estimate of 35 (95% likelihood interval 31-40) S-alleles. We performed crosses to confirm the association with SI function of 10 of the putative S-RNase allele sequences. Results in all cases were consistent with the expectation that these sequences represent functional alleles under single-locus gametophytic SI. We used the allele sequences to conduct an analysis of selection, as measured by the excess of nonsynonymous changes per site, and found evidence for adaptive changes both within the traditionally defined hypervariable regions and downstream, near the 3'-end of the molecule.

Population structure at the S-locus of Sorbus aucuparia L. (Rosaceae: Maloideae)

Low sequence divergence within functional alleles is predicted for the self-incompatibility locus because of strong negative frequency-dependent selection. Nevertheless, sequence variation within functional alleles is essential for current models of the evolution of new mating types. We genotyped the stylar self-incompatibility RNase of 20 Sorbus aucuparia from a population in the Pyrenees mountains of France in order to compare alleles found there to those previously sampled in a Belgian population. Both populations returned 20 different alleles from samples of 20 individuals, providing maximum-likelihood estimates of 24.4 (95% CI 20-34) alleles in each. Ten alleles occurred in both samples. The maximum likelihood (ML) estimate of the overlap in the alleles present in both populations was 16, meaning that an estimated eight alleles are private to each population, and a total of 32 alleles occur across the two populations examined. We used Fisher's (1961) missing plot method to estimate that 40 alleles occur in the species. In accord with population genetics theory, we observed at most one synonymous sequence difference between copies of alleles sampled from the different populations and no variation within populations. Phylogenetic analysis shows that nearly every allele in S. aucuparia arose prior to divergence of this species from members of three different genera of the Rosaceae subfamily, Maloideae. Lack of observable sequence variation within alleles, coupled with the slow pace of allelic relative to taxonomic diversification, implies that finding intermediate stages in the process of new allele creation will be difficult in this group.

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.

Analysis of S-RNase alleles of almond (Prunus dulcis): characterization of new sequences, resolution of synonyms and evidence of intragenic recombination

Cross-compatibility relationships in almond are controlled by a gametophytically expressed incompatibility system partly mediated by stylar RNases, of which 29 have been reported. To resolve possible synonyms and to provide data for phylogenetic analysis, 21 almond S-RNase alleles were cloned and sequenced from SP (signal peptide region) or C1 (first conserved region) to C5, except for the S29 allele, which could be cloned only from SP to C1. Nineteen sequences (S4, S6, S11-S22, S25-S29)) were potentially new whereas S10 and S24 had previously been published but with different labels. The sequences for S16 and S17 were identical to that for S1, published previously; likewise, S15 was identical to S5. In addition, S4 and S20 were identical, as were S13 and S19. A revised version of the standard table of almond incompatibility genotypes is presented. Several alleles had AT or GA tandem repeats in their introns. Sequences of the 23 distinct newly cloned or already published alleles were aligned. Sliding windows analysis of Ka/Ks identified regions where positive selection may operate; in contrast to the Maloideae, most of the region from the beginning of C3 to the beginning of RC4 appeared not to be under positive selection. Phylogenetic analysis indicated four pairs of alleles had "bootstrap" support > 80%: S5/S10, S4/S8, S11/S24, and S3/S6. Various motifs up to 19 residues long occurred in at least two alleles, and their distributions were consistent with intragenic recombination, as were separate phylogenetic analyses of the 5' and 3' sections. Sequence comparison of phylogenetically related alleles indicated the significance of the region between RC4 and C5 in defining specificity.

Gametophytic self-incompatibility in Lycium parishii (Solanaceae): allelic diversity, genealogical structure, and patterns of molecular evolution at the S-RNase locus

We characterized allelic diversity at the locus controlling self-incompatibility (SI) for a population of Lycium parishii (Solanaceae) from Organ Pipe National Monument, Arizona. Twenty-four partial sequences of S-RNase alleles were recovered from 25 individuals. Estimates of allelic diversity range from 23 to 27 alleles and, consistent with expectations for SI, individuals are heterozygous. We compare S-RNase diversity, patterns of molecular evolution, and the genealogical structure of alleles from L. parishii to a previously studied population of its congener L. andersonii. Gametophytic SI is well characterized for Solanaceae and although balancing selection is hypothesized to be responsible for high levels of allelic divergence, the pattern of selection varies depending on the portion of the gene considered. Site-specific models investigating patterns of selection for L. parishii and L. andersonii indicate that positive selection occurs in those regions of the S-RNase gene hypothesized as important to the recognition response, whereas positive selection was not detected for any position within regions previously characterized as conserved. A 10-species genealogy including S-RNases from a pair of congeners from each of five genera in Solanaceae reveals extensive transgeneric evolution of L. parishii S-RNases. Further, within Lycium, the Dn/Ds ratios for pairs of closely related alleles for intraspecific versus interspecific comparisons were not significantly different, suggesting that the S-RNase diversity recovered in these two species was present prior to the speciation event separating them. Despite this, two S-RNases from L. parishii are identical to two previously reported alleles for L. andersonii, suggesting gene flow between these species.

Ancient polymorphism reveals unidirectional breeding system shifts

Loss of complex characters is thought to be irreversible (Dollo's law). However, hypotheses of irreversible evolution are remarkably difficult to test, especially when character transitions are frequent. In such cases, inference of ancestral states, in the absence of fossil evidence, is uncertain and represents the single greatest constraint for reconstructing the evolutionary history of characters. Breeding system character transitions are of particular interest because they affect the amount and distribution of genetic variation within species. Transitions from obligate outcrossing to partial or predominant self-fertilization are thought to represent one of the most common trends in flowering plants. We use the unique molecular genetic properties (manifested as deep persistent polymorphisms) of the locus that enforces outcrossing to demonstrate that its loss is irreversible in the plant family Solanaceae. We argue that current phylogenetic methods of reconstruction are potentially inadequate in cases where ancestral state information is inferred by using only the phylogeny and the distribution of character states in extant taxa. This study shows in a statistical framework that a particular character transition is irreversible, consistent with Dollo's law.

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.

Variability patterns and positively selected sites at the gametophytic self-incompatibility pollen SFB gene in a wild self-incompatible Prunus spinosa (Rosaceae) population

Current models for the generation of new gametophytic self-incompatibility specificities require that neutral variability segregates within specificity classes. Furthermore, one of the models predicts greater ratios of nonsynonymous to synonymous substitutions in pollen than in pistil specificity genes. All models assume that new specificities arise by mutation only. To test these models, 21 SFB (the pollen S-locus) alleles from a wild Prunus spinosa (Rosaceae) population were obtained. For seven of these, the corresponding S-haplotype was also characterized. The SFB data set was also used to identify positively selected sites. Those sites are likely to be the ones responsible for defining pollen specificities. Of the 23 sites identified as being positively selected, 21 are located in the variable (including a new region described here) and hypervariable regions. Little variability is found within specificity classes. There is no evidence for selective sweeps being more frequent in pollen than in pistil specificity genes. The S-RNase and the SFB genes have only partially correlated evolutionary histories. None of the models is compatible with the variability patterns found in the SFB and the S-haplotype data.

Evolution of reproductive proteins from animals and plants

Sexual reproduction is a fundamental biological process common among eukaryotes. Because of the significance of reproductive proteins to fitness, the diversity and rapid divergence of proteins acting at many stages of reproduction is surprising and suggests a role of adaptive diversification in reproductive protein evolution. Here we review the evolution of reproductive proteins acting at different stages of reproduction among animals and plants, emphasizing common patterns. Although we are just beginning to understand these patterns, by making comparisons among stages of reproduction for diverse organisms we can begin to understand the selective forces driving reproductive protein diversity and the functional consequences of reproductive protein evolution.

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.

High-resolution mapping of theSandZloci ofPhalaris coerulescens

Self incompatibility (SI) in Phalaris coerulescens is gametophytically determined by two unlinked multi allelic loci (S and Z). Neither the S nor Z genes have yet been cloned. As part of a map-based cloning strategy, high-resolution maps of the S and Z regions were generated from distorted segregating populations using RFLP probes from wheat, barley, oat, and Phalaris. The S locus was delimited to 0.26 cM with two boundary markers (Xwg811 and Xpsr168) and cosegregated with Xbm2 and Xbcd762. Xbcd266 was the closest marker linked to Z (0.9 cM). A high level of colinearity in the S and Z regions was found in both self-incompatible and -compatible species. The S locus was localized to the subcentromere region of chromosome 1 and the Z locus to the long arm end of chromosome 2. Several rice BAC clones orthologous to the S and Z locus regions were identified. This opens the possibility of using the rice genome sequence data to generate more closely linked markers and identify SI candidate genes. These results add further support to the conservation of gene order in the S and Z regions of the grass genomes.Key words: Phalaris coerulescens, self-incompatibility, distorted segregation, mapping, map-based cloning, synteny mapping.

Genetic mapping and molecular characterization of the self-incompatibility (S) locus in Petunia inflata

Gametophytic self-incompatibility (SI) possessed by the Solanaceae is controlled by a highly polymorphic locus called the S locus. The S locus contains two linked genes, S-RNase, which determines female specificity, and the as yet unidentified pollen S gene, which determines male specificity in SI interactions. To identify the pollen S gene of Petunia inflata, we had previously used mRNA differential display and subtractive hybridization to identify 13 pollen-expressed genes that showed S -haplotype-specific RFLP. Here, we carried out recombination analysis of 1205 F2 plants to determine the genetic distance between each of these S -linked genes and S-RNase. Recombination was observed between four of the genes (3.16, G211, G212, and G221) and S-RNase, whereas no recombination was observed for the other nine genes (3.2, 3.15, A113, A134, A181, A301, G261, X9, and X11). A genetic map of the S locus was constructed, with 3.16 and G221 delimiting the outer limits. None of the observed crossovers disrupted SI, suggesting that all the genes required for SI are contained in the chromosomal region defined by 3.16 and G221. These results and our preliminary chromosome walking results suggest that the S locus is a huge multi-gene complex. Allelic sequence diversity of G221 and 3.16, as well as of 3.2, 3.15, A113, A134 and G261, was determined by comparing two or three alleles of their cDNA and/or genomic sequences. In contrast to S-RNase, all these genes showed very low degrees of allelic sequence diversity in the coding regions, introns, and flanking regions.

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.

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.

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

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

Evaluation of methods for detecting recombination from DNA sequences: empirical data

The performance of 14 different recombination detection methods was evaluated by analyzing several empirical data sets where the presence of recombination has been suggested or where recombination is assumed to be absent. In general, recombination methods seem to be more powerful with increasing levels of divergence, but different methods showed distinct performance. Substitution methods using summary statistics gave more accurate inferences than most phylogenetic methods. However, definitive conclusions about the presence of recombination should not be derived on the basis of a single method. Performance patterns observed from the analysis of real data sets coincided very well with previous computer simulation results. Previous recombination inferences from some of the data sets analyzed here should be reconsidered. In particular, recombination in HIV-1 seems to be much more widespread than previously thought. This finding might have serious implications on vaccine development and on the reliability of previous inferences of HIV-1 evolutionary history and dynamics.