Estimation of bumblebee queen dispersal distances using sibship reconstruction method

Limited introgression from non-native commercial strains and signatures of adaptation in the key pollinator Bombus terrestris

Insect pollination is fundamental for natural ecosystems and agricultural crops. The bumblebee species Bombus terrestris has become a popular choice for commercial crop pollination worldwide due to its effectiveness and ease of mass rearing. Bumblebee colonies are mass produced for the pollination of more than 20 crops and imported into over 50 countries including countries outside their native ranges, and the risk of invasion by commercial non-native bumblebees is considered an emerging issue for global conservation and biological diversity. Here, we use genome-wide data from seven wild populations close to and far from farms using commercial colonies, as well as commercial populations, to investigate the implications of utilizing commercial bumblebee subspecies in the UK. We find evidence for generally low levels of introgression between commercial and wild bees, with higher admixture proportions in the bees occurring close to farms. We identify genomic regions putatively involved in local and global adaptation, and genes in locally adaptive regions were found to be enriched for functions related to taste receptor activity, oxidoreductase activity, fatty acid and lipid biosynthetic processes. Despite more than 30 years of bumblebee colony importation into the UK, we observe low impact on the genetic integrity of local B. terrestris populations, but we highlight that even limited introgression might negatively affect locally adapted populations.

Occupancy of Bombus affinis (Hymenoptera: Apidae) in Minnesota is highest in developed areas when standardized surveys are employed

Mounting evidence of bumble bee declines and the listing of the rusty patched bumble bee (Bombus affinis Cresson) as federally endangered in the United States in 2017 and Canada in 2012 has stimulated an interest in monitoring and conservation. Understanding the influence of land use on occupancy patterns of imperiled species is crucial to successful recovery planning. Using detection data from community surveys, we assessed land use associations for 7 bumble bee species in Minnesota, USA, including B. affinis. We used multispecies occupancy models to assess the effect of 3 major land use types (developed, agricultural, and natural) within 0.5 and 1.5 km on occupancy of 7 Bombus (Hymenoptera: Apidae) species, while accounting for detection uncertainty. We found that B. affinis occupancy and detection were highest in developed landscapes and lowest in agricultural landscapes, representing an inverse relationship with the relative landcover ratios of these landscapes in Minnesota. Occupancy of 2 bumble bee species had strong positive associations with natural landscapes within 1.5 km and 2 species had strong negative associations with agricultural landscapes within 1.5 km. Our results suggest that best practices for imperiled Bombus monitoring and recovery planning depends upon the surrounding major land use patterns. We provide recommendations for urban planning to support B. affinis based on conservation success in the metropolitan areas of Minneapolis-St. Paul. We also encourage substantial survey effort be employed in agricultural and natural regions of the state historically occupied by B. affinis to determine the current occupancy state.

Dispersal inference from population genetic variation using a convolutional neural network

The geographic nature of biological dispersal shapes patterns of genetic variation over landscapes, making it possible to infer properties of dispersal from genetic variation data. Here, we present an inference tool that uses geographically distributed genotype data in combination with a convolutional neural network to estimate a critical population parameter: the mean per-generation dispersal distance. Using extensive simulation, we show that our deep learning approach is competitive with or outperforms state-of-the-art methods, particularly at small sample sizes. In addition, we evaluate varying nuisance parameters during training-including population density, demographic history, habitat size, and sampling area-and show that this strategy is effective for estimating dispersal distance when other model parameters are unknown. Whereas competing methods depend on information about local population density or accurate inference of identity-by-descent tracts, our method uses only single-nucleotide-polymorphism data and the spatial scale of sampling as input. Strikingly, and unlike other methods, our method does not use the geographic coordinates of the genotyped individuals. These features make our method, which we call "disperseNN," a potentially valuable new tool for estimating dispersal distance in nonmodel systems with whole genome data or reduced representation data. We apply disperseNN to 12 different species with publicly available data, yielding reasonable estimates for most species. Importantly, our method estimated consistently larger dispersal distances than mark-recapture calculations in the same species, which may be due to the limited geographic sampling area covered by some mark-recapture studies. Thus genetic tools like ours complement direct methods for improving our understanding of dispersal.

Pollinator conservation: Where will bees go in the Anthropocene future?

For pollinator conservation in human-transformed landscapes, it is crucial to know whether species can overcome gaps between fragments of natural habitat. A new study reveals why colony size, recruitment communication, and flower constancy increase the foraging ranges in social bees.

