Population Genomic Variation Research Articles

Extensive Admixture Among Karst-Obligate Salamanders Reveals Evidence of Recent Divergence and Gene Exchange Through Aquifers.

Karst ecosystems often contain extraordinary biodiversity, but the complex underground aquifers of karst regions present challenges for assessing and conserving stygobiont diversity and investigating their evolutionary history. We examined the karst-obligate salamanders of the Eurycea neotenes species complex in the Edwards Plateau region of central Texas using population genomics data to address questions about population connectivity and the potential for gene exchange within the underlying aquifer system. The E. neotenes species complex has historically been divided into three nominal species, but their status, and spatial extent of species ranges, have remained uncertain. We discovered evidence of extensive admixture among species within the complex and with adjacent lineages. We observed relatively low levels of differentiation among all sampling localities which supports the hypothesis of recent divergence. Nominal taxonomy, aquifer region and geography each accounted for a modest amount of the overall population genomic variation; however, these predictors were largely confounded and difficult to disentangle. Importantly, current taxonomy of the three nominal species does not reflect the admixture apparent in clustering analyses. Inference of migration events revealed a complex pattern of gene exchange, suggesting that Eurycea salamanders have a dynamic history of dispersal through the aquifer system. These results highlight the need for greater understanding of how stygobiont populations are connected via dispersal and gene exchange through karst aquifers. These results also highlight the applicability of population genomics data as a powerful lever for investigating connectivity among populations in systems where direct detection of dispersal paths is difficult, as in underground, aquatic systems.

Every Cog and Wheel: Identifying Biocomplexity at the Genomic and Phenotypic Level in a Population Complex of Chinook Salmon

Genetic diversity is the fundamental building block of biodiversity and the necessary ingredient for adaptation. Specifically, the intra-specific diversity (biocomplexity) comprised of phenotypic and genetic variation partitioned within and among populations can determine the ability of a species to respond to changing environmental conditions. Here, we explore the biocomplexity of California’s Central Valley Chinook salmon (Oncorhynchus tshawytscha) population complex at the genomic level by quantifying population genomic diversity among and within migration life-history phenotypes. Notably, despite apparent gene flow among populations with the same migration (life history) phenotypes inhabiting different tributaries, each group is characterized by a distinct component of unique genomic diversity. While enumerating biodiversity contained within individual hierarchical levels is informative, it is important to consider inter- and intra-specific diversity simultaneously as there may be emergent properties at higher levels due to presence of diversity at lower ones. Our results emphasize the importance of formulating conservation goals focused to maintain biocomplexity at both the phenotypic and genotypic level. Doing so will preserve the species’ adaptive potential and increase the probability of persistence of the population complex despite changing environmental pressures.

Comparative Population Genomic Diversity and Differentiation in Trapdoor Spiders and Relatives (Araneae, Mygalomorphae).

Although patterns of population genomic variation are well-studied in animals, there remains room for studies that focus on non-model taxa with unique biologies. Here we characterise and attempt to explain such patterns in mygalomorph spiders, which are generally sedentary, often occur as spatially clustered demes and show remarkable longevity. Genome-wide single nucleotide polymorphism (SNP) data were collected for 500 individuals across a phylogenetically representative sample of taxa. We inferred genetic populations within focal taxa using a phylogenetically informed clustering approach, and characterised patterns of diversity and differentiation within- and among these genetic populations, respectively. Using phylogenetic comparative methods we asked whether geographical range sizes and ecomorphological variables (behavioural niche and body size) significantly explain patterns of diversity and differentiation. Specifically, we predicted higher genetic diversity in genetic populations with larger geographical ranges, and in small-bodied taxa. We also predicted greater genetic differentiation in small-bodied taxa, and in burrowing taxa. We recovered several significant predictors of genetic diversity, but not genetic differentiation. However, we found generally high differentiation across genetic populations for all focal taxa, and a consistent signal for isolation-by-distance irrespective of behavioural niche or body size. We hypothesise that high population genetic structuring, likely reflecting combined dispersal limitation and microhabitat specificity, is a shared trait for all mygalomorphs. Few studies have found ubiquitous genetic structuring for an entire ancient and species-rich animal clade.

Fitness consequences of population bottlenecks in an invasive blowfly.

