Genetic divergence and isolation by thermal environment in geothermal populations of an aquatic invertebrate

Temperature is one of the most influential forces of natural selection impacting all biological levels. In the face of increasing global temperatures, studies over small geographic scales allowing investigations on the effects of gene flow are of great value for understanding thermal adaptation. Here, we investigated genetic population structure in the freshwater gastropod Radix balthica originating from contrasting thermal habitats in three areas of geothermal activity in Iceland. Snails from 32 sites were genotyped at 208 AFLP loci. Five AFLPs were identified as putatively under divergent selection in Lake Mývatn, a geothermal lake with an almost 20 °C difference in mean temperature across a distance of a few kilometres. In four of these loci, variation across all study populations was correlated with temperature. We found significant population structure in neutral markers both within and between the areas. Cluster analysis using neutral markers classified the sites mainly by geography, whereas analyses using markers under selection differentiated the sites based on temperature. Isolation by distance was stronger in the neutral than in the outlier loci. Pairwise differences based on outlier F_ST were significantly correlated with temperature at different spatial scales, even after correcting for geographic distance or neutral pairwise F_ST differences. In general, genetic variation decreased with increasing environmental temperature, possibly suggesting that natural selection had reduced the genetic diversity in the warm origin sites. Our results emphasize the influence of environmental temperature on the genetic structure of populations and suggest local thermal adaptation in these geothermal habitats.     

Ecophysiological responses of nine floodplain meadow species to changing hydrological conditions

Here, we investigated how species from different floodplain meadow plant communities differ in their ecophysiological responses to an abrupt change in hydrological conditions. We simulated two contrasting hydrological scenarios for 5 weeks under controlled conditions: the waterlogging of a mesophilic species mixture (flooding scenario) and the drying of a hygrophilic species mixture (drying scenario). The mesophilic mixture was composed of three characteristic species of dry habitats (dry species) and three indifferent species with regard to the hydrological conditions; the hygrophilic mixture was composed of three characteristic species of wet habitats (wet species) and the same three indifferent species. The flooding scenario induced a significant decrease in photochemical efficiency (F _v/F _m) of all species for at least one day, and four of these species did not readjust their F _v/F _m at the end of the experiment. The F _v/F _m of four species was significantly reduced by the drying scenario, but was readjusted within the time of the experiment. Significant leaf plastic responses were mainly detected under the drying scenario, through increasing specific leaf area (SLA) and decreasing leaf dry matter content (LDMC) for two indifferent species, and through decreasing SLA and increasing LDMC for two wet species. Finally, the flooding scenario significantly decreased the aboveground biomass of one dry species, but increased the aboveground biomass of the three indifferent species, suggesting an acclimation response to waterlogging through increasing shoot/root ratio. Since the dry species did not show ecophysiological acclimation to changing hydrological conditions (in contrast to the wet species and the indifferent species), this study demonstrated that mesophilic meadow communities should be particularly affected by a change from dry to wet conditions. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

MULTIVARIATE ADAPTATION BUT NO INCREASE IN COMPETITIVE ABILITY IN INVASIVE GERANIUM CAROLINIANUM L. (GERANIACEAE)

Adaptive evolution can affect the successful establishment of invasive species, but changes in selective pressures, loss of genetic variation in relevant traits, and/or altered trait correlations can make adaptation difficult to predict. We used a common‐garden experiment to assess trait correlations and patterns of adaptation in the invasive plant, Geranium carolinianum, sampled across 20 populations in its native (United States) and invasive (China) ranges. We used multivariate Q_ST – F_ST tests to determine if phenotypic differences between countries are attributable to adaptation. We also compared population‐level variation within each country to assess whether local adaptation resulted in similar multivariate phenotypes in the United States and China. Between countries, most phenotypic differences are indistinguishable from genetic drift, although we detected a signature of adaptation to the colder, drier winters in China. There was no evidence for increases in invasive traits in China. Within countries, strong multivariate adaptation appears to be driven by latitudinal climatic variation in the United States, but not in China. Additionally, adaptive trait combinations as well as their underlying correlations differ between the two countries, indicating that adaptation in invasive populations does not parallel patterns in native populations due to differences in selection pressures, genetic constraints, or both.     

