Genomics of local adaptation with gene flow

Ecological Interactions, Environmental Gradients, and Gene Flow in Local Adaptation

Despite long-standing interest in local adaptation of plants to their biotic and abiotic environment, existing theory, and many case studies, little work to date has addressed within-species evolution of concerted strategies and how these might contrast with patterns across species. Here we consider the interactions between pollinators, herbivores, and resource availability in shaping plant local adaptation, how these interactions impact plant phenotypes and gene flow, and the conditions where multiple traits align along major environmental gradients such as latitude and elevation. Continued work in emerging model systems will benefit from the melding of classic experimental approaches with novel population genetic analyses to reveal patterns and processes in plant local adaptation.

The role of neutral and adaptive genomic variation in population diversification and speciation in two ground squirrel species of conservation concern

Understanding the neutral (demographic) and adaptive processes leading to the differentiation of species and populations is a critical component of evolutionary and conservation biology. In this context, recently diverged taxa represent a unique opportunity to study the process of genetic differentiation. Northern and southern Idaho ground squirrels (Urocitellus brunneus - NIDGS, and U. endemicus - SIDGS, respectively) are a recently diverged pair of sister species that have undergone dramatic declines in the last 50 years and are currently found in metapopulations across restricted spatial areas with distinct environmental pressures. Here we genotyped single-nucleotide polymorphisms (SNPs) from buccal swabs with restriction site-associated DNA sequencing (RADseq). With these data we evaluated neutral genetic structure at both theinter- and intraspecific level, and identified putatively adaptive SNPs using population structure outlier detection and genotype-environment association (GEA) analyses. At the interspecific level, we detected a clear separation between NIDGS and SIDGS, and evidence for adaptive differentiation putatively linked to torpor patterns. At the intraspecific level, we found evidence of both neutral and adaptive differentiation. For NIDGS, elevation appears to be the main driver of adaptive differentiation, while neutral variation patterns match and expand information on the low connectivity between some populations identified in previous studies using microsatellite markers. For SIDGS, neutral substructure generally reflected natural geographic barriers, while adaptive variation reflected differences in land cover and temperature, as well as elevation. These results clearly highlight the roles of neutral and adaptive processes for understanding the complexity of the processes leading to species and population differentiation, which can have important conservation implications in susceptible and threatened species.

Multidimensional divergent selection, local adaptation, and speciation

Divergent selection applied to one or more traits drives local adaptation and may lead to ecological speciation. Divergent selection on many traits might be termed "multidimensional" divergent selection. There is a commonly held view that multidimensional divergent selection is likely to promote local adaptation and speciation to a greater extent than unidimensional divergent selection. We disentangle the core concepts underlying dimensionality as a property of the environment, phenotypes, and genome. In particular, we identify a need to separate the overall strength of selection and the number of loci affected from dimensionality per se, and to distinguish divergence dimensionality from dimensionality of stabilizing selection. We then critically scrutinize this commonly held view that multidimensional selection promotes speciation, re-examining the evidence base from theory, experiments, and nature. We conclude that the evidence base is currently weak and generally suffers from confounding of possible causal effects. Finally, we propose several mechanisms by which multidimensional divergent selection and related processes might influence divergence, both as a driver and as a barrier.

Genomics of altitude-associated wing shape in two tropical butterflies

Understanding how organisms adapt to their local environment is central to evolution. With new whole-genome sequencing technologies and the explosion of data, deciphering the genomic basis of complex traits that are ecologically relevant is becoming increasingly feasible. Here, we studied the genomic basis of wing shape in two Neotropical butterflies that inhabit large geographical ranges. Heliconius butterflies at high elevations have been shown to generally have rounder wings than those in the lowlands. We reared over 1,100 butterflies from 71 broods of H. erato and H. melpomene in common-garden conditions and showed that wing aspect ratio, that is, elongatedness, is highly heritable in both species and that elevation-associated wing aspect ratio differences are maintained. Genome-wide associations with a published data set of 666 whole genomes from across a hybrid zone, uncovered a highly polygenic basis to wing aspect ratio variation in the wild. We identified several genes that have roles in wing morphogenesis or wing aspect ratio variation in Drosophila flies, making them promising candidates for future studies. There was little evidence for molecular parallelism in the two species, with only one shared candidate gene, nor for a role of the four known colour pattern loci, except for optix in H. erato. Thus, we present the first insights into the heritability and genomic basis of within-species wing aspect ratio in two Heliconius species, adding to a growing body of evidence that polygenic adaptation may underlie many ecologically relevant traits.

Geographic patterns of seed trait variation in an invasive species: how much can close populations differ?

Seeds play a major role in plant species persistence and expansion, and therefore they are essential when modeling species dynamics. However, homogeneity in seed traits is generally assumed, underestimating intraspecific trait variability across the geographic space, which might bias species success models. The aim of this study was to evaluate the existence and consequences of interpopulation variability in seed traits of the invasive species Carpobrotus edulis at different geographical scales. We measured seed production, morphology, vigour and longevity of nine populations of C. edulis along the Catalan coast (NE Spain) from three differentiated zones with a human presence gradient. Geographic distances between populations were contrasted against individual and multivariate trait distances to explore trait variation along the territory, evaluating the role of bioclimatic variables and human density of the different zones. The analysis revealed high interpopulation variability that was not explained by geographic distance, as regardless of the little distance between some populations (< 0.5 km), significant differences were found in several seed traits. Seed production, germination, and persistence traits showed the strongest spatial variability up to 6000% of percent trait variability between populations, leading to differentiated C. edulis soil seed bank dynamics at small distances, which may demand differentiated strategies for a cost-effective species management. Seed trait variability was influenced by human density but also bioclimatic conditions, suggesting a potential impact of increased anthropogenic pressure and climate shifts. Geographic interpopulation trait variation should be included in ecological models and will be important for assessing species responses to environmental heterogeneity and change.

Genomic Adaptations to Salinity Resist Gene Flow in the Evolution of Floridian Watersnakes

The migration-selection balance often governs the evolution of lineages, and speciation with gene flow is now considered common across the tree of life. Ecological speciation is a process that can facilitate divergence despite gene flow due to strong selective pressures caused by ecological differences; however, the exact traits under selection are often unknown. The transition from freshwater to saltwater habitats provides strong selection targeting traits with osmoregulatory function. Several lineages of North American watersnakes (Nerodia spp.) are known to occur in saltwater habitat and represent a useful system for studying speciation by providing an opportunity to investigate gene flow and evaluate how species boundaries are maintained or degraded. We use double digest restriction-site associated DNA sequencing to characterize the migration-selection balance and test for evidence of ecological divergence within the Nerodia fasciata-clarkii complex in Florida. We find evidence of high intraspecific gene flow with a pattern of isolation-by-distance underlying subspecific lineages. However, we identify genetic structure indicative of reduced gene flow between inland and coastal lineages suggesting divergence due to isolation-by-environment. This pattern is consistent with observed environmental differences where the amount of admixture decreases with increased salinity. Furthermore, we identify significantly enriched terms related to osmoregulatory function among a set of candidate loci, including several genes that have been previously implicated in adaptation to salinity stress. Collectively, our results demonstrate that ecological differences, likely driven by salinity, cause strong divergent selection which promotes divergence in the N. fasciata-clarkii complex despite significant gene flow.

Adaptive divergence despite low effective population size in a peripherally isolated population of the pygmy rabbit, Brachylagus idahoensis

