Unpacking Conditional Neutrality: Genomic Signatures of Selection on Conditionally Beneficial and Conditionally Deleterious Mutations

A Historical Misstep: Niche Shift to Specialisation Is Pushing Insular Ginger Towards an Evolutionary Dead End

Niche specialisation is a double-edged sword as it aids species in adapting to a particular environment but makes them more susceptible to environmental change, which may result in species extinction. Although it has long been debated whether niche specialisation necessarily falls into an 'evolutionary dead end', empirical evidence from a population genetics perspective remains scant, especially when comparing both ecological generalists and specialists simultaneously. In this study, we scrutinised two Taiwan endemic gingers (Zingiber pleiostachyum and Z. shuanglongense) to evaluate how their contrasting patterns in niche breadth evolution have shaped their evolutionary trajectories. We utilised a genome-wide sequencing approach to investigate the demographic histories of each species, assess their maladaptation to future climate change, and estimate their mutational loads. Our results revealed distinct demographic histories between these two gingers. Z. shuanglongense, as the specialist, despite an initial increase during the Last Glacial Maximum (~22 Kya), has been subjected to a long-term decrease in effective population size (Ne), while Z. pleiostachyum is on the contrary increasing, leading to a significantly larger current Ne. Furthermore, ecological specialists are much more vulnerable to future climate change and exhibit greater drift-associated deleterious mutations compared to generalists, directly affecting species' fitness. This study strongly supports the idea that the transition in niche breadth towards specialisation will push Z. shuanglongense perilously close to extinction and also sheds light on species conservation within limited migratory space.

Evolving a plant-beneficial bacterium in soil vs. nutrient-rich liquid culture has contrasting effects on in-soil fitness

Inoculation of plant-beneficial microbes into agricultural soils can improve crop growth, but such outcomes depend on microbial survival. Here, we assessed how exposure to prior environmental conditions impacts microbial in-soil fitness, particularly focusing on incubation in liquid culture as an unavoidable phase of inoculant production and on pre-incubation in target soils as a potential method to improve performance. We conducted experimental evolution on a phosphorus-solubilizing bacterial species, Priestia megaterium, in (i) soil only, (ii) liquid media only, and (iii) soil followed by liquid media, using population metagenomic sequencing to track mutations over time. Several typical in vitro evolutionary phenomena were observed in liquid media-incubated populations, including clonal interference, genetic hitchhiking, and mutation parallelism between replicate populations, particularly in the sporulation transcription factor spo0A. Liquid media-incubated populations also developed a clear fitness reduction in soil compared to the ancestral isolate. However, soil-incubated populations grew slowly, experienced far fewer generations despite longer absolute time, and accumulated minimal mutational changes. Correspondingly, soil-incubated populations did not display improved survival compared to the ancestral isolate in their target soils, though there did appear to be minor fitness reductions in unfamiliar soil. This work demonstrates that adaptation to liquid media and/or a native soil can impact bacterial fitness in new soil and that bacterial evolution in more complex real-world habitats does not closely resemble bacterial evolution in liquid media.

IMPORTANCE

Innovative solutions are needed to address emerging challenges in agriculture while reducing its environmental footprint. Management of soil microbiomes could contribute to this effort, as plant growth-promoting microorganisms provide key ecosystem services that support crops. Yet, inoculating beneficial microbes into farm soils yields unreliable results. We require a greater knowledge of the ecology of these taxa to improve their functioning in sustainable agroecosystems. In this report, we demonstrate that exposure to laboratory media and lingering adaptation to another soil can negatively impact the in-soil survival of a phosphorus-solubilizing bacterial species. We go further to highlight the underlying mutations that give rise to these patterns. These insights can be leveraged to improve our understanding of how soil-dwelling beneficial microorganisms adapt to different evolutionary pressures.

Locally Specific Genome-Wide Signatures of Adaptation to Environmental Variation at High Resolution in an Alpine Plant

