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

How Have Self-Incompatibility Haplotypes Diversified? Generation of New Haplotypes during the Evolution of Self-Incompatibility from Self-Compatibility

I developed a gametophytic self-incompatibility (SI) model to study the conditions leading to diversification in SI haplotypes. In the model, the SI system is assumed to be incomplete, and the pollen expressing a given specificity is not fully rejected by the pistils expressing the same specificity. I also assumed that mutations can occur that enhance the rejection of pollen by pistils with the same haplotype variant and reduce rejection by pistils with other variants in the same haplotype. I found that if such mutations occur, the new haplotypes (mutant variants) can stably coexist with the ancestral haplotype in which the mutant arose. This is because pollen bearing the new haplotype is most strongly rejected by pistils bearing the same new haplotype among the pistils in the population; hence, negative frequency-dependent selection prevents their fixation. I also performed simulations and found that the nearly complete SI system evolves from completely self-compatible populations and that SI haplotypes can increase to about 40-50 within a few thousand generations. On the basis of my findings, I propose that diversification of SI haplotypes occurred during the evolution of SI from self-compatibility.

Conserved microstructure of the Brassica B Genome of Brassica nigra in relation to homologous regions of Arabidopsis thaliana, B. rapa and B. oleracea

BACKGROUND

The Brassica B genome is known to carry several important traits, yet there has been limited analyses of its underlying genome structure, especially in comparison to the closely related A and C genomes. A bacterial artificial chromosome (BAC) library of Brassica nigra was developed and screened with 17 genes from a 222 kb region of A. thaliana that had been well characterised in both the Brassica A and C genomes.

RESULTS

Fingerprinting of 483 apparently non-redundant clones defined physical contigs for the corresponding regions in B. nigra. The target region is duplicated in A. thaliana and six homologous contigs were found in B. nigra resulting from the whole genome triplication event shared by the Brassiceae tribe. BACs representative of each region were sequenced to elucidate the level of microscale rearrangements across the Brassica species divide.

CONCLUSIONS

Although the B genome species separated from the A/C lineage some 6 Mya, comparisons between the three paleopolyploid Brassica genomes revealed extensive conservation of gene content and sequence identity. The level of fractionation or gene loss varied across genomes and genomic regions; however, the greatest loss of genes was observed to be common to all three genomes. One large-scale chromosomal rearrangement differentiated the B genome suggesting such events could contribute to the lack of recombination observed between B genome species and those of the closely related A/C lineage.

The population genetics of sporophytic self-incompatibility in three hybridizing senecio (asteraceae) species with contrasting population histories

Hybridization generates evolutionary novelty and spreads adaptive variation. By promoting outcrossing, plant self-incompatibility (SI) systems also favor interspecific hybridization because the S locus is under strong negative frequency-dependent balancing selection. This study investigates the SI mating systems of three hybridizing Senecio species with contrasting population histories. Senecio aethnensis and S. chrysanthemifolius native to Sicily, form a hybrid zone at intermediate altitudes on Mount Etna, and their neo-homoploid hybrid species, S. squalidus, has colonized disturbed urban habitats in the UK during the last 150 years. We show that all three species express sporophytic SI (SSI), where pollen incompatibility is controlled by the diploid parental genome, and that SSI is inherited and functions normally in hybrids. Large-scale crossing studies of wild sampled populations allowed direct comparison of SSI between species and found that the main impacts of colonization in S. squalidus compared to Sicilian Senecio was a reduced number of S alleles, increased S allele frequencies, and increased interpopulation S allele sharing. In general, many S alleles were shared between species and the S locus showed reduced intra- and interspecific population genetic structure compared to molecular genetic markers, indicative of enhanced effective gene flow due to balancing selection.

Learning from rejection: the evolutionary biology of single-locus incompatibility

The self-incompatibility (S-) locus of flowering plants is among the most polymorphic known. PCR methods can now be used to estimate both the number of alleles in natural populations and their sequence diversity. The number of alleles provides an estimate of recent effective population size, thus the S-locus provides a tool for examining how species characteristics affect population size. Sequence relationships among alleles provide another estimate of population size extending millions of years into the past. Relationships between S-alleles and related genes provide a means of dating the age of origin of incompatibility systems and determining which, if any, angiosperm families share incompatibility by homology.

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.

Self-incompatibility is a common genetic system preventing selfing through recognition and rejection of self-pollen in hermaphroditic flowering plants. In the Brassicaceae family, this system is controlled by a single genomic region, called the S-locus, where many distinct specificities segregate in natural populations. In this study, we obtained genomic sequences comprising the S-locus in two closely related Brassicaceae species, Arabidopsis lyrata and A. halleri, and analyzed their diversity and patterns of molecular evolution. We report compelling evidence that the S-locus presents many similar properties with other genomic regions involved in the determination of mating-types in mammals, insects, plants, or fungi. In particular, in spite of their diversity, these genomic regions all show absence of similarity in intergenic sequences, large depth of genealogies, highly divergent organization, and accumulation of transposable elements. Moreover, some of these features were found to vary according to dominance of the S-locus specificities, suggesting that dominance/recessivity interactions are key drivers of the evolution of this 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.

