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.     

Interhaplotypic heterogeneity and heterochromatic features may contribute to recombination suppression at the S locus in apple (Malusxdomestica)

Gametophytic self-incompatibility (GSI) is controlled by a complex S locus containing the pistil determinant S-RNase and pollen determinant SFB/SLF. Tight linkage of the pistil and pollen determinants is necessary to guarantee the self-incompatibility (SI) function. However, multiple probable pollen determinants of apple and Japanese pear, SFBBs (S locus F-box brothers), exist in each S haplotype, and how these multiple genes maintain the SI function remains unclear. It is shown here by high-resolution fluorescence in situ hybridization (FISH) that SFBB genes of the apple S ( 9 ) haplotype are physically linked to the S ( 9 ) -RNase gene, and the S locus is located in the subtelomeric region. FISH analyses also determined the relative order of SFBB genes and S-RNase in the S ( 9 ) haplotype, and showed that gene order differs between the S ( 9 ) and S ( 3 ) haplotypes. Furthermore, it is shown that the apple S locus is located in a knob-like large heterochromatin block where DNA is highly methylated. It is proposed that interhaplotypic heterogeneity and the heterochromatic nature of the S locus help to suppress recombination at the S locus in apple.     

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

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

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

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

Patterns of variation within self-incompatibility loci

Diverse self-incompatibility (SI) mechanisms permit flowering plants to inhibit fertilization by pollen that express specificities in common with the pistil. Characteristic of at least two model systems is greatly reduced recombination across large genomic tracts surrounding the S-locus, which regulates SI. In three angiosperm families, including the Solanaceae, the gene that controls the expression of gametophytic SI in the pistil encodes a ribonuclease (S-RNase). The gene that controls pollen SI expression is currently unknown, although several candidates have recently been proposed. Although each candidate shows a high level of polymorphism and complete allelic disequilibrium with the S-RNase gene, such properties may merely reflect tight linkage to the S-locus, irrespective of any functional role in SI. We analyzed the magnitude and nature of nucleotide variation, with the objective of distinguishing likely candidates for regulators of SI from other genes embedded in the S-locus region. We studied the S-RNase gene of the Solanaceae and 48A, a candidate for the pollen gene in this system, and we also conducted a parallel analysis of the regulators of sporophytic SI in Brassica, a system in which both the pistil and pollen genes are known. Although the pattern of variation shown by the pollen gene of the Brassica system is consistent with its role as a determinant of pollen specificity, that of 48A departs from expectation. Our analysis further suggests that recombination between 48A and S-RNase may have occurred during the interval spanned by the gene genealogy, another indication that 48A may not regulate SI expression in pollen.     

The implications of intergenic polymorphism for major histocompatibility complex evolution

A systematic survey of six intergenic regions flanking the human HLA-B locus in eight haplotypes reveals the regions to be up to 20 times more polymorphic than the reported average degree of human neutral polymorphism. Furthermore, the extent of polymorphism is directly related to the proximity to the HLA-B locus. Apparently linkage to HLA-B locus alleles, which are under balancing selection, maintains the neutral polymorphism of adjacent regions. For these linked polymorphisms to persist, recombination in the 200-kb interval from HLA-B to TNF must occur at a low frequency. The high degree of polymorphism found distal to HLA-B suggests that recombination is uncommon on both sides of the HLA-B locus. The least-squares estimate is 0.15% per megabase with an estimated range from 0.02 to 0.54%. These findings place strong restrictions on possible recombinational mechanisms for the generation of diversity at the HLA-B.     

Likelihood analysis of disequilibrium mapping, and related problems.

In this paper a theory is developed that provides the sampling distribution of alleles at a diallelic marker locus closely linked to a low-frequency allele that arose as a single mutant. The sampling distribution provides a basis for maximum-likelihood estimation of either the recombination rate, the mutation rate, or the age of the allele, provided that the two other parameters are known. This theory is applied to (1) the data of Hästbacka et al., to estimate the recombination rate between a locus associated with diastrophic dysplasia and a linked RFLP marker; (2) the data of Risch et al., to estimate the age of a presumptive allele causing idiopathic distortion dystonia in Ashkenazi jews; and (3) the data of Tishkoff et al., to estimate the date at which, at the CD4 locus, non-African lineages diverged from African lineages. We conclude that the extent of linkage disequilibrium can lead to relatively accurate estimates of recombination and mutation rates and that those estimates are not very sensitive to parameters, such as the population age, whose values are not known with certainty. In contrast, we also conclude that, in many cases, linkage disequilibrium may not lead to useful estimates of allele age, because of the relatively large degree of uncertainly in those estimates.     

Diversity and linkage of genes in the self-incompatibility gene family in Arabidopsis lyrata

