Footprints of local adaptation span hundreds of linked genes in the Atlantic silverside genome

Signatures of adaptation at key insecticide resistance loci in Anopheles gambiae in Southern Ghana revealed by reduced-coverage WGS

Resistance to insecticides and adaptation to a diverse range of environments present challenges to Anopheles gambiae s.l. mosquito control efforts in sub-Saharan Africa. Whole-genome-sequencing is often employed for identifying the genomic basis underlying adaptation in Anopheles, but remains expensive for large-scale surveys. Reduced coverage whole-genome-sequencing can identify regions of the genome involved in adaptation at a lower cost, but is currently untested in Anopheles mosquitoes. Here, we use reduced coverage WGS to investigate population genetic structure and identify signatures of local adaptation in Anopheles mosquitoes across southern Ghana. In contrast to previous analyses, we find no structuring by ecoregion, with Anopheles coluzzii and Anopheles gambiae populations largely displaying the hallmarks of large, unstructured populations. However, we find signatures of selection at insecticide resistance loci that appear ubiquitous across ecoregions in An. coluzzii, and strongest in forest ecoregions in An. gambiae. Our study highlights resistance candidate genes in this region, and validates reduced coverage WGS, potentially to very low coverage levels, for population genomics and exploratory surveys for adaptation in Anopheles taxa.

Loci associated with cave-derived traits concentrate in specific regions of the Mexican cavefish genome

A major goal of modern evolutionary biology is connecting phenotypic evolution with its underlying genetic basis. The Mexican cavefish (Astyanax mexicanus), a characin fish species comprised of a surface ecotype and a cave-derived ecotype, is well suited as a model to study the genetic mechanisms underlying adaptation to extreme environments. Here we map 206 previously published quantitative trait loci (QTL) for cave-derived traits in A. mexicanus to the newest version of the surface fish genome assembly, AstMex3. This analysis revealed that QTL cluster in the genome more than expected by chance, and this clustering is not explained by the distribution of genes in the genome. To investigate whether certain characteristics of the genome facilitate phenotypic evolution, we tested whether genomic characteristics, such as highly mutagenic CpG sites, are reliable predictors of the sites of trait evolution but did not find any significant trends. Finally, we combined the QTL map with previously collected expression and selection data to identify a list of 36 candidate genes that may underlie the repeated evolution of cave phenotypes, including rgrb which is predicted to be involved in phototransduction. We found this gene has disrupted exons in all non-hybrid cave populations but intact reading frames in surface fish. Overall, our results suggest specific "evolutionary hotspots" in the genome may play significant roles in driving adaptation to the cave environment in Astyanax mexicanus and demonstrate how this compiled dataset can facilitate our understanding of the genetic basis of repeated evolution in the Mexican cavefish.

Predicting the future of our oceans-Evaluating genomic forecasting approaches in marine species

Climate change is restructuring biodiversity on multiple scales and there is a pressing need to understand the downstream ecological and genomic consequences of this change. Recent advancements in the field of eco-evolutionary genomics have sought to include evolutionary processes in forecasting species' responses to climate change (e.g., genomic offset), but to date, much of this work has focused on terrestrial species. Coastal and offshore species, and the fisheries they support, may be even more vulnerable to climate change than their terrestrial counterparts, warranting a critical appraisal of these approaches in marine systems. First, we synthesize knowledge about the genomic basis of adaptation in marine species, and then we discuss the few examples where genomic forecasting has been applied in marine systems. Next, we identify the key challenges in validating genomic offset estimates in marine species, and we advocate for the inclusion of historical sampling data and hindcasting in the validation phase. Lastly, we describe a workflow to guide marine managers in incorporating these predictions into the decision-making process.

Fine-scale environmentally associated spatial structure of lumpfish (Cyclopterus lumpus) across the Northwest Atlantic

