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

Evolutionary and structural aspects of Solanaceae RNases T2

Plant RNases T2 are involved in several physiological and developmental processes, including inorganic phosphate starvation, senescence, wounding, defense against pathogens, and the self-incompatibility system. Solanaceae RNases form three main clades, one composed exclusively of S-RNases and two that include S-like RNases. We identified several positively selected amino acids located in highly flexible regions of these molecules, mainly close to the B1 and B2 substrate-binding sites in S-like RNases and the hypervariable regions of S-RNases. These differences between S- and S-like RNases in the flexibility of amino acids in substrate-binding regions are essential to understand the RNA-binding process. For example, in the S-like RNase NT, two positively selected amino acid residues (Tyr156 and Asn134) are located at the most flexible sites on the molecular surface. RNase NT is induced in response to tobacco mosaic virus infection; these sites may thus be regions of interaction with pathogen proteins or viral RNA. Differential selective pressures acting on plant ribonucleases have increased amino acid variability and, consequently, structural differences within and among S-like RNases and S-RNases that seem to be essential for these proteins play different functions.

Phylogeography and Population Genetics Analyses Reveal Evolutionary History of the Desert Resource Plant Lycium ruthenicum (Solanaceae)

Climactic oscillations during the Quaternary played a significant role in the formation of genetic diversity and historical demography of numerous plant species in northwestern China. In this study, we used 11 simple sequence repeats derived from expressed sequence tag (EST-SSR), two chloroplast DNA (cpDNA) fragments, and ecological niche modeling (ENM) to investigate the population structure and the phylogeographic history of Lycium ruthenicum, a plant species adapted to the climate in northwestern China. We identified 20 chloroplast haplotypes of which two were dominant and widely distributed in almost all populations. The species has high haplotype diversity and low nucleotide diversity based on the cpDNA data. The EST-SSR results showed a high percentage of total genetic variation within populations. Both the cpDNA and EST-SSR results indicated no significant differentiation among populations. By combining the evidence from ENM and demographic analysis, we confirmed that both the last interglacial (LIG) and late-glacial maximum (LGM) climatic fluctuations, aridification might have substantially narrowed the distribution range of this desert species, the southern parts of the Junggar Basin, the Tarim Basin, and the eastern Pamir Plateau were the potential glacial refugia for L. ruthenicum during the late middle Pleistocene to late Pleistocene Period. During the early Holocene, the warm, and humid climate promoted its demographic expansion in northwestern China. This work may provide new insights into the mechanism of formation of plant diversity in this arid region.

Molecular signatures of long-distance oceanic dispersal and the colonization of Pacific islands in Lycium carolinianum

PREMISE

Long-distance dispersal has been important in explaining the present distributions of many plant species. Despite being infrequent, such dispersal events have considerable evolutionary consequences, because bottlenecks during colonization can result in reduced genetic diversity. We examined the phylogeographic history of Lycium carolinianum, a widespread taxon that ranges from southeastern North America to several Pacific islands, with intraspecific diversity in sexual and mating systems.

METHODS

We used Bayesian, likelihood, and coalescent approaches with nuclear and plastid sequence data and genome-wide single nucleotide polymorphisms to reconstruct the dispersal history of this species. We also compared patterns of genetic variation in mainland and island populations using single nucleotide polymorphisms and allelic diversity at the S-RNase mating system gene.

RESULTS

Lycium carolinianum is monophyletic and dispersed once from the North American mainland, colonizing the Pacific islands ca. 40,100 years ago. This dispersal was accompanied by a loss of genetic diversity in SNPs and the S-RNase locus due to a colonization bottleneck and the loss of self-incompatibility. Additionally, we documented at least two independent transitions to gynodioecy: once following the colonization of the Hawaiian Islands and loss of self-incompatibility, and a second time associated with polyploidy in the Yucatán region of Mexico.

CONCLUSIONS

Long-distance dispersal via fleshy, bird dispersed fruits best explains the unusually widespread distribution of L. carolinianum. The collapse of diversity at the S-RNase locus in island populations suggests that self-fertilization may have facilitated the subsequent colonization of Pacific islands following a single dispersal from mainland North America.

