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

A role for DNA methylation in bumblebee morphogenesis hints at female-specific developmental erasure

Epigenetic mechanisms, such as DNA methylation, are crucial factors in animal development. In some mammals, almost all DNA methylation is erased during embryo development and re-established in a sex- and cell-specific manner. This erasure and re-establishment is thought to primarily be a vertebrate-specific trait. Insects are particularly interesting in terms of development as many species often undergo remarkable morphological changes en route to maturity, that is, morphogenesis. However, little is known about the role of epigenetic mechanisms in this process across species. We have used whole-genome bisulfite sequencing to track genome-wide DNA methylation changes through the development of an economically and environmentally important pollinator species, the bumblebee Bombus terrestris (Hymenoptera:Apidae Linnaeus). We find overall levels of DNA methylation vary throughout development, and we find developmentally relevant differentially methylated genes throughout. Intriguingly, we have identified a depletion of DNA methylation in ovaries/eggs and an enrichment of highly methylated genes in sperm. We suggest this could represent a sex-specific DNA methylation erasure event. To our knowledge, this is the first suggestion of possible developmental DNA methylation erasure in an insect species. This study lays the required groundwork for functional experimental work to determine if there is a causal nature to the DNA methylation differences identified. Additionally, the application of single-cell methylation sequencing to this system will enable more accurate identification of if or when DNA methylation is erased during development.

Bumble Bees (Bombus terrestris) Use Time-Memory to Associate Reward with Color and Time of Day

Circadian clocks regulate ecologically important complex behaviors in honey bees, but it is not clear whether similar capacities exist in other species of bees. One key behavior influenced by circadian clocks is time-memory, which enables foraging bees to precisely time flower visitation to periods of maximal pollen or nectar availability and reduces the costs of visiting a non-rewarding flower patch. Bumble bees live in smaller societies and typically forage over shorter distances than honey bees, and it is therefore not clear whether they can similarly associate reward with time of day. We trained individually marked bumble bee (Bombus terrestris) workers to forage for sugar syrup in a flight cage with yellow or blue feeders rewarding either during the morning or evening. After training for over two weeks, we recorded all visitations to colored feeders filled with only water. We performed two experiments, each with a different colony. We found that bees tended to show higher foraging activity during the morning and evening training sessions compared to other times during the day. During the test day, the trained bees were more likely to visit the rewarding rather than the non-rewarding colored feeders at the same time of day during the test sessions, indicating that they associated time of day and color with the sugar syrup reward. These observations lend credence to the hypothesis that bumble bees have efficient time-memory, indicating that this complex behavior is not limited to honey bees that evolved sophisticated social foraging behaviors over large distances.

DNA methylation is associated with codon degeneracy in a species of bumblebee

Social insects display extreme phenotypic differences between sexes and castes even though the underlying genome can be almost identical. Epigenetic processes have been proposed as a possible mechanism for mediating these phenotypic differences. Using whole genome bisulfite sequencing of queens, males, and reproductive female workers we have characterised the sex- and caste-specific methylome of the bumblebee Bombus terrestris. We have identified a potential role for DNA methylation in histone modification processes which may influence sex and caste phenotypic differences. We also find differentially methylated genes generally show low levels of DNA methylation which may suggest a separate function for lowly methylated genes in mediating transcriptional plasticity, unlike highly methylated genes which are usually involved in housekeeping functions. We also examined the relationship between the underlying genome and the methylome using whole genome re-sequencing of the same queens and males. We find DNA methylation is enriched at zero-fold degenerate sites. We suggest DNA methylation may be acting as a targeted mutagen at these sites, providing substrate for selection via non-synonymous changes in the underlying genome. However, we did not see any relationship between DNA methylation and rates of positive selection in our samples. In order to fully assess a possible role for DNA methylation in adaptive processes a specifically designed study using natural population data is needed.

Infection Prevalence of Microsporidia Vairimorpha (Nosema) spp. in Japanese Bumblebees

Microsporidia are spore-forming intracellular parasites of various invertebrates and vertebrates. Vairimorpha bombi negatively affects the fitness of bumblebees and its prevalence correlates with declining bumblebee populations. The invasive alien species Bombus terrestris colonized Japan and possibly introduced new parasites. To assess the infection prevalence of V. bombi in Japanese bumblebees and B. terrestris, we investigated V. bombi infections using PCR and microscopy. The prevalence of sporulating V. bombi infections in three Bombus s. str. species/subspecies was low, whereas that of non/low-sporulating Vairimorpha sp. infections in three Diversobombus species/subspecies was high. Invasive B. terrestris showed low prevalence of non/low-sporulating V. bombi infections and shared the same V. bombi haplotype with B. hypocrita found in Hokkaido, where B. terrestris is present, and in Honshu, where B. terrestris is absent. Although V. bombi may have been introduced with B. terrestris colonies imported from Europe, it seems to be originally distributed in Japan. Furthermore, a new Vairimorpha sp. was found in Japanese bumblebee species. V. bombi and Vairimorpha sp. showed different organ and host specificities in bumblebees. There are no reports on the specific effects of other Vairimorpha spp. on bumblebees; further studies are needed to clarify the individual characteristics of Vairimorpha spp.

