Patterns of variation within self-incompatibility loci

Relaxed Selection and the Rapid Evolution of Reproductive Genes

Evolutionary genomic studies find that reproductive protein genes, those directly involved in reproductive processes, diversify more rapidly than most other gene categories. Strong postcopulatory sexual selection acting within species is the predominant hypothesis proposed to account for the observed pattern. Recently, relaxed selection due to sex-specific gene expression has also been put forward to explain the relatively rapid diversification. We contend that relaxed selection due to sex-limited gene expression is the correct null model for tests of molecular evolution of reproductive genes and argue that it may play a more significant role in the evolutionary diversification of reproductive genes than previously recognized. We advocate for a re-evaluation of adaptive explanations for the rapid diversification of reproductive genes.

Phylogenomics of the genus Populus reveals extensive interspecific gene flow and balancing selection

Phylogenetic analysis is complicated by interspecific gene flow and the presence of shared ancestral polymorphisms, particularly those maintained by balancing selection. In this study, we aimed to examine the prevalence of these factors during the diversification of Populus, a model tree genus in the Northern Hemisphere. We constructed phylogenetic trees of 29 Populus taxa using 80 individuals based on re-sequenced genomes. Our species tree analyses recovered four main clades in the genus based on consensus nuclear phylogenies, but in conflict with the plastome phylogeny. A few interspecific relationships remained unresolved within the multiple-species clade because of inconsistent gene trees. Our results indicated that gene flow has been widespread within each clade and also occurred among the four clades during their early divergence. We identified 45 candidate genes with ancient polymorphisms maintained by balancing selection. These genes were mainly associated with mating compatibility, growth and stress resistance. Both gene flow and selection-mediated ancient polymorphisms are prevalent in the genus Populus. These are potentially important contributors to adaptive variation. Our results provide a framework for the diversification of model tree genus that will facilitate future comparative studies.

CaDHN4, a Salt and Cold Stress-Responsive Dehydrin Gene from Pepper Decreases Abscisic Acid Sensitivity in Arabidopsis

Dehydrins play an important role in improving plant resistance to abiotic stresses. In this study, we isolated a dehydrin gene from pepper (Capsicum annuum L.) leaves, designated as CaDHN4. Sub-cellular localization of CaDHN4 was to be found in the nucleus and membrane. To authenticate the function of CaDHN4 in cold- and salt-stress responses and abscisic acid (ABA) sensitivity, we reduced the CaDHN4 expression using virus-induced gene silencing (VIGS), and overexpressed the CaDHN4 in Arabidopsis. We found that silencing of CaDHN4 reduced the growth of pepper seedlings and CaDHN4-silenced plants exhibited more serious wilting, higher electrolyte leakage, and more accumulation of ROS in the leaves compared to pTRV2:00 plants after cold stress, and lower chlorophyll contents and higher electrolyte leakage compared to pTRV2:00 plants under salt stress. However, CaDHN4-overexpressing Arabidopsis plants had higher seed germination rates and post-germination primary root growth, compared to WT plants under salt stress. In response to cold and salt stresses, the CaDHN4-overexpressed Arabidopsis exhibited lower MDA content, and lower relative electrolyte leakage compared to the WT plants. Under ABA treatments, the fresh weight and germination rates of transgenic plants were higher than WT plants. The transgenic Arabidopsis expressing a CaDHN4 promoter displayed a more intense GUS staining than the normal growth conditions under treatment with hormones including ABA, methyl jasmonate (MeJA), and salicylic acid (SA). Our results suggest that CaDHN4 can protect against cold and salt stresses and decrease ABA sensitivity in Arabidopsis.

Organization and evolution of mating-type genes in three Stagonosporopsis species causing gummy stem blight of cucurbits and leaf spot and dry rot of papaya

