Environmental Association Analyses Research Articles

Signatures of natural selection in a foundation tree along Mediterranean climatic gradients.

Temperature and precipitation regimes are rapidly changing, resulting in forest dieback and extinction events, particularly in Mediterranean‐type climates (MTC). Forest management that enhance forests’ resilience is urgently required, however adaptation to climates in heterogeneous landscapes with multiple selection pressures is complex. For widespread trees in MTC we hypothesized that: patterns of local adaptation are associated with climate; precipitation is a stronger factor of adaptation than temperature; functionally related genes show similar signatures of adaptation; and adaptive variants are independently sorting across the landscape. We sampled 28 populations across the geographic distribution of Eucalyptus marginata (jarrah), in South‐west Western Australia, and obtained 13,534 independent single nucleotide polymorphic (SNP) markers across the genome. Three genotype‐association analyses that employ different ways of correcting population structure were used to identify putatively adapted SNPs associated with independent climate variables. While overall levels of population differentiation were low (F ST = 0.04), environmental association analyses found a total of 2336 unique SNPs associated with temperature and precipitation variables, with 1440 SNPs annotated to genic regions. Considerable allelic turnover was identified for SNPs associated with temperature seasonality and mean precipitation of the warmest quarter, suggesting that both temperature and precipitation are important factors in adaptation. SNPs with similar gene functions had analogous allelic turnover along climate gradients, while SNPs among temperature and precipitation variables had uncorrelated patterns of adaptation. These contrasting patterns provide evidence that there may be standing genomic variation adapted to current climate gradients, providing the basis for adaptive management strategies to bolster forest resilience in the future.

Genetic divergence and local adaptation of Liriodendron driven by heterogeneous environments.

Ecological adaptive differentiation alters both the species diversity and intraspecific genetic diversity in forests, thus affecting the stability of forest ecosystems. Therefore, knowledge of the genetic underpinnings of the ecological adaptive differentiation of forest species is critical for effective species conservation. In this study, single-nucleotide polymorphisms (SNPs) from population transcriptomes were used to investigate the spatial distribution of genetic variation in Liriodendron to assess whether environmental variables can explain genetic divergence. We examined the contributions of environmental variables to population divergence and explored the genetic underpinnings of local adaptation using a landscape genomic approach. Niche models and statistical analyses showed significant niche divergence between L. chinense and L. tulipifera, suggesting that ecological adaptation may play a crucial role in driving interspecific divergence. We detected a new fine-scale genetic structure in L. chinense, and divergence of the six groups occurred during the late Pliocene to early Pleistocene. Redundancy analysis (RDA) revealed significant associations between genetic variation and multiple environmental variables. Environmental association analyses identified 67 environmental association loci (EALs; nonsynonymous SNPs) that underwent interspecific or intraspecific differentiation, 28 of which were associated with adaptive genes. These 28 candidate adaptive loci provide substantial evidence for local adaptation in Liriodendron. Our findings reveal ecological adaptive divergence pattern between Liriodendron species and provide novel insight into the role of heterogeneous environments in shaping genetic structure and driving local adaptation among populations, informing future L. chinense conservation efforts.

Adaptive divergence in bill morphology and other thermoregulatory traits is facilitated by restricted gene flow in song sparrows on the California Channel Islands.