Predicting Specificities Under the Non-self Gametophytic Self-Incompatibility Recognition Model

Non-self gametophytic self-incompatibility (GSI) recognition system is characterized by the presence of multiple F-box genes tandemly located in the S-locus, that regulate pollen specificity. This reproductive barrier is present in Solanaceae, Plantaginacea and Maleae (Rosaceae), but only in Petunia functional assays have been performed to get insight on how this recognition mechanism works. In this system, each of the encoded S-pollen proteins (called SLFs in Solanaceae and Plantaginaceae /SFBBs in Maleae) recognizes and interacts with a sub-set of non-self S-pistil proteins, called S-RNases, mediating their ubiquitination and degradation. In Petunia there are 17 SLF genes per S-haplotype, making impossible to determine experimentally each SLF specificity. Moreover, domain -swapping experiments are unlikely to be performed in large scale to determine S-pollen and S-pistil specificities. Phylogenetic analyses of the Petunia SLFs and those from two Solanum genomes, suggest that diversification of SLFs predate the two genera separation. Here we first identify putative SLF genes from nine Solanum and 10 Nicotiana genomes to determine how many gene lineages are present in the three genera, and the rate of origin of new SLF gene lineages. The use of multiple genomes per genera precludes the effect of incompleteness of the genome at the S-locus. The similar number of gene lineages in the three genera implies a comparable effective population size for these species, and number of specificities. The rate of origin of new specificities is one per 10 million years. Moreover, here we determine the amino acids positions under positive selection, those involved in SLF specificity recognition, using 10 Petunia S-haplotypes with more than 11 SLF genes. These 16 amino acid positions account for the differences of self-incompatible (SI) behavior described in the literature. When SLF and S-RNase proteins are divided according to the SI behavior, and the positively selected amino acids classified according to hydrophobicity, charge, polarity and size, we identified fixed differences between SI groups. According to the in silico 3D structure of the two proteins these amino acid positions interact. Therefore, this methodology can be used to infer SLF/S-RNase specificity recognition.

Colonization, Baker's law, and the evolution of gynodioecy in Hawaii: implications from a study of Lycium carolinianum

PREMISE

As Baker's law suggests, the successful colonization of oceanic islands is often associated with uniparental reproduction (self-fertility), but the high incidence of dimorphism (dioecy, gynodioecy) on islands complicates this idea. Lycium carolinianum is widespread, occurring on the North American mainland and the Hawaiian Islands. We examined Baker's ideas for mainland and island populations of L. carolinianum and examined inbreeding depression as a possible contributor to the evolution of gynodioecy on Maui.

METHODS

Controlled crosses were conducted in two mainland populations and two populations in Hawaii. Treatments included self and cross pollination, unmanipulated controls, and autogamy/agamospermy. Alleles from the self-incompatibility S-RNase gene were isolated and compared between mainland and island populations. Given self-compatibility in Hawaii, we germinated seeds from self- and cross- treatments and estimated inbreeding depression using seven traits and a measure of cumulative fitness.

RESULTS

Mainland populations of Lycium carolinianum are predominately self-incompatible with some polymorphism for self-fertility, whereas Hawaiian populations are self-compatible. Concordantly, S-RNase allelic diversity is reduced in Hawaii compared to the mainland. Hawaiian populations also exhibit significant inbreeding depression.

CONCLUSIONS

Self-compatibility in Hawaii and individual variation in self-fertility in mainland populations suggests that a colonization filter promoting uniparental reproduction may be acting in this system. Comparison of S-RNase variation suggests a collapse of allelic diversity and heterozygosity at the S-RNase locus in Hawaii, which likely contributed to mate limitation upon arrival to the Pacific. Inbreeding depression coupled with autonomous self-fertilization may have led to the evolution of gynodioecy on Maui.

Correlated polymorphism in cytotype and sexual system within a monophyletic species, Lycium californicum

BACKGROUND AND AIMS

