Christopher G. Oakley

Genetic mapping of adaptation reveals fitness tradeoffs in Arabidopsis thaliana

Significance Adaptation to local environmental conditions is common, but the genetic mechanisms of adaptation are poorly known. We produced recombinant inbred lines (RILs) of the model plant Arabidopsis thaliana by crossing populations that inhabit drastically different climates in Sweden and Italy, grew the RILs at the parental sites for 3 y, and genetically mapped quantitative trait loci (QTL) for fitness. The results demonstrate that surprisingly few QTL explain much of the adaptive divergence between the two plant populations. Moreover, we find strong evidence for tradeoffs (i.e., adaptation to one environment reduces performance elsewhere). The results shed light on processes governing the evolution of biological diversity and the potential for adaptive evolution in response to environmental change. Organisms inhabiting different environments are often locally adapted, and yet despite a considerable body of theory, the genetic basis of local adaptation is poorly understood. Unanswered questions include the number and effect sizes of adaptive loci, whether locally favored loci reduce fitness elsewhere (i.e., fitness tradeoffs), and whether a lack of genetic variation limits adaptation. To address these questions, we mapped quantitative trait loci (QTL) for total fitness in 398 recombinant inbred lines derived from a cross between locally adapted populations of the highly selfing plant Arabidopsis thaliana from Sweden and Italy and grown for 3 consecutive years at the parental sites (>40,000 plants monitored). We show that local adaptation is controlled by relatively few genomic regions of small to modest effect. A third of the 15 fitness QTL we detected showed evidence of tradeoffs, which contrasts with the minimal evidence for fitness tradeoffs found in previous studies. This difference may reflect the power of our multiyear study to distinguish conditionally neutral QTL from those that reflect fitness tradeoffs. In Sweden, but not in Italy, the local genotype underlying fitness QTL was often maladaptive, suggesting that adaptation there is constrained by a lack of adaptive genetic variation, attributable perhaps to genetic bottlenecks during postglacial colonization of Scandinavia or to recent changes in selection regime caused by climate change. Our results suggest that adaptation to markedly different environments can be achieved through changes in relatively few genomic regions, that fitness tradeoffs are common, and that lack of genetic variation can limit adaptation.

QTL mapping of freezing tolerance: links to fitness and adaptive trade-offs.

Local adaptation, defined as higher fitness of local vs. nonlocal genotypes, is commonly identified in reciprocal transplant experiments. Reciprocally adapted populations display fitness trade‐offs across environments, but little is known about the traits and genes underlying fitness trade‐offs in reciprocally adapted populations. We investigated the genetic basis and adaptive significance of freezing tolerance using locally adapted populations of Arabidopsis thaliana from Italy and Sweden. Previous reciprocal transplant studies of these populations indicated that subfreezing temperature is a major selective agent in Sweden. We used quantitative trait locus (QTL) mapping to identify the contribution of freezing tolerance to previously demonstrated local adaptation and genetic trade‐offs. First, we compared the genomic locations of freezing tolerance QTL to those for previously published QTL for survival in Sweden, and overall fitness in the field. Then, we estimated the contributions to survival and fitness across both field sites of genotypes at locally adaptive freezing tolerance QTL. In growth chamber studies, we found seven QTL for freezing tolerance, and the Swedish genotype increased freezing tolerance for five of these QTL. Three of these colocalized with locally adaptive survival QTL in Sweden and with trade‐off QTL for overall fitness. Two freezing tolerance QTL contribute to genetic trade‐offs across environments for both survival and overall fitness. A major regulator of freezing tolerance, CBF2, is implicated as a candidate gene for one of the trade‐off freezing tolerance QTL. Our study provides some of the first evidence of a trait and gene that mediate a fitness trade‐off in nature.

Flowering time QTL in natural populations of Arabidopsis thaliana and implications for their adaptive value