Disentangling the effects of neutral and adaptive processes in maintaining phenotypic variation across environmental gradients is challenging in natural populations. Song sparrows (Melospiza melodia) on the California Channel Islands occupy a pronounced east-west climate gradient within a small spatial scale, providing a unique opportunity to examine the interaction of genetic isolation (reduced gene flow) and the environment (selection) in driving variation. We used reduced representation genomic libraries to infer the role of neutral processes (drift and restricted gene flow) and divergent selection in driving variation in thermoregulatory traits with an emphasis on the mechanisms that maintain bill divergence among islands. Analyses of 22,029 neutral SNPs confirm distinct population structure by island with restricted gene flow and relatively large effective population sizes, suggesting bill differences are probably not a product of genetic drift. Instead, we found strong support for local adaptation using 3294 SNPs in differentiation-based and environmental association analyses coupled with genome-wide association tests. Specifically, we identified several putatively adaptive and candidate loci in or near genes involved in bill development pathways (e.g., BMP, CaM, Wnt), confirming the highly complex and polygenic architecture underlying bill morphology. Furthermore, we found divergence in genes associated with other thermoregulatory traits (i.e., feather structure, plumage colour, and physiology). Collectively, these results suggest strong divergent selection across an island archipelago results in genomic changes in a suite of traits associated with climate adaptation over small spatial scales. Future research should move beyond studying univariate traits to better understand multidimensional responses to complex environmental conditions.

Microbial-induced concrete corrosion under high-salt conditions: Microbial community composition and environmental multivariate association analysis

The corrosion of concrete materials in urban sewage pipes by microorganisms is a serious issue in wastewater networks around the world. There has been no systematic study conducted on material corrosion by microorganisms, particularly for high-salt environments in coastal cities. This study compared and analysed reactors under different salinities, and the 16S ribosomal ribonucleic acid gene sequence method was used to analyze the bacterial communities associated with concrete corrosion. Redundancy analysis demonstrated that chemical oxygen demand, volatile fatty acids, and sulphate altered the structure and distribution of the microbial community. The predominant bacteria, Proteobacteria, accounted for 41.85% of the seawater group. Among them, the sulphur functional microorganism, Desulfomicrobium, accounted for 4.14%. These bacteria can decompose macromolecular organic matter to provide energy for reproduction. Furthermore, they continue to provide sulphur for the eosinophilic sulphur-oxidising bacteria attached to the surface of the high alkaline concrete sample. The aggregated sulphur-oxidising bacteria produce biological sulphuric acid, leading to corrosion and damage to the concrete structure. Salinity promoted the aggregation of corrosion-inducing bacteria, accelerating the growth of corroding microorganisms on the concrete material of coastal urban sewage pipes.

Finding Needles in a Haystack: Using Geo-References to Enhance the Selection and Utilization of Landraces in Breeding for Climate-Resilient Cultivars of Upland Cotton (Gossypium hirsutum L.).

The genetic uniformity of cultivated cotton as a consequence of domestication and modern breeding makes it extremely vulnerable to abiotic challenges brought about by major climate shifts. To sustain productivity amidst worsening agro-environments, future breeding objectives need to seriously consider introducing new genetic variation from diverse resources into the current germplasm base of cotton. Landraces are genetically heterogeneous, population complexes that have been primarily selected for their adaptability to specific localized or regional environments. This makes them an invaluable genetic resource of novel allelic diversity that can be exploited to enhance the resilience of crops to marginal environments. The utilization of cotton landraces in breeding programs are constrained by the phenology of the plant and the lack of phenotypic information that can facilitate efficient selection of potential donor parents for breeding. In this review, the genetic value of cotton landraces and the major challenges in their utilization in breeding are discussed. Two strategies namely Focused Identification of Germplasm Strategy and Environmental Association Analysis that have been developed to effectively screen large germplasm collections for accessions with adaptive traits using geo-reference-based, mathematical modelling are highlighted. The potential applications of both approaches in mining available cotton landrace collections are also presented.

Looking for the needle in a downsized haystack: Whole-exome sequencing unravels genomic signals of climatic adaptation in Douglas-fir (Pseudotsuga menziesii).