Polyploidy has important effects on reproductive systems in plants and has been implicated in the evolution of dimorphic sexual systems. In particular, higher ploidy is associated with gender dimorphism across Lycium species (Solanaceae) and across populations within the species Lycium californicum. Previous research on the association of cytotype and sexual system within L. californicum sampled a limited portion of the species range, and did not investigate evolutionary transitions between sexual systems. Lycium californicum occurs in arid regions on offshore islands and mainland regions in the south-western United States and Mexico, motivating a more comprehensive analysis of intraspecific variation in sexual system and cytotype across the full range of this species.

METHODS

Sexual system (dimorphic vs. cosexual) was determined for 34 populations across the geographical range of L. californicum using field observations of pollen production, and was confirmed using morphological measurements and among-plant correlations of primary sexual traits. Ploidy was inferred using flow cytometry in 28 populations. DNA sequence data from four plastid and two nuclear regions were used to reconstruct relationships among populations and to map transitions in sexual system and ploidy.

KEY RESULTS

Lycium californicum is monophyletic, ancestrally diploid and cosexual, and the association of gender dimorphism and polyploidy appears to have two evolutionary origins in this species. Compared with cosexual populations, dimorphic populations had bimodal anther size distributions, negative correlations between male and female floral traits, and larger coefficients of variation for primary sexual traits. Flow cytometry confirmed tetraploidy in dimorphic populations, whereas cosexual populations were diploid.

CONCLUSIONS

Tetraploidy and gender dimorphism are perfectly correlated in L. californicum, and the distribution of tetraploid-dimorphic populations is restricted to populations in Arizona and the Baja California peninsula. The analysis suggests that tetraploidy and dimorphism likely established in Baja California and may have evolved multiple times.

Diversification and distinctive structural features of S-RNase alleles in the genus Solanum

The multigenic and multiallelic S-locus in plants is responsible for the gametophytic self-incompatibility system, which is important to prevent the detrimental effects of self-fertilization and inbreeding depression. Several studies have discussed the importance of punctual mutations, recombination, and natural selection in the generation of allelic diversity in the S-locus. However, there has been no wide-ranging study correlating the molecular evolution and structural aspects of the corresponding proteins in Solanum. Therefore, we evaluated the molecular evolution of one gene in this locus and generated a statistically well-supported phylogenetic tree, as well as evidence of positive selection, helping us to understand the diversification of S alleles in Solanum. The three-dimensional structures of some of the proteins corresponding to the major clusters of the phylogenetic tree were constructed and subsequently submitted to molecular dynamics to stabilize the folding and obtain the native structure. The positively selected amino acid residues were predominantly located in the hyper variable regions and on the surface of the protein, which appears to be fundamental for allele specificity. One of the positively selected residues was identified adjacent to a conserved strand that is crucial for enzymatic catalysis. Additionally, we have shown significant differences in the electrostatic potential among the predicted molecular surfaces in S-RNases. The structural results indicate that local changes in the three-dimensional structure are present in some regions of the molecule, although the general structure seems to be conserved. No previous study has described such structural variations in S-RNases.

Intraspecific variation in gender strategies in Lycium (Solanaceae): associations with ploidy and changes in floral form following the evolution of gender dimorphism

•

PREMISE OF THE STUDY

An association between polyploidy and gender dimorphism has been noted in several plant lineages. Whereas the majority of Lycium species are diploid and have hermaphroditic flowers in cosexual populations, gender dimorphism (gynodioecy, dioecy) has been shown to be uniformly associated with polyploidy in previous studies. Preliminary field observations suggested that some populations of Lycium carolinianum were dimorphic, providing a test of this association.•

METHODS

We assessed sexual systems and cytotype variation (to infer ploidy) across 17 populations of L. carolinianum. Comparison of flowers in cosexual and dimorphic populations were used to infer changes in reproductive morphology associated with the evolution of gynodioecy.•

KEY RESULTS

The majority of populations were cosexual in gender expression, but dimorphism was present in the Yucatán and in some populations in Hawaii. Populations varied in ploidy and were either diploid or tetraploid. Floral sexual dimorphism was present in all gynodioecious populations, though the magnitude differed and was cryptic in some cases. Our results are consistent with the hypothesis that following the evolution of gynodioecy, flowers on hermaphrodites increased in size.•

CONCLUSIONS

