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.     

Adaptive and neutral markers both show continent‐wide population structure of mountain pine beetle (Dendroctonus ponderosae)

Assessments of population genetic structure and demographic history have traditionally been based on neutral markers while explicitly excluding adaptive markers. In this study, we compared the utility of putatively adaptive and neutral single‐nucleotide polymorphisms (SNPs) for inferring mountain pine beetle population structure across its geographic range. Both adaptive and neutral SNPs, and their combination, allowed range‐wide structure to be distinguished and delimited a population that has recently undergone range expansion across northern British Columbia and Alberta. Using an equal number of both adaptive and neutral SNPs revealed that adaptive SNPs resulted in a stronger correlation between sampled populations and inferred clustering. Our results suggest that adaptive SNPs should not be excluded prior to analysis from neutral SNPs as a combination of both marker sets resulted in better resolution of genetic differentiation between populations than either marker set alone. These results demonstrate the utility of adaptive loci for resolving population genetic structure in a nonmodel organism.     

The population genomic signature of environmental association and gene flow in an ecologically divergent tree species Metrosideros polymorpha (Myrtaceae)

Genomewide markers enable us to study genetic differentiation within a species and the factors underlying it at a much higher resolution than before, which advances our understanding of adaptation in organisms. We investigated genomic divergence in Metrosideros polymorpha, a woody species that occupies a wide range of ecological habitats across the Hawaiian Islands and shows remarkable phenotypic variation. Using 1659 single nucleotide polymorphism (SNP) markers annotated with the genome assembly, we examined the population genetic structure and demographic history of nine populations across five elevations and two ages of substrates on Mauna Loa, the island of Hawaii. The nine populations were differentiated into two genetic clusters distributed on the lower and higher elevations and were largely admixed on the middle elevation. Demographic modelling revealed that the two genetic clusters have been maintained in the face of gene flow, and the effective population size of the high‐altitude cluster was much smaller. A F_ST‐based outlier search among the 1659 SNPs revealed that 34 SNPs (2.05%) were likely to be under divergent selection and the allele frequencies of 21 of them were associated with environmental changes along elevations, such as temperature and precipitation. This study shows a genomic mosaic of M. polymorpha, in which contrasting divergence patterns were found. While most genomic polymorphisms were shared among populations, a small fraction of the genome was significantly differentiated between populations in diverse environments and could be responsible for the dramatic adaptation to a wide range of environments.     

Coast to coast: High genomic connectivity in North American scoters

Dispersal shapes demographic processes and therefore is fundamental to understanding biological, ecological, and evolutionary processes acting within populations. However, assessing population connectivity in scoters (Melanitta sp.) is challenging as these species have large spatial distributions that span remote landscapes, have varying nesting distributions (disjunct vs. continuous), exhibit unknown levels of dispersal, and vary in the timing of the formation of pair bonds (winter vs. fall/spring migration) that may influence the distribution of genetic diversity. Here, we used double‐digest restriction‐associated DNA sequence (ddRAD) and microsatellite genotype data to assess population structure within the three North American species of scoter (black scoter, M. americana; white‐winged scoter, M. deglandi; surf scoter, M. perspicillata), and between their European congeners (common scoter, M. nigra; velvet scoter, M. fusca). We uncovered no or weak genomic structure (ddRAD Φ_ST < 0.019; microsatellite F_ST < 0.004) within North America but high levels of structure among European congeners (ddRAD Φ_ST > 0.155, microsatellite F_ST > 0.086). The pattern of limited genomic structure within North America is shared with other sea duck species and is often attributed to male‐biased dispersal. Further, migratory tendencies (east vs. west) of female surf and white‐winged scoters in central Canada are known to vary across years, providing additional opportunities for intracontinental dispersal and a mechanism for the maintenance of genomic connectivity across North America. In contrast, the black scoter had relatively elevated levels of divergence between Alaska and Atlantic sites and a second genetic cluster found in Alaska at ddRAD loci was concordant with its disjunct breeding distribution suggestive of a dispersal barrier (behavioral or physical). Although scoter populations appear to be connected through a dispersal network, a small percentage (<4%) of ddRAD loci had elevated divergence which may be useful in linking areas (nesting, molting, staging, and wintering) throughout the annual cycle.     