Conifers often occur along steep gradients of diverse climates throughout their natural ranges, which is expected to result in spatially varying selection to local climate conditions. However, signals of climatic adaptation can often be confounded, because unraveled clines covary with signals caused by neutral evolutionary processes such as gene flow and genetic drift. Consequently, our understanding of how selection and gene flow have shaped phenotypic and genotypic differentiation in trees is still limited.A 40‐year‐old common garden experiment comprising 16 Douglas‐fir (Pseudotsuga menziesii) provenances from a north‐to‐south gradient of approx. 1,000 km was analyzed, and genomic information was obtained from exome capture, which resulted in an initial genomic dataset of >90,000 single nucleotide polymorphisms. We used a restrictive and conservative filtering approach, which permitted us to include only SNPs and individuals in environmental association analysis (EAA) that were free of potentially confounding effects (LD, relatedness among trees, heterozygosity deficiency, and deviations from Hardy–Weinberg proportions). We used four conceptually different genome scan methods based on FST outlier detection and gene–environment association in order to disentangle truly adaptive SNPs from neutral SNPs.We found that a relatively small proportion of the exome showed a truly adaptive signal (0.01%–0.17%) when population substructuring and multiple testing was accounted for. Nevertheless, the unraveled SNP candidates showed significant relationships with climate at provenance origins, which strongly suggests that they have featured adaptation in Douglas‐fir along a climatic gradient. Two SNPs were independently found by three of the employed algorithms, and one of them is in close proximity to an annotated gene involved in circadian clock control and photoperiodism as was similarly found in Populus balsamifera. Synthesis. We conclude that despite neutral evolutionary processes, phenotypic and genomic signals of adaptation to climate are responsible for differentiation, which in particular explain disparity between the well‐known coastal and interior varieties of Douglas‐fir.

Signatures of local adaptation to climate in natural populations of sweet chestnut (Castanea sativa Mill.) from southern Europe

• Key message Understanding the adaptive mechanisms of forest species is vital to ensure their survival in a climate change scenario. This study aimed at uncovering the relationship between genetic variability and environmental variables in natural Castanea sativa populations, unveiling how different climate scenarios drove local adaption processes using a landscape genomics approach. Our findings provide useful data for future management of this species. • Context Temperate forest species, such as chestnut (Castanea sativa Mill.), are currently threatened by increasing temperature together with disruption and reduction of precipitation due to climate change. In this context, understanding the adaptation processes of species will help to manage and ensure the conservation of forests.• Aims We studied the relationship between genetic variability and climate variables in natural populations of C. sativa using a landscape genomics approach aimed to identify local adaption processes.• Methods Using five genomic SSRs and eight functional EST-SSRs markers, 268 individuals belonging to ten different natural European chestnut populations distributed in contrasting climatic sites were genotyped. In addition, associations between allelic variation and climatic variables (environmental association analyses approach) were performed using Samβada and LFMM.• Results Results highlighted a strong inter-relationship between climate variables and evolutionary processes resulting in adaptive variation. STRUCTURE analysis based on functional markers split the populations in three separate gene pools (K = 3), mostly in agreement with the different climatic conditions existing in the studied areas. Divergent spatial patterns of genetic variation between rainy and arid areas were found. We detected a total of 202 associations with climate among 22 different alleles, 9% of which related with the outlier locus FIR059, known to be implicated in regulatory mechanisms during water stress adaptation processes.• Conclusion Landscape genomics analyses revealed a pattern of adaptive variation, where specific climatic variables influenced the frequencies distribution and fixation of several alleles, resulting in local adaptation processes of the populations in the investigated areas. Our findings underline the close inter-relationship existing between climate and genetic variability and indicate how this approach could provide valuable information for the management of forest species in a rapidly changing environment.

Little parallelism in genomic signatures of local adaptation in two sympatric, cryptic sister species.

