A rapid return to normal: temporal gene expression patterns following cold exposure in the bumble bee Bombus impatiens.

Bumble bees are common in cooler climates and many species likely experience periodic exposure to very cold temperatures, but little is known about the temporal dynamics of cold response mechanisms following chill exposure, especially how persistent effects of cold exposure may facilitate tolerance of future events. To investigate molecular processes involved in the temporal response by bumble bees to acute cold exposure, we compared mRNA transcript abundance in Bombus impatiens workers exposed to 0°C for 75 min (inducing chill coma) and control bees maintained at a constant ambient temperature (28°C). We sequenced the 3' end of mRNA transcripts (TagSeq) to quantify gene expression in thoracic tissue of bees at several time points (0, 10, 30, 120 and 720 min) following cold exposure. Significant differences from control bees were only detectable within 30 min after the treatment, with most occurring at the 10 min recovery time point. Genes associated with gluconeogenesis and glycolysis were most notably upregulated, while genes related to lipid and purine metabolism were downregulated. The observed patterns of expression indicate a rapid recovery after chill coma, suggesting an acute differential transcriptional response during recovery from chill coma and return to baseline expression levels within an hour, with no long-term gene expression markers of this cold exposure. Our work highlights the functions and pathways important for acute cold recovery, provides an estimated time frame for recovery from cold exposure in bumble bees, and suggests that cold hardening may be less important for these heterothermic insects.

Whole genome analyses reveal weak signatures of population structure and environmentally associated local adaptation in an important North American pollinator, the bumble bee Bombus vosnesenskii.

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 in Bombus 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 that B. vosnesenskii exhibits 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 indicate B. vosnesenskii is 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.

Community-wide correlations between species richness, abundance and population genomic diversity in a freshwater biodiversity hotspot.

Understanding patterns of diversity across macro (e.g. species-level) and micro (e.g. molecular-level) scales can shed light on community function and stability by elucidating the abiotic and biotic drivers of diversity within ecological communities. We examined the relationships among taxonomic and genetic metrics of diversity in freshwater mussels (Bivalvia: Unionidae), an ecologically important and species-rich group in the southeastern United States. Using quantitative community surveys and reduced-representation genome sequencing across 22 sites in seven rivers and two river basins, we surveyed 68 mussel species and sequenced 23 of these species to characterize intrapopulation genetic variation. We tested for the presence of species diversity-abundance correlations (i.e. the more-individuals hypothesis, MIH), species-genetic diversity correlations (SGDCs) and abundance-genetic diversity correlations (AGDCs) across all sites to evaluate relationships between different metrics of diversity. Sites with greater cumulative multispecies density (a standardized metric of abundance) had a greater number of species, consistent with the MIH hypothesis. Intrapopulation genetic diversity was strongly associated with the density of most species, indicating the presence of AGDCs. However, there was no consistent evidence for SGDCs. Although sites with greater overall densities of mussels had greater species richness, sites with higher genetic diversity did not always exhibit positive correlations with species richness, suggesting that there are spatial and evolutionary scales at which the processes influencing community-level diversity and intraspecific diversity differ. Our work reveals the importance of local abundance as indicator (and possibly a driver) of intrapopulation genetic diversity.

It is buzziness time: rearing, mating, and overwintering Bombus vosnesenskii (Hymenoptera: Apidae).

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.

Local adaptation across a complex bioclimatic landscape in two montane bumble bee species.

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.

Test of the invasive pathogen hypothesis of bumble bee decline in North America.