Local adaptation can occur when spatially separated populations are subjected to contrasting environmental conditions. Historically, understanding the genetic basis of adaptation has been difficult, but increased availability of genome-wide markers facilitates studies of local adaptation in non-model organisms of conservation concern. The pygmy rabbit (Brachylagus idahoensis) is an imperiled lagomorph that relies on sagebrush for forage and cover. This reliance has led to widespread population declines following reductions in the distribution of sagebrush, leading to geographic separation between populations. In this study, we used >20,000 single nucleotide polymorphisms, genotype-environment association methods, and demographic modeling to examine neutral genetic variation and local adaptation in the pygmy rabbit in Nevada and California. We identified 308 loci as outliers, many of which had functional annotations related to metabolism of plant secondary compounds. Likewise, patterns of spatial variation in outlier loci were correlated with landscape and climatic variables including proximity to streams, sagebrush cover, and precipitation. We found that populations in the Mono Basin of California probably diverged from other Great Basin populations during late Pleistocene climate oscillations, and that this region is adaptively differentiated from other regions in the southern Great Basin despite limited gene flow and low effective population size. Our results demonstrate that peripherally isolated populations can maintain adaptive divergence.

Hybridization and introgression in sympatric and allopatric populations of four oak species

BACKGROUND

Hybridization and introgression are vital sources of novel genetic variation driving diversification during reticulated evolution. Quercus is an important model clade, having extraordinary diverse and abundant members in the Northern hemisphere, that are used to studying the introgression of species boundaries and adaptive processes. China is the second-largest distribution center of Quercus, but there are limited studies on introgressive hybridization.

RESULTS

Here, we screened 17 co-dominant nuclear microsatellite markers to investigate the hybridization and introgression of four oaks (Quercus acutissima, Quercus variabilis, Quercus fabri, and Quercus serrata) in 10 populations. We identified 361 alleles in the four-oak species across 17 loci, and all loci were characterized by high genetic variability (HE = 0.844-0.944) and moderate differentiation (FST = 0.037-0.156) levels. A population differentiation analysis revealed the following: allopatric homologous (FST = 0.064) < sympatric heterogeneous (FST = 0.071) < allopatric heterogeneous (FST = 0.084). A Bayesian admixture analysis determined four types of hybrids (Q. acutissima × Q. variabilis, Q. fabri × Q. serrata, Q. acutissima × Q. fabri, and Q. acutissima × Q. variabilis × Q. fabri) and their asymmetric introgression. Our results revealed that interspecific hybridization is commonly observed within the section Quercus, with members having tendency to hybridize.

CONCLUSIONS

Our study determined the basic hybridization and introgression states among the studied four oak species and extended our understanding of the evolutionary role of hybridization. The results provide useful theoretical data for formulating conservation strategies.

Chromosome-Level Assembly of the Atlantic Silverside Genome Reveals Extreme Levels of Sequence Diversity and Structural Genetic Variation

The levels and distribution of standing genetic variation in a genome can provide a wealth of insights about the adaptive potential, demographic history, and genome structure of a population or species. As structural variants are increasingly associated with traits important for adaptation and speciation, investigating both sequence and structural variation is essential for wholly tapping this potential. Using a combination of shotgun sequencing, 10x Genomics linked reads and proximity-ligation data (Chicago and Hi-C), we produced and annotated a chromosome-level genome assembly for the Atlantic silverside (Menidia menidia)-an established ecological model for studying the phenotypic effects of natural and artificial selection-and examined patterns of genomic variation across two individuals sampled from different populations with divergent local adaptations. Levels of diversity varied substantially across each chromosome, consistently being highly elevated near the ends (presumably near telomeric regions) and dipping to near zero around putative centromeres. Overall, our estimate of the genome-wide average heterozygosity in the Atlantic silverside is among the highest reported for a fish, or any vertebrate (1.32-1.76% depending on inference method and sample). Furthermore, we also found extreme levels of structural variation, affecting ∼23% of the total genome sequence, including multiple large inversions (> 1 Mb and up to 12.6 Mb) associated with previously identified haploblocks showing strong differentiation between locally adapted populations. These extreme levels of standing genetic variation are likely associated with large effective population sizes and may help explain the remarkable adaptive divergence among populations of the Atlantic silverside.

Limited Evidence for Parallel Evolution Among Desert-Adapted Peromyscus Deer Mice

Warming climate and increasing desertification urge the identification of genes involved in heat and dehydration tolerance to better inform and target biodiversity conservation efforts. Comparisons among extant desert-adapted species can highlight parallel or convergent patterns of genome evolution through the identification of shared signatures of selection. We generate a chromosome-level genome assembly for the canyon mouse (Peromyscus crinitus) and test for a signature of parallel evolution by comparing signatures of selective sweeps across population-level genomic resequencing data from another congeneric desert specialist (Peromyscus eremicus) and a widely distributed habitat generalist (Peromyscus maniculatus), that may be locally adapted to arid conditions. We identify few shared candidate loci involved in desert adaptation and do not find support for a shared pattern of parallel evolution. Instead, we hypothesize divergent molecular mechanisms of desert adaptation among deer mice, potentially tied to species-specific historical demography, which may limit or enhance adaptation. We identify a number of candidate loci experiencing selective sweeps in the P. crinitus genome that are implicated in osmoregulation (Trypsin, Prostasin) and metabolic tuning (Kallikrein, eIF2-alpha kinase GCN2, APPL1/2), which may be important for accommodating hot and dry environmental conditions.

Topography shapes the local coexistence of tree species within species complexes of Neotropical forests

Forest inventories in Amazonia include around 5000 described tree species belonging to more than 800 genera. Numerous species-rich genera share genetic variation among species because of recent speciation and/or recurrent hybridisation, forming species complexes. Despite the key role that tree species complexes play in understanding Neotropical diversification, and their need to exploit a diversity of niches, little is known about the mechanisms that allow local coexistence of tree species complexes and their species in sympatry. In this study, we explored the fine-scale distribution of five tree species complexes and 22 species within these complexes. Combining forest inventories, botanical determination, and LiDAR-derived topographic data over 120 ha of permanent plots in French Guiana, we used a Bayesian modelling framework to test the role of fine-scale topographic wetness and tree neighbourhood on the occurrence of species complexes and the relative distribution of species within complexes. Species complexes of Neotropical trees were widely spread across the topographic wetness gradient at the local scale. Species within complexes showed pervasive niche differentiation along with topographic wetness and competition gradients. Similar patterns of species-specific habitat preferences were observed within several species complexes: species more tolerant to competition for resources grow in drier and less fertile plateaus and slopes. If supported by partial reproductive isolation of species and adaptive introgression at the species complex level, our results suggest that both species-specific habitat specialisation within species complexes and the broad ecological distribution of species complexes might explain the success of these species complexes at the regional scale.

The maintenance of polygenic sex determination depends on the dominance of fitness effects which are predictive of the role of sexual antagonism

In species with polygenic sex determination (PSD), multiple male- and female-determining loci on different proto-sex chromosomes segregate as polymorphisms within populations. The extent to which these polymorphisms are at stable equilibria is not yet resolved. Previous work demonstrated that PSD is most likely to be maintained as a stable polymorphism when the proto-sex chromosomes have opposite (sexually antagonistic) fitness effects in males and females. However, these models usually consider PSD systems with only two proto-sex chromosomes, or they do not broadly consider the dominance of the alleles under selection. To address these shortcomings, I used forward population genetic simulations to identify selection pressures that can maintain PSD under different dominance scenarios in a system with more than two proto-sex chromosomes (modeled after the house fly). I found that overdominant fitness effects of male-determining proto-Y chromosomes are more likely to maintain PSD than dominant, recessive, or additive fitness effects. The overdominant fitness effects that maintain PSD tend to have proto-Y chromosomes with sexually antagonistic effects (male-beneficial and female-detrimental). In contrast, dominant fitness effects that maintain PSD tend to have sexually antagonistic multi-chromosomal genotypes, but the individual proto-sex chromosomes do not have sexually antagonistic effects. These results demonstrate that sexual antagonism can be an emergent property of the multi-chromosome genotype without individual sexually antagonistic chromosomes. My results further illustrate how the dominance of fitness effects has consequences for both the likelihood that PSD will be maintained as well as the role sexually antagonistic selection is expected to play in maintaining the polymorphism.