Microevolutionary processes shape adaptive responses to heterogeneous environments, where these effects vary both among and within species. However, it remains largely unknown to which degree signatures of adaptation to environmental drivers can be detected based on the choice of spatial scale and genomic marker. We studied signatures of local adaptation across two levels of spatial extents, investigating complementary types of genomic variants-single-nucleotide polymorphisms (SNPs) and polymorphic transposable elements (TEs)-in populations of the alpine model plant species Arabis alpina . We coupled environmental factors, derived from remote sensed digital elevation models (DEMs) at very high resolution (0.5 m), with whole-genome sequencing data of 304 individuals across four populations. By comparing putatively adaptive loci detected between each local population versus a regional assessment including all populations simultaneously, we demonstrate that responses of A. alpina to similar amounts of abiotic variation are largely governed by local evolutionary processes. Furthermore, we find minimally overlapping signatures of local adaptation between SNPs and polymorphic TEs. Notably, functional annotations of candidate genes for adaptation revealed several symbiosis-related genes associated with the abiotic factors studied, which could represent selective pressures from biotic agents. Our results highlight the importance of considering different spatial extents and types of genomic polymorphisms when searching for signatures of adaptation to environmental variation. Such insights provide key information on microevolutionary processes and could guide management decisions to mitigate negative impacts of climate change on alpine plant populations.

Spatial variation in genomic signatures of local adaptation during the cane toad invasion of Australia

Adaptive evolution can facilitate species' range expansions across environmentally heterogeneous landscapes. However, serial founder effects can limit the efficacy of selection, and the evolution of increased dispersal during range expansions may result in gene flow swamping local adaptation. Here, we study how genetic drift, gene flow and selection interact during the cane toad's (Rhinella marina) invasion across the heterogeneous landscape of Australia. Following its introduction in 1935, the cane toad colonised eastern Australia and established several stable range edges. The ongoing, more rapid range expansion in north-central Australia has occurred concomitant with an evolved increase in dispersal capacity. Using reduced representation genomic data of Australian cane toads from the expansion front and from two areas of their established range, we test the hypothesis that high gene flow constrains local adaptation at the expansion front relative to established areas. Genetic analyses indicate the three study areas are genetically distinct but show similar levels of allelic richness, heterozygosity and inbreeding. Markedly higher gene flow or recency of colonisation at the expansion front have likely hindered local adaptation at the time of sampling, as indicated by reduced slopes of genetic-environment associations (GEAs) estimated using a novel application of geographically weighted regression that accounts for allele surfing; GEA slopes are significantly steeper in established parts of the range. Our work bolsters evidence supporting adaptation of invasive species post-introduction and adds novel evidence for differing strengths of evolutionary forces among geographic areas with different invasion histories.

Conditionally Deleterious Mutation Load Accumulates in Genomic Islands of Local Adaptation but Can Be Purged with Sufficient Genotypic Redundancy

AbstractLocal adaptation frequently evolves in patches or environments that are connected via migration. In these cases, genomic regions that are linked to a locally adapted locus experience reduced effective migration rates. Via individual-based simulations of a two-patch system, we show that this reduced effective migration results in the accumulation of conditionally deleterious mutations, but not universally deleterious mutations, adjacent to adaptive loci. When there is redundancy in the genetic basis of local adaptation (i.e., genotypic redundancy), turnover of locally adapted polymorphisms allows conditionally deleterious mutation load to be purged. The amount of mutational load that accumulates adjacent to locally adapted loci is dependent on redundancy, recombination rate, migration rate, population size, strength of selection, and the phenotypic effect size of adaptive alleles. Our results highlight the need to be cautious when interpreting patterns of local adaptation at the level of phenotype or fitness, as the genetic basis of local adaptation can be transient, and evolution may confer a degree of maladaptation to nonlocal environments.

Genetic load and viability of a future restored northern white rhino population

