POPULATION STRUCTURE IN ARABIDOPSIS LYRATA: EVIDENCE FOR DIVERGENT SELECTION ON TRICHOME PRODUCTION

Intense browsing by sika deer (Cervus nippon) drives the genetic differentiation of hairy nettle (Urtica thunbergiana) populations

Many studies have inferred the way in which natural selection, genetic drift and gene flow shape the population genetic structures, but very few have quantified the population differentiation under spatially and temporally varying levels of selection pressure, population fluctuation and gene flow. In Nara Park (6.6 km²), central Japan, where several hundred sika deer (Cervus nippon) have been protected for more than 1,200 years, heavily- or moderately-haired nettle (Urtica thunbergiana) populations have evolved probably in response to intense deer browsing. Here, we analysed the genetic structure of two Nara Park populations and five surrounding populations using amplified fragment length polymorphism markers. A total of 546 marker loci were genotyped from 210 individuals. A Bayesian method estimated 5.5% of these loci to be outliers, which are putatively under natural selection. Neighbour-joining, principal coordinates and Bayesian clustering analyses using all-loci, non-outlier loci and outlier loci datasets showed that the Nara Park populations formed a cluster distinct from the surroundings. These results indicate the genome-wide differentiation of the Nara Park populations from the surroundings. Moreover, these imply the following: (1) gene flow is limited between these populations and thus genetic drift is a major factor causing the differentiation; and (2) natural selection imposed by intense deer browsing has contributed to some extent to the differentiation. In conclusion, sika deer seems to have counteracted genetic drift to drive the genetic differentiation of hairy nettles in Nara Park. This study suggests that a single herbivore species could lead to genetic differentiation among plant populations.

Leaf Trichome Distribution Pattern in Arabidopsis Reveals Gene Expression Variation Associated with Environmental Adaptation

Gene expression varies stochastically even in both heterogenous and homogeneous cell populations. This variation is not simply useless noise; rather, it is important for many biological processes. Unicellular organisms or cultured cell lines are useful for analyzing the variation in gene expression between cells; however, owing to technical challenges, the biological relevance of this variation in multicellular organisms such as higher plants remain unclear. Here, we addressed the biological relevance of this variation between cells by examining the genetic basis of trichome distribution patterns in Arabidopsis thaliana. The distribution pattern of a trichome on a leaf is stochastic and can be mathematically represented using Turing’s reaction-diffusion (RD) model. We analyzed simulations based on the RD model and found that the variability in the trichome distribution pattern increased with the increase in stochastic variation in a particular gene expression. Moreover, differences in heat-dependent variability of the trichome distribution pattern between the accessions showed a strong correlation with environmental factors to which each accession was adapted. Taken together, we successfully visualized variations in gene expression by quantifying the variability in the Arabidopsis trichome distribution pattern. Thus, our data provide evidence for the biological importance of variations in gene expression for environmental adaptation.

Ontogenetic changes in the targets of natural selection in three plant defenses

The evolution of plant defenses has traditionally been studied at single plant ontogenetic stages, overlooking the fact that natural selection acts continuously on organisms along their development, and that the adaptive value of phenotypes can change along ontogeny. We exposed 20 replicated genotypes of Turnera velutina to field conditions to evaluate whether the targets of natural selection on different defenses and their adaptative value change across plant development. We found that low chemical defense was favored in seedlings, which seems to be explained by the assimilation efficiency and the ability of the specialist herbivore to sequester cyanogenic glycosides. Whereas trichome density was unfavored in juvenile plants, it increased relative plant fitness in reproductive plants. At this stage we also found a positive correlative gradient between cyanogenic potential and sugar content in extrafloral nectar. We visualize this complex multi-trait combination as an ontogenetic defensive strategy. The inclusion of whole-plant ontogeny as a key source of variation in plant defense revealed that the targets and intensity of selection change along the development of plants, indicating that the influence of natural selection cannot be inferred without the assessment of ontogenetic strategies in the expression of multiple defenses.

Local adaptation and ecological differentiation under selection, migration, and drift in Arabidopsis lyrata