Multiple paths to the same destination: Influence of gene flow on convergent evolution

Elucidation of the genetic mechanisms of convergent evolution, the evolution of similar or the same phenotypes in phylogenetically independent lineages, helps predict how populations will respond to the same selective pressures. Convergent evolution can be caused by either the fixation of identical-by-descent alleles, independent mutations at the same gene, or mutations in different genes controlling the same trait. To what extent does the fixation of identical-by-descent alleles lead to convergent evolution in isolated populations where inflow of adaptive alleles from other populations is limited? In a From the Cover article in this issue of Molecular Ecology, Kemppainen et al. (2021) compared the genetic basis for the reduction of pelvic structures in three isolated freshwater populations of nine-spined stickleback (Pungitius pungitius) from Northern Europe. The authors used quantitative trait loci (QTL) mapping to reveal that the pelvic reduction in these three populations was caused by mutations at different genetic loci. In contrast to studies in three-spined stickleback (Gasterosteus aculeatus), where independently derived Pitx1 mutations were shown to be responsible for plate reduction across multiple freshwater populations, Kemppainen et al. (2021) found Pitx1 to be the candidate causative gene for only one population of P. pungitius. This study highlights the importance of genetic studies of convergent evolution, not only in the presence of gene flow but also in its absence for a better understanding of the genetic architecture of convergent evolution.

The river shapes the genetic diversity of common reed in the Yellow River Delta via hydrochory dispersal and habitat selection

Understanding the driving mechanisms of local genetic diversity is a fundamental challenge under the global environmental changes. Rivers provide an excellent study system to demonstrate the effects of hydrochory dispersal and habitat selection on genetic diversity of riparian flora. In this study, we focused on the genetic variation of common reed (Phragmites australis) in the Yellow River Delta, China. Firstly, samples were collected in the Yellow River Delta, its neighboring wetland and its upstream plain. The genetic variation of P. australis was investigated using two chloroplast DNA fragments and eleven nuclear microsatellites. The findings showed that the genetic variation of P. australis in the Yellow River Delta belonged to two distinct lineages (haplotype O and haplotype P), which were similar to the upstream, and to the neighboring populations, respectively. Moreover, the genetic results suggested the potential dispersal of haplotype O from upstream to downstream. Secondly, we surveyed the plant functional traits of common reed from the Yellow River Delta in the field and in the common garden. The results showed significant differences between riverine and non-riverine populations in plant functional traits (e.g. specific leaf area and leaf length), haplotype composition and genetic clustering, which implied natural selection by habitat conditions. Lastly, we re-analyzed the plant performance data from a salt manipulation experiment with different haplotypes, and the results supported that salinity is a significant selective stressor on P. australis lineages in the Yellow River Delta. Our study highlights the significance of hydrochory dispersal and habitat selection in the river effects on genetic diversity of riparian flora, and provides important information for biodiversity conservation and wetland management in the Yellow River Delta.

Stabilization and optimization of host-microbe-environment interactions as a potential reason for the behavior of natal philopatry

Many animals engage in a behavior known as natal philopatry, where after sexual maturity they return to their own birthplaces for subsequent reproduction. There are many proposed ultimate factors that may underlie the evolution of natal philopatry, such as genetic optimization, suitable living conditions, and friendly neighbors, which can improve the survival rates of offspring. However, here we propose that a key factor that has been overlooked could be the colonization of gut microbiota during early life and the effects these microorganisms have on host performance and fitness. In addition to the bacteria transmitted from the mother to offspring, microbes from the surrounding environment also account for a large proportion of the developing gut microbiome. While it was long believed that microbial species all have global distributions, we now know that there are substantial geographic differences and dispersal limitations to environmental microbes. The establishment of gut microbiota during early life has enormous impacts on animal development, including energy metabolism, training of the immune system, and cognitive development. Moreover, these microbial effects scale to influence animal performance and fitness, raising the possibility for natural selection to act on the integrated combination of gut microbial communities and host genetics (i.e. the holobiont). Therefore, in this paper, we propose a hypothesis: that optimization of host-microbe-environment interactions represents a potentially important yet overlooked reason for natal philopatry. Microbiota obtained by natal philopatry could help animals adapt to the environment and improve the survival rates of their young. We propose future directions to test these ideas, and the implications that this hypothesis has for our understanding of host-microbe interactions.

Adaptive Divergence under Gene Flow along an Environmental Gradient in Two Coexisting Stickleback Species

There is a general and solid theoretical framework to explain how the interplay between natural selection and gene flow affects local adaptation. Yet, to what extent coexisting closely related species evolve collectively or show distinctive evolutionary responses remains a fundamental question. To address this, we studied the population genetic structure and morphological differentiation of sympatric three-spined and nine-spined stickleback. We conducted genotyping-by-sequencing and morphological trait characterisation using 24 individuals of each species from four lowland brackish water (LBW), four lowland freshwater (LFW) and three upland freshwater (UFW) sites in Belgium and the Netherlands. This combination of sites allowed us to contrast populations from isolated but environmentally similar locations (LFW vs. UFW), isolated but environmentally heterogeneous locations (LBW vs. UFW), and well-connected but environmentally heterogenous locations (LBW vs. LFW). Overall, both species showed comparable levels of genetic diversity and neutral genetic differentiation. However, for all three spatial scales, signatures of morphological and genomic adaptive divergence were substantially stronger among populations of the three-spined stickleback than among populations of the nine-spined stickleback. Furthermore, most outlier SNPs in the two species were associated with local freshwater sites. The few outlier SNPs that were associated with the split between brackish water and freshwater populations were located on one linkage group in three-spined stickleback and two linkage groups in nine-spined stickleback. We conclude that while both species show congruent evolutionary and genomic patterns of divergent selection, both species differ in the magnitude of their response to selection regardless of the geographical and environmental context.

Persistent panmixia despite extreme habitat loss and population decline in the threatened tricolored blackbird (Agelaius tricolor)

Habitat loss and alteration has driven many species into decline, often to the point of requiring protection and intervention to avert extinction. Genomic data provide the opportunity to inform conservation and recovery efforts with details about vital evolutionary processes with a resolution far beyond that of traditional genetic approaches. The tricolored blackbird (Agelaius tricolor) has suffered severe losses during the previous century largely due to anthropogenic impacts on their habitat. Using a dataset composed of a whole genome paired with reduced representation libraries (RAD-Seq) from samples collected across the species' range, we find evidence for panmixia using multiple methods, including PCA (no geographic clustering), admixture analyses (ADMIXTURE and TESS conclude K = 1), and comparisons of genetic differentiation (average FST = 0.029). Demographic modeling approaches recovered an ancient decline that had a strong impact on genetic diversity but did not detect any effect from the known recent decline. We also did not detect any evidence for selection, and hence adaptive variation, at any site, either geographic or genomic. These results indicate that species continues to have high vagility across its range despite population decline and habitat loss and should be managed as a single unit.

Environmental correlates of genetic variation in the invasive European starling in North America

