Local adaptation: Causal agents of selection and adaptive trait divergence

Evolution of parasites in the Anthropocene: new pressures, new adaptive directions

The Anthropocene is seeing the human footprint rapidly spreading to all of Earth's ecosystems. The fast-changing biotic and abiotic conditions experienced by all organisms are exerting new and strong selective pressures, and there is a growing list of examples of human-induced evolution in response to anthropogenic impacts. No organism is exempt from these novel selective pressures. Here, we synthesise current knowledge on human-induced evolution in eukaryotic parasites of animals, and present a multidisciplinary framework for its study and monitoring. Parasites generally have short generation times and huge fecundity, features that predispose them for rapid evolution. We begin by reviewing evidence that parasites often have substantial standing genetic variation, and examples of their rapid evolution both under conditions of livestock production and in serial passage experiments. We then present a two-step conceptual overview of the causal chain linking anthropogenic impacts to parasite evolution. First, we review the major anthropogenic factors impacting parasites, and identify the selective pressures they exert on parasites through increased mortality of either infective stages or adult parasites, or through changes in host density, quality or immunity. Second, we discuss what new phenotypic traits are likely to be favoured by the new selective pressures resulting from altered parasite mortality or host changes; we focus mostly on parasite virulence and basic life-history traits, as these most directly influence the transmission success of parasites and the pathology they induce. To illustrate the kinds of evolutionary changes in parasites anticipated in the Anthropocene, we present a few scenarios, either already documented or hypothetical but plausible, involving parasite taxa in livestock, aquaculture and natural systems. Finally, we offer several approaches for investigations and real-time monitoring of rapid, human-induced evolution in parasites, ranging from controlled experiments to the use of state-of-the-art genomic tools. The implications of fast-evolving parasites in the Anthropocene for disease emergence and the dynamics of infections in domestic animals and wildlife are concerning. Broader recognition that it is not only the conditions for parasite transmission that are changing, but the parasites themselves, is needed to meet better the challenges ahead.

Biotic interactions promote local adaptation to soil in plants

Although different ecological factors shape adaptative evolution in natural habitats, we know little about how their interactions impact local adaptation. Here we used eight generations of experimental evolution with outcrossing Brassica rapa plants as a model system, in eight treatment groups that varied in soil type, herbivory (with/without aphids), and pollination mode (hand- or bumblebee-pollination), to study how biotic interactions affect local adaptation to soil. First, we show that several plant traits evolved in response to biotic interactions in a soil-specific way. Second, using a reciprocal transplant experiment, we demonstrate that significant local adaptation to soil-type evolved in the "number of open flowers", a trait used as a fitness proxy, but only in plants that evolved with herbivory and bee pollination. Whole genome re-sequencing of experimental lines revealed that biotic interactions caused a 10-fold increase in the number of SNPs across the genome with significant allele frequency change, and that alleles with opposite allele frequency change in different soil types (antagonistic pleiotropy) were most common in plants with an evolutionary history of herbivory and bee pollination. Our results demonstrate that the interaction with mutualists and antagonists can facilitate local adaptation to soil type through antagonistic pleiotropy.

Population origin and heritable effects mediate road salt toxicity and thermal stress in an amphibian

Human impacts on wild populations are numerous and extensive, degrading habitats and causing population declines across taxa. Though these impacts are often studied individually, wild populations typically face suites of stressors acting concomitantly, compromising the fitness of individuals and populations in ways poorly understood and not easily predicted by the effects of any single stressor. Developing understanding of the effects of multiple stressors and their potential interactions remains a critical challenge in environmental biology. Here, we focus on assessing the impacts of two prominent stressors associated with anthropogenic activities that affect many organisms across the planet - elevated salinity (e.g., from road de-icing salt) and temperature (e.g. from climate change). We examined a suite of physiological traits and components of fitness across populations of wood frogs originating from ponds that differ in their proximity to roads and thus their legacy of exposure to pollution from road salt. When experimentally exposed to road salt, wood frogs showed reduced survival (especially those from ponds adjacent to roads), divergent developmental rates, and reduced longevity. Family-level effects mediated these outcomes, but high salinity generally eroded family-level variance. When combined, exposure to both temperature and salt resulted in very low survival, and this effect was strongest in roadside populations. Taken together, these results suggest that temperature is an important stressor capable of exacerbating impacts from a prominent contaminant confronting many freshwater organisms in salinized habitats. More broadly, it appears likely that toxicity might often be underestimated in the absence of multi-stressor approaches.

Predicting species invasiveness with genomic data: Is genomic offset related to establishment probability?