How the balance between selection, migration, and drift influences the evolution of local adaptation has been under intense theoretical scrutiny. Yet, empirical studies that relate estimates of local adaptation to quantification of gene flow and effective population sizes have been rare. Here, we conducted a reciprocal transplant trial, a common garden trial, and a whole-genome-based demography analysis to examine these effects among Arabidopsis lyrata populations from two altitudinal gradients in Norway. Demography simulations indicated that populations within the two gradients are connected by gene flow (0.1 < 4Ne m < 11) and have small effective population sizes (Ne <  6000), suggesting that both migration and drift can counteract local selection. However, the three-year field experiments showed evidence of local adaptation at the level of hierarchical multiyear fitness, attesting to the strength of differential selection. In the lowland habitat, local superiority was associated with greater fecundity, while viability accounted for fitness differences in the alpine habitat. We also demonstrate that flowering time differentiation has contributed to adaptive divergence between these locally adapted populations. Our results show that despite the estimated potential of gene flow and drift to hinder differentiation, selection among these A. lyrata populations has resulted in local adaptation.

Role of seed germination in adaptation and reproductive isolation in Arabidopsis lyrata

Seed germination is an important developmental and life history stage. Yet, the evolutionary impact of germination has mainly been studied in the context of dormancy, or for its role in reproductive isolation between species. Here, we aim to examine multiple consequences of genetic divergence on germination traits between two Arabidopsis lyrata subspecies: ssp. petraea (Eurasia) and ssp. lyrata (North America). Postdormancy germination time, a potentially adaptive trait, showed differentiation between the populations, and quantitative trait loci (QTL) mapping revealed that the trait variation is mainly controlled by two antagonistic loci. These QTL areas contain several candidate genes with known function in postdormancy germination in A. thaliana. The sequence variation of three genes was consistent with differential selection, and they also included fixed nonsynonymous substitutions with potential to account for the phenotypic differentiation. We further show that the divergence between the subspecies has led to a slight but significant reduction in hybrid germination proportions, indicating incipient reproductive isolation. Comparison of reciprocal F₁ and F₂ progenies suggests that Bateson-Dobzhansky-Muller incompatibilities likely act through uniparentally inherited factors. Examination of genomewide transmission ratio distortion further revealed that cytonuclear interactions cause substantial pregermination inviability in the hybrids. These results confirm that seed germination has adaptive potential beyond the dormancy stage and that hybrid seed inviability can be one of the first reproductive barriers to arise during divergence.

Fine-scale frequency differentiation along a herbivory gradient in the trichome dimorphism of a wild Arabidopsis

Geographic variation is commonly observed in plant resistance traits, where plant species might experience different selection pressure across a heterogeneous landscape. Arabidopsis halleri subsp. gemmifera is dimorphic for trichome production, generating two morphs, trichome-producing (hairy) and trichomeless (glabrous) plants. Trichomes of A. halleri are known to confer resistance against the white butterfly, cabbage sawfly, and brassica leaf beetle, but not against flea beetles. We combined leaf damage, microclimate, and microsatellite loci data of 26 A. halleri populations in central Japan, to explore factors responsible for fine-scale geographic variation in the morph frequency. We found that hairy plants were less damaged than glabrous plants within populations, but the among-site variation was the most significant source of variation in the individual-level damage. Fixation index (Gst″) of a putative trichome locus exhibited a significant divergence along population-level damage with an exception of an outlier population, inferring the local adaptation to herbivory. Notably, this outlier was a population wherein our previous study reported a balancing role of the brassica leaf beetle Phaedon brassicae on the morph frequency. This differentiation of the trichome locus was unrelated to neutral genetic differentiation (evaluated by Gst″ of microsatellite loci) and meteorological factors (including temperature and solar radiation). The present findings, combined with those of our previous work, provide suggestive evidence that herbivore-driven divergence and occasional outbreak of a specific herbivore have jointly contributed to the ecogeographic pattern in the frequency of two morphs.

Difference in evolutionary patterns of strongly or weakly selected characters among ant populations

Despite being a central issue in evolutionary biology, few studies have examined the stasis of characters in populations with no gene flow. A possible mechanism of such stasis is stabilizing selection with similar peaks in each population. This study examined the evolutionary patterns of morphological characters with and without strong selection in ant populations. We show that compared to a character that seems to be less important, characters that are more important were less variable within and among populations. Microsatellite analyses showed significant genetic differences between populations, implying limited gene flow between them. The observed levels of genetic differentiation cannot be attributed to recent population separations. Thus, the observed differences in morphological variance seem to reflect the degree of selection on each character. The less important character changed proportionately with time, but such a pattern was not observed in more important characters. These results suggest that stabilizing selection maintains morphological stasis between populations of the same species with minimal gene flow independent of divergence times.

