Widespread parallel population adaptation to climate variation across a radiation: implications for adaptation to climate change

Global warming will impact species in a number of ways, and it is important to know the extent to which natural populations can adapt to anthropogenic climate change by natural selection. Parallel microevolution within separate species can demonstrate natural selection, but several studies of homoplasy have not yet revealed examples of widespread parallel evolution in a generic radiation. Taking into account primary phylogeographic divisions, we investigate numerous quantitative traits (size, shape, scalation, colour pattern and hue) in anole radiations from the mountainous Lesser Antillean islands. Adaptation to climatic differences can lead to very pronounced differences between spatially close populations with all studied traits showing some evidence of parallel evolution. Traits from shape, scalation, pattern and hue (particularly the latter) show widespread evolutionary parallels within these species in response to altitudinal climate variation greater than extreme anthropogenic climate change predicted for 2080. This gives strong evidence of the ability to adapt to climate variation by natural selection throughout this radiation. As anoles can evolve very rapidly, it suggests anthropogenic climate change is likely to be less of a conservation threat than other factors, such as habitat loss and invasive species, in this, Lesser Antillean, biodiversity hot spot.     

Evidence for parallel adaptation to climate across the natural range of Arabidopsis thaliana

How organisms adapt to different climate habitats is a key question in evolutionary ecology and biological conservation. Species distributions are often determined by climate suitability. Consequently, the anthropogenic impact on earth's climate is of key concern to conservation efforts because of our relatively poor understanding of the ability of populations to track and evolve to climate change. Here, we investigate the ability of Arabidopsis thaliana to occupy climate space by quantifying the extent to which different climate regimes are accessible to different A. thaliana genotypes using publicly available data from a large‐scale genotyping project and from a worldwide climate database. The genetic distance calculated from 149 single‐nucleotide polymorphisms (SNPs) among 60 lineages of A. thaliana was compared to the corresponding climate distance among collection localities calculated from nine different climatic factors. A. thaliana was found to be highly labile when adapting to novel climate space, suggesting that populations may experience few constraints when adapting to changing climates. Our results also provide evidence of a parallel or convergent evolution on the molecular level supporting recent generalizations regarding the genetics of adaptation.     

Lineage‐specific adaptation to climate involves flowering time in North American Arabidopsis lyrata

Adaptation to local climatic conditions is commonly found within species, but whether it involves the same intraspecific genomic variants is unknown. We studied this question in North American Arabidopsis lyrata, whose current distribution is shaped by post‐glacial range expansion from two refugia, resulting in two distinct genetic clusters covering comparable climatic gradients. Using pooled whole‐genome sequence data of 41 outcrossing populations, we identified loci associated with three niche‐determining climatic variables in the two clusters and compared these outliers. Little evidence was found for parallelism in climate adaptation for single nucleotide polymorphisms (SNPs) and for genes with an accumulation of outlier SNPs. Significantly increased selection coefficients supported them as candidates of climate adaptation. However, the fraction of gene ontology (GO) terms shared between clusters was higher compared to outlier SNPs and outlier genes, suggesting that selection acts on similar pathways but not necessarily the same genes. Enriched GO terms involved responses to abiotic and biotic stress, circadian rhythm and development, with flower development and reproduction being among the most frequently detected. In line with GO enrichment, regulators of flowering time were detected as outlier genes. Our results suggest that while adaptation to environmental gradients on the genomic level are lineage‐specific in A. lyrata, similar biological processes seem to be involved. Differential loss of standing genetic variation, probably driven by genetic drift, can in part account for the lack of parallel evolution on the genomic level.     

Evidence of genomic adaptation to climate in Eucalyptus microcarpa: Implications for adaptive potential to projected climate change