Prediction of white cabbage (Brassica oleracea var. capitata) self-incompatibility based on neural network and discriminant analysis of complex electrophoretic patterns

Self-incompatibility (SI) of white cabbage (Brassica oleracea var. capitata), a biological mechanism triggering the ability to pollinate plant ovules, is modelled based on the statistical analysis of electrophoretic patterns of stigma extracts. The patterns obtained for 72 plants grown in three different seasons (2003-2005) have been explored with discriminant analysis and non-linear neural networks. The NN models obtained perform a flawless classification of SI and provide a useful and fast technique for SI prediction before the end of the season. The discriminant analysis turns out to be less effective (74% of good predictions), but together with sensitivity analysis of NN, it points out the important markers of the SI process (peaks 1, 3, 5, 15 and 18).

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

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

The Prunus self-incompatibility locus (S locus) is seldom rearranged

Self-incompatibility enables flowering plants to discriminate between self- and non-selfpollen. In Prunus, the 2 genes determining specificity are the S-RNase (the female determinant that is a glycoprotein with ribonuclease activity) and the SFB (the male determinant, a protein with an F-box motif). In all Prunus S haplotypes characterized so far, with the exception of Prunus armeniaca S(2) haplotype, the 2 genes have opposite transcription orientations. Nevertheless, the relative transcription orientation observed in P. armeniaca S(2) haplotype has been postulated to be the one present in all S haplotypes from this species. We show that this is not the case by demonstrating that that the relative transcription orientation of the pollen and pistil genes of the P. armeniaca S(17) haplotype is that which is commonly found in Prunus. Using a phylogenetic approach, we show that the relative transcription orientation of the S-RNase and SFB genes is seldom changed (less than once every 380 million years). This contrasts with the Brassica sporophytic S locus where chromosomal rearrangements are often observed in the region between the pollen and pistil genes.

Molecular characterization of Lal2, an SRK-like gene linked to the S-locus in the wild mustard Leavenworthia alabamica

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

Genomic consequences of selection on self-incompatibility genes

Frequency-dependent selection at plant self-incompatibility systems is inherent and well understood theoretically. A self-incompatibility locus leads to a strong peak of diversity in the genome, to a unique distribution of diversity across the species and possibly to increased introgression between closely related species. We review recent empirical studies demonstrating these features and relate the empirical findings to theoretical predictions. We show how these features are being exploited in searches for other genes under multi-allelic balancing selection and for inference on recent breakdown of self-incompatibility.

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

Background

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

Results

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

Conclusion

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

Patch Aging and theS‐Allee Effect: Breeding System Effects on the Demographic Response of Plants to Habitat Fragmentation

We used empirical and modeling approaches to examine effects of plant breeding systems on demographic responses to habitat fragmentation. Empirically, we investigated effects of local flowering plant density on pollination and of population size on mate availability in a common, self-incompatible purple coneflower, Echinacea angustifolia, growing in fragmented prairie habitat. Pollination and recruitment increased with weighted local density around individual flowering plants. This positive density dependence is an Allee effect. In addition, mean mate compatibility between pairs of plants increased with population size. Based on this empirical study, we developed an individual-based, spatially explicit demographic model that incorporates autosomal loci and an S locus. We simulated habitat fragmentation in populations identical except for their breeding system, self-incompatible (SI) or self-compatible (SC). Both populations suffered reduced reproduction in small patches because of scarcity of plants within pollination distance (potential mates) and inbreeding depression. But SI species experienced an additional, genetic contribution to the Allee effect (S-Allee effect) caused by allele loss at the S locus, which reduces mate availability, thereby decreasing reproduction. The strength of the S-Allee effect increases through time (i.e., patches age) because random genetic drift reduces S-allele richness. We investigate how patch aging influences extinction and discuss how the S-Allee effect influences communities in fragmented habitat.

Genetic diversity and the reproductive system in related species of antirrhinum

BACKGROUND AND AIMS

Seven related species of Antirrhinum (A. siculum, A. majus, A. latifolium, A. linkianum, A. litigiosum, A. cirrhigherum and A. tortuosum) were studied in order to compare levels of genetic variation and its partitioning in them, and to check relationships between genetic patterns and the reproductive system.

METHODS

Eight hundred and fifty-one plants were screened for variability at 13 allozyme loci by means of horizontal starch gel electrophoresis. Parameters of genetic diversity and its partitioning, the inbreeding coefficient as well as an indirect estimate of gene flow based on the equation: Nm = (1 - G(ST))/4G(ST), were calculated.

KEY RESULTS

Genetic variability in A. siculum was found to be the lowest known in the genus. Mean values of F(IT) and F(IS) were mostly positive and not significantly different from zero. Population differentiation (F(ST)) ranged between 6.1 in A. tortuosum and 17.6 in A. linkianum. The inbreeding coefficient within populations ranged between F(IS) = -0.5 in A. tortuosum and F(IS) = 1 in A. siculum. Estimates of gene flow ranged between Nm = 15 in A. majus (considered as very high) to Nm = 0.42 in A. siculum (considered as low).

