Invoking adaptation to decipher the genetic legacy of past climate change

Genetic admixture drives climate adaptation in the bank vole

Genetic admixture introduces new variants at relatively high frequencies, potentially aiding rapid responses to environmental changes. Here, we evaluate its role in adaptive variation related to climatic conditions in bank voles (Clethrionomys glareolus) in Britain, using whole-genome data. Our results reveal loci showing excess ancestry from one of the two postglacial colonist populations inconsistent with overall admixture patterns. Notably, loci associated with climate adaptation exhibit disproportionate amounts of excess ancestry, highlighting the impact of admixture between colonist populations on local adaptation. The results suggest strong and localized selection on climate-adaptive loci, as indicated by steep clines and/or shifted cline centres, during population replacement. A subset, including a haemoglobin gene, is associated with oxidative stress responses, underscoring a role of oxidative stress in local adaptation. Our study highlights the important contribution of admixture during secondary contact between populations from distinct climatic refugia enriching adaptive diversity. Understanding these dynamics is crucial for predicting future adaptive capacity to anthropogenic climate change.

Future carbon sequestration potential in a widespread transcontinental boreal tree species: Standing genetic variation matters!

Climate change (CC) necessitates reforestation/afforestation programs to mitigate its impacts and maximize carbon sequestration. But comprehending how tree growth, a proxy for fitness and resilience, responds to CC is critical to maximize these programs' effectiveness. Variability in tree response to CC across populations can notably be influenced by the standing genetic variation encompassing both neutral and adaptive genetic diversity. Here, a framework is proposed to assess tree growth potential at the population scale while accounting for standing genetic variation. We applied this framework to black spruce (BS, Picea mariana [Mill] B.S.P.), with the objectives to (1) determine the key climate variables having impacted BS growth response from 1974 to 2019, (2) examine the relative roles of local adaptation and the phylogeographic structure in this response, and (3) project BS growth under two Shared Socioeconomic Pathways while taking standing genetic variation into account. We modeled growth using a machine learning algorithm trained with dendroecological and genetic data obtained from over 2600 trees (62 populations divided in three genetic clusters) in four 48-year-old common gardens, and simulated growth until year 2100 at the common garden locations. Our study revealed that high summer and autumn temperatures negatively impacted BS growth. As a consequence of warming, this species is projected to experience a decline in growth by the end of the century, suggesting maladaptation to anticipated CC and a potential threat to its carbon sequestration capacity. This being said, we observed a clear difference in response to CC within and among genetic clusters, with the western cluster being more impacted than the central and eastern clusters. Our results show that intraspecific genetic variation, notably associated with the phylogeographic structure, must be considered when estimating the response of widespread species to CC.

Environmental and Biogeographic Drivers behind Alpine Plant Thermal Tolerance and Genetic Variation

In alpine ecosystems, elevation broadly functions as a steep thermal gradient, with plant communities exposed to regular fluctuations in hot and cold temperatures. These conditions lead to selective filtering, potentially contributing to species-level variation in thermal tolerance and population-level genetic divergence. Few studies have explored the breadth of alpine plant thermal tolerances across a thermal gradient or the underlying genetic variation thereof. We measured photosystem heat (Tcrit-hot) and cold (Tcrit-cold) thresholds of ten Australian alpine species across elevation gradients and characterised their neutral genetic variation. To reveal the biogeographical drivers of present-day genetic signatures, we also reconstructed temporal changes in habitat suitability across potential distributional ranges. We found intraspecific variation in thermal thresholds, but this was not associated with elevation, nor underpinned by genetic differentiation on a local scale. Instead, regional population differentiation and considerable homozygosity within populations may, in part, be driven by distributional contractions, long-term persistence, and migrations following habitat suitability. Our habitat suitability models suggest that cool-climate-distributed alpine plants may be threatened by a warming climate. Yet, the observed wide thermal tolerances did not reflect this vulnerability. Conservation efforts should seek to understand variations in species-level thermal tolerance across alpine microclimates.

Climate-related naturally occurring epimutation and their roles in plant adaptation in A. thaliana

DNA methylation has been proposed to be an important mechanism that allows plants to respond to their environments sometimes entirely uncoupled from genetic variation. To understand the genetic basis, biological functions and climatic relationships of DNA methylation at a population scale in Arabidopsis thaliana, we performed a genome-wide association analysis with high-quality single nucleotide polymorphisms (SNPs), and found that ~56% on average, especially in the CHH sequence context (71%), of the differentially methylated regions (DMRs) are not tagged by SNPs. Among them, a total of 3235 DMRs are significantly associated with gene expressions and potentially heritable. 655 of the 3235 DMRs are associated with climatic variables, and we experimentally verified one of them, HEI10 (HUMAN ENHANCER OF CELL INVASION NO.10). Such epigenetic loci could be subjected to natural selection thereby affecting plant adaptation, and would be expected to be an indicator of accessions at risk. We therefore incorporated these climate-related DMRs into a gradient forest model, and found that the natural A. thaliana accessions in Southern Europe that may be most at risk under future climate change. Our findings highlight the importance of integrating DNA methylation that is independent of genetic variations, and climatic data to predict plants' vulnerability to future climate change.

Rolling with the punches: Organism-environment interactions shape spatial pattern of adaptive differentiation in the widespread mantis shrimp Oratosquilla oratoria

Investigating spatial pattern of adaptive variation and its underlying processes can inform the adaptive potential distributed within species ranges, which is increasingly important in the context of a changing climate. A correct interpretation of adaptive variation pattern requires that population history and the ensuing population genetic structure are taken into account. Here we carried out such a study by integrating population genomic analyses, demographic model testing and species distribution modeling to investigate patterns and causes of adaptive differentiation in a widespread mantis shrimp, Oratosquilla oratoria, along a replicated, broad-scale temperature gradient in the northwestern Pacific (NWP). Our results supported a strong hierarchical ecogeographic structure dominated by habitat-linked divergence among O. oratoria populations accompanied with introgressive hybridization. A combined FST outlier and environmental correlation analyses revealed remarkable temperature-associated clines in allele frequency across paired North-South populations on Chinese and Japanese coasts, and identified a suite of loci associated with temperature adaptation. Further demographic model testing revealed the observed clinal variation derived partly from Pleistocene divergence followed by recent secondary contact. More importantly, the likelihood of hybridization is predicted to increase as climate change progresses, which would break barriers to gene flow and enable the spread of adaptive genetic variation. These results support that not only is temperature-driven adaptive differentiation occurs in O. oratoria but that such pattern is likely attributed to ancient adaptive variation, sustained by contemporary ocean conditions and a semi-permeable barrier to gene flow maintained by selection. They moreover provide genomic insights into the distribution of adaptive potential across O. oratoria' s species range. This work can serve as a case study to characterize adaptive diversity of marine species in the NWP by integrating environmental and genetic data at temporal and spatial scales in a population genomic framework, which would improve management and conservation actions under climate change.