Understanding whether populations can adapt in situ or whether interventions are required is of key importance for biodiversity management under climate change. Landscape genomics is becoming an increasingly important and powerful tool for rapid assessments of climate adaptation, especially in long‐lived species such as trees. We investigated climate adaptation in Eucalyptus microcarpa using the DArTseq genomic approach. A combination of F_ST outlier and environmental association analyses were performed using >4200 genomewide single nucleotide polymorphisms (SNPs) from 26 populations spanning climate gradients in southeastern Australia. Eighty‐one SNPs were identified as putatively adaptive, based on significance in F_ST outlier tests and significant associations with one or more climate variables related to temperature (70/81), aridity (37/81) or precipitation (35/81). Adaptive SNPs were located on all 11 chromosomes, with no particular region associated with individual climate variables. Climate adaptation appeared to be characterized by subtle shifts in allele frequencies, with no consistent fixed differences identified. Based on these associations, we predict adaptation under projected changes in climate will include a suite of shifts in allele frequencies. Whether this can occur sufficiently rapidly through natural selection within populations, or would benefit from assisted gene migration, requires further evaluation. In some populations, the absence or predicted increases to near fixation of particular adaptive alleles hint at potential limits to adaptive capacity. Together, these results reinforce the importance of standing genetic variation at the geographic level for maintaining species’ evolutionary potential.     

Genomic studies on the nature of species: adaptation and speciation in Mimulus

Evolutionary biology is in an exciting era, in which powerful genomic tools make the answers accessible to long‐standing questions about variation, adaptation and speciation. The availability of a suite of genomic resources, a shared knowledge base and a long history of study have made the phenotypically diverse plant genus Mimulus an important system for understanding ecological and evolutionary processes. An international Mimulus Research Meeting was held at Duke University in June 2014 to discuss developments in ecological and evolutionary genetic studies in Mimulus. Here, we report major recent discoveries presented at the meeting that use genomic approaches to advance our understanding of three major themes: the parallel genetic basis of adaptation; the ecological genomics of speciation; and the evolutionary significance of structural genetic variation. We also suggest future research directions for studies of Mimulus and highlight challenges faced when developing new ecological and evolutionary model systems.     

Rapid adaptation to climate change

In recent years, amid growing concerns that changing climate is affecting species distributions and ecosystems, predicting responses to rapid environmental change has become a major goal. In this issue, Franks and colleagues take a first step towards this objective (Franks et al. 2016). They examine genomewide signatures of selection in populations of Brassica rapa after a severe multiyear drought. Together with other authors, Franks had previously shown that flowering time was reduced after this particular drought and that the reduction was genetically encoded. Now, the authors have sequenced previously stored samples to compare allele frequencies before and after the drought and identify the loci with the most extreme shifts in frequencies. The loci they identify largely differ between populations, suggesting that different genetic variants may be responsible for reduction in flowering time in the two populations.     

Parallel evolution of an alpine type ecomorph in a scorpionfly: Independent adaptation to high‐altitude environments in multiple mountain locations

Elucidation of the diversification process of organisms is one of the important tasks of biology. From the viewpoint of species diversity, insects are the most successful group among the diverse organisms on earth and evolutionary adaptation is one of the important factors driving this pattern. Evolutionary adaptation is one of the important factors in the diversification of insects. One of the representative examples of environmental adaptation in insects is the shortening and loss of wings in subalpine and alpine zones. In this study, we focused on the Japanese scorpionfly, Panorpodes paradoxus. In this species, individuals that inhabit mountainous regions and subalpine zones have long wings (the “general type”), and individuals that inhabit higher altitudinal ranges have short wings (the “alpine type”). We collected samples of all Japanese Panorpodes species and one Korean Panorpodes species, and conducted molecular phylogenetic analyses of the mtDNA COI (610 bp), COII (688 bp), and 16S rRNA (888 bp) and nuDNA EF1‐α (658 bp) and 28S rRNA (524 bp) regions in order to reveal the evolutionary history of the alpine type of P. paradoxus. As a result of molecular phylogenetic analyses, it was revealed that the alpine type of P. paradoxus was polyphyletic, and had evolved to become the alpine type at least twice independently at separate mountain locations. In addition, the result of divergence time estimation suggested that the alpine type is an “ecomorph”, having recently adapted to low temperature habitats following mountain uplift within the Japanese Archipelago and subsequent glacial‐interglacial cycles.     

Breathing life into climate change adaptation