Emergent fungal diseases are critical factors in global biodiversity declines. The fungal pathogenNosema bombiwas recently found to be widespread in declining species of North American bumble bees (Bombus), with circumstantial evidence suggesting an exotic introduction from Europe. This interpretation has been hampered by a lack of knowledge of global genetic variation, geographic origin, and changing prevalence patterns ofN. bombiin declining North American populations. Thus, the temporal and spatial emergence ofN. bombiand its potential role in bumble bee decline remain speculative. We analyzeNosemaprevalence and genetic variation in the United States and Europe from 1980, before an alleged introduction in the early 1990s, to 2011, extractingNosemaDNA fromBombusnatural history collection specimens from across this time period.Nosema bombiprevalence increased significantly from low detectable frequency in the 1980s to significantly higher frequency in the mid- to late-1990s, corresponding to a period of reported massive infectious outbreak ofN. bombiin commercial bumble bee rearing stocks in North America. Despite the increased frequency, we find no conclusive evidence of an exoticN. bombiorigin based on genetic analysis of globalNosemapopulations; the widespreadNosemastrain found currently in declining United States bumble bees was present in the United States before commercial colony trade. Notably, the USN. bombiis not detectably different from that found predominantly throughout Western Europe, with both regions characterized by low genetic diversity compared with high levels of diversity found in Asia, where commercial bee breeding activities are low or nonexistent.

Distance, elevation and environment as drivers of diversity and divergence in bumble bees across latitude and altitude.

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.

Phylogenetic Analysis of Molluscan Metallothioneins: Evolutionary Insight from Crassostrea virginica.

Mechanisms by which organisms genetically adapt to environmental conditions are of fundamental importance to studies of evolutionary biology and environmental physiology. Natural selection acts on existing genetic variation leading to adaptation through selection of new mutations that confer beneficial advantages to populations. The American oyster, Crassostrea virginica, is an excellent model to investigate interactions between environmental and ecological factors as driving forces for natural selection. A great example of this is represented by the diversity of C. virginica metallothioneins (CvMT), metal-binding proteins involved in homeostasis and tolerance, that have resulted from a series of duplication events to produce the greatest structural diversity of MT proteins found in a single species. We present phylogenetic evidence of two distinct ancestral ß-domain MTs that gave rise to a variety of ßß and αß CvMT proteins, as well as CvMT-II proteins consisting solely of one to four α-domains. Furthermore, we annotate the complete locus containing the paralogous CvMT-I, -II, and -IV genes, providing supporting evidence of a hypothesized series of exon and gene duplication events that gave rise to the various CvMT-I and -II isoforms. We also highlight unique MT expression profiles from four separate C. virginica populations to demonstrate differences in gene diversity and copy number which appear to be enriched in southeastern U.S. oyster populations. These observations contribute to a better understanding of the molecular mechanisms leading to adaptation in organisms that experience substantial environmental stress, with a specific focus on evolutionary adaptations of gene structure.

Molecular tools and bumble bees: revealing hidden details of ecology and evolution in a model system.

Bumble bees are a longstanding model system for studies on behaviour, ecology and evolution, due to their well-studied social lifestyle, invaluable role as wild and managed pollinators, and ubiquity and diversity across temperate ecosystems. Yet despite their importance, many aspects of bumble bee biology have remained enigmatic until the rise of the genetic and, more recently, genomic eras. Here, we review and synthesize new insights into the ecology, evolution and behaviour of bumble bees that have been gained using modern genetic and genomic techniques. Special emphasis is placed on four areas of bumble bee biology: the evolution of eusociality in this group, population-level processes, large-scale evolutionary relationships and patterns, and immunity and resistance to pesticides. We close with a prospective on the future of bumble bee genomics research, as this rapidly advancing field has the potential to further revolutionize our understanding of bumble bees, particularly in regard to adaptation and resilience. Worldwide, many bumble bee populations are in decline. As such, throughout the review, connections are drawn between new molecular insights into bumble bees and our understanding of the causal factors involved in their decline. Ongoing and potential applications to bumble bee management and conservation are also included to demonstrate how genetics- and genomics-enabled research aids in the preservation of this threatened group.

Patterns of widespread decline in North American bumble bees.