Populations of invasive species that colonize and spread in novel environments may differentiate both through demographic processes and local selection. European starlings (Sturnus vulgaris) were introduced to New York in 1890 and subsequently spread throughout North America, becoming one of the most widespread and numerous bird species on the continent. Genome-wide comparisons across starling individuals and populations can identify demographic and/or selective factors that facilitated this rapid and successful expansion. We investigated patterns of genomic diversity and differentiation using reduced-representation genome sequencing of 17 winter-season sampling sites. Consistent with this species' high dispersal rate and rapid expansion history, we found low geographical differentiation and few F_ST outliers even at a continental scale. Despite starting from a founding population of ~180 individuals, North American starlings show only a moderate genetic bottleneck, and models suggest a dramatic increase in effective population size since introduction. In genotype-environment associations we found that ~200 single-nucleotide polymorphisms are correlated with temperature and/or precipitation against a background of negligible genome- and range-wide divergence. Given this evidence, we suggest that local adaptation in North American starlings may have evolved rapidly even in this wide-ranging and evolutionarily young system. This survey of genomic signatures of expansion in North American starlings is the most comprehensive to date and complements ongoing studies of world-wide local adaptation in these highly dispersive and invasive birds.

Genomic variation, population history and within-archipelago adaptation between island bird populations

Oceanic island archipelagos provide excellent models to understand evolutionary processes. Colonization events and gene flow can interact with selection to shape genetic variation at different spatial scales. Landscape-scale variation in biotic and abiotic factors may drive fine-scale selection within islands, while long-term evolutionary processes may drive divergence between distantly related populations. Here, we examine patterns of population history and selection between recently diverged populations of the Berthelot's pipit (Anthus berthelotii), a passerine endemic to three North Atlantic archipelagos. First, we use demographic trees and f3 statistics to show that genome-wide divergence across the species range is largely shaped by colonization and bottlenecks, with evidence of very weak gene flow between populations. Then, using a genome scan approach, we identify signatures of divergent selection within archipelagos at single nucleotide polymorphisms (SNPs) in genes potentially associated with craniofacial development and DNA repair. We did not detect within-archipelago selection at the same SNPs as were detected previously at broader spatial scales between archipelagos, but did identify signatures of selection at loci associated with similar biological functions. These findings suggest that similar ecological factors may repeatedly drive selection between recently separated populations, as well as at broad spatial scales across varied landscapes.

Genome-wide diversity and habitat underlie fine-scale phenotypic differentiation in the rainbow darter (Etheostoma caeruleum)

Adaptation to environmental change requires that populations harbor the necessary genetic variation to respond to selection. However, dispersal-limited species with fragmented populations and reduced genetic diversity may lack this variation and are at an increased risk of local extinction. In freshwater fish species, environmental change in the form of increased stream temperatures places many cold-water species at-risk. We present a study of rainbow darters (Etheostoma caeruleum) in which we evaluated the importance of genetic variation on adaptive potential and determined responses to extreme thermal stress. We compared fine-scale patterns of morphological and thermal tolerance differentiation across eight sites, including a unique lake habitat. We also inferred contemporary population structure using genomic data and characterized the relationship between individual genetic diversity and stress tolerance. We found site-specific variation in thermal tolerance that generally matched local conditions and morphological differences associated with lake-stream divergence. We detected patterns of population structure on a highly local spatial scale that could not be explained by isolation by distance or stream connectivity. Finally, we showed that individual thermal tolerance was positively correlated with genetic variation, suggesting that sites with increased genetic diversity may be better at tolerating novel stress. Our results highlight the importance of considering intraspecific variation in understanding population vulnerability and stress response.

The Effects of Weather on Avian Growth and Implications for Adaptation to Climate Change

Climate change is forecasted to generate a range of evolutionary changes and plastic responses. One important aspect of avian responses to climate change is how weather conditions may change nestling growth and development. Early life growth is sensitive to environmental effects and can potentially have long-lasting effects on adult phenotypes and fitness. A detailed understanding of both how and when weather conditions affect the entire growth trajectory of a nestling may help predict population changes in phenotypes and demography under climate change. This review covers three main topics on the impacts of weather variation (air temperature, rainfall, wind speed, solar radiation) on nestling growth. Firstly, we highlight why understanding the effects of weather on nestling growth might be important in understanding adaptation to, and population persistence in, environments altered by climate change. Secondly, we review the documented effects of weather variation on nestling growth curves. We investigate both altricial and precocial species, but we find a limited number of studies on precocial species in the wild. Increasing temperatures and rainfall have mixed effects on nestling growth, while increasing windspeeds tend to have negative impacts on the growth rate of open cup nesting species. Thirdly, we discuss how weather variation might affect the evolution of nestling growth traits and suggest that more estimates of the inheritance of and selection acting on growth traits in natural settings are needed to make evolutionary predictions. We suggest that predictions will be improved by considering concurrently changing selection pressures like urbanization. The importance of adaptive plastic or evolutionary changes in growth may depend on where a species or population is located geographically and the species’ life-history. Detailed characterization of the effects of weather on growth patterns will help answer whether variation in avian growth frequently plays a role in adaption to climate change.

Harnessing landscape genomics to identify future climate resilient genotypes in a desert annual

Local adaptation features critically in shaping species responses to changing environments, complicating efforts to revegetate degraded areas. Rapid climate change poses an additional challenge that could reduce fitness of even locally sourced seeds in restoration. Predictive restoration strategies that apply seeds with favourable adaptations to future climate may promote long-term resilience. Landscape genomics is increasingly used to assess spatial patterns in local adaption and may represent a cost-efficient approach for identifying future-adapted genotypes. To demonstrate such an approach, we genotyped 760 plants from 64 Mojave Desert populations of the desert annual Plantago ovata. Genome scans on 5,960 SNPs identified 184 potentially adaptive loci related to climate and satellite vegetation metrics. Causal modelling indicated that variation in potentially adaptive loci was not confounded by isolation by distance or isolation by habitat resistance. A generalized dissimilarity model (GDM) attributed spatial turnover in potentially adaptive loci to temperature, precipitation and NDVI amplitude, a measure of vegetation green-up potential. By integrating a species distribution model (SDM), we find evidence that summer maximum temperature may both constrain the range of P. ovata and drive adaptive divergence in populations exposed to higher temperatures. Within the species' current range, warm-adapted genotypes are predicted to experience a fivefold expansion in climate niche by midcentury and could harbour key adaptations to cope with future climate. We recommend eight seed transfer zones and project each zone into its relative position in future climate. Prioritizing seed collection efforts on genotypes with expanding future habitat represents a promising strategy for restoration practitioners to address rapidly changing climates.

Adaptation, ancestral variation and gene flow in a 'Sky Island' Drosophila species

Over time, populations of species can expand, contract, fragment and become isolated, creating subpopulations that must adapt to local conditions. Understanding how species maintain variation after divergence as well as adapt to these changes in the face of gene flow is of great interest, especially as the current climate crisis has caused range shifts and frequent migrations for many species. Here, we characterize how a mycophageous fly species, Drosophila innubila, came to inhabit and adapt to its current range which includes mountain forests in south-western USA separated by large expanses of desert. Using population genomic data from more than 300 wild-caught individuals, we examine four populations to determine their population history in these mountain forests, looking for signatures of local adaptation. In this first extensive study, establishing D. innubila as a key genomic "Sky Island" model, we find D. innubila spread northwards during the previous glaciation period (30-100 KYA) and have recently expanded even further (0.2-2 KYA). D. innubila shows little evidence of population structure, consistent with a recent establishment and genetic variation maintained since before geographic stratification. We also find some signatures of recent selective sweeps in chorion proteins and population differentiation in antifungal immune genes suggesting differences in the environments to which flies are adapting. However, we find little support for long-term recurrent selection in these genes. In contrast, we find evidence of long-term recurrent positive selection in immune pathways such as the Toll signalling system and the Toll-regulated antimicrobial peptides.