Population divergence and speciation of closely related lineages can result from reproductive differences leading to genetic isolation. An increasing number of fungal diseases of plants and animals have been determined to be caused by morphologically indistinguishable species that are genetically distinct, thereby representing cryptic species. We were interested in identifying if mating systems among three Stagonosporopsis species (S. citrulli, S. cucurbitacearum, and S. caricae) causing gummy stem blight (GSB) of cucurbits or leaf spot and dry rot of papaya differed, possibly underlying species divergence. Additionally, we were interested in identifying evolutionary pressures acting on the genes controlling mating in these fungi. The mating-type loci (MAT1) of three isolates from each of the three species were identified in draft genome sequences. For the three species, MAT1 was structurally identical and contained both mating-type genes necessary for sexual reproduction, which suggests that all three species are homothallic. However, both MAT1-1-1 and MAT1-2-1 were divergent among species showing rapid evolution with a much greater number of amino acid-changing substitutions detected for the reproductive genes compared with genes flanking MAT1. Positive selection was detected in MAT1-2-1, especially in the highly conserved high mobility group (MATA_HMG-box) domain. Thus, the mating-type genes are rapidly evolving in GSB fungi, but a difference in mating systems among the three species does not underlie their divergence.

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.

Evidence for selection on a chordate histocompatibility locus

Allorecognition is the ability of an organism to differentiate self or close relatives from unrelated individuals. The best known applications of allorecognition are the prevention of inbreeding in hermaphroditic species (e.g., the self-incompatibility [SI] systems in plants), the vertebrate immune response to foreign antigens mediated by MHC loci, and somatic fusion, where two genetically independent individuals physically join to become a chimera. In the few model systems where the loci governing allorecognition outcomes have been identified, the corresponding proteins have exhibited exceptional polymorphism. But information about the evolution of this polymorphism outside MHC is limited. We address this subject in the ascidian Botryllus schlosseri, where allorecognition outcomes are determined by a single locus, called FuHC (Fusion/HistoCompatibility). Molecular variation in FuHC is distributed almost entirely within populations, with very little evidence for differentiation among different populations. Mutation plays a larger role than recombination in the creation of FuHC polymorphism. A selection statistic, neutrality tests, and distribution of variation within and among different populations all provide evidence for selection acting on FuHC, but are not in agreement as to whether the selection is balancing or directional.

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.

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.

Temporally variable selection on proteolysis-related reproductive tract proteins in Drosophila

In order to gain further insight into the processes underlying rapid reproductive protein evolution, we have conducted a population genetic survey of 44 reproductive tract-expressed proteases, protease inhibitors, and targets of proteolysis in Drosophila melanogaster and Drosophila simulans. Our findings suggest that positive selection on this group of genes is temporally heterogeneous, with different patterns of selection inferred using tests sensitive at different time scales. Such variation in the strength and targets of selection through time may be expected under models of sexual conflict and/or host-pathogen interaction. Moreover, available functional information concerning the genes that show evidence of selection suggests that both sexual selection and immune processes have been important in the evolutionary history of this group of molecules.

Differential strengths of selection on S-RNases from Physalis and Solanum (Solanaceae)

Background

The S-RNases of the Solanaceae are highly polymorphic self-incompatibility (S-) alleles subject to strong balancing selection. Relatively recent diversification of S-alleles has occurred in the genus Physalis following a historical restriction of S-allele diversity. In contrast, the genus Solanum did not undergo a restriction of S-locus diversity and its S-alleles are generally much older. Because recovery from reduced S-locus diversity should involve increased selection, we employ a statistical framework to ask whether S-locus selection intensities are higher in Physalis than Solanum. Because different S-RNase lineages diversify in Physalis and Solanum, we also ask whether different sites are under selection in different lineages.

Results