We report studies of seven members of the S-domain gene family in Arabidopsis lyrata, a member of the Brassicaceae that has a sporophytic self-incompatibility (SI) system. Orthologs for five loci are identifiable in the self-compatible relative A. thaliana. Like the Brassica stigmatic incompatibility protein locus (SRK), some of these genes have kinase domains. We show that several of these genes are unlinked to the putative A. lyrata SRK, Aly13. These genes have much lower nonsynonymous and synonymous polymorphism than Aly13 in the S-domains within natural populations, and differentiation between populations is higher, consistent with balancing selection at the Aly13 locus. One gene (Aly8) is linked to Aly13 and has high diversity. No departures from neutrality were detected for any of the loci. Comparing different loci within A. lyrata, sites corresponding to hypervariable regions in the Brassica S-loci (SLG and SRK) and in comparable regions of Aly13 have greater replacement site divergence than the rest of the S-domain. This suggests that the high polymorphism in these regions of incompatibility loci is due to balancing selection acting on sites within or near these regions, combined with low selective constraints.     

The use of intraallelic variability for testing neutrality and estimating population growth rate

To better understand the forces affecting individual alleles, we introduce a method for finding the joint distribution of the frequency of a neutral allele and the extent of variability at closely linked marker loci (the intraallelic variability). We model three types of intraallelic variability: (a) the number of nonrecombinants at a linked biallelic marker locus, (b) the length of a conserved haplotype, and (c) the number of mutations at a linked marker locus. If the population growth rate is known, the joint distribution provides the basis for a test of neutrality by testing whether the observed level of intraallelic variability is consistent with the observed allele frequency. If the population growth rate is unknown but neutrality can be assumed, the joint distribution provides the likelihood of the growth rate and leads to a maximum-likelihood estimate. We apply the method to data from published data sets for four loci in humans. We conclude that the Delta32 allele at CCR5 and a disease-associated allele at MLH1 arose recently and have been subject to strong selection. Alleles at PAH appear to be neutral and we estimate the recent growth rate of the European population to be approximately 0.027 per generation with a support interval of (0.017-0.037). Four of the relatively common alleles at CFTR also appear to be neutral but DeltaF508 appears to be significantly advantageous to heterozygous carriers.     

Recent selection for self‐compatibility in a population of Leavenworthia alabamica

The evolution of self‐compatibility (SC) is the first step in the evolutionary transition in plants from outcrossing enforced by self‐incompatibility (SI) to self‐fertilization. In the Brassicaceae, SI is controlled by alleles of two tightly linked genes at the S‐locus. Despite permitting inbreeding, mutations at the S‐locus leading to SC may be selected if they provide reproductive assurance and/or gain a transmission advantage in a population when SC plants self‐ and outcross. Positive selection can leave a genomic signature in the regions physically linked to the focus of selection when selection has occurred recently. From an SC population of Leavenworthia alabamica with a known nonfunctional mutation at the S‐locus, we collected sequence data from a ～690 Kb region surrounding the S‐locus, as well as from regions not linked to the S‐locus. To test for recent positive selection acting at the S‐locus, we examined polymorphism and the site‐frequency spectra. Using forward simulations, we demonstrate that recent selection of the strength expected for SC at a locus formerly under balancing selection can generate patterns similar to those seen in our empirical data.     

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

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

Disentangling the effects of breakdown of self-incompatibility and transition to selfing in North American Arabidopsis lyrata

A breakdown of self‐incompatibility (SI) followed by a shift to selfing is commonly observed in the evolution of flowering plants. Both are expected to reduce the levels of heterozygosity and genetic diversity. However, breakdown of SI should most strongly affect the region of the SI locus (S‐locus) because of the relaxation of balancing selection that operates on a functional S‐locus, and a potential selective sweep. In contrast, a transition to selfing should affect the whole genome. We set out to disentangle the effects of breakdown of SI and transition to selfing on the level and distribution of genetic diversity in North American populations of Arabidopsis lyrata. Specifically, we compared sequence diversity of loci linked and unlinked to the S‐locus for populations ranging from complete selfing to fully outcrossing. Regardless of linkage to the S‐locus, heterozygosity and genetic diversity increased with population outcrossing rate. High heterozygosity of self‐compatible individuals in outcrossing populations suggests that SI is not the only factor preventing the evolution of self‐fertilization in those populations. There was a strong loss of diversity in selfing populations, which was more pronounced at the S‐locus. In addition, selfing populations showed an accumulation of derived mutations at the S‐locus. Our results provide evidence that beyond the genome‐wide consequences of the population bottleneck associated with the shift to selfing, the S‐locus of A. lyrata shows a specific signal either reflecting the relaxation of balancing selection or positive selection.     

