Revisiting comparisons of genetic diversity in stable and declining species: assessing genome-wide polymorphism in North American bumble bees using RAD sequencing

Effects of landscape, resource use, and body size on genetic structure in bee populations

Quantifying genetic structure and levels of genetic variation are fundamentally important to predicting the ability of populations to persist in human-altered landscapes and adapt to future environmental changes. Genetic structure reflects the dispersal of individuals over generations, which can be mediated by species-level traits or environmental factors. Dispersal distances are commonly positively associated with body size and negatively associated with the amount of degraded habitat between sites, motivating the investigation of these potential drivers of dispersal concomitantly. We quantified genetic structure and genetic variability within populations of seven bee species from the genus Euglossa across fragmented landscapes. We genotyped bees at SNP loci and tested the following predictions: (1) deforested areas restrict gene flow; (2) larger species have lower genetic structure; (3) species with greater resource specialization have higher genetic structure; and (4) sites surrounded by more intact habitat have higher genetic diversity. Contrasting with previous work on bees, we found no associations between body size and genetic structure. Genetic structure was higher for species with greater resource specialization, and the amount of intact habitat between or surrounding sites was positively associated with parameters reflecting gene flow and genetic diversity. These results challenge the dominant paradigm that individuals of larger species disperse farther, and they suggest that landscape and resource requirements are important factors mediating dispersal.

Population Genomics of the Critically Endangered Brazilian Merganser

The Brazilian merganser (Mergus octosetaceus) is one of the most endangered bird species in South America and comprises less than 250 mature individuals in wild environments. This is a species extremely sensitive to environmental disturbances and restricted to a few "pristine" freshwater habitats in Brazil, and it has been classified as Critically Endangered on the IUCN Red List since 1994. Thus, biological conservation studies are vital to promote adequate management strategies and to avoid the decline of merganser populations. In this context, to understand the evolutionary dynamics and the current genetic diversity of remaining Brazilian merganser populations, we used the "Genotyping by Sequencing" approach to genotype 923 SNPs in 30 individuals from all known areas of occurrence. These populations revealed a low genetic diversity and high inbreeding levels, likely due to the recent population decline associated with habitat loss. Furthermore, it showed a moderate level of genetic differentiation between all populations located in four separated areas of the highly threatened Cerrado biome. The results indicate that urgent actions for the conservation of the species should be accompanied by careful genetic monitoring to allow appropriate in situ and ex situ management to increase the long-term species' survival in its natural environment.

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.

Development of EST-SSR primers and genetic diversity analysis of the southern blight pathogen Sclerotium rolfsii using transcriptome data

Introduction

Sclerotium rolfsii Sacc. is a globally dispersed pathogenic fungus that causes southern blight disease in many crops and Chinese herbal medicine. The high degree of variation and diversity in the fungi altered population genetic structure. Therefore, the important factors of variation within the pathogen population should be considered during the development of management strategies for the disease.

Methods

In this study, S. rolfsii isolates from 13 hosts in 7 provinces of China were collected and analyzed to identify their morphological features and perform molecular characterization. To develop EST-SSR primers, transcriptome sequencing was performed on isolated CB1, and its SSR loci were comprehensively analyzed. In addition, we analyzed the polymorphisms among different populations based on screened EST-SSR primers.

Results

The results showed that all of these clean reads with total 36,165,475 assembled bases were clustered into 28,158 unigenes, ranged from 201 bp to 16,402 bp on the length, of which the average length was 1,284 bp. Of these, the SSR sequence appeared at an average interval of 15.43 kB, and the frequency of SSR was 0.0648 SSR/kB. Polymorphism of 9 primers was observed among 22 populations, and was verified by the Shannon's index (average = 1.414) and polymorphic information index (> 0.50). The genetic diversity analysis revealed diversity in all host populations and geographical populations. Further, molecular variance analysis (AMOVA) showed that the differences between groups were mainly related to geographical location. Based on cluster analysis, the 7 populations were roughly divided into 3 groups, and the results were highly consistent with those based on the geographical location, ultimately aligning with the results of STRUCTURE analysis.

Discussion

The findings build on current knowledge of the distribution of S. rolfsii in the southwest area of China, adding value to current knowledge base on the population structure and genetic diversity of S. rolfsii, specifically in the context of Chinese herbal medicine cultivation in China. Overall, our findings may provide valuable information for breeding of crops with enhanced resistance toward S. rolfsii.

Whole genome demographic models indicate divergent effective population size histories shape contemporary genetic diversity gradients in a montane bumble bee

Understanding historical range shifts and population size variation provides an important context for interpreting contemporary genetic diversity. Methods to predict changes in species distributions and model changes in effective population size (N e) using whole genomes make it feasible to examine how temporal dynamics influence diversity across populations. We investigate N e variation and climate-associated range shifts to examine the origins of a previously observed latitudinal heterozygosity gradient in the bumble bee Bombus vancouverensis Cresson (Hymenoptera: Apidae: Bombus Latreille) in western North America. We analyze whole genomes from a latitude-elevation cline using sequentially Markovian coalescent models of N e through time to test whether relatively low diversity in southern high-elevation populations is a result of long-term differences in N e. We use Maxent models of the species range over the last 130,000 years to evaluate range shifts and stability. N e fluctuates with climate across populations, but more genetically diverse northern populations have maintained greater N e over the late Pleistocene and experienced larger expansions with climatically favorable time periods. Northern populations also experienced larger bottlenecks during the last glacial period, which matched the loss of range area near these sites; however, bottlenecks were not sufficient to erode diversity maintained during periods of large N e. A genome sampled from an island population indicated a severe postglacial bottleneck, indicating that large recent postglacial declines are detectable if they have occurred. Genetic diversity was not related to niche stability or glacial-period bottleneck size. Instead, spatial expansions and increased connectivity during favorable climates likely maintain diversity in the north while restriction to high elevations maintains relatively low diversity despite greater stability in southern regions. Results suggest genetic diversity gradients reflect long-term differences in N e dynamics and also emphasize the unique effects of isolation on insular habitats for bumble bees. Patterns are discussed in the context of conservation under climate change.