Genetic diversity of Couratari multiflora and Couratari guianensis (Lecythidaceae): consequences of two types of rarity in central Amazonia

We quantified the within-population genetic variation of Couratari multiflora and C. guianensis, two tree species found in terra firme forests of central Amazonia. Both species have some ecological features in common, but they differ in population abundance across their geographic ranges. While C. multiflora has been found only in low-density populations in all sites studied to date, C. guianensis is relatively common in some sites and very scarce in others. In a 400-ha plot, we found 41 and 29 adults of C. multiflora and C. guianensis, respectively. Twenty-two saplings of C. guianensis and 103 seedlings of C. multiflora were also examined. The mean expected heterozygosities (H_em) of seedlings and adults of C. multiflora were 0.431 and 0.436, and the mean fixation indices (F_m), 0.114 and 0.176, respectively. For C. guianensis, saplings and adults presented H_em equal to 0.425 and 0.429, and the F_m were 0.393 and 0.527, respectively. These low-density populations of two congeneric species did not differ in terms of genetic diversity, but rather they differed in terms of mean observed heterozygosity (H_om), and therefore F_m. The species with variable population density had lower H_om and greater F_m relative to the species that is always found in low-density.     

Long-Term Waterlogging: Nutrient,Gas Exchange,Photochemical and Pigment Characteristics of Eucalyptus nitens Saplings

Imposition of waterlogging for eight months induced morphological adaptation in Eucalyptus nitens in the form of adventitious and aerenchymatous roots and hypertrophy of stems. Foliar calcium (1.3-fold), potassium (2-fold) and phosphorus (2.4-fold) were lower and iron (5.6-fold) was higher in waterlogged than control saplings. Stem Ca (1.7-fold) was lower, whereas Mn (1.8-fold) and Fe (117-fold) were higher in waterlogged than control saplings. Distinct purple pigmentation was observed in xylem tissues of waterlogged saplings. A significant reduction in the maximum photosynthetic rate and photochemical efficiency was observed in waterlogged compared to control saplings. Although chlorophyll levels were similar, the xanthophyll cycle pool size was significantly greater in waterlogged saplings and may have contributed to the relatively greater capacity for light energy dissipation observed. Predawn xanthophyll cycle engagement was significantly greater in waterlogged than control saplings, despite relatively mild temperatures. Foliar anthocyanin concentration was higher in waterlogged than control saplings.     

Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.     

Climate‐related genetic variation in drought‐resistance of Douglas‐fir (Pseudotsuga menziesii)

There is a general assumption that intraspecific populations originating from relatively arid climates will be better adapted to cope with the expected increase in drought from climate change. For ecologically and economically important species, more comprehensive, genecological studies that utilize large distributions of populations and direct measures of traits associated with drought‐resistance are needed to empirically support this assumption because of the implications for the natural or assisted regeneration of species. We conducted a space‐for‐time substitution, common garden experiment with 35 populations of coast Douglas‐fir (Pseudotsuga menziesii var. menziesii) growing at three test sites with distinct summer temperature and precipitation (referred to as ‘cool/moist’, ‘moderate’, or ‘warm/dry’) to test the hypotheses that (i) there is large genetic variation among populations and regions in traits associated with drought‐resistance, (ii) the patterns of genetic variation are related to the native source‐climate of each population, in particular with summer temperature and precipitation, (iii) the differences among populations and relationships with climate are stronger at the warm/dry test site owing to greater expression of drought‐resistance traits (i.e., a genotype × environment interaction). During midsummer 2012, we measured the rate of water loss after stomatal closure (transpiration_min), water deficit (% below turgid saturation), and specific leaf area (SLA, cm² g^?1) on new growth of sapling branches. There was significant genetic variation in all plant traits, with populations originating from warmer and drier climates having greater drought‐resistance (i.e., lower transpiration_min, water deficit and SLA), but these trends were most clearly expressed only at the warm/dry test site. Contrary to expectations, populations from cooler climates also had greater drought‐resistance across all test sites. Multiple regression analysis indicated that Douglas‐fir populations from regions with relatively cool winters and arid summers may be most adapted to cope with drought conditions that are expected in the future.     

