Search

Ecological theory predicts that species turnover among communities (e.g. β-diversity) and genetic turnover among populations within species (e.g. F ST ) should be positively correlated if similar processes influence colonization and occupancy of species and gene flow and genetic drift of populations within a metacommunity. Using recently published population genomic data from multiple populations of 15 freshwater mussel (Unionidae) species across seven rivers in the Mobile and Tennessee River basins of the south-eastern USA, we conducted novel analyses examining the relationship between taxonomic turnover (β-diversity) among communities and genetic differentiation ( F ST ) within these species. F ST and β-diversity were both hierarchically structured, and strong basin effects and isolation-by-distance were observed for β-diversity and for F ST among populations within most species. Furthermore, β-diversity and F ST were directly correlated for the overall community and among sites for individual species, indicating that factors shaping turnover among mussel assemblages are similar at the species and genetic levels. The widespread associations between turnover metrics at the community and population genetic levels of biological organization suggest that parallel processes govern species composition and intraspecific connectivity in freshwater mussel metacommunities.

The gut microbiome is influenced by host species and the environment, but how the environment influences the microbiome of animals introduced into a new ecosystem has rarely been investigated. Freshwater mussels are aquatic fauna, with some threatened or endangered species propagated in hatcheries and introduced into natural systems as part of conservation efforts. The effects of the environment on the freshwater mussel gut microbiome were assessed for two hatchery-propagated species ( Lampsilis ovata , Lampsilis ornata ) introduced into rivers within their natural range. Mussels were placed in rivers for 8 weeks, after which one subset was collected, another subset remained in that river, and a third subset was reciprocally transplanted to another river in the same river basin for a further 8 weeks. Gut microbiome composition and diversity were characterized for all mussels. After the initial 8 weeks, mussels showed increased gut bacterial species richness and distinct community composition compared to hatchery mussels, but gut microbiome diversity then decreased for mussels that remained in the same river for all 16 weeks. The gut bacterial community of mussels transplanted between rivers shifted to resemble that of mussels placed initially into the recipient river and that remained there for the whole study. All mussels showed high proportions of Firmicutes in their gut microbiome after 8 weeks, suggesting an essential role of this phylum in the gut of Lampsilis species. These findings show that the mussel gut microbiome shifts in response to new environments and provide insights into conservation strategies that involve species reintroductions., Competing Interests: The authors report that they have no conflicts of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of the International Society for Microbial Ecology.)

Understanding the myriad avenues through which spatial and environmental factors shape evolution is a major focus in biological research. From a molecular perspective, much work has been focused on genomic sequence variation; however, recently there has been increased interest in how epigenetic variation may be shaped by different variables across the landscape. DNA methylation has been of particular interest given that it is dynamic and can alter gene expression, potentially offering a path for a rapid response to environmental change. We utilized whole genome enzymatic methyl sequencing to evaluate the distribution of CpG methylation across the genome and to analyze patterns of spatial and environmental association in the methylomes of two broadly distributed montane bumble bees (Bombus vancouverensis Cresson and Bombus vosnesenskii Radoszkowski) across elevational gradients in the western US. Methylation patterns in both species are similar at the genomic scale with ∼1% of CpGs being methylated and most methylation being found in exons. At the landscape scale, neither species exhibited strong spatial or population structuring in patterns of methylation, although some weak relationships between methylation and distance or environmental variables were detected. Differential methylation analysis suggests a stronger environment association in B. vancouverensis given the larger number of differentially methylated CpG's compared to B. vosnesenskii. We also observed only a handful of genes with both differentially methylated CpGs and previously detected environmentally associated outlier SNPs. Overall results reveal a weak but present pattern in variation in methylation over the landscape in both species., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