The exploration of evolutionary biology and biological adaptation can inform society's adaptation to climate change, particularly the mechanisms that bring about adaptability, such as phenotypic plasticity, epigenetics, and horizontal gene transfer. Learning from unplanned autonomous biological adaptation may be considered undesirable and incompatible with human endeavor. However, it is argued that there is no need for agency, and planned adaptation is not necessarily preferable over autonomous adaptation. What matters is the efficacy of adaptive mechanisms and their capacity to increase societal resilience to current and future impacts. In addition, there is great scope for industrial ecology (IE) to contribute approaches to climate change adaptation that generate system models and baseline data to inform decision making. The problem of “uncertainty” was chosen as an example of a challenge that is shared by biological systems, IE, and climate change adaptation to show how biological adaptation might contribute solutions. Finally, the Coastal Climate Adaptation Decision Support tool was used to demonstrate how IE and biological adaptation approaches may be mainstreamed in climate change adaptation planning and practice. In conclusion, there is close conceptual alignment between evolutionary biology and IE. The integration of biological adaptation thinking can enrich IE, add new perspectives to climate change adaptation science, and support IE's engagement with climate change adaptation. There should be no major obstacles regarding the collaboration of industrial ecologists with the climate change adaptation community, but mainstreaming of biological adaptation solutions depends greatly on successful knowledge transfer and the engagement of open‐minded and informed adaptation stakeholders.     

From nature to the laboratory: the impact of founder effects on adaptation

Most founding events entail a reduction in population size, which in turn leads to genetic drift effects that can deplete alleles. Besides reducing neutral genetic variability, founder effects can in principle shift additive genetic variance for phenotypes that underlie fitness. This could then lead to different rates of adaptation among populations that have undergone a population size bottleneck as well as an environmental change, even when these populations have a common evolutionary history. Thus, theory suggests that there should be an association between observable genetic variability for both neutral markers and phenotypes related to fitness. Here, we test this scenario by monitoring the early evolutionary dynamics of six laboratory foundations derived from founders taken from the same source natural population of Drosophila subobscura. Each foundation was in turn three‐fold replicated. During their first few generations, these six foundations showed an abrupt increase in their genetic differentiation, within and between foundations. The eighteen populations that were monitored also differed in their patterns of phenotypic adaptation according to their immediately ancestral founding sample. Differences in early genetic variability and in effective population size were found to predict differences in the rate of adaptation during the first 21 generations of laboratory evolution. We show that evolution in a novel environment is strongly contingent not only on the initial composition of a newly founded population but also on the stochastic changes that occur during the first generations of colonization. Such effects make laboratory populations poor guides to the evolutionary genetic properties of their ancestral wild populations.     

Improved genome assembly of Chinese shrimp (Fenneropenaeus chinensis) suggests adaptation to the environment during evolution and domestication

A high-quality reference genome is necessary to determine the molecular mechanisms underlying important biological phenomena; therefore, in the present study, a chromosome-level genome assembly of the Chinese shrimp Fenneropenaeus chinensis was performed. Muscle of a male shrimp was sequenced using PacBio platform, and assembled by Hi-C technology. The assembled F. chinensis genome was 1.47 Gb with contig N50 of 472.84 Kb, including 57.73% repetitive sequences, and was anchored to 43 pseudochromosomes, with scaffold N50 of 36.87 Mb. In total, 25,026 protein-coding genes were predicted. The genome size of F. chinensis showed significant contraction in comparison with that of other penaeid species, which is likely related to migration observed in this species. However, the F. chinensis genome included several expanded gene families related to cellular processes and metabolic processes, and the contracted gene families were associated with virus infection process. The findings signify the adaptation of F. chinensis to the selection pressure of migration and cold environment. Furthermore, the selection signature analysis identified genes associated with metabolism, phototransduction, and nervous system in cultured shrimps when compared with wild population, indicating targeted, artificial selection of growth, vision, and behavior during domestication. The construction of the genome of F. chinensis provided valuable information for the further genetic mechanism analysis of important biological processes, and will facilitate the research of genetic changes during evolution.     

Parallel genotypic adaptation: when evolution repeats itself