Genetic diversity and structure of the 4th cycle breeding population of Chinese fir (Cunninghamia lanceolata (lamb.) hook)

Studying population genetic structure and diversity is crucial for the marker-assisted selection and breeding of coniferous tree species. In this study, using RAD-seq technology, we developed 343,644 high-quality single nucleotide polymorphism (SNP) markers to resolve the genetic diversity and population genetic structure of 233 Chinese fir selected individuals from the 4^th cycle breeding program, representing different breeding generations and provenances. The genetic diversity of the 4^th cycle breeding population was high with nucleotide diversity (P_i ) of 0.003, and H_o and H_e of 0.215 and 0.233, respectively, indicating that the breeding population has a broad genetic base. The genetic differentiation level between the different breeding generations and different provenances was low (F_st < 0.05), with population structure analysis results dividing the 233 individuals into four subgroups. Each subgroup has a mixed branch with interpenetration and weak population structure, which might be related to breeding rather than provenance, with aggregation from the same source only being in the local branches. Our results provide a reference for further research on the marker-assisted selective breeding of Chinese fir and other coniferous trees.

Genomic divergence, local adaptation, and complex demographic history may inform management of a popular sportfish species complex

The Neosho Bass (Micropterus velox), a former subspecies of the keystone top-predator and globally popular Smallmouth Bass (M. dolomieu), is endemic and narrowly restricted to small, clear streams of the Arkansas River Basin in the Central Interior Highlands (CIH) ecoregion, USA. Previous studies have detected some morphological, genetic, and genomic differentiation between the Neosho and Smallmouth Basses; however, the extent of neutral and adaptive divergence and patterns of intraspecific diversity are poorly understood. Furthermore, lineage diversification has likely been impacted by gene flow in some Neosho populations, which may be due to a combination of natural biogeographic processes and anthropogenic introductions. We assessed: (1) lineage divergence, (2) local directional selection (adaptive divergence), and (3) demographic history among Smallmouth Bass populations in the CIH using population genomic analyses of 50,828 single-nucleotide polymorphisms (SNPs) obtained through ddRAD-seq. Neosho and Smallmouth Bass formed monophyletic clades with 100% bootstrap support. We identified two major lineages within each species. We discovered six Neosho Bass populations (two nonadmixed and four admixed) and three nonadmixed Smallmouth Bass populations. We detected 29 SNPs putatively under directional selection in the Neosho range, suggesting populations may be locally adapted. Two populations were admixed via recent asymmetric secondary contact, perhaps after anthropogenic introduction. Two other populations were likely admixed via combinations of ancient and recent processes. These species comprise independently evolving lineages, some having experienced historical and natural admixture. These results may be critical for management of Neosho Bass as a distinct species and may aid in the conservation of other species with complex biogeographic histories.

Demographic inference provides insights into the extirpation and ecological dominance of eusocial snapping shrimps

Although eusocial animals often achieve ecological dominance in the ecosystems where they occur, many populations are unstable, resulting in local extinction. Both patterns may be linked to the characteristic demography of eusocial species-high reproductive skew and reproductive division of labor support stable effective population sizes that make eusocial groups more competitive in some species, but also lower effective population sizes that increase susceptibility to population collapse in others. Here, we examine the relationship between demography and social organization in Synalpheus snapping shrimps, a group in which eusociality has evolved recently and repeatedly. We show using coalescent demographic modelling that eusocial species have had lower but more stable effective population sizes across 100,000 generations. Our results are consistent with the idea that stable population sizes may enable competitive dominance in eusocial shrimps, but they also suggest that recent population declines are likely caused by eusocial shrimps' heightened sensitivity to environmental changes, perhaps as a result of their low effective population sizes and localized dispersal. Thus, although the unique life histories and demography of eusocial shrimps have likely contributed to their persistence and ecological dominance over evolutionary timescales, these social traits may also make them vulnerable to contemporary environmental change.

Genomic divergence and a lack of recent introgression between commercial and wild bumblebees (Bombus terrestris)

The global movement of bees for agricultural pollination services can affect local pollinator populations via hybridization. When commercial bumblebees are of the same species but of different geographic origin, intraspecific hybridization may result in beneficial integration of new genetic variation, or alternatively may disrupt locally adapted gene complexes. However, neither the existence nor the extent of genomic introgression and evolutionary divergence between wild and commercial bumblebees is fully understood. We obtained whole-genome sequencing data from wild and commercial Bombus terrestris collected from sites in Southern Sweden with and without long-term use of commercially imported B. terrestris. We search for evidence of introgression, dispersal and genome-wide differentiation in a comparative genomic analysis of wild and commercial bumblebees. Commercial B. terrestris were found in natural environments near sites where commercial bumblebees were used, as well as drifting wild B. terrestris in commercial bumblebee colonies. However, we found no evidence for widespread, recent genomic introgression of commercial B. terrestris into local wild conspecific populations. We found that wild B. terrestris had significantly higher nucleotide diversity (Nei's pi, π), while the number of segregating sites (Watterson's theta, θw) was higher in commercial B. terrestris. A highly divergent region on chromosome 11 was identified in commercial B. terrestris and found to be enriched with structural variants. The genes present in this region are involved in flight muscle contraction and structure and pathogen immune response, providing evidence for differing evolutionary processes operating in wild and commercial B. terrestris. We did not find evidence for recent introgression, suggesting that co-occurring commercial B. terrestris have not disrupted evolutionary processes in wild B. terrestris populations.

Integrating genomic and epidemiologic data to accelerate progress toward schistosomiasis elimination

