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.     

Uneven segregation of sporophytic self-incompatibility alleles in Arabidopsis lyrata

Self-incompatibility in Arabidopsis lyrata is sporophytically controlled by the multi-allelic S-locus. Self-incompatibility alleles (S-alleles) are under strong negative frequency dependent selection because pollen carrying common S-alleles have fewer mating opportunities. Population genetics theory predicts that deleterious alleles can accumulate if linked to the S-locus. This was tested by studying segregation of S-alleles in 11 large full sib families in A. lyrata. Significant segregation distortion leading to an up to fourfold difference in transmission rates was found in six families. Differences in transmission rates were not significantly different in reciprocal crosses and the distortions observed were compatible with selection acting at the gametic stage alone. The S-allele with the largest segregation advantage is also the most recessive, and is very common in natural populations concordant with its apparent segregation advantage. These results imply that frequencies of S-alleles in populations of A. lyrata cannot be predicted based on simple models of frequency-dependent selection alone.     

Inheritance and dominance of self-incompatibility alleles in polyploid Arabidopsis lyrata

Natural populations of diploid Arabidopsis lyrata exhibit the sporophytic type of self-incompatibility system characteristic of Brassicaceae, in which complicated dominance interactions among alleles in the diploid parent determine self-recognition phenotypes of both pollen and stigma. The purpose of this study was to investigate how polyploidy affects this already complex system. One tetraploid population (Arabidopsis lyrata ssp kawasakiana from Japan) showed complete self-compatibility and produced viable selfed progeny for at least three generations subsequent to field collection. In contrast, individuals from a second tetraploid population (A. lyrata ssp petraea from Austria) were strongly self-incompatible (SI). Segregation of SI genotypes in this population followed Mendelian patterns based on a tetrasomic model of inheritance, with two to four alleles per individual, independent segregation of alleles, and little evidence of dosage effects of alleles found in multiple copies. Similar to results from diploids, anomalous compatibility patterns involving particular combinations of individuals occurred at a low frequency in the tetraploids, suggesting altered dominance in certain genetic backgrounds that could be due to the influence of a modifier locus. Overall, dominance relationships among S-alleles in self-incompatible tetraploid families were remarkably similar to those in related diploids, suggesting that this very important and complicated locus has not undergone extensive modification subsequent to polyploidization.     

Linkage disequilibrium and recombination rate estimates in the self-incompatibility region of Arabidopsis lyrata

Genetic diversity is unusually high at loci in the S-locus region of the self-incompatible species of the flowering plant, Arabidopsis lyrata, not just in the S loci themselves, but also at two nearby loci. In a previous study of a single natural population from Iceland, we attributed this elevated polymorphism to linkage disequilibrium (LD) between variants at loci close to the S locus and the S alleles, which are maintained in the population by balancing selection. With the four S-flanking loci whose diversity we previously studied, we could not determine the extent of the region linked to the S loci in which neutral sites are affected. We also could not exclude the possibility of a population bottleneck, or of admixture, as causes of the LD. We have now studied four more distant loci flanking the S-locus region, and more populations, and we analyze the results using a theoretical model of the effect of balancing selection on diversity at linked neutral sites within and between different functional S-allelic classes. In the model, diversity is a function of the number of selectively maintained alleles and the recombination distances from the selectively maintained sites. We use the model to estimate the number of different functional S alleles, their turnover rate, and recombination rates between the S-locus region and other loci. Our estimates suggest that there is a small region of very low recombination surrounding the S-locus region.     

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

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

Testing for effects of recombination rate on nucleotide diversity in natural populations of Arabidopsis lyrata

We investigated DNA sequence diversity for loci on chromosomes 1 and 2 in six natural populations of Arabidopsis lyrata and tested for the role of natural selection in structuring genomewide patterns of variability, specifically examining the effects of recombination rate on levels of silent polymorphism. In contrast with theoretical predictions from models of genetic hitchhiking, maximum-likelihood-based analyses of diversity and divergence do not suggest reduction of diversity in the region of suppressed recombination near the centromere of chromosome 1, except in a single population from Russia, in which the pericentromeric region may have undergone a local selective sweep or demographic process that reduced variability. We discuss various possibilities that might explain why nucleotide diversity in most A. lyrata populations is not related to recombination rate, including genic recombination hotspots, and low gene density in the low recombination rate region.     