CONCLUSIONS

Correlation was found between levels of diversity and differentiation on one hand, and the reproductive system of the studied taxa on the other. Striking differences among species in the inbreeding coefficient (F(IS)) show different reproductive systems, which mostly support previous reports. Strategies for the conservation of A. siculum are recommended, such as preservation of natural populations as well as ex situ preservation of seeds from different populations.

Nonneutral evolution of organelle genes in Silene vulgaris

Knowledge of mitochondrial gene evolution in angiosperms has taken a dramatic shift within the past decade, from universal slow rates of nucleotide change to a growing realization of high variation in rates among lineages. Additionally, evidence of paternal inheritance of plant mitochondria and recombination among mitochondrial genomes within heteroplasmic individuals has led to speculation about the potential for independent evolution of organellar genes. We report intraspecific mitochondrial and chloroplast sequence variation in a cosmopolitan sample of 42 Silene vulgaris individuals. There was remarkably high variation in two mitochondrial genes (atp1, atp9) and additional variation within a third gene (cob). Tests for patterns of nonneutral evolution were significant for atp1 and atp9, indicative of the maintenance of balanced polymorphisms. Two chloroplast genes (matK, ndhF) possessed less, but still high, variation and no divergence from neutral expectations. Phylogenetic patterns of organelle genes in both the chloroplast and mitochondria were incongruent, indicating the potential for independent evolutionary trajectories. Evidence indicated reassociation among cytoplasmic genomes and recombination between mitochondrial genes and within atp1, implying transient heteroplasmy in ancestral lineages. Although the mechanisms for long-term maintenance of mitochondrial polymorphism are currently unknown, frequency-dependent selection on linked cytoplasmic male sterility genes is a potential candidate.

Analyses of synteny between Arabidopsis thaliana and species in the Asteraceae reveal a complex network of small syntenic segments and major chromosomal rearrangements

Comparative genomic studies among highly divergent species have been problematic because reduced gene similarities make orthologous gene pairs difficult to identify and because colinearity is expected to be low with greater time since divergence from the last common ancestor. Nevertheless, synteny between divergent taxa in several lineages has been detected over short chromosomal segments. We have examined the level of synteny between the model species Arabidopsis thaliana and species in the Compositae, one of the largest and most diverse plant families. While macrosyntenic patterns covering large segments of the chromosomes are not evident, significant levels of local synteny are detected at a fine scale covering segments of 1-Mb regions of A. thaliana and regions of <5 cM in lettuce and sunflower. These syntenic patches are often not colinear, however, and form a network of regions that have likely evolved by duplications followed by differential gene loss.

Molecular clocks: four decades of evolution

During the past four decades, the molecular-clock hypothesis has provided an invaluable tool for building evolutionary timescales, and has served as a null model for testing evolutionary and mutation rates in different species. Molecular clocks have also influenced the development of theories of molecular evolution. As DNA-sequencing technologies have progressed, the use of molecular clocks has increased, with a profound effect on our understanding of the temporal diversification of species and genomes.

Detecting site-specific physicochemical selective pressures: applications to the Class I HLA of the human major histocompatibility complex and the SRK of the plant sporophytic self-incompatibility system

Models of codon substitution are developed that incorporate physicochemical properties of amino acids. When amino acid sites are inferred to be under positive selection, these models suggest the nature and extent of the physicochemical properties under selection. This is accomplished by first partitioning the codons on the basis of some property of the encoded amino acids. This partition is used to parametrize the rates of property-conserving and property-altering base substitutions at the codon level by means of finite mixtures of Markov models that also account for codon and transition:transversion biases. Here, we apply this method to two positively selected receptors involved in ligand-recognition: the class I alleles of the human major histocompatibility complex (MHC) of known structure and the S-locus receptor kinase (SRK) of the sporophytic self-incompatibility system (SSI) in cruciferous plants (Brassicaceae), whose structure is unknown. Through likelihood ratio tests we demonstrate that at some sites, the positively selected MHC and SRK proteins are under physicochemical selective pressures to alter polarity, volume, polarity and/or volume, and charge to various extents. An empirical Bayes approach is used to identify sites that may be important for ligand recognition in these proteins.

Chromosome triplication found across the tribe Brassiceae

We have used an approximately 8.7-Mb BAC contig of Arabidopsis thaliana Chromosome 4 to trace homeologous chromosome regions in 21 species of the family Brassicaceae. Homeologs of this segment could be identified in all tested species. Painting of pachytene chromosomes of Calepina, Conringia, and Sisymbrium species (2n = 14, 16), traditionally placed in tribe Brassiceae, showed one homeologous copy of the Arabidopsis contig, while the remaining taxa of the tribe (2n = 14-30) revealed three, and three Brassica species (2n = 34, 36, and 38) and Erucastrum gallicum (2n = 30) had six copies corresponding to the 8.7-Mb segment. The multiple homeologous copies corresponded structurally to the Arabidopsis segment or were rearranged by inversions and translocations within the diploidized genomes. These chromosome rearrangements accompanied by chromosome fusions/fissions led to the present-day chromosome number variation within the Brassiceae. Phylogenetic relationships based on the chloroplast 5'-trnL (UAA)-trnF(GAA) region and estimated divergence times based on sequence data of the chalcone synthase gene are congruent with comparative painting data and place Calepina, Conringia, and Sisymbrium outside the clade of Brassiceae species with triplicated genomes. Most likely, species containing three or six copy pairs descended from a common hexaploid ancestor with basic genomes similar to that of Arabidopsis. The presumed hexaploidization event occurred after the Arabidopsis-Brassiceae split, between 7.9 and 14.6 Mya.

