Population genomics of divergence among extreme and intermediate color forms in a polymorphic insect

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.

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.

Signals of adaptation to agricultural stress in the genomes of two European bumblebees

Human-induced environmental impacts on wildlife are widespread, causing major biodiversity losses. One major threat is agricultural intensification, typically characterised by large areas of monoculture, mechanical tillage, and the use of agrochemicals. Intensification leads to the fragmentation and loss of natural habitats, native vegetation, and nesting and breeding sites. Understanding the adaptability of insects to these changing environmental conditions is critical to predicting their survival. Bumblebees, key pollinators of wild and cultivated plants, are used as model species to assess insect adaptation to anthropogenic stressors. We investigated the effects of agricultural pressures on two common European bumblebees, Bombus pascuorum and B. lapidarius. Restriction-site Associated DNA Sequencing was used to identify loci under selective pressure across agricultural-natural gradients over 97 locations in Europe. 191 unique loci in B. pascuorum and 260 in B. lapidarius were identified as under selective pressure, and associated with agricultural stressors. Further investigation suggested several candidate proteins including several neurodevelopment, muscle, and detoxification proteins, but these have yet to be validated. These results provide insights into agriculture as a stressor for bumblebees, and signal for conservation action in light of ongoing anthropogenic changes.

Facilitating population genomics of non-model organisms through optimized experimental design for reduced representation sequencing

BACKGROUND

Genome-wide data are invaluable to characterize differentiation and adaptation of natural populations. Reduced representation sequencing (RRS) subsamples a genome repeatedly across many individuals. However, RRS requires careful optimization and fine-tuning to deliver high marker density while being cost-efficient. The number of genomic fragments created through restriction enzyme digestion and the sequencing library setup must match to achieve sufficient sequencing coverage per locus. Here, we present a workflow based on published information and computational and experimental procedures to investigate and streamline the applicability of RRS.

RESULTS

In an iterative process genome size estimates, restriction enzymes and size selection windows were tested and scaled in six classes of Antarctic animals (Ostracoda, Malacostraca, Bivalvia, Asteroidea, Actinopterygii, Aves). Achieving high marker density would be expensive in amphipods, the malacostracan target taxon, due to the large genome size. We propose alternative approaches such as mitogenome or target capture sequencing for this group. Pilot libraries were sequenced for all other target taxa. Ostracods, bivalves, sea stars, and fish showed overall good coverage and marker numbers for downstream population genomic analyses. In contrast, the bird test library produced low coverage and few polymorphic loci, likely due to degraded DNA.

CONCLUSIONS

Prior testing and optimization are important to identify which groups are amenable for RRS and where alternative methods may currently offer better cost-benefit ratios. The steps outlined here are easy to follow for other non-model taxa with little genomic resources, thus stimulating efficient resource use for the many pressing research questions in molecular ecology.

De Novo Genome Assemblies for Three North American Bumble Bee Species: Bombus bifarius, Bombus vancouverensis, and Bombus vosnesenskii

Bumble bees are ecologically and economically important insect pollinators. Three abundant and widespread species in western North America, Bombus bifarius, Bombus vancouverensis, and Bombus vosnesenskii, have been the focus of substantial research relating to diverse aspects of bumble bee ecology and evolutionary biology. We present de novo genome assemblies for each of the three species using hybrid assembly of Illumina and Oxford Nanopore Technologies sequences. All three assemblies are of high quality with large N50s (> 2.2 Mb), BUSCO scores indicating > 98% complete genes, and annotations producing 13,325 - 13,687 genes, comparing favorably with other bee genomes. Analysis of synteny against the most complete bumble bee genome, Bombus terrestris, reveals a high degree of collinearity. These genomes should provide a valuable resource for addressing questions relating to functional genomics and evolutionary biology in these species.

Bayesian Divergence-Time Estimation with Genome-Wide Single-Nucleotide Polymorphism Data of Sea Catfishes (Ariidae) Supports Miocene Closure of the Panamanian Isthmus