The global community has adopted ambitious goals to eliminate schistosomiasis as a public health problem, and new tools are needed to achieve them. Mass drug administration programs, for example, have reduced the burden of schistosomiasis, but the identification of hotspots of persistent and reemergent transmission threaten progress toward elimination and underscore the need to couple treatment with interventions that reduce transmission. Recent advances in DNA sequencing technologies make whole-genome sequencing a valuable and increasingly feasible option for population-based studies of complex parasites such as schistosomes. Here, we focus on leveraging genomic data to tailor interventions to distinct social and ecological circumstances. We consider two priority questions that can be addressed by integrating epidemiological, ecological, and genomic information: (1) how often do non-human host species contribute to human schistosome infection? and (2) what is the importance of locally acquired versus imported infections in driving transmission at different stages of elimination? These questions address processes that can undermine control programs, especially those that rely heavily on treatment with praziquantel. Until recently, these questions were difficult to answer with sufficient precision to inform public health decision-making. We review the literature related to these questions and discuss how whole-genome approaches can identify the geographic and taxonomic sources of infection, and how such information can inform context-specific efforts that advance schistosomiasis control efforts and minimize the risk of reemergence.

Maine's Bumble Bees (Hymenoptera: Apidae)-Part 2: Comparisons of a Common (Bombus ternarius) and a Rare (Bombus terricola) Species

As part of a quantitative survey of Maine's bumble bee fauna (Butler et al. 2021), we compared and contrasted genetic diversity, parasite and pathogen burdens, and pesticide exposure of the relatively common Bombus ternarius Say, 1937 and the spatially rare Bombus terricola Kirby, 1837. We recorded 11 Bombus species at 40 survey sites across three Maine ecoregions, and B. ternarius was the most common species, while B. terricola was spatially rare. Nonmetric multidimensional scaling indicated that B. terricola was associated with higher elevation sites in Maine, while B. ternarius was more broadly distributed in the state. Pollinator networks constructed for each bee indicated B. ternarius foraged on more plant species than B. terricola, but that there was considerable overlap (73%) in plant species visited. Genetic diversity was greater in the spatially restricted B. terricola, whereas the widely distributed B. ternarius was characterized by greater genetic differentiation among regions. Bombus terricola had higher molecular marker levels of the microsporidian fungi Nosema spp. and the trypanosome Crithidia spp., and both species had high levels of Trypanosoma spp. exposure. No Western Honey Bee (Apis mellifera, Linnaeus, 1758) viruses were detected in either species. Pesticides were not detected in pollen samples collected from workers of either species, and B. ternarius worker tissue samples exhibited only trace levels of diflubenzuron.

Temporal Analysis of Effective Population Size and Mating System in a Social Wasp

Highly social species are successful because they cooperate in obligately integrated societies. We examined temporal genetic variation in the eusocial wasp Vespula maculifrons to gain a greater understanding of evolution in highly social taxa. First, we wished to test if effective population sizes of eusocial species were relatively low due to the reproductive division of labor that characterizes eusocial taxa. We thus estimated the effective population size of V. maculifrons by examining temporal changes in population allele frequencies. We sampled the genetic composition of a V. maculifrons population at 3 separate timepoints spanning a 13-year period. We found that effective population size ranged in the hundreds of individuals, which is similar to estimates in other, non-eusocial taxa. Second, we estimated levels of polyandry in V. maculifrons in different years to determine if queen mating system varied over time. We found no significant change in the number or skew of males mated to queens. In addition, mating skew was not significant within V. maculifrons colonies. Therefore, our data suggest that queen mate number may be subject to stabilizing selection in this taxon. Overall, our study provides novel insight into the selective processes operating in eusocial species by analyzing temporal genetic changes within populations.

Evidence for genetic isolation and local adaptation in the field cricket Gryllus campestris

Understanding how species can thrive in a range of environments is a central challenge for evolutionary ecology. There is strong evidence for local adaptation along large-scale ecological clines in insects. However, potential adaptation among neighbouring populations differing in their environment has been studied much less. We used RAD sequencing to quantify genetic divergence and clustering of ten populations of the field cricket Gryllus campestris in the Cantabrian Mountains of northern Spain, and an outgroup on the inland plain. Our populations were chosen to represent replicate high and low altitude habitats. We identified genetic clusters that include both high and low altitude populations indicating that the two habitat types do not hold ancestrally distinct lineages. Using common-garden rearing experiments to remove environmental effects, we found evidence for differences between high and low altitude populations in physiological and life-history traits. As predicted by the local adaptation hypothesis, crickets with parents from cooler (high altitude) populations recovered from periods of extreme cooling more rapidly than those with parents from warmer (low altitude) populations. Growth rates also differed between offspring from high and low altitude populations. However, contrary to our prediction that crickets from high altitudes would grow faster, the most striking difference was that at high temperatures, growth was fastest in individuals from low altitudes. Our findings reveal that populations a few tens of kilometres apart have independently evolved adaptations to their environment. This suggests that local adaptation in a range of traits may be commonplace even in mobile invertebrates at scales of a small fraction of species' distributions.

Conservation insights from wild bee genetic studies: Geographic differences, susceptibility to inbreeding, and signs of local adaptation

Conserving bees are critical both ecologically and economically. Genetic tools are valuable for monitoring these vital pollinators since tracking these small, fast-flying insects by traditional means is difficult. By surveying the current state of the literature, this review discusses how recent advances in landscape genetic and genomic research are elucidating how wild bees respond to anthropogenic threats. Current literature suggests that there may be geographic differences in the vulnerability of bee species to landscape changes. Populations of temperate bee species are becoming more isolated and more genetically depauperate as their landscape becomes more fragmented, but tropical bee species appear unaffected. These differences may be an artifact of historical differences in land-use, or it suggests that different management plans are needed for temperate and tropical bee species. Encouragingly, genetic studies on invasive bee species indicate that low levels of genetic diversity may not lead to rapid extinction in bees as once predicted. Additionally, next-generation sequencing has given researchers the power to identify potential genes under selection, which are likely critical to species' survival in their rapidly changing environment. While genetic studies provide insights into wild bee biology, more studies focusing on a greater phylogenetic and life-history breadth of species are needed. Therefore, caution should be taken when making broad conservation decisions based on the currently few species examined.

A reference-free approach to analyse RADseq data using standard next generation sequencing toolkits