Landscape genomics reveals adaptive divergence of indigenous sheep in different ecological environments of Xinjiang, China

Sheep are important livestock animals that have evolved under various ecological pressures. Xinjiang is a region with diverse and harsh environments that have shaped many local sheep breeds with unique characteristics and environmental adaptability. However, these breeds are losing ecological flexibility due to the promotion of intensive farming practices. Here we sequenced 14 local sheep breeds from Xinjiang and analyzed their genetic structure and gene flow with other sheep breeds from neighboring regions. The Tibetan Plateau was the geographic origin of Xinjiang native sheep evolution. We performed genome-environment association analysis and identified Bio9: Mean Temperature of Driest Quarter and Bio15: Precipitation Seasonality as the key environmental factors affecting Xinjiang local sheep and the key genes involved in their survival and adaptation. We classified Xinjiang native sheep breeds into six groups based on their differential genes by pairwise selective sweep analysis and Community Network Analysis. We analyzed transcriptome expression data of 832 sheep tissues and detected tissue-specific enrichment of six group-specific genes in different biological systems. Our results revealed the genetic basis of year-round estrus, drought tolerance, hypoxia resistance, and cold tolerance traits of Xinjiang sheep breeds. Moreover, we proposed conservation strategies for Xinjiang local sheep breeds and provided theoretical guidance for breeding new sheep breeds under global extreme environments.

Local adaptation shapes functional traits and resource allocation in black spruce

Climate change is rapidly altering weather patterns, resulting in shifts in climatic zones. The survival of trees in specific locations depends on their functional traits. Local populations exhibit trait adaptations that ensure their survival and accomplishment of growth and reproduction processes during the growing season. Studying these traits offers valuable insights into species responses to present and future environmental conditions, aiding the implementation of measures to ensure forest resilience and productivity. This study investigates the variability in functional traits among five black spruce (Picea mariana (Mill.) B.S.P.) provenances originating from a latitudinal gradient along the boreal forest, and planted in a common garden in Quebec, Canada. We examined differences in bud phenology, growth performance, lifetime first reproduction, and the impact of a late-frost event on tree growth and phenological adjustments. The findings revealed that trees from northern sites exhibit earlier budbreak, lower growth increments, and reach reproductive maturity earlier than those from southern sites. Late-frost damage affected growth performance, but no phenological adjustment was observed in the successive year. Local adaptation in the functional traits may lead to maladaptation of black spruce under future climate conditions or serve as a potent evolutionary force promoting rapid adaptation under changing environmental conditions.

Local adaptation and future climate vulnerability in a wild rodent

As climate change continues, species pushed outside their physiological tolerance limits must adapt or face extinction. When change is rapid, adaptation will largely harness ancestral variation, making the availability and characteristics of that variation of critical importance. Here, we used whole-genome sequencing and genetic-environment association analyses to identify adaptive variation and its significance in the context of future climates in a small Palearctic mammal, the bank vole (Clethrionomys glareolus). We found that peripheral populations of bank vole in Britain are already at the extreme bounds of potential genetic adaptation and may require an influx of adaptive variation in order to respond. Analyses of adaptive loci suggest regional differences in climate variables select for variants that influence patterns of population adaptive resilience, including genes associated with antioxidant defense, and support a pattern of thermal/hypoxic cross-adaptation. Our findings indicate that understanding potential shifts in genomic composition in response to climate change may be key to predicting species' fate under future climates.

Reforestation of high elevation pines: Direct seeding success depends on seed source and sowing environment

Forest persistence in regions impacted by increasing water and temperature stress will depend upon species' ability to either rapidly adjust to novel conditions or migrate to track ecological niches. Predicted, rapid climate change is likely to outpace the adaptive and migratory capacity of long-lived isolated tree species, and reforestation may be critical to species' persistence. Facilitating persistence both within and beyond a species' range requires identification of seed lots best adapted to the current and future conditions predicted with rapid climate change. We evaluate variation in emergent seedling performance that leads to differential survival among species and populations for three high elevation five-needle pines. We paired a fully reciprocal field common garden experiment with a greenhouse common garden study to (1) quantify variation in seedling emergence and functional traits, (2) ask how functional traits affect performance under different establishment conditions, and (3) evaluate whether trait and performance variation demonstrates local adaptation and plasticity. Among study species-limber, Great Basin bristlecone, and whitebark pines-we found divergence in emergence and functional traits, though soil moisture was the strongest driver of seedling emergence and abundance across all species. Generalist limber pine had a clear emergence advantage as well as traits associated with drought adaptation, while edaphic specialist bristlecone pine was characterized by low emergence yet high early survival once established. Despite evidence for edaphic specialization, soil characteristics alone did not explain bristlecone success. Across species, trait-environment relationships provided some evidence for local adaptation in drought-adapted traits, but we found no evidence of local adaptation in emergence or survival at this early life stage. For managers looking to promote persistence, sourcing seed from drier environments is likely to impart greater drought resistance into reforestation efforts through strategies such as greater root investment, increasing the probability of early seedling survival. This research demonstrates, through a rigorous reciprocal transplant experimental design, that it may be possible to select climate- and soil-appropriate seed sources for reforestation. However, planting success will ultimately rely on a suitable establishment environment, requiring careful consideration of interannual climate variability for management interventions in these climate and disturbance-impacted tree species.

Invasion success of a Lessepsian symbiont-bearing foraminifera linked to high dispersal ability, preadaptation and suppression of sexual reproduction

Among the most successful Lessepsian invaders is the symbiont-bearing benthic foraminifera Amphistegina lobifera. In its newly conquered habitat, this prolific calcifier and ecosystem engineer is exposed to environmental conditions that exceed the range of its native habitat. To disentangle which processes facilitated the invasion success of A. lobifera into the Mediterranean Sea we analyzed a ~ 1400 bp sequence fragment covering the SSU and ITS gene markers to compare the populations from its native regions and along the invasion gradient. The genetic variability was studied at four levels: intra-genomic, population, regional and geographical. We observed that the invasion is not associated with genetic differentiation, but the invasive populations show a distinct suppression of intra-genomic variability among the multiple copies of the rRNA gene. A reduced genetic diversity compared to the Indopacific is observed already in the Red Sea populations and their high dispersal potential into the Mediterranean appears consistent with a bridgehead effect resulting from the postglacial expansion from the Indian Ocean into the Red Sea. We conclude that the genetic structure of the invasive populations reflects two processes: high dispersal ability of the Red Sea source population pre-adapted to Mediterranean conditions and a likely suppression of sexual reproduction in the invader. This discovery provides a new perspective on the cost of invasion in marine protists: The success of the invasive A. lobifera in the Mediterranean Sea comes at the cost of abandonment of sexual reproduction.