Species living in sympatry and sharing a similar niche often express parallel phenotypes as a response to similar selection pressures. The degree of parallelism within underlying genomic levels is often unexplored, but can give insight into the mechanisms of natural selection and adaptation. Here, we use multi-dimensional genomic associations to assess the basis of local and climate adaptation in two sympatric, cryptic Crematogaster levior ant species along a climate gradient. Additionally, we investigate the genomic basis of chemical communication in both species. Communication in insects is mainly mediated by cuticular hydrocarbons (CHCs), which also protect against water loss and, hence, are subject to changes via environmental acclimation or adaptation. The combination of environmental and chemical association analyses based on genome-wide Pool-Seq data allowed us to identify single nucleotide polymorphisms (SNPs) associated with climate and with chemical differences. Within species, CHC changes as a response to climate seem to be driven by phenotypic plasticity, since there is no overlap between climate- and CHC-associated SNPs. The only exception is the odorant receptor OR22c, which may be a candidate for population-specific CHC recognition in one of the species. Within both species, climate is significantly correlated with CHC differences, as well as to allele frequency differences. However, associated candidate SNPs, genes and functions are largely species-specific and we find evidence for minimal parallel evolution only on the level of genomic regions (J=0.04). This highlights that even closely related species may follow divergent evolutionary trajectories when expressing similar adaptive phenotypes.

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-2200m), and one narrower, K. tristis (1600-2000m). We examine signatures of selection with respect to environmental and morphological variables in two mountain regions using FST 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.

Spatial and environmental influences on selection in a clock gene coding trinucleotide repeat in Canada lynx (Lynx canadensis).

Clock genes exhibit substantial control over gene expression and ultimately life-histories using external cues such as photoperiod, and are thus likely to be critical for adaptation to shifting seasonal conditions and novel environments as species redistribute their ranges under climate change. Coding trinucleotide repeats (cTNRs) are found within several clock genes, and may be interesting targets of selection due to their containment within exonic regions and elevated mutation rates. Here, we conduct inter-specific characterization of the NR1D1 cTNR between Canada lynx and bobcat, and intra-specific spatial and environmental association analyses of neutral microsatellites and our functional cTNR marker, to investigate the role of selection on this locus in Canada lynx. We report signatures of divergent selection between lynx and bobcat, with the potential for hybrid-mediated gene flow in the area of range overlap. We also provide evidence that this locus is under selection across Canada lynx in eastern Canada, with both spatial and environmental variables significantly contributing to the explained variation, after controlling for neutral population structure. These results suggest that cTNRs may play an important role in the generation of functional diversity within some mammal species, and allow for contemporary rates of adaptation in wild populations in response to environmental change. We encourage continued investment into the study of cTNR markers to better understand their broader relevance to the evolution and adaptation of mammals.

Individual-based population genomics reveal different drivers of adaptation in sympatric fish

Connectivity and local adaptation are two contrasting evolutionary forces highly influencing population structure. To evaluate the impact of early-life traits and environmental conditions on genetic structuring and adaptation, we studied two sympatric fish species in the Western Mediterranean Sea: Symphodus tinca and S. ocellatus. We followed an individual-based approach and measured early-life history traits from otolith readings, gathered information on environmental variables and obtained genome-wide markers from genotyping-by-sequencing (GBS). The two species presented contrasting population structure across the same geographic gradient, with high and significant population differentiation in S. ocellatus, mostly determined by oceanographic fronts, and low differentiation and no front effect in S. tinca. Despite their different levels of genetic differentiation, we identified in both species candidate regions for local adaptation by combining outlier analysis with environmental and phenotypic association analyses. Most candidate loci were associated to temperature and productivity in S. ocellatus and to temperature and turbulence in S. tinca suggesting that different drivers may determine genomic diversity and differentiation in each species. Globally, our study highlights that individual-based approach combining genomic, environmental and phenotypic information is key to identify signals of selection and the processes mediating them.

Disentangling the effects of geographic peripherality and habitat suitability on neutral and adaptive genetic variation in Swiss stone pine.