Conditional neutrality at two adjacent NBS-LRR disease resistance loci in natural populations of Arabidopsis lyrata

We examined patterns of nucleotide diversity at a genomic region containing two linked candidate disease resistance (NBS-LRR) genes in seven populations of the outcrossing plant Arabidopsis lyrata. In comparison with two adjacent control genes and neutral reference genes across the genome, the NBS-LRR genes exhibited elevated nonsynonymous variation and a large number of major-effect polymorphisms causing early stop codons and/or frameshift mutations. In contrast, analysis of synonymous diversity provided no evidence that the region was subject to long-term balancing selection or recent selective sweeps in any of the seven populations surveyed. Also in contrast with earlier surveys of one of these R genes, there was no evidence that the resistance genes or the major-effect mutations were subject to elevated differentiation between populations. We suggest that conditional neutrality in the absence of the corresponding pathogen, rather than long-term balancing selection or local adaptation, may in some circumstances be a significant cause of elevated functional polymorphism at R genes. In contrast with the R genes, analysis of diversity and differentiation at the flanking FERONIA locus showed high population divergence, suggesting local adaptation on this locus controlling male-female signalling during fertilization.     

Population genetic structure of Arabidopsis lyrata in Europe

Population genetic theory predicts that the self-incompatible and perennial herb, Arabidopsis lyrata, will have a genetic structure that differs from the self-fertilizing, annual Arabidopsis thaliana. We quantified the genetic structure for eight populations of A. lyrata ssp. petraea in historically nonglaciated regions of central Europe. Analysis of 20 microsatellite loci for 344 individuals demonstrated that, in accordance with predictions, diploid populations had high genome-wide heterozygosity (H(O) = 0.48; H(E) = 0.52), high within-population diversity (83% of total) compatible with mutation-drift equilibrium, and moderate differentiation among populations (F(ST) = 0.17). Within a single population, the vast majority of genetic variability (92%) was found at the smallest spatial scale (< 3 m). Although there was no evidence of biparental inbreeding or clonal propagation at this scale (F(IS) = 0.003), significant fine-scale spatial autocorrelation indicated localized gene flow presumably due to gravity dispersed seeds (Sp = 0.018). Limited gene flow between isolated population clusters (regions) separated by hundreds of kilometres has given rise to an isolation by distance pattern of diversification, with low, but significant, differentiation among regions (F(ST) = 0.05). The maintenance of geographically widespread polymorphisms and uniformly high diversity throughout central Europe is consistent with periglacial survival of A. lyrata ssp. petraea north of the Alps in steppe-tundra habitats during the last glacial maximum. As expected of northern and previously glaciated localities, A. lyrata in Iceland was genetically less diverse and highly differentiated from central Europe (H(E) = 0.37; F(ST) = 0.27).     

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.     

Local adaptation in European populations of Arabidopsis lyrata (Brassicaceae)

We studied local adaptation to contrasting environments using an organism that is emerging as a model for evolutionary plant biology-the outcrossing, perennial herb Arabidopsis lyrata subsp. petraea (Brassicaceae). With reciprocal transplant experiments, we found variation in cumulative fitness, indicating adaptive differentiation among populations. Nonlocal populations did not have significantly higher fitness than the local population. Experimental sites were located in Norway (alpine), Sweden (coastal), and Germany (continental). At all sites after one year, the local population had higher cumulative fitness, as quantified by survival combined with rosette area, than at least one of the nonlocal populations. At the Norwegian site, measurements were done for two additional years, and fitness differences persisted. The fitness components that contributed most to differences in cumulative fitness varied among sites. Relatively small rosette area combined with a large number of inflorescences produced by German plants may reflect differentiation in life history. The results of the current study demonstrate adaptive population differentiation in A. lyrata along a climatic gradient in Europe. The studied populations harbor considerable variation in several characters contributing to adaptive population differentiation. The wealth of genetic information available makes A. lyrata a highly attractive system also for examining the functional and genetic basis of local adaptation in plants.     

Gene, phenotype and function: GLABROUS1 and resistance to herbivory in natural populations of Arabidopsis lyrata