Using physiology to better support wild bee conservation

There is accumulating evidence that wild bees are experiencing a decline in terms of species diversity, abundance or distribution, which leads to major concerns about the sustainability of both pollination services and intrinsic biodiversity. There is therefore an urgent need to better understand the drivers of their decline, as well as design conservation strategies. In this context, the current approach consists of linking observed occurrence and distribution data of species to environmental features. While useful, a highly complementary approach would be the use of new biological metrics that can link individual bee responses to environmental alteration with population-level responses, which could communicate the actual bee sensitivity to environmental changes and act as early warning signals of bee population decline or sustainability. We discuss here through several examples how the measurement of bee physiological traits or performance can play this role not only in better assessing the impact of anthropogenic pressures on bees, but also in guiding conservation practices with the help of the documentation of species' physiological needs. Last but not least, because physiological changes generally occur well in advance of demographic changes, we argue that physiological traits can help in predicting and anticipating future population trends, which would represent a more proactive approach to conservation. In conclusion, we believe that future efforts to combine physiological, ecological and population-level knowledge will provide meaningful contributions to wild bee conservation-based research.

Anthropogenic effects on the body size of two neotropical orchid bees

To accommodate an ever-increasing human population, agriculture is rapidly intensifying at the expense of natural habitat, with negative and widely reported effects on biodiversity in general and on wild bee abundance and diversity in particular. Cities are similarly increasing in area, though the impact of urbanisation on wild bees is more equivocal and potentially positive in northern temperate regions. Yet agriculture and urbanisation both lead to the loss and alteration of natural habitat, its fragmentation, a potential reduction in floral availability, and warmer temperatures, factors thought to be drivers of wild bee decline. They have also been shown to be factors to which wild bee populations respond through morphological change. Body size is one such trait that, because of its relation to individual fitness, has received growing attention as a morphological feature that responds to human induced modification in land use. Here, we investigated the change in body size of two sympatric orchid bee species on the Yucatan Peninsula of Mexico in response to urbanization and agricultural intensification. By measuring 540 male individuals sampled from overall 24 sites, we found that Euglossa dilemma and Euglossa viridissima were on average smaller in urban and agricultural habitats than in natural ones. We discuss the potential role of reduced availability of resources in driving the observed body size shifts. Agricultural and urban land management in tropical regions might benefit wild bees if it encompassed the planting of flowering herbs and trees to enhance their conservation.

Pollinator movement activity influences genetic diversity and differentiation of spatially isolated populations of clonal forest herbs

In agricultural landscapes, forest herbs live in small, spatially isolated forest patches. For their long-term survival, their populations depend on animals as genetic linkers that provide pollen- or seed-mediated gene flow among different forest patches. However, whether insect pollinators serve as genetic linkers among spatially isolated forest herb populations in agricultural landscapes remains to be shown. Here, we used population genetic methods to analyze: (A) the genetic diversity and genetic differentiation of populations of two common, slow-colonizing temperate forest herb species [Polygonatum multiflorum (L.) All. and Anemone nemorosa L.] in spatially isolated populations within three agricultural landscapes in Germany and Sweden and (B) the movement activity of their most relevant associated pollinator species, i.e., the bumblebee Bombus pascuorum (Scopoli, 1,763) and the hoverfly Melanostoma scalare (Fabricus, 1,794), respectively, which differ in their mobility. We tested whether the indicated pollinator movement activity affected the genetic diversity and genetic differentiation of the forest herb populations. Bumblebee movement indicators that solely indicated movement activity between the forest patches affected both genetic diversity and genetic differentiation of the associated forest herb P. multiflorum in a way that can be explained by pollen-mediated gene flow among the forest herb populations. In contrast, movement indicators reflecting the total movement activity at a forest patch (including within-forest patch movement activity) showed unexpected effects for both plant-pollinator pairs that might be explained by accelerated genetic drift due to enhanced sexual reproduction. Our integrated approach revealed that bumblebees serve as genetic linkers of associated forest herb populations, even if they are more than 2 km apart from each other. No such evidence was found for the forest associated hoverfly species which showed significant genetic differentiation among forest patches itself. Our approach also indicated that a higher within-forest patch movement activity of both pollinator species might enhance sexual recruitment and thus diminishes the temporal buffer that clonal growth provides against habitat fragmentation effects.

Bumble bees exhibit body size clines across an urban gradient despite low genetic differentiation

