Multiple losses of self‐incompatibility in North‐American Arabidopsis lyrata?: Phylogeographic context and population genetic consequences

Arabidopsis lyrata is mostly outcrossing due to a sporophytic self‐incompatibility (SI) system but around the Great Lakes of North America some populations have experienced a loss of SI. We researched the loss of SI in a phylogeographic context. We used cpDNA and microsatellite markers to test if populations of North‐American A. lyrata around the Great Lakes have experienced different (recent) histories, and linked this with individually established selfing phenotype and population level realized outcrossing rates calculated based on variation in progeny arrays at multi‐locus microsatellite markers. We found three chloroplast haplotypes, in two of which the loss of self‐incompatibility had occurred independently. Shifts to high rates of inbreeding were most apparent in one of these lineages but individuals showing loss of SI occurred in all three. Self‐compatible individuals usually showed a reduction of observed heterozygosity (H_O) compared to outcrossing individuals. In the lineage that included the populations with the highest levels of inbreeding, this reduction was more substantial. This may indicate that the loss of SI in this lineage did not occur as recently as in the other lineage. The geographic distribution of the haplotypes suggested that there had been at least two independent colonization routes to the north of the Great Lakes following the last glaciation. This is consistent with postglacial migration patterns that have been suggested for other organisms with limited dispersal, such as reptiles and amphibians.     

The effect of mating system on growth of Arabidopsis lyrata in response to inoculation with the biotrophic parasite Albugo candida

The effects of variation in host reproductive systems on response to pathogens are not well understood. We inoculated individuals from outcrossing and inbreeding populations of North American Arabidopsis lyrata with Albugo candida (white blister rust) to test the effect of mating system and heterozygosity on disease response. We observed three host infection phenotypes, classified as fully resistant, partially resistant and fully susceptible. Overall, inbreeding populations had more susceptible and fewer partially resistant individuals than outcrossing populations, but the highest proportion of resistant individuals was found in two of the inbreeding populations. Mating system did not affect relative growth rate of inoculated plants, but there were strong effects of population and infection phenotype. We conclude that mating system per se does not determine the resistance of natural A. lyrata populations to infection by Albugo, but that the increased variability in responses among inbreeding populations may be due to reduced effective population size.     

Breakdown of self-incompatibility in the perennial Arabidopsis lyrata (Brassicaceae) and its genetic consequences

Mating systems in plants are known to be highly labile traits, with frequent transitions from outcrossing to selfing. The genetic basis for breakdown in self-incompatibility (SI) systems has been studied, but data on variation in selfing rates in species for which the molecular basis of SI is known are rare. This study surveyed such variation in Arabidopsis lyrata (Brassicaceae), which is often considered an obligately outcrossing species, to examine the causes and genetic consequences of changes in its breeding system. Based on controlled self-pollinations in the greenhouse, three populations from the Great Lakes region of North America included a minority of self-compatible (SC) individuals, while two showed larger proportions of SC individuals and all populations contained some individuals capable of setting selfed seeds. Loss of SI was not associated with particular haplotypes at the S-locus (as estimated by alleles amplified at the SRK locus, the gene controlling female specificity) and all populations contained similar numbers of SRK alleles, suggesting that some other genetic factor is responsible for modifying the SI reaction. The loss of SI has resulted in an effective shift in the mating system, as the two populations with a high frequency of SC individuals showed significantly lower microsatellite-based multilocus outcrossing rates and higher inbreeding coefficients than the other populations. Based on microsatellites, observed heterozygosities and genetic diversity were also significantly depressed in these populations. These findings provide the unique opportunity to examine in detail the consequences of mating system changes within a species with a well-characterized SI system.     

Limited phenological and pollinator-mediated isolation among selfing and outcrossing Arabidopsis lyrata populations

Transitions from outcrossing to selfing have been a frequent evolutionary shift in plants and clearly play a role in species divergence. However, many questions remain about the initial mechanistic basis of reproductive isolation during the evolution of selfing. For instance, how important are pre-zygotic pre-pollination mechanisms (e.g. changes in phenology and pollinator visitation) in maintaining reproductive isolation between newly arisen selfing populations and their outcrossing ancestors? To test whether changes in phenology and pollinator visitation isolate selfing populations of Arabidopsis lyrata from outcrossing populations, we conducted a common garden experiment with plants from selfing and outcrossing populations as well as their between-population hybrids. Specifically, we asked whether there was isolation between outcrossing and selfing plants and their between-population hybrids through differences in (1) the timing or intensity of flowering; and/or (2) pollinator visitation. We found that phenology largely overlapped between plants from outcrossing and selfing populations. There were also no differences in pollinator preference related to mating system. Additionally, pollinators preferred to visit flowers on the same plant rather than exploring nearby plants, creating a large opportunity for self-fertilization. Overall, this suggests that pre-zygotic pre-pollination mechanisms do not strongly reproductively isolate plants from selfing and outcrossing populations of Arabidopsis lyrata.     