Until recently, parallel genotypic adaptation was considered unlikely because phenotypic differences were thought to be controlled by many genes. There is increasing evidence, however, that phenotypic variation sometimes has a simple genetic basis and that parallel adaptation at the genotypic level may be more frequent than previously believed. Here, we review evidence for parallel genotypic adaptation derived from a survey of the experimental evolution, phylogenetic, and quantitative genetic literature. The most convincing evidence of parallel genotypic adaptation comes from artificial selection experiments involving microbial populations. In some experiments, up to half of the nucleotide substitutions found in independent lineages under uniform selection are the same. Phylogenetic studies provide a means for studying parallel genotypic adaptation in non-experimental systems, but conclusive evidence may be difficult to obtain because homoplasy can arise for other reasons. Nonetheless, phylogenetic approaches have provided evidence of parallel genotypic adaptation across all taxonomic levels, not just microbes. Quantitative genetic approaches also suggest parallel genotypic evolution across both closely and distantly related taxa, but it is important to note that this approach cannot distinguish between parallel changes at homologous loci versus convergent changes at closely linked non-homologous loci. The finding that parallel genotypic adaptation appears to be frequent and occurs at all taxonomic levels has important implications for phylogenetic and evolutionary studies. With respect to phylogenetic analyses, parallel genotypic changes, if common, may result in faulty estimates of phylogenetic relationships. From an evolutionary perspective, the occurrence of parallel genotypic adaptation provides increasing support for determinism in evolution and may provide a partial explanation for how species with low levels of gene flow are held together.     

A genome‐wide search for local adaptation in a terrestrial‐breeding frog reveals vulnerability to climate change

Terrestrial‐breeding amphibians are likely to be vulnerable to warming and drying climates, as their embryos require consistent moisture for successful development. Adaptation to environmental change will depend on sufficient genetic variation existing within or between connected populations. Here, we use Single Nucleotide Polymorphism (SNP) data to investigate genome‐wide patterns in genetic diversity, gene flow and local adaptation in a terrestrial‐breeding frog (Pseudophryne guentheri) subject to a rapidly drying climate and recent habitat fragmentation. The species was sampled across 12 central and range‐edge populations (192 samples), and strong genetic structure was apparent, as were high inbreeding coefficients. Populations showed differences in genetic diversity, and one population lost significant genetic diversity in a decade. More than 500 SNP loci were putatively under directional selection, and 413 of these loci were correlated with environmental variables such as temperature, rainfall, evaporation and soil moisture. One locus showed homology to a gene involved in the activation of maturation in Xenopus oocytes, which may facilitate rapid development of embryos in drier climates. The low genetic diversity, strong population structuring and presence of local adaptation revealed in this study shows why management strategies such as targeted gene flow may be necessary to assist isolated populations to adapt to future climates.     

Genomic and phenotypic signatures of climate adaptation in an Anolis lizard

Integrated knowledge on phenotype, physiology, and genomic adaptations is required to understand the effects of climate on evolution. The functional genomic basis of organismal adaptation to changes in the abiotic environment, its phenotypic consequences, and its possible convergence across vertebrates are still understudied. In this study, we use a comparative approach to verify predicted gene functions for vertebrate thermal adaptation with observed functions underlying repeated genomic adaptations in response to elevation in the lizard Anolis cybotes. We establish a direct link between recurrently evolved phenotypes and functional genomics of altitude‐related climate adaptation in three highland and lowland populations in the Dominican Republic. We show that across vertebrates, genes contained in this interactome are expressed within the brain, the endocrine system, and during development. These results are relevant to elucidate the effect of global climate change across vertebrates and might aid in furthering insight into gene–environment relationships under disturbances to homeostasis.     

Experimental evidence for accelerated adaptation to desiccation through sexual selection on males

The impact of sexual selection on the adaptive process remains unclear. On the one hand, sexual selection might hinder adaptation by favouring costly traits and preferences that reduce nonsexual fitness. On the other hand, condition dependence of success in sexual selection may accelerate adaptation. Here, we used replicate populations of Drosophila melanogaster to artificially select on male desiccation resistance while manipulating the opportunity for precopulatory sexual selection in a factorial design. Following five generations of artificial selection, we measured the desiccation resistance of males and females to test whether the addition of sexual selection accelerated adaptation. We found a significant interaction between the effects of natural selection and sexual selection: desiccation resistance was highest in populations where sexual selection was allowed to operate. Despite only selecting on males, we also found a correlated response in females. These results provide empirical support for the idea that sexual selection can accelerate the rate of adaptation.     