Invasive species often undergo demographic bottlenecks that cause a decrease in genetic diversity and associated reductions in population fitness. Despite this, they manage to thrive in novel environments. Investigating the effects of inbreeding and genetic bottlenecks on population fitness for invasive species is, therefore, key to understanding how they may survive in new environments. We used the blowfly Calliphora vicina (Sciences, Mathématiques et Physique, 1830, 2, 1), which is native to Europe and was introduced to Australia and New Zealand, to examine the effects of genetic diversity on population fitness. We first collected 59 samples from 15 populations across New Zealand and one in Australia, and used 20,501 biallelic SNPs to investigate population genomic diversity, structure and admixture. We then explored the impacts of repeated experimental bottlenecks on population fitness by creating inbred and outbred lines of C. vicina and measuring a variety of fitness traits. In wild-caught samples, we found low overall genetic diversity, signals of genetic admixture and limited (<3%) genetic differentiation between North and South Island populations, with genetic links between the South Island and Australia. Following experimental bottlenecks, we found significant reductions in fitness for inbred lines. However, fitness effects were not felt equally across all phenotypic traits. Moreover, they were not enough to cause population collapse in any experimental line, suggesting that C. vicina (when under relaxed selection, as in laboratory settings) may be able to compensate for population bottlenecks even when highly inbred. Our results demonstrate the value of a tractable experimental system for investigating processes that may facilitate or hamper biological invasion.

The 'queen of the Andes' (Puya raimondii) is genetically fragile and fragmented: a consequence of long generation time and semelparity?

Understanding how life history shapes genetic diversity is a fundamental issue in evolutionary biology, with important consequences for conservation. However, we still have an incomplete picture of the impact of life history on genome-wide patterns of diversity, especially in long-lived semelparous plants. Puya raimondii is a high-altitude semelparous species from the Andes that flowers at 40-100 years of age. We sequenced the whole genome and estimated the nucleotide diversity of 200 individuals sampled from nine populations. Coalescent-based approaches were then used to infer past population dynamics. Finally, these results were compared with results obtained for the iteroparous species, Puya macrura. The nine populations of P. raimondii were highly divergent, highly inbred, and carried an exceptionally high genetic load. They are genetically depauperate, although, locally in the genome, balancing selection contributed to the maintenance of genetic polymorphism. While both P. raimondii and P. macrura went through a severe bottleneck during the Pleistocene, P. raimondii did not recover from it and continuously declined, while P. macrura managed to bounce back. Our results demonstrate the importance of life history, in particular generation time and reproductive strategy, in affecting population dynamics and genomic variation, and illustrate the genetic fragility of long-lived semelparous plants.

Patterns of genomic variation reveal a single evolutionary origin of the wild allotetraploid Mimulus sookensis.

Polyploidy occurs across the tree of life and is especially common in plants. Because newly formed cytotypes are often incompatible with their progenitors, polyploidy is also said to trigger "instantaneous" speciation. If a polyploid can self-fertilize or reproduce asexually, it is even possible for one individual to produce an entirely new lineage, but how often this scenario occurs is unclear. Here, we investigate the evolutionary history of the wild allotetraploid Mimulus sookensis, which was formed through hybridization between self-compatible, diploid species in the Mimulus guttatus complex. We generate a chromosome-scale reference assembly for M. sookensis and define its distinct subgenomes. Despite previous reports suggesting multiple origins of this highly selfing polyploid, we discover patterns of population genomic variation that provide unambiguous support for a single origin. One M. sookensis subgenome is clearly derived from the selfer Mimulus nasutus, which organellar variation suggests is the maternal progenitor. The ancestor of the other subgenome is less certain, but it shares variation with both Mimulus decorus and M. guttatus, two outcrossing diploids with geographic ranges that overlap broadly with M. sookensis. This study establishes M. sookensis as an example of instantaneous speciation, likely facilitated by the polyploid's predisposition to self-fertilize.

Shared patterns of population genomic variation and phenotypic response across rapid range expansions in two invasive lady beetle species

Non-native lady beetle species have often been introduced, with variable success, into North America for biological control of aphids, scales, whiteflies, and other agricultural pests. Two predatory lady beetle species, Propylea quatuordecimpunctata and Hippodamia variegata, both originating from Eurasia, were first discovered near Montreal, Quebec, in North America in 1968 and 1984, respectively, and have since expanded into northeastern North America and the midwestern United States. In this study, we estimate the range-wide population structure, establishment and range-expansion, and recent evolutionary history of these lady beetle species using reduced-representation genotyping-by-sequencing via ddRADseq. In addition, we quantified the responses to a key abiotic factor, photoperiod, that regulates adult reproductive diapause in these two species and may influence their geographical range in North America. Our analyses detect: (1) non-significant genetic differentiation and divergence among North American populations that likely originated from a single accidental introduction, (2) evidence of reduced contemporary gene flow within the continental US, and (3) minor phenotypic differences in diapause induction between populations of the two species from the same location.