Maximum-likelihood estimation of rates of recombination within mating-type regions

Features common to many mating-type regions include recombination suppression over large genomic tracts and cosegregation of genes of various functions, not necessarily related to reproduction. Model systems for homomorphic self-incompatibility (SI) in flowering plants share these characteristics. We introduce a method for the exact computation of the joint probability of numbers of neutral mutations segregating at the determinant of mating type and at a linked marker locus. The underlying Markov model incorporates strong balancing selection into a two-locus coalescent. We apply the method to obtain a maximum-likelihood estimate of the rate of recombination between a marker locus, 48A, and S-RNase, the determinant of SI specificity in pistils of Nicotiana alata. Even though the sampled haplotypes show complete allelic linkage disequilibrium and recombinants have never been detected, a highly significant deficiency of synonymous substitutions at 48A compared to S-RNase suggests a history of recombination. Our maximum-likelihood estimate indicates a rate of recombination of perhaps 3 orders of magnitude greater than the rate of synonymous mutation. This approach may facilitate the construction of genetic maps of regions tightly linked to targets of strong balancing selection.

Molecular mechanisms of self-incompatibility in flowering plants and fungi - different means to the same end

Hermaphrodite flowering plants and fungi face the same sexual dilemma - how to avoid self-fertilization. Both have evolved ingenious recognition systems that reduce or eliminate the possibility of selfing. These self-incompatibility (SI) systems offer unique opportunities to study recognition and signalling in non-animal cells and also represent model systems for studying the evolution of breeding systems at a molecular level. In this review, the authors discuss recent molecular data that predict an astonishing diversity in the cellular mechanisms of SI operating in flowering plants and fungi.

Evolution under tight linkage to mating type

Recent large-scale sequencing studies of mating type loci in a number of organisms offer insight into the origin and evolution of these genomic regions. Extensive tracts containing genes with a wide diversity of functions typically cosegregate with mating type. Cases in which mating type determination entails complementarity between distinct transcription units may descend from systems in which close physical linkage facilitated the coordinated expression and cosegregation of the interacting genes. In response to the particular selection pressures associated with the maintenance of more than one mating type, this nucleus of low recombination may expand over evolutionary time, engulfing neighboring tracts bearing genes with no direct role in reproduction. This scenario is consistent with the present-day structure of some mating type loci, including regulators of homomorphic self-incompatibility in angiosperms (S-loci). Recombination suppression and enforced S-locus heterozygosity accelerate the accumulation of genetic load and promote genetic associations between S-alleles and degenerating genes in cosegregating tracts. This S-allele-specific load may influence the evolution of self-incompatibility systems.

Molecular evidence on plant divergence times

Estimation of divergence times from sequence data has become increasingly feasible in recent years. Conflicts between fossil evidence and molecular dates have sparked the development of new methods for inferring divergence times, further encouraging these efforts. In this paper, available methods for estimating divergence times are reviewed, especially those geared toward handling the widespread variation in rates of molecular evolution observed among lineages. The assumptions, strengths, and weaknesses of local clock, Bayesian, and rate smoothing methods are described. The rapidly growing literature applying these methods to key divergence times in plant evolutionary history is also reviewed. These include the crown group ages of green plants, land plants, seed plants, angiosperms, and major subclades of angiosperms. Finally, attempts to infer divergence times are described in the context of two very different temporal settings: recent adaptive radiations and much more ancient biogeographic patterns.

Plant self-incompatibility in natural populations: a critical assessment of recent theoretical and empirical advances

Self-incompatibility systems in plants are genetic systems that prevent self-fertilization in hermaphrodites through recognition and rejection of pollen expressing the same allelic specificity as that expressed in the pistils. The evolutionary properties of these self-recognition systems have been revealed through a fascinating interplay between empirical advances and theoretical developments. In 1939, Wright suggested that the main evolutionary force driving the genetic and molecular properties of these systems was strong negative frequency-dependent selection acting on pollination success. The empirical observation of high allelic diversity at the self-incompatibility locus in several species, followed by the discovery of very high molecular divergence among alleles in all plant families where the locus has been identified, supported Wright's initial theoretical predictions as well as many of its later developments. In the last decade, however, advances in the molecular characterization of the incompatibility reaction and in the analysis of allelic frequencies and allelic divergence from natural populations have stimulated new theoretical investigations that challenged some important assumptions of Wright's model of gametophytic self-incompatibility. We here review some of these recent empirical and theoretical advances that investigated: (i) the hypothesis that S-alleles are selectively equivalent, and the evolutionary consequences of genetic interactions between alleles; (ii) the occurrence of frequency-dependent selection in female fertility; (iii) the evolutionary genetics of self-incompatibility systems in subdivided populations; (iv) the evolutionary implications of the self-incompatibility locus's genetic architecture; and (v) of its interactions with the genomic environment.