The genetics of extreme microgeographic adaptation: an integrated approach identifies a major gene underlying leaf trichome divergence in Yellowstone Mimulus guttatus

Microgeographic adaptation provides a particularly interesting context for understanding the genetic basis of phenotypic divergence and may also present unique empirical challenges. In particular, plant adaptation to extreme soil mosaics may generate barriers to gene flow or shifts in mating system that confound simple genomic scans for adaptive loci. Here, we combine three approaches - quantitative trait locus (QTL) mapping of candidate intervals in controlled crosses, population resequencing (PoolSeq) and analyses of wild recombinant individuals - to investigate one trait associated with Mimulus guttatus (yellow monkeyflower) adaptation to geothermal soils in Yellowstone National Park. We mapped a major QTL causing dense leaf trichomes in thermally adapted plants to a <50-kb region of linkage Group 14 (Tr14) previously implicated in trichome divergence between independent M. guttatus populations. A PoolSeq scan of Tr14 region revealed a cluster of six genes, coincident with the inferred QTL peak, with high allele frequency differences sufficient to explain observed phenotypic differentiation. One of these, the R2R3 MYB transcription factor Migut.N02661, is a plausible functional candidate and was also strongly associated (r²  = 0.27) with trichome phenotype in analyses of wild-collected admixed individuals. Although functional analyses will be necessary to definitively link molecular variants in Tr14 with trichome divergence, our analyses are a major step in that direction. They point to a simple, and parallel, genetic basis for one axis of Mimulus guttatus adaptation to an extreme habitat, suggest a broadly conserved genetic basis for trichome variation across flowering plants and pave the way for further investigations of this challenging case of microgeographic incipient speciation.

Contrasting patterns of genetic structuring in natural populations of Arabidopsis lyrata Subsp. petraea across different regions in northern Europe

Level and partitioning of genetic diversity is expected to vary between contrasting habitats, reflecting differences in strength of ecological and evolutionary processes. Therefore, it is necessary to consider processes acting on different time scales when trying to explain diversity patterns in different parts of species' distributions. To explore how historical and contemporary factors jointly may influence patterns of genetic diversity and population differentiation, we compared genetic composition in the perennial herb Arabidopsis lyrata ssp. petraea from the northernmost parts of its distribution range on Iceland to that previously documented in Scandinavia. Leaf tissue and soil were sampled from ten Icelandic populations of A. lyrata. Seedlings were grown from soil samples, and tissue from above-ground and seed bank individuals were genotyped with 21 microsatellite markers. Seed bank density in Icelandic populations was low but not significantly different from that observed in Norwegian populations. While within-population genetic diversity was relatively high on Iceland (H(E) = 0.35), among-population differentiation was low (F(ST) = 0.10) compared to Norwegian and Swedish populations. Population differentiation was positively associated with geographical distance in both Iceland and Scandinavia, but the strength of this relationship varied between regions. Although topography and a larger distribution range may explain the higher differentiation between mountainous Norwegian relative to lowland populations in Sweden, these factors cannot explain the lower differentiation in Icelandic compared to Swedish populations. We propose that low genetic differentiation among Icelandic populations is not caused by differences in connectivity, but is rather due to large historical effective population sizes. Thus, rather than contemporary processes, historical factors such as survival of Icelandic lineages in northern refugia during the last glacial period may have contributed to the observed pattern.

The genetic architecture of constitutive and induced trichome density in two new recombinant inbred line populations of Arabidopsis thaliana: phenotypic plasticity, epistasis, and bidirectional leaf damage response

BACKGROUND

Herbivory imposes an important selective pressure on plants. In Arabidopsis thaliana leaf trichomes provide a key defense against insect herbivory; however, trichome production incurs a fitness cost in the absence of herbivory. Previous work on A. thaliana has shown an increase in trichome density in response to leaf damage, suggesting a mechanism by which the cost associated with constitutively high trichome density might be mitigated; however, the genetic basis of trichome density induction has not been studied.

RESULTS