Signals of adaptation to agricultural stress in the genomes of two European bumblebees

Human-induced environmental impacts on wildlife are widespread, causing major biodiversity losses. One major threat is agricultural intensification, typically characterised by large areas of monoculture, mechanical tillage, and the use of agrochemicals. Intensification leads to the fragmentation and loss of natural habitats, native vegetation, and nesting and breeding sites. Understanding the adaptability of insects to these changing environmental conditions is critical to predicting their survival. Bumblebees, key pollinators of wild and cultivated plants, are used as model species to assess insect adaptation to anthropogenic stressors. We investigated the effects of agricultural pressures on two common European bumblebees, Bombus pascuorum and B. lapidarius. Restriction-site Associated DNA Sequencing was used to identify loci under selective pressure across agricultural-natural gradients over 97 locations in Europe. 191 unique loci in B. pascuorum and 260 in B. lapidarius were identified as under selective pressure, and associated with agricultural stressors. Further investigation suggested several candidate proteins including several neurodevelopment, muscle, and detoxification proteins, but these have yet to be validated. These results provide insights into agriculture as a stressor for bumblebees, and signal for conservation action in light of ongoing anthropogenic changes.

Biogeography and Diversification of Bumblebees (Hymenoptera: Apidae), with Emphasis on Neotropical Species

A detailed phylogeny of bumblebees is urgently needed to understand speciation and biogeographic diversification in the Neotropical region. We sequenced autosomal and mtDNA loci from nine Brazilian bumblebee species and compiled it with the data already available to obtain highly resolved phylogenetic trees with fossil-calibrated dates. The ancestral Bombus lineage was estimated to diversify between 47.08 and 34.27 million years ago (Ma) in the Holarctic region, but largely restricted to the eastern Old World. The Neotropical region was initially colonized in the Late Miocene, where bumblebee diversification was shown to be consistent with geologic and climatic events of the Late Cenozoic. Neotropical bumblebees likely originated from Nearctic lineages, which dispersed towards South America after 29 Ma.

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.

Environmental differences explain subtle yet detectable genetic structure in a widespread pollinator

BACKGROUND

The environment is a strong driver of genetic structure in many natural populations, yet often neglected in population genetic studies. This may be a particular problem in vagile species, where subtle structure cannot be explained by limitations to dispersal. Consequently, these species might falsely be considered quasi-panmictic and hence potentially mismanaged. A species this might apply to, is the buff-tailed bumble bee (Bombus terrestris), an economically important and widespread pollinator, which is considered to be quasi-panmictic at mainland continental scales. Here we aimed to (i) quantify genetic structure in 21+ populations of the buff-tailed bumble bee, sampled throughout two Eastern European countries, and (ii) analyse the degree to which structure is explained by environmental differences, habitat permeability and geographic distance. Using 12 microsatellite loci, we characterised populations of this species with Fst analyses, complemented by discriminant analysis of principal components and Bayesian clustering approaches. We then applied generalized dissimilarity modelling to simultaneously assess the informativeness of geographic distance, habitat permeability and environmental differences among populations in explaining divergence.

RESULTS

Genetic structure of the buff-tailed bumble bee quantified by means of Fst was subtle and not detected by Bayesian clustering. Discriminant analysis of principal components suggested insignificant but still noticeable structure that slightly exceeded estimates obtained through Fst analyses. As expected, geographic distance and habitat permeability were not informative in explaining the spatial pattern of genetic divergence. Yet, environmental variables related to temperature, vegetation and topography were highly informative, explaining between 33 and 39% of the genetic variation observed.

CONCLUSIONS

In contrast to previous studies reporting quasi-panmixia in continental populations of this species, we demonstrated the presence of subtle population structure related to environmental heterogeneity. Environmental data proved to be highly useful in unravelling the drivers of genetic structure in this vagile and opportunistic species. We highlight the potential of including these data to obtain a better understanding of population structure and the processes driving it in species considered to be quasi-panmictic.

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.

Genetic and Usurpation Data Support High Incidence of Bumble Bee Nest Invasion by Socially Parasitic Bumble Bee, Bombus insularis