As biodiversity loss outpaces recovery, conservationists are increasingly turning to novel tools for preventing extinction, including cloning and in vitro gametogenesis of biobanked cells. However, restoration of populations can be hindered by low genetic diversity and deleterious genetic load. The persistence of the northern white rhino (Ceratotherium simum cottoni) now depends on the cryopreserved cells of 12 individuals. These banked genomes have higher genetic diversity than southern white rhinos (C. s. simum), a sister subspecies that successfully recovered from a severe bottleneck, but the potential impact of genetic load is unknown. We estimated how demographic history has shaped genome-wide genetic load in nine northern and 13 southern white rhinos. The bottleneck left southern white rhinos with more fixed and homozygous deleterious alleles and longer runs of homozygosity, whereas northern white rhinos retained more deleterious alleles masked in heterozygosity. To gauge the impact of genetic load on the fitness of a northern white rhino population restored from biobanked cells, we simulated recovery using fitness of southern white rhinos as a benchmark for a viable population. Unlike traditional restoration, cell-derived founders can be reintroduced in subsequent generations to boost lost genetic diversity and relieve inbreeding. In simulations with repeated reintroduction of founders into a restored population, the fitness cost of genetic load remained lower than that borne by southern white rhinos. Without reintroductions, rapid growth of the restored population (>20-30% per generation) would be needed to maintain comparable fitness. Our results suggest that inbreeding depression from genetic load is not necessarily a barrier to recovery of the northern white rhino and demonstrate how restoration from biobanked cells relieves some constraints of conventional restoration from a limited founder pool. Established conservation methods that protect healthy populations will remain paramount, but emerging technologies hold promise to bolster these tools to combat the extinction crisis.

Functional genomic diversity is correlated with neutral genomic diversity in populations of an endangered rattlesnake

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.

The Contribution of Mutation to Variation in Temperature-Dependent Sprint Speed in Zebrafish, Danio rerio

AbstractThe contribution of new mutations to phenotypic variation and the consequences of this variation for individual fitness are fundamental concepts for understanding genetic variation and adaptation. Here, we investigated how mutation influenced variation in a complex trait in zebrafish, Danio rerio. Typical of many ecologically relevant traits in ectotherms, swimming speed in fish is temperature dependent, with evidence of adaptive evolution of thermal performance. We chemically induced novel germline point mutations in males and measured sprint speed in their sons at six temperatures (between 16°C and 34°C). Heterozygous mutational effects on speed were strongly positively correlated among temperatures, resulting in statistical support for only a single axis of mutational variation, reflecting temperature-independent variation in speed (faster-slower mode). These results suggest pleiotropic effects on speed across different temperatures; however, spurious correlations arise via linkage or heterogeneity in mutation number when mutations have consistent directional effects on each trait. Here, mutation did not change mean speed, indicating no directional bias in mutational effects. The results contribute to emerging evidence that mutations may predominantly have synergistic cross-environment effects, in contrast to conditionally neutral or antagonistic effects that underpin thermal adaptation. We discuss several aspects of experimental design that may affect resolution of mutations with nonsynergistic effects.

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.

The contribution of historical processes to contemporary extinction risk in placental mammals

Species persistence can be influenced by the amount, type, and distribution of diversity across the genome, suggesting a potential relationship between historical demography and resilience. In this study, we surveyed genetic variation across single genomes of 240 mammals that compose the Zoonomia alignment to evaluate how historical effective population size (Ne) affects heterozygosity and deleterious genetic load and how these factors may contribute to extinction risk. We find that species with smaller historical Ne carry a proportionally larger burden of deleterious alleles owing to long-term accumulation and fixation of genetic load and have a higher risk of extinction. This suggests that historical demography can inform contemporary resilience. Models that included genomic data were predictive of species' conservation status, suggesting that, in the absence of adequate census or ecological data, genomic information may provide an initial risk assessment.

Genome-Wide Association Studies across Environmental and Genetic Contexts Reveal Complex Genetic Architecture of Symbiotic Extended Phenotypes