Here, we describe the mapping of quantitative trait loci (QTL) for constitutive and damage induced trichome density in two new recombinant inbred line populations of A. thaliana; mapping for constitutive and induced trichome density also allowed for the investigation of damage response (plasticity) QTL. Both novel and previously identified QTL for constitutive trichome density and the first QTL for induced trichome density and response are identified. Interestingly, two of the four parental accessions and multiple RILs in each population exhibited lower trichome density following leaf damage, a response not previously described in A. thaliana. Importantly, a single QTL was mapped for the response phenotype and allelic variation at this locus appears to determine response trajectory in RILs. The data also show that epistatic interactions are a significant component of the genetic architecture of trichome density.

CONCLUSIONS

Together, our results provide further insights into the genetic architecture of constitutive trichome density and new insights into induced trichome density in A. thaliana specifically and to our understanding of the genetic underpinnings of natural variation generally.

Involvement of ASR genes in aluminium tolerance mechanisms in rice

Among cereal crops, rice is considered the most tolerant to aluminium (Al). However, variability among rice genotypes leads to remarkable differences in the degree of Al tolerance for distinct cultivars. A number of studies have demonstrated that rice plants achieve Al tolerance through an unknown mechanism that is independent of root tip Al exclusion. We have analysed expression changes of the rice ASR gene family as a function of Al treatment. The gene ASR5 was differentially regulated in the Al-tolerant rice ssp. Japonica cv. Nipponbare. However, ASR5 expression did not respond to Al exposure in Indica cv. Taim rice roots, which are highly Al sensitive. Transgenic plants carrying RNAi constructs that targeted the ASR genes were obtained, and increased Al susceptibility was observed in T1 plants. Embryogenic calli of transgenic rice carrying an ASR5-green fluorescent protein fusion revealed that ASR5 was localized in both the nucleus and cytoplasm. Using a proteomic approach to compare non-transformed and ASR-RNAi plants, a total of 41 proteins with contrasting expression patterns were identified. We suggest that the ASR5 protein acts as a transcription factor to regulate the expression of different genes that collectively protect rice cells from Al-induced stress responses.

Genetic basis of local adaptation and flowering time variation in Arabidopsis lyrata

Understanding how genetic variation at individual loci contributes to adaptation of populations to different local environments is an important topic in modern evolutionary biology. To date, most evidence has pointed to conditionally neutral quantitative trait loci (QTL) showing fitness effects only in some environments, while there has been less evidence for single-locus fitness trade-offs. At QTL underlying local adaptation, alleles from the local population are expected to show a fitness advantage. Cytoplasmic genomes also can have a role in local adaptation, but the role of cytonuclear interactions in adaptive differentiation has remained largely unknown. We mapped genomic regions underlying adaptive differentiation in multiple fitness components and flowering time in diverged populations of a perennial plant Arabidopsis lyrata. Experimental hybrids for this purpose were grown in natural field conditions of the parental populations in Norway and North Carolina (NC), USA, and in the greenhouse. We found QTL where high fitness and early flowering were associated with local alleles, indicating a role of different selection pressures in phenotypic differentiation. At two QTL regions, a fitness component showing local adaptation between the parental populations also showed signs of putative fitness trade-offs. Beneficial dominance effects of conditionally neutral QTL for different fitness components resulted in hybrid vigour at the Norwegian site in the F(2) hybrids. We also found that cytoplasmic genomes contributed to local adaptation and hybrid vigour by interacting with nuclear QTL, but these interactions did not show evidence for cytonuclear coadaptation (high fitness of local alleles combined with the local cytoplasm).

Natural variation in GL1 and its effects on trichome density in Arabidopsis thaliana

The ultimate understanding of how biological diversity arises, is maintained, and lost depends on identifying the genes responsible. Although a good deal has been discovered about gene function over the past few decades, far less is understood about gene effects, that is, how natural variation in a gene contributes to natural variation in phenotypes. Trichome density in Arabidopsis thaliana is an ideal trait for studies of natural molecular and phenotypic variation, as trichome initiation is genetically well-characterized and trichome density is highly variable in and among natural populations. Here, we show that variation at GLABRA1 (GL1), an R2R3 MYB transcription factor gene, which has a role in trichome initiation, has qualitative and likely quantitative effects on trichome density in natural accessions of A. thaliana. Specifically, we characterize four independent loss-of-function alleles for GL1, each of which yields a glabrous phenotype. Further, we find that a pattern of common polymorphisms confined to the GL1 locus is associated with quantitative variation for trichome density. While mutations resulting in a glabrous phenotype are primarily coding changes, the pattern resulting in quantitative variation spans both coding and regulatory regions. These data show that GL1 is an important source of trichome density variation within A. thaliana and, along with recent reports, suggest that the TTG1 epidermal cell fate pathway generally may be the key molecular genetic source of natural trichome density variation and an important model for the study of molecular evolution.