The role of selection and historical factors in driving population differentiation along an elevational gradient in an island bird

Adaptation to local environmental conditions and the range dynamics of populations can influence evolutionary divergence along environmental gradients. Thus, it is important to investigate patterns of both phenotypic and genetic variations among populations to reveal the respective roles of these two types of factors in driving population differentiation. Here, we test for evidence of phenotypic and genetic structure across populations of a passerine bird (Zosterops borbonicus) distributed along a steep elevational gradient on the island of Réunion. Using 11 microsatellite loci screened in 401 individuals from 18 localities distributed along the gradient, we found that genetic differentiation occurred at two spatial levels: (i) between two main population groups corresponding to highland and lowland areas, respectively, and (ii) within each of these two groups. In contrast, several morphological traits varied gradually along the gradient. Comparison of neutral genetic differentiation (F_ST) and phenotypic differentiation (P_ST) showed that P_ST largely exceeds F_ST at several morphological traits, which is consistent with a role for local adaptation in driving morphological divergence along the gradient. Overall, our results revealed an area of secondary contact midway up the gradient between two major, cryptic, population groups likely diverged in allopatry. Remarkably, local adaptation has shaped phenotypic differentiation irrespective of population history, resulting in different patterns of variation along the elevational gradient. Our findings underscore the importance of understanding both historical and selective factors when trying to explain variation along environmental gradients.     

Selection on life‐history traits and genetic population divergence in rotifers

A combination of founder effects and local adaptation – the Monopolization hypothesis – has been proposed to reconcile the strong population differentiation of zooplankton dwelling in ponds and lakes and their high dispersal abilities. The role genetic drift plays in genetic differentiation of zooplankton is well documented, but the impact of natural selection has received less attention. Here, we compare differentiation in neutral genetic markers (F_ST) and in quantitative traits (Q_ST) in six natural populations of the rotifer Brachionus plicatilis to assess the importance of natural selection in explaining genetic differentiation of life‐history traits. Five life‐history traits were measured in four temperature × salinity combinations in common‐garden experiments. Population differentiation for neutral genetic markers – 11 microsatellite loci – was very high (F_ST = 0.482). Differentiation in life‐history traits was higher in traits related to sexual reproduction than in those related to asexual reproduction. Q_ST values for diapausing egg production (a trait related to sexual reproduction) were higher than their corresponding F_ST in some pairs of populations. Our results indicate the importance of divergent natural selection in these populations and suggest local adaptation to the unpredictability of B. plicatilis habitats.     

Climate change shifts natural selection and the adaptive potential of the perennial forb Boechera stricta in the Rocky Mountains

Heritable genetic variation is necessary for populations to evolve in response to anthropogenic climate change. However, antagonistic genetic correlations among traits may constrain the rate of adaptation, even if substantial genetic variation exists. We examine potential genetic responses to selection by comparing multivariate genetic variance–covariances of traits and fitness (multivariate Robertson–Price identities) across different environments in a reciprocal transplant experiment of the forb Boechera stricta in the Rocky Mountains. By transplanting populations into four common gardens arrayed along an elevational gradient, and exposing populations to control and snow removal treatments, we simulated future and current climates and snowmelt regimes. Genetic variation in flowering and germination phenology declined in plants moved downslope to warmer, drier sites, suggesting that these traits may have a limited ability to evolve under future climates. Simulated climate change via snow removal altered the strength of selection on flowering traits, but we found little evidence that genetic correlations among traits are likely to affect the rate of adaptation to climate change. Overall, our results suggest that climate change may alter the evolutionary potential of B. stricta, but reduced expression of genetic variation may be a larger impediment to adaptation than constraints imposed by antagonistic genetic correlations.     

Is local selection so widespread in river organisms? Fractal geometry of river networks leads to high bias in outlier detection