Climate‐associated genetic variation in Fagus sylvatica and potential responses to climate change in the French Alps

Local adaptation patterns have been found in many plants and animals, highlighting the genetic heterogeneity of species along their range of distribution. In the next decades, global warming is predicted to induce a change in the selective pressures that drive this adaptive variation, forcing a reshuffling of the underlying adaptive allele distributions. For species with low dispersion capacity and long generation time such as trees, the rapidity of the change could impede the migration of beneficial alleles and lower their capacity to track the changing environment. Identifying the main selective pressures driving the adaptive genetic variation is thus necessary when investigating species capacity to respond to global warming. In this study, we investigate the adaptive landscape of Fagus sylvatica along a gradient of populations in the French Alps. Using a double‐digest restriction‐site‐associated DNA (ddRAD) sequencing approach, we identified 7,000 SNPs from 570 individuals across 36 different sites. A redundancy analysis (RDA)‐derived method allowed us to identify several SNPs that were strongly associated with climatic gradients; moreover, we defined the primary selective gradients along the natural populations of F. sylvatica in the Alps. Strong effects of elevation and humidity, which contrast north‐western and south‐eastern site, were found and were believed to be important drivers of genetic adaptation. Finally, simulations of future genetic landscapes that used these findings allowed identifying populations at risk for F. sylvatica in the Alps, which could be helpful for future management plans.     

Assessing polar bear (Ursus maritimus) population structure in the Hudson Bay region using SNPs

Defining subpopulations using genetics has traditionally used data from microsatellite markers to investigate population structure; however, single‐nucleotide polymorphisms (SNPs) have emerged as a tool for detection of fine‐scale structure. In Hudson Bay, Canada, three polar bear (Ursus maritimus) subpopulations (Foxe Basin (FB), Southern Hudson Bay (SH), and Western Hudson Bay (WH)) have been delineated based on mark–recapture studies, radiotelemetry and satellite telemetry, return of marked animals in the subsistence harvest, and population genetics using microsatellites. We used SNPs to detect fine‐scale population structure in polar bears from the Hudson Bay region and compared our results to the current designations using 414 individuals genotyped at 2,603 SNPs. Analyses based on discriminant analysis of principal components (DAPC) and STRUCTURE support the presence of four genetic clusters: (i) Western—including individuals sampled in WH, SH (excluding Akimiski Island in James Bay), and southern FB (south of Southampton Island); (ii) Northern—individuals sampled in northern FB (Baffin Island) and Davis Strait (DS) (Labrador coast); (iii) Southeast—individuals from SH (Akimiski Island in James Bay); and (iv) Northeast—individuals from DS (Baffin Island). Population structure differed from microsatellite studies and current management designations demonstrating the value of using SNPs for fine‐scale population delineation in polar bears.     

Radiation of the polymorphic Little Devil poison frog (Oophaga sylvatica) in Ecuador

Some South American poison frogs (Dendrobatidae) are chemically defended and use bright aposematic colors to warn potential predators of their unpalatability. Aposematic signals are often frequency‐dependent where individuals deviating from a local model are at a higher risk of predation. However, extreme diversity in the aposematic signal has been documented in poison frogs, especially in Oophaga. Here, we explore the phylogeographic pattern among color‐divergent populations of the Little Devil poison frog Oophaga sylvatica by analyzing population structure and genetic differentiation to evaluate which processes could account for color diversity within and among populations. With a combination of PCR amplicons (three mitochondrial and three nuclear markers) and genome‐wide markers from a double‐digested RAD (ddRAD) approach, we characterized the phylogenetic and genetic structure of 199 individuals from 13 populations (12 monomorphic and 1 polymorphic) across the O. sylvatica distribution. Individuals segregated into two main lineages by their northern or southern latitudinal distribution. A high level of genetic and phenotypic polymorphism within the northern lineage suggests ongoing gene flow. In contrast, low levels of genetic differentiation were detected among the southern lineage populations and support recent range expansions from populations in the northern lineage. We propose that a combination of climatic gradients and structured landscapes might be promoting gene flow and phylogenetic diversification. Alternatively, we cannot rule out that the observed phenotypic and genomic variations are the result of genetic drift on near or neutral alleles in a small number of genes.     

