A hot topic: the genetics of adaptation to geothermal vents in Mimulus guttatus

Identifying the individual loci and mutations that underlie adaptation to extreme environments has long been a goal of evolutionary biology. However, finding the genes that underlie adaptive traits is difficult for several reasons. First, because many traits and genes evolve simultaneously as populations diverge, it is difficult to disentangle adaptation from neutral demographic processes. Second, finding the individual loci involved in any trait is challenging given the respective limitations of quantitative and population genetic methods. In this issue of Molecular Ecology, Hendrick et al. (2016) overcome these difficulties and determine the genetic basis of microgeographic adaptation between geothermal vent and nonthermal populations of Mimulus guttatus in Yellowstone National Park. The authors accomplish this by combining population and quantitative genetic techniques, a powerful, but labour‐intensive, strategy for identifying individual causative adaptive loci that few studies have used (Stinchcombe & Hoekstra 2008 ). In a previous common garden experiment (Lekberg et al. 2012), thermal M. guttatus populations were found to differ from their closely related nonthermal neighbours in various adaptive phenotypes including trichome density. Hendrick et al. (2016) combine quantitative trait loci (QTL) mapping, population genomic scans for selection and admixture mapping to identify a single genetic locus underlying differences in trichome density between thermal and nonthermal M. guttatus. The candidate gene, R2R3 MYB, is homologous to genes involved in trichome development across flowering plants. The major trichome QTL, Tr14, is also involved in trichome density differences in an independent M. guttatus population comparison (Holeski et al. 2010) making this an example of parallel genetic evolution.     

The candidate gene approach in plant genetics: a review