Environmental heterogeneity resulting from human-modified landscapes can increase intraspecific trait variation. However, less known is whether such phenotypic variation is driven by plastic or adaptive responses to local environments. Here, we study five bumble bee (Apidae: Bombus) species across an urban gradient in the greater Saint Louis, Missouri region in the North American Midwest and ask: (1) Can urban environments induce intraspecific spatial structuring of body size, an ecologically consequential functional trait? And, if so, (2) is this body size structure the result of plasticity or adaptation? We additionally estimate genetic diversity, inbreeding, and colony density of these species—three factors that affect extinction risk. Using ≥ 10 polymorphic microsatellite loci per species and measurements of body size, we find that two of these species (Bombus impatiens, Bombus pensylvanicus) exhibit body size clines across the urban gradient, despite a lack of population genetic structure. We also reaffirm reports of low genetic diversity in B. pensylvanicus and find evidence that Bombus griseocollis, a species thought to be thriving in North America, is inbred in the greater Saint Louis region. Collectively, our results have implications for conservation in urban environments and suggest that plasticity can cause phenotypic clines across human-modified landscapes.

Sunflower pollen reduces a gut pathogen in the model bee species, Bombus impatiens, but has weaker effects in three wild congeners

Commercial bumblebees have become popular models to understand stressors and solutions for pollinator health, but few studies test whether results translate to other pollinators. Consuming sunflower pollen dramatically reduces infection by the gut parasite Crithidia bombi in commercially reared Bombus impatiens. We assessed the effect of sunflower pollen on infection in wild B. impatiens, Bombus griseocollis, Bombus bimaculatus and Bombus vagans. We also asked how pollen diet (50% sunflower pollen versus wildflower pollen) and infection (yes/no) affected performance in wild B. impatiens microcolonies. Compared to controls, sunflower pollen dramatically reduced Crithidia infection in commercial and wild B. impatiens, had similar but less dramatic effects in B. bimaculatus and B. vagans, and no effect in B. griseocollis. Bombus impatiens, B. bimaculatus and B. vagans are in the same subgenus, suggesting that responses to sunflower pollen may be phylogenetically conserved. In microcolonies, 50% sunflower pollen reduced infection compared to wildflower pollen, but also reduced reproduction. Sunflower pollen could control Crithidia infections in B. impatiens and potentially close relatives, but may hinder reproduction if other resources are scarce. We caution that research using managed bee species, such as B. impatiens, be interpreted carefully as findings may not relate to all bee species.

Maine's Bumble Bees (Hymenoptera: Apidae)-Part 2: Comparisons of a Common (Bombus ternarius) and a Rare (Bombus terricola) Species

As part of a quantitative survey of Maine's bumble bee fauna (Butler et al. 2021), we compared and contrasted genetic diversity, parasite and pathogen burdens, and pesticide exposure of the relatively common Bombus ternarius Say, 1937 and the spatially rare Bombus terricola Kirby, 1837. We recorded 11 Bombus species at 40 survey sites across three Maine ecoregions, and B. ternarius was the most common species, while B. terricola was spatially rare. Nonmetric multidimensional scaling indicated that B. terricola was associated with higher elevation sites in Maine, while B. ternarius was more broadly distributed in the state. Pollinator networks constructed for each bee indicated B. ternarius foraged on more plant species than B. terricola, but that there was considerable overlap (73%) in plant species visited. Genetic diversity was greater in the spatially restricted B. terricola, whereas the widely distributed B. ternarius was characterized by greater genetic differentiation among regions. Bombus terricola had higher molecular marker levels of the microsporidian fungi Nosema spp. and the trypanosome Crithidia spp., and both species had high levels of Trypanosoma spp. exposure. No Western Honey Bee (Apis mellifera, Linnaeus, 1758) viruses were detected in either species. Pesticides were not detected in pollen samples collected from workers of either species, and B. ternarius worker tissue samples exhibited only trace levels of diflubenzuron.

A genomic and morphometric analysis of alpine bumblebees: Ongoing reductions in tongue length but no clear genetic component

Over the last six decades, populations of the bumblebees Bombus sylvicola and Bombus balteatus in Colorado have experienced decreases in tongue length, a trait important for plant-pollinator mutualisms. It has been hypothesized that this observation reflects selection resulting from shifts in floral composition under climate change. Here we used morphometrics and population genomics to determine whether morphological change is ongoing, investigate the genetic basis of morphological variation, and analyse population structure in these populations. We generated a genome assembly of B. balteatus. We then analysed whole-genome sequencing data and morphometric measurements of 580 samples of both species from seven high-altitude localities. Out of 281 samples originally identified as B. sylvicola, 67 formed a separate genetic cluster comprising a newly-discovered cryptic species ("incognitus"). However, an absence of genetic structure within species suggests that gene flow is common between mountains. We found a significant decrease in tongue length between bees collected between 2012-2014 and in 2017, indicating that morphological shifts are ongoing. We did not discover any genetic associations with tongue length, but a SNP related to production of a proteolytic digestive enzyme was implicated in body size variation. We identified evidence of covariance between kinship and both tongue length and body size, which is suggestive of a genetic component of these traits, although it is possible that shared environmental effects between colonies are responsible. Our results provide evidence for ongoing modification of a morphological trait important for pollination and indicate that this trait probably has a complex genetic and environmental basis.