SNP-SVant: A Computational Workflow to Predict and Annotate Genomic Variants in Organisms Lacking Benchmarked Variants.

Whole-genome sequencing is widely used to investigate population genomic variation in organisms of interest. Assorted tools have been independently developed to call variants from short-read sequencing data aligned to a reference genome, including single nucleotide polymorphisms (SNPs) and structural variations (SVs). We developed SNP-SVant, an integrated, flexible, and computationally efficient bioinformatic workflow that predicts high-confidence SNPs and SVs in organisms without benchmarked variants, which are traditionally used for distinguishing sequencing errors from real variants. In the absence of these benchmarked datasets, we leverage multiple rounds of statistical recalibration to increase the precision of variant prediction. The SNP-SVant workflow is flexible, with user options to tradeoff accuracy for sensitivity. The workflow predicts SNPs and small insertions and deletions using the Genome Analysis ToolKit (GATK) and predicts SVs using the Genome Rearrangement IDentification Software Suite (GRIDSS), and it culminates in variant annotation using custom scripts. A key utility of SNP-SVant is its scalability. Variant calling is a computationally expensive procedure, and thus, SNP-SVant uses a workflow management system with intermediary checkpoint steps to ensure efficient use of resources by minimizing redundant computations and omitting steps where dependent files are available. SNP-SVant also provides metrics to assess the quality of called variants and converts between VCF and aligned FASTA format outputs to ensure compatibility with downstream tools to calculate selection statistics, which are commonplace in population genomics studies. By accounting for both small and large structural variants, users of this workflow can obtain a wide-ranging view of genomic alterations in an organism of interest. Overall, this workflow advances our capabilities in assessing the functional consequences of different types of genomic alterations, ultimately improving our ability to associate genotypes with phenotypes. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Predicting single nucleotide polymorphisms and structural variations Support Protocol 1: Downloading publicly available sequencing data Support Protocol 2: Visualizing variant loci using Integrated Genome Viewer Support Protocol 3: Converting between VCF and aligned FASTA formats.

Unveiling recent and ongoing adaptive selection in human populations.

Genome-wide scans for signals of selection have become a routine part of the analysis of population genomic variation datasets and have resulted in compelling evidence of selection during recent human evolution. This Essay spotlights methodological innovations that have enabled the detection of selection over very recent timescales, even in contemporary human populations. By harnessing large-scale genomic and phenotypic datasets, these new methods use different strategies to uncover connections between genotype, phenotype, and fitness. This Essay outlines the rationale and key findings of each strategy, discusses challenges in interpretation, and describes opportunities to improve detection and understanding of ongoing selection in human populations.

Disentangling the contemporary and historical effects of landscape on the population genomic variation of two bird species restricted to the highland forest enclaves of northeastern Brazil.

Investigating the impact of landscape features on patterns of genetic variation is crucial to understand spatially dependent evolutionary processes. Here, we assess the population genomic variation of two bird species (Conopophaga cearae and Sclerurus cearensis) through the Caatinga moist forest enclaves in northeastern Brazil. To infer the evolutionary dynamics of bird populations through the Late Quaternary, we used genome-wide polymorphism data obtained from double-digestion restriction-site-associated DNA sequencing (ddRADseq), and integrated population structure analyses, historical demography models, paleodistribution modeling, and landscape genetics analyses. We found the population differentiation among enclaves to be significantly related to the geographic distance and historical resistance across the rugged landscape. The climate changes at the end of the Pleistocene to the Holocene likely triggered synchronic population decline in all enclaves for both species. Our findings revealed that both geographic distance and historical connectivity through highlands are important factors that can explain the current patterns of genetic variation. Our results further suggest that levels of population differentiation and connectivity cannot be explained purely on the basis of contemporary environmental conditions. By combining historical demographic analyses and niche modeling predictions in a historical framework, we provide strong evidence that climate fluctuations of the Quaternary promoted population differentiation and a high degree of temporal synchrony among population size changes in both species.

Population genomic diversity and structure in the golden bandicoot: a history of isolation, extirpation, and conservation