Cuckoo bumble bees (Psithyrus) (Lepeletier, 1832) (Hymenoptera: Apidae) are a unique lineage of bees that depend exclusively on a host bumble bee species to provide nesting material, nutritional resources, and labor to rear offspring. In this study, we document usurpation incidence and population genetic data of Bombus insularis (Smith, 1861) (Hymenoptera: Apidae), a bumble bee species in the Psithyrus subgenus, on field-deployed B. huntii colonies in northern Utah, United States. Within 12 d of deploying B. huntii Greene, 1860 (Hymenoptera: Apidae) colonies at two field sites, 13 of the 16 colonies contained at least one established B. insularis female. Although our results demonstrate that field-deployed bumble bee colonies are highly susceptible to B. insularis usurpation, applying a fabricated excluder to prevent the inquiline from invading a colony was 100% effective. Sibship analysis using microsatellite genotype data of 59 B. insularis females estimates that they originated from at least 49 unique colonies. Furthermore, sibship analysis found siblings distributed between the field sites that were 7.04 km apart. Our result suggests that B. insularis females have the capacity to disperse across the landscape in search of host colonies at distances of at least 3.52 km and up to 7.04 km. Our study underscores the detrimental impact B. insularis usurpation has on the host bumble bee colony. As B. insularis significantly impacts the success of bumble bee colonies, we briefly discuss how the utilization of excluders may be useful for commercial bumble bee colonies that are used to pollinate open field crops.

Mating systems and predictors of relative reproductive success in a Cutthroat Trout subspecies of conservation concern

Mating systems and patterns of reproductive success in fishes play an important role in ecology and evolution. While information on the reproductive ecology of many anadromous salmonids (Oncorhynchus spp.) is well detailed, there is less information for nonanadromous species including the Yellowstone Cutthroat Trout (O. clarkii bouvieri), a subspecies of recreational angling importance and conservation concern. Using data from a parentage-based tagging study, we described the genetic mating system of a migratory population of Yellowstone Cutthroat Trout, tested for evidence of sexual selection, and identified predictors of mating and reproductive success. The standardized variance in mating success (i.e., opportunity for sexual selection) was significantly greater for males relative to females, and while the relationship between mating success and reproductive success (i.e., Bateman gradient) was significantly positive for both sexes, a greater proportion of reproductive success was explained by mating success for males (r 2 = 0.80) than females (r 2 = 0.59). Overall, the population displayed a polygynandrous mating system, whereby both sexes experienced variation in mating success due to multiple mating, and sexual selection was variable across sexes. Tests for evidence of sexual selection indicated the interaction between mating success and total length best-predicted relative reproductive success. We failed to detect a signal of inbreeding avoidance among breeding adults, but the group of parents that produced progeny were on average slightly less related than adults that did not produce progeny. Lastly, we estimated the effective number of breeders (N b) and effective population size (N e) and identified while N b was lower than N e, both are sufficiently high to suggest Yellowstone Cutthroat Trout in Burns Creek represent a genetically stable and diverse population.

Living with relatives offsets the harm caused by pathogens in natural populations

Living with relatives can be highly beneficial, enhancing reproduction and survival. High relatedness can, however, increase susceptibility to pathogens. Here, we examine whether the benefits of living with relatives offset the harm caused by pathogens, and if this depends on whether species typically live with kin. Using comparative meta-analysis of plants, animals, and a bacterium (nspecies = 56), we show that high within-group relatedness increases mortality when pathogens are present. In contrast, mortality decreased with relatedness when pathogens were rare, particularly in species that live with kin. Furthermore, across groups variation in mortality was lower when relatedness was high, but abundances of pathogens were more variable. The effects of within-group relatedness were only evident when pathogens were experimentally manipulated, suggesting that the harm caused by pathogens is masked by the benefits of living with relatives in nature. These results highlight the importance of kin selection for understanding disease spread in natural populations.

Flexibility in the Critical Period of Nutrient Sequestration in Bumble Bee Queens

Synopsis

Bumble bee queens undergo a nutrient storage period prior to entering diapause wherein they sequester glycogen and lipids that are metabolized during overwintering. In the laboratory under optimal food availability conditions, the majority of nutrients are sequestered during the first few days of adulthood. However, if food resources are scarce during this narrow window of time, wild queen bumble bees might be limited in their ability to obtain adequate food resources for overwintering. Here we used a laboratory experiment to examine whether queen bumble bees exhibit flexibility in the timing of pre-overwintering nutrient sequestration, by limiting their access to either nectar (artificial) or pollen, the two primary foods for bumble bees, for varying periods of time. In response to these treatments, we quantified queen survival, changes in weight, and glycogen and lipids levels. We found evidence that queens are able to recuperate almost entirely from food resource limitation, with respect to nutrient storage, especially when it is experienced for shorter durations (up to 6 days). This study sheds light on how bumble bee queens are impacted by food resource availability at a critical life stage.