Predicting the risk of establishment and spread of populations outside their native range represents a major challenge in evolutionary biology. Various methods have recently been developed to estimate population (mal)adaptation to a new environment with genomic data via so-called Genomic Offset (GO) statistics. These approaches are particularly promising for studying invasive species but have still rarely been used in this context. Here, we evaluated the relationship between GO and the establishment probability of a population in a new environment using both in silico and empirical data. First, we designed invasion simulations to evaluate the ability to predict establishment probability of two GO computation methods (Geometric GO and Gradient Forest) under several conditions. Additionally, we aimed to evaluate the interpretability of absolute Geometric GO values, which theoretically represent the adaptive genetic distance between populations from distinct environments. Second, utilizing public empirical data from the crop pest species Bactrocera tryoni, a fruit fly native from Northern Australia, we computed GO between "source" populations and a diverse range of locations within invaded areas. This practical application of GO within the context of a biological invasion underscores its potential in providing insights and guiding recommendations for future invasion risk assessment. Overall, our results suggest that GO statistics represent good predictors of the establishment probability and may thus inform invasion risk, although the influence of several factors on prediction performance (e.g., propagule pressure or admixture) will need further investigation.

Natural genetic variation in dynamic photosynthesis is correlated with stomatal anatomical traits in diverse tomato species across geographical habitats

Plants grown under field conditions experience fluctuating light. Understanding the natural genetic variations for a similarly dynamic photosynthetic response among untapped germplasm resources, as well as the underlying mechanisms, may offer breeding strategies to improve production using molecular approaches. Here, we measured gas exchange under fluctuating light, along with stomatal density and size, in eight wild tomato species and two tomato cultivars. The photosynthetic induction response showed significant diversity, with some wild species having faster induction rates than the two cultivars. Species with faster photosynthetic induction rates had higher daily integrated photosynthesis, but lower average water use efficiency because of high stomatal conductance under natural fluctuating light. The variation in photosynthetic induction was closely associated with the speed of stomatal responses, highlighting its critical role in maximizing photosynthesis under fluctuating light conditions. Moreover, stomatal size was negatively correlated with stomatal density within a species, and plants with smaller stomata at a higher density had a quicker photosynthetic response than those with larger stomata at lower density. Our findings show that the response of stomatal conductance plays a pivotal role in photosynthetic induction, with smaller stomata at higher density proving advantageous for photosynthesis under fluctuating light in tomato species. The interspecific variation in the rate of stomatal responses could offer an untapped resource for optimizing dynamic photosynthetic responses under field conditions.

The evolutionary genomics of adaptation to stress in wild rhizobium bacteria

Microbiota comprise the bulk of life's diversity, yet we know little about how populations of microbes accumulate adaptive diversity across natural landscapes. Adaptation to stressful soil conditions in plants provides seminal examples of adaptation in response to natural selection via allelic substitution. For microbes symbiotic with plants however, horizontal gene transfer allows for adaptation via gene gain and loss, which could generate fundamentally different evolutionary dynamics. We use comparative genomics and genetics to elucidate the evolutionary mechanisms of adaptation to physiologically stressful serpentine soils in rhizobial bacteria in western North American grasslands. In vitro experiments demonstrate that the presence of a locus of major effect, the nre operon, is necessary and sufficient to confer adaptation to nickel, a heavy metal enriched to toxic levels in serpentine soil, and a major axis of environmental soil chemistry variation. We find discordance between inferred evolutionary histories of the core genome and nreAXY genes, which often reside in putative genomic islands. This suggests that the evolutionary history of this adaptive variant is marked by frequent losses, and/or gains via horizontal acquisition across divergent rhizobium clades. However, different nre alleles confer distinct levels of nickel resistance, suggesting allelic substitution could also play a role in rhizobium adaptation to serpentine soil. These results illustrate that the interplay between evolution via gene gain and loss and evolution via allelic substitution may underlie adaptation in wild soil microbiota. Both processes are important to consider for understanding adaptive diversity in microbes and improving stress-adapted microbial inocula for human use.

Weak local adaptation to drought in seedlings of a widespread conifer