Using genetic information to develop and implement conservation programs is vital for maintaining biodiversity and ecosystem resilience. Evaluation of the genetic variability within and among remnant populations can inform management of both natural and translocated populations to maximise species’ adaptive potential, mitigate negative impacts of inbreeding, and subsequently minimise risk of extinction. Here we use reduced representation sequencing to undertake a genetic assessment of the golden bandicoot (Isoodon auratus), a threatened marsupial endemic to Australia. The currently recognised taxon consists of three subspecies distributed among multiple natural and translocated populations. After confirming the genetic distinctiveness of I. auratus from two closely related taxa, I. fusciventer and I. macrourus, we identified four genetic clusters within I. auratus. These clusters exhibited substantial genetic differentiation (pairwise FST values ranging from 0.18 to 0.65, pairwise DXY ranging from 0.1 to 0.168), reflecting long-term isolation of some populations on offshore islands and the influence of genetic drift. Mainland natural populations in the Kimberley region had the highest genetic diversity and the largest contribution to overall allelic and gene diversity compared to both natural and translocated island populations. A population translocated to Guluwuru Island in the Northern Territory had the lowest genetic diversity. Our data suggest that island populations can appear genetically unique due to genetic drift and this needs to be taken into account when considering genetic diversity in conservation efforts to maintain overall genetic diversity of the species. We effectively demonstrate how genomic information can guide practical conservation planning, especially when declining species are represented by multiple isolated populations.

Chronosequence of invasion reveals minimal losses of population genomic diversity, niche expansion, and trait divergence in the polyploid, leafy spurge.

Rapid evolution may play an important role in the range expansion of invasive species and modify forecasts of invasion, which are the backbone of land management strategies. However, losses of genetic variation associated with colonization bottlenecks may constrain trait and niche divergence at leading range edges, thereby impacting management decisions that anticipate future range expansion. The spatial and temporal scales over which adaptation contributes to invasion dynamics remain unresolved. We leveraged detailed records of the ~130-year invasion history of the invasive polyploid plant, leafy spurge (Euphorbia virgata), across ~500 km in Minnesota, U.S.A. We examined the consequences of range expansion for population genomic diversity, niche breadth, and the evolution of germination behavior. Using genotyping-by-sequencing, we found some population structure in the range core, where introduction occurred, but panmixia among all other populations. Range expansion was accompanied by only modest losses in sequence diversity, with small, isolated populations at the leading edge harboring similar levels of diversity to those in the range core. The climatic niche expanded during most of the range expansion, and the niche of the range core was largely non-overlapping with the invasion front. Ecological niche models indicated that mean temperature of the warmest quarter was the strongest determinant of habitat suitability and that populations at the leading edge had the lowest habitat suitability. Guided by these findings, we tested for rapid evolution in germination behavior over the time course of range expansion using a common garden experiment and temperature manipulations. Germination behavior diverged from the early to late phases of the invasion, with populations from later phases having higher dormancy at lower temperatures. Our results suggest that trait evolution may have contributed to niche expansion during invasion and that distribution models, which inform future management planning, may underestimate invasion potential without accounting for evolution.

Evidence of introgression, ecological divergence and adaptation in Asterias sea stars.

Hybrid zones are important windows into the evolutionary dynamics of populations, revealing how processes like introgression and adaptation structure population genomic variation. Importantly, they are useful for understanding speciation and how species respond to their environments. Here, we investigate two closely related sea star species, Asterias rubens and A. forbesi, distributed along rocky European and North American coastlines of the North Atlantic, and use genome-wide molecular markers to infer the distribution of genomic variation within and between species in this group. Using genomic data and environmental niche modelling, we document hybridization occurring between northern New England and the southern Canadian Maritimes. We investigate the factors that maintain this hybrid zone, as well as the environmental variables that putatively drive selection within and between species. We find that the two species differ in their environmental niche breadth; Asterias forbesi displays a relatively narrow environmental niche while conversely, A. rubens has a wider niche breadth. Species distribution models accurately predict hybrids to occur within environmental niche overlap, thereby suggesting environmental selection plays an important role in the maintenance of the hybrid zone. Our results imply that the distribution of genomic variation in North Atlantic sea stars is influenced by the environment, which will be crucial to consider as the climate changes.

Community-wide correlations between species richness, abundance and population genomic diversity in a freshwater biodiversity hotspot.