The molecular genetic basis of adaptive variation is of fundamental importance for evolutionary dynamics, but is still poorly known. Only in very few cases has the relationship between genetic variation at the molecular level, phenotype and function been established in natural populations. We examined the functional significance and genetic basis of a polymorphism in production of leaf hairs, trichomes, in the perennial herb Arabidopsis lyrata. Earlier studies suggested that trichome production is subject to divergent selection. Here we show that the production of trichomes is correlated with reduced damage from insect herbivores in natural populations, and using statistical methods developed for medical genetics we document an association between loss of trichome production and mutations in the regulatory gene GLABROUS1. Sequence data suggest that independent mutations in this regulatory gene have provided the basis for parallel evolution of reduced resistance to insect herbivores in different populations of A. lyrata and in the closely related Arabidopsis thaliana. The results show that candidate genes identified in model organisms provide a valuable starting point for analysis of the genetic basis of phenotypic variation in natural populations.     

Time scales of divergence and speciation among natural populations and subspecies of Arabidopsis lyrata (Brassicaceae)

? Premise of the study: Plant populations that face new environments adapt and diverge simultaneously, and both processes leave footprints in their genetic diversity. Arabidopsis lyrata is an excellent species for studying these processes. Pairs of populations and subspecies of A. lyrata represent different stages of divergence. These populations are also known to be locally adapted and display various stages of emerging reproductive isolation. ? Methods: We used nucleotide diversity data from 19 loci to estimate divergence times and levels of diversity among nine A. lyrata populations. Traditional distance-based methods and model-based clustering analysis were used to supplement pairwise coalescence-based analysis of divergence. ? Key results: Estimated divergence times varied from 130000 generations between North American and European subspecies to 39000 generations between central European and Scandinavian populations. In concordance with previous studies, the highest level of diversity was found in Central Europe and the lowest in North America and a diverged Russian Karhum?ki population. Local adaptation among Northern and central European populations has emerged during the last 39000 generations. Populations that are reproductively isolated by prezygotic mechanisms have been separated for a longer time period of ～70000 generations but still have shared nucleotide polymorphism. ? Conclusions: In A. lyrata, reproductively isolated populations started to diverge ～70000 generations ago and more closely related, locally adapted populations have been separate lineages for ～39000 generations. However, based on the posterior distribution of divergence times, the processes leading to reproductive isolation and local adaptation are likely to temporally coincide.     

Patterns of genetic diversity in outcrossing and selfing populations of Arabidopsis lyrata

Arabidopsis lyrata is normally considered an obligately outcrossing species with a strong self-incompatibility system, but a shift in mating system towards inbreeding has been found in some North American populations (subspecies A. lyrata ssp. lyrata). This study provides a survey of the Great Lakes region of Canada to determine the extent of this mating system variation and how outcrossing rates are related to current population density, geographical distribution, and genetic diversity. Based on variation at microsatellite markers (progeny arrays to estimate multilocus outcrossing rates and population samples to estimate diversity measures) and controlled greenhouse pollinations, populations can be divided into two groups: (i) group A, consisting of individuals capable of setting selfed seed (including autogamous fruit set in the absence of pollinators), showing depressed outcrossing rates (T(m) = 0.2-0.6), heterozygosity (H(O) = 0.02-0.06) and genetic diversity (H(E) = 0.08-0.10); and (ii) group B, consisting of individuals that are predominantly self-incompatible (T(m) > 0.8), require pollinators for seeds set, and showing higher levels of heterozygosity (H(O) = 0.13-0.31) and diversity (H(E) = 0.19-0.410). Current population density is not related to the shift in mating system but does vary with latitude. Restricted gene flow among populations was evident among all but two populations (F(ST) = 0.11-0.8). Group A populations were more differentiated from one another (F(ST) = 0.78) than they were from group B populations (F(ST) = 0.59), with 41% of the variation partitioned within populations, 47% between populations, and 12% between groups. No significant relationship was found between genetic and geographical distance. Results are discussed in the context of possible postglacial expansion scenarios in relation to loss of self-incompatibility.     

Estimating the number,frequency, and dominance of S-alleles in a natural population of Arabidopsis lyrata(Brassicaceae) with sporophytic control of self-incompatibility