Lumpfish, Cyclopterus lumpus, have historically been harvested throughout Atlantic Canada and are increasingly in demand as a solution to controlling sea lice in Atlantic salmon farms-a process which involves both the domestication and the transfer of lumpfish between geographic regions. At present, little is known regarding population structure and diversity of wild lumpfish in Atlantic Canada, limiting attempts to assess the potential impacts of escaped lumpfish individuals from salmon pens on currently at-risk wild populations. Here, we characterize the spatial population structure and genomic-environmental associations of wild populations of lumpfish throughout the Northwest Atlantic using both 70K SNP array data and whole-genome re-sequencing data (WGS). At broad spatial scales, our results reveal a large environmentally associated genetic break between the southern populations (Gulf of Maine and Bay of Fundy) and northern populations (Newfoundland and the Gulf of St. Lawrence), linked to variation in ocean temperature and ice cover. At finer spatial scales, evidence of population structure was also evident in a distinct coastal group in Newfoundland and significant isolation by distance across the northern region. Both evidence of consistent environmental associations and elevated genome-wide variation in F ST values among these three regional groups supports their biological relevance. This study represents the first extensive description of population structure of lumpfish in Atlantic Canada, revealing evidence of broad and fine geographic scale environmentally associated genomic diversity. Our results will facilitate the commercial use of lumpfish as a cleaner fish in Atlantic salmon aquaculture, the identification of lumpfish escapees, and the delineation of conservation units of this at-risk species throughout Atlantic Canada.

Importance of timely metadata curation to the global surveillance of genetic diversity

Genetic diversity within species represents a fundamental yet underappreciated level of biodiversity. Because genetic diversity can indicate species resilience to changing climate, its measurement is relevant to many national and global conservation policy targets. Many studies produce large amounts of genome-scale genetic diversity data for wild populations, but most (87%) do not include the associated spatial and temporal metadata necessary for them to be reused in monitoring programs or for acknowledging the sovereignty of nations or Indigenous peoples. We undertook a distributed datathon to quantify the availability of these missing metadata and to test the hypothesis that their availability decays with time. We also worked to remediate missing metadata by extracting them from associated published papers, online repositories, and direct communication with authors. Starting with 848 candidate genomic data sets (reduced representation and whole genome) from the International Nucleotide Sequence Database Collaboration, we determined that 561 contained mostly samples from wild populations. We successfully restored spatiotemporal metadata for 78% of these 561 data sets (n = 440 data sets with data on 45,105 individuals from 762 species in 17 phyla). Examining papers and online repositories was much more fruitful than contacting 351 authors, who replied to our email requests 45% of the time. Overall, 23% of our email queries to authors unearthed useful metadata. The probability of retrieving spatiotemporal metadata declined significantly as age of the data set increased. There was a 13.5% yearly decrease in metadata associated with published papers or online repositories and up to a 22% yearly decrease in metadata that were only available from authors. This rapid decay in metadata availability, mirrored in studies of other types of biological data, should motivate swift updates to data-sharing policies and researcher practices to ensure that the valuable context provided by metadata is not lost to conservation science forever.

Detecting parallel polygenic adaptation to novel evolutionary pressure in wild populations: a case study in Atlantic cod (Gadus morhua)

Populations can adapt to novel selection pressures through dramatic frequency changes in a few genes of large effect or subtle shifts in many genes of small effect. The latter (polygenic adaptation) is expected to be the primary mode of evolution for many life-history traits but tends to be more difficult to detect than changes in genes of large effect. Atlantic cod (Gadus morhua) were subjected to intense fishing pressure over the twentieth century, leading to abundance crashes and a phenotypic shift toward earlier maturation across many populations. Here, we use spatially replicated temporal genomic data to test for a shared polygenic adaptive response to fishing using methods previously applied to evolve-and-resequence experiments. Cod populations on either side of the Atlantic show covariance in allele frequency change across the genome that are characteristic of recent polygenic adaptation. Using simulations, we demonstrate that the degree of covariance in allele frequency change observed in cod is unlikely to be explained by neutral processes or background selection. As human pressures on wild populations continue to increase, understanding and attributing modes of adaptation using methods similar to those demonstrated here will be important in identifying the capacity for adaptive responses and evolutionary rescue. This article is part of the theme issue 'Detecting and attributing the causes of biodiversity change: needs, gaps and solutions'.