Limited potential for adaptation to climate change in a broadly distributed marine crustacean

The extent to which acclimation and genetic adaptation might buffer natural populations against climate change is largely unknown. Most models predicting biological responses to environmental change assume that species' climatic envelopes are homogeneous both in space and time. Although recent discussions have questioned this assumption, few empirical studies have characterized intraspecific patterns of genetic variation in traits directly related to environmental tolerance limits. We test the extent of such variation in the broadly distributed tidepool copepod Tigriopus californicus using laboratory rearing and selection experiments to quantify thermal tolerance and scope for adaptation in eight populations spanning more than 17° of latitude. Tigriopus californicus exhibit striking local adaptation to temperature, with less than 1 per cent of the total quantitative variance for thermal tolerance partitioned within populations. Moreover, heat-tolerant phenotypes observed in low-latitude populations cannot be achieved in high-latitude populations, either through acclimation or 10 generations of strong selection. Finally, in four populations there was no increase in thermal tolerance between generations 5 and 10 of selection, suggesting that standing variation had already been depleted. Thus, plasticity and adaptation appear to have limited capacity to buffer these isolated populations against further increases in temperature. Our results suggest that models assuming a uniform climatic envelope may greatly underestimate extinction risk in species with strong local adaptation.     

All roads lead to weediness: Patterns of genomic divergence reveal extensive recurrent weedy rice origins from South Asian Oryza

Weedy rice (Oryza spp.), a weedy relative of cultivated rice (O. sativa), infests and persists in cultivated rice fields worldwide. Many weedy rice populations have evolved similar adaptive traits, considered part of the ‘agricultural weed syndrome’, making this an ideal model to study the genetic basis of parallel evolution. Understanding parallel evolution hinges on accurate knowledge of the genetic background and origins of existing weedy rice groups. Using population structure analyses of South Asian and US weedy rice, we show that weeds in South Asia have highly heterogeneous genetic backgrounds, with ancestry contributions both from cultivated varieties (aus and indica) and wild rice. Moreover, the two main groups of weedy rice in the USA, which are also related to aus and indica cultivars, constitute a separate origin from that of Asian weeds. Weedy rice populations in South Asia largely converge on presence of red pericarps and awns and on ease of shattering. Genomewide divergence scans between weed groups from the USA and South Asia, and their crop relatives are enriched for loci involved in metabolic processes. Some candidate genes related to iconic weedy traits and competitiveness are highly divergent between some weed‐crop pairs, but are not shared among all weed‐crop comparisons. Our results show that weedy rice is an extreme example of recurrent evolution, and suggest that most populations are evolving their weedy traits through different genetic mechanisms.     

Genomic evidence for parallel adaptation to cities

Urban evolutionary biology is the study of rapid evolutionary change in response to humans and our uses of land to support city dwellers. Because cities are relatively modern additions to the natural world, research on urban evolution tends to focus on microevolutionary change that has happened across a few to many hundreds of generations. These questions still fall under the broad purview of evolutionary ecology. However, the severity, rapidity and replication of environmental changes that drive evolution in this context make it worthy of specific attention. Urban evolution provides the opportunity to study the earliest stages of evolution in a context that is scientifically interesting and societally important. The newness of urban populations and their proximity to natural populations also creates challenges when trying to detect population genetic change. In a From the Cover article in this issue of Molecular Ecology, Mueller et al. use whole genome resequencing data to address some of these challenges while exploring genetic changes associated with urbanization in three replicate urban‐rural burrowing owl (Athene cunicularia) populations. Combining multiple approaches across these sample sites Mueller et al. find evidence for selection on genes whose function is related to synapses, neuron projections, brain connectivity and cognitive function in general. That selection was parallel suggests that phenotypes related to brain processes were probably particularly important for urban adaptation.     

Crassulacean acid metabolism (CAM) offers sustainable bioenergy production and resilience to climate change