Microgeographical adaptation corresponds to elevational distributions of congeneric montane grasshoppers

Local adaptation can occur at small spatial scales relative to the dispersal capacity of species. Alpine ecosystems have sharp environmental clines that offer an opportunity to investigate the effects of fine-scale shifts in species' niche breadth on adaptive genetic processes. Here we examine two grasshopper species endemic to the Australian Alps (Kosciuscola spp.) that differ in elevational niche breadth: one broader, K. usitatus (1400-2200 m), and one narrower, K. tristis (1600-2000 m). We examine signatures of selection with respect to environmental and morphological variables in two mountain regions using F_ST outlier tests and environmental association analyses (EAAs) applied to single nucleotide polymorphism (SNP) data (K. usitatus: 9017 SNPs, n = 130; K. tristis: 7363 SNPs, n = 135). Stronger genetic structure was found in the more narrowly distributed K. tristis, which showed almost twice the number of SNPs under putative selection (10.8%) compared with K. usitatus (5.3%). When examining SNPs in common across species (n = 3058), 260 SNPs (8.5%) were outliers shared across species, and these were mostly associated with elevation, a proxy for temperature, suggesting parallel adaptive processes in response to climatic drivers. Additive polygenic scores (an estimate of the cumulative signal of selection across all candidate loci) were nonlinearly and positively correlated with elevation in both species. However, a steeper correlation in K. tristis indicated a stronger signal of spatially varying selection towards higher elevations. Our study illustrates that the niche breadth of co-occurring and related species distributed along the same environmental cline is associated with differences in patterns of microgeographical adaptation.

Rapid Adaptation to Temperature via a Potential Genomic Island of Divergence in the Invasive Green Crab, Carcinus maenas

Widespread species often adapt easily to novel conditions – both those found in new habitats and those generated by climate change. However, rapid adaptation may be hindered in the marine realm, where long-distance dispersal and consequently high gene flow are predicted to limit potential for local adaptation. Here, we use a highly dispersive invasive marine crab to test the nature and speed of adaptation to temperature in the sea. Using single nucleotide polymorphisms (SNPs) generated from cardiac transcriptome sequencing, we characterized six populations of the European green crab (Carcinus maenas) located across parallel thermal gradients in their native and invasive ranges. We compared SNP frequencies with local temperatures and previously generated data on cardiac heat and cold tolerance to identify candidate markers associated with population-level differences in thermal physiology. Of 10,790 SNPs, 104 were identified as frequency outliers, a signal that was strongly driven by association with temperature and/or cold tolerance. Seventy-two of these outlier markers, representing 28 different genes, were in a cluster of SNPs identified as a potential inversion polymorphism using linkage disequilibrium network analysis. This SNP cluster was unique in the data set, which was otherwise characterized by low levels of linkage disequilibrium, and markers in this cluster showed a significant enrichment of coding substitutions relative to the full SNP set. These 72 outlier SNPs appear to be transmitted as a unit, and represent a putative genomic island of divergence which varied in frequency with organismal cold tolerance. This relationship was strikingly similar across both native and invasive populations, all of which showed a very strong correlation with cold tolerance (R2 = 0.96 over all six populations). Notably, three of these populations have diverged recently (&lt;100 years) and show little to no neutral divergence, suggesting that this genomic region may be responding to temperature on a relatively short time scale. This relationship indicates adaptation to temperature based on the action of a putative genomic island of divergence, perhaps partially explaining the extraordinary invasive ability of this species.

Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal, poorly dispersing snail

Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab-adapted ecotypes that exhibit partial reproductive isolation from wave-adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co-association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis. In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes-consistent with models of adaptation with gene flow-but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co-association of loci during ecological speciation with ongoing gene flow.

Genomic Prediction of (Mal)Adaptation Across Current and Future Climatic Landscapes

Signals of local adaptation have been found in many plants and animals, highlighting the heterogeneity in the distribution of adaptive genetic variation throughout species ranges. In the coming decades, global climate change is expected to induce shifts in the selective pressures that shape this adaptive variation. These changes in selective pressures will likely result in varying degrees of local climate maladaptation and spatial reshuffling of the underlying distributions of adaptive alleles. There is a growing interest in using population genomic data to help predict future disruptions to locally adaptive gene-environment associations. One motivation behind such work is to better understand how the effects of changing climate on populations’ short-term fitness could vary spatially across species ranges. Here we review the current use of genomic data to predict the disruption of local adaptation across current and future climates. After assessing goals and motivationsunderlying the approach, we review the main steps and associated statistical methods currently in use and explore our current understanding of the limits and future potential of using genomics to predict climate change (mal)adaptation.

Gene flow relates to evolutionary divergence among populations at the range margin

Background

Morphological differentiation between populations resulting from local adaptations to environmental conditions is likely to be more pronounced in populations with increasing genetic isolation. In a previous study a positive clinal variation in body size was observed in isolated Roesel’s bush-cricket, Metrioptera roeselii, populations, but were absent from populations within a continuous distribution at the same latitudinal range. This observational study inferred that there was a phenotypic effect of gene flow on climate-induced selection in this species.

Methods

To disentangle genetic versus environmental drivers of population differences in morphology, we measured the size of four different body traits in wild-caught individuals from the two most distinct latitudinally-matched pairs of populations occurring at about 60°N latitude in northern Europe, characterised by either restricted or continuous gene flow, and corresponding individuals raised under laboratory conditions.

Results

Individuals that originated from the genetically isolated populations were always bigger (femur, pronotum and genital appendages) when compared to individuals from latitudinally-matched areas characterised by continuous gene flow between populations. The magnitude of this effect was similar for wild-caught and laboratory-reared individuals. We found that previously observed size cline variation in both male and female crickets was likely to be the result of local genetic adaptation rather than phenotypic plasticity.

Conclusions

This strongly suggests that restricted gene flow is of major importance for frequencies of alleles that participate in climate-induced selection acting to favour larger phenotypes in isolated populations towards colder latitudes.

A population genomics approach to uncover the CNVs, and their evolutionary significance, hidden in reduced-representation sequencing data sets

The importance of structural variation in adaptation and speciation is becoming increasingly evident in the literature. Among SVs, copy number variants (CNVs) are known to affect phenotypes through changes in gene expression and can potentially reduce recombination between alleles with different copy numbers. However, little is known about their abundance, distribution and frequency in natural populations. In a "From the Cover" article in this issue of Molecular Ecology, Dorant et al. (2020) present a new cost-effective approach to genotype copy number variants (CNVs) from large reduced-representation sequencing (RRS) data sets in nonmodel organisms, and thus to analyse sequence and structural variation jointly. They show that in American lobsters (Homarus americanus), CNVs exhibit strong population structure and several significant associations with annual variance in sea surface temperature, while SNPs fail to uncover any population structure or genotype-environment associations. Their results clearly illustrate that structural variants like CNVs can potentially store important information on differentiation and adaptive differences that cannot be retrieved from the analysis of sequence variation alone. To better understand the factors affecting the evolution of CNVs and their role in adaptation and speciation, we need to compare and synthesize data from a wide variety of species with different demographic histories and genome structure. The approach developed by Dorant et al. (2020) now allows to gain crucial knowledge on CNVs in a cost-effective way, even in species with limited genomic resources.

Genomic signatures of convergent adaptation to Alpine environments in three Brassicaceae species