Hill-Robertson interference is a minor determinant of variations in codon bias across Drosophila melanogaster and Caenorhabditis elegans genomes

According to population genetics models, genomic regions with lower crossing-over rates are expected to experience less effective selection because of Hill-Robertson interference (HRi). The effect of genetic linkage is thought to be particularly important for a selection of weak intensity such as selection affecting codon usage. Consistent with this model, codon bias correlates positively with recombination rate in Drosophila melanogaster and Caenorhabditis elegans. However, in these species, the G+C content of both noncoding DNA and synonymous sites correlates positively with recombination, which suggests that mutation patterns and recombination are associated. To remove this effect of mutation patterns on codon bias, we used the synonymous sites of lowly expressed genes that are expected to be effectively neutral sites. We measured the differences between codon biases of highly expressed genes and their lowly expressed neighbors. In D. melanogaster we find that HRi weakly reduces selection on codon usage of genes located in regions of very low recombination; but these genes only comprise 4% of the total. In C. elegans we do not find any evidence for the effect of recombination on selection for codon bias. Computer simulations indicate that HRi poorly enhances codon bias if the local recombination rate is greater than the mutation rate. This prediction of the model is consistent with our data and with the current estimate of the mutation rate in D. melanogaster. The case of C. elegans, which is highly self-fertilizing, is discussed. Our results suggest that HRi is a minor determinant of variations in codon bias across the genome.     

Self-incompatibility of Prunus tenella and evidence that reproductively isolated species of Prunus have different SFB alleles coupled with an identical S-RNase allele

Many species of Prunus display an S-RNase-based gametophytic self-incompatibility (SI), controlled by a single highly polymorphic multigene complex termed the S-locus. This comprises tightly linked stylar- and pollen-expressed genes that determine the specificity of the SI response. We investigated SI of Prunus tenella, a wild species found in small, isolated populations on the Balkan peninsula, initially by pollination experiments and identifying stylar-expressed RNase alleles. Nine P. tenella S-RNase alleles (S(1)-S(9)) were cloned; their sequence analysis showed a very high ratio of non-synonymous to synonymous nucleotide substitutions (K(a)/K(s)) and revealed that S-RNase alleles of P. tenella, unlike those of Prunus dulcis, show positive selection in all regions except the conserved regions and that between C2 and RHV. Remarkably, S(8)-RNase, was found to be identical to S(1)-RNase from Prunus avium, a species that does not interbreed with P. tenella and, except for just one amino acid, to S(11) of P. dulcis. However, the corresponding introns and S-RNase-SFB intergenic regions showed considerable differences. Moreover, protein sequences of the pollen-expressed SFB alleles were not identical, harbouring 12 amino-acid replacements between those of P. tenella SFB(8) and P. avium SFB(1). Implications of this finding for hypotheses about the evolution of new S-specificities are discussed.     

Gene duplication and evidence for balancing selection acting on MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis)

Allelic polymorphism and evolution mechanism of major histocompatibility complex (MHC) genes has been investigated in many mammals, however, much less is known in teleost. In order to investigate the mechanisms creating and maintaining variability at the MHC class II DAA locus, we examined the polymorphism, gene duplication and balancing selection of MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis). We described 33 alleles in the C. semilaevis, recombination and gene duplication seems to play more important roles in the origin of new alleles. The rate of non-synonymous substitutions (d(N)) occurred at a significantly higher frequency than that of synonymous substitutions (d(S)) in peptide-binding region (PBR) and non-PBR, suggesting balancing selection for maintaining polymorphisms at the MHC II DAA locus. Many positive selection sites were found to act very intensively on antigen-binding sites. Our founding suggests a snapshot in an evolutionary process of MHC-DAA gene evolution of the C. semilaevis.     

Detecting the signature of selection on immune genes in highly structured populations of wild sheep (Ovis dalli)

The confounding effects of population structure complicate efforts to identify regions of the genome under the influence of selection in natural populations. Here we test for evidence of selection in three genes involved in vertebrate immune function - the major histocompatibility complex (MHC), interferon gamma (IFNG) and natural resistance associated macrophage polymorphism (NRAMP) - in highly structured populations of wild thinhorn sheep (Ovis dalli). We examined patterns of variation at microsatellite loci linked to these gene regions and at the DNA sequence level. Simple Watterson's tests indicated balancing selection at all three gene regions. However, evidence for selection was confounded by population structure, as the Watterson's test statistics from linked markers were not outside of the range of values from unlinked and presumably neutral microsatellites. The translated coding sequences of thinhorn IFNG and NRAMP are fixed and identical to those of domestic sheep (Ovis aries). In contrast, the thinhorn MHC DRB locus shows significant evidence of overdominance through both an excess of nonsynonymous substitution and trans-species polymorphism. The failure to detect balancing selection at microsatellite loci linked to the MHC is likely the result of recombination between the markers and expressed gene regions.     