Additive and non-additive effects of simulated leaf and inflorescence damage on survival, growth and reproduction of the perennial herb Arabidopsis lyrata

Herbivores may damage both leaves and reproductive structures, and although such combined damage may affect plant fitness non-additively, this has received little attention. We conducted a 2-year field experiment with a factorial design to examine the effects of simulated leaf (0, 12.5, 25, or 50% of leaf area removed) and inflorescence damage (0 vs. 50% of inflorescences removed) on survival, growth and reproduction in the perennial herb Arabidopsis lyrata. Leaf and inflorescence damage negatively and independently reduced flower, fruit and seed production in the year of damage; leaf damage also reduced rosette size by the end of the first season and flower production in the second year. Leaf damage alone reduced the proportion of flowers forming a fruit and fruit production per plant the second year, but when combined with inflorescence damage no such effect was observed (significant leaf × inflorescence damage interaction). Damage to leaves (sources) caused a greater reduction in future reproduction than did simultaneous damage to leaves and inflorescences (sinks). This demonstrates that a full understanding of the effects of herbivore damage on plant fitness requires that consequences of damage to vegetative and reproductive structures are evaluated over more than 1 year and that non-additive effects are considered.

Extremely low genetic variability and highly structured local populations of Arabidopsis thaliana at higher latitudes

The genetic diversity and population structure of Arabidopsis thaliana populations from Norway were studied and compared to a worldwide sample of A. thaliana to investigate the demographic history and elucidate possible colonization routes of populations at the northernmost species limit. We genotyped 282 individuals from 31 local populations using 149 single nucleotide polymorphism markers. A high level of population subdivision (F(ST) = 0.85 ± 0.007) was found indicating that A. thaliana is highly structured at the regional level. Significant relationships between genetic and geographical distances were found, suggesting an isolation by distance mode of evolution. Genetic diversity was much lower, and the level of linkage disequilibrium was higher in populations from the north (65-68°N) compared to populations from the south (59-62°N); this is consistent with a northward expansion pattern. A neighbour-joining tree showed that populations from northern Norway form a separate cluster, while the remaining populations are distributed over a few minor clusters. Minimal gene flow seems to have occurred between populations in different regions, especially between the geographically distant northern and southern populations. Our data suggest that northern populations represent a homogenous group that may have been established from a few founders during northward expansions, while populations in the central part of Norway constitute an admixed group established by founders of different origins, most probably as a result of human-mediated gene flow. Moreover, Norwegian populations appeared to be homogenous and isolated compared to a worldwide sample of A. thaliana, but they are still grouped with Swedish populations, which may indicate common colonization histories.

Heterogeneous selection on trichome production in Alaskan Arabidopsis kamchatica (Brassicaceae)

PREMISE OF THE STUDY

Environmental heterogeneity is thought to be one of the primary factors in the evolutionary maintenance of morphological variation. Here, we explore the role of environmental heterogeneity in the maintenance of variation in leaf hair (trichome) production in Arabidopsis kamchatica. •

METHODS

We investigate abiotic correlates of trichome production in A. kamchatica via surveys of both herbarium specimens and wild populations. In addition, we examine patterns of phenotypic selection on trichome production among populations that differ in environmental characteristics. •

KEY RESULTS

Trichome-producing herbarium specimens were more likely to occur at lower latitudes and in locations with lower mean annual precipitation and less annual variation in temperature than glabrous specimens. In surveys of wild populations, frequencies of trichome-producing plants were higher in drier habitats than in wetter environments. Using phenotypic selection analysis, we found divergent selection through female fitness (fruit production) on trichome number in populations that differ in environmental characteristics; there was selection for reduced trichome number in one population and selection for increased trichome number in another population. In a population containing both glabrous and trichome-producing plants, glabrous plants produced significantly more fruits than trichome-producing individuals, which indicates selection against the trichome morph. •

CONCLUSIONS