The closure of the Isthmus of Panama has long been considered to be one of the best defined biogeographic calibration points for molecular divergence-time estimation. However, geological and biological evidence has recently cast doubt on the presumed timing of the initial isthmus closure around 3 Ma but has instead suggested the existence of temporary land bridges as early as the Middle or Late Miocene. The biological evidence supporting these earlier land bridges was based either on only few molecular markers or on concatenation of genome-wide sequence data, an approach that is known to result in potentially misleading branch lengths and divergence times, which could compromise the reliability of this evidence. To allow divergence-time estimation with genomic data using the more appropriate multispecies coalescent (MSC) model, we here develop a new method combining the single-nucleotide polymorphism-based Bayesian species-tree inference of the software SNAPP with a molecular clock model that can be calibrated with fossil or biogeographic constraints. We validate our approach with simulations and use our method to reanalyze genomic data of Neotropical army ants (Dorylinae) that previously supported divergence times of Central and South American populations before the isthmus closure around 3 Ma. Our reanalysis with the MSC model shifts all of these divergence times to ages younger than 3 Ma, suggesting that the older estimates supporting the earlier existence of temporary land bridges were artifacts resulting at least partially from the use of concatenation. We then apply our method to a new restriction-site associated DNA-sequencing data set of Neotropical sea catfishes (Ariidae) and calibrate their species tree with extensive information from the fossil record. We identify a series of divergences between groups of Caribbean and Pacific sea catfishes around 10 Ma, indicating that processes related to the emergence of the isthmus led to vicariant speciation already in the Late Miocene, millions of years before the final isthmus closure.

Quaternary climate instability is correlated with patterns of population genetic variability in Bombus huntii

Climate oscillations have left a significant impact on the patterns of genetic diversity observed in numerous taxa. In this study, we examine the effect of Quaternary climate instability on population genetic variability of a bumble bee pollinator species, Bombus huntii in western North America. Pleistocene and contemporary B. huntii habitat suitability (HS) was estimated with an environmental niche model (ENM) by associating 1,035 locality records with 10 bioclimatic variables. To estimate genetic variability, we genotyped 380 individuals from 33 localities at 13 microsatellite loci. Bayesian inference was used to examine population structure with and without a priori specification of geographic locality. We compared isolation by distance (IBD) and isolation by resistance (IBR) models to examine population differentiation within and among the Bayesian inferred genetic clusters. Furthermore, we tested for the effect of environmental niche stability (ENS) on population genetic diversity with linear regression. As predicted, high-latitude B. huntii habitats exhibit low ENS when compared to low-latitude habitats. Two major genetic clusters of B. huntii inhabit western North America: (a) a north genetic cluster predominantly distributed north of 28°N and (b) a south genetic cluster distributed south of 28°N. In the south genetic cluser, both IBD and IBR models are significant. However, in the north genetic cluster, IBD is significant but not IBR. Furthermore, the IBR models suggest that low-latitude montane populations are surrounded by habitat with low HS, possibly limiting dispersal, and ultimately gene flow between populations. Finally, we detected high genetic diversity across populations in regions that have been climatically unstable since the last glacial maximum (LGM), and low genetic diversity across populations in regions that have been climatically stable since the LGM. Understanding how species have responded to climate change has the potential to inform management and conservation decisions of both ecological and economic concerns.

Loci under selection and markers associated with host plant and host-related strains shape the genetic structure of Brazilian populations of Spodoptera frugiperda (Lepidoptera, Noctuidae)

We applied the ddRAD genotyping-by-sequencing technique to investigate the genetic distinctiveness of Brazilian populations of the noctuid moth Spodoptera frugiperda, the fall armyworm (FAW), and the role of host-plant association as a source of genetic diversification. By strain-genotyping all field-collected individuals we found that populations collected from corn were composed primarily of corn-strain individuals, while the population collected from rice was composed almost entirely of rice-strain individuals. Outlier analyses indicated 1,184 loci putatively under selection (ca. 15% of the total) related to 194 different Gene Ontologies (GOs); the most numerous GOs were nucleotide binding, ATP binding, metal-ion binding and nucleic-acid binding. The association analyses indicated 326 loci associated with the host plant, and 216 loci associated with the individual strain, including functions related to Bacillus thuringiensis and insecticide resistance. The genetic-structure analyses indicated a moderate level of differentiation among all populations, and lower genetic structure among populations collected exclusively from corn, which suggests that the population collected from rice has a strong influence on the overall genetic structure. Populations of S. frugiperda are structured partially due to the host plant, and pairs of populations using the same host plant are more genetically similar than pairs using different hosts. Loci putatively under selection are the main factors responsible for the genetic structure of these populations, which indicates that adaptive selection on important traits, including the response to control tactics, is acting in the genetic differentiation of FAW populations in Brazil.

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.

Assessing genotype-phenotype associations in three dorsal colour morphs in the meadow spittlebug Philaenus spumarius (L.) (Hemiptera: Aphrophoridae) using genomic and transcriptomic resources

Background

Colour polymorphisms are common among animal species. When combined with genetic and ecological data, these polymorphisms can be excellent systems in which to understand adaptation and the molecular changes underlying phenotypic evolution. The meadow spittlebug, Philaenus spumarius (L.) (Hemiptera, Aphrophoridae), a widespread insect species in the Holarctic region, exhibits a striking dorsal colour/pattern balanced polymorphism. Although experimental crosses have revealed the Mendelian inheritance of this trait, its genetic basis remains unknown. In this study we aimed to identify candidate genomic regions associated with the colour balanced polymorphism in this species.

