Temporal Variation in Selection Influences Microgeographic Local Adaptation.

AbstractEcological heterogeneity can lead to local adaptation when populations exhibit fitness trade-offs among habitats. However, the degree to which local adaptation is affected by the spatial and temporal scale of environmental variation is poorly understood. A multiyear reciprocal transplant experiment was performed with populations of the annual plant Leptosiphon parviflorus living on adjacent serpentine and nonserpentine soil. Local adaptation over this small geographic scale was observed, but there were differences in the temporal variability of selection across habitats. On serpentine soil, the local population had a consistently large survival advantage, presumably as a result of the temporal stability in selection imposed by soil cation content. In contrast, a fecundity advantage was observed for the sandstone population on its native soil type but only in the two study years with the highest rainfall. A manipulative greenhouse experiment demonstrated that the fitness advantage of the sandstone population in its native soil type depends critically on water availability. The temporal variability in local adaptation driven by variation in precipitation suggests that continued drought conditions have the potential to erode local adaptation in these populations. These results show how different selective factors can influence spatial and temporal patterns of variation in fitness trade-offs.

Ecological genetics of local adaptation in Arabidopsis: An 8-year field experiment.

There is considerable evidence for local adaptation in nature, yet important questions remain regarding its genetic basis. How many loci are involved? What are their effect sizes? What is the relative importance of conditional neutrality versus genetic trade-offs? Here we address these questions in the self-pollinating, annual plant Arabidopsis thaliana. We used 400 recombinant inbred lines (RILs) derived from two locally adapted populations in Italy and Sweden, grew the RILs and parents at the parental locations, and mapped quantitative trait loci (QTL) for mean fitness (fruits/seedling planted). We previously published results from the first 3 years of the study, and here add five additional years, providing a unique opportunity to assess how temporal variation in selection might affect QTL detection and classification. We found 10 adaptive and one maladaptive QTL in Italy, and six adaptive and four maladaptive QTL in Sweden. The discovery of maladaptive QTL at both sites suggests that even locally adapted populations are not always at their genotypic optimum. Mean effect sizes for adaptive QTL, 0.97 and 0.55 fruits in Italy and Sweden, respectively, were large relative to the mean fitness of the RILs (approximately 8 fruits/seedling planted at both sites). Both genetic trade-offs (four cases) and conditional neutrality (seven cases) contribute to local adaptation in this system. The 8-year dataset provided greater power to detect QTL and to estimate their locations compared to our previous 3-year study, identifying one new genetic trade-off and resolving one genetic trade-off into two conditionally adaptive QTL.

Genome assemblies and comparison of two Neotropical spiral gingers: Costus pulverulentus and C. lasius.

The spiral gingers (Costus L.) are a pantropical genus of herbaceous perennial monocots; the Neotropical clade of Costus radiated rapidly in the past few million years into over 60 species. The Neotropical spiral gingers have a rich history of evolutionary and ecological research that can motivate and inform modern genetic investigations. Here, we present the first 2 chromosome-level genome assemblies in the genus, for C. pulverulentus and C. lasius, and briefly compare their synteny. We assembled the C. pulverulentus genome from a combination of short-read data, Chicago and Dovetail Hi-C chromatin-proximity sequencing, and alignment with a linkage map. We annotated the genome by mapping a C. pulverulentus transcriptome and querying mapped transcripts against a protein database. We assembled the C. lasius genome with Pacific Biosciences HiFi long reads and alignment to the C. pulverulentus genome. These 2 assemblies are the first published genomes for non-cultivated tropical plants. These genomes solidify the spiral gingers as a model system and will facilitate research on the poorly understood genetic basis of tropical plant diversification.

Genetic and physiological mechanisms of freezing tolerance in locally adapted populations of a winter annual.