Genetic signatures of natural selection in response to air pollution in red spruce (Picea rubens,Pinaceae)

One of the most important drivers of local adaptation for forest trees is climate. Coupled to these patterns, however, are human‐induced disturbances through habitat modification and pollution. The confounded effects of climate and disturbance have rarely been investigated with regard to selective pressure on forest trees. Here, we have developed and used a population genetic approach to search for signals of selection within a set of 36 candidate genes chosen for their putative effects on adaptation to climate and human‐induced air pollution within five populations of red spruce (Picea rubens Sarg.), distributed across its natural range and air pollution gradient in eastern North America. Specifically, we used F_ST outlier and environmental correlation analyses to highlight a set of seven single nucleotide polymorphisms (SNPs) that were overly correlated with climate and levels of sulphate pollution after correcting for the confounding effects of population history. Use of three age cohorts within each population allowed the effects of climate and pollution to be separated temporally, as climate‐related SNPs (n = 7) showed the strongest signals in the oldest cohort, while pollution‐related SNPs (n = 3) showed the strongest signals in the youngest cohorts. These results highlight the usefulness of population genetic scans for the identification of putatively nonneutral evolution within genomes of nonmodel forest tree species, but also highlight the need for the development and application of robust methodologies to deal with the inherent multivariate nature of the genetic and ecological data used in these types of analyses.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.     

Low genetic diversity and strong but shallow population differentiation suggests genetic homogenization by metapopulation dynamics in a social spider

Mating systems and population dynamics influence genetic diversity and structure. Species that experience inbreeding and limited gene flow are expected to evolve isolated, divergent genetic lineages. Metapopulation dynamics with frequent extinctions and colonizations may, on the other hand, deplete and homogenize genetic variation, if extinction rate is sufficiently high compared to the effect of drift in local demes. We investigated these theoretical predictions empirically in social spiders that are highly inbred. Social spiders show intranest mating, female‐biased sex ratio, and frequent extinction and colonization events, factors that deplete genetic diversity within nests and populations and limit gene flow. We characterized population genetic structure in Stegodyphus sarasinorum, a social spider distributed across the Indian subcontinent. Species‐wide genetic diversity was estimated over approximately 2800 km from Sri Lanka to Himalayas, by sequencing 16 protein‐coding nuclear loci. We found 13 SNPs in 6592 bp (π = 0.00045) indicating low species‐wide nucleotide diversity. Three genetic lineages were strongly differentiated; however, only one fixed difference among them suggests recent divergence. This is consistent with a scenario of metapopulation dynamics that homogenizes genetic diversity across the species' range. Ultimately, low standing genetic variation may hamper a species' ability to track environmental change and render social inbreeding spiders ‘evolutionary dead‐ends’.     

netview p: a network visualization tool to unravel complex population structure using genome‐wide SNPs

Network‐based approaches are emerging as valuable tools for the analysis of complex genetic structure in wild and captive populations. netview p combines data quality control with the construction of population networks through mutual k‐nearest neighbours thresholds applied to genome‐wide SNPs. The program is cross‐platform compatible, open‐source and efficiently operates on data ranging from hundreds to hundreds of thousands of SNPs. The pipeline was used for the analysis of pedigree data from simulated (n = 750, SNPs = 1279) and captive silver‐lipped pearl oysters (n = 415, SNPs = 1107), wild populations of the European hake from the Atlantic and Mediterranean (n = 834, SNPs = 380) and grey wolves from North America (n = 239, SNPs = 78 255). The population networks effectively visualize large‐ and fine‐scale genetic structure within and between populations, including family‐level structure and relationships. netview p comprises a network‐based addition to other population analysis tools and provides user‐friendly access to a complex network analysis pipeline through implementation in python .     

Population structure of a vector of human diseases: Aedes aegypti in its ancestral range,Africa