Genotyping-by-sequencing methods such as RADseq are popular for generating genomic and population-scale data sets from a diverse range of organisms. These often lack a usable reference genome, restricting users to RADseq specific software for processing. However, these come with limitations compared to generic next generation sequencing (NGS) toolkits. Here, we describe and test a simple pipeline for reference-free RADseq data processing that blends de novo elements from STACKS with the full suite of state-of-the art NGS tools. Specifically, we use the de novo RADseq assembly employed by STACKS to create a catalogue of RAD loci that serves as a reference for read mapping, variant calling and site filters. Using RADseq data from 28 zebra sequenced to ~8x depth-of-coverage we evaluate our approach by comparing the site frequency spectra (SFS) to those from alternative pipelines. Most pipelines yielded similar SFS at 8x depth, but only a genotype likelihood based pipeline performed similarly at low sequencing depth (2-4x). We compared the RADseq SFS with medium-depth (~13x) shotgun sequencing of eight overlapping samples, revealing that the RADseq SFS was persistently slightly skewed towards rare and invariant alleles. Using simulations and human data we confirm that this is expected when there is allelic dropout (AD) in the RADseq data. AD in the RADseq data caused a heterozygosity deficit of ~16%, which dropped to ~5% after filtering AD. Hence, AD was the most important source of bias in our RADseq data.

Genetic Differentiation of a New World Screwworm Fly Population from Uruguay Detected by SNPs, Mitochondrial DNA and Microsatellites in Two Consecutive Years

The New World screwworm (NWS) fly, Cochliomyia hominivorax (Diptera: Calliphoridae), is an economically important ectoparasite currently distributed in South America and in the Caribbean basin. The successful eradication of this species in USA, Mexico and continental Central America was achieved by a control program based on the sterile insect technique (SIT). In order to implement a genetic control strategy over the NWS fly's current area of occurrence, first, it is necessary to understand the species dynamics and population structure. In order to address this objective, the spatial genetic structure of the NWS fly was previously reported in South America based on different genetic markers; however, to date, no study has investigated temporal changes in the genetic composition of its populations. In the current study, the temporal genetic structure of a NWS fly population from Uruguay was investigated through two consecutive samplings from the same locality over an interval of approximately 18 generations. The genetic structure was accessed with neutral and under selection SNPs obtained with genotyping-by-sequencing. The results gathered with these data were compared to estimates achieved with mitochondrial DNA sequences and eight microsatellite markers. Temporal changes in the genetic composition were revealed by all three molecular markers, which may be attributed to seasonal changes in the NWS fly's southern distribution. SNPs were employed for the first time for estimating the genetic structure in a NWS fly population; these results provide new clues and perspectives on its population genetic structure. This approach could have significant implications for the planning and implementation of management programs.

Cryptic Patterns of Speciation in Cryptic Primates: Microendemic Mouse Lemurs and the Multispecies Coalescent

Mouse lemurs (Microcebus) are a radiation of morphologically cryptic primates distributed throughout Madagascar for which the number of recognized species has exploded in the past two decades. This taxonomic revision has prompted understandable concern that there has been substantial oversplitting in the mouse lemur clade. Here, we investigate mouse lemur diversity in a region in northeastern Madagascar with high levels of microendemism and predicted habitat loss. We analyzed RADseq data with multispecies coalescent (MSC) species delimitation methods for two pairs of sister lineages that include three named species and an undescribed lineage previously identified to have divergent mtDNA. Marked differences in effective population sizes, levels of gene flow, patterns of isolation-by-distance, and species delimitation results were found among the two pairs of lineages. Whereas all tests support the recognition of the presently undescribed lineage as a separate species, the species-level distinction of two previously described species, M. mittermeieri and M. lehilahytsara is not supported-a result that is particularly striking when using the genealogical discordance index (gdi). Nonsister lineages occur sympatrically in two of the localities sampled for this study, despite an estimated divergence time of less than 1 Ma. This suggests rapid evolution of reproductive isolation in the focal lineages and in the mouse lemur clade generally. The divergence time estimates reported here are based on the MSC calibrated with pedigree-based mutation rates and are considerably more recent than previously published fossil-calibrated relaxed-clock estimates. We discuss the possible explanations for this discrepancy, noting that there are theoretical justifications for preferring the MSC estimates in this case. [Cryptic species; effective population size; microendemism; multispecies coalescent; speciation; species delimitation.].

Sex-specific impact of inbreeding on pathogen load in the striped dolphin

The impact of inbreeding on fitness has been widely studied and provides consequential inference about adaptive potential and the impact on survival for reduced and fragmented natural populations. Correlations between heterozygosity and fitness are common in the literature, but they rarely inform about the likely mechanisms. Here, we investigate a pathology with a clear impact on health in striped dolphin hosts (a nematode infection that compromises lung function). Dolphins varied with respect to their parasite burden of this highly pathogenic lung nematode (Skrjabinalius guevarai). Genetic diversity revealed by high-resolution restriction-associated DNA (43 018 RADseq single nucleotide polymorphisms) analyses showed a clear association between heterozygosity and pathogen load, but only for female dolphins, for which the more heterozygous individuals had lower Sk. guevarai burden. One locus identified by RADseq was a strong outlier in association with parasite load (heterozygous in all uninfected females, homozygous for 94% of infected females), found in an intron of the citron rho-interacting serine/threonine kinase locus (associated with milk production in mammals). Allelic variation at the Class II major histocompatability complex DQB locus was also assessed and found to be associated with both regional variation and with pathogen load. Both sex specificity and the identification of associating functional loci provide insight into the mechanisms by which more inbred individuals may be more susceptible to the infection of this parasite. This provides important insight towards our understanding of the impact of inbreeding in natural populations, relevant to both evolutionary and practical conservation considerations.

Comparative phylogeography of mainland and insular species of Neotropical molossid bats (Molossus)