Portuguese

As abelhas rainhas do gênero Bombus armazenam nutrientes antes de entrarem em diapausa, sequestrando o glicogênio e os lipídios que serão metabolizados durante o inverno. Em condições ideais de disponibilidade de alimento no laboratório, a maioria dos nutrientes é sequestrada nos primeiros dias de vida adulta. No entanto, em condições de escassez de alimento na natureza, as rainhas podem sofrer limitações em sua capacidade de obter recursos para o inverno. Nesse contexto, em condições controladas, examinamos se as rainhas exibem variações no sequestro de nutrientes, limitando o acesso ao néctar (artificial) ou pólen, seus principais alimentos, em diferentes intervalos de tempo. Em resposta a esses tratamentos, quantificamos a taxa de sobrevivência das rainhas, as mudanças no peso e os níveis de glicogênio e lipídios. Encontramos evidências de que as rainhas são capazes de recuperar a capacidade de armazenar nutrientes quase inteiramente, especialmente em períodos mais curtos de escassez de alimento (até 6 dias). Este estudo lança luz sobre como as rainhas são afetadas pela variação na disponibilidade de recursos alimentares em um estágio crítico da vida.

Spanish

Las abejas reinas de generó Bombus, mejor conocidas como reinas de abejorro se someten a un período de almacenamiento de nutrientes antes de entrar en diapausa, en el cual secuestran glucógeno y lípidos que se metabolizan durante el invierno. En el laboratorio, en condiciones óptimas de disponibilidad de alimentos, la mayoría de los nutrientes se secuestran durante los primeros días de la edad adulta. Sin embargo, si los recursos alimenticios son escasos durante esta estrecha ventana de tiempo, las abejas reinas silvestres podrían verse limitadas en su capacidad para obtener recursos alimenticios adecuados para pasar el invierno. Aquí utilizamos un experimento de laboratorio para examinar si las abejas reinas exhiben flexibilidad en el momento del secuestro de nutrientes antes de la hibernación, al limitar su acceso al néctar (artificial) o al polen, los dos alimentos principales de los abejorros, durante períodos variables. En respuesta a estos tratamientos, cuantificamos la supervivencia de la reina, los cambios de peso y los niveles de glucógeno y lípidos. Encontramos evidencia de que las reinas pueden recuperarse casi por completo de la limitación de los recursos alimenticios, con respecto al almacenamiento de nutrientes, especialmente cuando se experimenta por períodos más cortos (hasta 6 días). Este estudio arroja luz sobre cómo las abejas reinas se ven afectadas por la disponibilidad de recursos alimenticios en una etapa crítica de la vida.

Corridors through time: Does resource continuity impact pollinator communities, populations, and individuals?

Spatial aspects of connectivity have received considerable attention from ecologists and conservationists, yet temporal connectivity, the periodic linking of habitats, plays an equally important, but largely overlooked role. Different biological and biophysical attributes of ecosystems underpin temporal connectivity, but here we focus on resource continuity, the uninterrupted availability of foraging sites. We test the response of pollinators to resource continuity at community, population, and individual levels using a novel natural experiment consisting of farms with either single or sequential cropping systems. We found significant effects at the population level; colony density of an important crop pollinator (Bombus impatiens L.) was greater when crop floral resources were continuously available. However, we did not find significant effects at the community or individual level; wild bee abundance, diversity and body size did not respond to resource continuity. Raspberry farms with greater early season resources provided by blueberry had greater bumble bee populations, suggesting beneficial effects on resource availability due to crop diversity. Better understanding the impact of resource continuity via crop diversity on broader patterns of biodiversity is essential for the co-management of biodiversity and ecosystem services.

Connection, isolation and reconnection: Quaternary climatic oscillations and the Andes shaped the phylogeographical patterns of the Patagonian bee Centris cineraria (Apidae)