Aedes aegypti, the major vector of dengue, yellow fever, chikungunya, and Zika viruses, remains of great medical and public health concern. There is little doubt that the ancestral home of the species is Africa. This mosquito invaded the New World 400‐500 years ago and later, Asia. However, little is known about the genetic structure and history of Ae. aegypti across Africa, as well as the possible origin(s) of the New World invasion. Here, we use ~17,000 genome‐wide single nucleotide polymorphisms (SNPs) to characterize a heretofore undocumented complex picture of this mosquito across its ancestral range in Africa. We find signatures of human‐assisted migrations, connectivity across long distances in sylvan populations, and of local admixture between domestic and sylvan populations. Finally, through a phylogenetic analysis combined with the genetic structure analyses, we suggest West Africa and especially Angola as the source of the New World's invasion, a scenario that fits well with the historic record of 16th‐century slave trade between Africa and Americas.     

Using RAD‐seq to develop sex‐linked markers in a haplodiplontic alga

For many taxa, including isomorphic haplodiplontic macroalgae, determining sex and ploidy is challenging, thereby limiting the scope of some population demographic and genetic studies. Here, we used double‐digest restriction site‐associated DNA sequencing (ddRAD‐seq) to identify sex‐linked molecular markers in the widespread red alga Agarophyton vermiculophyllum. In the ddRAD‐seq library, we included 10 female gametophytes, 10 male gametophytes, and 16 tetrasporophytes from one native and one non‐native site (N = 40 gametophytes and N = 32 tetrasporophytes total). We identified seven putatively female‐linked and 19 putatively male‐linked sequences. Four female‐ and eight male‐linked markers amplified in all three life cycle stages. Using one female‐ and one male‐linked marker that were sex‐specific, we developed a duplex PCR and tested the efficacy of this assay on a subset of thalli sampled at two sites in the non‐native range. We confirmed ploidy based on the visual observation of reproductive structures and previous microsatellite genotyping at 10 polymorphic loci. For 32 vegetative thalli, we were able to assign sex and confirm ploidy in these previously genotyped thalli. These markers will be integral to ongoing studies of A. vermiculophyllum invasion. We discuss the utility of RAD‐seq over other approaches previously used, such as RAPDs (random amplified polymorphic DNA), for future work designing sex‐linked markers in other haplodiplontic macroalgae for which genomes are lacking.     

Phylogeography and population genetics of pine butterflies: Sky islands increase genetic divergence

The sky islands of southeastern Arizona (AZ) mark a major transition zone between tropical and temperate biota and are considered a neglected biodiversity hotspot. Dispersal ability and host plant specificity are thought to impact the history and diversity of insect populations across the sky islands. We aimed to investigate the population structure and phylogeography of two pine‐feeding pierid butterflies, the pine white (Neophasia menapia) and the Mexican pine white (Neophasia terlooii), restricted to these “islands” at this transition zone. Given their dependence on pines as the larval hosts, we hypothesized that habitat connectivity affects population structure and is at least in part responsible for their allopatry. We sampled DNA from freshly collected butterflies from 17 sites in the sky islands and adjacent high‐elevation habitats and sequenced these samples using ddRADSeq. Up to 15,399 SNPs were discovered and analyzed in population genetic and phylogenetic contexts with Stacks and pyRAD pipelines. Low genetic differentiation in N. menapia suggests that it is panmictic. Conversely, there is strong evidence for population structure within N. terlooii. Each sky island likely contains a population of N. terlooii, and clustering is hierarchical, with populations on proximal mountains being more related to each other. The N. menapia habitat, which is largely contiguous, facilitates panmixia, while the N. terlooii habitat, restricted to the higher elevations on each sky island, creates distinct population structure. Phylogenetic results corroborate those from population genetic analyses. The historical climate‐driven fluxes in forest habitat connectivity have implications for understanding the biodiversity of fragmented habitats.     

An ecological genetic delineation of local seed‐source provenance for ecological restoration