Historical biogeography and local adaptation explain population genetic structure in a widespread terrestrial orchid

BACKGROUND AND AIMS

Historical changes in environmental conditions and colonization-extinction dynamics have a direct impact on the genetic structure of plant populations. However, understanding how past environmental conditions influenced the evolution of species with high gene flow is challenging when signals for genetic isolation and adaptation are swamped by gene flow. We investigated the spatial distribution and genetic structure of the widespread terrestrial orchid Epipactis helleborine to identify glacial refugia, characterize postglacial population dynamics and assess its adaptive potential.

METHODS

Ecological niche modelling was used to locate possible glacial refugia and postglacial recolonization opportunities of E. helleborine. A large single-nucleotide polymorphism (SNP) dataset obtained through genotyping by sequencing was used to define population genetic diversity and structure and to identify sources of postglacial gene flow. Outlier analyses were used to elucidate how adaptation to the local environment contributed to population divergence.

KEY RESULTS

The distribution of climatically suitable areas was restricted during the Last Glacial Maximum to the Mediterranean, south-western Europe and small areas in the Alps and Carpathians. Within-population genetic diversity was high in E. helleborine (mean expected heterozygosity, 0.373 ± 0.006; observed heterozygosity, 0.571 ± 0.012; allelic richness, 1.387 ± 0.007). Italy and central Europe are likely to have acted as important genetic sources during postglacial recolonization. Adaptive SNPs were associated with temperature, elevation and precipitation.

CONCLUSIONS

Forests in the Mediterranean and Carpathians are likely to have acted as glacial refugia for Epipactis helleborine. Postglacial migration northwards and to higher elevations resulted in the dispersal and diversification of E. helleborine in central Europe and Italy, and to geographical isolation and divergent adaptation in Greek and Italian populations. Distinguishing adaptive from neutral genetic diversity allowed us to conclude that E. helleborine has a high adaptive potential to climate change and demonstrates that signals of adaptation and historical isolation can be identified even in species with high gene flow.

Climate-induced range shifts drive adaptive response via spatio-temporal sieving of alleles

Quaternary climate fluctuations drove many species to shift their geographic ranges, in turn shaping their genetic structures. Recently, it has been argued that adaptation may have accompanied species range shifts via the "sieving" of genotypes during colonisation and establishment. However, this has not been directly demonstrated, and knowledge remains limited on how different evolutionary forces, which are typically investigated separately, interacted to jointly mediate species responses to past climatic change. Here, through whole-genome re-sequencing of over 1200 individuals of the carnation Dianthus sylvestris coupled with integrated population genomic and gene-environment models, we reconstruct the past neutral and adaptive landscape of this species as it was shaped by the Quaternary glacial cycles. We show that adaptive responses emerged concomitantly with the post-glacial range shifts and expansions of this species in the last 20 thousand years. This was due to the heterogenous sieving of adaptive alleles across space and time, as populations expanded out of restrictive glacial refugia into the broader and more heterogeneous range of habitats available in the present-day inter-glacial. Our findings reveal a tightly-linked interplay of migration and adaptation under past climate-induced range shifts, which we show is key to understanding the spatial patterns of adaptive variation we see in species today.

Neutral and adaptive genetic diversity in plants: An overview

Genetic diversity is a prerequisite for evolutionary change in all kinds of organisms. It is generally acknowledged that populations lacking genetic variation are unable to evolve in response to new environmental conditions (e.g., climate change) and thus may face an increased risk of extinction. Although the importance of incorporating genetic diversity into the design of conservation measures is now well understood, less attention has been paid to the distinction between neutral (NGV) and adaptive (AGV) genetic variation. In this review, we first focus on the utility of NGV by examining the ways to quantify it, reviewing applications of NGV to infer ecological and evolutionary processes, and by exploring its utility in designing conservation measures for plant populations and species. Against this background, we then summarize the ways to identify and estimate AGV and discuss its potential use in plant conservation. After comparing NGV and AGV and considering their pros and cons in a conservation context, we conclude that there is an urgent need for a better understanding of AGV and its role in climate change adaptation. To date, however, there are only a few AGV studies on non-model plant species aimed at deciphering the genetic and genomic basis of complex trait variation. Therefore, conservation researchers and practitioners should keep utilizing NGV to develop relevant strategies for rare and endangered plant species until more estimates of AGV are available.

Jack pine of all trades: Deciphering intraspecific variability of a key adaptive trait at the rear edge of a widespread fire-embracing North American conifer

PREMISE

Understanding mechanisms fostering long-term persistence of marginal populations should provide key insights about species resilience facing climate change. Cone serotiny is a key adaptive trait in Pinus banksiana (jack pine), which shows phenotypic variation according to the fire regime. Compared to range-core populations within the fire-prone boreal forest, low and variable serotiny in rear-edge populations suggest local adaptation to uncommon and unpredictable wildfire regime. We assessed environmental/physiological factors that might modulate intraspecific variation in cone serotiny.

METHODS

We experimentally subjected closed cones to incrementing temperatures, then tested seed germination to determine whether and how various ecological factors (cone age, branch height, tree size, tree age) are related to cone dehiscence and seed viability in jack pines from rear-edge and range-core populations in eastern Canada.

RESULTS

Cones from rear-edge populations dehisce at a lower opening temperature, which increases with cone age. Cones from range-core stands open at a more constant, yet higher temperature. Cones from rear-edge stands take between 13 and 27 years to reach the level of serotiny achieved at the range core. At the rear edge, seed viability is steady (51%), whereas it decreases from 70% to 30% in 20 years at the range core.

CONCLUSIONS

We inferred the mechanisms of a bet-hedging strategy in rear-edge populations, which ensures steady recruitment during fire-free intervals and successful postfire regeneration. This capacity to cope with infrequent and unpredictable fire regime should increase the resilience of jack pine populations as global changes alter fire dynamics of the boreal forest.

Phylogeny of Saxifraga section Saxifraga subsection Arachnoideae (Saxifragaceae) and the origin of low elevation shade-dwelling species