The joint effect of the Andes as a geographical barrier and the Quaternary glaciations as promoters of genetic divergence remains virtually unexplored in southern South America. To help fill this knowledge gap, in this study we investigated the demographic history of Centris cineraria, a solitary bee mainly distributed in Patagonia. We used mitochondrial and nuclear markers and performed phylogeographical and dating analyses, adjusted spatio-temporal diffusion and species distribution models, and used Approximate Bayesian Computation to identify likely historical demographic scenarios. Our results revealed that during glacial periods the Andes represented a barrier due to the extent of the ice-sheets and the occurrence of unsuitable habitats, while interglacials allowed for gene flow across the Andes. Secondary contact between previously isolated lineages was evident across at least two low-altitude Andean areas, the northern one being a putative glacial refugium. Our findings also suggest that C. cineraria has persisted in situ in four periglacial refugia located along a north–south transect, congruent with the maximum extent of the ice sheet during the Greatest Patagonian Glaciation. As the first phylogeographical study of Patagonian insects, our work reveals that the interaction between Quaternary climatic oscillations and the Andes as a barrier was the main driver of the spatial and demographic history of C. cineraria.

Using spatial genetics to quantify mosquito dispersal for control programs

Background

Hundreds of millions of people get a mosquito-borne disease every year and nearly one million die. Transmission of these infections is primarily tackled through the control of mosquito vectors. The accurate quantification of mosquito dispersal is critical for the design and optimization of vector control programs, yet the measurement of dispersal using traditional mark-release-recapture (MRR) methods is logistically challenging and often unrepresentative of an insect’s true behavior. Using Aedes aegypti (a major arboviral vector) as a model and two study sites in Singapore, we show how mosquito dispersal can be characterized by the spatial analyses of genetic relatedness among individuals sampled over a short time span without interruption of their natural behaviors.

Results

Using simple oviposition traps, we captured adult female Ae. aegypti across high-rise apartment blocks and genotyped them using genome-wide SNP markers. We developed a methodology that produces a dispersal kernel for distance which results from one generation of successful breeding (effective dispersal), using the distance separating full siblings and 2nd- and 3rd-degree relatives (close kin). The estimated dispersal distance kernel was exponential (Laplacian), with a mean dispersal distance (and dispersal kernel spread σ) of 45.2 m (95% CI 39.7–51.3 m), and 10% probability of a dispersal > 100 m (95% CI 92–117 m). Our genetically derived estimates matched the parametrized dispersal kernels from previous MRR experiments. If few close kin are captured, a conventional genetic isolation-by-distance analysis can be used, as it can produce σ estimates congruent with the close-kin method if effective population density is accurately estimated. Genetic patch size, estimated by spatial autocorrelation analysis, reflects the spatial extent of the dispersal kernel “tail” that influences, for example, the critical radii of release zones and the speed of Wolbachia spread in mosquito replacement programs.

Conclusions

We demonstrate that spatial genetics can provide a robust characterization of mosquito dispersal. With the decreasing cost of next-generation sequencing, the production of spatial genetic data is increasingly accessible. Given the challenges of conventional MRR methods, and the importance of quantified dispersal in operational vector control decisions, we recommend genetic-based dispersal characterization as the more desirable means of parameterization.

Wildfire reveals transient changes to individual traits and population responses of a native bumble bee Bombus vosnesenskii

Fire-induced changes in the abundance and distribution of organisms, especially plants, can alter resource landscapes for mobile consumers driving bottom-up effects on their population sizes, morphologies and reproductive potential. We expect these impacts to be most striking for obligate visitors of plants, like bees and other pollinators, but these impacts can be difficult to interpret due to the limited information provided by forager counts in the absence of survival or fitness proxies. Increased bumble bee worker abundance is often coincident with the pulses of flowers that follow recent fire. However, it is unknown if observed postfire activity is due to underlying population growth or a stable pool of colonies recruiting more foragers to abundant resource patches. This distinction is necessary for determining the net impact of disturbance on bumble bees: are there population-wide responses or do just a few colonies reap the rewards? We estimated colony abundance before and after fire in burned and unburned areas using a genetic mark-recapture framework. We paired colony abundance estimates with measures of body size, counts of queens, and estimates of foraging and dispersal to assess changes in worker size, reproductive output, and landscape-scale movements. Higher floral abundance following fire not only increased forager abundance but also the number of colonies from which those foragers came. Importantly, despite a larger population size, we also observed increased mean worker size. Two years following fire, queen abundance was higher in both burned and unburned sites, potentially due to the dispersal of queens from burned into unburned areas. The effects of fire were transient; within two growing seasons, worker abundance was substantially reduced across the entire sampling area and body sizes were similar between burned and unburned sites. Our results reveal how disturbance can temporarily release populations from resource limitation, boosting the genetic diversity, body size, and reproductive output of populations. Given that the effects of fire on bumble bees acted indirectly through pulsed resource availability, it is likely our results are generalizable to other situations, such as habitat restorations, where resource density is enhanced within the landscape.