It is generally accepted that the spatial distribution of neutral genetic diversity within a species' native range mostly depends on effective population size, demographic history, and geographic position. However, it is unclear how genetic diversity at adaptive loci correlates with geographic peripherality or with habitat suitability within the ecological niche. Using exome-wide genomic data and distribution maps of the Alpine range, we first tested whether geographic peripherality correlates with four measures of population genetic diversity at >17,000 SNP loci in 24 Alpine populations (480 individuals) of Swiss stone pine (Pinus cembra) from Switzerland. To distinguish between neutral and adaptive SNP sets, we used four approaches (two gene diversity estimates, FST outlier test, and environmental association analysis) that search for signatures of selection. Second, we established ecological niche models for P. cembra in the study range and investigated how habitat suitability correlates with genetic diversity at neutral and adaptive loci. All estimates of neutral genetic diversity decreased with geographic peripherality, but were uncorrelated with habitat suitability. However, heterozygosity (He ) at adaptive loci based on Tajima's D declined significantly with increasingly suitable conditions. No other diversity estimates at adaptive loci were correlated with habitat suitability. Our findings suggest that populations at the edge of a species' geographic distribution harbour limited neutral genetic diversity due to demographic properties. Moreover, we argue that populations from suitable habitats went through strong selection processes, are thus well adapted to local conditions, and therefore exhibit reduced genetic diversity at adaptive loci compared to populations at niche margins.

Insights into the neutral and adaptive processes shaping the spatial distribution of genomic variation in the economically important Moroccan locust (Dociostaurus maroccanus).

Understanding the processes that shape neutral and adaptive genomic variation is a fundamental step to determine the demographic and evolutionary dynamics of pest species. Here, we use genomic data obtained via restriction site‐associated DNA sequencing to investigate the genetic structure of Moroccan locust (Dociostaurus maroccanus) populations from the westernmost portion of the species distribution (Iberian Peninsula and Canary Islands), infer demographic trends, and determine the role of neutral versus selective processes in shaping spatial patterns of genomic variation in this pest species of great economic importance. Our analyses showed that Iberian populations are characterized by high gene flow, whereas the highly isolated Canarian populations have experienced strong genetic drift and loss of genetic diversity. Historical demographic reconstructions revealed that all populations have passed through a substantial genetic bottleneck around the last glacial maximum (~21 ka BP) followed by a sharp demographic expansion at the onset of the Holocene, indicating increased effective population sizes during warm periods as expected from the thermophilic nature of the species. Genome scans and environmental association analyses identified several loci putatively under selection, suggesting that local adaptation processes in certain populations might not be impeded by widespread gene flow. Finally, all analyses showed few differences between outbreak and nonoutbreak populations. Integrated pest management practices should consider high population connectivity and the potential importance of local adaptation processes on population persistence.

Local adaptation across a complex bioclimatic landscape in two montane bumble bee species.

Understanding evolutionary responses to variation in temperature and precipitation across species ranges is of fundamental interest given ongoing climate change. The importance of temperature and precipitation for multiple aspects of bumble bee (Bombus) biology, combined with large geographic ranges that expose populations to diverse environmental pressures, make these insects well-suited for studying local adaptation. Here, we analyzed genome-wide sequence data from two widespread bumble bees, Bombus vosnesenskii and Bombus vancouverensis, using multiple environmental association analysis methods to investigate climate adaptation across latitude and altitude. The strongest signatures of selection were observed in B.vancouverensis, but despite unique responses between species for most loci, we detected several shared responses. Genes relating to neural and neuromuscular function and ion transport were especially evident with respect to temperature variables, while genes relating to cuticle formation, tracheal and respiratory system development, and homeostasis were associated with precipitation variables. Our data thus suggest that adaptive responses for tolerating abiotic variation are likely to be complex, but that several parallels among species can emerge even for these complex traits and landscapes. Results provide the framework for future work into mechanisms of thermal and desiccation tolerance in bumble bees and a set of genomic targets that might be monitored for future conservation efforts.

Fine-scale seascape genomics of an exploited marine species, the common cockle Cerastoderma edule, using a multimodelling approach.