Maximum-likelihood and Bayesian coalescent methods found higher intensities of selection and more sites under significant positive selection in the 48 Physalis S-RNase alleles than the 49 from Solanum. Highest posterior densities of dN/dS (ω) estimates show that the strength of selection is greater for Physalis at 36 codons. A nested maximum likelihood method was more conservative, but still found 16 sites with greater selection in Physalis. Neither method found any codons under significantly greater selection in Solanum. A random effects likelihood method that examines data from both taxa jointly confirmed higher selection intensities in Physalis, but did not find different proportions of sites under selection in the two datasets. The greatest differences in strengths of selection were found in the most variable regions of the S-RNases, as expected if these regions encode self-recognition specificities. Clade-specific likelihood models indicated some codons were under greater selection in background Solanum lineages than in specific lineages of Physalis implying that selection on sites may differ among lineages.

Conclusions

Likelihood and Bayesian methods provide a statistical approach to testing differential selection across populations or species. These tests appear robust to the levels of polymorphism found in diverse S-allele collections subject to strong balancing selection. As predicted, the intensity of selection at the S-locus was higher in the taxon with more recent S-locus diversification. This is the first confirmation by statistical test of differing selection intensities among self-incompatibility alleles from different populations or species.

The molecular evolution of animal reproductive tract proteins: what have we learned from mating-system comparisons?

Postcopulatory sexual selection is thought to drive the rapid evolution of reproductive tract genes in many animals. Recently, a number of studies have sought to test this hypothesis by examining the effects of mating system variation on the evolutionary rates of reproductive tract genes. Perhaps surprisingly, there is relatively little evidence that reproductive proteins evolve more rapidly in species subject to strong postcopulatory sexual selection. This emerging trend may suggest that other processes, such as host-pathogen interactions, are the main engines of rapid reproductive gene evolution. I suggest that such a conclusion is as yet unwarranted; instead, I propose that more rigorous analytical techniques, as well as multigene and population-based approaches, are required for a full understanding of the consequences of mating system variation for the evolution of reproductive tract genes.

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.

Molecular population genetics of the SRK and SCR self-incompatibility genes in the wild plant species Brassica cretica (Brassicaceae)

Self-incompatibility (SI) in plants is a classic example of a trait evolving under strong frequency-dependent selection. As a consequence, population genetic theory predicts that the S locus, which controls SI, should maintain numerous alleles, display a high level of nucleotide diversity, and, in structured populations, show a lower level of among-population differentiation compared to neutral loci. Population-level investigations of DNA sequence variation at the S locus have recently been carried out in the genus Arabidopsis, largely confirming results from theoretical models of S-locus evolutionary dynamics, but no comparable studies have been done in wild Brassica species. In this study, we sequenced parts of the S-locus genes SRK and SCR, two tightly linked genes that are directly involved in the determination of SI specificity in samples from four natural populations of the wild species Brassica cretica. The amount and distribution of nucleotide diversity, as well as the frequency spectrum of putative functional haplotypes, observed at the S locus in B. cretica fit very well with expectations from theoretical models, providing strong evidence for frequency-dependent selection acting on the S locus in a wild Brassica species.

RNase-based self-incompatibility: puzzled by pollen S

Many plants have a genetically determined self-incompatibility system in which the rejection of self pollen grains is controlled by alleles of an S locus. A common feature of these S loci is that separate pollen- and style-expressed genes (pollen S and style S, respectively) determine S allele identity. The long-held view has been that pollen S and style S must be a coevolving gene pair in order for allelic recognition to be maintained as new S alleles arise. In at least three plant families, the Solanaceae, Rosaceae, and Plantaginaceae, the style S gene has long been known to encode an extracellular ribonuclease called the S-RNase. Pollen S in these families has more recently been identified and encodes an F-box protein known as either SLF or SFB. In this perspective, we describe the puzzling evolutionary relationship that exists between the SLF/SFB and S-RNase genes and show that in most cases cognate pairs of genes are not coevolving in the expected manner. Because some pollen S genes appear to have arisen much more recently than their style S cognates, we conclude that either some pollen S genes have been falsely identified or that there is a major problem with our understanding of how the S locus evolves.

Unequal segregation of SRK alleles at the S locus in Brassica cretica