Molecular evolution of the s locus controlling mating in the brassicaceae

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

Specificity determinants and diversification of the Brassica self-incompatibility pollen ligand

Self-incompatibility in crucifers is effected by allele-specific interactions between the highly polymorphic stigmatic S locus receptor kinase (SRK) and its pollen ligand, the S locus cysteine-rich protein (SCR). Here we show that specificity in SCR function is determined by four contiguous amino acids in one variant, indicating that the minimum sequence requirement for gaining a new specificity can be low. We also provide evidence for an extraordinarily high degree of evolutionary flexibility in SCR, whereby SCR can tolerate extensive amino acid changes within the limits of maintaining the same predicted overall structure. This remarkable adaptability suggests a hypothesis for generation of new self-incompatibility specificities by gradual modification of SRK-SCR affinities and, more generally, for functional specialization within families of homologous ligands and receptors.

Effects of inbreeding on the genetic diversity of populations

The study of variability within species is important to all biologists who use genetic markers. Since the discovery of molecular variability among normal individuals, data have been collected from a wide range of organisms, and it is important to understand the major factors affecting diversity levels and patterns. Comparisons of inbreeding and outcrossing populations can contribute to this understanding, and therefore studying plant populations is important, because related species often have different breeding systems. DNA sequence data are now starting to become available from suitable plant and animal populations, to measure and compare variability levels and test predictions.

Genealogy-dependent variation in viability among self-incompatibility genotypes

Many hermaphroditic plants avoid self-fertilization by rejecting pollen that express genetically determined specificities in common with the pistil. The S-locus, comprising the determinants of pistil and pollen specificity, typically shows extremely high polymorphism, with dozens to hundreds of specificities maintained within species. This article explores a conjecture, motivated by empirical findings, that the expression of recessive deleterious factors at sites closely linked to the S-locus may cause greater declines in the viability of zygotes constituted from more closely related S-alleles. Diffusion approximation models incorporating variation in viability among S-locus genotypes and antagonistic interactions between a new specificity and its immediate parent specificity are constructed and analyzed. Results indicate that variation in viability tends to reduce the number of specificities maintained in a population at stochastic steady state, and that genealogy-based antagonism reduces the rate of bifurcation of S-allele lineages. These effects may account for some of the unusual features observed in empirical studies of S-allele genealogies.

Genomic clocks and evolutionary timescales

For decades, molecular clocks have helped to illuminate the evolutionary timescale of life, but now genomic data pose a challenge for time estimation methods. It is unclear how to integrate data from many genes, each potentially evolving under a different model of substitution and at a different rate. Current methods can be grouped by the way the data are handled (genes considered separately or combined into a 'supergene') and the way gene-specific rate models are applied (global versus local clock). There are advantages and disadvantages to each of these approaches, and the optimal method has not yet emerged. Fortunately, time estimates inferred using many genes or proteins have greater precision and appear to be robust to different approaches.

Just how complex is the Brassica S-receptor complex?

Of the plant self-incompatibility (SI) systems investigated to date, that possessed by members of the Brassicaceae is currently the best understood. Whilst the recent demonstrations of interactions between the male determinant (S-locus cysteine rich protein, SCR) and the female determinant (S-locus receptor kinase, SRK) indicate the minimal requirement for SI in Brassica, no consensus exists as to the nature of these molecules in vivo and the potential involvement of accessory molecules in establishing the active S-receptor complex. Variation between S haplotypes appears to be present in the molecular composition of the receptor complex, the regulation of downstream signalling and the requirement for accessory molecules. This review discusses what constitutes an active receptor complex and highlights potential differences between haplotypes. The role of accessory molecules, in particular SLG (S-locus glycoprotein) and low molecular weight pollen coat proteins (PCPs), in pollination are discussed, as is the link between SI and unilateral incompatibility (UI).

Mathematical consequences of the genealogical species concept

A genealogical species is defined as a basal group of organisms whose members are all more closely related to each other than they are to any organisms outside the group ("exclusivity"), and which contains no exclusive group within it. In practice, a pair of species is so defined when phylogenies of alleles from a sample of loci shows them to be reciprocally monophyletic at all or some specified fraction of the loci. We investigate the length of time it takes to attain this status when an ancestral population divides into two descendant populations of equal size with no gene exchange, and when genetic drift and mutation are the only evolutionary forces operating. The number of loci used has a substantial effect on the probability of observing reciprocal monophyly at different times after population separation, with very long times needed to observe complete reciprocal monophyly for a large number of loci. In contrast, the number of alleles sampled per locus has a relatively small effect on the probability of reciprocal monophyly. Because a single mitochondrial or chloroplast locus becomes reciprocally monophyletic much faster than does a single nuclear locus, it is not advisable to use mitochondrial and chloroplast DNA to recognize genealogical species for long periods after population divergence. Using a weaker criterion of assigning genealogical species status when more than 50% of sampled nuclear loci show reciprocal monophyly, genealogical species status depends much less on the number of sampled loci, and is attained at roughly 4-7 N generations after populations are isolated, where N is the historically effective population size of each descendant. If genealogical species status is defined as more than 95% of sampled nuclear loci showing reciprocal monophyly, this status is attained after roughly 9-12 N generations.