A goal of modern biology is to develop the genotype-phenotype (G→P) map, a predictive understanding of how genomic information generates trait variation that forms the basis of both natural and managed communities. As microbiome research advances, however, it has become clear that many of these traits are symbiotic extended phenotypes, being governed by genetic variation encoded not only by the host's own genome, but also by the genomes of myriad cryptic symbionts. Building a reliable G→P map therefore requires accounting for the multitude of interacting genes and even genomes involved in symbiosis. Here, we use naturally occurring genetic variation in 191 strains of the model microbial symbiont Sinorhizobium meliloti paired with two genotypes of the host Medicago truncatula in four genome-wide association studies (GWAS) to determine the genomic architecture of a key symbiotic extended phenotype-partner quality, or the fitness benefit conferred to a host by a particular symbiont genotype, within and across environmental contexts and host genotypes. We define three novel categories of loci in rhizobium genomes that must be accounted for if we want to build a reliable G→P map of partner quality; namely, (i) loci whose identities depend on the environment, (ii) those that depend on the host genotype with which rhizobia interact, and (iii) universal loci that are likely important in all or most environments. IMPORTANCE Given the rapid rise of research on how microbiomes can be harnessed to improve host health, understanding the contribution of microbial genetic variation to host phenotypic variation is pressing, and will better enable us to predict the evolution of (and select more precisely for) symbiotic extended phenotypes that impact host health. We uncover extensive context-dependency in both the identity and functions of symbiont loci that control host growth, which makes predicting the genes and pathways important for determining symbiotic outcomes under different conditions more challenging. Despite this context-dependency, we also resolve a core set of universal loci that are likely important in all or most environments, and thus, serve as excellent targets both for genetic engineering and future coevolutionary studies of symbiosis.

A Pleiotropic Flowering Time QTL Exhibits Gene-by-Environment Interaction for Fitness in a Perennial Grass

Appropriate flowering time is a crucial adaptation impacting fitness in natural plant populations. Although the genetic basis of flowering variation has been extensively studied, its mechanisms in nonmodel organisms and its adaptive value in the field are still poorly understood. Here, we report new insights into the genetic basis of flowering time and its effect on fitness in Panicum hallii, a native perennial grass. Genetic mapping in populations derived from inland and coastal ecotypes identified flowering time quantitative trait loci (QTL) and many exhibited extensive QTL-by-environment interactions. Patterns of segregation within recombinant hybrids provide strong support for directional selection driving ecotypic divergence in flowering time. A major QTL on chromosome 5 (q-FT5) was detected in all experiments. Fine-mapping and expression studies identified a gene with orthology to a rice FLOWERING LOCUS T-like 9 (PhFTL9) as the candidate underlying q-FT5. We used a reciprocal transplant experiment to test for local adaptation and the specific impact of q-FT5 on performance. We did not observe local adaptation in terms of fitness tradeoffs when contrasting ecotypes in home versus away habitats. However, we observed that the coastal allele of q-FT5 conferred a fitness advantage only in its local habitat but not at the inland site. Sequence analyses identified an excess of low-frequency polymorphisms at the PhFTL9 promoter in the inland lineage, suggesting a role for either selection or population expansion on promoter evolution. Together, our findings demonstrate the genetic basis of flowering variation in a perennial grass and provide evidence for conditional neutrality underlying flowering time divergence.

Standing Genetic Diversity and Transmission Bottleneck Size Drive Adaptation in Bacteriophage Qβ

A critical issue to understanding how populations adapt to new selective pressures is the relative contribution of the initial standing genetic diversity versus that generated de novo. RNA viruses are an excellent model to study this question, as they form highly heterogeneous populations whose genetic diversity can be modulated by factors such as the number of generations, the size of population bottlenecks, or exposure to new environment conditions. In this work, we propagated at nonoptimal temperature (43 °C) two bacteriophage Qβ populations differing in their degree of heterogeneity. Deep sequencing analysis showed that, prior to the temperature change, the most heterogeneous population contained some low-frequency mutations that had previously been detected in the consensus sequences of other Qβ populations adapted to 43 °C. Evolved populations with origin in this ancestor reached similar growth rates, but the adaptive pathways depended on the frequency of these standing mutations and the transmission bottleneck size. In contrast, the growth rate achieved by populations with origin in the less heterogeneous ancestor did depend on the transmission bottleneck size. The conclusion is that viral diversification in a particular environment may lead to the emergence of mutants capable of accelerating adaptation when the environment changes.