Parasite Prevalence May Drive the Biotic Impoverishment of New England (USA) Bumble Bee Communities

Simple Summary

Here we discuss widespread changes in the community structure of bumble bees (Bombus spp.) found in the coastal-zone community of New England. One species in particular, Bombus impatiens Cresson, 1863, has increased in relative abundance nearly 45% since the 1990s to become the dominant species in the region, representing nearly 75% of all Bombus individuals collected in our studies. These changes in abundance may be, in part, due to differences in infection rates by microparasites, with B. impatiens having significantly fewer microparasites than several other less common and declining Bombus species. We discuss the possible role of microparasites in influencing the community composition of Bombus species in our region, and how these infections might be compounding declines in conjunction with habitat loss and climate change.

Abstract

Numerous studies have reported a diversity of stressors that may explain continental-scale declines in populations of native pollinators, particularly those in the genus Bombus. However, there has been little focus on the identification of the local-scale dynamics that may structure currently impoverished Bombus communities. For example, the historically diverse coastal-zone communities of New England (USA) now comprise only a few species and are primarily dominated by a single species, B. impatiens. To better understand the local-scale factors that might be influencing this change in community structure, we examined differences in the presence of parasites in different species of Bombus collected in coastal-zone communities. Our results indicate that Bombus species that are in decline in this region were more likely to harbor parasites than are B. impatiens populations, which were more likely to be parasite-free and to harbor fewer intense infections or co-infections. The contrasting parasite burden between co-occurring winners and losers in this community may impact the endgame of asymmetric contests among species competing for dwindling resources. We suggest that under changing climate and landscape conditions, increasing domination of communities by healthy, synanthropic Bombus species (such as B. impatiens) may be another factor hastening the further erosion of bumble bee diversity.

Urbanization is associated with shifts in bumblebee body size, with cascading effects on pollination

Urbanization is a global phenomenon with major effects on species, the structure of community functional traits and ecological interactions. Body size is a key species trait linked to metabolism, life-history and dispersal as well as a major determinant of ecological networks. Here, using a well-replicated urban-rural sampling design in Central Europe, we investigate the direction of change of body size in response to urbanization in three common bumblebee species, Bombus lapidarius, Bombus pascuorum and Bombus terrestris, and potential knock-on effects on pollination service provision. We found foragers of B. terrestris to be larger in cities and the body size of all species to be positively correlated with road density (albeit at different, species-specific scales); these are expected consequences of habitat fragmentation resulting from urbanization. High ambient temperature at sampling was associated with both a small body size and an increase in variation of body size in all three species. At the community level, the community-weighted mean body size and its variation increased with urbanization. Urbanization had an indirect positive effect on pollination services through its effects not only on flower visitation rate but also on community-weighted mean body size and its variation. We discuss the eco-evolutionary implications of the effect of urbanization on body size, and the relevance of these findings for the key ecosystem service of pollination.

Urban evolution comes into its own: Emerging themes and future directions of a burgeoning field

Urbanization has recently emerged as an exciting new direction for evolutionary research founded on our growing understanding of rapid evolution paired with the expansion of novel urban habitats. Urbanization can influence adaptive and nonadaptive evolution in urban-dwelling species, but generalized patterns and the predictability of urban evolutionary responses within populations remain unclear. This editorial introduces the special feature "Evolution in Urban Environments" and addresses four major emerging themes, which include: (a) adaptive evolution and phenotypic plasticity via physiological responses to urban climate, (b) adaptive evolution via phenotype-environment relationships in urban habitats, (c) population connectivity and genetic drift in urban landscapes, and (d) human-wildlife interactions in urban spaces. Here, we present the 16 articles (12 empirical, 3 review, 1 capstone) within this issue and how they represent each of these four emerging themes in urban evolutionary biology. Finally, we discuss how these articles address previous questions and have now raised new ones, highlighting important new directions for the field.

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.

Beyond the Decline of Wild Bees: Optimizing Conservation Measures and Bringing Together the Actors