In the Brassicaceae plant family, which includes the Arabidopsis and Brassica genera, self-incompatibility (SI) is controlled by genes at the S locus. Using experimental crosses, we studied the pattern of inheritance of S-locus alleles in the wild species Brassica cretica. Four full-sib families were established and unequal segregation of alleles at the SRK SI gene was found in one family. The segregation distortion acted in favour of a recessive (class II) allele and was best explained by some form of gametic-level selection. Our findings are discussed in the light of theoretical predictions of differential accumulation of deleterious mutations among S-locus alleles.

Reproductive protein evolution within and between species: maintenance of divergent ZP3 alleles in Peromyscus

In a variety of animal taxa, proteins involved in reproduction evolve more rapidly than nonreproductive proteins. Most studies of reproductive protein evolution, however, focus on divergence between species, and little is known about differentiation among populations within a species. Here we investigate the molecular population genetics of the protein ZP3 within two Peromyscus species. ZP3 is an egg coat protein involved in primary binding of egg and sperm and is essential for fertilization. We find that amino acid polymorphism in the sperm-combining region of ZP3 is high relative to silent polymorphism in both species of Peromyscus. In addition, while there is geographical structure at a mitochondrial gene (Cytb), a nuclear gene (Lcat) and eight microsatellite loci, we find no evidence for geographical structure at Zp3 in Peromyscus truei. These patterns are consistent with the maintenance of ZP3 alleles by balancing selection, possibly due to sexual conflict or pathogen resistance. However, we do not find evidence that reinforcement promotes ZP3 diversification; allelic variation in P. truei is similar among populations, including populations allopatric and sympatric with sibling species. In fact, most alleles are present in all populations sampled across P. truei's range. While additional data are needed to identify the precise evolutionary forces responsible for sequence variation in ZP3, our results suggest that in Peromyscus, selection to maintain divergent alleles within species contributes to the pattern of rapid amino acid substitution observed among species.

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.

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.

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

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

Expression of 10 S-class SLF-like genes in Nicotiana alata pollen and its implications for understanding the pollen factor of the S locus

The S locus of Nicotiana alata encodes a polymorphic series of ribonucleases (S-RNases) that determine the self-incompatibility (SI) phenotype of the style. The pollen product of the S locus (pollen S) in N. alata is unknown, but in species from the related genus Petunia and in self-incompatible members of the Plantaginaceae and Rosaceae, this function has been assigned to an F-box protein known as SLF or SFB. Here we describe the identification of 10 genes (designated DD1-10) encoding SLF-related proteins that are expressed in N. alata pollen. Because our approach to cloning the DD genes was based on sequences of SLFs from other species, we presume that one of the DD genes encodes the N. alata SLF ortholog. Seven of the DD genes were exclusively expressed in pollen and a low level of sequence variation was found in alleles of each DD gene. Mapping studies confirmed that all 10 DD genes were linked to the S locus and that at least three were located in the same chromosomal segment as pollen S. Finally, the different topologies of the phylogenetic trees produced using available SLF-related sequences and those produced using S-RNase sequences suggests that pollen S and the S-RNase have different evolutionary histories.

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.

Evolution under strong balancing selection: how many codons determine specificity at the female self-incompatibility gene SRK in Brassicaceae?

Background

Molecular lock-and-key systems are common among reproductive proteins, yet their evolution remains a major puzzle in evolutionary biology. In the Brassicaceae, the genes encoding self-incompatibility have been identified, but technical challenges currently prevent detailed analyses of the molecular interaction between the male and female components. In the present study, we investigate sequence polymorphism in the female specificity determinant SRK of Arabidopsis halleri from throughout Europe. Using a comparative approach based on published SRK sequences in A. lyrata and Brassica, we track the signature of frequency-dependent selection acting on these genes at the codon level. Using simulations, we evaluate power and accuracy of our approach and estimate the proportion of codon sites involved in the molecular interaction.