Our results demonstrate that there is heterogeneity in selection among populations, which could be responsible for the maintenance of trichome variation in Alaskan populations of A. kamchatica.

A population genomic approach to map recent positive selection in model species

Based on nearly complete genome sequences from a variety of organisms data on naturally occurring genetic variation on the scale of hundreds of loci to entire genomes have been collected in recent years. In parallel, new statistical tests have been developed to infer evidence of recent positive selection from these data and to localize the target regions of selection in the genome. These methods have now been successfully applied to Drosophila melanogaster, humans, mice and a few plant species. In genomic regions of normal recombination rates, the targets of positive selection have been mapped down to the level of individual genes.

Patterns of polymorphism and demographic history in natural populations of Arabidopsis lyrata

BACKGROUND

Many of the processes affecting genetic diversity act on local populations. However, studies of plant nucleotide diversity have largely ignored local sampling, making it difficult to infer the demographic history of populations and to assess the importance of local adaptation. Arabidopsis lyrata, a self-incompatible, perennial species with a circumpolar distribution, is an excellent model system in which to study the roles of demographic history and local adaptation in patterning genetic variation.

PRINCIPAL FINDINGS

We studied nucleotide diversity in six natural populations of Arabidopsis lyrata, using 77 loci sampled from 140 chromosomes. The six populations were highly differentiated, with a median FST of 0.52, and structure analysis revealed no evidence of admixed individuals. Average within-population diversity varied among populations, with the highest diversity found in a German population; this population harbors 3-fold higher levels of silent diversity than worldwide samples of A. thaliana. All A. lyrata populations also yielded positive values of Tajima's D. We estimated a demographic model for these populations, finding evidence of population divergence over the past 19,000 to 47,000 years involving non-equilibrium demographic events that reduced the effective size of most populations. Finally, we used the inferred demographic model to perform an initial test for local adaptation and identified several genes, including the flowering time gene FCA and a disease resistance locus, as candidates for local adaptation events.

CONCLUSIONS

Our results underscore the importance of population-specific, non-equilibrium demographic processes in patterning diversity within A. lyrata. Moreover, our extensive dataset provides an important resource for future molecular population genetic studies of local adaptation in A. lyrata.

Genetic divergence of climatically marginal populations of Vicia pisiformis on the Scandinavian Peninsula

Vicia pisiformis L. is a perennial leguminous plant with a main distribution in broadleaved forest-steppes of eastern Europe. The species is classified as endangered (EN) according to the IUCN red-lists in both Norway and Sweden, due to severe fragmentation, small population sizes and continuing population decline. The populations on the Scandinavian Peninsula constitute the northern limit of the species distribution and are mostly restricted to warm stony slopes with predominantly southern aspects. In this study we used the AFLP method, which is a high-resolution genetic fingerprint method. Samples were collected from 22 Scandinavian populations. The overall genetic structure was analysed in an AMOVA, in a Mantel test and through constrained correspondence analysis (CCA). The ordination scores representing non-geographic genetic divergence were extracted from the CCA and analysed in a linear model using habitat variables and population size as explanatory variables. We found (i) a strong geographic structure, (ii) significant genetic divergence between populations, (iii) that this genetic divergence remained significant even after removing the effect of geography in a partial CCA and (iv) that the remaining non-geographic part of genetic divergence (distance from the ordination centre) was associated with aspect, populations with a northern aspect were more genetically divergent. Aspect explains more variation than population size and is the only variable retained in the minimal adequate model. We suggest that local adaptation has caused this divergence from an expected geographical pattern of genetic variation. This explanation is further supported by the association between aspect and specific AFLP fragments.

Genome-wide effects of postglacial colonization in Arabidopsis lyrata

The perennial outcrossing Arabidopsis lyrata is becoming a plant model species for molecular ecology and evolution. However, its evolutionary history, and especially the impact of the climatic oscillations of the Pleistocene on its genetic diversity and population structure, is not well known. We analyzed the broad-scale population structure of the species based on microsatellite variation at 22 loci. A wide sample in Europe revealed that glaciations and postglacial colonization have caused high divergence and high variation in variability between populations. Colonization from Central Europe to Iceland and Scandinavia was associated with a strong decrease of genetic diversity from South to North. On the other hand, the Russian population included in our data set may originate from a different refugium probably located more to the East. These genome-wide patterns must be taken into account in studies aiming at elucidating the genetic basis of local adaptation. As shown by sequence data, most of the loci used in this study do not evolve like typical microsatellite loci and show variable levels of homoplasy: this mode of evolution makes these markers less suitable to investigate the between-continent divergence and more generally the worldwide evolution of the species. Finally, a strong negative correlation was detected between levels of within-population diversity and indices of differentiation such as F(ST). We discuss the causes of this correlation as well as the potential bias it induces on the quantification and interpretation of population structure.

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.