Wild bees are facing a global decline mostly induced by numerous human factors for the last decades. In parallel, public interest for their conservation increased considerably, namely through numerous scientific studies relayed in the media. In spite of this broad interest, a lack of knowledge and understanding of the subject is blatant and reveals a gap between awareness and understanding. While their decline is extensively studied, information on conservation measures is often scattered in the literature. We are now beyond the precautionary principle and experts are calling for effective actions to promote wild bee diversity and the enhancement of environment quality. In this review, we draw a general and up-to-date assessment of the conservation methods, as well as their efficiency and the current projects that try to fill the gaps and optimize the conservation measures. Targeting bees, we focused our attention on (i) the protection and restoration of wild bee habitats, (ii) the conservation measures in anthropogenic habitats, (iii) the implementation of human made tools, (iv) how to deal with invasive alien species, and finally (v) how to communicate efficiently and accurately. This review can be considered as a needed catalyst to implement concrete and qualitative conversation actions for 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.

Utilizing field collected insects for next generation sequencing: Effects of sampling, storage, and DNA extraction methods

DNA sequencing technologies continue to advance the biological sciences, expanding opportunities for genomic studies of non-model organisms for basic and applied questions. Despite these opportunities, many next generation sequencing protocols have been developed assuming a substantial quantity of high molecular weight DNA (>100 ng), which can be difficult to obtain for many study systems. In particular, the ability to sequence field-collected specimens that exhibit varying levels of DNA degradation remains largely unexplored. In this study we investigate the influence of five traditional insect capture and curation methods on Double-Digest Restriction Enzyme Associated DNA (ddRAD) sequencing success for three wild bee species. We sequenced a total of 105 specimens (between 7-13 specimens per species and treatment). We additionally investigated how different DNA quality metrics (including pre-sequence concentration and contamination) predicted downstream sequencing success, and also compared two DNA extraction methods. We report successful library preparation for all specimens, with all treatments and extraction methods producing enough highly reliable loci for population genetic analyses. Although results varied between species, we found that specimens collected by net sampling directly into 100% EtOH, or by passive trapping followed by 100% EtOH storage before pinning tended to produce higher quality ddRAD assemblies, likely as a result of rapid specimen desiccation. Surprisingly, we found that specimens preserved in propylene glycol during field sampling exhibited lower-quality assemblies. We provide recommendations for each treatment, extraction method, and DNA quality assessment, and further encourage researchers to consider utilizing a wider variety of specimens for genomic analyses.

Population genomics of Bombus terrestris reveals high but unstructured genetic diversity in a potential glacial refugium

Ongoing climate change is expected to cause an increase in temperature and a reduction of precipitation levels in the Mediterranean region, which might cause changes in many species distributions. These effects negatively influence species gene pools, decreasing genetic variability and adaptive potential. Here, we use mitochondrial DNA and RADseq to analyse population genetic structure and genetic diversity of the bumblebee species Bombus terrestris (subspecies Bombus terrestris lusitanicus), in the Iberian Peninsula. Although this subspecies shows a panmictic pattern of population structure across Iberia and beyond, we found differentiation between subspecies B. t. lusitanicus and B. t. africanus, probably caused by the existence of barriers to gene flow between Iberia and North Africa. Furthermore, the results revealed that the Iberian Peninsula harbours a large fraction of B. terrestris intraspecific genetic variation, with the highest number of mitochondrial haplotypes found when compared with any other region in Europe studied so far, suggesting a potential role for the Iberian Peninsula as a glacial refugium. Our findings strengthen the idea that Iberia is a very important source of diversity for the global genetic pool of this species, because rare alleles might play a role in population resilience against human- or climate-mediated changes.

Methylation and gene expression differences between reproductive and sterile bumblebee workers

Phenotypic plasticity is the production of multiple phenotypes from a single genome and is notably observed in social insects. Multiple epigenetic mechanisms have been associated with social insect plasticity, with DNA methylation being explored to the greatest extent. DNA methylation is thought to play a role in caste determination in Apis mellifera, and other social insects, but there is limited knowledge on its role in other bee species. In this study, we analyzed whole genome bisulfite sequencing and RNA-seq data sets from head tissue of reproductive and sterile castes of the eusocial bumblebee Bombus terrestris. We found that genome-wide methylation in B. terrestris is similar to other holometabolous insects and does not differ between reproductive castes. We did, however, find differentially methylated genes between castes, which are enriched for multiple biological processes including reproduction. However, we found no relationship between differential methylation and differential gene expression or differential exon usage between castes. Our results also indicate high intercolony variation in methylation. These findings suggest that methylation is associated with caste differences but may serve an alternate function, other than direct caste determination in this species. This study provides the first insights into the nature of a bumblebee caste-specific methylome as well as its interaction with gene expression and caste-specific alternative splicing, providing greater understanding of the role of methylation in phenotypic plasticity within social bee species. Future experimental work is needed to determine the function of methylation and other epigenetic mechanisms in insects.

Identification and functional characterisation of a novel N-cyanoamidine neonicotinoid metabolising cytochrome P450, CYP9Q6, from the buff-tailed bumblebee Bombus terrestris