Saxifraga section Saxifraga subsection Arachnoideae is a lineage of 12 species distributed mainly in the European Alps. It is unusual in terms of ecological diversification by containing both high elevation species from exposed alpine habitats and low elevation species from shady habitats such as overhanging rocks and cave entrances. Our aims are to explore which of these habitat types is ancestral, and to identify the possible drivers of this remarkable ecological diversification. Using a Hybseq DNA-sequencing approach and a complete species sample we reconstructed and dated the phylogeny of subsection Arachnoideae. Using Landolt indicator values, this phylogenetic tree was used for the reconstruction of the evolution of temperature, light and soil pH requirements in this lineage. Diversification of subsection Arachnoideae started in the late Pliocene and continued through the Pleistocene. Both diversification among and within clades was largely allopatric, and species from shady habitats with low light requirements are distributed in well-known refugia. We hypothesize that low light requirements evolved when species persisting in cold-stage refugia were forced into marginal habitats by more competitive warm-stage vegetation. While we do not claim that such competition resulted in speciation, it very likely resulted in adaptive evolution.

Gene-by-environment interactions in plants: Molecular mechanisms, environmental drivers, and adaptive plasticity

Plants demonstrate a broad range of responses to environmental shifts. One of the most remarkable responses is plasticity, which is the ability of a single plant genotype to produce different phenotypes in response to environmental stimuli. As with all traits, the ability of plasticity to evolve depends on the presence of underlying genetic diversity within a population. A common approach for evaluating the role of genetic variation in driving differences in plasticity has been to study genotype-by-environment interactions (G × E). G × E occurs when genotypes produce different phenotypic trait values in response to different environments. In this review, we highlight progress and promising methods for identifying the key environmental and genetic drivers of G × E. Specifically, methodological advances in using algorithmic and multivariate approaches to understand key environmental drivers combined with new genomic innovations can greatly increase our understanding about molecular responses to environmental stimuli. These developing approaches can be applied to proliferating common garden networks that capture broad natural environmental gradients to unravel the underlying mechanisms of G × E. An increased understanding of G × E can be used to enhance the resilience and productivity of agronomic systems.

Genomic insights into local adaptation and future climate-induced vulnerability of a keystone forest tree in East Asia

Rapid global climate change is posing a substantial threat to biodiversity. The assessment of population vulnerability and adaptive capacity under climate change is crucial for informing conservation and mitigation strategies. Here we generate a chromosome-scale genome assembly and re-sequence genomes of 230 individuals collected from 24 populations for Populus koreana, a pioneer and keystone tree species in temperate forests of East Asia. We integrate population genomics and environmental variables to reveal a set of climate-associated single-nucleotide polymorphisms, insertion/deletions and structural variations, especially numerous adaptive non-coding variants distributed across the genome. We incorporate these variants into an environmental modeling scheme to predict a highly spatiotemporal shift of this species in response to future climate change. We further identify the most vulnerable populations that need conservation priority and many candidate genes and variants that may be useful for forest tree breeding with special aims. Our findings highlight the importance of integrating genomic and environmental data to predict adaptive capacity of a key forest to rapid climate change in the future.

Chromosomal inversions can limit adaptation to new environments

Chromosomal inversions are often thought to facilitate local adaptation and population divergence because they can link multiple adaptive alleles into non-recombining genomic blocks. Selection should thus be more efficient in driving inversion-linked adaptive alleles to high frequency in a population, particularly in the face of maladaptive gene flow. But what if ecological conditions and hence selection on inversion-linked alleles change? Reduced recombination within inversions could then constrain the formation of optimal combinations of pre-existing alleles under these new ecological conditions. Here, we outline this idea of inversions limiting adaptation and divergence when ecological conditions change across time or space. We reason and use simulations to illustrate that the benefit of inversions for local adaptation and divergence under one set of ecological conditions can come with a concomitant constraint for adaptation to novel sets of ecological conditions. This limitation of inversions to adaptation may contribute to the maintenance of polymorphism within species.

Tracking habitat or testing its suitability? Similar distributional patterns can hide very different histories of persistence versus nonequilibrium dynamics

The expansions and contractions of a species' range in response to temporal changes in selective filters leave genetic signatures that can inform a more accurate reconstruction of their evolutionary history across the landscape. After a long period of continental decline, Australian rainforests settled into localized patterns of contraction or expansion during the climatic fluctuations of the Quaternary. The environmental impacts of recurring glacial and interglacial periods also intensified the arrival of new lineages from the Sunda shelf, and it can be expected that immigrant versus locally persistent taxa responded to environmental challenges in quantifiably different manner. To investigate how such differences impact on species' distribution, we contrast landscape genomic patterns and changes in habitat availability between a species with a long continental history on Doryphora sassafras and a Sunda-derived species (Toona ciliata), across a distributional overlap. Extensive landscape-level homogeneity across chloroplast and nuclear genomes for the Sunda-derived T. ciliata, characterize the genetic signature of a very recent invasion and a rapid southern "exploratory" expansion that had not been previously recorded in the Australian flora (i.e., of Gondwanan origin or Sahul-derived). In contrast, D. sassafras is consistent with other Sahul-derived species characterized by strong geographical divergence and regional differentiation. Interestingly, our findings suggest that admixture between genetically divergent populations during expansion events might be a contributing factor to the successful colonization of novel habitats. Overall, this study identifies some of the mechanisms regulating the rearrangements in species distributions and assemblage composition that follow major environmental shifts, and reminds us how a species' current range might not necessarily define species' habitat preference, with the consequence that estimates of past or future range might not always be reliable.

The Bank Vole (Clethrionomys glareolus) as a Model System for Adaptive Phylogeography in the European Theater

The legacy of climatic changes during the Pleistocene glaciations allows inferences to be made about the patterns and processes associated with range expansion/colonization, including evolutionary adaptation. With the increasing availability of population genomic data, we have the opportunity to examine these questions in detail and in a variety of non-traditional model species. As an exemplar, here we review more than two decades of work by our group and others that illustrate the potential of a single “non-model model” mammal species - the bank vole (Clethrionomys glareolus), which is particularly well suited to illustrate the complexities that may be associated with range expansion and the power of genomics (and other datasets) to uncover them. We first summarize early phylogeographic work using mitochondrial DNA and then describe new phylogeographic insights gained from population genomic analysis of genome-wide SNP data to highlight the bank vole as one of the most compelling examples of a forest mammal, that survived in cryptic extra-Mediterranean (“northern”) glacial refugia in Europe, and as one of the species in which substantial replacement and mixing of lineages originating from different refugia occurred during end-glacial colonization. Our studies of bank vole hemoglobin structure and function, as well as our recent ecological niche modeling study examining differences among bank vole lineages, led us to develop the idea of “adaptive phylogeography.” This is what we call the study of the role of adaptive differences among populations in shaping phylogeographic patterns. Adaptive phylogeography provides a link between past population history and adaptation that can ultimately help predict the potential of future species responses to climate change. Because the bank vole is part of a community of organisms whose range has repeatedly contracted and then expanded in the past, what we learn from the bank vole will be useful for our understanding of a broad range of species.