Bumble bees (Bombus) are vitally important pollinators of wild plants and agricultural crops worldwide. Fragmentary observations, however, have suggested population declines in several North American species. Despite rising concern over these observations in the United States, highlighted in a recent National Academy of Sciences report, a national assessment of the geographic scope and possible causal factors of bumble bee decline is lacking. Here, we report results of a 3-y interdisciplinary study of changing distributions, population genetic structure, and levels of pathogen infection in bumble bee populations across the United States. We compare current and historical distributions of eight species, compiling a database of >73,000 museum records for comparison with data from intensive nationwide surveys of >16,000 specimens. We show that the relative abundances of four species have declined by up to 96% and that their surveyed geographic ranges have contracted by 23-87%, some within the last 20 y. We also show that declining populations have significantly higher infection levels of the microsporidian pathogen Nosema bombi and lower genetic diversity compared with co-occurring populations of the stable (nondeclining) species. Higher pathogen prevalence and reduced genetic diversity are, thus, realistic predictors of these alarming patterns of decline in North America, although cause and effect remain uncertain.

Genome-wide DNA methylation patterns in bumble bee (Bombus vosnesenskii) populations from spatial-environmental range extremes.

Unraveling molecular mechanisms of adaptation to complex environments is crucial to understanding tolerance of abiotic pressures and responses to climatic change. Epigenetic variation is increasingly recognized as a mechanism that can facilitate rapid responses to changing environmental cues. To investigate variation in genetic and epigenetic diversity at spatial and thermal extremes, we use whole genome and methylome sequencing to generate a high-resolution map of DNA methylation in the bumble bee Bombus vosnesenskii. We sample two populations representing spatial and environmental range extremes (a warm southern low-elevation site and a cold northern high-elevation site) previously shown to exhibit differences in thermal tolerance and determine positions in the genome that are consistently and variably methylated across samples. Bisulfite sequencing reveals methylation characteristics similar to other arthropods, with low global CpG methylation but high methylation concentrated in gene bodies and in genome regions with low nucleotide diversity. Differentially methylated sites (n = 2066) were largely hypomethylated in the northern high-elevation population but not related to local sequence differentiation. The concentration of methylated and differentially methylated sites in exons and putative promoter regions suggests a possible role in gene regulation, and this high-resolution analysis of intraspecific epigenetic variation in wild Bombus suggests that the function of methylation in niche adaptation would be worth further investigation.

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.

Metabolomes of bumble bees reared in common garden conditions suggest constitutive differences in energy and toxin metabolism across populations.

Cold tolerance of ectotherms can vary strikingly among species and populations. Variation in cold tolerance can reflect differences in genomes and transcriptomes that confer cellular-level protection from cold; additionally, shifts in protein function and abundance can be altered by other cellular constituents as cold-exposed insects often have shifts in their metabolomes. Even without a cold challenge, insects from different populations may vary in cellular composition that could alter cold tolerance, but investigations of constitutive differences in metabolomes across wild populations remain rare. To address this gap, we reared Bombus vosnesenskii queens collected from Oregon and California (USA) that differ in cold tolerance (CTmin = -6 °C and 0 °C, respectively) in common garden conditions, and measured offspring metabolomes using untargeted LC-MS/MS. Oregon bees had higher levels of metabolites associated with carbohydrate (sorbitol, lactitol, maltitol, and sorbitol-6-phosphate) and amino acid (hydroxyproline, ornithine, and histamine) metabolism. Exogenous metabolites, likely derived from the diet, also varied between Oregon and California bees, suggesting population-level differences in toxin metabolism. Overall, our results reveal constitutive differences in metabolomes for bumble bees reared in common garden conditions from queens collected in different locations despite no previous cold exposure.

A trait dataset for freshwater mussels of the United States of America.

The United States of America has a diverse collection of freshwater mussels comprising 301 species distributed among 59 genera and two families (Margaritiferidae and Unionidae), each having a unique suite of traits. Mussels are among the most imperilled animals and are critical components of their ecosystems, and successful management, conservation and research requires a cohesive and widely accessible data source. Although trait-based analysis for mussels has increased, only a small proportion of traits reflecting mussel diversity in this region has been collated. Decentralized and non-standardized trait information impedes large-scale analysis. Assembling trait data in a synthetic dataset enables comparison across species and lineages and identification of data gaps. We collated data from the primary literature, books, state and federal reports, theses and dissertations, and museum collections into a centralized dataset covering information on taxonomy, morphology, reproductive ecology and life history, fish hosts, habitats, thermal tolerance, geographic distribution, available genetic information, and conservation status. By collating these traits, we aid researchers in assessing variation in mussel traits and modelling ecosystem change.