Seed dispersal and realized gene flow of two forest orchids in a fragmented landscape

Species with vast production of dust-like windborne seeds, such as orchids, should not be limited by seed dispersal. This paradigm, however, does not fit recent studies showing that many sites suitable for orchids are unoccupied and most seeds land close to their maternal plant. To explore this issue, we studied seed dispersal and gene flow of two forest orchid species, Epipactis atrorubens and Cephalanthera rubra, growing in a fragmented landscape of forested limestone hills in southwest Bohemia, Czech Republic. We used a combination of seed trapping and plant genotyping methods (microsatellite DNA markers) to quantify short- and long-distance dispersal, respectively. In addition, seed production of both species was estimated. We found that most seeds landed very close to maternal plants (95% of captured seeds were within 7.2 m) in both species, and dispersal distance was influenced by forest type in E. atrorubens. In addition, C. rubra showed clonal reproduction (20% of plants were of clonal origin) and very low fruiting success (only 1.6% of plants were fruiting) in comparison with E. atrorubens (25.7%). Gene flow was frequent up to 2 km in C. rubra and up to 125 km in E. atrorubens, and we detected a relatively high dispersal rate among regions in both species. Although both species occupy similar habitats and have similar seed dispersal abilities, C. rubra is notably rarer in the study area. Considerably low fruiting success in this species likely limits its gene flow to longer distances and designates it more sensitive to habitat loss and fragmentation.

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.

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.

Harmonic radar tracking reveals random dispersal pattern of bumblebee (Bombus terrestris) queens after hibernation

The dispersal of animals from their birth place has profound effects on the immediate survival and longer-term persistence of populations. Molecular studies have estimated that bumblebee colonies can be established many kilometers from their queens’ natal nest site. However, little is known about when and how queens disperse during their lifespan. One possible life stage when dispersal may occur, is directly after emerging from hibernation. Here, harmonic radar tracking of artificially over-wintered Bombus terrestris queens shows that they spend most of their time resting on the ground with intermittent very short flights (duration and distance). We corroborate these behaviors with observations of wild queen bees, which show similar prolonged resting periods between short flights, indicating that the behavior of our radar-monitored bees was not due to the attachment of transponders nor an artifact of the bees being commercially reared. Radar-monitored flights were not continuously directed away from the origin, suggesting that bees were not intentionally trying to disperse from their artificial emergence site. Flights did not loop back to the origin suggesting bees were not trying to remember or get back to the original release site. Most individuals dispersed from the range of the harmonic radar within less than two days and did not return. Flight directions were not different from a uniform distribution and flight lengths followed an exponential distribution, both suggesting random dispersal. A random walk model based on our observed data estimates a positive net dispersal from the origin over many flights, indicating a biased random dispersal, and estimates the net displacement of queens to be within the range of those estimated in genetic studies. We suggest that a distinct post-hibernation life history stage consisting mostly of rest with intermittent short flights and infrequent foraging fulfils the dual purpose of ovary development and dispersal prior to nest searching.

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.

Climate change-driven range losses among bumblebee species are poised to accelerate

Climate change has shaped bee distributions over the past century. Here, we conducted the first species-specific assessment of future climate change impacts on North American bumblebee distributions, using the most recent global change scenarios developed in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC). We assessed potential shifts in bumblebee species distributions with models generated using Maxent. We tested different assumptions about bumblebee species’ dispersal capacities, drawing on observed patterns of range shifts to date, dispersal rates observed for bumblebee queens, and, lastly, assuming unlimited dispersal. Models show significant contractions of current ranges even under scenarios in which dispersal rates were high. Results suggest that dispersal rates may not suffice for bumblebees to track climate change as rapidly as required under any IPCC scenario for future climate change. Areas where species losses are projected overlap for many species and climate scenarios, and are concentrated in eastern parts of the continent. Models also show overlap for range expansions across many species, suggesting the presence of “hotspots” where management activities could benefit many species, across all climate scenarios. Broad-scale strategies are likely to be necessary to improve bumblebee conservation prospects under climate change.