Historical events, habitat preferences, and geographic barriers might result in distinct genetic patterns in insular versus mainland populations. Comparison between these two biogeographic systems provides an opportunity to investigate the relative role of isolation in phylogeographic patterns and to elucidate the importance of evolution and demographic history in population structure. Herein, we use a genotype-by-sequencing approach (GBS) to explore population structure within three species of mastiff bats (Molossus molossus, M. coibensis, and M. milleri), which represent different ecological histories and geographical distributions in the genus. We tested the hypotheses that oceanic straits serve as barriers to dispersal in Caribbean bats and that isolated island populations are more likely to experience genetic drift and bottlenecks in comparison with highly connected ones, thus leading to different phylogeographic patterns. We show that population structures vary according to general habitat preferences, levels of population isolation, and historical fluctuations in climate. In our dataset, mainland geographic barriers played only a small role in isolation of lineages. However, oceanic straits posed a partial barrier to the dispersal for some populations within some species (M. milleri), but do not seem to disrupt gene flow in others (M. molossus). Lineages on distant islands undergo genetic bottlenecks more frequently than island lineages closer to the mainland, which have a greater exchange of haplotypes.

Population genomics of Bombus terrestris reveals high but unstructured genetic diversity in a potential glacial refugium

Ongoing climate change is expected to cause an increase in temperature and a reduction of precipitation levels in the Mediterranean region, which might cause changes in many species distributions. These effects negatively influence species gene pools, decreasing genetic variability and adaptive potential. Here, we use mitochondrial DNA and RADseq to analyse population genetic structure and genetic diversity of the bumblebee species Bombus terrestris (subspecies Bombus terrestris lusitanicus), in the Iberian Peninsula. Although this subspecies shows a panmictic pattern of population structure across Iberia and beyond, we found differentiation between subspecies B. t. lusitanicus and B. t. africanus, probably caused by the existence of barriers to gene flow between Iberia and North Africa. Furthermore, the results revealed that the Iberian Peninsula harbours a large fraction of B. terrestris intraspecific genetic variation, with the highest number of mitochondrial haplotypes found when compared with any other region in Europe studied so far, suggesting a potential role for the Iberian Peninsula as a glacial refugium. Our findings strengthen the idea that Iberia is a very important source of diversity for the global genetic pool of this species, because rare alleles might play a role in population resilience against human- or climate-mediated changes.

Defining Pollinator Health: A Holistic Approach Based on Ecological, Genetic, and Physiological Factors

Evidence for global bee population declines has catalyzed a rapidly evolving area of research that aims to identify the causal factors and to effectively assess the status of pollinator populations. The term pollinator health emerged through efforts to understand causes of bee decline and colony losses, but it lacks a formal definition. In this review, we propose a definition for pollinator health and synthesize the available literature on the application of standardized biomarkers to assess health at the individual, colony, and population levels. We focus on biomarkers in honey bees, a model species, but extrapolate the potential application of these approaches to monitor the health status of wild bee populations. Biomarker-guided health measures can inform beekeeper management decisions, wild bee conservation efforts, and environmental policies. We conclude by addressing challenges to pollinator health from a One Health perspective that emphasizes the interplay between environmental quality and human, animal, and bee health.

Macrogeographic genetic structure of Lutzomyia longipalpis complex populations using Next Generation Sequencing

Lutzomyia longipalpis is the main vector of Leishmania infantum, the causative agent of visceral leishmaniasis in the Neotropical realm. Its taxonomic status has been widely discussed once it encompasses a complex of species. The knowledge about the genetic structure of insect vector populations helps the elucidation of components and interactions of the disease ecoepidemiology. Thus, the objective of this study was to genotypically analyze populations of the Lu. longipalpis complex from a macrogeographic perspective using Next Generation Sequencing. Polymorphism analysis of three molecular markers was used to access the levels of population genetic structure among nine different populations of sand flies. Illumina Amplicon Sequencing Protocol^® was used to identify possible polymorphic sites. The library was sequenced on paired-end Illumina MiSeq platform. Significant macrogeographical population differentiation was observed among Lu. longipalpis populations via PCA and DAPC analyses. Our results revealed that populations of Lu. longipalpis from the nine municipalities were grouped into three clusters. In addition, it was observed that the levels of Lu. longipalpis population structure could be associated with distance isolation. This new sequencing method allowed us to study different molecular markers after a single sequencing run, and to evaluate population and inter-species differences on a macrogeographic scale.

Pattern of population structuring between Belgian and Estonian bumblebees

Several population genetic studies investigated the extent of gene flow and population connectivity in bumblebees. In general, no restriction in gene flow is considered for mainland populations of common bumblebee species. Whether this assumption holds true for all species is not known. An assessment of bumblebee genetic structure in the context of their geographic distribution is needed to prioritize conservation and management needs. Here, we conducted a genetic study on seven bumblebee species occurring in Belgium and Estonia. Using 16 microsatellite markers, we investigated genetic diversity and population structuring in each species. This is the first study investigating population structuring of both declining and stable bumblebee species on both small and large geographic scales. Our results showed no or only low population structuring between the populations of the restricted and declining bumblebee species on both scales, while significant structuring was found for populations of the common species on the larger scale. The latter result, which may be due to human or environmental changes in the landscape, implies the need for the conservation of also widespread bumblebee species. Conservation strategies to improve gene flow and connectivity of populations could avoid the isolation and future losses of populations of these important species.

Influence of a climatic gradient on genetic exchange between two oak species

PREMISE

In plant groups with limited intrinsic barriers to gene flow, it is thought that environmental conditions can modulate interspecific genetic exchange. Oaks are known for limited barriers to gene flow among closely related species. Here, we use Quercus as a living laboratory in which to pursue a fundamental question in plant evolution: Do environmental gradients restrict or promote genetic exchange between species?

METHODS

We focused on two North American oaks, the rare Quercus dumosa and the widespread Q. berberidifolia. We sampled intensively along a contact zone in California, USA. We sequenced restriction site-associated DNA markers and measured vegetative phenotype. We tested for genetic exchange, the association with climate, and the effect on phenotype.

RESULTS

There is evidence for genetic exchange between the species. Admixed plants are found in areas of intermediate climate, while less admixed plants are found at the extremes of the climatic gradient. Genetic and phenotypic patterns are out of phase in the contact zone; some plants display the phenotype of one species but are genetically associated with another.

CONCLUSIONS

Our results support the hypothesis that a strong climatic gradient can promote genetic exchange between species. The overall weak correlation between genotype and phenotype in the contact zone between the species suggests that genetic exchange can lead to the breakdown of trait combinations used to define species. This incongruency predicts ongoing problems for conservation of Q. dumosa, with implications for conservation of other oaks.