The genetic basis of phenotypic traits is of great interest to evolutionary biologists, but their contribution to adaptation in nature is often unknown. To determine the genetic architecture of flowering time in ecologically relevant conditions, we used a recombinant inbred line population created from two locally adapted populations of Arabidopsis thaliana from Sweden and Italy. Using these RILs, we identified flowering time QTL in growth chambers that mimicked the natural temperature and photoperiod variation across the growing season in each native environment. We also compared the genomic locations of flowering time QTL to those of fitness (total fruit number) QTL from a previous three‐year field study. Ten total flowering time QTL were found, and in all cases, the Italy genotype caused early flowering regardless of the conditions. Two QTL were consistent across chamber environments, and these had the largest effects on flowering time. Five of the fitness QTL colocalized with flowering time QTL found in the Italy conditions, and in each case, the local genotype was favoured. In contrast, just two flowering time QTL found in the Sweden conditions colocalized with fitness QTL and in only one case was the local genotype favoured. This implies that flowering time may be more important for adaptation in Italy than Sweden. Two candidate genes (FLC and VIN3) underlying the major flowering time QTL found in the current study are implicated in local adaptation.

Effects of population size and isolation on heterosis, mean fitness, and inbreeding depression in a perennial plant

• In small isolated populations, genetic drift is expected to increase chance fixation of partly recessive, mildly deleterious mutations, reducing mean fitness and inbreeding depression within populations and increasing heterosis in outcrosses between populations. • We estimated relative effective sizes and migration among populations and compared mean fitness, heterosis, and inbreeding depression for eight large and eight small populations of a perennial plant on the basis of fitness of progeny produced by hand pollinations within and between populations. • Migration was limited, and, consistent with expectations for drift, mean fitness was 68% lower in small populations; heterosis was significantly greater for small (mean = 70%, SE = 14) than for large populations (mean = 7%, SE = 27); and inbreeding depression was lower, although not significantly so, in small (mean = -0.29%, SE = 28) than in large (mean = 0.28%, SE = 23) populations. • Genetic drift promotes fixation of deleterious mutations in small populations, which could threaten their persistence. Limited migration will exacerbate drift, but data on migration and effective population sizes in natural populations are scarce. Theory incorporating realistic variation in population size and patterns of migration could better predict genetic threats to small population persistence.

Adaptive divergence in flowering time among natural populations of Arabidopsis thaliana: Estimates of selection and QTL mapping.

To identify the ecological and genetic mechanisms of local adaptation requires estimating selection on traits, identifying their genetic basis, and evaluating whether divergence in adaptive traits is due to conditional neutrality or genetic trade‐offs. To this end, we conducted field experiments for three years using recombinant inbred lines (RILs) derived from two ecotypes of Arabidopsis thaliana (Italy, Sweden), and at each parental site examined selection on flowering time and mapped quantitative trait loci (QTL). There was strong selection for early flowering in Italy, but weak selection in Sweden. Eleven distinct flowering time QTL were detected, and for each the Italian genotype caused earlier flowering. Twenty‐seven candidate genes were identified, two of which (FLC and VIN3) appear under major flowering time QTL in Italy. Seven of eight QTL in Italy with narrow credible intervals colocalized with previously reported fitness QTL, in comparison to three of four in Sweden. The results demonstrate that the magnitude of selection on flowering time differs strikingly between our study populations, that the genetic basis of flowering time variation is multigenic with some QTL of large effect, and suggest that divergence in flowering time between ecotypes is due mainly to conditional neutrality.

Factors influencing the effect size distribution of adaptive substitutions

The distribution of effect sizes of adaptive substitutions has been central to evolutionary biology since the modern synthesis. Early theory proposed that because large-effect mutations have negative pleiotropic consequences, only small-effect mutations contribute to adaptation. More recent theory suggested instead that large-effect mutations could be favoured when populations are far from their adaptive peak. Here we suggest that the distributions of effect sizes are expected to differ among study systems, reflecting the wide variation in evolutionary forces and ecological conditions experienced in nature. These include selection, mutation, genetic drift, gene flow, and other factors such as the degree of pleiotropy, the distance to the phenotypic optimum, whether the optimum is stable or moving, and whether new mutation or standing genetic variation provides the source of adaptive alleles. Our goal is to review how these factors might affect the distribution of effect sizes and to identify new research directions. Until more theory and empirical work is available, we feel that it is premature to make broad generalizations about the effect size distribution of adaptive substitutions important in nature.