The candidate gene (CG) approach has been applied in plant genetics in the past decade for the characterisation and cloning of Mendelian and quantitative trait loci (QTLs). It constitutes a complementary strategy to map-based cloning and insertional mutagenesis. The goal of this paper is to present an overview of CG analyses in plant genetics. CG analysis is based on the hypothesis that known-function genes (the candidate genes) could correspond to loci controlling traits of interest. CGs refer either to cloned genes presumed to affect a given trait (`functional CGs') or to genes suggested by their close proximity on linkage maps to loci controlling the trait (`positional CGs'). In plant genetics, the most common way to identify a CG is to look for map co-segregation between CGs and loci affecting the trait. Statistical association analyses between molecular polymorphisms of the CG and variation in the trait of interest have also been carried out in a few studies. The final validation of a CG will be provided through physiological analyses, genetic transformation and/or sexual complementation. Theoretical and practical applications of validated CGs in plant genetics and breeding are discussed.     

The genomics of adaptation,divergence and speciation: a congealing theory

In this issue, Flaxman et al. ( 2014 ) report the results of sophisticated whole‐genome simulations of speciation with gene flow, enhancing our understanding of the process by building on previous single‐locus, multilocus and analytical works. Their findings provide us with new insights about how genomes can diverge and the importance of statistical and chromosomal linkage in facilitating reproductive isolation. The authors characterize the conditions under which, even with high gene flow and weak divergent selection, reproductive isolation between populations can occur due to the emergent stochastic process of genomewide congealing, where numerous statistically or physically linked loci of small effect allow selection to limit effective migration rates. The initial congealing event can occur within a broad range conditions, and once initiated, the self‐reinforcing process leads to rapid divergence and ultimately two reproductively isolated populations. Flaxman et al.'s ( 2014 ) work is a valuable contribution to our understanding of speciation with gene flow and in making a more predictive field of evolutionary genomics and speciation.     

In the presence of population structure: From genomics to candidate genes underlying local adaptation

Understanding the genomic signatures, genes, and traits underlying local adaptation of organisms to heterogeneous environments is of central importance to the field evolutionary biology. To identify loci underlying local adaptation, models that combine allelic and environmental variation while controlling for the effects of population structure have emerged as the method of choice. Despite being evaluated in simulation studies, there has not been a thorough investigation of empirical evidence supporting local adaptation across these alleles. To evaluate these methods, we use 875 Arabidopsis thaliana Eurasian accessions and two mixed models (GEMMA and LFMM) to identify candidate SNPs underlying local adaptation to climate. Subsequently, to assess evidence of local adaptation and function among significant SNPs, we examine allele frequency differentiation and recent selection across Eurasian populations, in addition to their distribution along quantitative trait loci (QTL) explaining fitness variation between Italy and Sweden populations and cis‐regulatory/nonsynonymous sites showing significant selective constraint. Our results indicate that significant LFMM/GEMMA SNPs show low allele frequency differentiation and linkage disequilibrium across locally adapted Italy and Sweden populations, in addition to a poor association with fitness QTL peaks (highest logarithm of odds score). Furthermore, when examining derived allele frequencies across the Eurasian range, we find that these SNPs are enriched in low‐frequency variants that show very large climatic differentiation but low levels of linkage disequilibrium. These results suggest that their enrichment along putative functional sites most likely represents deleterious variation that is independent of local adaptation. Among all the genomic signatures examined, only SNPs showing high absolute allele frequency differentiation (AFD) and linkage disequilibrium (LD) between Italy and Sweden populations showed a strong association with fitness QTL peaks and were enriched along selectively constrained cis‐regulatory/nonsynonymous sites. Using these SNPs, we find strong evidence linking flowering time, freezing tolerance, and the abscisic‐acid pathway to local adaptation.     

Dissecting the role of a large chromosomal inversion in life history divergence throughout the Mimulus guttatus species complex

Chromosomal inversions can play an important role in adaptation, but the mechanism of their action in many natural populations remains unclear. An inversion could suppress recombination between locally beneficial alleles, thereby preventing maladaptive reshuffling with less‐fit, migrant alleles. The recombination suppression hypothesis has gained much theoretical support but empirical tests are lacking. Here, we evaluated the evolutionary history and phenotypic effects of a chromosomal inversion which differentiates annual and perennial forms of Mimulus guttatus. We found that perennials likely possess the derived orientation of the inversion. In addition, this perennial orientation occurs in a second perennial species, M. decorus, where it is strongly associated with life history differences between co‐occurring M. decorus and annual M. guttatus. One prediction of the recombination suppression hypothesis is that loci contributing to local adaptation will predate the inversion. To test whether the loci influencing perenniality pre‐date this inversion, we mapped QTLs for life history traits that differ between annual M. guttatus and a more distantly related, collinear perennial species, M. tilingii. Consistent with the recombination suppression hypothesis, we found that this region is associated with life history in the absence of the inversion, and this association can be broken into at least two QTLs. However, the absolute phenotypic effect of the LG8 inversion region on life history is weaker in M. tilingii than in perennials which possess the inversion. Thus, while we find support for the recombination suppression hypothesis, the contribution of this inversion to life history divergence in this group is likely complex.     

Genetic basis for rapidly evolved tolerance in the wild: adaptation to toxic pollutants by an estuarine fish species

Atlantic killifish (Fundulus heteroclitus) residing in some urban and industrialized estuaries of the US eastern seaboard demonstrate recently evolved and extreme tolerance to toxic aryl hydrocarbon pollutants, characterized as dioxin‐like compounds (DLCs). Here, we provide an unusually comprehensive accounting (69%) through quantitative trait locus (QTL) analysis of the genetic basis for DLC tolerance in killifish inhabiting an urban estuary contaminated with PCB congeners, the most toxic of which are DLCs. Consistent with mechanistic knowledge of DLC toxicity in fish and other vertebrates, the aryl hydrocarbon receptor (ahr2) region accounts for 17% of trait variation; however, QTL on independent linkage groups and their interactions have even greater explanatory power (44%). QTL interpreted within the context of recently available Fundulus genomic resources and shared synteny among fish species suggest adaptation via interacting components of a complex stress response network. Some QTL were also enriched in other killifish populations characterized as DLC‐tolerant and residing in distant urban estuaries contaminated with unique mixtures of pollutants. Together, our results suggest that DLC tolerance in killifish represents an emerging example of parallel contemporary evolution that has been driven by intense human‐mediated selection on natural populations.     

Adaptation genomics: next generation sequencing reveals a shared haplotype for rapid early development in geographically and genetically distant populations of rainbow trout

Local adaptation occurs when a population evolves a phenotype that confers a selective advantage in its local environment, but which may not be advantageous in other habitats. Restricted gene flow and strong selection pressures are prerequisites for local adaptation. Fishes in the family Salmonidae are predicted to provide numerous examples of local adaptation because of the high fidelity of returning to spawn in their natal streams, which results in highly structured populations, and the wide diversity of environments that salmonids have colonized. These conditions are ideally suited for producing a set of specialist phenotypes, whose fitness is maximized for one specific habitat, rather than a generalist phenotype similarly viable in several environments. Understanding patterns and processes leading to local adaptations has long been a goal of evolutionary biology, but it is only recently that identifying the molecular basis for local adaptation has become feasible because of advances in genomic technologies. The study of shared adaptive phenotypes in populations that are both geographically distant and genetically distinct should reveal some of the fundamental molecular mechanisms associated with local adaptation. In this issue of Molecular Ecology, Miller et al. (2012) make a significant contribution to the development of adaptation genomics. This study suggests that salmonids use standing genetic variation to select beneficial alleles for local adaptations rather than de novo mutations in the same gene or alternative physiological pathways. Identifying the genetic basis for local adaptation has major implications for the management, conservation and potential restoration of salmonid populations.     

Local adaptation of sex induction in a facultative sexual crustacean: insights from QTL mapping and natural populations of Daphnia magna

Dormancy is a common adaptation in invertebrates to survive harsh conditions. Triggered by environmental cues, populations produce resting eggs that allow them to survive temporally unsuitable conditions. Daphnia magna is a crustacean that reproduces by cyclical parthenogenesis, alternating between the production of asexual offspring and the sexual reproduction of diapausing eggs (ephippia). Prior to ephippia production, males (necessary to ensure ephippia fertilization) are produced parthenogenetically. Both the production of ephippia and the parthenogenetic production of males are induced by environmental factors. Here, we test the hypothesis that the induction of D. magna resting egg production shows a signature of local adaptation. We postulated that Daphnia from permanent ponds would produce fewer ephippia and males than Daphnia from intermittent ponds and that the frequency and season of habitat deterioration would correlate with the timing and amount of male and ephippia production. To test this, we quantified the production of males and ephippia in clonal D. magna populations in several different controlled environments. We found that the production of both ephippia and males varies strongly among populations in a way that suggests local adaptation. By performing quantitative trait locus mapping with parent clones from contrasting pond environments, we identified nonoverlapping genomic regions associated with male and ephippia production. As the traits are influenced by two different genomic regions, and both are necessary for successful resting egg production, we suggest that the genes for their induction co‐evolve.     

Genome-wide association mapping identifies the genetic basis of discrete and quantitative variation in sexual weaponry in a wild sheep population

Understanding the genetic architecture of phenotypic variation in natural populations is a fundamental goal of evolutionary genetics. Wild Soay sheep (Ovis aries) have an inherited polymorphism for horn morphology in both sexes, controlled by a single autosomal locus, Horns. The majority of males have large normal horns, but a small number have vestigial, deformed horns, known as scurs; females have either normal horns, scurs or no horns (polled). Given that scurred males and polled females have reduced fitness within each sex, it is counterintuitive that the polymorphism persists within the population. Therefore, identifying the genetic basis of horn type will provide a vital foundation for understanding why the different morphs are maintained in the face of natural selection. We conducted a genome-wide association study using ～36000 single nucleotide polymorphisms (SNPs) and determined the main candidate for Horns as RXFP2, an autosomal gene with a known involvement in determining primary sex characters in humans and mice. Evidence from additional SNPs in and around RXFP2 supports a new model of horn-type inheritance in Soay sheep, and for the first time, sheep with the same horn phenotype but different underlying genotypes can be identified. In addition, RXFP2 was shown to be an additive quantitative trait locus (QTL) for horn size in normal-horned males, accounting for up to 76% of additive genetic variation in this trait. This finding contrasts markedly from genome-wide association studies of quantitative traits in humans and some model species, where it is often observed that mapped loci only explain a modest proportion of the overall genetic variation.     

Genomic signatures of host‐associated divergence and adaptation in a coral‐eating snail,Coralliophila violacea (Kiener, 1836)

The fluid nature of the ocean, combined with planktonic dispersal of marine larvae, lowers physical barriers to gene flow. However, divergence can still occur despite gene flow if strong selection acts on populations occupying different ecological niches. Here, we examined the population genomics of an ectoparasitic snail, Coralliophila violacea (Kiener 1836), that specializes on Porites corals in the Indo‐Pacific. Previous genetic analyses revealed two sympatric lineages associated with different coral hosts. In this study, we examined the mechanisms promoting and maintaining the snails’ adaptation to their coral hosts. Genome‐wide single nucleotide polymorphism (SNP) data from type II restriction site‐associated DNA (2b‐RAD) sequencing revealed two differentiated clusters of C. violacea that were largely concordant with coral host, consistent with previous genetic results. However, the presence of some admixed genotypes indicates gene flow from one lineage to the other. Combined, these results suggest that differentiation between host‐associated lineages of C. violacea is occurring in the face of ongoing gene flow, requiring strong selection. Indeed, 2.7% of all SNP loci were outlier loci (73/2,718), indicative of divergence with gene flow, driven by adaptation of each C. violacea lineage to their specific coral hosts.     

Landscape genetics and hierarchical genetic structure in Atlantic salmon: the interaction of gene flow and local adaptation

Disentangling evolutionary forces that may interact to determine the patterns of genetic differentiation within and among wild populations is a major challenge in evolutionary biology. The objective of this study was to assess the genetic structure and the potential influence of several ecological variables on the extent of genetic differentiation at multiple spatial scales in a widely distributed species, the Atlantic salmon, Salmo salar . A total of 2775 anadromous fish were sampled from 51 rivers along the North American Atlantic coast and were genotyped using 13 microsatellites. A Bayesian analysis clustered these populations into seven genetically and geographically distinct groups, characterized by different environmental and ecological factors, mainly temperature. These groups were also characterized by different extent of genetic differentiation among populations. Dispersal was relatively high and of the same magnitude within compared to among regional groups, which contrasted with the maintenance of a regional genetic structure. However, genetic differentiation was lower among populations exchanging similar rates of local as opposed to inter-regional migrants, over the same geographical scale. This raised the hypothesis that gene flow could be constrained by local adaptation at the regional scale. Both coastal distance and temperature regime were found to influence the observed genetic structure according to landscape genetic analyses. The influence of other factors such as latitude, river length and altitude, migration tactic, and stocking was not significant at any spatial scale. Overall, these results suggested that the interaction between gene flow and thermal regime adaptation mainly explained the hierarchical genetic structure observed among Atlantic salmon populations.     

Toward the evolutionary genomics of gametophytic divergence: patterns of transmission ratio distortion in monkeyflower (Mimulus) hybrids reveal a complex genetic basis for conspecific pollen precedence

Conspecific pollen precedence (CPP) is a major component of reproductive isolation between many flowering plant taxa and may reveal mechanisms of gametophytic evolution within species, but little is known about the genetic basis and evolutionary history of CPP. We systematically investigated the genetic architecture of CPP using patterns of transmission ratio distortion (TRD) in F2 and backcross hybrids between closely related species of Mimulus (Phrymaceae) with divergent mating systems. We found that CPP in Mimulus hybrids was polygenic and was the majority source of interspecific TRD genome-wide, with at least eight genomic regions contributing to the transmission advantage of M. guttatus pollen grains on M. guttatus styles. In aggregate, these male-specific transmission ratio distorting loci (TRDLs) were more than sufficient to account for the 100% precedence of pure M. guttatus pollen over M. nasutus pollen in mixed pollinations of M. guttatus. All but one of these pollen TRDLs were style-dependent; that is, we observed pollen TRD in F(1) and/or M. guttatus styles, but not in M. nasutus styles. These findings suggest that species-specific differences in pollen tube performance accumulate gradually and may have been driven by coevolution between pollen and style in the predominantly outcrossing M. guttatus.     

Ready on arrival: Standing variation at a chromosomal inversion contributes to rapid adaptation in an invasive marine crab

In this issue of Molecular Ecology, Tepolt et al. (2021) illustrate how the genetic architecture of adaptation and life history influence invasive success. A marvel of many invasive species is that they are incredibly successful despite evolutionary expectations that they will have low adaptive potential and suffer inbreeding depression due to initially small founding population sizes. Determining the combinations of ecoevolutionary factors that permit this apparent “genetic paradox of invasions” is an ongoing endeavour of invasive species research. Tepolt et al. (2021) study the European green crab in its invasive range on the North American west coast. Following a single introduction into California, this crab quickly spread across a wide latitude gradient, despite low diversity in the original founding population. Adaptation of this crab to clinal variation in temperature appeared largely driven by an inferred chromosomal inversion. This inversion exists as a balanced polymorphism in the European home range of green crabs and is associated with thermal tolerance. Tepolt et al. (2021) therefore demonstrate that adaptive evolution post introduction need not be impeded by bottlenecks if variation at key parts of the genome is available and can be maintained in introduced populations. Moreover, Tepolt et al. (2021) show how chromosomal inversions acting as large-effect loci might facilitate rapid responses to selection in introduced populations.     

Bypassing of a polygenic Microcyclus ulei resistance in rubber tree, analyzed by QTL detection

Genetic resistance components of the Hevea brasiliensis x H. benthamiana RO 38 cultivar to Microcyclus ulei disease were investigated by inoculating isolates which succeeded in partially or completely infecting genotypes of a mapping population. Progeny of a cross between RO 38 and a susceptible cultivar was inoculated under controlled conditions with three isolates and scored for two resistance traits. Interval mapping and a nonparametric test were used to detect resistance quantitative trait loci (QTLs). Eight significant QTLs were detected, all of them inherited from the interspecific parent. Among these QTLs, only one contributed to the partial resistance against a highly pathogenic isolate, and no QTL was detected for resistance against the most pathogenic isolate. As an unexpected result, a single isolate can thus completely bypass this polygenic resistance. This complex situation, where no clear relationship can be established between number of resistance factors and qualitative vs quantitative or partial vs complete resistance, is discussed and suggestions as to the detection of new and sustainable resistance sources are proposed.     

Interval mapping of genes for quantitative resistance of maize to Setosphaeria turcica,cause of northern leaf blight,in a tropical environment

Quantitative trait loci (QTL) involved in the resistance of maize to Setosphaeria turcica, the causal agent of northern leaf blight, were located by interval mapping analysis of 121 F_2:3 lines derived from a cross between Mo17 (moderately resistant) and B52 (susceptible). A linkage map spanning 112 RFLP loci with 15 cM mean interval length was constructed, based on marker data recorded in a previous study. Field tests with artificial inoculation were conducted at three sites in tropical mid- to high-altitude regions of Kenya, East Africa. Host-plant response was measured in terms of incubation period, disease severity (five scoring dates), and the area under the disease progress curve (AUDPC). Heritability of all traits was high (around 0.75). QTL associated with the incubation period were located on chromosomes 2S and 8L. For disease severity and AUDPC, significant QTL were detected in the putative centromeric region of chromosome 1 and on 2S, 3L, 5S, 6L, 7L, 8L and 9S. On 2S the same marker interval which carried a gene enhancing latent period was also associated with reduced disease severity of juvenile plants. QTL on chromosomes 3L, 5S, 7L and 8L were significant across environments but all other QTL were affected by a large genotype x environment interaction. Partially dominant gene action for resistance as well as for susceptibility was prevailing. Single QTL explained 10 to 38% of the phenotypic variation of the traits. All but the QTL on chromosomes 1, 6 and 9 were contributed by the resistant parent Mo17. On chromosome 8L a QTL mapped to the same region as the major race-specific gene Ht2, supporting the hypothesis that some qualitative and quantitative resistance genes may be allelic.Abbreviations AUDPC area under the disease progress curve - CIMMYT International Maize and Wheat Improvement Center - KARI Kenya Agricultural Research Institute - NCLB northern corn leaf blight - QTL quantitative trait locus/loci     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.     

The genetic basis of floral traits associated with mating system evolution in Leptosiphon (Polemoniaceae): an analysis of quantitative trait loci

Mapping of quantitative trait loci (QTL) was used to investigate the genetic architecture of divergence in floral characters associated with the mating system, an important adaptive trait in angiosperms. Two species of Leptosiphon (Polemoniaceae), one strongly self-fertilizing (L. bicolor) and the other partially outcrossing (L. jepsonii), were crossed to produce F2 and both backcross progenies. For each crossing population, a linkage map was created using amplified fragment length polymorphism markers, and QTL were identified for several dimensions of floral size. For each of the five traits examined, three to seven QTL were detected, with independent datasets yielding congruent results in some but not all cases. The phenotypic effect of individual QTL was generally moderate. We estimated that many of the QTL were additive or showed dominance toward L. bicolor, whereas comparison of mean trait values for parental and cross progenies showed apparent overall dominance of L. jepsonii traits. Colocalization of QTL for different dimensions of floral size was consistent with high phenotypic correlations between floral traits. Substantial segregation distortion was observed in marker loci, the majority favoring alleles from the large-flowered parent. A low frequency of male sterility in the F2 population is consistent with the Dobzhansky-Muller model for the evolution of reproductive isolation.     

The selfing syndrome: a model for studying the genetic and evolutionary basis of morphological adaptation in plants

Background

In angiosperm evolution, autogamously selfing lineages have been derived from outbreeding ancestors multiple times, and this transition is regarded as one of the most common evolutionary tendencies in flowering plants. In most cases, it is accompanied by a characteristic set of morphological and functional changes to the flowers, together termed the selfing syndrome. Two major areas that have changed during evolution of the selfing syndrome are sex allocation to male vs. female function and flower morphology, in particular flower (mainly petal) size and the distance between anthers and stigma.

Scope

A rich body of theoretical, taxonomic, ecological and genetic studies have addressed the evolutionary modification of these two trait complexes during or after the transition to selfing. Here, we review our current knowledge about the genetics and evolution of the selfing syndrome.

Conclusions

We argue that because of its frequent parallel evolution, the selfing syndrome represents an ideal model for addressing basic questions about morphological evolution and adaptation in flowering plants, but that realizing this potential will require the molecular identification of more of the causal genes underlying relevant trait variation.     

The genetic basis of adaptive population differentiation: a quantitative trait locus analysis of fitness traits in two wild barley populations from contrasting habitats

We used a quantitative trait locus (QTL) approach to study the genetic basis of population differentiation in wild barley, Hordeum spontaneum. Several ecotypes are recognized in this model species, and population genetic studies and reciprocal transplant experiments have indicated the role of local adaptation in shaping population differences. We derived a mapping population from a cross between a coastal Mediterranean population and a steppe inland population from Israel and assessed F3 progeny fitness in the natural growing environments of the two parental populations. Dilution of the local gene pool, estimated as the proportion of native alleles at 96 marker loci in the recombinant lines, negatively affected fitness traits at both sites. QTLs for fitness traits tended to differ in the magnitude but not in the direction of their effects across sites, with beneficial alleles generally conferring a greater fitness advantage at their native site. Several QTLs showed fitness effects at one site only, but no opposite selection on individual QTLs was observed across the sites. In a common-garden experiment, we explored the hypothesis that the two populations have adapted to divergent nutrient availabilities. In the different nutrient environments of this experiment, but not under field conditions, fitness of the F3 progeny lines increased with the number of heterozygous marker loci. Comparison of QTL-effects that underlie genotype x nutrient interaction in the common-garden experiment and genotype x site interaction in the field suggested that population differentiation at the field sites may have been driven by divergent nutrient availabilities to a limited extent. Also in this experiment no QTLs were observed with opposite fitness effects in contrasting environments. Our data are consistent with the view that adaptive differentiation can be based on selection on multiple traits changing gradually along ecological gradients. This can occur without QTLs showing opposite fitness effects in the different environments, that is, in the absence of genetic trade-offs in performance between environments.