A polygenic architecture with habitat-dependent effects underlies ecological differentiation in Silene

Ecological differentiation can drive speciation but it is unclear how the genetic architecture of habitat-dependent fitness contributes to lineage divergence. We investigated the genetic architecture of cumulative flowering, a fitness component, in second-generation hybrids between Silene dioica and Silene latifolia transplanted into the natural habitat of each species. We used reduced-representation sequencing and Bayesian sparse linear mixed models (BSLMMs) to analyze the genetic control of cumulative flowering in each habitat. Our results point to a polygenic architecture of cumulative flowering. Allelic effects were mostly beneficial or deleterious in one habitat and neutral in the other. Positive-effect alleles often were derived from the native species, whereas negative-effect alleles, at other loci, tended to originate from the non-native species. We conclude that ecological differentiation is governed and maintained by many loci with small, habitat-dependent effects consistent with conditional neutrality. This pattern may result from differences in selection targets in the two habitats and from environmentally dependent deleterious load. Our results further suggest that selection for native alleles and against non-native alleles acts as a barrier to gene flow between species.

Genetic load: genomic estimates and applications in non-model animals

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

The genetic architecture underlying prey-dependent performance in a microbial predator

Natural selection should favour generalist predators that outperform specialists across all prey types. Two genetic solutions could explain why intraspecific variation in predatory performance is, nonetheless, widespread: mutations beneficial on one prey type are costly on another (antagonistic pleiotropy), or mutational effects are prey-specific, which weakens selection, allowing variation to persist (relaxed selection). To understand the relative importance of these alternatives, we characterised natural variation in predatory performance in the microbial predator Dictyostelium discoideum. We found widespread nontransitive differences among strains in predatory success across different bacterial prey, which can facilitate stain coexistence in multi-prey environments. To understand the genetic basis, we developed methods for high throughput experimental evolution on different prey (REMI-seq). Most mutations (~77%) had prey-specific effects, with very few (~4%) showing antagonistic pleiotropy. This highlights the potential for prey-specific effects to dilute selection, which would inhibit the purging of variation and prevent the emergence of an optimal generalist predator.

What prevents a generalist predator from evolving and outperforming specialist predators? By combing analyses of natural variation with experimental evolution, Stewart et al. suggest that predator variation persists because most mutations have prey-specific effects, which results in relaxed selection

Physical geography, isolation by distance and environmental variables shape genomic variation of wild barley (Hordeum vulgare L. ssp. spontaneum) in the Southern Levant

Determining the extent of genetic variation that reflects local adaptation in crop-wild relatives is of interest for the purpose of identifying useful genetic diversity for plant breeding. We investigated the association of genomic variation with geographical and environmental factors in wild barley (Hordeum vulgare L. ssp. spontaneum) populations of the Southern Levant using genotyping by sequencing (GBS) of 244 accessions in the Barley 1K+ collection. The inference of population structure resulted in four genetic clusters that corresponded to eco-geographical habitats and a significant association between lower gene flow rates and geographical barriers, e.g. the Judaean Mountains and the Sea of Galilee. Redundancy analysis (RDA) revealed that spatial autocorrelation explained 45% and environmental variables explained 15% of total genomic variation. Only 4.5% of genomic variation was solely attributed to environmental variation if the component confounded with spatial autocorrelation was excluded. A synthetic environmental variable combining latitude, solar radiation, and accumulated precipitation explained the highest proportion of genomic variation (3.9%). When conditioned on population structure, soil water capacity was the most important environmental variable explaining 1.18% of genomic variation. Genome scans with outlier analysis and genome-environment association studies were conducted to identify adaptation signatures. RDA and outlier methods jointly detected selection signatures in the pericentromeric regions, which have reduced recombination, of the chromosomes 3H, 4H, and 5H. However, selection signatures mostly disappeared after correction for population structure. In conclusion, adaptation to the highly diverse environments of the Southern Levant over short geographical ranges had a limited effect on the genomic diversity of wild barley. This highlighted the importance of nonselective forces in genetic differentiation.

