skip to main content


Title: 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.  more » « less
Award ID(s):
1921585
NSF-PAR ID:
10171881
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
G3
ISSN:
2160-1836
Page Range / eLocation ID:
g3.401437.2020
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Studies of species that experience environmental heterogeneity across their distributions have become an important tool for understanding mechanisms of adaptation and predicting responses to climate change. We examine population structure, demographic history and environmentally associated genomic variation inBombus vosnesenskii, a common bumble bee in the western USA, using whole genome resequencing of populations distributed across a broad range of latitudes and elevations. We find thatB. vosnesenskiiexhibits minimal population structure and weak isolation by distance, confirming results from previous studies using other molecular marker types. Similarly, demographic analyses with Sequentially Markovian Coalescent models suggest that minimal population structure may have persisted since the last interglacial period, with genomes from different parts of the species range showing similar historical effective population size trajectories and relatively small fluctuations through time. Redundancy analysis revealed a small amount of genomic variation explained by bioclimatic variables. Environmental association analysis with latent factor mixed modelling (LFMM2) identified few outlier loci that were sparsely distributed throughout the genome and although a few putative signatures of selective sweeps were identified, none encompassed particularly large numbers of loci. Some outlier loci were in genes with known regulatory relationships, suggesting the possibility of weak selection, although compared with other species examined with similar approaches, evidence for extensive local adaptation signatures in the genome was relatively weak. Overall, results indicateB. vosnesenskiiis an example of a generalist with a high degree of flexibility in its environmental requirements that may ultimately benefit the species under periods of climate change.

     
    more » « less
  2. Abstract

    Understanding evolutionary responses to variation in temperature and precipitation across species ranges is of fundamental interest given ongoing climate change. The importance of temperature and precipitation for multiple aspects of bumble bee (Bombus) biology, combined with large geographic ranges that expose populations to diverse environmental pressures, make these insects well‐suited for studying local adaptation. Here, we analyzed genome‐wide sequence data from two widespread bumble bees,Bombus vosnesenskiiandBombus vancouverensis, using multiple environmental association analysis methods to investigate climate adaptation across latitude and altitude. The strongest signatures of selection were observed inB. vancouverensis, but despite unique responses between species for most loci, we detected several shared responses. Genes relating to neural and neuromuscular function and ion transport were especially evident with respect to temperature variables, while genes relating to cuticle formation, tracheal and respiratory system development, and homeostasis were associated with precipitation variables. Our data thus suggest that adaptive responses for tolerating abiotic variation are likely to be complex, but that several parallels among species can emerge even for these complex traits and landscapes. Results provide the framework for future work into mechanisms of thermal and desiccation tolerance in bumble bees and a set of genomic targets that might be monitored for future conservation efforts.

     
    more » « less
  3. Ware, Jessica (Ed.)
    Abstract Broadly distributed species experience divergent abiotic conditions across their ranges that may drive local adaptation. Montane systems where populations are distributed across both latitudinal and elevational gradients are especially likely to produce local adaptation due to spatial variation in multiple abiotic factors, including temperature, oxygen availability, and air density. We use whole-genome resequencing to evaluate the landscape genomics of Bombus vancouverensis Cresson (Hymenoptera: Apidae), a common montane bumble bee that is distributed throughout the western part of North America. Combined statistical approaches revealed several large windows of outlier SNPs with unusual levels of differentiation across the region and indicated that isothermality and elevation were the environmental features most strongly associated with these variants. Genes found within these regions had diverse biological functions, but included neuromuscular function, ion homeostasis, oxidative stress, and hypoxia that could be associated with tolerance of temperature, desiccation, or high elevation conditions. The whole-genome sequencing approach revealed outliers occurred in genome regions with elevated linkage disequilibrium, elevated mean FST, and low intrapopulation nucleotide diversity. Other kinds of structural variations were not widely associated with environmental predictors but did broadly match geographic separation. Results are consistent with other studies suggesting that regions of low recombination may harbor adaptive variation in bumble bees within as well as between species and refine our understanding of candidate genes that could be further investigated as possible targets of selection across the B. vancouverensis range. 
    more » « less
  4. Leppla, Norman (Ed.)
    Abstract

    Bombus vosnesenskii Radowszkowski, 1862 is one of three bumble bee species commercially available for pollination services in North America; however, little is documented about B. vosnesenskii colony life cycle or the establishment of ex situ rearing, mating, and overwintering practices. In this study, we documented nest success, colony size, and gyne production; recorded the duration of mating events; assessed overwintering survival of mated gynes; and evaluated second-generation nest success for colonies established from low- and high-elevation wild-caught B. vosnesenskii gynes. Of the 125 gynes installed, 62.4% produced brood cells (nest initiation) and 43.2% had at least 1 worker eclose (nest establishment). High-elevation B. vosnesenskii gynes had significantly higher nest initiation and establishment success than low-elevation gynes. However, low-elevation colonies were significantly larger with queens producing more gynes on average. Mating was recorded for 200 low-elevation and 37 high-elevation gynes, resulting in a mean duration of 62 and 51 min, respectively. Mated gynes were then placed into cold storage for 54 days to simulate overwintering, which resulted in 59.1% of low-elevation gynes surviving and 91.9% of high-elevation gynes surviving. For second-generation low-elevation gynes, 26.4% initiated nesting and 14.3% established nesting. Second-generation high-elevation gynes did not initiate nesting despite CO2 narcosis treatments. Overall, these results increase our understanding of B. vosnesenskii nesting, mating, and overwintering biology from 2 elevations. Furthermore, this study provides information on successful husbandry practices that can be used by researchers and conservationists to address knowledge gaps and enhance the captive rearing of bumble bees.

     
    more » « less
  5. Abstract

    As hybrid zones exhibit selective patterns of gene flow between otherwise distinct lineages, they can be especially valuable for informing processes of microevolution and speciation. The bumble bee,Bombus melanopygus, displays two distinct color forms generated by Müllerian mimicry: a northern “Rocky Mountain'’ color form with ferruginous mid‐abdominal segments (B.m.melanopygus) and a southern “Pacific'’ form with black mid‐abdominal segments (B.m.edwardsii). These morphs meet in a mimetic transition zone in northern California and southern Oregon that is more narrow and transitions further west than comimetic bumble bee species. To understand the historical formation of this mimicry zone, we assessed color distribution data forB.melanopygusfrom the last 100 years. We then examined gene flow among the color forms in the transition zone by comparing sequences from mitochondrial COI barcode sequences, color‐controlling loci, and the rest of the nuclear genome. These data support two geographically distinct mitochondrial haplogroups aligned to the ancestrally ferruginous and black forms that meet within the color transition zone. This clustering is also supported by the nuclear genome, which, while showing strong admixture across individuals, distinguishes individuals most by their mitochondrial haplotype, followed by geography. These data suggest the two lineages most likely were historically isolated, acquired fixed color differences, and then came into secondary contact with ongoing gene flow. The transition zone, however, exhibits asymmetries: mitochondrial haplotypes transition further south than color pattern, and both transition over shorter distances in the south. This system thus demonstrates alternative patterns of gene flow that occur in contact zones, presenting another example of mito‐nuclear discordance. Discordant gene flow is inferred to most likely be driven by a combination of mimetic selection, dominance effects, and assortative mating.

     
    more » « less