It has long been discussed to what extent related species develop similar genetic mechanisms to adapt to similar environments. Most studies documenting such convergence have either used different lineages within species or surveyed only a limited portion of the genome. Here, we investigated whether similar or different sets of orthologous genes were involved in genetic adaptation of natural populations of three related plant species to similar environmental gradients in the Alps. We used whole-genome pooled population sequencing to study genome-wide SNP variation in 18 natural populations of three Brassicaceae (Arabis alpina, Arabidopsis halleri, and Cardamine resedifolia) from the Swiss Alps. We first de novo assembled draft reference genomes for all three species. We then ran population and landscape genomic analyses with ~3 million SNPs per species to look for shared genomic signatures of selection and adaptation in response to similar environmental gradients acting on these species. Genes with a signature of convergent adaptation were found at significantly higher numbers than expected by chance. The most closely related species pair showed the highest relative over-representation of shared adaptation signatures. Moreover, the identified genes of convergent adaptation were enriched for nonsynonymous mutations, suggesting functional relevance of these genes, even though many of the identified candidate genes have hitherto unknown or poorly described functions based on comparison with Arabidopsis thaliana. We conclude that adaptation to heterogeneous Alpine environments in related species is partly driven by convergent evolution, but that most of the genomic signatures of adaptation remain species-specific.

Genome-wide genetic diversity yields insights into genomic responses of candidate climate-selected loci in an Andean wetland plant

Assessing population evolutionary potential has become a central tenet of conservation biology. Since adaptive responses require allelic variation at functional genes, consensus has grown that genetic variation at genes under selection is a better surrogate for adaptive evolutionary potential than neutral genetic diversity. Although consistent with prevailing theory, this argument lacks empirical support and ignores recent theoretical advances questioning the very concept of neutral genetic diversity. In this study, we quantified genome-wide responses of single nucleotide polymorphism loci linked to climatic factors over a strong latitudinal gradient in natural populations of the high Andean wetland plant, Carex gayana, and then assessed whether genetic variation of candidate climate-selected loci better predicted their genome-wide responses than genetic variation of non-candidate loci. Contrary to this expectation, genomic responses of climate-linked loci only related significantly to environmental variables and genetic diversity of non-candidate loci. The effects of genome-wide genetic diversity detected in this study may be a result of either the combined influence of small effect variants or neutral and demographic factors altering the adaptive evolutionary potential of C. gayana populations. Regardless of the processes involved, our results redeem genome-wide genetic diversity as a potentially useful indicator of population adaptive evolutionary potential.

Small-scale population divergence is driven by local larval environment in a temperate amphibian

Genomic variation within and among populations is shaped by the interplay between natural selection and the effects of genetic drift and gene flow. Adaptive divergence can be found in small-scale natural systems even when population sizes are small, and the potential for gene flow is high, suggesting that local environments exert selection pressures strong enough to counteract the opposing effects of drift and gene flow. Here, we investigated genomic differentiation in nine moor frog (Rana arvalis) populations in a small-scale network of local wetlands using 16,707 ddRAD-seq SNPs, relating levels of differentiation with local environments, as well as with properties of the surrounding landscape. We characterized population structure and differentiation, and partitioned the effects of geographic distance, local larval environment, and landscape features on total genomic variation. We also conducted gene-environment association studies using univariate and multivariate approaches. We found small-scale population structure corresponding to 6-8 clusters. Local larval environment was the most influential component explaining 2.3% of the total genetic variation followed by landscape features (1.8%) and geographic distance (0.8%), indicative of isolation-by-environment, -by-landscape, and -by-distance, respectively. We identified 1000 potential candidate SNPs putatively under divergent selection mediated by the local larval environment. The candidate SNPs were involved in, among other biological functions, immune system function and development. Our results suggest that small-scale environmental differences can exert selection pressures strong enough to counteract homogenizing effects of gene flow and drift in this small-scale system, leading to observable population differentiation.

Strong Gene Flow Undermines Local Adaptations in a Host Parasite System

Simple Summary

The co-evolution of hosts and parasites depends on their ability to adapt to each other’s defense and counter-defense mechanisms. The strength of selection on those mechanisms may vary among populations, resulting in a geographical mosaic of co-evolution. The boreo-montane paper wasp Polistes biglumis and its parasite Polistes atrimandibularis exemplify this type of co-evolutionary system. Here, we used genetic markers to examine the genetic population structures of these wasps in the western Alps. We found that both host and parasite populations displayed similar levels of genetic variation. In the host species, populations located near to each other were genetically similar; in both the host and the parasite species populations farther apart were significantly different. Thus, apparent dispersal barriers (i.e., high mountains) did not seem to restrict gene flow across populations as expected. Furthermore, there were no major differences in gene flow between the two species, perhaps because P. atrimandibularis parasitizes both alpine and lowland host species and annually migrates between alpine and lowland populations. The presence of strong gene flow in a system where local populations experience variable levels of selection pressure challenges the classical hypothesis that restricted gene flow is required for local adaptations to evolve.

Abstract

The co-evolutionary pathways followed by hosts and parasites strongly depend on the adaptive potential of antagonists and its underlying genetic architecture. Geographically structured populations of interacting species often experience local differences in the strength of reciprocal selection pressures, which can result in a geographic mosaic of co-evolution. One example of such a system is the boreo-montane social wasp Polistes biglumis and its social parasite Polistes atrimandibularis, which have evolved local defense and counter-defense mechanisms to match their antagonist. In this work, we study spatial genetic structure of P. biglumis and P. atrimandibularis populations at local and regional scales in the Alps, by using nuclear markers (DNA microsatellites, AFLP) and mitochondrial sequences. Both the host and the parasite populations harbored similar amounts of genetic variation. Host populations were not genetically structured at the local scale, but geographic regions were significantly differentiated from each other in both the host and the parasite in all markers. The net dispersal inferred from genetic differentiation was similar in the host and the parasite, which may be due to the annual migration pattern of the parasites between alpine and lowland populations. Thus, the apparent dispersal barriers (i.e., high mountains) do not restrict gene flow as expected and there are no important gene flow differences between the species, which contradict the hypothesis that restricted gene flow is required for local adaptations to evolve.

Spatial grain of adaptation is much finer than ecoregional-scale common gardens reveal

Adaptive variation among plant populations must be known for effective conservation and restoration of imperiled species and predicting their responses to a changing climate. Common-garden experiments, in which plants sourced from geographically distant populations are grown together such that genetic differences may be expressed, have provided much insight on adaptive variation. Common-garden experiments also form the foundation for climate-based seed-transfer guidelines. However, the spatial scale at which population differentiation occurs is rarely addressed, leaving a critical information gap for parameterizing seed-transfer guidelines and assessing species' climate vulnerability. We asked whether adaptation was evident among populations of a foundational perennial within a single "empirical" seed-transfer zone (based on previous common-garden findings evaluating very distant populations) but different "provisional" seed zones (groupings of areas of similar climate and are not parameterized from common-garden data). Seedlings from three populations originating from similar conditions within an intermediate elevation were planted into gardens nearby at the same elevation, or 250-450 m higher or lower in elevation and 0.4-25 km away. Substantial variation was observed between gardens in survival (ranging 2%-99%), foliar crown volume (7.8-22.6 dm3), and reproductive effort (0%-65%), but not among the three transplanted populations. The between garden variation was inversely related to climatic differences between the gardens and seed-source populations, specifically the site differences in maximum-minimum annual temperatures. Results suggest that substantial site-specificity in adaptation can occur at finer scales than is accounted for in empirical seed-transfer guidance when the guidance is derived from broadscale common-garden studies. Being within the same empirical seed zone, geographic unit, and even within 10 km distance may not qualify as "local" in the context of seed transfer. Moving forward, designing common-garden experiments so that they allow for testing the scale of adaptation will help in translating the resulting seed-transfer guidance to restoration projects.

