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-associated DNA sequencing (RADseq) in two bumble bee species, Bombus vosnesenskii and Bombus 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 while B. 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, with B. vosnesenskii exhibiting relatively consistent levels of genetic diversity across its range, while B. bifarius has 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.

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 vosnesenskii and Bombus vancouverensis, using multiple environmental association analysis methods to investigate climate adaptation across latitude and altitude. The strongest signatures of selection were observed in B. 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.

Flies transport specific bacteria with their larvae that provide a wider range of nutrients for those bacteria. Our hypothesis was that this symbiotic interaction may depend on interkingdom signaling. We obtained Proteus mirabilis from the salivary glands of the blow fly Lucilia sericata; this strain swarmed significantly and produced a strong odor that attracts blow flies. To identify the putative interkingdom signals for the bacterium and flies, we reasoned that as swarming is used by this bacterium to cover the food resource and requires bacterial signaling, the same bacterial signals used for swarming may be used to communicate with blow flies. Using transposon mutagenesis, we identified six novel genes for swarming (ureR, fis, hybG, zapB, fadE and PROSTU_03490), then, confirming our hypothesis, we discovered that fly attractants, lactic acid, phenol, NaOH, KOH and ammonia, restore swarming for cells with the swarming mutations. Hence, compounds produced by the bacterium that attract flies also are utilized for swarming. In addition, bacteria with the swarming mutation rfaL attracted fewer blow flies and reduced the number of eggs laid by the flies. Therefore, we have identified several interkingdom signals between P. mirabilis and blow flies.

Global temperature changes have emphasized the need to understand how species adapt to thermal stress across their ranges. Genetic mechanisms may contribute to variation in thermal tolerance, providing evidence for how organisms adapt to local environments. We determine physiological thermal limits and characterize genome-wide transcriptional changes at these limits in bumble bees using laboratory-reared Bombus vosnesenskii workers. We analyze bees reared from latitudinal (35.7–45.7°N) and altitudinal (7–2154 m) extremes of the species’ range to correlate thermal tolerance and gene expression among populations from different climates. We find that critical thermal minima (CTMIN) exhibit strong associations with local minimums at the location of queen origin, while critical thermal maximum (CTMAX) was invariant among populations. Concordant patterns are apparent in gene expression data, with regional differentiation following cold exposure, and expression shifts invariant among populations under high temperatures. Furthermore, we identify several modules of co-expressed genes that tightly correlate with critical thermal limits and temperature at the region of origin. Our results reveal that local adaptation in thermal limits and gene expression may facilitate cold tolerance across a species range, whereas high temperature responses are likely constrained, both of which may have implications for climate change responses of bumble bees.

There can be substantial negative consequences for insects colonizing a resource in the presence of competitors. We hypothesized that bacteria, associated with an oviposition resource and the insect eggs deposited on that resource, serve as a mechanism regulating subsequent insect attraction, colonization, and potentially succession of insect species. We isolated and identified bacterial species associated with insects associated with vertebrate carrion and used these bacteria to measure their influence on the oviposition preference of adult black soldier flies which utilizes animal carcasses and is an important species in waste management and forensics. We also ascertained that utilizing a mixture of bacteria, rather than a single species, differentially influenced behavioral responses of the flies, as did bacterial concentration and the species of fly from which the bacteria originated. These studies provide insight into interkingdom interactions commonly occurring during decomposition, but not commonly studied.

The Calliphoridae or blow flies are a family of insects that occupy diverse habitats and perform important ecological roles, particularly the decomposition of animal remains. Some Calliphoridae species are also important in the forensic sciences, in agriculture (e.g. as livestock pests) and in medicine (e.g. maggot therapy). Calliphoridae provide striking examples in support of the hypothesis that sex determination regulatory gene hierarchies evolve in the reverse order, with the gene at the top being the most recently added. Unlike the model fly Drosophila melanogaster, where sex is determined by the number of X chromosomes, in the Australian sheep blow fly (Lucilia cuprina) sex is determined by a Y-linked male-determining gene (M). A different regulatory system appears to operate in the hairy maggot blow fly (Chrysomya rufifacies) where the maternal genotype determines sex. It is hypothesized that females heterozygous for a dominant female-determining factor (F/f) produce only female offspring and homozygous f/f females produce only sons. The bottom of the regulatory hierarchy appears to be the same in D. melanogaster and L. cuprina, with sex-specific splicing of doublesex transcripts being controlled by the female-specific Transformer (TRA) protein. We discuss a model that has been proposed for how tra transcripts are sex-specifically spliced in calliphorids, which is very different from D. melanogaster.

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.

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.