Recognition and rejection of self in plant reproduction

Plant self-incompatibility (SI) systems are unique among self/nonself recognition systems in being based on the recognition of self rather than nonself. SI in crucifer species is controlled by highly polymorphic and co-evolving genes linked in a complex. Self recognition is based on allele-specific interactions between stigma receptors and pollen ligands that result in the arrest of pollen tube development. Commonalities and differences between SI and other self/nonself discrimination systems are discussed.

Molecular mechanisms underlying the breakdown of gametophytic self-incompatibility

The breakdown of self-incompatibility has occurred repeatedly throughout the evolution of flowering plants and has profound impacts on the genetic structure of populations. Recent advances in understanding of the molecular basis of self-incompatibility have provided insights into the mechanisms of its loss in natural populations, especially in the tomato family, the Solanaceae. In the Solanaceae, the gene that controls self-incompatibility in the style codes for a ribonuclease that causes the degradation of RNA in pollen tubes bearing an allele at the S-locus that matches either of the two alleles held by the maternal plant. The pollen component of the S-locus has yet to be identified. Loss of self-incompatibility can be attributed to three types of causes: duplication of the S-locus, mutations that cause loss of S-RNase activity, and mutations that do not cause loss of S-RNase activity. Duplication of the S-locus has been well studied in radiation-induced mutants but may be a relatively rare cause of the breakdown of self-incompatibility in nature. Point mutations within the S-locus that disrupt the production of S-RNase have been documented in natural populations. There are also a number of mutants in which S-RNase production is unimpaired, yet self-incompatibility is disrupted. The identity and function of these mutations is not well understood. Careful work on a handful of model organisms will enable population biologists to better understand the breakdown of self-incompatibility in nature.

The dominance of alleles controlling self-incompatibility in Brassica pollen is regulated at the RNA level

Self-incompatibility (SI) in Brassica is controlled sporophytically by the multiallelic S-locus. The SI phenotype of pollen in an S-heterozygote is determined by the relationship between the two S-haplotypes it carries, and dominant/recessive relationships often are observed between the two S-haplotypes. The S-locus protein 11 (SP11, also known as the S-locus cysteine-rich protein) gene has been cloned from many pollen-dominant S-haplotypes (class I) and shown to encode the pollen S-determinant. However, SP11 from pollen-recessive S-haplotypes (class II) has never been identified by homology-based cloning strategies, and how the dominant/recessive interactions between the two classes occur was not known. We report here the identification and molecular characterization of SP11s from six class II S-haplotypes of B. rapa and B. oleracea. Phylogenetic analysis revealed that the class II SP11s form a distinct group separated from class I SP11s. The promoter sequences and expression patterns of SP11s also were different between the two classes. The mRNA of class II SP11, which was detected predominantly in the anther tapetum in homozygotes, was not detected in the heterozygotes of class I and class II S-haplotypes, suggesting that the dominant/recessive relationships of pollen are regulated at the mRNA level of SP11s.

Recognition specificity of self-incompatibility maintained after the divergence of Brassica oleracea and Brassica rapa

The determinants of recognition specificity of self-incompatibility in Brassica are SRK in the stigma and SP11/SCR in the pollen, respectively. In the pair of S haplotypes BrS46 (S46 in B. rapa) and BoS7 (S7 in B. oleracea), which have highly similar SRK alleles, the SP11 alleles were found to be similar, with 96.1% identity in the deduced amino acid sequence. Two other pairs of S haplotypes, BrS47 and BoS12, and BrS8 and BoS32, having highly similar SRK and SP11 alleles between the two species were also found. The haplotypes in each pair are considered to have been derived from a single S haplotype in the ancestral species. The allotetraploid produced by interspecific hybridization between homozygotes of BrS46 and BoS15 showed incompatibility with a BoS7 homozygote and compatibility with other B. oleracea S haplotypes in reciprocal crossings. This result indicates that BrS46 and BoS7 have maintained the same recognition specificity after the divergence of the two species and that amino acid substitutions found in such cases in both SRK alleles and SP11 alleles do not alter the recognition specificity. DNA blot analysis of SRK, SP11, SLG and other S-locus genes showed different DNA fragment sizes between the interspecific pairs of S haplotypes. A much lower level of sequence similarity was observed outside the genes of SRK and SP11 between BrS46 and BoS7. These results suggest that the DNA sequences of the regions intervening between the S-locus genes were diversified after or at the time of speciation. This is the first report demonstrating the presence of common S haplotypes in different plant species and presenting definite evidence of the trans-specific evolution of self-incompatibility genes.

Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach

Rates of molecular evolution vary widely between lineages, but quantification of how rates change has proven difficult. Recently proposed estimation procedures have mainly adopted highly parametric approaches that model rate evolution explicitly. In this study, a semiparametric smoothing method is developed using penalized likelihood. A saturated model in which every lineage has a separate rate is combined with a roughness penalty that discourages rates from varying too much across a phylogeny. A data-driven cross-validation criterion is then used to determine an optimal level of smoothing. This criterion is based on an estimate of the average prediction error associated with pruning lineages from the tree. The methods are applied to three data sets of six genes across a sample of land plants. Optimally smoothed estimates of absolute rates entailed 2- to 10-fold variation across lineages.

Intrahaplotype polymorphism at the Brassica S locus

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

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

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

On the origin of self-incompatibility haplotypes: transition through self-compatible intermediates

Self-incompatibility (SI) in flowering plants entails the inhibition of fertilization by pollen that express specificities in common with the pistil. In species of the Solanaceae, Rosaceae, and Scrophulariaceae, the inhibiting factor is an extracellular ribonuclease (S-RNase) secreted by stylar tissue. A distinct but as yet unknown gene (provisionally called pollen-S) appears to determine the specific S-RNase from which a pollen tube accepts inhibition. The S-RNase gene and pollen-S segregate with the classically defined S-locus. The origin of a new specificity appears to require, at minimum, mutations in both genes. We explore the conditions under which new specificities may arise from an intermediate state of loss of self-recognition. Our evolutionary analysis of mutations that affect either pistil or pollen specificity indicates that natural selection favors mutations in pollen-S that reduce the set of pistils from which the pollen accepts inhibition and disfavors mutations in the S-RNase gene that cause the nonreciprocal acceptance of pollen specificities. We describe the range of parameters (rate of receipt of self-pollen and relative viability of inbred offspring) that permits the generation of a succession of new specificities. This evolutionary pathway begins with the partial breakdown of SI upon the appearance of a mutation in pollen-S that frees pollen from inhibition by any S-RNase presently in the population and ends with the restoration of SI by a mutation in the S-RNase gene that enables pistils to reject the new pollen type.

Self-incompatibility in the genus Arabidopsis: characterization of the S locus in the outcrossing A. lyrata and its autogamous relative A. thaliana

As a starting point for a phylogenetic study of self-incompatibility (SI) in crucifers and to elucidate the genetic basis of transitions between outcrossing and self-fertilizing mating systems in this family, we investigated the SI system of Arabidopsis lyrata. A. lyrata is an outcrossing close relative of the self-fertile A. thaliana and is thought to have diverged from A. thaliana approximately 5 million years ago and from Brassica spp 15 to 20 million years ago. Analysis of two S (sterility) locus haplotypes demonstrates that the A. lyrata S locus contains tightly linked orthologs of the S locus receptor kinase (SRK) gene and the S locus cysteine-rich protein (SCR) gene, which are the determinants of SI specificity in stigma and pollen, respectively, but lacks an S locus glycoprotein gene. As described previously in Brassica, the S haplotypes of A. lyrata differ by the rearranged order of their genes and by their variable physical sizes. Comparative mapping of the A. lyrata and Brassica S loci indicates that the S locus of crucifers is a dynamic locus that has undergone several duplication events since the Arabidopsis--Brassica split and was translocated as a unit between two distant chromosomal locations during diversification of the two taxa. Furthermore, comparative analysis of the S locus region of A. lyrata and its homeolog in self-fertile A. thaliana identified orthologs of the SRK and SCR genes and demonstrated that self-compatibility in this species is associated with inactivation of SI specificity genes.

Self-incompatibility in the genus Arabidopsis: characterization of the S locus in the outcrossing A. lyrata and its autogamous relative A. thaliana

As a starting point for a phylogenetic study of self-incompatibility (SI) in crucifers and to elucidate the genetic basis of transitions between outcrossing and self-fertilizing mating systems in this family, we investigated the SI system of Arabidopsis lyrata. A. lyrata is an outcrossing close relative of the self-fertile A. thaliana and is thought to have diverged from A. thaliana approximately 5 million years ago and from Brassica spp 15 to 20 million years ago. Analysis of two S (sterility) locus haplotypes demonstrates that the A. lyrata S locus contains tightly linked orthologs of the S locus receptor kinase (SRK) gene and the S locus cysteine-rich protein (SCR) gene, which are the determinants of SI specificity in stigma and pollen, respectively, but lacks an S locus glycoprotein gene. As described previously in Brassica, the S haplotypes of A. lyrata differ by the rearranged order of their genes and by their variable physical sizes. Comparative mapping of the A. lyrata and Brassica S loci indicates that the S locus of crucifers is a dynamic locus that has undergone several duplication events since the Arabidopsis--Brassica split and was translocated as a unit between two distant chromosomal locations during diversification of the two taxa. Furthermore, comparative analysis of the S locus region of A. lyrata and its homeolog in self-fertile A. thaliana identified orthologs of the SRK and SCR genes and demonstrated that self-compatibility in this species is associated with inactivation of SI specificity genes.