PREMISE: Despite myriad examples of local adaptation, the phenotypes and genetic variants underlying such adaptive differentiation are seldom known. Recent work on freezing tolerance and local adaptation in ecotypes of Arabidopsis thaliana from Italy and Sweden provides an essential foundation for uncovering the genotype-phenotype-fitness map for an adaptive response to a key environmental stress. METHODS: We examined the consequences of a naturally occurring loss-of-function (LOF) mutation in an Italian allele of the gene that encodes the transcription factor CBF2, which underlies a major freezing-tolerance locus. We used four lines with a Swedish genetic background, each containing a LOF CBF2 allele. Two lines had introgression segments containing the Italian CBF2 allele, and two contained deletions created using CRISPR-Cas9. We used a growth chamber experiment to quantify freezing tolerance and gene expression before and after cold acclimation. RESULTS: Freezing tolerance was lower in the Italian (11%) compared to the Swedish (72%) ecotype, and all four experimental CBF2 LOF lines had reduced freezing tolerance compared to the Swedish ecotype. Differential expression analyses identified 10 genes for which all CBF2 LOF lines, and the IT ecotype had similar patterns of reduced cold responsive expression compared to the SW ecotype. CONCLUSIONS: We identified 10 genes that are at least partially regulated by CBF2 that may contribute to the differences in cold-acclimated freezing tolerance between the Italian and Swedish ecotypes. These results provide novel insight into the molecular and physiological mechanisms connecting a naturally occurring sequence polymorphism to an adaptive response to freezing conditions.

Heterosis is common and inbreeding depression absent in natural populations of Arabidopsis thaliana.

The importance of genetic drift in shaping patterns of adaptive genetic variation in nature is poorly known. Genetic drift should drive partially recessive deleterious mutations to high frequency, and inter-population crosses may therefore exhibit heterosis (increased fitness relative to intra-population crosses). Low genetic diversity and greater genetic distance between populations should increase the magnitude of heterosis. Moreover, drift and selection should remove strongly deleterious recessive alleles from individual populations, resulting in reduced inbreeding depression. To estimate heterosis, we crossed 90 independent line pairs of Arabidopsis thaliana from 15 pairs of natural populations sampled across Fennoscandia and crossed an additional 41 line pairs from a subset of four of these populations to estimate inbreeding depression. We measured lifetime fitness of crosses relative to parents in a large outdoor common garden (8,448 plants in total) in central Sweden. To examine the effects of genetic diversity and genetic distance on heterosis, we genotyped parental lines for 869 SNPs. Overall, genetic variation within populations was low (median expected heterozygosity = 0.02), and genetic differentiation was high (median FST  = 0.82). Crosses between 10 of 15 population pairs exhibited significant heterosis, with magnitudes of heterosis as high as 117%. We found no significant inbreeding depression, suggesting that the observed heterosis is due to fixation of mildly deleterious alleles within populations. Widespread and substantial heterosis indicates an important role for drift in shaping genetic variation, but there was no significant relationship between fitness of crosses relative to parents and genetic diversity or genetic distance between populations.

Plant-herbivore coevolution and plant speciation.

More than five decades ago, Ehrlich and Raven proposed a revolutionary idea-that the evolution of novel plant defense could spur adaptive radiation in plants. Despite motivating much work on plant-herbivore coevolution and defense theory, Ehrlich and Raven never proposed a mechanism for their "escape and radiate" model. Recent intriguing mechanisms proposed by Marquis et al. include sympatric divergence, pleiotropic effects of plant defense traits on reproductive isolation, and strong postzygotic isolation, but these may not be general features of herbivore-mediated speciation. An alternate view is that herbivores may impose strong divergent selection on defenses in allopatric plant populations, with plant-herbivore coevolution driving local adaptation resulting in plant speciation. Building on these ideas, we propose three scenarios that consider the role of herbivores in plant speciation. These include (1) vicariance, subsequent coevolution within populations and adaptive divergence between geographically isolated populations, (2) colonization of a new habitat lacking effective herbivores followed by loss of defense and then re-evolution and coevolution of defense in response to novel herbivores, and (3) evolution of a new defense followed by range expansion, vicariance, and coevolution. We discuss the general role of coevolution in plant speciation and consider outstanding issues related to understanding: (1) the mechanisms behind cospeciation of plants and insects, (2) geographic variation in defense phenotypes, (3) how defensive traits and geography map on plant phylogenies, and (4) the role of herbivores in driving character displacement in defense phenotypes of related species in sympatry.