Glucosinolate and trichome defenses in a natural Arabidopsis lyrata population

Glucosinolates (GS) and trichomes contribute to plant resistance against insect herbivores in the model Arabidopsis thaliana. The functional and genetic characteristics of herbivore defense, however, can differ even between closely related species. In a quantitative genetic experiment with the out-crossing perennial Arabidopsis lyrata spp. petraea, we measured constitutive GS composition, trichome density, leaf thickness, and plant resistance in four different herbivore interactions. In a single population of A. lyrata, we found heritable variation for trichome density as well as GS amount and carbon side-chain elongation ratios associated with activity in methylthioalkylmalate synthase (MAM). Unexpectedly, heritabilities for indole GS in A. lyrata were high and less affected by differences in plant age and environment than aliphatic GS. We found significant heritability in plant resistance to the specialist Plutella xylostella and generalist Trichoplusia ni, but not to the specialists Pieris brassicae and Phyllotreta cruciferae. Analyses of phenotypic and genetic correlations between candidate defense traits and insect resistance suggested that A. lyrata resistance was conferred by a combination of indole GS amount and trichome density, and, to a lesser extent, aliphatic GS ratios and leaf thickness. Variation in the most abundant compound, the aliphatic 3-hydroxypropyl GS, had little impact on A. lyrata herbivore resistance. The contribution of defense traits to resistance depended on the experimental herbivory context, and resistances were weakly correlated. A diversified defense strategy is likely to be important for long-lived individuals of A. lyrata that are subject to attack by many different herbivores in nature.

Selection on flowering time and floral display in an alpine and a lowland population of Arabidopsis lyrata

To determine whether population differentiation in flowering time is consistent with differences in current selection, we quantified phenotypic selection acting through female reproductive success on flowering phenology and floral display in two Scandinavian populations of the outcrossing, perennial herb Arabidopsis lyrata in two years. One population was located in an alpine environment strongly affected by grazing, whereas the other was close to sea level and only moderately affected by herbivory. Multiple regression models indicated directional selection for early end of flowering in one year in the lowland population, and directional selection for early start of flowering in one year in the alpine population. As expected, there was selection for more inflorescences in the lowland population. However, in the alpine population, plants with many inflorescences were selectively grazed and the number of inflorescences produced was negatively related to female fitness in one year and not significantly related to female fitness in the second year. The results are consistent with the hypothesis that genetic differentiation in flowering phenology between the study populations is adaptive, and indicate that interactions with selective grazers may strongly influence selection on floral display in A. lyrata.

Evolution and genetic differentiation among relatives of Arabidopsis thaliana

Arabidopsis thaliana is one of the most intensively studied plant species. More recently, information is accumulating about its closest relatives, the former genus Cardaminopsis. A. thaliana diverged from these relatives, actually treated within three major lineages (Arabidopsis lyrata, Arabidopsis halleri, and Arabidopsis arenosa), approximately 5 mya. Significant karyotype evolution in A. thaliana with base chromosome number reduction from x=8 to x=5 might indicate and favor effective genetic isolation from these other species, although hybrids are occurring naturally and have been also constituted under controlled conditions. We tested the evolutionary significance to separate the x=5 from the x=8 lineage using DNA sequence data from the plastome and the nuclear ribosomal DNA based on an extensive, representative worldwide sampling of nearly all taxonomic entities. We conclude that (i) A. thaliana is clearly separated phylogenetically from the x=8 lineage, (ii) five major lineages outside A. thaliana can be identified (A. lyrata, A. arenosa, A. halleri, Arabidopsis croatica, and Arabidopsis pedemontana) together with Arabidopsis cebennensis, and (iii) centers of genetic and morphological diversity are mostly in congruence and are located close to the Balkans in Austria and Slovakia outside glaciated and permafrost regions with few notable exceptions.

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.

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