Evolution of polygenic traits under global vs local adaptation

Observations about the number, frequency, effect size, and genomic distribution of alleles associated with complex traits must be interpreted in light of evolutionary process. These characteristics, which constitute a trait's genetic architecture, can dramatically affect evolutionary outcomes in applications from agriculture to medicine, and can provide a window into how evolution works. Here, I review theoretical predictions about the evolution of genetic architecture under spatially homogeneous, global adaptation as compared with spatially heterogeneous, local adaptation. Due to the tension between divergent selection and migration, local adaptation can favor "concentrated" genetic architectures that are enriched for alleles of larger effect, clustered in a smaller number of genomic regions, relative to expectations under global adaptation. However, the evolution of such architectures may be limited by many factors, including the genotypic redundancy of the trait, mutation rate, and temporal variability of environment. I review the circumstances in which predictions differ for global vs local adaptation and discuss where progress can be made in testing hypotheses using data from natural populations and lab experiments. As the field of comparative population genomics expands in scope, differences in architecture among traits and species will provide insights into how evolution works, and such differences must be interpreted in light of which kind of selection has been operating.

Genetic trade-offs underlie divergent life history strategies for local adaptation in white clover

Local adaptation is common in plants, yet characterization of its underlying genetic basis is rare in herbaceous perennials. Moreover, while many plant species exhibit intraspecific chemical defence polymorphisms, their importance for local adaptation remains poorly understood. We examined the genetic architecture of local adaptation in a perennial, obligately-outcrossing herbaceous legume, white clover (Trifolium repens). This widespread species displays a well-studied chemical defence polymorphism for cyanogenesis (HCN release following tissue damage) and has evolved climate-associated cyanogenesis clines throughout its range. Two biparental F₂  mapping populations, derived from three parents collected in environments spanning the U.S. latitudinal species range (Duluth, MN, St. Louis, MO and Gainesville, FL), were grown in triplicate for two years in reciprocal common garden experiments in the parental environments (6,012 total plants). Vegetative growth and reproductive fitness traits displayed trade-offs across reciprocal environments, indicating local adaptation. Genetic mapping of fitness traits revealed a genetic architecture characterized by allelic trade-offs between environments, with 100% and 80% of fitness QTL in the two mapping populations showing significant QTL×E interactions, consistent with antagonistic pleiotropy. Across the genome there were three hotspots of QTL colocalization. Unexpectedly, we found little evidence that the cyanogenesis polymorphism contributes to local adaptation. Instead, divergent life history strategies in reciprocal environments were major fitness determinants: selection favoured early investment in flowering at the cost of multiyear survival in the southernmost site versus delayed flowering and multiyear persistence in the northern environments. Our findings demonstrate that multilocus genetic trade-offs contribute to contrasting life history characteristics that allow for local adaptation in this outcrossing herbaceous perennial.

The New Is Old: Novel Germination Strategy Evolved From Standing Genetic Variation in Weedy Rice

Feralization of crop plants has aroused an increasing interest in recent years, not only for the reduced yield and quality of crop production caused by feral plants but also for the rapid evolution of novel traits that facilitate the evolution and persistence of weedy forms. Weedy rice (Oryza sativa f. spontanea) is a conspecific weed of cultivated rice, with separate and independent origins. The weedy rice distributed in eastern and northeastern China did not diverge from their cultivated ancestors by reverting to the pre-domestication trait of seed dormancy during feralization. Instead, they developed a temperature-sensing mechanism to control the timing of seed germination. Subsequent divergence in the minimum critical temperature for germination has been detected between northeastern and eastern populations. An integrative analysis was conducted using combinations of phenotypic, genomic and transcriptomic data to investigate the genetic mechanism underlying local adaptation and feralization. A dozen genes were identified, which showed extreme allele frequency differences between eastern and northeastern populations, and high correlations between allele-specific gene expression and feral phenotypes. Trancing the origin of potential adaptive alleles based on genomic sequences revealed the presence of most selected alleles in wild and cultivated rice genomes, indicating that weedy rice drew upon pre-existing, "conditionally neutral" alleles to respond to the feral selection regimes. The cryptic phenotype was exposed by activating formerly silent alleles to facilitate the transition from cultivation to wild existence, promoting the evolution and persistence of weedy forms.