Tree seedlings from populations native to drier regions are often assumed to be more drought tolerant than those from wetter provenances. However, intraspecific variation in drought tolerance has not been well-characterized despite being critical for developing climate change mitigation and adaptation strategies, and for predicting the effects of drought on forests. We used a large-scale common garden drought-to-death experiment to assess range-wide variation in drought tolerance, measured by decline of photosynthetic efficiency, growth, and plastic responses to extreme summer drought in seedlings of 73 natural populations of the two main varieties of Douglas-fir (Pseudotsuga menziesii var. menziesii and var. glauca). Local adaptation to drought was weak in var. glauca and nearly absent in menziesii. Var. glauca showed higher tolerance to drought but slower growth than var. menziesii. Clinal variation in drought tolerance and growth species-wide was mainly associated with temperature rather than precipitation. A higher degree of plasticity for growth was observed in var. menziesii in response to extreme drought. Genetic variation for drought tolerance in seedlings within varieties is maintained primarily within populations. Selective breeding within populations may facilitate adaptation to drought more than assisted gene flow.

A large-effect fitness trade-off across environments is explained by a single mutation affecting cold acclimation

Identifying the genetic basis of local adaptation and fitness trade-offs across environments is a central goal of evolutionary biology. Cold acclimation is an adaptive plastic response for surviving seasonal freezing, and costs of acclimation may be a general mechanism for fitness trade-offs across environments in temperate zone species. Starting with locally adapted ecotypes of Arabidopsis thaliana from Italy and Sweden, we examined the fitness consequences of a naturally occurring functional polymorphism in CBF2. This gene encodes a transcription factor that is a major regulator of cold-acclimated freezing tolerance and resides within a locus responsible for a genetic trade-off for long-term mean fitness. We estimated the consequences of alternate genotypes of CBF2 on 5-y mean fitness and fitness components at the native field sites by comparing near-isogenic lines with alternate genotypes of CBF2 to their genetic background ecotypes. The effects of CBF2 were validated at the nucleotide level using gene-edited lines in the native genetic backgrounds grown in simulated parental environments. The foreign CBF2 genotype in the local genetic background reduced long-term mean fitness in Sweden by more than 10%, primarily via effects on survival. In Italy, fitness was reduced by more than 20%, primarily via effects on fecundity. At both sites, the effects were temporally variable and much stronger in some years. The gene-edited lines confirmed that CBF2 encodes the causal variant underlying this genetic trade-off. Additionally, we demonstrated a substantial fitness cost of cold acclimation, which has broad implications for potential maladaptive responses to climate change.

The genomic consequences of selection across development

Understanding how natural selection drives diversification in nature has been at the forefront of biological research for over a century. The main idea is simple: natural selection favours individuals best suited to pass on their genes. However, the journey from birth to reproduction is complex as organisms experience multiple developmental stages, each influenced by genetic and environmental factors (Orr, 2009). These complexities compound even further as each stage of development might be governed by a unique underlying set of alleles and genes. In this issue of Molecular Ecology, Goebl et al. (2022) examine the role of natural selection in driving ecotypic divergence across different life history stages of the prairie sunflower Helianthus petiolaris. The authors used reciprocal transplant experiments, demographic models, and genomic sequencing to explore fitness variation across developmental stages. They show how natural selection impacts population divergence across multiple life history stages and evaluate the resulting allele frequency changes. Goebl et al. link these results to the role of chromosomal inversions, thus furthering our understanding of how ecological divergence proceeds in the face of gene flow. Below, we explore these results in detail and complement their interpretation by considering the evolution of genetic correlations amongst traits governing fitness.

Ecophysiological adaptations shape distributions of closely related trees along a climatic moisture gradient

Tradeoffs between the energetic benefits and costs of traits can shape species and trait distributions along environmental gradients. Here we test predictions based on such tradeoffs using survival, growth, and 50 photosynthetic, hydraulic, and allocational traits of ten Eucalyptus species grown in four common gardens along an 8-fold gradient in precipitation/pan evaporation (P/Ep) in Victoria, Australia. Phylogenetically structured tests show that most trait-environment relationships accord qualitatively with theory. Most traits appear adaptive across species within gardens (indicating fixed genetic differences) and within species across gardens (indicating plasticity). However, species from moister climates have lower stomatal conductance than others grown under the same conditions. Responses in stomatal conductance and five related traits appear to reflect greater mesophyll photosynthetic sensitivity of mesic species to lower leaf water potential. Our data support adaptive cross-over, with realized height growth of most species exceeding that of others in climates they dominate. Our findings show that pervasive physiological, hydraulic, and allocational adaptations shape the distributions of dominant Eucalyptus species along a subcontinental climatic moisture gradient, driven by rapid divergence in species P/Ep and associated adaptations.

A few essential genetic loci distinguish Penstemon species with flowers adapted to pollination by bees or hummingbirds