Landscape Genomics in Tree Conservation Under a Changing Environment

Understanding the genetic basis of how species respond to changing environments is essential to the conservation of species. However, the molecular mechanisms of adaptation remain largely unknown for long-lived tree species which always have large population sizes, long generation time, and extensive gene flow. Recent advances in landscape genomics can reveal the signals of adaptive selection linking genetic variations and landscape characteristics and therefore have created novel insights into tree conservation strategies. In this review article, we first summarized the methods of landscape genomics used in tree conservation and elucidated the advantages and disadvantages of these methods. We then highlighted the newly developed method "Risk of Non-adaptedness," which can predict the genetic offset or genomic vulnerability of species via allele frequency change under multiple scenarios of climate change. Finally, we provided prospects concerning how our introduced approaches of landscape genomics can assist policymaking and improve the existing conservation strategies for tree species under the ongoing global changes.

Rules of Plant Species Ranges: Applications for Conservation Strategies

Earth is changing rapidly and so are many plant species’ ranges. Here, we synthesize eco-evolutionary patterns found in plant range studies and how knowledge of species ranges can inform our understanding of species conservation in the face of global change. We discuss whether general biogeographic “rules” are reliable and how they can be used to develop adaptive conservation strategies of native plant species across their ranges. Rules considered include (1) factors that set species range limits and promote range shifts; (2) the impact of biotic interactions on species range limits; (3) patterns of abundance and adaptive properties across species ranges; (4) patterns of gene flow and their implications for genetic rescue, and (5) the relationship between range size and conservation risk. We conclude by summarizing and evaluating potential species range rules to inform future conservation and management decisions. We also outline areas of research to better understand the adaptive capacity of plants under environmental change and the properties that govern species ranges. We advise conservationists to extend their work to specifically consider peripheral and novel populations, with a particular emphasis on small ranges. Finally, we call for a global effort to identify, synthesize, and analyze prevailing patterns or rules in ecology to help speed conservation efforts.

Persistent yet vulnerable: resurvey of an Abies ecotone reveals few differences but vulnerability to climate change

Climate change is shifting forest tree species distributions across elevational and latitudinal gradients, and these changes are often pronounced at ecotones where species meet their climatic bounds and are replaced by other species. Using an extensive ecotone composed of lower-montane white fir (Abies concolor var. lowiana) and upper-montane red fir (Abies magnifica var. magnifica) in the central Sierra Nevada range of California, USA, we (1) examined how the demographics of the ecotone have responded to recent climate using a field observational study and a historical dataset, (2) quantified climate drivers across species life stages using contemporary demographic data, and (3) tested the potential impacts of future climate on species-specific seedling survival and growth in a fully factorial growth chamber experiment that varied temperature, growing season length, and water availability. A re-examination of the ecotone midpoint after 35 yr suggested a reduction in A. concolor sapling and tree densities and a rise in A. magnifica proportional dominance between surveys. Seedling abundances across the ecotone indicated that A. magnifica tends to dominate the regeneration layer and currently forms an important component of the seedling community at elevations below those where A. magnifica saplings or trees begin to co-dominate stands. Observational and experimental assessments suggest that temperature and precipitation serve as important drivers, differentiating A. concolor vs. A. magnifica distributions, and are primary stressors at the seedling stage. Seedlings of both species were adversely affected by experimental climate treatments, although A. concolor exhibited greater survival and a more conservative growth strategy under extreme climatic stress than A. magnifica. Projections indicate that historical climate conditions will rise by an amount greater than the ecotone's current elevational extent by the end of the 21st century. Differential drivers of species abundances suggest that the projected climate will expand conditions that promote A. concolor abundance and impede A. magnifica abundance across the ecotone; however, disturbance activity and microclimatic conditions will also influence regeneration and overstory tree dynamics. Our study demonstrates the importance of quantifying species-specific responses to climate and indicates that widespread regeneration failure may be one possible consequence in which species exhibit strong sensitivity to projected climate conditions.

Genetic Differentiation and Demographic History of Three Cerris Oak Species in China Based on Nuclear Microsatellite Makers

Knowledge of interspecific divergence and population expansions/contractions of dominant forest trees in response to geological events and climatic oscillations is of major importance to understand their evolution and demography. However, the interspecific patterns of genetic differentiation and spatiotemporal population dynamics of three deciduous Cerris oak species (Q. acutissima, Q. variabilis and Q. chenii) that are widely distributed in China remain poorly understood. In this study, we genotyped 16 nuclear loci in 759 individuals sampled from 44 natural populations of these three sibling species to evaluate the plausible demographical scenarios of the closely related species. We also tested the hypothesis that macro- and microevolutionary processes of the three species had been triggered and molded by Miocene–Pliocene geological events and Quaternary climatic change. The Bayesian cluster analysis showed that Q. acutissima and Q. chenii were clustered in the same group, whereas Q. variabilis formed a different genetic cluster. An approximate Bayesian computation (ABC) analyses suggested that Q. variabilis and Q. acutissima diverged from their most common ancestor around 19.84 Ma, and subsequently Q. chenii diverged from Q. acutissima at about 9.6 Ma, which was significantly associated with the episodes of the Qinghai–Tibetan Plateau (QTP). In addition, ecological niche modeling and population history analysis showed that these three Cerris oak species repeatedly underwent considerable ‘expansion–contraction’ during the interglacial and glacial periods of the Pleistocene, although they have varying degrees of tolerance for the climatic change. Overall, these findings indicated geological and climatic changes during the Miocene–Pliocene and Pleistocene as causes of species divergence and range shifts of dominant tree species in the subtropical and warm temperature areas in China.

Genetic signatures of divergent selection in European beech (Fagus sylvatica L.) are associated with the variation in temperature and precipitation across its distribution range