Geographical patterns of genetic divergence in the widespread Mesoamerican bumble bee Bombus ephippiatus (Hymenoptera: Apidae).

Bumble bees (Bombus Latreille) are an important group of social insects, well recognized throughout northern temperate regions as important pollinators of wild and agricultural plants. Little is known about the biology of this group in southern portions of the Americas, especially in Mesoamerica, a region of geological and ecological complexity from Mexico through Central America. One ubiquitous Mesoamerican species, Bombus ephippiatus, is enigmatic. Like many other Bombus, this species is homogeneous in body structure yet exhibits striking intraspecific color pattern polymorphism across its range, leading to uncertainty about its genealogical boundaries. It has been grouped taxonomically with B. wilmattae, a species narrowly restricted to southern Mexico and northern Guatamala. Furthermore, the relationships between these two taxa and a third species, B. impatiens, found only in America north of Mexico, have been controversial. Our phylogenetic analysis of DNA sequences from mitochondrial COI and nuclear PEPCK and CAD resolves the phylogeny of these three taxa as (B. impatiens, (B. ephippiatus, B. wilmattae)). Additional data from eight nuclear microsatellite markers reveal complex patterns of genetic divergence and isolation among populations of B. ephippiatus across its extensive geographic range, providing evidence for multiple independent evolutionary lineages. These lineages correspond not only to geographic and habitat variation across their range, but also to distinct color pattern groups present in the species. Knowledge of the phylogeny and genetic divergence of the B. ephippiatus group will provide a framework for understanding evolutionary and ecological origins of color pattern polymorphism in bumble bees, as well as providing insight into geographical factors enhancing speciation in Mesoamerica.

Interspecific geographic distribution and variation of the pathogens Nosema bombi and Crithidia species in United States bumble bee populations.

Several bumble bee (Bombus) species in North America have undergone range reductions and rapid declines in relative abundance. Pathogens have been suggested as causal factors, however, baseline data on pathogen distributions in a large number of bumble bee species have not been available to test this hypothesis. In a nationwide survey of the US, nearly 10,000 specimens of 36 bumble bee species collected at 284 sites were evaluated for the presence and prevalence of two known Bombus pathogens, the microsporidium Nosema bombi and trypanosomes in the genus Crithidia. Prevalence of Crithidia was ≤10% for all host species examined but was recorded from 21% of surveyed sites. Crithidia was isolated from 15 of the 36 Bombus species screened, and were most commonly recovered from Bombus bifarius, Bombus bimaculatus, Bombus impatiens and Bombus mixtus. Nosema bombi was isolated from 22 of the 36 US Bombus species collected. Only one species with more than 50 sampled bees, Bombus appositus, was free of the pathogen; whereas, prevalence was highest in Bombus occidentalis and Bombus pensylvanicus, two species that are reportedly undergoing population declines in North America. A variant of a tetranucleotide repeat in the internal transcribed spacer (ITS) of the N. bombi rRNA gene, thus far not reported from European isolates, was isolated from ten US Bombus hosts, appearing in varying ratios in different host species. Given the genetic similarity of the rRNA gene of N. bombi sampled in Europe and North America to date, the presence of a unique isolate in US bumble could reveal one or more native North American strains and indicate that N. bombi is enzootic across the Holarctic Region, exhibiting some genetic isolation.

RAD-tag and mitochondrial DNA sequencing reveal the genetic structure of a widespread and regionally imperiled freshwater mussel, Obovaria olivaria (Bivalvia: Unionidae).