Lower bumblebee colony reproductive success in agricultural compared with urban environments

Urbanization represents a rapidly growing driver of land-use change. While it is clear that urbanization impacts species abundance and diversity, direct effects of urban land use on animal reproductive success are rarely documented. Here, we show that urban land use is linked to long-term colony reproductive output in a key pollinator. We reared colonies from wild-caught bumblebee (Bombus terrestris) queens, placed them at sites characterized by varying degrees of urbanization from inner city to rural farmland and monitored the production of sexual offspring across the entire colony cycle. Our land-use cluster analysis identified three site categories, and this categorization was a strong predictor of colony performance. Crucially, colonies in the two clusters characterized by urban development produced more sexual offspring than those in the cluster dominated by agricultural land. These colonies also reached higher peak size, had more food stores, encountered fewer parasite invasions and survived for longer. Our results show a link between urbanization and bumblebee colony reproductive success, supporting the theory that urban areas provide a refuge for pollinator populations in an otherwise barren agricultural landscape.

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.

Pollinator population size and pollination ecosystem service responses to enhancing floral and nesting resources

Modeling pollination ecosystem services requires a spatially explicit, process-based approach because they depend on both the behavioral responses of pollinators to the amount and spatial arrangement of habitat and on the within- and between-season dynamics of pollinator populations in response to land use. We describe a novel pollinator model predicting flower visitation rates by wild central-place foragers (e.g., nesting bees) in spatially explicit landscapes. The model goes beyond existing approaches by: (1) integrating preferential use of more rewarding floral and nesting resources; (2) considering population growth over time; (3) allowing different dispersal distances for workers and reproductives; (4) providing visitation rates for use in crop pollination models. We use the model to estimate the effect of establishing grassy field margins offering nesting resources and a low quantity of flower resources, and/or late-flowering flower strips offering no nesting resources but abundant flowers, on bumble bee populations and visitation rates to flowers in landscapes that differ in amounts of linear seminatural habitats and early mass-flowering crops. Flower strips were three times more effective in increasing pollinator populations and visitation rates than field margins, and this effect increased over time. Late-blooming flower strips increased early-season visitation rates, but decreased visitation rates in other late-season flowers. Increases in population size over time in response to flower strips and amounts of linear seminatural habitats reduced this apparent competition for pollinators. Our spatially explicit, process-based model generates emergent patterns reflecting empirical observations, such that adding flower resources may have contrasting short- and long-term effects due to apparent competition for pollinators and pollinator population size increase. It allows exploring these effects and comparing effect sizes in ways not possible with other existing models. Future applications include species comparisons, analysis of the sensitivity of predictions to life-history traits, as well as large-scale management intervention and policy assessment.

Comparative phylogeography in the Atlantic forest and Brazilian savannas: pleistocene fluctuations and dispersal shape spatial patterns in two bumblebees

Background

Bombus morio and B. pauloensis are sympatric widespread bumblebee species that occupy two major Brazilian biomes, the Atlantic forest and the savannas of the Cerrado. Differences in dispersion capacity, which is greater in B. morio, likely influence their phylogeographic patterns. This study asks which processes best explain the patterns of genetic variation observed in B. morio and B. pauloensis, shedding light on the phenomena that shaped the range of local populations and the spatial distribution of intra-specific lineages.

Results

Results suggest that Pleistocene climatic oscillations directly influenced the population structure of both species. Correlative species distribution models predict that the warmer conditions of the Last Interglacial contributed to population contraction, while demographic expansion happened during the Last Glacial Maximum. These results are consistent with physiological data suggesting that bumblebees are well adapted to colder conditions. Intra-specific mitochondrial genealogies are not congruent between the two species, which may be explained by their documented differences in dispersal ability.

Conclusions

While populations of the high-dispersal B. morio are morphologically and genetically homogeneous across the species range, B. pauloensis encompasses multiple (three) mitochondrial lineages, and show clear genetic, geographic, and morphological differences. Because the lineages of B. pauloensis are currently exposed to distinct climatic conditions (and elevations), parapatric diversification may occur within this taxon. The eastern portion of the state of São Paulo, the most urbanized area in Brazil, represents the center of genetic diversity for B. pauloensis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0803-0) contains supplementary material, which is available to authorized users.