Landscape genomics to the rescue of a tropical bee threatened by habitat loss and climate change

Habitat degradation and climate change are currently threatening wild pollinators, compromising their ability to provide pollination services to wild and cultivated plants. Landscape genomics offers powerful tools to assess the influence of landscape modifications on genetic diversity and functional connectivity, and to identify adaptations to local environmental conditions that could facilitate future bee survival. Here, we assessed range-wide patterns of genetic structure, genetic diversity, gene flow, and local adaptation in the stingless bee Melipona subnitida, a tropical pollinator of key biological and economic importance inhabiting one of the driest and hottest regions of South America. Our results reveal four genetic clusters across the species' full distribution range. All populations were found to be under a mutation-drift equilibrium, and genetic diversity was not influenced by the amount of reminiscent natural habitats. However, genetic relatedness was spatially autocorrelated and isolation by landscape resistance explained range-wide relatedness patterns better than isolation by geographic distance, contradicting earlier findings for stingless bees. Specifically, gene flow was enhanced by increased thermal stability, higher forest cover, lower elevations, and less corrugated terrains. Finally, we detected genomic signatures of adaptation to temperature, precipitation, and forest cover, spatially distributed in latitudinal and altitudinal patterns. Taken together, our findings shed important light on the life history of M. subnitida and highlight the role of regions with large thermal fluctuations, deforested areas, and mountain ranges as dispersal barriers. Conservation actions such as restricting long-distance colony transportation, preserving local adaptations, and improving the connectivity between highlands and lowlands are likely to assure future pollination services.

Is genomic diversity a useful proxy for census population size? Evidence from a species-rich community of desert lizards

Species abundance data are critical for testing ecological theory, but obtaining accurate empirical estimates for many taxa is challenging. Proxies for species abundance can help researchers circumvent time and cost constraints that are prohibitive for long-term sampling. Under simple demographic models, genetic diversity is expected to correlate with census size, such that genome-wide heterozygosity may provide a surrogate measure of species abundance. We tested whether nucleotide diversity is correlated with long-term estimates of abundance, occupancy and degree of ecological specialization in a diverse lizard community from arid Australia. Using targeted sequence capture, we obtained estimates of genomic diversity from 30 species of lizards, recovering an average of 5,066 loci covering 3.6 Mb of DNA sequence per individual. We compared measures of individual heterozygosity to a metric of habitat specialization to investigate whether ecological preference exerts a measurable effect on genetic diversity. We find that heterozygosity is significantly correlated with species abundance and occupancy, but not habitat specialization. Demonstrating the power of genomic sampling, the correlation between heterozygosity and abundance/occupancy emerged from considering just one or two individuals per species. However, genetic diversity does no better at predicting abundance than a single day of traditional sampling in this community. We conclude that genetic diversity is a useful proxy for regional-scale species abundance and occupancy, but a large amount of unexplained variation in heterozygosity suggests additional constraints or a failure of ecological sampling to adequately capture variation in true population size.

Conservation Genomics of the Declining North American Bumblebee Bombus terricola Reveals Inbreeding and Selection on Immune Genes

The yellow-banded bumblebee Bombus terricola was common in North America but has recently declined and is now on the IUCN Red List of threatened species. The causes of B. terricola's decline are not well understood. Our objectives were to create a partial genome and then use this to estimate population data of conservation interest, and to determine whether genes showing signs of recent selection suggest a specific cause of decline. First, we generated a draft partial genome (contig set) for B. terricola, sequenced using Pacific Biosciences RS II at an average depth of 35×. Second, we sequenced the individual genomes of 22 bumblebee gynes from Ontario and Quebec using Illumina HiSeq 2500, each at an average depth of 20×, which were used to improve the PacBio genome calls and for population genetic analyses. The latter revealed that several samples had long runs of homozygosity, and individuals had high inbreeding coefficient F, consistent with low effective population size. Our data suggest that B. terricola's effective population size has decreased orders of magnitude from pre-Holocene levels. We carried out tests of selection to identify genes that may have played a role in ameliorating environmental stressors underlying B. terricola's decline. Several immune-related genes have signatures of recent positive selection, which is consistent with the pathogen-spillover hypothesis for B. terricola's decline. The new B. terricola contig set can help solve the mystery of bumblebee decline by enabling functional genomics research to directly assess the health of pollinators and identify the stressors causing declines.

Genomic signatures of population bottleneck and recovery in Northwest Atlantic pinnipeds

Population increases over the past several decades provide natural settings in which to study the evolutionary processes that occur during bottleneck, growth, and spatial expansion. We used parallel natural experiments of historical decline and subsequent recovery in two sympatric pinniped species in the Northwest Atlantic, the gray seal (Halichoerus grypus atlantica) and harbor seal (Phoca vitulina vitulina), to study the impact of recent demographic change in genomic diversity. Using restriction site-associated DNA sequencing, we assessed genomic diversity at over 8,700 polymorphic gray seal loci and 3,700 polymorphic harbor seal loci in samples from multiple cohorts collected throughout recovery over the past half-century. Despite significant differences in the degree of genetic diversity assessed in the two species, we found signatures of historical bottlenecks in the contemporary genomes of both gray and harbor seals. We evaluated temporal trends in diversity across cohorts, as well as compared samples from sites at both the center and edge of a recent gray seal range expansion, but found no significant change in genomewide diversity following recovery. We did, however, find that the variance and degree of allele frequency change measured over the past several decades were significantly different from neutral expectations of drift under population growth. These two cases of well-described demographic history provide opportunities for critical evaluation of current approaches to simulating and understanding the genetic effects of historical demographic change in natural populations.

Genomic associations with bill length and disease reveal drift and selection across island bird populations