Latitudinal patterns of herbivore pressure in a temperate herb support the biotic interactions hypothesis.

The longstanding biotic interactions hypothesis predicts that herbivore pressure declines with latitude, but the evidence is mixed. To address gaps in previous studies, we measured herbivory and defence in the same system, quantified defence with bioassays, and considered effects of leaf age. We quantified herbivory and defence of young and mature leaves along a continental gradient in eastern North America in the native herb Phytolacca americana L. Herbivory in the field declined with latitude and was strongly correlated with lepidopteran abundance. Laboratory bioassays revealed that leaf palatability was positively correlated with latitude of origin. Young leaves were more damaged than mature leaves at lower latitudes in the field, but less palatable in bioassays. Both defence and palatability displayed non-linear latitudinal patterns, suggesting potential mechanisms based on biological or climatic thresholds. In sum, observational and experimental studies find patterns consistent with high herbivore pressure and stronger plant defences at lower latitudes.

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

Local adaptation is common, but the traits and genes involved are often unknown. Physiological responses to cold probably contribute to local adaptation in wide-ranging species, but the genetic basis underlying natural variation in these traits has rarely been studied. Using a recombinant inbred (495 lines) mapping population from locally adapted populations of Arabidopsis thaliana from Sweden and Italy, we grew plants at low temperature and mapped quantitative trait loci (QTLs) for traits related to photosynthesis: maximal quantum efficiency (Fv/Fm), rapidly reversible photoprotection (NPQfast), and photoinhibition of PSII (NPQslow) using high-throughput, whole-plant measures of chlorophyll fluorescence. In response to cold, the Swedish line had greater values for all traits, and for every trait, large effect QTLs contributed to parental differences. We found one major QTL affecting all traits, as well as unique major QTLs for each trait. Six trait QTLs overlapped with previously published locally adaptive QTLs based on fitness measured in the native environments over 3 years. Our results demonstrate that photosynthetic responses to cold can vary dramatically within a species, and may predominantly be caused by a few QTLs of large effect. Some photosynthesis traits and QTLs probably contribute to local adaptation in this system.

The Edaphic Environment Mediates Flowering-Time Differentiation Between Adjacent Populations of Leptosiphon Parviflorus.

Flowering time is an important life history trait in plants that often affects fitness. The optimal time to flower may be influenced by trade-offs between flowering time and growth-related traits and is thus likely to differ among habitats. Because flowering-time differences between populations can also reduce gene flow, understanding the factors that contribute to variation in flowering time among closely adjacent populations that experience gene flow is of particular interest. Plant adaptation to different edaphic environments provides some of the best examples of adaptive divergence at small spatial scales, and often coincides with flowering-time shifts. The current study addresses the causes of flowering-time differences in two populations of Leptosiphon parviflorus that are locally adapted to adjacent serpentine and sandstone soils despite moderate levels of gene flow and close geographic proximity. Field reciprocal-transplant studies and watering manipulations in the greenhouse demonstrate the contribution of both the genotype and the environment to observed flowering-time differences. The plasticity of flowering time in response to soil type appears to be driven by differences in soil moisture. In addition, selection on flowering time was measured in both soil types across 4 years of study using a set of F5 advanced-generation hybrids and found to differ between the habitats. Therefore, both selection and plasticity contribute to flowering-time differences between these populations and thus have likely played an important role in the initiation and/or maintenance of adaptive divergence in this system.

Speciation and the Latitudinal Diversity Gradient: Insights from the Global Distribution of Endemic Fish.