Recent work has shown that two bumblebee (Bombus terrestris) cytochrome P450s of the CYP9Q subfamily, CYP9Q4 and CYP9Q5, are important biochemical determinants of sensitivity to neonicotinoid insecticides. Here, we report the characterisation of a third P450 gene CYP9Q6, previously mis-annotated in the genome of B. terrestris, encoding an enzyme that metabolises the N-cyanoamidine neonicotinoids thiacloprid and acetamiprid with high efficiency. The genomic location and complete ORF of CYP9Q6 was corroborated by PCR and its metabolic activity characterised in vitro by expression in an insect cell line. CYP9Q6 metabolises both thiacloprid and acetamiprid more rapidly than the previously reported CYP9Q4 and CYP9Q5. We further demonstrate a direct, in vivo correlation between the expression of the CYP9Q6 enzyme in transgenic Drosophila melanogaster and an increased tolerance to thiacloprid and acetamiprid. We conclude that CYP9Q6 is an efficient metaboliser of N-cyanoamidine neonicotinoids and likely plays a key role in the high tolerance of B. terrestris to these insecticides.

Pattern of population structuring between Belgian and Estonian bumblebees

Several population genetic studies investigated the extent of gene flow and population connectivity in bumblebees. In general, no restriction in gene flow is considered for mainland populations of common bumblebee species. Whether this assumption holds true for all species is not known. An assessment of bumblebee genetic structure in the context of their geographic distribution is needed to prioritize conservation and management needs. Here, we conducted a genetic study on seven bumblebee species occurring in Belgium and Estonia. Using 16 microsatellite markers, we investigated genetic diversity and population structuring in each species. This is the first study investigating population structuring of both declining and stable bumblebee species on both small and large geographic scales. Our results showed no or only low population structuring between the populations of the restricted and declining bumblebee species on both scales, while significant structuring was found for populations of the common species on the larger scale. The latter result, which may be due to human or environmental changes in the landscape, implies the need for the conservation of also widespread bumblebee species. Conservation strategies to improve gene flow and connectivity of populations could avoid the isolation and future losses of populations of these important species.

A trait-based approach to predict population genetic structure in bees

Understanding population genetic structure is key to developing predictions about species susceptibility to environmental change, such as habitat fragmentation and climate change. It has been theorized that life-history traits may constrain some species in their dispersal and lead to greater signatures of population genetic structure. In this study, we use a quantitative comparative approach to assess if patterns of population genetic structure in bees are driven by three key species-level life-history traits: body size, sociality, and diet breadth. Specifically, we reviewed the current literature on bee population genetic structure, as measured by the differentiation indices Nei's G_ST, Hedrick's G'_ST , and Jost's D. We then used phylogenetic generalised linear models to estimate the correlation between the evolution of these traits and patterns of genetic differentiation. Our analyses revealed a negative and significant effect of body size on genetic structure, regardless of differentiation index utilized. For Hedrick's G'_ST and Jost's D, we also found a significant impact of sociality, where social species exhibited lower levels of differentiation than solitary species. We did not find an effect of diet specialization on population genetic structure. Overall, our results suggest that physical dispersal or other functions related to body size are among the most critical for mediating population structure for bees. We further highlight the importance of standardizing population genetic measures to more easily compare studies and to identify the most susceptible species to landscape and climatic changes.

Spatial ecology of a range-expanding bumble bee pollinator

Molecular methods have greatly increased our understanding of the previously cryptic spatial ecology of bumble bees (Bombus spp.), with knowledge of the spatial ecology of these bees being central to conserving their essential pollination services. Bombus hypnorum, the Tree Bumble Bee, is unusual in that it has recently rapidly expanded its range, having colonized much of the UK mainland since 2001. However, the spatial ecology of B. hypnorum has not previously been investigated. To address this issue, and to investigate whether specific features of the spatial ecology of B. hypnorum are associated with its rapid range expansion, we used 14 microsatellite markers to estimate worker foraging distance, nest density, between-year lineage survival rate and isolation by distance in a representative UK B. hypnorum population. After assigning workers to colonies based on full or half sibship, we estimated the mean colony-specific worker foraging distance as 103.6 m, considerably less than values reported from most other bumble bee populations. Estimated nest density was notably high (2.56 and 0.72 colonies ha-1 in 2014 and 2015, respectively), estimated between-year lineage survival rate was 0.07, and there was no evidence of fine-scale isolation by distance. In addition, genotyping stored sperm dissected from sampled queens confirmed polyandry in this population (mean minimum mating frequency of 1.7 males per queen). Overall, our findings establish critical spatial ecological parameters and the mating system of this unusual bumble bee population and suggest that short worker foraging distances and high nest densities are associated with its rapid range expansion.

Diet and nutritional status during early adult life have immediate and persistent effects on queen bumble bees