The Hunt bumble bee, Bombus huntii, is a widely distributed pollinator in western North America. The species produces large colony sizes in captive rearing conditions, experiences low parasite and pathogen loads, and has been demonstrated to be an effective pollinator of tomatoes grown in controlled environment agriculture systems. These desirable traits have galvanized producer efforts to develop commercial Bombus huntii colonies for growers to deliver pollination services to crops. To better understand Bombus huntii biology and support population genetic studies and breeding decisions, we sequenced and assembled the Bombus huntii genome from a single haploid male. High-fidelity sequencing of the entire genome using PacBio, along with HiC sequencing, led to a comprehensive contig assembly of high continuity. This assembly was further organized into a chromosomal arrangement, successfully identifying 18 chromosomes spread across the 317.4 Mb assembly with a BUSCO score indicating 97.6% completeness. Synteny analysis demonstrates shared chromosome number (n = 18) with Bombus terrestris, a species belonging to a different subgenus, matching the expectation that presence of 18 haploid chromosomes is an ancestral trait at least between the subgenera Pyrobombus and Bombus sensu stricto. In conclusion, the assembly outcome, alongside the minimal tissue sampled destructively, showcases efficient techniques for producing a comprehensive, highly contiguous genome., Competing Interests: Conflicts of interests The authors declare no conflicts of interest., (Published by Oxford University Press on behalf of The Genetics Society of America 2024.)

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., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2024. Published by The Company of Biologists Ltd.)

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 (CT min  = -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., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

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., (© 2023 John Wiley & Sons Ltd.)

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., (© 2023. The Author(s).)

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., (© 2023 John Wiley & Sons Ltd.)

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., (© 2023. Springer Nature Limited.)

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., (© The Author(s) 2023. Published by Oxford University Press on behalf of Entomological Society of America.)

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., Competing Interests: Neither I nor my co‐authors have a conflict of interest. No materials should require permission to reproduce., (© 2023 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Chiarello, Bucholz, McCauley, Vaughn, Hopper, Sánchez González, Atkinson, Lozier and Jackson.)

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., Competing Interests: None declared., (© 2022 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd.)

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., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

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 (CT MIN ) exhibit strong associations with local minimums at the location of queen origin, while critical thermal maximum (CT MAX ) 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.

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., (Copyright © 2020 Heraghty et al.)

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., (© 2020 John Wiley & Sons Ltd.)

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., (Copyright © 2019 Elsevier Inc. All rights reserved.)

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., (© 2018 John Wiley & Sons Ltd.)

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.

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.

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.

Geographic variation in insect coloration is among the most intriguing examples of rapid phenotypic evolution and provides opportunities to study mechanisms of phenotypic change and diversification in closely related lineages. The bumble bee Bombus bifarius comprises two geographically disparate color groups characterized by red-banded and black-banded abdominal pigmentation, but with a range of spatially and phenotypically intermediate populations across western North America. Microsatellite analyses have revealed that B. bifarius in the USA are structured into two major groups concordant with geography and color pattern, but also suggest ongoing gene flow among regional populations. In this study, we better resolve the relationships among major color groups to better understand evolutionary mechanisms promoting and maintaining such polymorphism. We analyze >90,000 and >25,000 single-nucleotide polymorphisms derived from transcriptome (RNAseq) and double digest restriction site associated DNA sequencing (ddRAD), respectively, in representative samples from spatial and color pattern extremes in B. bifarius as well as phenotypic and geographic intermediates. Both ddRAD and RNAseq data illustrate substantial genome-wide differentiation of the red-banded (eastern) color form from both black-banded (western) and intermediate (central) phenotypes and negligible differentiation among the latter populations, with no obvious admixture among bees from the two major lineages. Results thus indicate much stronger background differentiation among B. bifarius lineages than expected, highlighting potential challenges for revealing loci underlying color polymorphism from population genetic data alone. These findings will have significance for resolving taxonomic confusion in this species and in future efforts to investigate color-pattern evolution in B. bifarius and other polymorphic bumble bee species.

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., (© 2015 John Wiley & Sons Ltd.)

Background: The shift from solitary to social behavior is one of the major evolutionary transitions. Primitively eusocial bumblebees are uniquely placed to illuminate the evolution of highly eusocial insect societies. Bumblebees are also invaluable natural and agricultural pollinators, and there is widespread concern over recent population declines in some species. High-quality genomic data will inform key aspects of bumblebee biology, including susceptibility to implicated population viability threats., Results: We report the high quality draft genome sequences of Bombus terrestris and Bombus impatiens, two ecologically dominant bumblebees and widely utilized study species. Comparing these new genomes to those of the highly eusocial honeybee Apis mellifera and other Hymenoptera, we identify deeply conserved similarities, as well as novelties key to the biology of these organisms. Some honeybee genome features thought to underpin advanced eusociality are also present in bumblebees, indicating an earlier evolution in the bee lineage. Xenobiotic detoxification and immune genes are similarly depauperate in bumblebees and honeybees, and multiple categories of genes linked to social organization, including development and behavior, show high conservation. Key differences identified include a bias in bumblebee chemoreception towards gustation from olfaction, and striking differences in microRNAs, potentially responsible for gene regulation underlying social and other traits., Conclusions: These two bumblebee genomes provide a foundation for post-genomic research on these key pollinators and insect societies. Overall, gene repertoires suggest that the route to advanced eusociality in bees was mediated by many small changes in many genes and processes, and not by notable expansion or depauperation.

Genetic variation is of key importance for a species' evolutionary potential, and its estimation is a major component of conservation studies. New DNA sequencing technologies have enabled the analysis of large portions of the genome in nonmodel species, promising highly accurate estimates of such population genetic parameters. Restriction site-associated DNA sequencing (RADseq) is used to analyse thousands of variants in the bumble bee species Bombus impatiens, which is common, and Bombus pensylvanicus, which is in decline. Previous microsatellite-based analyses have shown that gene diversity is lower in the declining B. pensylvanicus than in B. impatiens. RADseq nucleotide diversities appear much more similar in the two species. Both species exhibit allele frequencies consistent with historical population expansions. Differences in diversity observed at microsatellites thus do not appear to have arisen from long-term differences in population size and are either recent in origin or may result from mutational processes. Additional research is needed to explain these discrepancies and to investigate the best ways to integrate next-generation sequencing data and more traditional molecular markers in studies of genetic diversity., (© 2013 John Wiley & Sons Ltd.)

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., (Copyright © 2012 Elsevier Inc. All rights reserved.)

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., (Copyright © 2011 Elsevier Inc. All rights reserved.)

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., (© 2011 Blackwell Publishing Ltd.)

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.

Background: The integration of multiple complementary approaches is a powerful way to understand the processes of diversification and speciation. The parasitoid wasp Aphidius transcaspicus Telenga (Hymenoptera: Braconidae) is a parasitoid of Hyalopterus aphids across a wide geographic range. This species shows a remarkable degree of genetic structure among western, central, and eastern Mediterranean population clusters. In this paper we attempt to better characterize this genetic structure., Methodology/principal Findings: We use a Bayesian coalescent analysis of gene flow under the Isolation with Migration model using mitochondrial and microsatellite markers together with climate-based ecological niche models to better understand the genetic structure of A. transcaspicus in the Mediterranean. The coalescent analysis revealed low levels of migration among western and eastern Mediterranean populations (Nm<1) that were not statistically distinguishable from zero. Niche models showed that localities within population clusters each occupy areas of continuously high environmental suitability, but are separated from each other by large regions of completely unsuitable habitat that could limit dispersal. Overall, environmental characteristics were similar among the population clusters, though significant differences did emerge., Conclusions/significance: These results support contemporary allopatric isolation of Mediterranean populations of A. transcaspicus, which together with previous analyses indicating partial behaviorally mediated reproductive isolation, suggest that the early stages of cryptic speciation may be in progress.

Direct comparison of genetic patterns between museum specimens and contemporary collections can be a powerful approach for detecting recent demographic changes. Using microsatellite markers, we examined historical and contemporary genetic variation from an apparently declining bumble bee species, Bombus pensylvanicus, and from a stable species, Bombus impatiens, in central Illinois. For each species, we genotyped specimens from the Illinois Natural History Survey collected from three populations between 1969-1972 and from a resurvey of the same areas conducted in 2008. Population structure in B. pensylvanicus increased markedly over the last four decades (from theta(ST) = 0.001 to 0.027) while no structure was detected in B. impatiens for either time period (theta(ST) = -0.006 to -0.003). Changes in genetic diversity were not significant for either species, although small reductions were observed for B. pensylvanicus in all three populations. Coalescent simulations incorporating both contemporary and historical samples suggest that this small change is not surprising for recent population declines, as large reductions in genetic diversity were only apparent under the most severe bottleneck scenarios. These results demonstrate how comparisons of genetic patterns between temporal periods and species can help elucidate potential threats to population health and suggest several strategies that might be useful in the conservation of B. pensylvanicus in the Midwestern USA.

Host plant associated genetic differentiation is a common phenomenon in phytophagous insects, but the degree to which such associations sequentially drive diversification at higher trophic levels is not as well understood. A recent study examining neutral molecular markers in Hyalopterus aphids revealed that genetic structure in this genus is strongly determined by primary host plant use (Prunus spp.). In this paper, we take a similar approach to determine whether this host plant specificity has affected genetic structure in the parasitoid Aphidius transcaspicus, an important natural enemy of Hyalopterus spp. in the Mediterranean. Mitochondrial DNA (428 bp) and seven microsatellite loci were examined in parasitoids collected from aphid populations on almond, apricot, peach and plum trees from Spain and Greece. In contrast to the previous findings for Hyalopterus from the same regions, here we find no evidence for host associated diversification in A. transcaspicus at the species level or below, though geographic structure between regional populations is exceptionally high. These findings have several implications for our understanding of the ecology and evolution of A. transcaspicus as well as for its use as a biological control agent for Hyalopterus, suggesting that a consideration of host plant specificity may be less critical than factors such as climatic suitability or geographic origins of invasive populations.

Over the past several decades biologists' fascination with plant-herbivore interactions has generated intensive research into the implications of these interactions for insect diversification. The study of closely related phytophagous insect species or populations from an evolutionary perspective can help illuminate ecological and selective forces that drive these interactions. Here we present such an analysis for aphids in the genus Hyalopterus (Hemiptera: Aphididae), a cosmopolitan group that feeds on plants in the genus Prunus (Rosaceae). Hyalopterus currently contains two recognized species associated with different Prunus species, although the taxonomy and evolutionary history of the group is poorly understood. Using mitochondrial COI sequences, 16S rDNA sequences from the aphid endosymbiont Buchnera aphidicola, and nine microsatellite loci we investigated population structure in Hyalopterus from the most commonly used Prunus host species throughout the Mediterranean as well as in California, where the species H. pruni is an invasive pest. We found three deeply divergent lineages structured in large part by specific associations with plum, almond, and peach trees. There was no evidence that geographic or temporal barriers could explain the overall diversity in the genus. Levels of genetic differentiation are consistent with that typically attributed to aphid species and indicate divergence times older than the domestication of Prunus for agriculture. Interestingly, in addition to their typical hosts, aphids from each of the three lineages were frequently found on apricot trees. Apricot also appears to act as a resource mediated hybrid zone for plum and almond associated lineages. Together, results suggest that host plants have played a role in maintaining host-associated differentiation in Hyalopterus for as long as several million years, despite worldwide movement of host plants and the potential for ongoing hybridization.