Spatial gradients of introgressed ancestry reveal cryptic connectivity patterns in a high gene flow marine fish

Assessing genetic connectivity among populations in high gene flow species is sometimes insufficient to evaluate demographic connectivity. Genetic differentiation quickly becomes zero as soon as a few dozen migrants are exchanged per generation. This provides little information to determine whether migration can ensure demographic coupling. The resulting difficulties in delineating conservation units for the management of commercially exploited marine fish species are well illustrated in the case of the European sea bass (Dicentrarchus labrax). Previous attempts to assess connectivity patterns in the northeast Atlantic have been hampered by a lack of spatial genetic structure. In contrast, mark-recapture data suggested low migration rates between regional spawning areas. Here, we show how a spatial gradient of introgressed Mediterranean ancestry across the northeast Atlantic reflects cryptic patterns of genetic and demographic connectivity. Using a 1K SNP chip data set in 827 individuals sampled from Portugal to the North Sea, we found null overall genetic differentiation across the northeast Atlantic. We however detected a subtle latitudinal admixture gradient originating at the edge of the contact zone with the Mediterranean sea bass lineage. Two significant breaks in the ancestry gradient at the tip of Galicia and northern Brittany indicated barriers to effective dispersal between demographically distinct units. Moreover, a northward expansion signal in Irish and North Seas was revealed by the surfing of rare Mediterranean alleles at the edge of the species range. Our results show that introgressed ancestry gradients offer a powerful alternative to assess genetic and demographic connectivity when the neutral migration-drift balance is not informative.

Footprints of local adaptation span hundreds of linked genes in the Atlantic silverside genome

The study of local adaptation in the presence of ongoing gene flow is the study of natural selection in action, revealing the functional genetic diversity most relevant to contemporary pressures. In addition to individual genes, genome‐wide architecture can itself evolve to enable adaptation. Distributed across a steep thermal gradient along the east coast of North America, Atlantic silversides (Menidia menidia) exhibit an extraordinary degree of local adaptation in a suite of traits, and the capacity for rapid adaptation from standing genetic variation, but we know little about the patterns of genomic variation across the species range that enable this remarkable adaptability. Here, we use low‐coverage, whole‐transcriptome sequencing of Atlantic silversides sampled along an environmental cline to show marked signatures of divergent selection across a gradient of neutral differentiation. Atlantic silversides sampled across 1371 km of the southern section of its distribution have very low genome‐wide differentiation (median F_ST = 0.006 across 1.9 million variants), consistent with historical connectivity and observations of recent migrants. Yet almost 14,000 single nucleotide polymorphisms (SNPs) are nearly fixed (F_ST > 0.95) for alternate alleles. Highly differentiated SNPs cluster into four tight linkage disequilibrium (LD) blocks that span hundreds of genes and several megabases. Variants in these LD blocks are disproportionately nonsynonymous and concentrated in genes enriched for multiple functions related to known adaptations in silversides, including variation in lipid storage, metabolic rate, and spawning behavior. Elevated levels of absolute divergence and demographic modeling suggest selection maintaining divergence across these blocks under gene flow. These findings represent an extreme case of heterogeneity in levels of differentiation across the genome, and highlight how gene flow shapes genomic architecture in continuous populations. Locally adapted alleles may be common features of populations distributed along environmental gradients, and will likely be key to conserving variation to enable future responses to environmental change.

Consequences of Single-Locus and Tightly Linked Genomic Architectures for Evolutionary Responses to Environmental Change

Genetic and genomic architectures of traits under selection are key factors influencing evolutionary responses. Yet, knowledge of their impacts has been limited by a widespread assumption that most traits are controlled by unlinked polygenic architectures. Recent advances in genome sequencing and eco-evolutionary modeling are unlocking the potential for integrating genomic information into predictions of population responses to environmental change. Using eco-evolutionary simulations, we demonstrate that hypothetical single-locus control of a life history trait produces highly variable and unpredictable harvesting-induced evolution relative to the classically applied multilocus model. Single-locus control of complex traits is thought to be uncommon, yet blocks of linked genes, such as those associated with some types of structural genomic variation, have emerged as taxonomically widespread phenomena. Inheritance of linked architectures resembles that of single loci, thus enabling single-locus-like modeling of polygenic adaptation. Yet, the number of loci, their effect sizes, and the degree of linkage among them all occur along a continuum. We review how linked architectures are often associated, directly or indirectly, with traits expected to be under selection from anthropogenic stressors and are likely to play a large role in adaptation to environmental disturbance. We suggest using single-locus models to explore evolutionary extremes and uncertainties when the trait architecture is unknown, refining parameters as genomic information becomes available, and explicitly incorporating linkage among loci when possible. By overestimating the complexity (e.g., number of independent loci) of the genomic architecture of traits under selection, we risk underestimating the complexity (e.g., nonlinearity) of their evolutionary dynamics.

Spatial predictors of genomic and phenotypic variation differ in a lowland Middle American bird (Icterus gularis)

Spatial patterns of intraspecific variation are shaped by geographical distance among populations, historical changes in gene flow and interactions with local environments. Although these factors are not mutually exclusive and operate on both genomic and phenotypic variation, it is unclear how they affect these two axes of variation. We address this question by exploring the predictors of genomic and phenotypic divergence in Icterus gularis, a broadly distributed Middle American bird that exhibits marked geographical variation in body size across its range. We combined a comprehensive single nucleotide polymorphism and phenotypic data set to test whether genome-wide genetic and phenotypic differentiation are best explained by (i) isolation by distance, (ii) isolation by history or (iii) isolation by environment. We find that the pronounced genetic and phenotypic variation in I. gularis are only partially correlated and differ regarding spatial predictors. Whereas genomic variation is largely explained by historical barriers to gene flow, phenotypic diversity can be best predicted by contemporary environmental heterogeneity. Our genomic analyses reveal strong phylogeographical structure coinciding with the Chivela Pass at the Isthmus of Tehuantepec that was formed during the Pleistocene, when populations were isolated in north-south refugia. In contrast, we found a strong association between body size and environmental variables, such as temperature and precipitation. The relationship between body size and local climate is consistent with a pattern produced by either natural selection or environmental plasticity. Overall, these results provide empirical evidence for why phenotypic and genomic data are often in conflict in taxonomic and phylogeographical studies.

Replicated Landscape Genomics Identifies Evidence of Local Adaptation to Urbanization in Wood Frogs

Native species that persist in urban environments may benefit from local adaptation to novel selection factors. We used double-digest restriction-side associated DNA (RAD) sequencing to evaluate shifts in genome-wide genetic diversity and investigate the presence of parallel evolution associated with urban-specific selection factors in wood frogs (Lithobates sylvaticus). Our replicated paired study design involved 12 individuals from each of 4 rural and urban populations to improve our confidence that detected signals of selection are indeed associated with urbanization. Genetic diversity measures were less for urban populations; however, the effect size was small, suggesting little biological consequence. Using an FST outlier approach, we identified 37 of 8344 genotyped single nucleotide polymorphisms with consistent evidence of directional selection across replicates. A genome-wide association study analysis detected modest support for an association between environment type and 12 of the 37 FST outlier loci. Discriminant analysis of principal components using the 37 FST outlier loci produced correct reassignment for 87.5% of rural samples and 93.8% of urban samples. Eighteen of the 37 FST outlier loci mapped to the American bullfrog (Rana [Lithobates] catesbeiana) genome, although none were in coding regions. This evidence of parallel evolution to urban environments provides a powerful example of the ability of urban landscapes to direct evolutionary processes.

Genetic architecture influences when and how hybridization contributes to colonization