Many insects sequester nutrients during developmentally programmed periods, which they metabolize during subsequent life history stages. During these periods, failure to store adequate nutrients can have persistent effects on fitness. Here, we examined a critical but under-studied nutrient storage period in queen bumble bees: the first days of adult life, which are followed by a diapause period typically coinciding with winter. We experimentally manipulated availability of pollen (the primary dietary source of lipids and protein) and the sugar concentration of artificial nectar (the primary source of carbohydrates) for laboratory-reared queens during this period and examined three nutritional phenomena: (i) diet impacts on nutritional status, (ii) the timescale upon which nutrient sequestration occurs and (iii) the fitness consequences of nutrient sequestration, specifically related to survival across the life cycle. We found evidence that pollen and nectar starvation negatively impact lipid storage, whereas nectar sugar concentration impacts stored carbohydrates. The majority of nutrients were stored during the first ~ 3 days of adult life. Nutrients derived from pollen during this period appear to be more critical for surviving earlier life stages, whereas nutrients sequestered from nectar become more important for surviving the diapause and post-diapause periods. Negative impacts of a poor diet during early life persisted in our experiment, even when pollen and a relatively high (50%) nectar sugar concentration were provided post-diapause. Based on these findings, we posit that the nutritional environment during the early adult life of queens has both immediate and persistent impacts on fitness. These findings underscore the importance of examining effects of stage-specific nutritional limitations on physiology and life history traits in this social insect group. Moreover, the findings may shed light on how declining food resources are contributing to the decline of wild bumble bee populations.

One size does not fit all: Caste and sex differences in the response of bumblebees (Bombus impatiens) to chronic oral neonicotinoid exposure

Neonicotinoid insecticides have been implicated in the rapid global decline of bumblebees over recent years, particularly in agricultural and urban areas. While there is much known about neonicotinoid toxicity effects at the colony stage of the bumblebee annual cycle, far less is known about such effects at other stages critical for the maintenance of wild populations. In the present work, individual-based feeding assays were used to show that chronic consumption of the widely used neonicotinoid clothianidin at a field-realistic average rate of 3.6 and 4.0 ng/g·bee/day reduces survival of queen and male bumblebees, respectively, within a 7-day period. In contrast, worker survival was unaffected at a similar consumption rate of 3.9 ng/g·bee/day. To test the hypothesis that males have a lower tolerance for oral clothianidin exposure than workers due to their haploid genetic status, RNAseq analysis was used to compare the transcriptomic responses of workers and males to chronic intake of clothianidin at a sub-lethal dose of 0.37ng/bee/day for 5 days. Surprisingly, clothianidin consumption only altered the expression of 19 putative detoxification genes in a sex-specific manner, with 11/19 genes showing increased expression in workers. Sub-lethal clothianidin exposure also altered the expression of 40 genes associated with other major biological functions, including locomotion, reproduction, and immunity. Collectively, these results suggest that chronic oral toxicity effects of neonicotinoids are greatest during mating and nest establishment phases of the bumblebee life cycle. Chronic oral toxicity testing on males and queens is therefore required in order to fully assess the impact of neonicotinoids on wild bumblebee populations.

Quaternary climate instability is correlated with patterns of population genetic variability in Bombus huntii

Climate oscillations have left a significant impact on the patterns of genetic diversity observed in numerous taxa. In this study, we examine the effect of Quaternary climate instability on population genetic variability of a bumble bee pollinator species, Bombus huntii in western North America. Pleistocene and contemporary B. huntii habitat suitability (HS) was estimated with an environmental niche model (ENM) by associating 1,035 locality records with 10 bioclimatic variables. To estimate genetic variability, we genotyped 380 individuals from 33 localities at 13 microsatellite loci. Bayesian inference was used to examine population structure with and without a priori specification of geographic locality. We compared isolation by distance (IBD) and isolation by resistance (IBR) models to examine population differentiation within and among the Bayesian inferred genetic clusters. Furthermore, we tested for the effect of environmental niche stability (ENS) on population genetic diversity with linear regression. As predicted, high-latitude B. huntii habitats exhibit low ENS when compared to low-latitude habitats. Two major genetic clusters of B. huntii inhabit western North America: (a) a north genetic cluster predominantly distributed north of 28°N and (b) a south genetic cluster distributed south of 28°N. In the south genetic cluser, both IBD and IBR models are significant. However, in the north genetic cluster, IBD is significant but not IBR. Furthermore, the IBR models suggest that low-latitude montane populations are surrounded by habitat with low HS, possibly limiting dispersal, and ultimately gene flow between populations. Finally, we detected high genetic diversity across populations in regions that have been climatically unstable since the last glacial maximum (LGM), and low genetic diversity across populations in regions that have been climatically stable since the LGM. Understanding how species have responded to climate change has the potential to inform management and conservation decisions of both ecological and economic concerns.

Fabrication of an Effective Avermectin Nanoemulsion Using a Cleavable Succinic Ester Emulsifier