Understanding patterns of diversity across macro (e.g. species-level) and micro (e.g. molecular-level) scales can shed light on community function and stability by elucidating the abiotic and biotic drivers of diversity within ecological communities. We examined the relationships among taxonomic and genetic metrics of diversity in freshwater mussels (Bivalvia: Unionidae), an ecologically important and species-rich group in the southeastern United States. Using quantitative community surveys and reduced-representation genome sequencing across 22 sites in seven rivers and two river basins, we surveyed 68 mussel species and sequenced 23 of these species to characterize intrapopulation genetic variation. We tested for the presence of species diversity-abundance correlations (i.e. the more-individuals hypothesis, MIH), species-genetic diversity correlations (SGDCs) and abundance-genetic diversity correlations (AGDCs) across all sites to evaluate relationships between different metrics of diversity. Sites with greater cumulative multispecies density (a standardized metric of abundance) had a greater number of species, consistent with the MIH hypothesis. Intrapopulation genetic diversity was strongly associated with the density of most species, indicating the presence of AGDCs. However, there was no consistent evidence for SGDCs. Although sites with greater overall densities of mussels had greater species richness, sites with higher genetic diversity did not always exhibit positive correlations with species richness, suggesting that there are spatial and evolutionary scales at which the processes influencing community-level diversity and intraspecific diversity differ. Our work reveals the importance of local abundance as indicator (and possibly a driver) of intrapopulation genetic diversity.

Growth, drought response, and climate-associated genomic structure in whitebark pine in the Sierra Nevada of California.

Whitebark pine (Pinus albicaulis Engelm.) has experienced rapid population declines and is listed as threatened under the Endangered Species Act in the United States. Whitebark pine in the Sierra Nevada of California represents the southernmost end of the species' distribution and, like other portions of its range, faces threats from an introduced pathogen, native bark beetles, and a rapidly warming climate. Beyond these chronic stressors, there is also concern about how this species will respond to acute stressors, such as drought. We present patterns of stem growth from 766 large (average diameter at breast height >25 cm), disease-free whitebark pine across the Sierra Nevada before and during a recent period of drought. We contextualize growth patterns using population genomic diversity and structure from a subset of 327 trees. Sampled whitebark pine generally had positive to neutral stem growth trends from 1970 to 2011, which was positively correlated with minimum temperature and precipitation. Indices of stem growth during drought years (2012 to 2015) relative to a predrought interval were mostly positive to neutral at our sampled sites. Individual tree growth response phenotypes appeared to be linked to genotypic variation in climate-associated loci, suggesting that some genotypes can take better advantage of local climatic conditions than others. We speculate that reduced snowpack during the 2012 to 2015 drought years may have lengthened the growing season while retaining sufficient moisture to maintain growth at most study sites. Growth responses may differ under future warming, however, particularly if drought severity increases and modifies interactions with pests and pathogens.

Population-based genomic study of Plasmodium vivax malaria in seven Brazilian states and across South America

Brazil is a unique and understudied setting for malaria, with complex foci of transmission associated with human and environmental conditions. An understanding of the population genomic diversity of P.vivax parasites across Brazil can support malaria control strategies. Through whole genome sequencing of P.vivax isolates across 7 Brazilian states, we use population genomic approaches to compare genetic diversity within country (n=123), continent (6 countries, n=315) and globally (26countries, n=885). We confirm that South American isolates are distinct, have more ancestral populations than the other global regions, with differentiating mutations in genes under selective pressure linked to antimalarial drugs (pvmdr1, pvdhfr-ts) and mosquito vectors (pvcrmp3, pvP45/48, pvP47). We demonstrate Brazil as a distinct parasite population, with signals of selection including ABC transporter (PvABCI3) and PHIST exported proteins. Brazil has a complex population structure, with evidence of P.simium infections and Amazonian parasites separating into multiple clusters. Overall, our work provides the first Brazil-wide analysis of P.vivax population structure and identifies important mutations, which can inform future research and control measures. AI is funded by an MRC LiD PhD studentship. TGC is funded by the Medical Research Council (Grant no. MR/M01360X/1, MR/N010469/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1). SC is funded by Medical Research Council UK grants (MR/M01360X/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1) and Bloomsbury SET (ref. CCF17-7779). FN is funded by The Shloklo Malaria Research Unit - part of the Mahidol Oxford Research Unit, supported by the Wellcome Trust (Grant no. 220211). ARSB is funded by São Paulo Research Foundation - FAPESP (Grant no. 2002/09546-1). RLDM is funded by Brazilian National Council for Scientific and Technological Development - CNPq (Grant no. 302353/2003-8 and 471605/2011-5); CRFM is funded by FAPESP (Grant no. 2020/06747-4) and CNPq (Grant no. 302917/2019-5 and 408636/2018-1); JGD is funded by FAPESP fellowships (2016/13465-0 and 2019/12068-5) and CNPq (Grant no. 409216/2018-6).