The transition to self-compatibility in Arabidopsis thaliana and evolution within S-haplotypes over 10 Myr

A recent investigation found evidence that the transition of Arabidopsis thaliana from ancestral self-incompatibility (SI) to full self-compatibility occurred very recently and suggested that this occurred through a selective fixation of a nonfunctional allele (PsiSCR1) at the SCR gene, which determines pollen specificity in the incompatibility response. The main evidence is the lack of polymorphism at the SCR locus in A. thaliana. However, the nearby SRK gene, which determines stigma specificity in self-incompatible Brassicaceae species, has extremely high sequence diversity, with 3 very divergent SRK haplotypes, 2 of them present in multiple strains. Such high diversity is extremely unusual in this species, and it suggests the possibility that multiple, different SRK haplotypes may have been preserved from A. thaliana's self-incompatible ancestor. To study the evolution of S-haplotypes in the A. thaliana lineage, we searched the 2 most closely related Arabidopsis species Arabidopsis lyrata and Arabidopsis halleri, in which most populations have retained SI, and found SRK sequences corresponding to all 3 A. thaliana haplogroup sequences. Our molecular evolutionary analyses of these 3 S-haplotypes provide an independent estimate of the timing of the breakdown of SI and again exclude an ancient transition to selfing in A. thaliana. Comparing sequences of each of the 3 haplogroups between species, we find that 2 of the 3 SRK sequences (haplogroups A and B) are similar throughout their length, suggesting that little or no recombination with other SRK alleles has occurred since these species diverged. The diversity difference between the SCR and SRK loci in A. thaliana, however, suggests crossing-over, either within SRK or between the SCR and SRK loci. If the loss of SI involved fixation of the PsiSCR1 sequence, the exchange must have occurred during its fixation. Divergence between the species is much lower at the S-locus, compared with reference loci, and we discuss two contributory possibilities. Introgression may have occurred between A. lyrata and A. halleri and between their ancestral lineage and A. thaliana, at least for some period after their split. In addition, the coalescence times of sequences of individual S-haplogroups are expected to be less than those of alleles at non-S-loci.

Trans-specificity at loci near the self-incompatibility loci in Arabidopsis

We compared allele sequences of two loci near the Arabidopsis lyrata self-incompatibility (S) loci with sequences of A. thaliana orthologs and found high numbers of shared polymorphisms, even excluding singletons and sites likely to be highly mutable. This suggests maintenance of entire S-haplotypes for long evolutionary times and extreme recombination suppression in the region.

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.

Selectionism and neutralism in molecular evolution

Charles Darwin proposed that evolution occurs primarily by natural selection, but this view has been controversial from the beginning. Two of the major opposing views have been mutationism and neutralism. Early molecular studies suggested that most amino acid substitutions in proteins are neutral or nearly neutral and the functional change of proteins occurs by a few key amino acid substitutions. This suggestion generated an intense controversy over selectionism and neutralism. This controversy is partially caused by Kimura's definition of neutrality, which was too strict (|2Ns|< or =1). If we define neutral mutations as the mutations that do not change the function of gene products appreciably, many controversies disappear because slightly deleterious and slightly advantageous mutations are engulfed by neutral mutations. The ratio of the rate of nonsynonymous nucleotide substitution to that of synonymous substitution is a useful quantity to study positive Darwinian selection operating at highly variable genetic loci, but it does not necessarily detect adaptively important codons. Previously, multigene families were thought to evolve following the model of concerted evolution, but new evidence indicates that most of them evolve by a birth-and-death process of duplicate genes. It is now clear that most phenotypic characters or genetic systems such as the adaptive immune system in vertebrates are controlled by the interaction of a number of multigene families, which are often evolutionarily related and are subject to birth-and-death evolution. Therefore, it is important to study the mechanisms of gene family interaction for understanding phenotypic evolution. Because gene duplication occurs more or less at random, phenotypic evolution contains some fortuitous elements, though the environmental factors also play an important role. The randomness of phenotypic evolution is qualitatively different from allele frequency changes by random genetic drift. However, there is some similarity between phenotypic and molecular evolution with respect to functional or environmental constraints and evolutionary rate. It appears that mutation (including gene duplication and other DNA changes) is the driving force of evolution at both the genic and the phenotypic levels.