In this study, a new emulsifier precursor was prepared via esterification of avermectin with succinic anhydride. The chemical structure of the product was confirmed to be monosubstituted avermectin. After neutralization with triethanolamine, it exhibited adequate emulsification ability for avermectin. Avermectin was then encapsulated in nanoparticles in the nanoemulsion with a high drug loading up to 60 wt % and high stability. The nanoemulsion of nanoparticles that serves as a carrier of avermectin shows highly efficient pesticide characteristics, including low surface tension, high affinity to leaves, and improved photostability. In the presence of esterase or under strongly basic conditions, the ester bonds of the emulsifier can be hydrolyzed, and the encapsulated avermectin molecules can be released in an accelerated manner. The nanoemulsion exhibited improved insecticidal effect compared with commercial emulsifiable concentrate, which was attributed to the cleavage of ester bonds of the emulsifier by esterase in vivo.

Fabrication of Novel Avermectin Nanoemulsion Using a Polyurethane Emulsifier with Cleavable Disulfide Bonds

In this study, a polyurethane emulsifer with various functional groups was prepared from isophorone diisocyanate, avermectin, 2,2-dimethylol propionic acid, and bis(2-hydroxyethyl) disulfide. The chemical structure of the polymer was confirmed by Fourier transform infrared spectroscopy, proton nuclear magnetic resonance, and element analysis. The polymer exhibited adequate emulsification ability for avermectin after neutralization with triethylamine. A satisfaying nanoemulsion was obtained, in which avermectin was encapsulated in nanoparticles with 50 wt % drug loading, low organic solvent content, and high stability under dilution and centrifuging treatment in addition to low surface tension, high affinity to crop leaf, and improved avermectin photostability. The resulting nanoparticles showed degradability in the presence of dl-dithiothreitol or inside the insect as a result of the disulfide bonds, promoting the release of avermectin. As a result, the avermectin nanoparticles showed higher insecticidal ability compared to both the avermectin nanoparticles without a disulfide group and the avermectin emulsifiable concentrate.

Genome-wide single nucleotide polymorphism scan suggests adaptation to urbanization in an important pollinator, the red-tailed bumblebee (Bombus lapidarius L.)

Urbanization is considered a global threat to biodiversity; the growth of cities results in an increase in impervious surfaces, soil and air pollution, fragmentation of natural vegetation and invasion of non-native species, along with numerous environmental changes, including the heat island phenomenon. The combination of these effects constitutes a challenge for both the survival and persistence of many native species, while also imposing altered selective regimes. Here, using 110 314 single nucleotide polymorphisms generated by restriction-site-associated DNA sequencing, we investigated the genome-wide effects of urbanization on putative neutral and adaptive genomic diversity in a major insect pollinator, Bombus lapidarius, collected from nine German cities and nine paired rural sites. Overall, genetic differentiation among sites was low and there was no obvious genome-wide genetic structuring, suggesting the absence of strong effects of urbanization on gene flow. We nevertheless identified several loci under directional selection, a subset of which was associated with urban land use, including the percentage of impervious surface surrounding each sampling site. Overall, our results provide evidence of local adaptation to urbanization in the face of gene flow in a highly mobile insect pollinator.

Population genetic structure of Diaphorina citri Kuwayama (Hemiptera: Liviidae): host-driven genetic differentiation in China

The Asian citrus psyllid Diaphorina citri Kuwayama is a major pest in citrus production, transmitting Candidatus Liberibacter asiaticus. It has spread widely across eastern and southern China. Unfortunately, little is known about the genetic diversity and population structure of D. citri, making pest control difficult. In this study, nine specifically developed SSR markers and three known mitochondrial DNA were used for population genetics study of D. citri using 225 samples collected from all 7 distribution regions in China. Based on the SSR data, D. citri was found highly diverse with a mean observed heterozygosity of 0.50, and three subgroups were structured by host plant: (i) Shatangju, NF mandarin and Ponkan; (ii) Murraya paniculata and Lemon; (iii) Citrus unshiu, Bingtangcheng, Summer orange and Navel. No significant genetic differences were found with mtDNA data. We suggested the host-associated divergence is likely to have occurred very recently. A unimodal distribution of paired differences, the negative and significant Tajima’s D and Fu’s F_S parameters among mtDNA suggested a recent demographic expansion. The extensive citrus cultivation and increased suitable living habitat was recommended as a key for this expansion event.

Bumble bee ecophysiology: integrating the changing environment and the organism

Bumble bees are among the most ecologically and economically important pollinators worldwide, yet many of their populations are being threatened by a suite of interrelated, human-mediated environmental changes. Here, I discuss recent progress in our understanding of bumble bee ecophysiology, including advances related to thermal biology in light of global warming; nutritional biology in the context of declining food resources; and the capacity for bumble bees to exhibit physiological plasticity or adaptations to novel or extreme environments, with reference to their evolutionary history and current biogeography.

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.