Where to draw the line? Expanding the delineation of conservation units to highly mobile taxa

Conservation units (CUs) are an essential tool for maximizing evolutionary potential and prioritizing areas across a species' range for protection when implementing conservation and management measures. However, current workflows for identifying CUs on the basis of neutral and adaptive genomic variation largely ignore information contained in patterns of isolation by distance (IBD), frequently the primary signal of population structure in highly mobile taxa, such as birds, bats, and marine organisms with pelagic larval stages. While individuals located on either end of a species' distribution may exhibit clear genetic, phenotypic, and ecological differences, IBD produces subtle changes in allele frequencies across space, making it difficult to draw clear boundaries for conservation purposes in the absence of discrete population structure. Here, we highlight potential pitfalls that arise when applying common methods for delineating CUs to continuously distributed organisms and review existing methods for detecting subtle breakpoints in patterns of IBD that can indicate barriers to gene flow in highly mobile taxa. In addition, we propose a new framework for identifying CUs in all organisms, including those characterized by continuous genomic differentiation, and suggest several possible ways to harness the information contained in patterns of IBD to guide conservation and management decisions.

Latitudinal variation and plasticity in response to temperature in Geukensia demissa

As global temperatures warm, species must adapt to a changing climate or transition to a different location suitable for their survival. Understanding the extent to which species are able to do so, particularly keystone species, is imperative to ensuring the survival of key ecosystems. The ribbed mussel Geukensia demissa is an integral part of salt marshes along the Atlantic coast of North America. Spatial patterns of genomic and phenotypic divergence have been previously documented, although their link with coastal environmental variation is unknown. Here, we study how populations of G. demissa in the northern (Massachusetts) and southern (Georgia) portions of the species range respond to changes in temperature. We combine assays of variation in oxygen consumption and RNA transcriptomic data with genomic divergence analyses to identify how separate populations of G. demissa may vary in distinct thermal environments. Our results show differences in constitutive oxygen consumption between mussels from Georgia and Massachusetts, as well as shared and disparate patterns of gene expression across temperature profiles. We also find that metabolic genes seem to be a strong component of divergence between these two populations. Our analysis highlights the importance of studying integrative patterns of genomic and phenotypic variation in species that are key for particular ecosystems, and how they might respond to further changes in climate.

Karyotype differentiation in the Nothobranchiusugandensis species group (Teleostei, Cyprinodontiformes), seasonal fishes from the east African inland plateau, in the context of phylogeny and biogeography

The karyotype differentiation of the twelve known members of the Nothobranchiusugandensis Wildekamp, 1994 species group is reviewed and the karyotype composition of seven of its species is described herein for the first time using a conventional cytogenetic protocol. Changes in the architecture of eukaryotic genomes often have a major impact on processes underlying reproductive isolation, adaptation and diversification. African annual killifishes of the genus Nothobranchius Peters, 1868 (Teleostei: Nothobranchiidae), which are adapted to an extreme environment of ephemeral wetland pools in African savannahs, feature extensive karyotype evolution in small, isolated populations and thus are suitable models for studying the interplay between karyotype change and species evolution. The present investigation reveals a highly conserved diploid chromosome number (2n = 36) but a variable number of chromosomal arms (46-64) among members of the N.ugandensis species group, implying a significant role of pericentric inversions and/or other types of centromeric shift in the karyotype evolution of the group. When superimposed onto a phylogenetic tree based on molecular analyses of two mitochondrial genes the cytogenetic characteristics did not show any correlation with the phylogenetic relationships within the lineage. While karyotypes of many other Nothobranchius spp. studied to date diversified mainly via chromosome fusions and fissions, the N.ugandensis species group maintains stable 2n and the karyotype differentiation seems to be constrained to intrachromosomal rearrangements. Possible reasons for this difference in the trajectory of karyotype differentiation are discussed. While genetic drift seems to be a major factor in the fixation of chromosome rearrangements in Nothobranchius, future studies are needed to assess the impact of predicted multiple inversions on the genome evolution and species diversification within the N.ugandensis species group.