The role of genetic architecture in adaptation to novel environments has received considerable attention when the source of adaptive variation is de novo mutation. Relatively less is known when the source of adaptive variation is inter- or intraspecific hybridization. We model hybridization between divergent source populations and subsequent colonization of an unoccupied novel environment using individual-based simulations to understand the influence of genetic architecture on the timing of colonization and the mode of adaptation. We find that two distinct categories of genetic architecture facilitate rapid colonization but that they do so in qualitatively different ways. For few and/or tightly linked loci, the mode of adaptation is via the recovery of adaptive parental genotypes. With many unlinked loci, the mode of adaptation is via the generation of novel hybrid genotypes. The first category results in the shortest colonization lag phases across the widest range of parameter space, but further adaptation is mutation limited. The second category takes longer and is more sensitive to genetic variance and dispersal rate, but can facilitate adaptation to environmental conditions that exceed the tolerance of parental populations. These findings have implications for understanding the origins of biological invasions and the success of hybrid populations.

Rapid local adaptation linked with phenotypic plasticity

Models of "plasticity-first" evolution are attractive because they explain the rapid evolution of new complex adaptations. Nevertheless, it is unclear whether plasticity can facilitate rapid microevolutionary change between diverging populations. Here, we show how plasticity may have generated adaptive differences in fecundity between neighboring wild populations of burying beetles Nicrophorus vespilloides. These populations occupy distinct Cambridgeshire woodlands that are just 2.5 km apart and that probably originated from a common ancestral population about 1000-4000 years ago. We find that populations are divergently adapted to breed on differently sized carrion. Adaptive differences in clutch size and egg size are associated with divergence at loci connected with oogenesis. The populations differ specifically in the elevation of the reaction norm linking clutch size to carrion size (i.e., genetic accommodation), and in the likelihood that surplus offspring will be lost after hatching. We suggest that these two processes may have facilitated rapid local adaptation on a fine-grained spatial scale.

A Roadmap for Understanding the Evolutionary Significance of Structural Genomic Variation

Structural genomic variants (SVs) are ubiquitous and play a major role in adaptation and speciation. Yet, comparative and population genomics have focused predominantly on gene duplications and large-effect inversions. The lack of a common framework for studying all SVs is hampering progress towards a more systematic assessment of their evolutionary significance. Here we (i) review how different types of SVs affect ecological and evolutionary processes; (ii) suggest unifying definitions and recommendations for future studies; and (iii) provide a roadmap for the integration of SVs in ecoevolutionary studies. In doing so, we lay the foundation for population genomics, theoretical, and experimental approaches to understand how the full spectrum of SVs impacts ecological and evolutionary processes.

Genome-scale data reveal the role of hybridization in lichen-forming fungi

Advancements in molecular genetics have revealed that hybridization may be common among plants, animals, and fungi, playing a role in evolutionary dynamics and speciation. While hybridization has been well-documented in pathogenic fungi, the effects of these processes on speciation in fungal lineages with different life histories and ecological niches are largely unexplored. Here we investigated the potential influence of hybridization on the emergence of morphologically and reproductively distinct asexual lichens. We focused on vagrant forms (growing obligately unattached to substrates) within a clade of rock-dwelling, sexually reproducing species in the Rhizoplaca melanophthalma (Lecanoraceae, Ascomycota) species complex. We used phylogenomic data from both mitochondrial and nuclear genomes to infer evolutionary relationships and potential patterns of introgression. We observed multiple instances of discordance between the mitochondrial and nuclear trees, including the clade comprising the asexual vagrant species R. arbuscula, R. haydenii, R. idahoensis, and a closely related rock-dwelling lineage. Despite well-supported phylogenies, we recovered strong evidence of a reticulated evolutionary history using a network approach that incorporates both incomplete lineage sorting and hybridization. These data suggest that the rock-dwelling western North American subalpine endemic R. shushanii is potentially the result of a hybrid speciation event, and introgression may have also played a role in other taxa, including vagrant species R. arbuscula, R. haydenii and R. idahoensis. We discuss the potential roles of hybridization in terms of generating asexuality and novel morphological traits in lichens. Furthermore, our results highlight the need for additional study of reticulate phylogenies when investigating species boundaries and evolutionary history, even in cases with well-supported topologies inferred from genome-scale data.

Evolutionary genomics can improve prediction of species' responses to climate change

Global climate change (GCC) increasingly threatens biodiversity through the loss of species, and the transformation of entire ecosystems. Many species are challenged by the pace of GCC because they might not be able to respond fast enough to changing biotic and abiotic conditions. Species can respond either by shifting their range, or by persisting in their local habitat. If populations persist, they can tolerate climatic changes through phenotypic plasticity, or genetically adapt to changing conditions depending on their genetic variability and census population size to allow for de novo mutations. Otherwise, populations will experience demographic collapses and species may go extinct. Current approaches to predicting species responses to GCC begin to combine ecological and evolutionary information for species distribution modelling. Including an evolutionary dimension will substantially improve species distribution projections which have not accounted for key processes such as dispersal, adaptive genetic change, demography, or species interactions. However, eco‐evolutionary models require new data and methods for the estimation of a species' adaptive potential, which have so far only been available for a small number of model species. To represent global biodiversity, we need to devise large‐scale data collection strategies to define the ecology and evolutionary potential of a broad range of species, especially of keystone species of ecosystems. We also need standardized and replicable modelling approaches that integrate these new data to account for eco‐evolutionary processes when predicting the impact of GCC on species' survival. Here, we discuss different genomic approaches that can be used to investigate and predict species responses to GCC. This can serve as guidance for researchers looking for the appropriate experimental setup for their particular system. We furthermore highlight future directions for moving forward in the field and allocating available resources more effectively, to implement mitigation measures before species go extinct and ecosystems lose important functions.

Intraspecific diversification of the crop wild relative Brassica cretica Lam. using demographic model selection

BACKGROUND

Crop wild relatives (CWRs) contain genetic diversity, representing an invaluable resource for crop improvement. Many of their traits have the potential to help crops to adapt to changing conditions that they experience due to climate change. An impressive global effort for the conservation of various CWR will facilitate their use in crop breeding for food security. The genus Brassica is listed in Annex I of the International Treaty on Plant Genetic Resources for Food and Agriculture. Brassica oleracea (or wild cabbage), a species native to southern and western Europe, has become established as an important human food crop plant because of its large reserves stored over the winter in its leaves. Brassica cretica Lam. (Bc) is a CWR in the brassica group and B. cretica subsp. nivea (Bcn) has been suggested as a separate subspecies. The species Bc has been proposed as a potential gene donor to brassica crops, including broccoli, cabbage, cauliflower, oilseed rape, etc. RESULTS: We sequenced genomes of four Bc individuals, including two Bcn and two Bc. Demographic analysis based on our whole-genome sequence data suggests that populations of Bc are not isolated. Classification of the Bc into distinct subspecies is not supported by the data. Using only the non-coding part of the data (thus, the parts of the genome that has evolved nearly neutrally), we find the gene flow between different Bc population is recent and its genomic diversity is high.

CONCLUSIONS

Despite predictions on the disruptive effect of gene flow in adaptation, when selection is not strong enough to prevent the loss of locally adapted alleles, studies show that gene flow can promote adaptation, that local adaptations can be maintained despite high gene flow, and that genetic architecture plays a fundamental role in the origin and maintenance of local adaptation with gene flow. Thus, in the genomic era it is important to link the selected demographic models with the underlying processes of genomic variation because, if this variation is largely selectively neutral, we cannot assume that a diverse population of crop wild relatives will necessarily exhibit the wide-ranging adaptive diversity required for further crop improvement.