skip to main content


Title: Divergence in Body Mass, Wing Loading, and Population Structure Reveals Species-Specific and Potentially Adaptive Trait Variation Across Elevations in Montane Bumble Bees
Abstract Biogeographic clines in morphology along environmental gradients can illuminate forces influencing trait evolution within and between species. Latitude has long been studied as a driver of morphological clines, with a focus on body size and temperature. However, counteracting environmental pressures may impose constraints on body size. In montane landscapes, declines in air density with elevation can negatively impact flight performance in volant species, which may contribute to selection for reduced body mass despite declining temperatures. We examine morphology in two bumble bee (Hymenoptera: Apidae: Bombus Latreille) species, Bombus vancouverensis Cresson and Bombus vosnesenskii Radoszkowski, across mountainous regions of California, Oregon, and Washington, United States. We incorporate population genomic data to investigate the relationship between genomic ancestry and morphological divergence. We find that B. vancouverensis, which tends to be more specialized for high elevations, exhibits stronger spatial-environmental variation, being smaller in the southern and higher elevation parts of its range and having reduced wing loading (mass relative to wing area) at high elevations. Bombus vosnesenskii, which is more of an elevational generalist, has substantial trait variation, but spatial-environmental correlations are weak. Population structure is stronger in the smaller B. vancouverensis, and we find a significant association between elevation and wing loading after accounting for genetic structure, suggesting the possibility of local adaptation for this flight performance trait. Our findings suggest that some conflicting results for body size trends may stem from distinct environmental pressures that impact different aspects of bumble bee ecology, and that different species show different morphological clines in the same region.  more » « less
Award ID(s):
1921562 1826834
NSF-PAR ID:
10298268
Author(s) / Creator(s):
; ; ; ; ; ; ;
Editor(s):
Hines, Heather
Date Published:
Journal Name:
Insect Systematics and Diversity
Volume:
5
Issue:
5
ISSN:
2399-3421
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. 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
  4. Abstract

    Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site‐associatedDNAsequencing (RADseq) in two bumble bee species,Bombus vosnesenskiiandBombus bifarius,across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A.Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure whileB. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, withB. vosnesenskiiexhibiting relatively consistent levels of genetic diversity across its range, whileB. bifariushas reduced genetic diversity at low latitudes, where it is restricted to high‐elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.

     
    more » « less
  5. Abstract

    Phenotypic polymorphism can constitute an inherent challenge for species delimitation. This issue is exemplified in bumble bees (Bombus), where species can exhibit high colour variation across their range, but otherwise exhibit little morphological variation to distinguish them from close relatives. We examine the species status of one of the most abundant North American bumble bees,Bombus bifariusCresson, which historically comprised two major taxa,bifariuss.s. andnearcticus. These lineages are recognized primarily by red and black variation in their mid‐abdominal coloration; however, a continuum from black (nearcticus) to red (bifariuss.s.) variation has led to their historic synonymization. Integrating mitochondrial and nuclear data and whole‐genome sequencing, we reveal a high level of both mitochondrial and nuclear divergence delimiting two morphologically cryptic species – the redbifariuss.s. and the colour‐variable (black to red)nearcticus. Population genomic analysis supports an absence of recent genomic admixture and a strong population structure between the two clades, even in sympatry. Species distribution models predict partially differentiated niches between the genetically inferred clades with annual precipitation being a leading differentiating variable. Thebifariuss.s. lineage also occupies significantly higher elevations, with regions of sympatry being among the highest elevations innearcticus. Our data also support a subspecies‐level divergence between the broadly distributednearcticusand the island populationvancouverensis. In this paper, we formally recognize the two species,Bombus bifariusCresson andBombus vancouverensisCresson, the latter including the subspeciesB. vancouverensis vancouverensiscomb.n.andB. vancouverensis nearcticuscomb.n., withvancouverensisthe name bearer due to year priority.

     
    more » « less