Island species provide excellent models for investigating how selection and drift operate in wild populations, and for determining how these processes act to influence local adaptation and speciation. Here, we examine the role of selection and drift in shaping genomic and phenotypic variation across recently separated populations of Berthelot's pipit (Anthus berthelotii), a passerine bird endemic to three archipelagos in the Atlantic. We first characterized genetic diversity and population structuring that supported previous inferences of a history of recent colonizations and bottlenecks. We then tested for regions of the genome associated with the ecologically important traits of bill length and malaria infection, both of which vary substantially across populations in this species. We identified a SNP associated with variation in bill length among individuals, islands, and archipelagos; patterns of variation at this SNP suggest that both phenotypic and genotypic variation in bill length is largely shaped by founder effects. Malaria was associated with SNPs near/within genes involved in the immune response, but this relationship was not consistent among archipelagos, supporting the view that disease resistance is complex and rapidly evolving. Although we found little evidence for divergent selection at candidate loci for bill length and malaria resistance, genome scan analyses pointed to several genes related to immunity and metabolism as having important roles in divergence and adaptation. Our findings highlight the utility and challenges involved with combining association mapping and population genetic analysis in nonequilibrium populations, to disentangle the effects of drift and selection on shaping genotypes and phenotypes.

Genome-wide analyses of the Bemisia tabaci species complex reveal contrasting patterns of admixture and complex demographic histories

Once considered a single species, the whitefly, Bemisia tabaci, is a complex of numerous morphologically indistinguishable species. Within the last three decades, two of its members (MED and MEAM1) have become some of the world's most damaging agricultural pests invading countries across Europe, Africa, Asia and the Americas and affecting a vast range of agriculturally important food and fiber crops through both feeding-related damage and the transmission of numerous plant viruses. For some time now, researchers have relied on a single mitochondrial gene and/or a handful of nuclear markers to study this species complex. Here, we move beyond this by using 38,041 genome-wide Single Nucleotide Polymorphisms, and show that the two invasive members of the complex are closely related species with signatures of introgression with a third species (IO). Gene flow patterns were traced between contemporary invasive populations within MED and MEAM1 species and these were best explained by recent international trade. These findings have profound implications for delineating the B. tabaci species status and will impact quarantine measures and future management strategies of this global pest.

Genomic data reveal a loss of diversity in two species of tuco-tucos (genus Ctenomys) following a volcanic eruption

Marked reductions in population size can trigger corresponding declines in genetic variation. Understanding the precise genetic consequences of such reductions, however, is often challenging due to the absence of robust pre- and post-reduction datasets. Here, we use heterochronous genomic data from samples obtained before and immediately after the 2011 eruption of the Puyehue-Cordón Caulle volcanic complex in Patagonia to explore the genetic impacts of this event on two parapatric species of rodents, the colonial tuco-tuco (Ctenomys sociabilis) and the Patagonian tuco-tuco (C. haigi). Previous analyses using microsatellites revealed no post-eruption changes in genetic variation in C. haigi, but an unexpected increase in variation in C. sociabilis. To explore this outcome further, we used targeted gene capture to sequence over 2,000 putatively neutral regions for both species. Our data revealed that, contrary to the microsatellite analyses, the eruption was associated with a small but significant decrease in genetic variation in both species. We suggest that genome-level analyses provide greater power than traditional molecular markers to detect the genetic consequences of population size changes, particularly changes that are recent, short-term, or modest in size. Consequently, genomic analyses promise to generate important new insights into the effects of specific environmental events on demography and genetic variation.

Population genomics reveals a candidate gene involved in bumble bee pigmentation

Variation in bumble bee color patterns is well‐documented within and between species. Identifying the genetic mechanisms underlying such variation may be useful in revealing evolutionary forces shaping rapid phenotypic diversification. The widespread North American species Bombus bifarius exhibits regional variation in abdominal color forms, ranging from red‐banded to black‐banded phenotypes and including geographically and phenotypically intermediate forms. Identifying genomic regions linked to this variation has been complicated by strong, near species level, genome‐wide differentiation between red‐ and black‐banded forms. Here, we instead focus on the closely related black‐banded and intermediate forms that both belong to the subspecies B. bifarius nearcticus. We analyze an RNA sequencing (RNAseq) data set and identify a cluster of single nucleotide polymorphisms (SNPs) within one gene, Xanthine dehydrogenase/oxidase‐like, that exhibit highly unusual differentiation compared to the rest of the sequenced genome. Homologs of this gene contribute to pigmentation in other insects, and results thus represent a strong candidate for investigating the genetic basis of pigment variation in B. bifarius and other bumble bee mimicry complexes.

Genomic data detect corresponding signatures of population size change on an ecological time scale in two salamander species

Understanding the demography of species over recent history (e.g. <100 years) is critical in studies of ecology and evolution, but records of population history are rarely available. Surveying genetic variation is a potential alternative to census-based estimates of population size, and can yield insight into the demography of a population. However, to assess the performance of genetic methods, it is important to compare their estimates of population history to known demography. Here, we leveraged the exceptional resources from a wetland with 37 years of amphibian mark-recapture data to study the utility of genetically based demographic inference on salamander species with documented population declines (Ambystoma talpoideum) and expansions (A. opacum), patterns that have been shown to be correlated with changes in wetland hydroperiod. We generated ddRAD data from two temporally sampled populations of A. opacum (1993, 2013) and A. talpoideum (1984, 2011) and used coalescent-based demographic inference to compare alternate evolutionary models. For both species, demographic model inference supported population size changes that corroborated mark-recapture data. Parameter estimation in A. talpoideum was robust to our variations in analytical approach, while estimates for A. opacum were highly inconsistent, tempering our confidence in detecting a demographic trend in this species. Overall, our robust results in A. talpoideum suggest that genome-based demographic inference has utility on an ecological scale, but researchers should also be cognizant that these methods may not work in all systems and evolutionary scenarios. Demographic inference may be an important tool for population monitoring and conservation management planning.

Estimating genomic diversity and population differentiation - an empirical comparison of microsatellite and SNP variation in Arabidopsis halleri

BACKGROUND