High genetic variation and extensive gene flow may help forest trees with adapting to ongoing climate change, yet the genetic bases underlying their adaptive potential remain largely unknown. We investigated range-wide patterns of potentially adaptive genetic variation in 64 populations of European beech (Fagus sylvatica L.) using 270 SNPs from 139 candidate genes involved either in phenology or in stress responses. We inferred neutral genetic structure and processes (drift and gene flow) and performed differentiation outlier analyses and gene-environment association (GEA) analyses to detect signatures of divergent selection. Beech range-wide genetic structure was consistent with the species' previously identified postglacial expansion scenario and recolonization routes. Populations showed high diversity and low differentiation along the major expansion routes. A total of 52 loci were found to be putatively under selection and 15 of them turned up in multiple GEA analyses. Temperature and precipitation related variables were equally represented in significant genotype-climate associations. Signatures of divergent selection were detected in the same proportion for stress response and phenology-related genes. The range-wide adaptive genetic structure of beech appears highly integrated, suggesting a balanced contribution of phenology and stress-related genes to local adaptation, and of temperature and precipitation regimes to genetic clines. Our results imply a best-case scenario for the maintenance of high genetic diversity during range shifts in beech (and putatively other forest trees) with a combination of gene flow maintaining within-population neutral diversity and selection maintaining between-population adaptive differentiation.

Genetic diversity and genetic differentiation of Megalobrama populations inferred by mitochondrial markers

BACKGROUND

Megalobrama is economically one of the most important freshwater fish genera in China. However, phylogenetic relationships among M. amblycephala, M. skolkovii, M. hoffmanni and M. pellegrini remain unresolved.

OBJECTIVE

To explore the genetic diversity and phylogenetic relationship of Megalobrama populations belonging to all four species.

METHODS

The concatenated sequences of mitochondrial cytochrome b (Cytb) and control region (CR) were used to analyze the genetic variation, genetic differentiation and population expansion of 15 Megalobrama populations.

RESULTS

The study showed that haplotype diversity and nucleotide diversity of M. hoffmanni and M. skolkovii were high, and that M. hoffmanni was the most genetically divergent of the four species. Haplotype network analysis revealed that M. hoffmanni and M. amblycephala formed a monophyletic group each, while M. skolkovii and M. pellegrini clustered together. There was a high genetic differentiation among the four Megalobrama species, and genetic distance among populations was not affected by geographical distance. Additionally, the results indicated that there was gene flow between the Liangzi Lake (LZL) population and Jinsha River Reservoir (JS) population. Also, Zhaoqing (ZQ) population of M. hoffmanni might have experienced a population expansion.

CONCLUSION

Our study verifies genetic diversity and differentiation of Megalobrama populations, and these findings will represent a significant contribution to the conservation and utilization of germplasm resources of Megalobrama.

How is adaptive potential distributed within species ranges?

Quantitative genetic variation (QGV) represents a major component of adaptive potential and, if reduced toward range-edge populations, could prevent a species' expansion or adaptive response to rapid ecological change. It has been hypothesized that QGV will be lower at the range edge due to small populations-often the result of poor habitat quality-and potentially decreased gene flow. However, whether central populations are higher in QGV is unknown. We used a meta-analytic approach to test for a general QGV-range position relationship, including geographic and climatic distance from range centers. We identified 35 studies meeting our criteria, yielding nearly 1000 estimates of QGV (including broad-sense heritability, narrow-sense heritability, and evolvability) from 34 species. The relationship between QGV and distance from the geographic range or climatic niche center depended on the focal trait and how QGV was estimated. We found some evidence that QGV declines from geographic centers but that it increases toward niche edges; niche and geographic distances were uncorrelated. Nevertheless, few studies have compared QGV in both central and marginal regions or environments within the same species. We call for more research in this area and discuss potential research avenues related to adaptive potential in the context of global change.

Niche differentiation in a postglacial colonizer, the bank vole Clethrionomys glareolus

Species-level environmental niche modeling has been crucial in efforts to understand how species respond to climate variation and change. However, species often exhibit local adaptation and intraspecific niche differences that may be important to consider in predicting responses to climate. Here, we explore whether phylogeographic lineages of the bank vole originating from different glacial refugia (Carpathian, Western, Eastern, and Southern) show niche differentiation, which would suggest a role for local adaptation in biogeography of this widespread Eurasian small mammal. We first model the environmental requirements for the bank vole using species-wide occurrences (210 filtered records) and then model each lineage separately to examine niche overlap and test for niche differentiation in geographic and environmental space. We then use the models to estimate past [Last Glacial Maximum (LGM) and mid-Holocene] habitat suitability to compare with previously hypothesized glacial refugia for this species. Environmental niches are statistically significantly different from each other for all pairs of lineages in geographic and environmental space, and these differences cannot be explained by habitat availability within their respective ranges. Together with the inability of most of the lineages to correctly predict the distributions of other lineages, these results support intraspecific ecological differentiation in the bank vole. Model projections of habitat suitability during the LGM support glacial survival of the bank vole in the Mediterranean region and in central and western Europe. Niche differences between lineages and the resulting spatial segregation of habitat suitability suggest ecological differentiation has played a role in determining the present phylogeographic patterns in the bank vole. Our study illustrates that models pooling lineages within a species may obscure the potential for different responses to climate change among populations.

Adaptation across geographic ranges is consistent with strong selection in marginal climates and legacies of range expansion

Every species experiences limits to its geographic distribution. Some evolutionary models predict that populations at range edges are less well adapted to their local environments due to drift, expansion load, or swamping gene flow from the range interior. Alternatively, populations near range edges might be uniquely adapted to marginal environments. In this study, we use a database of transplant studies that quantify performance at broad geographic scales to test how local adaptation, site quality, and population quality change from spatial and climatic range centers toward edges. We find that populations from poleward edges perform relatively poorly, both on average across all sites (15% lower population quality) and when compared to other populations at home (31% relative fitness disadvantage), consistent with these populations harboring high genetic load. Populations from equatorial edges also perform poorly on average (18% lower population quality) but, in contrast, outperform foreign populations (16% relative fitness advantage), suggesting that populations from equatorial edges have strongly adapted to unique environments. Finally, we find that populations from sites that are thermally extreme relative to the species' niche demonstrate strong local adaptation, regardless of their geographic position. Our findings indicate that both nonadaptive processes and adaptive evolution contribute to variation in adaptation across species' ranges.

Adaptation to drought is coupled with slow growth, but independent from phenology in marginal silver fir (Abies alba Mill.) populations