Consequences of mutation accumulation for growth performance are more likely to be resource-dependent at higher temperatures

Background

Mutation accumulation (MA) has profound ecological and evolutionary consequences. One example is that accumulation of conditionally neutral mutations leads to fitness trade-offs among heterogenous habitats which cause population divergence. Here we suggest that temperature, which controls the rates of all biochemical and biophysical processes, should play a crucial role for determining mutational effects. Particularly, warmer temperatures may mitigate the effects of some, not all, deleterious mutations and cause stronger environmental dependence in MA effects.

Results

We experimentally tested the above hypothesis by measuring the growth performance of ten Escherichia coli genotypes on six carbon resources across ten temperatures, where the ten genotypes were derived from a single ancestral strain and accumulated spontaneous mutations. We analyzed resource dependence of MA consequences for growth yields. The MA genotypes typically showed reduced growth yields relative to the ancestral type; and the magnitude of reduction was smaller at intermediate temperatures. Stronger resource dependence in MA consequences for growth performance was observed at higher temperatures. Specifically, the MA genotypes were more likely to show impaired growth performance on all the six carbon resources when grown at lower temperatures; but suffered growth performance loss only on some, not all the six, carbon substrates at higher temperatures.

Conclusions

Higher temperatures increase the chance that MA causes conditionally neutral fitness effects while MA is more likely to cause fitness loss regardless of available resources at lower temperatures. This finding has implications for understanding how geographic patterns in population divergence may emerge, and how conservation practices, particularly protection of diverse microhabitats, may mitigate the impacts of global warming.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01846-1.

Does geographic variation in thermal tolerance in Daphnia represent trade-offs or conditional neutrality?

Geographic variation in thermal tolerance in Daphnia seems to represent genetic load at the loci specifically responsible for heat tolerance resulting from conditional neutrality. We see no evidence of trade-offs between fitness-related traits at 25 °C vs. 10 °C or between two algal diets across Daphnia magna clones from a variety of locations representing the opposite ends of the distribution of long-term heat tolerance. Likewise, we found no evidence of within-environment trade-offs between heat tolerance and fitness-related traits in any of the environments. Neither short-term and long-term heat tolerance shows any consistent relationship with lipid fluorescence polarization and lipid peroxidation across clones or environments. Pervasive positive correlations between fitness-related traits indicate differences in genetic load rather than trade-off based local adaptation or thermal specialization. For heat tolerance such differences may be caused by either relaxation of stabilizing selection due to lower exposure to high temperature extremes, i.e., conditional neutrality, or by small effective population size followed by the recent range expansion.

Adaptation and competition in deteriorating environments

Evolution might rescue populations from extinction in changing environments. Using experimental evolution with microalgae, we investigated if competition influences adaptation to an abiotic stressor, and vice versa, if adaptation to abiotic change influences competition. In a first set of experiments, we propagated monocultures of five species with and without increasing salt stress for approximately 180 generations. When assayed in monoculture, two of the five species showed signatures of adaptation, that is, lines with a history of salt stress had higher population growth rates at high salt than lines without prior exposure to salt. When assayed in mixtures of species, however, only one of these two species had increased population size at high salt, indicating that competition can alter how adaptation to abiotic change influences population dynamics. In a second experiment, we cultivated two species in monocultures and in pairs, with and without increasing salt. While we found no effect of competition on adaptation to salt, our experiment revealed that evolutionary responses to salt can influence competition. Specifically, one of the two species had reduced competitive ability in the no-salt environment after long-term exposure to salt stress. Collectively, our results highlight the complex interplay of adaptation to abiotic change and competitive interactions.