Whole-genome resequencing reveals signature of local adaptation and divergence in wild soybean.

Global climate change has threatened world crop production and food security. Decoding the adaptive genetic basis of wild relatives provides an invaluable genomic resource for climate-smart crop breedinG.Here, we performed whole-genome sequencing of 185 diverse wild soybean (Glycine soja) accessions collected from three major agro-ecological zones in China to parse the genomic basis of local adaptation in wild soybean. The population genomic diversity pattern exhibited clear agro-ecological zone-based population structure, and multiple environmental factors were observed to contribute to the genetic divergence. Demographic analysis shows that wild soybeans from the three ecological zones diverged about 1 × 105years ago, and then the effective population sizes have undergone different degrees of expansions. Genome-environment association identified multiple genes involved in the local adaptation, such as flowering time and temperature-related genes. A locus containing two adjacent MADS-box transcription factors on chromosome 19 was identified for multiple environmental factors, and it experienced positive selection that enables the adaptation to high-latitude environment. This study provides insights into the genetic mechanism of ecological adaptation in wild soybean that may facilitate climate-resilient soybean breeding.

A range-wide analysis of population structure and genomic variation within the critically endangered spiny daisy (Acanthocladium dockeri)

Understanding population structure and genetic diversity is important for designing effective conservation strategies. As a critically endangered shrub, the six remaining extant populations of spiny daisy (Acanthocladium dockeri) are restricted to country roadsides in the mid-north of South Australia, where the species faces many ongoing abiotic and biotic threats to survival. Currently the spiny daisy is managed by selecting individuals from the extant populations and translocating them to establish insurance populations. However, there is little information available on the genetic differentiation between populations and diversity within source populations, which are essential components of planning translocations. To help fill this knowledge gap, we analysed population structure within and among all six of its known wild populations using 7,742 SNPs generated by a genotyping-by-sequencing approach. Results indicated that each population was strongly differentiated, had low levels of genetic diversity, and there was no evidence of inter-population gene flow. Individuals within each population were generally closely related, however, the Melrose population consisted entirely of clones. Our results suggest genetic rescue should be applied to wild spiny daisy populations to increase genetic diversity that will subsequently lead to greater intra-population fitness and adaptability. As a starting point, we suggest focussing on improving seed viability via inter-population crosses such as through hand pollination experiments to experimentally assess their sexual compatibility with the hope of increasing spiny daisy sexual reproduction and long-term reproductive fitness.

Museum Genomics of an Agricultural Super-Pest, the Colorado Potato Beetle, Leptinotarsa decemlineata (Chrysomelidae), Provides Evidence of Adaptation from Standing Variation.

Despite extensive research on agricultural pests, our knowledge about their evolutionary history is often limited. A mechanistic understanding of the demographic changes and modes of adaptation remains an important goal, as it improves our understanding of organismal responses to environmental change and our ability to sustainably manage pest populations. Emerging genomic datasets now allow for characterization of demographic and adaptive processes, but face limits when they are drawn from contemporary samples, especially in the context of strong demographic change, repeated selection, or adaptation involving modest shifts in allele frequency at many loci. Temporal sampling, however, can improve our ability to reconstruct evolutionary events. Here, we leverage museum samples to examine whether population genomic diversity and structure has changed over time, and to identify genomic regions that appear to be under selection. We focus on the Colorado potato beetle (CPB), Leptinotarsa decemlineata (Say 1824; Coleoptera: Chrysomelidae), which is widely regarded as a super-pest due to its rapid, and repeated, evolution to insecticides. By combining whole genome resequencing data from 78 museum samples with modern sampling, we demonstrate that CPB expanded rapidly in the 19th century, leading to a reduction in diversity and limited genetic structure from the Midwest to Northeast United States. Temporal genome scans provide extensive evidence for selection acting in resistant field populations in Wisconsin and New York, including numerous known insecticide resistance genes. We also validate these results by showing that known selective sweeps in modern populations are identified by our genome scan. Perhaps most importantly, temporal analysis indicates selection on standing genetic variation, as we find evidence for parallel evolution in the two geographical regions. Parallel evolution involves a range of phenotypic traits not previously identified as under selection in CPB, such as reproductive and morphological functional pathways that might be important for adaptation to agricultural habitats.

Islands in the desert: assessing fine scale population genomic variation of a group of imperiled desert fishes

Islands in the desert: assessing fine scale population genomic variation of a group of imperiled desert fishes