Global, regional, and cryptic population structure in a high gene-flow transatlantic fish

Lumpfish (Cyclopterus lumpus) is a transatlantic marine fish displaying large population sizes and a high potential for dispersal and gene-flow. These features are expected to result in weak population structure. Here, we investigated population genetic structure of lumpfish throughout its natural distribution in the North Atlantic using two approaches: I) 4,393 genome wide SNPs and 95 individuals from 10 locations, and II) 139 discriminatory SNPs and 1,669 individuals from 40 locations. Both approaches identified extensive population genetic structuring with a major split between the East and West Atlantic and a distinct Baltic Sea population, as well as further differentiation of lumpfish from the English Channel, Iceland, and Greenland. The discriminatory loci displayed ~2-5 times higher divergence than the genome wide approach, revealing further evidence of local population substructures. Lumpfish from Isfjorden in Svalbard were highly distinct but resembled most fish from Greenland. The Kattegat area in the Baltic transition zone, formed a previously undescribed distinct genetic group. Also, further subdivision was detected within North America, Iceland, West Greenland, Barents Sea, and Norway. Although lumpfish have considerable potential for dispersal and gene-flow, the observed high levels of population structuring throughout the Atlantic suggests that this species may have a natal homing behavior and local populations with adaptive differences. This fine-scale population structure calls for consideration when defining management units for exploitation of lumpfish stocks and in decisions related to sourcing and moving lumpfish for cleaner fish use in salmonid aquaculture.

Seascape genomics of common dolphins (Delphinus delphis) reveals adaptive diversity linked to regional and local oceanography

Background

High levels of standing genomic variation in wide-ranging marine species may enhance prospects for their long-term persistence. Patterns of connectivity and adaptation in such species are often thought to be influenced by spatial factors, environmental heterogeneity, and oceanographic and geomorphological features. Population-level studies that analytically integrate genome-wide data with environmental information (i.e., seascape genomics) have the potential to inform the spatial distribution of adaptive diversity in wide-ranging marine species, such as many marine mammals. We assessed genotype-environment associations (GEAs) in 214 common dolphins (Delphinus delphis) along > 3000 km of the southern coast of Australia.

Results

We identified 747 candidate adaptive SNPs out of a filtered panel of 17,327 SNPs, and five putatively locally-adapted populations with high levels of standing genomic variation were disclosed along environmentally heterogeneous coasts. Current velocity, sea surface temperature, salinity, and primary productivity were the key environmental variables associated with genomic variation. These environmental variables are in turn related to three main oceanographic phenomena that are likely affecting the dispersal of common dolphins: (1) regional oceanographic circulation, (2) localised and seasonal upwellings, and (3) seasonal on-shelf circulation in protected coastal habitats. Signals of selection at exonic gene regions suggest that adaptive divergence is related to important metabolic traits.

Conclusion

To the best of our knowledge, this represents the first seascape genomics study for common dolphins (genus Delphinus). Information from the associations between populations and their environment can assist population management in forecasting the adaptive capacity of common dolphins to climate change and other anthropogenic impacts.

How does genetic architecture affect eco-evolutionary dynamics? A theoretical perspective

Recent studies have revealed the importance of feedbacks between contemporary rapid evolution (i.e. evolution that occurs through changes in allele frequencies) and ecological dynamics. Despite its inherent interdisciplinary nature, however, studies on eco-evolutionary feedbacks have been mostly ecological and tended to focus on adaptation at the phenotypic level without considering the genetic architecture of evolutionary processes. In empirical studies, researchers have often compared ecological dynamics when the focal species under selection has a single genotype with dynamics when it has multiple genotypes. In theoretical studies, common approaches are models of quantitative traits where mean trait values change adaptively along the fitness gradient and Mendelian traits with two alleles at a single locus. On the other hand, it is well known that genetic architecture can affect short-term evolutionary dynamics in population genetics. Indeed, recent theoretical studies have demonstrated that genetic architecture (e.g. the number of loci, linkage disequilibrium and ploidy) matters in eco-evolutionary dynamics (e.g. evolutionary rescue where rapid evolution prevents extinction and population cycles driven by (co)evolution). I propose that theoretical approaches will promote the synthesis of functional genomics and eco-evolutionary dynamics through models that combine population genetics and ecology as well as nonlinear time-series analyses using emerging big data.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Repeatability of Voluntary Thermal Maximum and Covariance with Water Loss Reveal Potential for Adaptation to Changing Climates