Identifying local adaptation is crucial in conservation biology to define ecotypes and establish management guidelines. Local adaptation is often inferred from the detection of loci showing a high differentiation between populations, the so‐called F_ST outliers. Methods of detection of loci under selection are reputed to be robust in most spatial population models. However, using simulations we showed that F_ST outlier tests provided a high rate of false‐positives (up to 60%) in fractal environments such as river networks. Surprisingly, the number of sampled demes was correlated with parameters of population genetic structure, such as the variance of F_STs, and hence strongly influenced the rate of outliers. This unappreciated property of river networks therefore needs to be accounted for in genetic studies on adaptation and conservation of river organisms.     

LOCAL ADAPTATION,PHENOTYPIC DIFFERENTIATION,AND HYBRID FITNESS IN DIVERGED NATURAL POPULATIONS OF ARABIDOPSIS LYRATA

Selection for local adaptation results in genetic differentiation in ecologically important traits. In a perennial, outcrossing model plant Arabidopsis lyrata, several differentiated phenotypic traits contribute to local adaptation, as demonstrated by fitness advantage of the local population at each site in reciprocal transplant experiments. Here we compared fitness components, hierarchical total fitness and differentiation in putatively ecologically important traits of plants from two diverged parental populations from different continents in the native climate conditions of the populations in Norway and in North Carolina (NC, U.S.A.). Survival and number of fruits per inflorescence indicated local advantage at both sites and aster life‐history models provided additional evidence for local adaptation also at the level of hierarchical total fitness. Populations were also differentiated in flowering start date and floral display. We also included reciprocal experimental F₁ and F₂ hybrids to examine the genetic basis of adaptation. Surprisingly, the F₂ hybrids showed heterosis at the study site in Norway, likely because of a combination of beneficial dominance effects from different traits. At the NC site, hybrid fitness was mostly intermediate relative to the parental populations. Local cytoplasmic origin was associated with higher fitness, indicating that cytoplasmic genomes also may contribute to the evolution of local adaptation.     

Waterlogging responses in dune,swale and marsh populations of Spartina patens under field conditions

Summary Soil waterlogging responses were examined in three Spartina patens populations along a steep flooding gradient in coastal Louisiana. Root anatomy and physiological indicators of anaerobic metabolism were examined to identify and compare flooding responses in dune, swale and marsh populations, while soil physicochemical factors were measured to characterize the three habitats. Soil waterlogging increased along the gradient from dune to marsh habitats and was accompanied by increases in root porosity (aerenchyma). Aerenchyma in marsh roots was apparently insufficient to provide enough oxygen for aerobic respiratory demand, as indicated by high root alcohol dehydrogenase activities and low energy charge ratios. Patterns of root metabolic indicators suggest that dune and swale roots generally respired aerobically, while anaerobic metabolism was important in marsh roots. However, in each population, relatively greater soil waterloging was accompanied by differences in enzyme activities leading to malate accumulation. In dune and swale roots under these circumstances, depressed adenylate energy charge ratios may have been the result of an absence of increased ethanol fermentation. These trends suggest that: 1) Aerenchyma formation was an important, albeit incomplete, long-term adaptation to the prevalent degree of soil waterlogging. 2) All populations adjusted root metabolism in response to a relative (short-term) increase in soil waterlogging.     

Local adaptation with high gene flow: temperature parameters drive adaptation to altitude in the common frog (Rana temporaria)

Both environmental and genetic influences can result in phenotypic variation. Quantifying the relative contributions of local adaptation and phenotypic plasticity to phenotypes is key to understanding the effect of environmental variation on populations. Identifying the selective pressures that drive divergence is an important, but often lacking, next step. High gene flow between high‐ and low‐altitude common frog (Rana temporaria) breeding sites has previously been demonstrated in Scotland. The aim of this study was to assess whether local adaptation occurs in the face of high gene flow and to identify potential environmental selection pressures that drive adaptation. Phenotypic variation in larval traits was quantified in R. temporaria from paired high‐ and low‐altitude sites using three common temperature treatments. Local adaptation was assessed using Q_ST–F_ST analyses, and quantitative phenotypic divergence was related to environmental parameters using Mantel tests. Although evidence of local adaptation was found for all traits measured, only variation in larval period and growth rate was consistent with adaptation to altitude. Moreover, this was only evident in the three mountains with the highest high‐altitude sites. This variation was correlated with mean summer and winter temperatures, suggesting that temperature parameters are potentially strong selective pressures maintaining local adaptation, despite high gene flow.     