Drought is one of the most important selection pressures for forest trees in the context of climate change. Yet, the different evolutionary mechanisms, and their environmental drivers, by which certain populations become more drought tolerant than others is still little understood. We studied adaptation to drought in 16 silver fir (Abies alba Mill.) populations from the French Mediterranean Alps by combining observations on seedlings from a greenhouse experiment (N = 8,199) and on adult tress in situ (N = 315). In the greenhouse, we followed half-sib families for four growing seasons for growth and phenology traits, and tested their water stress response in a "drought until death" experiment. Adult trees in the field were assessed for δ 13C, a proxy for water use efficiency, and genotyped at 357 SNP loci. SNP data was used to generate a null expectation for seedling trait divergence between populations in order to detect the signature of selection, and 31 environmental variables were used to identify the selective environment. We found that seedlings originating from populations with low soil water capacity grew more slowly, attained a smaller stature, and resisted water stress for a longer period of time in the greenhouse. Additionally, adult trees of these populations exhibited a higher water use efficiency as evidenced by their δ 13C. These results suggest a correlated evolution of the growth-drought tolerance trait complex. Population divergence in bud break phenology was adaptive only in the second growing season, and evolved independently from the growth-drought tolerance trait complex. Adaptive divergence in bud break phenology was principally driven by the inter- and intra-annual variation in temperature at the geographic origin of the population. Our results illustrate the different evolutionary strategies used by populations to cope with drought stress at the range limits across a highly heterogeneous landscape, and can be used to inform assisted migration programs.

Contrasted patterns of local adaptation to climate change across the range of an evergreen oak, Quercus aquifolioides

Long-lived tree species are genetically differentiated and locally adapted with respect to fitness-related traits, but the genetic basis of local adaptation remains largely unresolved. Recent advances in population genetics and landscape genomic analyses enable identification of putative adaptive loci and specific selective pressures acting on local adaptation. Here, we sampled 60 evergreen oak (Quercus aquifolioides) populations throughout the species' range and pool-sequenced 587 individuals at drought-stress candidate genes. We analyzed patterns of genetic diversity and differentiation for 381 single nucleotide polymorphisms (SNPs) from 65 candidate genes and eight microsatellites. Outlier loci were identified by genetic differentiation analysis and genome-environment associations. The response pattern of genetic variation to environmental gradient was assessed by linear isolation-by-distance/environment tests, redundancy analysis, and nonlinear methods. SNPs and microsatellites revealed two genetic lineages: Tibet and Hengduan Mountains-Western Sichuan Plateau (HDM-WSP), with reduced genetic diversity in Tibet lineage. More outlier loci were detected in HDM-WSP lineage than Tibet lineage. Among these, three SNPs in two genes responded to dry season precipitation in the HDM-WSP lineage but not in Tibet. By contrast, genetic variation in the Tibet lineage was related to geographic distance instead of the environment. Furthermore, risk of nonadaptedness (RONA) analyses suggested HDM-WSP lineage will have a better capacity to adapt in the predicted future climate compared with the Tibet lineage. We detected genetic imprints consistent with natural selection and molecular adaptation to drought on the Qinghai-Tibet Plateau (QTP) over a range of long-lived and widely distributed oak species in a changing environment. Our results suggest that different within-species adaptation processes occur in species occurring in heterogeneous environments.

Exploring genomic variation associated with drought stress in Picea mariana populations

Predicted increases in drought and heat stress will likely induce shifts in species bioclimatic envelopes. Genetic variants adapted to water limitation may prove pivotal for species response under scenarios of increasing drought. In this study, we aimed to explore this hypothesis by investigating genetic variation in 16 populations of black spruce (Picea mariana) in relation to climate variables in Alaska. A total of 520 single nucleotide polymorphisms (SNPs) were genotyped for 158 trees sampled from areas of contrasting climate regimes. We used multivariate and univariate genotype-by-environment approaches along with available gene annotations to investigate the relationship between climate and genetic variation among sampled populations. Nine SNPs were identified as having a significant association with climate, of which five were related to drought stress response. Outlier SNPs with respect to the overall environment were significantly overrepresented for several biological functions relevant for coping with variable hydric regimes, including osmotic stress response. This genomic imprint is consistent with local adaptation of black spruce to drought stress. These results suggest that natural selection acting on standing variation prompts local adaptation in forest stands facing water limitation. Improved understanding of possible adaptive responses could inform our projections about future forest dynamics and help prioritize populations that harbor valuable genetic diversity for conservation.

Antimicrobial peptide and sequence variation along a latitudinal gradient in two anurans

BACKGROUND

While there is evidence of both purifying and balancing selection in immune defense genes, large-scale genetic diversity in antimicrobial peptides (AMPs), an important part of the innate immune system released from dermal glands in the skin, has remained uninvestigated. Here we describe genetic diversity at three AMP loci (Temporin, Brevinin and Palustrin) in two ranid frogs (Rana arvalis and R. temporaria) along a 2000 km latitudinal gradient. We amplified and sequenced part of the Acidic Propiece domain and the hypervariable Mature Peptide domain (~ 150-200 bp) in the three genes using Illumina Miseq and expected to find decreased AMP genetic variation towards the northern distribution limit of the species similarly to studies on MHC genetic patterns.

RESULTS

We found multiple loci for each AMP and relatively high gene diversity, but no clear pattern of geographic genetic structure along the latitudinal gradient. We found evidence of trans-specific polymorphism in the two species, indicating a common evolutionary origin of the alleles. Temporin and Brevinin did not form monophyletic clades suggesting that they belong to the same gene family. By implementing codon evolution models we found evidence of strong positive selection acting on the Mature Peptide. We also found evidence of diversifying selection as indicated by divergent allele frequencies among populations and high Theta k values.

CONCLUSION

Our results suggest that AMPs are an important source of adaptive diversity, minimizing the chance of microorganisms developing resistance to individual peptides.

Ice-age persistence and genetic isolation of the disjunct distribution of larch in Alaska

Larix laricina (eastern larch, tamarack) is a transcontinental North American conifer with a prominent disjunction in the Yukon isolating the Alaskan distribution from the rest of its range. We investigate whether in situ persistence during the last glacial maximum (LGM) or long-distance postglacial migration from south of the ice sheets resulted in the modern-day Alaskan distribution. We analyzed variation in three chloroplast DNA regions of 840 trees from a total of 69 populations (24 new sampling sites situated on both sides of the Yukon range disjunction pooled with 45 populations from a published source) and conducted ensemble species distribution modeling (SDM) throughout Canada and United States to hindcast the potential range of L. laricina during the LGM. We uncovered the genetic signature of a long-term isolation of larch populations in Alaska, identifying three endemic chlorotypes and low levels of genetic diversity. Range-wide analysis across North America revealed the presence of a distinct Alaskan lineage. Postglacial gene flow across the Yukon divide was unidirectional, from Alaska toward previously glaciated Canadian regions, and with no evidence of immigration into Alaska. Hindcast SDM indicates one of the broadest areas of past climate suitability for L. laricina existed in central Alaska, suggesting possible in situ persistence of larch in Alaska during the LGM. Our results provide the first unambiguous evidence for the long-term isolation of L. laricina in Alaska that extends beyond the last glacial period and into the present interglacial period. The lack of gene flow into Alaska along with the overall probability of larch occurrence in Alaska being currently lower than during the LGM suggests that modern-day Alaskan larch populations are isolated climate relicts of broader glacial distributions, and so are particularly vulnerable to current warming trends.