AbstractAlthough climate warming poses a grave threat to amphibians, little is known about the capacity of this group to evolve in response to warming. The capacity of key traits to evolve depends on the presence of genetic variation on which selection can act. Here, we use repeatability estimates to estimate the potential upper bounds of heritable genetic variation in voluntary and critical thermal maxima of gray-cheeked salamanders (Plethodon metcalfi). Increases in thermal tolerance may also require concordant increases in resistance to water loss because hotter temperatures incur greater evaporative risk. Therefore, we also tested for a correlation between voluntary thermal maxima and resistance to water loss and conducted an acclimation study to test for covariation between these traits in response to warming. Voluntary thermal maxima exhibited low to moderate levels of repeatability (R=0.32, P=0.045), while critical thermal maxima exhibited no statistically significant repeatability (R=0.10, P=0.57). Voluntary thermal maxima also correlated positively with resistance to water loss (R=0.31, P=0.025) but only when controlling for body mass. Voluntary thermal maxima and resistance to water loss also exhibited different acclimatory responses across control (12°C-18°C) and warm (18°C-24°C) temperature regimes, indicating a potential decoupling of traits in different thermal environments. By addressing the repeatability of thermal tolerance and the potential for covariation with resistance to water loss, we begin to address some of the key requirements of amphibians to evolve in warming climates.

Invasion genomics uncover contrasting scenarios of genetic diversity in a widespread marine invader

Invasion rates have increased in the past 100 y irrespective of international conventions. What characterizes a successful invasion event? And how does genetic diversity translate into invasion success? Employing a whole-genome perspective using one of the most successful marine invasive species world-wide as a model, we resolve temporal invasion dynamics during independent invasion events in Eurasia. We reveal complex regionally independent invasion histories including cases of recurrent translocations, time-limited translocations, and stepping-stone range expansions with severe bottlenecks within the same species. Irrespective of these different invasion dynamics, which lead to contrasting patterns of genetic diversity, all nonindigenous populations are similarly successful. This illustrates that genetic diversity, per se, is not necessarily the driving force behind invasion success. Other factors such as propagule pressure and repeated introductions are an important contribution to facilitate successful invasions. This calls into question the dominant paradigm of the genetic paradox of invasions, i.e., the successful establishment of nonindigenous populations with low levels of genetic diversity.

Gene flow influences the genomic architecture of local adaptation in six riverine fish species

Understanding how gene flow influences adaptive divergence is important for predicting adaptive responses. Theoretical studies suggest that when gene flow is high, clustering of adaptive genes in fewer genomic regions would protect adaptive alleles from recombination and thus be selected for, but few studies have tested it with empirical data. Here, we used restriction site-associated sequencing to generate genomic data for six fish species with contrasting life histories from six reaches of the Upper Mississippi River System, USA. We used four differentiation-based outlier tests and three genotype-environment association analyses to define neutral single nucleotide polymorphisms (SNPs) and outlier SNPs that were putatively under selection. We then examined the distribution of outlier SNPs along the genome and investigated whether these SNPs were found in genomic islands of differentiation and inversions. We found that gene flow varied among species, and outlier SNPs were clustered more tightly in species with higher gene flow. The two species with the highest overall F_ST (0.0303-0.0720) and therefore lowest gene flow showed little evidence of clusters of outlier SNPs, with outlier SNPs in these species spreading uniformly across the genome. In contrast, nearly all outlier SNPs in the species with the lowest F_ST (0.0003) were found in a single large putative inversion. Two other species with intermediate gene flow (F_ST  ~ 0.0025-0.0050) also showed clustered genomic architectures, with most islands of differentiation clustered on a few chromosomes. Our results provide important empirical evidence to support the hypothesis that increasingly clustered architecture of local adaptation is associated with high gene flow.