Microsatellite markers are widely used for estimating genetic diversity within and differentiation among populations. However, it has rarely been tested whether such estimates are useful proxies for genome-wide patterns of variation and differentiation. Here, we compared microsatellite variation with genome-wide single nucleotide polymorphisms (SNPs) to assess and quantify potential marker-specific biases and derive recommendations for future studies. Overall, we genotyped 180 Arabidopsis halleri individuals from nine populations using 20 microsatellite markers. Twelve of these markers were originally developed for Arabidopsis thaliana (cross-species markers) and eight for A. halleri (species-specific markers). We further characterized 2 million SNPs across the genome with a pooled whole-genome re-sequencing approach (Pool-Seq).

RESULTS

Our analyses revealed that estimates of genetic diversity and differentiation derived from cross-species and species-specific microsatellites differed substantially and that expected microsatellite heterozygosity (SSR-H e) was not significantly correlated with genome-wide SNP diversity estimates (SNP-H e and θ Watterson) in A. halleri. Instead, microsatellite allelic richness (A r) was a better proxy for genome-wide SNP diversity. Estimates of genetic differentiation among populations (F ST) based on both marker types were correlated, but microsatellite-based estimates were significantly larger than those from SNPs. Possible causes include the limited number of microsatellite markers used, marker ascertainment bias, as well as the high variance in microsatellite-derived estimates. In contrast, genome-wide SNP data provided unbiased estimates of genetic diversity independent of whether genome- or only exome-wide SNPs were used. Further, we inferred that a few thousand random SNPs are sufficient to reliably estimate genome-wide diversity and to distinguish among populations differing in genetic variation.

CONCLUSIONS

We recommend that future analyses of genetic diversity within and differentiation among populations use randomly selected high-throughput sequencing-based SNP data to draw conclusions on genome-wide diversity patterns. In species comparable to A. halleri, a few thousand SNPs are sufficient to achieve this goal.

Application of RAD Sequencing for Evaluating the Genetic Diversity of Domesticated Panax notoginseng (Araliaceae)

Panax notoginseng, a traditional Chinese medicinal plant, has been cultivated and domesticated for approximately 400 years, mainly in Yunnan and Guangxi, two provinces in southwest China. This species was named according to cultivated rather than wild individuals, and no wild populations had been found until now. The genetic resources available on farms are important for both breeding practices and resource conservation. In the present study, the recently developed technology RADseq, which is based on next-generation sequencing, was used to analyze the genetic variation and differentiation of P. notoginseng. The nucleotide diversity and heterozygosity results indicated that P. notoginseng had low genetic diversity at both the species and population levels. Almost no genetic differentiation has been detected, and all populations were genetically similar due to strong gene flow and insufficient splitting time. Although the genetic diversity of P. notoginseng was low at both species and population levels, several traditional plantations had relatively high genetic diversity, as revealed by the He and π values and by the private allele numbers. These valuable genetic resources should be protected as soon as possible to facilitate future breeding projects. The possible geographical origin of Sanqi domestication was discussed based on the results of the genetic diversity analysis.

Population genetic inferences using immune gene SNPs mirror patterns inferred by microsatellites

Single nucleotide polymorphisms (SNPs) are replacing microsatellites for population genetic analyses, but it is not apparent how many SNPs are needed or how well SNPs correlate with microsatellites. We used data from the gopher tortoise, Gopherus polyphemus-a species with small populations, to compare SNPs and microsatellites to estimate population genetic parameters. Specifically, we compared one SNP data set (16 tortoises from four populations sequenced at 17 901 SNPs) to two microsatellite data sets, a full data set of 101 tortoises and a partial data set of 16 tortoises previously genotyped at 10 microsatellites. For the full microsatellite data set, observed heterozygosity, expected heterozygosity and F_ST were correlated between SNPs and microsatellites; however, allelic richness was not. The same was true for the partial microsatellite data set, except that allelic richness, but not observed heterozygosity, was correlated. The number of clusters estimated by structure differed for each data set (SNPs = 2; partial microsatellite = 3; full microsatellite = 4). Principle component analyses (PCA) showed four clusters for all data sets. More than 800 SNPs were needed to correlate with allelic richness, observed heterozygosity and expected heterozygosity, but only 100 were needed for F_ST . The number of SNPs typically obtained from next-generation sequencing (NGS) far exceeds the number needed to correlate with microsatellite parameter estimates. Our study illustrates that diversity, F_ST and PCA results from microsatellites can mirror those obtained with SNPs. These results may be generally applicable to small populations, a defining feature of endangered and threatened species, because theory predicts that genetic drift will tend to outweigh selection in small populations.

Genetic variability, local selection and demographic history: genomic evidence of evolving towards allopatric speciation in Asian seabass

Genomewide analysis of genetic divergence is critically important in understanding the genetic processes of allopatric speciation. We sequenced RAD tags of 131 Asian seabass individuals of six populations from South-East Asia and Australia/Papua New Guinea. Using 32 433 SNPs, we examined the genetic diversity and patterns of population differentiation across all the populations. We found significant evidence of genetic heterogeneity between South-East Asian and Australian/Papua New Guinean populations. The Australian/Papua New Guinean populations showed a rather lower level of genetic diversity. FST and principal components analysis revealed striking divergence between South-East Asian and Australian/Papua New Guinean populations. Interestingly, no evidence of contemporary gene flow was observed. The demographic history was further tested based on the folded joint site frequency spectrum. The scenario of ancient migration with historical population size changes was suggested to be the best fit model to explain the genetic divergence of Asian seabass between South-East Asia and Australia/Papua New Guinea. This scenario also revealed that Australian/Papua New Guinean populations were founded by ancestors from South-East Asia during mid-Pleistocene and were completely isolated from the ancestral population after the last glacial retreat. We also detected footprints of local selection, which might be related to differential ecological adaptation. The ancient gene flow was examined and deemed likely insufficient to counteract the genetic differentiation caused by genetic drift. The observed genomic pattern of divergence conflicted with the 'genomic islands' scenario. Altogether, Asian seabass have likely been evolving towards allopatric speciation since the split from the ancestral population during mid-Pleistocene.

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction-site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small Ne (range = 2.1-89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome-wide divergence. Nonetheless, outlier tests identified 3.6-6.6% of loci as high FST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high FST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small Ne in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.