The Evolution of Paleoecology

While the interplay between migration and adaptation dictates species response to climate change, technological limitations have obfuscated explicit tests on past adaptive responses. However, a surge in technology-driven advances in paleoecological methods coincides with breakthroughs in processing ancient DNA, providing the first opportunity to assess adaptation to past climate shifts.

A multiscale approach to detect selection in nonmodel tree species: Widespread adaptation despite population decline in Taxus baccata L

Detecting the molecular basis of local adaptation and identifying selective drivers is still challenging in nonmodel species. The use of purely population genetic approaches is limited by some characteristics of genetic systems, such as pleiotropy and polygenic control, and parallel evidence from phenotypic-based experimental comparisons is required. In long-lived organisms, the detection of selective pressures might also be precluded by evolutionary lag times in response to the environment. Here, we used the English yew to showcase an example of a multiscale integrative approach in a nonmodel species with limited plant and genomic resources. We combined information from two independent sources, phenotypes in a common environment and genomic data in natural populations, to investigate the signature of selection. Growth differences among populations in a common environment, and phenological patterns of both shoot elongation and male strobili maturation, were associated with climate clines, providing evidence for local adaptation and guiding us in the selection of populations for genomic analyses. We used information on over 25,000 SNPs from c. 1,200 genes to infer the demographic history and to test for molecular signatures of selection at different levels: SNP, gene, and biological pathway. Our results confirmed an overall demographic history of population decline, but we also found evidence for putative local adaptation at the molecular level. We identified or confirmed several candidate genes for positive and negative selection in forest trees, including the pseudo-response regulator 7 (PRR7), an essential component of the circadian clock in plants. In addition, we successfully tested an approach to detect polygenic adaptation in biological pathways, allowing us to identify the flavonoid biosynthesis pathway as a candidate stress-response pathway that deserves further attention in other plants. Finally, our study contributes to the emerging view that explaining contemporary standing genetic variation requires considering adaptation to past climates, especially for long-lived trees.

Spatial genetic structure and diversity of natural populations of Aesculus hippocastanum L. in Greece

Horse-chestnut (Aesculus hippocastanum L.) is an endemic and relict species from the Mediterranean biodiversity hotspot and a popular ornamental tree. Knowledge about the evolutionary history of this species remains scarce. Here, we ask what historical and ecological factors shaped the pattern of genetic diversity and differentiation of this species. We genotyped 717 individuals from nine natural populations using microsatellite markers. The influence of distance, topography and habitat variables on spatial genetic structure was tested within the approaches of isolation-by-distance and isolation-by-ecology. Species niche modeling was used to project the species theoretical range through time and space. The species showed high genetic diversity and moderate differentiation for which topography, progressive range contraction through the species’ history and long-term persistence in stable climatic refugia are likely responsible. A strong geographic component was revealed among five genetic clusters that are connected with very limited gene flow. The environmental variables were a significant factor in the spatial genetic structure. Modeling results indicated that future reduction of the species range may affect its survival. The possible impact of climate changes and high need of in situ conservation are discussed.

Terpene Synthase Genes in Quercus robur - Gene Characterization, Expression and Resulting Terpenes Due to Cockchafer Feeding

Root herbivory caused by larvae of the forest cockchafer (Melolontha hippocastani) enhances the impact of drought on trees, particularly in oak forest rejuvenations. In Germany, geographically distant oak stands show differences in infestation strength by the forest cockchafer. While in Southwestern Germany this insect causes severe damage, oak forests in northern Germany are rarely infested. It is known that root-released volatile organic compounds (VOCs) are perceived by soil herbivores, thus guiding the larvae toward the host roots. In this work, we exposed seedlings of two distant oak provenances to forest cockchafer larvae and studied their population genetic properties, their root-based VOC chemotypes, their attraction for larvae and terpene synthase gene expression. Based on nuclear and chloroplast marker analysis, we found both oak populations to be genetically highly variable while showing typical patterns of migration from different refugial regions. However, no clear association between genetic constitution of the different provenances and the abundance of cockchafer populations on site was observed. In contrast to observations in the field, bioassays revealed a preference of the larvae for the northeastern oak provenance. The behavior of larvae was most likely related to root-released volatile terpenes and benzenoids since their composition and quantity differed between oak populations. We assume repellent effects of these compounds because the populations attractive to insects showed low abundance of these compounds. Five different oak terpene synthase (TPS) genes were identified at the genomic level which can be responsible for biosynthesis of the released terpenes. TPS gene expression patterns in response to larval feeding revealed geographic variation rather than genotypic variation. Our results support the assumption that root-released VOC are influencing the perception of roots by herbivores.

Playing Hide-and-Seek in Beta-Globin Genes: Gene Conversion Transferring a Beneficial Mutation between Differentially Expressed Gene Duplicates

Increasing evidence suggests that adaptation to diverse environments often involves selection on existing variation rather than new mutations. A previous study identified a nonsynonymous single nucleotide polymorphism (SNP) in exon 2 of two paralogous β-globin genes of the bank vole (Clethrionomys glareolus) in Britain in which the ancestral serine (Ser) and the derived cysteine (Cys) allele represent geographically partitioned functional variation affecting the erythrocyte antioxidative capacity. Here we studied the geographical pattern of the two-locus Ser/Cys polymorphism throughout Europe and tested for the geographic correlation between environmental variables and allele frequency, expected if the polymorphism was under spatially heterogeneous environment-related selection. Although bank vole population history clearly is important in shaping the dispersal of the oxidative stress protective Cys allele, analyses correcting for population structure suggest the Europe-wide pattern is affected by geographical variation in environmental conditions. The β-globin phenotype is encoded by the major paralog HBB-T1 but we found evidence of bidirectional gene conversion of exon 2 with the low-expression paralog HBB-T2. Our data support the model where gene conversion reshuffling genotypes between high- and low- expressed paralogs enables tuning of erythrocyte thiol levels, which may help maintain intracellular redox balance under fluctuating environmental conditions. Therefore, our study suggests a possible role for gene conversion between differentially expressed gene duplicates as a mechanism of physiological adaptation of populations to new or changing environments.