Results

By using restriction site-associated DNA (RAD) sequencing we were able to obtain a set of 1,837 markers across 33 individuals to test for associations with three dorsal colour phenotypes (typicus, marginellus, and trilineatus). Single and multi-association analyses identified a total of 60 SNPs associated with dorsal colour morphs. The genome size of P. spumarius was estimated by flow cytometry, revealing a 5.3 Gb genome, amongst the largest found in insects. A partial genome assembly, representing 24% of the total size, and an 81.4 Mb transcriptome, were also obtained. From the SNPs found to be associated with colour, 35% aligned to the genome and 10% to the transcriptome. Our data suggested that major loci, consisting of multi-genomic regions, may be involved in dorsal colour variation among the three dorsal colour morphs analysed. However, no homology was found between the associated loci and candidate genes known to be responsible for coloration pattern in other insect species. The associated markers showed stronger differentiation of the trilineatus colour phenotype, which has been shown previously to be more differentiated in several life-history and physiological characteristics as well. It is possible that colour variation and these traits are linked in a complex genetic architecture.

Conclusions

The loci detected to have an association with colour and the genomic and transcriptomic resources developed here constitute a basis for further research on the genetic basis of colour pattern in the meadow spittlebug P. spumarius.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0455-5) contains supplementary material, which is available to authorized users.

Genotyping of whole genome amplified reduced representation libraries reveals a cryptic population of Culicoides brevitarsis in the Northern Territory, Australia

BACKGROUND

The advent of genotyping by Next Generation Sequencing has enabled rapid discovery of thousands of single nucleotide polymorphism (SNP) markers and high throughput genotyping of large populations at an affordable cost. Genotyping by sequencing (GBS), a reduced representation library sequencing method, allows highly multiplexed sequencing of genomic subsets. This method has limitations for small organisms with low amounts of genomic DNA, such as the bluetongue virus (BTV) vectors, Culicoides midges.

RESULTS

This study employed the GBS method to isolate SNP markers de novo from whole genome amplified Culicoides brevitarsis genomic DNA. The individuals were collected from regions representing two different Australian patterns of BTV strain distribution: the Northern Territory (NT) and the east coast. We isolated 8145 SNPs using GBS. Phylogenetic analysis conducted using the filtered 3263 SNPs revealed the presence of a distinct C. brevitarsis sub-population in the NT and this was confirmed by analysis of mitochondrial DNA. Two loci showed a very strong signal for selection and were unique to the NT population. Bayesian analysis with STRUCTURE indicated a possible two-population cluster.

CONCLUSIONS

The results suggest that genotyping vectors with high density markers in combination with biological and environmental data is useful. However, more extensive sampling over a wider spatial and temporal range is needed. The presence of sub-structure in populations and loci under natural selection indicates the need for further investigation of the role of vectors in shaping the two Australian systems of BTV transmission. The described workflow is transferable to genotyping of small, non-model organisms, including arthropod vectors of pathogens of economic and medical importance.

Altitudinal variation in bumble bee (Bombus) critical thermal limits

Organism critical thermal limits are often tightly linked to current geographic distribution and can therefore help predict future range shifts driven by changing environmental temperatures. Thermal tolerance of diverse organisms often varies predictably with latitude, with upper thermal limits changing little and lower thermal limits decreasing with latitude. Despite similarly steep gradients in environmental temperatures across altitude, few studies have investigated altitudinal variation in critical thermal limits. We estimated critical thermal minimum (CTmin), critical thermal maximum (CTmax) and recovery temperature (Trec) by tracking righting response of three bumble bee species during thermal ramps: Bombus huntii collected from 2180m asl, and Bombus bifarius and Bombus sylvicola collected from 3290m asl in Wyoming, USA. Overall, larger bees could tolerate more extreme temperatures, likely due to a thermal inertia driven lag between core body temperatures and air temperatures. Despite their smaller size, high altitude bumble bees tolerated colder air temperatures: they had ~1°C lower CTmin and recovered from cold exposure at ~3-4°C lower air temperatures. Conversely, low altitude bees tolerated ~5°C hotter air temperatures. These altitudinal differences in thermal tolerance parallel differences in average daily minimum (1.2°C) and maximum (7.5°C) temperatures between these sites. These results provide one of the few measurements of organism thermal tolerance across altitude and the first evidence for geographical differences in tolerance of temperature extremes in heterothermic bumble bees.