Detecting selection in low-coverage high-throughput sequencing data using principal component analysis

BACKGROUND

Identification of selection signatures between populations is often an important part of a population genetic study. Leveraging high-throughput DNA sequencing larger sample sizes of populations with similar ancestries has become increasingly common. This has led to the need of methods capable of identifying signals of selection in populations with a continuous cline of genetic differentiation. Individuals from continuous populations are inherently challenging to group into meaningful units which is why existing methods rely on principal components analysis for inference of the selection signals. These existing methods require called genotypes as input which is problematic for studies based on low-coverage sequencing data.

MATERIALS AND METHODS

We have extended two principal component analysis based selection statistics to genotype likelihood data and applied them to low-coverage sequencing data from the 1000 Genomes Project for populations with European and East Asian ancestry to detect signals of selection in samples with continuous population structure.

RESULTS

Here, we present two selections statistics which we have implemented in the PCAngsd framework. These methods account for genotype uncertainty, opening for the opportunity to conduct selection scans in continuous populations from low and/or variable coverage sequencing data. To illustrate their use, we applied the methods to low-coverage sequencing data from human populations of East Asian and European ancestries and show that the implemented selection statistics can control the false positive rate and that they identify the same signatures of selection from low-coverage sequencing data as state-of-the-art software using high quality called genotypes.

CONCLUSION

We show that selection scans of low-coverage sequencing data of populations with similar ancestry perform on par with that obtained from high quality genotype data. Moreover, we demonstrate that PCAngsd outperform selection statistics obtained from called genotypes from low-coverage sequencing data without the need for ad-hoc filtering.

Locally Adaptive Inversions Modulate Genetic Variation at Different Geographic Scales in a Seaweed Fly

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales.

Expanding the conservation genomics toolbox: Incorporating structural variants to enhance genomic studies for species of conservation concern

Structural variants (SVs) are large rearrangements (>50 bp) within the genome that impact gene function and the content and structure of chromosomes. As a result, SVs are a significant source of functional genomic variation, that is, variation at genomic regions underpinning phenotype differences, that can have large effects on individual and population fitness. While there are increasing opportunities to investigate functional genomic variation in threatened species via single nucleotide polymorphism (SNP) data sets, SVs remain understudied despite their potential influence on fitness traits of conservation interest. In this future-focused Opinion, we contend that characterizing SVs offers the conservation genomics community an exciting opportunity to complement SNP-based approaches to enhance species recovery. We also leverage the existing literature-predominantly in human health, agriculture and ecoevolutionary biology-to identify approaches for readily characterizing SVs and consider how integrating these into the conservation genomics toolbox may transform the way we manage some of the world's most threatened species.

Local Adaptation in Marine Foundation Species at Microgeographic Scales

AbstractNearshore foundation species in coastal and estuarine systems (e.g., salt marsh grasses, mangroves, seagrasses, corals) drive the ecological functions of ecosystems and entire biomes by creating physical structure that alters local abiotic conditions and influences species interactions and composition. The resilience of foundation species and the ecosystem functions they provide depends on their phenotypic and genetic responses to spatial and temporal shifts in environmental conditions. In this review, we explore what is known about the causes and consequences of adaptive genetic differentiation in marine foundation species over spatial scales shorter than dispersal capabilities (i.e., microgeographic scales). We describe the strength of coupling field and laboratory experiments with population genetic techniques to illuminate patterns of local adaptation, and we illustrate this approach by using several foundation species. Among the major themes that emerge from our review include (1) adaptive differentiation of marine foundation species repeatedly evolves along vertical (i.e., elevation or depth) gradients, and (2) mating system and phenology may facilitate this differentiation. Microgeographic adaptation is an understudied mechanism potentially underpinning the resilience of many sessile marine species, and this evolutionary mechanism likely has particularly important consequences for the ecosystem functions provided by foundation species.

Adaptive and nonadaptive causes of heterogeneity in genetic differentiation across the Mimulus guttatus genome