Genome scan in the mosquito Aedes rusticus: population structure and detection of positive selection after insecticide treatment

Identification of genes involved in local adaptation is particularly challenging for species functioning as a network of interconnected populations undergoing frequent extinctions–recolonizations, because populations are submitted to contrasted evolutionary pressures. Using amplified fragment length polymorphism markers, population genetic structure of the mosquito Aedes rusticus was analysed in five geographical areas of the French Rhône‐Alpes region. We included a number of sites that were treated with the bio‐insecticide Bacillus thuringiensis israelensis (Bti) for more than 15 years. Analysis of molecular variance revealed that most of the genetic variability was found within populations (96%), with no significant variation among geographical areas, although variation among populations within areas (4%) was significant. The global genetic differentiation index F_ST was low (0.0366 ± 0.167). However, pairwise F_ST values were significant and no isolation‐by‐distance at the regional level was observed, suggesting a metapopulation structure in this species. Bti‐treatment had no effect on genetic structure and on within‐population genetic diversity. Potential signatures of positive selection associated with Bti‐treatment were detected for five loci, even though toxicological bioassays performed on field‐collected larvae showed no significant difference in mortality between Bti‐treated and nontreated sites. The difficulty of detecting moderate resistance in field‐collected larvae together with possible differential persistence of toxins in the environment may explain our inability to detect a toxicological response to Bti in treated sites. The evidence for positive selection occurring at several genomic regions suggests a first step towards Bti resistance in natural mosquito populations treated with this bio‐insecticide. Furthermore, this signal was detectable using genomic tools before any toxicological evidence for resistance could be identified.     

GENETIC FINGERPRINT-INFERRED POPULATION SUBDIVISION AND SPATIAL GENETIC TESTS FOR ISOLATION BY DISTANCE AND ADAPTATION IN THE COASTAL PLANT LIMONIUM CAROLINIANUM

This paper examines two wild populations of Limonium carolinianum for population genetic subdivision and spatial patterns of genetic variation in an attempt to simultaneously test for both the action of local adaptation to tidal gradients and isolation by distance (IBD). A VNTR (variable number of tandem repeats) genetic “fingerprinting” marker was used to infer relatedness among mapped plants in two populations. Band sharing within and between populations estimated F'_ST, an approximate measure of F_ST. Regression models were used to analyze the relationship between band sharing and spatial separation in tidal elevation and horizontal distance, as well as the relationship of fecundity differences with band sharing and spatial distance. Populations differed in band size frequency distributions and mean number of bands per profile and, therefore, likely differed in effective population size. F'_ST was estimated at 0.0678 and was significantly greater than F'_ST among randomly constructed subpopulations. Band sharing decreased 0.13% per meter in one population but showed no significant relation to distance in the other. In the population with significant IBD band sharing increased with increasingly different tidal elevation, contrary to an adaptive hypothesis, possibly due to directional gene flow or drift. Deme sizes were approximately 25 meters and greater than 100 meters, spanning larger areas than the entire environmental gradient. Fecundity differences were not associated with spatial parameters or band sharing. Unequal potential maternal fecundity measured as variance in number of seeds per maternal family was a significant source of genetic sampling variance. The VNTR marker employed is capable of detecting adaptation as identity by descent in ecological time and is an appropriate method for estimating the net evolutionary fate of polygenic traits. The results show that the net balance between selection along an environmental gradient and the effects of IBD and unequal maternal fecundity favor genetic differentiation by random processes in populations of Limonium.     

Inferring demographic processes from the genetic structure of a metapopulation of <Emphasis Type="Italic">Boltonia decurrens</Emphasis>(Asteraceae)