Evidence for positive selection at the pantophysin (Pan I) locus in walleye pollock, Theragra chalcogramma

Nucleotide polymorphism at the pantophysin (Pan I) locus in walleye pollock, Theragra chalcogramma, was examined using DNA sequence data. Two distinct allelic lineages were detected in pollock, resulting from three amino acid replacement mutations in the first intravesicular domain of the protein. The common Pan I allelic group, comprising 94% of the samples, was less polymorphic (pi = 0.005) than the uncommon group (pi = 0.008), and nucleotide diversity in both was higher than for two allelic lineages in the related Atlantic cod, Gadus morhua. Phylogenetic analyses of Pan I sequences from these two species did not clearly resolve orthology among allelic groups, in part because of recombination that has occurred between the two pollock lineages. Conventional tests of neutrality comparing polymorphisms within and between homologous regions of the Pan I locus in walleye pollock and Atlantic cod did not detect the effects of selection. This result is likely attributed to low levels of synonymous divergence among allelic lineages and a lack of mutation-drift equilibrium inferred from nucleotide mismatch frequency distributions. However, the ratio of nonsynonymous to synonymous substitutions per site (dN/dS) exceeded unity in two intravesicular domains of the protein and the influence of positive selection at multiple codon sites was strongly inferred through the use of maximum-likelihood analyses. In addition, the frequency spectrum of linked neutral variation showed indirect effects of adaptive hitchhiking in pollock resulting from a selective sweep of the common allelic lineage. Recombination between the two allelic classes may have prevented complete loss of the older, more polymorphic lineage. The results suggest that recurrent sweeps driven by positive selection is the principle mode of evolution at the Pan I locus in gadid fishes.     

Chromosome Walking in the Petunia Inflata Self-Incompatibility (S-) Locus and Gene Identification in an 881-kb Contig Containing S 2 -RNase

Self-incompatibility (SI) in the Solanaceae, Rosaceae and Scrophulariaceae is controlled by the polymorphic S locus, which contains two separate genes encoding pollen and pistil determinants in SI interactions. The S-RNase gene encodes the pistil determinant, whereas the pollen determinant gene, named the pollen S gene, has not yet been identified. Here, we set out to construct an integrated genetic and physical map of the S locus of Petunia inflata and identify any additional genes located at this locus. We first conducted chromosome walking at the S2 locus using BAC clones that contained either S2-RNase or one of the nine markers tightly linked to the S locus. Ten separate contigs were constructed, which collectively spanned 4.4 Mb. To identify additional genes located at the S2 locus, a 328-kb region (part of an 881-kb BAC contig) containing S2-RNase was completely sequenced. Approximately 76% of the region contained repetitive sequences, including transposon-like sequences. Other than S2-RNase, an F-box gene, named PiSLF2 (S2-allele of P. inflata S-locus F-box gene), was the only predicted gene whose deduced amino acid sequence was similar to the sequences of known proteins in the database. Two different cDNA selection methods were used to identify additional genes in the 881-kb contig; 11 groups of cDNA clones were identified in addition to those for S2-RNase and PiSLF2. RT-PCR analysis of expression profiles and PCR analysis of BAC clones and genomic DNA confirmed that seven of these 11 newly identified genes were located in the 881-kb contig.     

Microsatellite variation and recombination rate in the human genome

Background (purifying) selection on deleterious mutations is expected to remove linked neutral mutations from a population, resulting in a positive correlation between recombination rate and levels of neutral genetic variation, even for markers with high mutation rates. We tested this prediction of the background selection model by comparing recombination rate and levels of microsatellite polymorphism in humans. Published data for 28 unrelated Europeans were used to estimate microsatellite polymorphism (number of alleles, heterozygosity, and variance in allele size) for loci throughout the genome. Recombination rates were estimated from comparisons of genetic and physical maps. First, we analyzed 61 loci from chromosome 22, using the complete sequence of this chromosome to provide exact physical locations. These 61 microsatellites showed no correlation between levels of variation and recombination rate. We then used radiation-hybrid and cytogenetic maps to calculate recombination rates throughout the genome. Recombination rates varied by more than one order of magnitude, and most chromosomes showed significant suppression of recombination near the centromere. Genome-wide analyses provided no evidence for a strong positive correlation between recombination rate and polymorphism, although analyses of loci with at least 20 repeats suggested a weak positive correlation. Comparisons of microsatellites in lowest-recombination and highest-recombination regions also revealed no difference in levels of polymorphism. Together, these results indicate that background selection is not a major determinant of microsatellite variation in humans.