Obovaria olivaria is a species of freshwater mussel native to the Mississippi River and Laurentian Great Lakes-St. Lawrence River drainages of North America. This mussel has experienced population declines across large parts of its distribution and is imperiled in many jurisdictions. Obovaria olivaria uses the similarly imperiled Acipenser fulvescens (Lake Sturgeon) as a host for its glochidia. We employed mitochondrial DNA sequencing and restriction site-associated DNA sequencing (RAD-seq) to assess patterns of genetic diversity and population structure of O. olivaria from 19 collection locations including the St. Lawrence River drainage, the Great Lakes drainage, the Upper Mississippi River drainage, the Ohioan River drainage, and the Mississippi Embayment. Heterozygosity was highest in Upper Mississippi and Great Lakes populations, followed by a reduction in diversity and relative effective population size in the St. Lawrence populations. Pairwise F ST ranged from 0.00 to 0.20, and analyses of genetic structure revealed two major ancestral populations, one including all St. Lawrence River/Ottawa River sites and the other including remaining sites; however, significant admixture and isolation by river distance across the range were evident. The genetic diversity and structure of O. olivaria is consistent with the existing literature on Acipenser fulvescens and suggests that, although northern and southern O. olivaria populations are genetically distinct, genetic structure in O. olivaria is largely clinal rather than discrete across its range. Conservation and restoration efforts of O. olivaria should prioritize the maintenance and restoration of locations where O. olivaria remain, especially in northern rivers, and to ensure connectivity that will facilitate dispersal of Acipenser fulvescens and movement of encysted glochidia.

Environment and Co-occurring Native Mussel Species, but Not Host Genetics, Impact the Microbiome of a Freshwater Invasive Species (Corbicula fluminea).

The Asian clam Corbicula fluminea (Family: Cyneridae) has aggressively invaded freshwater habitats worldwide, resulting in dramatic ecological changes and declines of native bivalves such as freshwater mussels (Family: Unionidae), one of the most imperiled faunal groups. Despite increases in our knowledge of invasive C. fluminea biology, little is known of how intrinsic and extrinsic factors, including co-occurring native species, influence its microbiome. We investigated the gut bacterial microbiome across genetically differentiated populations of C. fluminea in the Tennessee and Mobile River Basins in the Southeastern United States and compared them to those of six co-occurring species of native freshwater mussels. The gut microbiome of C. fluminea was diverse, differed with environmental conditions and varied spatially among rivers, but was unrelated to host genetic variation. Microbial source tracking suggested that the gut microbiome of C. fluminea may be influenced by the presence of co-occurring native mussels. Inferred functions from 16S rRNA gene data using PICRUST2 predicted a high prevalence and diversity of degradation functions in the C. fluminea microbiome, especially the degradation of carbohydrates and aromatic compounds. Such modularity and functional diversity of the microbiome of C. fluminea may be an asset, allowing to acclimate to an extensive range of nutritional sources in invaded habitats, which could play a vital role in its invasive success.

Patterns of range-wide genetic variation in six North American bumble bee (Apidae: Bombus) species.

The increasing evidence for population declines in bumble bee (Bombus) species worldwide has accelerated research efforts to explain losses in these important pollinators. In North America, a number of once widespread Bombus species have suffered serious reductions in range and abundance, although other species remain healthy. To examine whether declining and stable species exhibit different levels of genetic diversity or population fragmentation, we used microsatellite markers to genotype populations sampled across the geographic distributions of two declining (Bombus occidentalis and Bombus pensylvanicus) and four stable (Bombus bifarius; Bombus vosnesenskii; Bombus impatiens and Bombus bimaculatus) Bombus species. Populations of declining species generally have reduced levels of genetic diversity throughout their range compared to codistributed stable species. Genetic diversity can be affected by overall range size and degree of isolation of local populations, potentially confounding comparisons among species in some cases. We find no evidence for consistent differences in gene flow among stable and declining species, with all species exhibiting weak genetic differentiation over large distances (e.g. >1000 km). Populations on islands and at high elevations experience relatively strong genetic drift, suggesting that some conditions lead to genetic isolation in otherwise weakly differentiated species. B. occidentalis and B. bifarius exhibit stronger genetic differentiation than the other species, indicating greater phylogeographic structure consistent with their broader geographic distributions across topographically complex regions of western North America. Screening genetic diversity in North American Bombus should prove useful for identifying species that warrant monitoring, and developing management strategies that promote high levels of gene flow will be a key component in efforts to maintain healthy populations.