Subdivision and haplotype structure in natural populations of Arabidopsis lyrata

High levels of inbreeding are expected to cause a strong reduction in levels of genetic variability, effective recombination rates and in adaptation compared with related outcrossing populations. We examined patterns of DNA polymorphism at five nuclear loci and one chloroplast locus within and between four populations of the outcrossing plant Arabidopsis lyrata, a close relative of the highly self-fertilizing model species A. thaliana. The observed patterns are compared with species-wide polymorphism at orthologous loci, as well as within- and between-population patterns at other studied loci in A. thaliana. In addition to evidence for much higher average within-population diversity, species-wide levels of silent polymorphism are generally higher in A. lyrata than in A. thaliana, unlike the results from a previous study of the ADH locus. However, polymorphism is also low in the North American A. lyrata subspecies lyrata compared with the European subspecies petraea, suggesting either a population bottleneck in North American populations or recent admixture involving diverged European populations. Differentiation between the two subspecies is strong, although there are few fixed differences, suggesting that their isolation is recent. Estimates of intralocus recombination rates and analysis of haplotype structure in European A. lyrata populations indicate lower recombination than predicted based on the variability together with physical recombination rates estimated from A. thaliana. This may be due to strong population subdivision, or to recent departures from demographic equilibrium such as a bottleneck or population admixture. Alternatively, there may be consistently lower recombination rates in the outcrossing species. In contrast, estimates of recombination rates from species-wide samples of A. thaliana are close to the values expected assuming a high rate of self-fertilization. Complex population histories in both A. thaliana and A. lyrata complicate theoretical predictions and empirical tests of the effects of inbreeding on polymorphism and molecular evolution.     

Relatively weak inbreeding depression in selfing but also in outcrossing populations of North American Arabidopsis lyrata

Hermaphroditic plants can potentially self‐fertilize, but most possess adaptations that promote outcrossing. However, evolutionary transitions to higher selfing rates are frequent. Selfing comes with a transmission advantage over outcrossing, but self‐progeny may suffer from inbreeding depression, which forms the main barrier to the evolution of higher selfing rates. Here, we assessed inbreeding depression in the North American herb Arabidopsis lyrata, which is normally self‐incompatible, with a low frequency of self‐compatible plants. However, a few populations have become fixed for self‐compatibility and have high selfing rates. Under greenhouse conditions, we estimated mean inbreeding depression per seed (based on cumulative vegetative performance calculated as the product of germination, survival and aboveground biomass) to be 0.34 for six outcrossing populations, and 0.26 for five selfing populations. Exposing plants to drought and inducing defences with jasmonic acid did not magnify these estimates. For outcrossing populations, however, inbreeding depression per seed may underestimate true levels of inbreeding depression, because self‐incompatible plants showed strong reductions in seed set after (enforced) selfing. Inbreeding‐depression estimates incorporating seed set averaged 0.63 for outcrossing populations (compared to 0.30 for selfing populations). However, this is likely an overestimate because exposing plants to 5% CO₂ to circumvent self‐incompatibility to produce selfed seed might leave residual effects of self‐incompatibility that contribute to reduced seed set. Nevertheless, our estimates of inbreeding depression were clearly lower than previous estimates based on the same performance traits in outcrossing European populations of A. lyrata, which may help explain why selfing could evolve in North American A. lyrata.     

Effect of inbreeding depression on outcrossing rates among populations of a tropical pine

Inbreeding depression is common among plants and may distort mating system estimates. Mating system studies traditionally ignore this effect, nonetheless an assessment of inbreeding depression that may have occurred before progeny evaluation could be necessary. In the neotropical Pinus chiapensis inbreeding depression was evaluated using regression analysis relating progeny F-values with seed germinability, the mating system was analysed in three populations with contrasting size, using isozymes, obtained a corrected outcrossing rate. Selfing decreased seed viability by 19%, relative to an outcrossed plant. Multilocus outcrossing rates, t(m), varied widely among populations. In the two smallest populations t(m) congruent with 1. Therefore, inbreeding depression did not affect the estimates, but overestimated t(m) by 10% in the third population, which has a true mixed mating system (selfing was the major source of inbreeding), and an unusually low t(m) for pines (t(m) = 0.54, uncorrected, t(m) = 0.49, corrected). Inbreeding depression may be an uneven source of bias for outcrossing estimates even at the infraspecific level. Accuracy [corrected] but not precision [corrected] may be gained by including inbreeding depression in outcrossing estimates. Therefore, caution should be taken when comparing t(m) among species or even populations within the same 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.     