The nearly universal pattern that species richness increases from the poles to the equator (the latitudinal diversity gradient [LDG]) has been of intense interest since its discovery by early natural-history explorers. Among the many hypotheses proposed to explain the LDG, latitudinal variation in (1) productivity, (2) time and area available for diversification, and (3) speciation and/or extinction rates have recently received the most attention. Because tropical regions are older and were formerly more widespread, these factors are often intertwined, hampering efforts to distinguish their relative contributions to the LDG. Here we examine the global distribution of endemic lake fishes to determine how lake age, area, and latitude each affect the probability of speciation and the extent of diversification occurring within a lake. We analyzed the distribution of endemic fishes worldwide (1,933 species and subspecies from 47 families in 2,746 lakes) and find that the probability of a lake containing an endemic species and the total number of endemics per lake increase with lake age and area and decrease with latitude. Moreover, the geographic locations of endemics in 34 of 41 families are found at lower latitudes than those of nonendemics. We propose that the greater diversification of fish at low latitudes may be driven in part by ecological opportunities promoted by tropical climates and by the coevolution of species interactions.

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.

Increased heterosis in selfing populations of a perennial forb.

Quantifying the importance of random genetic drift in natural populations is central to understanding the potential limits to natural selection. One approach is to estimate the magnitude of heterosis, the increased fitness of progeny derived from crosses between populations relative to crosses within populations caused by the heterozygous masking of deleterious recessive or nearly recessive alleles that have been fixed by drift within populations. Self-fertilization is expected to reduce the effective population size by half relative to outcrossing, and population bottlenecks may be common during the transition to selfing. Therefore, chance fixation of deleterious alleles due to drift in selfing populations should increase heterosis between populations. Increased homozygosity due to fixation or loss of alleles should also decrease inbreeding depression within populations. Most populations of the perennial herb Arabidopsis lyrata ssp. lyrata are self-incompatible (SI), but several have evolved self-compatibility and are highly selfing. We quantified heterosis and inbreeding depression in two predominantly self-compatible (SC) and seven SI populations in a field common garden experiment within the species' native range and examined the correlation between these metrics to gauge the similarity in their genetic basis. We measured proportion germination in the lab, and survival and fecundity (flower and seed production) for 2 years in the field, and calculated estimates of cumulative fitness. We found 7.2-fold greater heterosis in SC compared with SI populations, despite substantial heterosis in SI populations (56 %). Inbreeding depression was >61 %, and not significantly different between SC and SI populations. There was no correlation between population estimates of heterosis and inbreeding depression, suggesting that they have somewhat different genetic bases. Combined with other sources of information, our results suggest a history of bottlenecks in all of these populations. The bottlenecks in SC populations may have been severe, but their strong inbreeding depression remains enigmatic.

Ecological differentiation and local adaptation in two sister species of Neotropical Costus (Costaceae).

Reciprocal transplant experiments have often provided evidence of local adaptation in temperate plants, but few such studies have been conducted in the tropics. To enhance our knowledge of local adaptation in tropical plants, we studied natural populations of two recently diverged Neotropical plant species, Costus allenii and C. villosissimus, in central Panama. We found that these species display a parapatric distribution that reflects local environmental differences on a fine geographic scale: C. allenii is found along ravines in the understory of primary forest, while C. villosissimus is found along forest edges. Light availability was lower in C. allenii habitats, while precipitation and soil moisture were lower in C. villosissimus habitats. We carried out reciprocal transplant experiments with seeds and clones of mature plants to test the hypothesis that the parapatric distribution of these species is due to divergent adaptation to their local habitats. We found strong evidence of local adaptation, i.e., when grown in their "home" sites, each species outperformed the species from an "away" site. Our finding that C. allenii and C. villosissimus are mainly isolated by their microhabitats provides a first step toward understanding the mechanisms of adaptation and speciation in the tropics.

Ecological and evolutionary perspectives on community assembly.

Ecologists often view community assembly as a process involving the dispersal of species from a static regional species pool followed by environmental filtering to establish the local community. This conceptual framework ignores the dynamic nature of species pools and fails to recognize that communities are assembled by processes operating over a vast range of temporal and spatial scales. Species pool richness and composition are influenced by metacommunity dynamics over short timescales and by speciation, extinction, and dispersal over long timescales. We suggest that a stronger focus on the geography of speciation, the formation of secondary sympatry, and the feedback between local and regional processes is needed to fully understand community assembly and the importance of dynamic species pools.

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.

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.

Genetic mapping of adaptation reveals fitness tradeoffs in Arabidopsis thaliana.

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.