In the formation of species, adaptation by natural selection generates distinct combinations of traits that function well together. The maintenance of adaptive trait combinations in the face of gene flow depends on the strength and nature of selection acting on the underlying genetic loci. Floral pollination syndromes exemplify the evolution of trait combinations adaptive for particular pollinators. The North American wildflower genus Penstemon displays remarkable floral syndrome convergence, with at least 20 separate lineages that have evolved from ancestral bee pollination syndrome (wide blue-purple flowers that present a landing platform for bees and small amounts of nectar) to hummingbird pollination syndrome (bright red narrowly tubular flowers offering copious nectar). Related taxa that differ in floral syndrome offer an attractive opportunity to examine the genomic basis of complex trait divergence. In this study, we characterized genomic divergence among 229 individuals from a Penstemon species complex that includes both bee and hummingbird floral syndromes. Field plants are easily classified into species based on phenotypic differences and hybrids displaying intermediate floral syndromes are rare. Despite unambiguous phenotypic differences, genome-wide differentiation between species is minimal. Hummingbird-adapted populations are more genetically similar to nearby bee-adapted populations than to geographically distant hummingbird-adapted populations, in terms of genome-wide dXY. However, a small number of genetic loci are strongly differentiated between species. These approximately 20 "species-diagnostic loci," which appear to have nearly fixed differences between pollination syndromes, are sprinkled throughout the genome in high recombination regions. Several map closely to previously established floral trait quantitative trait loci (QTLs). The striking difference between the diagnostic loci and the genome as whole suggests strong selection to maintain distinct combinations of traits, but with sufficient gene flow to homogenize the genomic background. A surprisingly small number of alleles confer phenotypic differences that form the basis of species identity in this species complex.

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.

Standing genetic variation and chromosome differences drove rapid ecotype formation in a major malaria mosquito

Species distributed across heterogeneous environments often evolve locally adapted ecotypes, but understanding of the genetic mechanisms involved in their formation and maintenance in the face of gene flow is incomplete. In Burkina Faso, the major African malaria mosquito Anopheles funestus comprises two strictly sympatric and morphologically indistinguishable yet karyotypically differentiated forms reported to differ in ecology and behavior. However, knowledge of the genetic basis and environmental determinants of An. funestus diversification was impeded by lack of modern genomic resources. Here, we applied deep whole-genome sequencing and analysis to test the hypothesis that these two forms are ecotypes differentially adapted to breeding in natural swamps versus irrigated rice fields. We demonstrate genome-wide differentiation despite extensive microsympatry, synchronicity, and ongoing hybridization. Demographic inference supports a split only ~1,300 y ago, closely following the massive expansion of domesticated African rice cultivation ~1,850 y ago. Regions of highest divergence, concentrated in chromosomal inversions, were under selection during lineage splitting, consistent with local adaptation. The origin of nearly all variations implicated in adaptation, including chromosomal inversions, substantially predates the ecotype split, suggesting that rapid adaptation was fueled mainly by standing genetic variation. Sharp inversion frequency differences likely facilitated adaptive divergence between ecotypes by suppressing recombination between opposing chromosomal orientations of the two ecotypes, while permitting free recombination within the structurally monomorphic rice ecotype. Our results align with growing evidence from diverse taxa that rapid ecological diversification can arise from evolutionarily old structural genetic variants that modify genetic recombination.

Breeding for Higher Yields of Wheat and Rice through Modifying Nitrogen Metabolism

Wheat and rice produce nutritious grains that provide 32% of the protein in the human diet globally. Here, we examine how genetic modifications to improve assimilation of the inorganic nitrogen forms ammonium and nitrate into protein influence grain yield of these crops. Successful breeding for modified nitrogen metabolism has focused on genes that coordinate nitrogen and carbon metabolism, including those that regulate tillering, heading date, and ammonium assimilation. Gaps in our current understanding include (1) species differences among candidate genes in nitrogen metabolism pathways, (2) the extent to which relative abundance of these nitrogen forms across natural soil environments shape crop responses, and (3) natural variation and genetic architecture of nitrogen-mediated yield improvement. Despite extensive research on the genetics of nitrogen metabolism since the rise of synthetic fertilizers, only a few projects targeting nitrogen pathways have resulted in development of cultivars with higher yields. To continue improving grain yield and quality, breeding strategies need to focus concurrently on both carbon and nitrogen assimilation and consider manipulating genes with smaller effects or that underlie regulatory networks as well as genes directly associated with nitrogen metabolism.

Taking quantitative genomics into the wild

We organized this special issue to highlight new work and review recent advances at the cutting edge of 'wild quantitative genomics'. In this editorial, we will present some history of wild quantitative genetic and genomic studies, before discussing the main themes in the papers published in this special issue and highlighting the future outlook of this dynamic field.