Boltonia decurrens(Asteraceae), a federally listed, threatened floodplain species, requires regular flooding for suitable habitat and seed dispersal. Flood suppression and habitat destruction have resulted in fewer than 25 populations remaining throughout its 400 km range. Because individual populations are widely interspaced (>10 km) and subject to frequent extinction and colonization, seed dispersal along the river, not pollen flow, is likely the primary determinant of population genetic structure. We used neutral genetic markers (isozymes) assayed for fourteen populations to determine which demographic processes contribute to the genetic structure of B. decurrens. Significant genetic differentiation was detected among populations (F _ST=0.098, P< 0.05) but not among regions (F _RT=0.013, P> 0.05), suggesting that long-distance dispersal events occur and involve seed from a small number of populations. Correspondingly, we found no evidence of isolation by distance, and admixture analyses indicate that colonization events involve seed from 3 to 5 source populations. Individual populations exhibited high levels of fixation (mean F _IS=0.192, P< 0.05), yet mean population outcrossing rates were high (t _m=0.87–0.95) and spatial autocorrelation analyses revealed no fine-scale within population structure, indicating that inbreeding alone cannot explain the observed fixation. Rather, genetic bottlenecks, detected for 12 of 14 populations, and admixture at population founding may be important sources of fixation. These observations are consistent with a metapopulation model and confirm the importance of regular flooding events, capable of producing suitable habitat and dispersing seed long distances, to the long-term persistence of B. decurrens.     

Postzygotic isolation drives genomic speciation between highly cryptic Hypocnemis antbirds from Amazonia

How species evolve reproductive isolation in the species-rich Amazon basin is poorly understood in vertebrates. Here, we sequenced a reference genome and used a genome-wide sample of SNPs to analyze a hybrid zone between two highly cryptic species of Hypocnemis warbling-antbirds—the Rondonia warbling-antbird (H. ochrogyna) and Spix's warbling-antbird (H. striata)—in a headwater region of southern Amazonia. We found that both species commonly hybridize, producing F₁s and a variety of backcrosses with each species but we detected only one F₂-like hybrid. Patterns of heterozygosity, hybrid index, and interchromosomal linkage disequilibrium in hybrid populations closely match expectations under strong postzygotic isolation. Hybrid zone width (15.4 km) was much narrower than expected (211 km) indicating strong selection against hybrids. A remarkably high degree of concordance in cline centers and widths across loci, and a lack of reduced interspecific F_st between populations close to versus far from the contact zone, suggest that genetic incompatibilities have rendered most of the genome immune to introgression. These results support intrinsic postzygotic isolation as a driver of speciation in a moderately young cryptic species pair from the Amazon and suggest that species richness of the Amazon may be grossly underestimated.     

Diversity and genetic connectivity among populations of a threatened tree (<Emphasis Type="Italic">Dalbergia nigra</Emphasis>) in a recently fragmented landscape of the Brazilian Atlantic Forest

In this study we evaluated the influence of recent landscape fragmentation on the dynamics of remnant fragments from the Brazilian Atlantic Forest. This biome is one of the richest in the world and has been extensively deforested and fragmented. We sampled five populations of the threatened Dalbergia nigra, a tree endemic to the Brazilian Atlantic Forest, two located in a large reserve of continuous forest and three in fragments of different sizes and levels of disturbance. In order to assess historical changes, considering the longevity of the analyzed species, 119 adults and 116 saplings were genotyped for six microsatellite loci. Lower levels of genetic diversity were found in the most impacted fragments when compared to the most preserved population located inside the reserve, and there was significant genetic structure among the populations studied (pairwise F _ST = 0.031–0.152; pairwise D _EST = 0.039–0.301). However, genetic structure among saplings (F _ST = 0.056; D _EST = 0.231) was significantly lower than among adults (F _ST = 0.088; D _EST = 0.275). Estimates of contemporary gene flow based on assignment tests corroborated this result, suggesting that fragmentation led to an increase in gene flow. This connectivity among remnant fragments could mitigate the loss of genetic diversity through a metapopulation dynamic, but the high rate of habitat loss and the unknown long-term genetic effects add uncertainty. These results, taken together with the presence of private alleles in disturbed populations, highlight the importance of preserving the extant fragments.     

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 F_ST = 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, F_ST 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.