Genetic diversity becomes structured among populations over time due to genetic drift and divergent selection. Although population structure is often treated as a uniform underlying factor, recent resequencing studies of wild populations have demonstrated that diversity in many regions of the genome may be structured quite dissimilar to the genome-wide pattern. Here, we explored the adaptive and nonadaptive causes of such genomic heterogeneity using population-level, whole genome resequencing data obtained from annual Mimulus guttatus individuals collected across a rugged environment landscape. We found substantial variation in how genetic differentiation is structured both within and between chromosomes, although, in contrast to other studies, known inversion polymorphisms appear to serve only minor roles in this heterogeneity. In addition, much of the genome can be clustered into eight among-population genetic differentiation patterns, but only two of these clusters are particularly consistent with patterns of isolation by distance. By performing genotype-environment association analysis, we also identified genomic intervals where local adaptation to specific climate factors has accentuated genetic differentiation among populations, and candidate genes in these windows indicate climate adaptation may proceed through changes affecting specialized metabolism, drought resistance, and development. Finally, by integrating our findings with previous studies, we show that multiple aspects of plant reproductive biology may be common targets of balancing selection and that variants historically involved in climate adaptation among populations have probably also fuelled rapid adaptation to microgeographic environmental variation within sites.

A beginner's guide to low-coverage whole genome sequencing for population genomics

Low-coverage whole genome sequencing (lcWGS) has emerged as a powerful and cost-effective approach for population genomic studies in both model and non-model species. However, with read depths too low to confidently call individual genotypes, lcWGS requires specialized analysis tools that explicitly account for genotype uncertainty. A growing number of such tools have become available, but it can be difficult to get an overview of what types of analyses can be performed reliably with lcWGS data, and how the distribution of sequencing effort between the number of samples analyzed and per-sample sequencing depths affects inference accuracy. In this introductory guide to lcWGS, we first illustrate how the per-sample cost for lcWGS is now comparable to RAD-seq and Pool-seq in many systems. We then provide an overview of software packages that explicitly account for genotype uncertainty in different types of population genomic inference. Next, we use both simulated and empirical data to assess the accuracy of allele frequency and genetic diversity estimation, detection of population structure, and selection scans under different sequencing strategies. Our results show that spreading a given amount of sequencing effort across more samples with lower depth per sample consistently improves the accuracy of most types of inference, with a few notable exceptions. Finally, we assess the potential for using imputation to bolster inference from lcWGS data in non-model species, and discuss current limitations and future perspectives for lcWGS-based population genomics research. With this overview, we hope to make lcWGS more approachable and stimulate its broader adoption.

Chromosome-Level Assembly of the Atlantic Silverside Genome Reveals Extreme Levels of Sequence Diversity and Structural Genetic Variation

The levels and distribution of standing genetic variation in a genome can provide a wealth of insights about the adaptive potential, demographic history, and genome structure of a population or species. As structural variants are increasingly associated with traits important for adaptation and speciation, investigating both sequence and structural variation is essential for wholly tapping this potential. Using a combination of shotgun sequencing, 10x Genomics linked reads and proximity-ligation data (Chicago and Hi-C), we produced and annotated a chromosome-level genome assembly for the Atlantic silverside (Menidia menidia)-an established ecological model for studying the phenotypic effects of natural and artificial selection-and examined patterns of genomic variation across two individuals sampled from different populations with divergent local adaptations. Levels of diversity varied substantially across each chromosome, consistently being highly elevated near the ends (presumably near telomeric regions) and dipping to near zero around putative centromeres. Overall, our estimate of the genome-wide average heterozygosity in the Atlantic silverside is among the highest reported for a fish, or any vertebrate (1.32-1.76% depending on inference method and sample). Furthermore, we also found extreme levels of structural variation, affecting ∼23% of the total genome sequence, including multiple large inversions (> 1 Mb and up to 12.6 Mb) associated with previously identified haploblocks showing strong differentiation between locally adapted populations. These extreme levels of standing genetic variation are likely associated with large effective population sizes and may help explain the remarkable adaptive divergence among populations of the Atlantic silverside.