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.

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.

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.

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 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.

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.

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.

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.

Body mass and sex, not local climate, drive differences in chill coma recovery times in common garden reared bumble bees.

The time required to recover from cold exposure (chill coma recovery time) may represent an important metric of performance and has been linked to geographic distributions of diverse species. Chill coma recovery time (CCRT) has rarely been measured in bumble bees (genus Bombus) but may provide insights regarding recent changes in their distributions. We measured CCRT of Bombus vosnesenskii workers reared in common garden laboratory conditions from queens collected across altitude and latitude in the Western United States. We also compared CCRTs of male and female bumble bees because males are often overlooked in studies of bumble bee ecology and physiology and may differ in their ability to respond to cold temperatures. We found no relationship between CCRT and local climate at the queen collection sites, but CCRT varied significantly with sex and body mass. Because differences in the ability to recover from cold temperatures have been shown in wild-caught Bombus, we predict that variability in CCRT may be strongly influenced by plasticity.

Biogeographic parallels in thermal tolerance and gene expression variation under temperature stress in a widespread bumble bee.

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.

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.

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.

Adaptation to the abiotic environment in insects: the influence of variability on ecophysiology and evolutionary genomics.

Advances in tools to gather environmental, phenotypic, and molecular data have accelerated our ability to detect abiotic drivers of variation across the genome-to-phenome spectrum in model and non-model insects. However, differences in the spatial and temporal resolution of these data sets may create gaps in our understanding of linkages between environment, genotype, and phenotype that yield missed or misleading results about adaptive variation. In this review we highlight sources of variability that might impact studies of phenotypic and 'omic environmental adaptation, challenges to collecting data at relevant scales, and possible solutions that link intensive fine-scale reductionist studies of mechanisms to large-scale biogeographic patterns.

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.

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.

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.

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.