Plant self-incompatibility systems: a molecular evolutionary perspective

Incompatibility recognition systems preventing self-fertilization have evolved several times in independent lineages of Angiosperm plants, and three main model systems are well characterized at the molecular level [the gametophytic self-incompatibility (SI) systems of Solanaceae, Rosaceae and Anthirrhinum, the very different system of poppy, and the system in Brassicaceae with sporophytic control of pollen SI reactions]. In two of these systems, the genes encoding both components of pollen-pistil recognition are now known, showing clearly that these two proteins are distinct, that is, SI is a lock-and-key mechanism. Here, we review recent findings in the three well-studied systems in the light of these results and analyse their implications for understanding polymorphism and coevolution of the two SI genes, in the context of a tightly linked genome region.

High-diversity genes in the Arabidopsis genome

High-diversity genes represent an important class of loci in organismal genomes. Since elevated levels of nucleotide variation are a key component of the molecular signature for balancing selection or local adaptation, high-diversity genes may represent loci whose alleles are selectively maintained as balanced polymorphisms. Comparison of 4300 random shotgun sequence fragments of the Arabidopsis thaliana Ler ecotype genome with the whole genomic sequence of the Col-0 ecotype identified 60 genes with putatively high levels of intraspecific variability. Eleven of these genes were sequenced in multiple A. thaliana accessions, 3 of which were found to display elevated levels of nucleotide polymorphism. These genes encode the myb-like transcription factor MYB103, a putative soluble starch synthase I, and a homeodomain-leucine zipper transcription factor. Analysis of these genes and 4-7 flanking genes in 14-20 A. thaliana ecotypes revealed that two of these loci show other characteristics of balanced polymorphisms, including broad peaks of nucleotide diversity spanning multiple linked genes and an excess of intermediate-frequency polymorphisms. Scanning genomes for high-diversity genomic regions may be useful in approaches to adaptive trait locus mapping for uncovering candidate balanced polymorphisms.

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.

Barley Dhn13 encodes a KS-type dehydrin with constitutive and stress responsive expression

Dehydrins (DHNs) compose a family of intrinsically unstructured proteins that have high water solubility and accumulate during late seed development, low temperature or water deficit conditions, and are thought to play a protective role in freezing and drought tolerance in plants. Twelve Dhn genes were previously described in the barley genome. Here, we report an additional member of this multigene family, Dhn13. The Dhn13 gene is located in chromosome 4 near marker MWG634 and encodes a 107-amino acid KS-type DHN. Semi-quantitative reverse transcriptase PCR data indicated that Dhn13 is constitutively expressed in seedling tissues and embryos of developing seeds. Microarray data were consistent with these results and showed a considerable increase of Dhn13 transcripts when plants were subjected to chilling and freezing temperatures. The highest transcript levels where observed in anthers. The presence of ABRE, MYC, DRE, and POLLEN1LELAT52 regulatory elements in the putative Dhn13 promoter region is in agreement with expression data.

Evolution under tight linkage to mating type

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

Molecular Population Genetics and the Search for Adaptive Evolution in Plants

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

A simple method for computing exact probabilities of mutation numbers

We describe a method for the recursive computation of exact probability distributions for the number of neutral mutations segregating in samples of arbitrary size and configuration. Construction of the recursions requires only characterization of evolutionary changes as a Markov process and determination of one-step transition matrices. We address the pattern of nucleotide diversity at a neutral marker locus linked to a determinant of mating type. Under a reformulation of parameters, the method also applies directly to metapopulation models with island migration among demes. Characterization of complete probability distributions facilitates parameter estimation and hypothesis testing by likelihood- as well as moment-based approaches.