Mutational origin of new mating type specificities in flowering plants

Many hermaphroditic plants avoid self-fertilization by rejecting pollen that express genetically-determined specificities in common with the pistil. Self-incompatibility systems typically show extremely high genetic diversity, some maintaining hundreds of specificities. This article addresses the genetic and evolutionary mechanisms through which new mating specificities arise. Recent investigations of the genetic and physiological basis of self-incompatibility are reviewed. Two evolutionary pathways are considered: one which requires full expression of self-incompatibility in all intermediates and one in which new mating specificities arise through episodes of partial breakdown and restoration of self-incompatibility.

Estimating divergence times in the presence of an overdispersed molecular clock

Molecular loci that fail relative-rate tests are said to be "overdispersed." Traditional molecular-clock approaches to estimating divergence times do not take this into account. In this study, a method was developed to estimate divergence times using loci that may be overdispersed. The approach was to replace the traditional Poisson process assumption with a more general stationary process assumption. A probability model was developed, and an accompanying computer program was written to find maximum-likelihood estimates of divergence times under both the Poisson process and the stationary process assumptions. In simulation, it was shown that confidence intervals under the traditional Poisson assumptions often vastly underestimate the true confidence limits for overdispersed loci. Both models were applied to two data sets: one from land plants, the other from the higher metazoans. In both cases, the traditional Poisson process model could be rejected with high confidence. Maximum-likelihood analysis of the metazoan data set under the more general stationary process suggested that their radiation occurred well over a billion years ago, but confidence intervals were extremely wide. It was also shown that a model consistent with a Cambrian (or nearly Cambrian) origination of the animal phyla, although significantly less likely than a much older divergence, fitted the data well. It is argued that without an a priori understanding of the variance in the time between substitutions, molecular data sets may be incapable of ever establishing the age of the metazoan radiation.

Evolutionary dynamics of self-incompatibility alleles in Brassica

Self-incompatibility in Brassica entails the rejection of pollen grains that express specificities held in common with the seed parent. In Brassica, pollen specificity is encoded at the multipartite S-locus, a complex region comprising many expressed genes. A number of species within the Brassicaceae express sporophytic self-incompatibility, under which individual pollen grains bear specificities determined by one or both S-haplotypes of the pollen parent. Classical genetic and nucleotide-level analyses of the S-locus have revealed a dichotomy in sequence and function among S-haplotypes; in particular, all class I haplotypes show dominance over all class II haplotypes in determination of pollen specificity. Analysis of an evolutionary model that explicitly incorporates features of the Brassica system, including the class dichotomy, indicates that class II haplotypes may invade populations at lower rates and decline to extinction at higher rates than class I haplotypes. This analysis suggests convergence to an evolutionarily persistent state characterized by the maintenance in high frequency of a single class II haplotype together with many class I haplotypes, each in low frequency. This expectation appears to be consistent with empirical observations of high frequencies of relatively few distinct recessive haplotypes.

Cellular and molecular mechanisms of sexual incompatibility in plants and fungi

Plants and fungi show an astonishing diversity of mechanisms to promote outbreeding, the most widespread of which is sexual incompatibility. Sexual incompatibility involves molecular recognition between mating partners. In fungi and algae, highly polymorphic mating-type loci mediate mating through complementary interactions between molecules encoded or regulated by different mating-type haplotypes, whereas in flowering plants polymorphic self-incompatibility loci regulate mate recognition through oppositional interactions between molecules encoded by the same self-incompatibility haplotypes. This subtle mechanistic difference is a consequence of the different life cycles of fungi, algae, and flowering plants. Recent molecular and biochemical studies have provided fascinating insights into the mechanisms of mate recognition and are beginning to shed light on evolution and population genetics of these extraordinarily polymorphic genetic systems of incompatibility.

The signature of balancing selection: fungal mating compatibility gene evolution

A key problem in evolutionary biology has been distinguishing the contributions of current and historical processes to the maintenance of genetic variation. Because alleles at self-recognition genes are under balancing selection, they exhibit extended residence times in populations and thus may provide unique insight into population demographic history. However, evidence for balancing selection and extended residence times has almost exclusively depended on identification of transspecific polymorphisms; polymorphisms retained in populations through speciation events. We present a broadly applicable approach for detecting balancing selection and apply it to the b1 mating type gene in the mushroom fungus Coprinus cinereus. The comparison of neutral molecular variation within and between allelic classes was used to directly estimate the strength of balancing selection. Different allelic classes are defined as encoding different mating compatibility types and are thus potentially subject to balancing selection. Variation within an allelic class, where all alleles have the same mating compatibility type, provided an internal standard of neutral evolution. Mating compatibility in this organism is determined by the complex A mating type locus, and b1 is one of several redundantly functioning genes. Consequently, we conducted numerical simulations of a model with two subloci and varying levels of recombination to show that balancing selection should operate at each sublocus. Empirical data show that strong balancing selection has indeed occurred at the b1 locus. The widespread geographic distribution of identical b1 alleles suggests that their association with differing A mating types is the result of recent recombination events.