An increasingly important practical application of the analysis of spatial genetic structure within plant species is to help define the extent of local provenance seed collection zones that minimize negative impacts in ecological restoration programs. Here, we derive seed sourcing guidelines from a novel range‐wide assessment of spatial genetic structure of 24 populations of Banksia menziesii (Proteaceae), a widely distributed Western Australian tree of significance in local ecological restoration programs. An analysis of molecular variance (AMOVA) of 100 amplified fragment length polymorphism (AFLP) markers revealed significant genetic differentiation among populations (Φ_PT = 0.18). Pairwise population genetic dissimilarity was correlated with geographic distance, but not environmental distance derived from 15 climate variables, suggesting overall neutrality of these markers with regard to these climate variables. Nevertheless, Bayesian outlier analysis identified four markers potentially under selection, although these were not correlated with the climate variables. We calculated a global R‐statistic using analysis of similarities (ANOSIM) to test the statistical significance of population differentiation and to infer a threshold seed collection zone distance of ~60 km (all markers) and 100 km (outlier markers) when genetic distance was regressed against geographic distance. Population pairs separated by >60 km were, on average, twice as likely to be significantly genetically differentiated than population pairs separated by <60 km, suggesting that habitat‐matched sites within a 30‐km radius around a restoration site genetically defines a local provenance seed collection zone for B. menziesii. Our approach is a novel probability‐based practical solution for the delineation of a local seed collection zone to minimize negative genetic impacts in ecological restoration.     

Extensive sampling of polar bears (Ursus maritimus) in the Northwest Passage (Canadian Arctic Archipelago) reveals population differentiation across multiple spatial and temporal scales

As global warming accelerates the melting of Arctic sea ice, polar bears (Ursus maritimus) must adapt to a rapidly changing landscape. This process will necessarily alter the species distribution together with population dynamics and structure. Detailed knowledge of these changes is crucial to delineating conservation priorities. Here, we sampled 361 polar bears from across the center of the Canadian Arctic Archipelago spanning the Gulf of Boothia (GB) and M'Clintock Channel (MC). We use DNA microsatellites and mitochondrial control region sequences to quantify genetic differentiation, estimate gene flow, and infer population history. Two populations, roughly coincident with GB and MC, are significantly differentiated at both nuclear (F_ST = 0.01) and mitochondrial (Φ_ST = 0.47; F_ST = 0.29) loci, allowing Bayesian clustering analyses to assign individuals to either group. Our data imply that the causes of the mitochondrial and nuclear genetic patterns differ. Analysis of mtDNA reveals the matrilineal structure dates at least to the Holocene, and is common to individuals throughout the species’ range. These mtDNA differences probably reflect both genetic drift and historical colonization dynamics. In contrast, the differentiation inferred from microsatellites is only on the scale of hundreds of years, possibly reflecting contemporary impediments to gene flow. Taken together, our data suggest that gene flow is insufficient to homogenize the GB and MC populations and support the designation of GB and MC as separate polar bear conservation units. Our study also provide a striking example of how nuclear DNA and mtDNA capture different aspects of a species demographic history.     

Distance,elevation and environment as drivers of diversity and divergence in bumble bees across latitude and altitude

Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site‐associated DNA sequencing (RADseq) in two bumble bee species, Bombus vosnesenskii and Bombus bifarius, across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A. Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure while B. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, with B. vosnesenskii exhibiting relatively consistent levels of genetic diversity across its range, while B. bifarius has reduced genetic diversity at low latitudes, where it is restricted to high‐elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.     

From microsatellites to single nucleotide polymorphisms for the genetic monitoring of a critically endangered sturgeon

The use of genetic information is crucial in conservation programs for the establishment of breeding plans and for the evaluation of restocking success. Short tandem repeats (STRs) have been the most widely used molecular markers in such programs, but next‐generation sequencing approaches have prompted the transition to genome‐wide markers such as single nucleotide polymorphisms (SNPs). Until now, most sturgeon species have been monitored using STRs. The low diversity found in the critically endangered European sturgeon (Acipenser sturio), however, makes its future genetic monitoring challenging, and the current resolution needs to be increased. Here, we describe the discovery of a highly informative set of 79 SNPs using double‐digest restriction‐associated DNA (ddRAD) sequencing and its validation by genotyping using the MassARRAY system. Comparing with STRs, the SNP panel proved to be highly efficient and reproducible, allowing for more accurate parentage and kinship assignments' on 192 juveniles of known pedigree and 40 wild‐born adults. We explore the effectiveness of both markers to estimated relatedness and inbreeding, using simulated and empirical datasets. Interestingly, we found significant correlations between STRs and SNPs at individual heterozygosity and inbreeding that give support to a reasonable representation of whole genome diversity for both markers. These results are useful for the conservation program of A. sturio in building a comprehensive studbook, which will optimize conservation strategies. This approach also proves suitable for other case studies in which highly discriminatory genetic markers are needed to assess parentage and kinship.