Natural variation in Arabidopsis lyrata vernalization requirement conferred by a FRIGIDA indel polymorphism

Species share homologous genes to a large extent, but it isnot yet known to what degree the same loci have been targetsfor natural selection in different species. Natural variationin flowering time is determined to a large degree by 2 genes,FLOWERING LOCUS C and FRIGIDA, in Arabidopsis thaliana. Here,we examine whether FRIGIDA has a role in differences in floweringtime between and within natural populations of Arabidopsis lyrata,a close outcrossing perennial relative of A. thaliana. We found2 FRIGIDA sequence variants producing potentially functionalproteins but with a length difference of 14 amino acids. Thesevariants conferred a 15-day difference in flowering time inan association experiment in 2 Scandinavian populations. Thedifference in flowering time between alleles was confirmed withtransformation to A. thaliana. Because the north European late-floweringpopulations harbor both late- and early sequence variants atintermediate frequencies and the late-flowering variant is mostfrequent in the southern early flowering European population,other genetic factors must be responsible for the floweringtime differences between the populations. The length polymorphismoccurs at high frequencies also in several North American populations.The occurrence of functional variants at intermediate frequenciesin several populations suggests that the variation may be maintainedby balancing selection. This is in contrast to A. thaliana,where independent loss-of-function mutations at the FRIGIDAgene are responsible for differences between populations andlocal adaptation.     

Genetic diversity and population structure in the outcrossing populations of Equisetum arvense and E. hyemale (Equisetaceae)

An allozyme examination was conducted to study the mating systems and genetic differentiation of populations of Equisetum arvense and E. hyemale. The study revealed that the rate of intragametophytic selfing in these homosporous pteridophytes is very low, i.e., on average 0.020 and 0.019, respectively, despite the potential hermaproditism and selfing of the gametophytes. Most populations consisted of numerous genotypes, and the average heterozygosities of E. arvense and E. hyemale equalled 0.092 and 0.134, respectively. The commonly observed excess of the heterozygote genotypes indicates that there are interclonal differences in the frequency of vegetative reproduction. The level of genetic divergence among populations was considerable even within a limited geographic area. It is suggested that the life history of Equisetum, characterized by the inefficiency of spore germination and gametophyte reproduction in noncolonizing situations, limits the level of gene flow and leads to a great genetic divergence between populations.     

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.     

The allopolyploid Arabidopsis kamchatica originated from multiple individuals of Arabidopsis lyrata and Arabidopsis halleri

Polyploidization, or genome duplication, has played a critical role in the diversification of animals, fungi and plants. Little is known about the population structure and multiple origins of polyploid species because of the difficulty in identifying multiple homeologous nuclear genes. The allotetraploid species Arabidopsis kamchatica is closely related to the model species Arabidopsis thaliana and is distributed in a broader climatic niche than its parental species. Here, we performed direct sequencing of homeologous pairs of the low-copy nuclear genes WER and CHS by designing homeolog-specific primers, and obtained also chloroplast and ribosomal internal transcribed spacer sequences. Phylogenetic analysis showed that 50 individuals covering the distribution range including North America are allopolyploids derived from Arabidopsis lyrata and Arabidopsis halleri . Three major clusters within A. kamchatica were detected using Bayesian clustering. One cluster has widespread distribution. The other two are restricted to the southern part of the distribution range including Japan, where the parent A. lyrata is not currently distributed. This suggests that the mountains in Central Honshu and surrounding areas in Japan served as refugia during glacial–interglacial cycles and retained this diversity. We also found that multiple haplotypes of nuclear and chloroplast sequences of A. kamchatica are identical to those of their parental species. This indicates that multiple diploid individuals contributed to the origin of A. kamchatica . The haplotypes of low-copy nuclear genes in Japan suggest independent polyploidization events rather than introgression. Our findings suggest that self-compatibility and gene silencing occurred independently in different origins.     

High selfing and high inbreeding depression in peripheral populations of Juncus atratus