Biomass production on low‐grade land is needed to meet future energy demands and minimize resource conflicts. This, however, requires improvements in plant water‐use efficiency (WUE) that are beyond conventional C3 and C4 dedicated bioenergy crops. Here we present the first global‐scale geographic information system (GIS)‐based productivity model of two highly water‐efficient crassulacean acid metabolism (CAM) candidates: Agave tequilana and Opuntia ficus‐indica. Features of these plants that translate to WUE advantages over C3 and C4 bioenergy crops include nocturnal stomatal opening, rapid rectifier‐like root hydraulic conductivity responses to fluctuating soil water potential and the capacity to buffer against periods of drought. Yield simulations for the year 2070 were performed under the four representative concentration pathway (RCPs) scenarios presented in the IPCC's 5th Assessment Report. Simulations on low‐grade land suggest that O. ficus‐indica alone has the capacity to meet ‘extreme’ bioenergy demand scenarios (>600 EJ yr^?1) and is highly resilient to climate change (?1%). Agave tequilana is moderately impacted (?11%). These results are significant because bioenergy demand scenarios >600 EJ yr^?1 could be met without significantly increasing conflicts with food production and contributing to deforestation. Both CAM candidates outperformed the C4 bioenergy crop, Panicum virgatum L. (switchgrass) in arid zones in the latitudinal range 30°S–30°N.     

The role of Bh4 in parallel evolution of hull colour in domesticated and weedy rice

The two independent domestication events in the genus Oryza that led to African and Asian rice offer an extremely useful system for studying the genetic basis of parallel evolution. This system is also characterized by parallel de‐domestication events, with two genetically distinct weedy rice biotypes in the US derived from the Asian domesticate. One important trait that has been altered by rice domestication and de‐domestication is hull colour. The wild progenitors of the two cultivated rice species have predominantly black‐coloured hulls, as does one of the two U.S. weed biotypes; both cultivated species and one of the US weedy biotypes are characterized by straw‐coloured hulls. Using Black hull 4 (Bh4) as a hull colour candidate gene, we examined DNA sequence variation at this locus to study the parallel evolution of hull colour variation in the domesticated and weedy rice system. We find that independent Bh4‐coding mutations have arisen in African and Asian rice that are correlated with the straw hull phenotype, suggesting that the same gene is responsible for parallel trait evolution. For the U.S. weeds, Bh4 haplotype sequences support current hypotheses on the phylogenetic relationship between the two biotypes and domesticated Asian rice; straw hull weeds are most similar to indica crops, and black hull weeds are most similar to aus crops. Tests for selection indicate that Asian crops and straw hull weeds deviate from neutrality at this gene, suggesting possible selection on Bh4 during both rice domestication and de‐domestication.     

Landscape genomics of Colorado potato beetle provides evidence of polygenic adaptation to insecticides

The ability of insect pests to rapidly and repeatedly adapt to insecticides has long challenged entomologists and evolutionary biologists. Since Crow's seminal paper on insecticide resistance in 1957, new data and insights continue to emerge that are relevant to the old questions about how insecticide resistance evolves: such as whether it is predominantly mono‐ or polygenic, and evolving from standing vs. de novo genetic variation. Many studies support the monogenic hypothesis, and current management recommendations assume single‐ or two‐locus models. But inferences could be improved by integrating data from a broader sample of pest populations and genomes. Here, we generate evidence relevant to these questions by applying a landscape genomics framework to the study of insecticide resistance in a major agricultural pest, Colorado potato beetle, Leptinotarsa decemlineata (Say). Genome–environment association tests using genomic variation from 16 populations spanning gradients of landscape variables associated with insecticide exposure over time revealed 42 strong candidate insecticide resistance genes, with potentially overlapping roles in multiple resistance mechanisms. Measurements of resistance to a widely used insecticide, imidacloprid, among 47 L. decemlineata populations revealed heterogeneity at a small (2 km) scale and no spatial signature of origin or spread throughout the landscape. Analysis of nucleotide diversity suggested candidate resistance loci have undergone varying degrees of selective sweeps, often maintaining similar levels of nucleotide diversity to neutral loci. This study suggests that many genes are involved in insecticide resistance in L. decemlineata and that resistance likely evolves from both de novo and standing genetic variation.