Population dynamics of marine species that are sessile as adults are driven by oceanographic dispersal of larvae from spawning to nursery grounds. This is mediated by life‐history traits such as the timing and frequency of spawning, larval behaviour and duration, and settlement success. Here, we use 1725 single nucleotide polymorphisms (SNPs) to study the fine‐scale spatial genetic structure in the commercially important cockle species Cerastoderma edule and compare it to environmental variables and current‐mediated larval dispersal within a modelling framework. Hydrodynamic modelling employing the NEMO Atlantic Margin Model (AMM15) was used to simulate larval transport and estimate connectivity between populations during spawning months (April–September), factoring in larval duration and interannual variability of ocean currents. Results at neutral loci reveal the existence of three separate genetic clusters (mean F ST = 0.021) within a relatively fine spatial scale in the north‐west Atlantic. Environmental association analysis indicates that oceanographic currents and geographic proximity explain over 20% of the variance observed at neutral loci, while genetic variance (71%) at outlier loci was explained by sea surface temperature extremes. These results fill an important knowledge gap in the management of a commercially important and overexploited species, bringing us closer to understanding the role of larval dispersal in connecting populations at a fine geographic scale.

SLAF-seq Uncovers the Genetic Diversity and Adaptation of Chinese Elm (Ulmus parvifolia) in Eastern China

The Chinese elm is an important tree ecologically; however, little is known about its genetic diversity and adaptation mechanisms. In this study, a total of 107 individuals collected from seven natural populations in eastern China were investigated by specific locus amplified fragment sequencing (SLAF-seq). Based on the single nucleotide polymorphisms (SNPs) detected by SLAF-seq, genetic diversity and markers associated with climate variables were identified. All seven populations showed medium genetic diversity, with PIC values ranging from 0.2632 to 0.2761. AMOVA and Fst indicated that a low genetic differentiation existed among populations. Environmental association analyses with three climate variables (annual rainfall, annual average temperature, and altitude) resulted in, altogether, 43 and 30 putative adaptive loci by Bayenv2 and LFMM, respectively. Five adaptive genes were annotated, which were related to the functions of glycosylation, peroxisome synthesis, nucleic acid metabolism, energy metabolism, and signaling. This study was the first on the genetic diversity and local adaptation in Chinese elms, and the results will be helpful in future work on molecular breeding.

Minimal genetic differentiation of the malaria vector Nyssorhynchus darlingi associated with forest cover level in Amazonian Brazil.

The relationship between deforestation and malaria in Amazonian Brazil is complex, and a deeper understanding of this relationship is required to inform effective control measures in this region. Here, we are particularly interested in characterizing the impact of land use and land cover change on the genetics of the major regional vector of malaria, Nyssorhynchus darlingi (Root). We used nextera-tagmented, Reductively Amplified DNA (nextRAD) genotyping-by-sequencing to genotype 164 Ny. darlingi collected from 16 collection sites with divergent forest cover levels in seven municipalities in four municipality groups that span the state of Amazonas in northwestern Amazonian Brazil: São Gabriel da Cachoeira, Presidente Figueiredo, four municipalities in the area around Cruzeiro do Sul, and Lábrea. Using a dataset of 5,561 Single Nucleotide Polymorphisms (SNPs), we investigated the genetic structure of these Ny. darlingi populations with a combination of model- and non-model-based analyses. We identified weak to moderate genetic differentiation among the four municipality groups. There was no evidence for microgeographic genetic structure of Ny. darlingi among forest cover levels within the municipality groups, indicating that there may be gene flow across areas of these municipalities with different degrees of deforestation. Additionally, we conducted an environmental association analysis using two outlier detection methods to determine whether individual SNPs were associated with forest cover level without affecting overall population genetic structure. We identified 14 outlier SNPs, and investigated functions associated with their proximal genes, which could be further characterized in future studies.

Detection of environmental and morphological adaptation despite high landscape genetic connectivity in a pest grasshopper (Phaulacridium vittatum).