POPULATION-LEVEL AND FAMILY-LEVEL INBREEDING DEPRESSION IN A CLEISTOGAMOUS PERENNIAL

Cleistogamous species constitute strong evidence for the stability of mixed selfing and outcrossing but have been largely overlooked as theoretical and empirical systems for studying the forces thought to influence mating‐system evolution. Inbreeding depression is expected to play a major role in the maintenance of outcrossing, yet reliable estimates of inbreeding depression for cleistogamous species are virtually nonexistent. We estimated both population‐ and family‐level inbreeding depression for all life cycle stages for a population of the cleistogamous perennial Viola septemloba. Our estimate of cumulative inbreeding depression of 43% is substantially greater than two previous estimates for cleistogamous species, but it is not sufficient by itself to explain the maintenance of potentially outcrossed chasmogamous flowers.

Heterosis and outbreeding depression in crosses between natural populations of Arabidopsis thaliana.

Understanding the causes and architecture of genetic differentiation between natural populations is of central importance in evolutionary biology. Crosses between natural populations can result in heterosis if recessive or nearly recessive deleterious mutations have become fixed within populations because of genetic drift. Divergence between populations can also result in outbreeding depression because of genetic incompatibilities. The net fitness consequences of between-population crosses will be a balance between heterosis and outbreeding depression. We estimated the magnitude of heterosis and outbreeding depression in the highly selfing model plant Arabidopsis thaliana, by crossing replicate line pairs from two sets of natural populations (C↔R, B↔S) separated by similar geographic distances (Italy↔Sweden). We examined the contribution of different modes of gene action to overall differences in estimates of lifetime fitness and fitness components using joint scaling tests with parental, reciprocal F1 and F2, and backcross lines. One of these population pairs (C↔R) was previously demonstrated to be locally adapted, but locally maladaptive quantitative trait loci were also found, suggesting a role for genetic drift in shaping adaptive variation. We found markedly different genetic architectures for fitness and fitness components in the two sets of populations. In one (C↔R), there were consistently positive effects of dominance, indicating the masking of recessive or nearly recessive deleterious mutations that had become fixed by genetic drift. The other set (B↔S) exhibited outbreeding depression because of negative dominance effects. Additional studies are needed to explore the molecular genetic basis of heterosis and outbreeding depression, and how their magnitudes vary across environments.

Genetic basis of photosynthetic responses to cold in two locally adapted populations of Arabidopsis thaliana

Dramatic intraspecific natural variation in photosynthetic responses to cold is controlled mainly by a few large effect QTLs, and some cold-responsive photosynthesis QTLs may contribute to local adaptation.

Genetics of water use physiology in locally adapted Arabidopsis thaliana

Identifying the genetic basis of adaptation to climate has long been a goal in evolutionary biology and has applications in agriculture. Adaptation to drought represents one important aspect of local adaptation, and drought is the major factor limiting agricultural yield. We examined local adaptation between Sweden and Italy Arabidopsis thaliana ecotypes, which show contrasting levels of water availability in their local environments. To identify quantitative trait loci (QTL) controlling water use physiology traits and adaptive trait QTL (genomic regions where trait QTL and fitness QTL colocalize), we performed QTL mapping on 374F9 recombinant inbred lines in well-watered and terminal drought conditions. We found 72 QTL (32 in well-watered, 31 in drought, 9 for plasticity) across five water use physiology traits: δ(13)C, rosette area, dry rosette weight, leaf water content and percent leaf nitrogen. Some of these genomic regions colocalize with fitness QTL and with other physiology QTL in defined hotspots. In addition, we found evidence of both constitutive and inducible water use physiology QTL. Finally, we identified highly divergent candidate genes, in silico. Our results suggest that many genes with minor effects may influence adaptation through water use physiology and that pleiotropic water use physiology QTL have fitness consequences.