Dimorphic sexual systems have likely evolved convergently in L. carolinianum. In contrast to previous studies, dimorphism is not perfectly associated with polyploidy. Although our sample from the Yucatán was both tetraploid and dimorphic, all populations in Hawaii were diploid regardless of sexual system. Ongoing phylogeographic and mating system studies will contribute to our understanding of reproductive evolution in this widespread, polymorphic species.

The timetable for allopolyploidy in flowering plants

BACKGROUND

Our understanding of the processes and dynamics of allopolyploid speciation, the long-term consequences of ploidal change, and the genetic and chromosomal changes in new emerged allopolyploids has substantially increased during the past few decades. Yet we remain uncertain about the time since lineage divergence when two taxa are capable of spawning such entities. Indeed, the matter has seemed intractable. Knowledge of the window of opportunity for allopolyploid production is very important because it provides temporal insight into a key evolutionary process, and a temporal reference against which other modes of speciation may be measured.

SCOPE

This Viewpoint paper reviews and integrates published information on the crossability of herbaceous species and the fertility of their hybrids in relation to species' divergence times. Despite limitations in methodology and sampling, the estimated times to hybrid sterility are somewhat congruent across disparate lineages. Whereas the waiting time for hybrid sterility is roughly 4-5 million years, the waiting time for cross-incompatibility is roughly 8-10 million years, sometimes considerably more. Strict allopolyploids may be formed in the intervening time window. The progenitors of several allopolyploids diverged between 4 and 6 million years before allopolyploid synthesis, as expected. This is the first study to propose a general temporal framework for strict allopolyploidy. This Viewpoint paper hopefully will stimulate interest in studying the tempo of speciation and the tempo of reproductive isolation in general.

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.

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 and genetic characterization of novel S-RNases from a natural population of Nicotiana alata

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

The utility of nuclear conserved ortholog set II (COSII) genomic regions for species-level phylogenetic inference in Lycium (Solanaceae)

The identification of genomic regions with sufficient variation to elucidate fine-scale relationships among closely related species is a major goal of phylogenetic systematics. However, the accumulation of such multi-locus data sets brings its own challenges, given that gene trees do not necessarily represent the true species tree. Using genomic tools developed for Solanum (Solanaceae), we have evaluated the utility of nuclear conserved ortholog set II (COSII) regions for phylogenetic inference in tribe Lycieae (Solanaceae). Five COSII regions, with intronic contents ranging from 68% to 91%, were sequenced in 10 species. Their phylogenetic utility was assessed and compared with data from more commonly used nuclear (GBSSI, nrITS) and cpDNA spacer data. We compared the effectiveness of a traditional total evidence concatenation approach versus the recently developed Bayesian estimation of species trees (BEST) method to infer species trees given multiple independent gene trees. All of the sampled COSII regions had high numbers of parsimony-informative (PI) characters, and two of the COSII regions had more PI characters than the GBSSI, ITS, and cpDNA spacer data sets combined. COSII markers are a promising new tool for phylogenetic inference in Solanaceae, and should be explored in related groups. Both the concatenation and BEST approaches yielded similar topologies; however, when multiple individuals with polyphyletic alleles were included, BEST was clearly the more robust approach for inferring species trees in the presence of gene tree incongruence.

Positive Darwinian selection at single amino acid sites conferring plant virus resistance

Explicit evaluation of the accuracy and power of maximum likelihood and Bayesian methods for detecting site-specific positive Darwinian selection presents a challenge because selective consequences of single amino acid changes are generally unknown. We exploited extensive molecular and functional characterization of amino acid substitutions in the plant gene eIF4E to evaluate the performance of these methods in detecting site-specific positive selection. We documented for the first time a molecular signature of positive selection within a recessive resistance gene in plants. We then used two statistical platforms, Phylogenetic Analysis Using Maximum Likelihood and Hypothesis Testing Using Phylogenies (HyPhy), to look for site-specific positive selection. Their relative power and accuracy are assessed by comparing the sites they identify as being positively selected with those of resistance-determining amino acids. Our results indicate that although both methods are surprisingly accurate in their identification of resistance sites, HyPhy appears to more accurately identify biologically significant amino acids using our data set.

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.

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.