The mating system of a plant is the prime determinant of its population genetic structure. However, mating system effects may be modified by postzygotic mechanisms like inbreeding depression. Furthermore, historical as well as contemporary ecological factors and population characteristics, like the location within the species range can contribute to genetic variability. Using microsatellite markers we assessed the population genetic structure of the wind-pollinated Juncus atratus in 16 populations from peripheral and nearly central areas of the distribution range and studied the mating system of the species. In three peripheral populations, outcrossing rates at seeds stage were low (mean t(m) = 5.6%), suggesting a highly autogamous mating system. Despite this fact, on adult stage both individual heterozygosity (mean H(O) = 0.48) and gene diversity (mean H(E) = 0.58) were high even in small populations. Inbreeding coefficients were consistently low among all populations (mean F(IS) = 0.15). Within the three peripheral populations indirect estimates of lifetime inbreeding depression were surprisingly high (delta(eq) = 0.96) and inbreeding depression could be shown to act mostly on early seedling establishment. Similar conditions of autogamy combined with high inbreeding depression are typical for plants with a large lifetime genomic mutation rate that cannot avoid selfing by geitonogamy. However, the results presented here are unexpected for small-statured, herbaceous plants. Substantial genetic differentiation among all populations was found (mean F(ST) = 0.24). An isolation-by-distance pattern was apparent on large scale but not on local scale suggesting that the overall pattern was largely influenced by historical factors, e.g. colonization, whereas locally genetic drift was of greater importance than gene flow. Peripheral populations exhibited lower genetic diversity and higher inbreeding coefficients when compared with subcentral populations.     

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).     

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.     

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.     

Evolutionary dynamics of mating system shifts in Arabidopsis lyrata

Outcrossing is the prevalent mode of reproduction in plants and animals despite its substantial costs, while selfing and mixed mating occur at much lower frequency. Comparative research on plants has demonstrated the lability of self‐incompatibility, but there is little information about the transition on a within‐species level from self‐incompatibility to predominant selfing. We studied variation in mating system among 18 populations of Arabidopsis lyrata within a phylogenetic context to shed light on the evolution of selfing. Realized and potential mating systems were assessed by genetic analysis with microsatellite markers and hand‐self‐pollinations on 30 plants from each population. The fraction of self‐incompatible plants in a population was highly correlated with the outcrossing rate, showing that the spread of self‐compatibility is accompanied by or soon followed by an increase in the rate of selfing. The four predominantly selfing populations (outcrossing rates < 0.25) fell into more than one phylogenetic cluster, suggesting that the transition to selfing occurred more than once independently. Hence, A. lyrata offers an opportunity for the comparative analysis of outcrossing as a predominant mode of reproduction in plants and of the causes of the shift to selfing.     

Development of EST‐derived microsatellite markers for Arabidopsis lyrata subspecies petraea (L.)

A set of expressed sequence tag–simple sequence repeat (EST‐SSR) loci has been developed for Arabidopsis lyrata ssp. petraea. From 768 root cDNA clones, 126 microsatellites, including di‐, tri‐, tetra‐ and pentanucleotide repeat motifs were identified and primers were designed to 24 EST‐SSRs. Eleven loci were subsequently screened on 150 individuals sampled from five natural populations, which revealed three to nine alleles per locus (mean 5.36) and expected heterozygosity (H_E) estimates ranging from 0.046 to 0.698. Significant deviations from random mating were observed at 10 EST‐SSR loci, likely due to inbreeding (global F_IS = 0.151) and population structure (global F_ST = 0.246).     

用遗传标记监测不同环境和经营措施下的加勒比松森林群体遗传动态

对古巴加勒比松的6个群体（包括天然林、采伐林、母树林和种子园）进行了同功酶分析,根据5个酶系统8个位点的同功酶数据,对各群 交配系统以及群体遗传变异和结构进行了分析,天然林、种子园和母树林的多位点异交率和绝大多数单位点异交率都和完全异交无显著差异,过渡采伐的松树岛群体多位点异交率显著小于完全异交,而只有一半单位点异交率显著小于完全异交,而且该群体单位点平均异交率和多位点异交率均低于其它3个群体的估计值。采伐群体中同功酶变异和基因多样性与天然林群体JAG的相似,但低于其它群体,其近交系数较大,但小于天然林MAN和中国栽培群体的近交系数。中国引种栽培群体无论是同功酶变异还是基因多样性都显著高于古巴群体,与所有古巴群体的遗传距离都显著大于古巴群体之间的遗传距离。结果表明过度采伐导致群体自交程度增加,营建种子园可有效减少近交。自然分布区以外的引种栽培群体遗传变化量大,无论遗传变异和基因多样性都比参试其它群体大。