Interspecific hybrids show a reduced adaptive potential under DNA damaging conditions

Hybridization may increase the probability of adaptation to extreme stresses. This advantage could be caused by an increased genome plasticity in hybrids, which could accelerate the search for adaptive mutations. High ultraviolet (UV) radiation is a particular challenge in terms of adaptation because it affects the viability of organisms by directly damaging DNA, while also challenging future generations by increasing mutation rate. Here we test whether hybridization accelerates adaptive evolution in response to DNA damage, using yeast as a model. We exposed 180 populations of hybrids between species (Saccharomyces cerevisiae and Saccharomyces paradoxus) and their parental strains to UV mimetic and control conditions for approximately 100 generations. Although we found that adaptation occurs in both hybrids and parents, hybrids achieved a lower rate of adaptation, contrary to our expectations. Adaptation to DNA damage conditions comes with a large and similar cost for parents and hybrids, suggesting that this cost is not responsible for the lower adaptability of hybrids. We suggest that the lower adaptive potential of hybrids in this condition may result from the interaction between DNA damage and the inherent genetic instability of hybrids.

Global adaptation complicates the interpretation of genome scans for local adaptation

Spatially varying selection promotes variance in allele frequencies, increasing genetic differentiation between the demes of a metapopulation. For that reason, outliers in the genome‐wide distribution of summary statistics measuring genetic differentiation, such as F_ST, are often interpreted as evidence for alleles that contribute to local adaptation. However, theoretical studies have shown that in spatially structured populations the spread of beneficial mutations with spatially uniform fitness effects can also induce transient genetic differentiation. In recent years, numerous empirical studies have suggested that such species‐wide, or global, adaptation makes a substantial contribution to molecular evolution. In this perspective, we discuss how commonly such global adaptation may influence the genome‐wide distribution of F_ST and generate genetic differentiation patterns, which could be mistaken for local adaptation. To illustrate this, we use forward‐in‐time population genetic simulations assuming parameters for the rate and strength of beneficial mutations consistent with estimates from natural populations. We demonstrate that the spread of globally beneficial mutations in parapatric populations may frequently generate F_ST outliers, which could be misinterpreted as evidence for local adaptation. The spread of beneficial mutations causes selective sweeps at flanking sites, so in some cases, the effects of global versus local adaptation may be distinguished by examining patterns of nucleotide diversity within and between populations in addition to F_ST. However, when local adaptation has been only recently established, it may be much more difficult to distinguish from global adaptation, due to less accumulation of linkage disequilibrium at flanking sites. Through our discussion, we conclude that a large fraction of F_ST outliers that are presumed to arise from local adaptation may instead be due to global adaptation.

Understanding Maladaptation by Uniting Ecological and Evolutionary Perspectives

Evolutionary biologists have long trained their sights on adaptation, focusing on the power of natural selection to produce relative fitness advantages while often ignoring changes in absolute fitness. Ecologists generally have taken a different tack, focusing on changes in abundance and ranges that reflect absolute fitness while often ignoring relative fitness. Uniting these perspectives, we articulate various causes of relative and absolute maladaptation and review numerous examples of their occurrence. This review indicates that maladaptation is reasonably common from both perspectives, yet often in contrasting ways. That is, maladaptation can appear strong from a relative fitness perspective, yet populations can be growing in abundance. Conversely, resident individuals can appear locally adapted (relative to nonresident individuals) yet be declining in abundance. Understanding and interpreting these disconnects between relative and absolute maladaptation, as well as the cases of agreement, is increasingly critical in the face of accelerating human-mediated environmental change. We therefore present a framework for studying maladaptation, focusing in particular on the relationship between absolute and relative fitness, thereby drawing together evolutionary and ecological perspectives. The unification of these ecological and evolutionary perspectives has the potential to bring together previously disjunct research areas while addressing key conceptual issues and specific practical problems.