Widespread species that exhibit both high gene flow and the capacity to occupy heterogeneous environments make excellent models for examining local selection processes along environmental gradients. Here we evaluate the influence of temperature and landscape variables on genetic connectivity and signatures of local adaptation in Phaulacridium vittatum, a widespread agricultural pest grasshopper, endemic to Australia. With sampling across a 900-km latitudinal gradient, we genotyped 185 P.vittatum from 19 sites at 11,408 single nucleotide polymorphisms (SNPs) using ddRAD sequencing. Despite high gene flow across sites (pairwise FST =0.0003-0.08), landscape genetic resistance modelling identified a positive nonlinear effect of mean annual temperature on genetic connectivity. Urban areas and water bodies had a greater influence on genetic distance among sites than pasture, agricultural areas and forest. Together, FST outlier tests and environmental association analysis (EAA) detected 242 unique SNPs under putative selection, with the highest numbers associated with latitude, mean annual temperature and body size. A combination of landscape genetic connectivity analysis together with EAA identified mean annual temperature as a key driver of both neutral gene flow and environmental selection processes. Gene annotation of putatively adaptive SNPs matched with gene functions for olfaction, metabolic detoxification and ultraviolet light shielding. Our results imply that this widespread agricultural pest has the potential to spread and adapt under shifting temperature regimes and land cover change.

Convergent evolution in Arabidopsis halleri and Arabidopsis arenosa on calamine metalliferous soils.

It is a plausible hypothesis that parallel adaptation events to the same environmental challenge should result in genetic changes of similar or identical effects, depending on the underlying fitness landscapes. However, systematic testing of this is scarce. Here we examine this hypothesis in two closely related plant species, Arabidopsis halleri and Arabidopsis arenosa, which co-occur at two calamine metalliferous (M) sites harbouring toxic levels of the heavy metals zinc and cadmium. We conduct individual genome resequencing alongside soil elemental analysis for 64 plants from eight populations on M and non-metalliferous (NM) soils, and identify genomic footprints of selection and local adaptation. Selective sweep and environmental association analyses indicate a modest degree of gene as well as functional network convergence, whereby the proximal molecular factors mediating this convergence mostly differ between site pairs and species. Notably, we observe repeated selection on identical single nucleotide polymorphisms in several A. halleri genes at two independently colonized M sites. Our data suggest that species-specific metal handling and other biological features could explain a low degree of convergence between species. The parallel establishment of plant populations on calamine M soils involves convergent evolution, which will probably be more pervasive across sites purposely chosen for maximal similarity in soil composition.This article is part of the theme issue ‘Convergent evolution in the genomics era: new insights and directions’.

Detecting the genetic basis of local adaptation in loblolly pine (Pinus taeda L.) using whole exome-wide genotyping and an integrative landscape genomics analysis approach.

In the Southern United States, the widely distributed loblolly pine contributes greatly to lumber and pulp production, as well as providing many important ecosystem services. Climate change may affect the productivity and range of loblolly pine. Nevertheless, we have insufficient knowledge of the adaptive potential and the genetics underlying the adaptability of loblolly pine. To address this, we tested the association of 2.8 million whole exome‐based single nucleotide polymorphisms (SNPs) with climate and geographic variables, including temperature, precipitation, latitude, longitude, and elevation data. Using an integrative landscape genomics approach by combining multiple environmental association and outlier detection analyses, we identified 611 SNPs associated with 56 climate and geographic variables. Longitude, maximum temperature of the warm months and monthly precipitation associated with most SNPs, indicating their importance and complexity in shaping the genetic variation in loblolly pine. Functions of candidate genes related to terpenoid synthesis, pathogen defense, transcription factors, and abiotic stress response. We provided evidence that environment‐associated SNPs also composed the genetic structure of adaptive phenotypic traits including height, diameter, metabolite levels, and gene transcript abundance. Our study promotes understanding of the genetic basis of local adaptation in loblolly pine and provides promising tools for selecting genotypes adapted to local environments in a changing climate.