In homomorphic plant self-incompatibility (SI) systems, large numbers of alleles may be maintained at a single Mendelian locus. Most estimators of the number of alleles present in natural populations are designed for gametophytic self-incompatibility systems (GSI) in which the recognition phenotype of the pollen is determined by its own haploid genotype. In sporophytic systems (SSI), the recognition phenotype of the pollen is determined by the diploid genotype of its parent, and dominance differs among alleles. We describe research aimed at estimates of S-allele numbers in a natural population of Arabidopsis lyrata (Brassicaceae), whose SSI system has recently been described. Using a combination of pollination studies and PCR-based identification of alleles at a locus equivalent to the Brassica SRK gene, we identified and sequenced 11 putative alleles in a sample of 20 individuals from different maternal seed sets. The pollination results indicate that we have not amplified all alleles that must be present. Extensive partial incompatibility, nonreciprocal compatibility differences, and evidence of weakened expression of SI in some genotypes, prevent us from determining the exact number of missing alleles based only on cross-pollination data. Although we show that none of the theoretical models currently proposed is completely appropriate for estimating the number of alleles in this system, we estimate that there are between 13 and 16 different S-alleles in our sample, probably between 16 and 25 alleles in the population, and discuss the relative frequency of alleles in relation to dominance.     

Effects of gene expression on molecular evolution in Arabidopsis thaliana and Arabidopsis lyrata

We analyzed the complete genome sequence of Arabidopsis thaliana and sequence data from 83 genes in the outcrossing A. lyrata, to better understand the role of gene expression on the strength of natural selection on synonymous and replacement sites in Arabidopsis. From data on tRNA gene abundance, we find a good concordance between codon preferences and the relative abundance of isoaccepting tRNAs in the complete A. thaliana genome, consistent with models of translational selection. Both EST-based and new quantitative measures of gene expression (MPSS) suggest that codon preferences derived from information on tRNA abundance are more strongly associated with gene expression than those obtained from multivariate analysis, which provides further support for the hypothesis that codon bias in Arabidopsis is under selection mediated by tRNA abundance. Consistent with previous results, analysis of protein evolution reveals a significant correlation between gene expression level and amino acid substitution rate. Analysis by MPSS estimates of gene expression suggests that this effect is primarily the result of a correlation between the number of tissues in which a gene is expressed and the rate of amino acid substitution, which indicates that the degree of tissue specialization may be an important determinant of the rate of protein evolution in Arabidopsis.     

Landscape structure, clonal propagation, and genetic diversity in Scandinavian populations of Arabidopsis lyrata (Brassicaceae)

Colonization history, landscape structure, and environmental conditions may influence patterns of neutral genetic variation because of their effects on gene flow and reproductive mode. We compared variation at microsatellite loci within and among 26 Arabidopsis lyrata populations in two disjunct areas of its distribution in northern Europe (Norway and Sweden). The two areas probably share a common colonization history but differ in size (Norwegian range markedly larger than Swedish range), landscape structure (mountains vs. coast), and habitat conditions likely to affect patterns of gene flow and opportunities for sexual reproduction. Within-population genetic diversity was not related to latitude but was higher in Sweden than in Norway. Population differentiation was stronger among Norwegian than among Swedish populations (F(ST) = 0.23 vs. F(ST) = 0.18). The frequency of clonal propagation (proportion of identical multilocus genotypes) increased with decreasing population size, was higher in Norwegian than in Swedish populations, but was not related to altitude or substrate. Differences in genetic structure are discussed in relation to population characteristics and range size in the two areas. The results demonstrate that the possibility of clonal propagation should be considered when developing strategies for sampling and analyzing data in ecological and genetic studies of this emerging model species.     

Highly structured nucleotide variation within and among Arabidopsis lyrata populations at the FAH1 and DFR gene regions

Nucleotide variation at the FAH1 and DFR gene regions was surveyed in four populations of Arabidopsis lyrata (two European A. l. petraea and two North American A. l. lyrata populations). In contrast to previous results, levels of variation were not consistently lower in A. l. lyrata than in A. l. petraea, and similar degrees of genetic differentiation were detected between and within subspecies. These observations and the significant genetic differentiation detected among populations suggest population substructure and no real subdivision between subspecies. For each gene studied, genotypic data were obtained, which allowed comparing nucleotide diversity within individuals (between sequences from the same individual) and within populations (between sequences from the same population). The generally lower level of variation within than among individuals detected in each population yielded a significant deviation from panmixia within populations. In three of the four populations studied, two highly divergent alleles were detected within populations at the highly variable DFR locus. This pattern and the significant excess of derived variants detected in most populations suggest that most variation segregating within populations results from rare migration events between relatively small and isolated populations exhibiting reduced panmixia.