Substantial genetic divergence and lack of recent gene flow support cryptic speciation in a colour polymorphic bumble bee (Bombus bifarius) species complex

Phenotypic polymorphism can constitute an inherent challenge for species delimitation. This issue is exemplified in bumble bees (Bombus), where species can exhibit high colour variation across their range, but otherwise exhibit little morphological variation to distinguish them from close relatives. We examine the species status of one of the most abundant North American bumble bees, Bombus bifarius Cresson, which historically comprised two major taxa, bifarius s.s. and nearcticus. These lineages are recognized primarily by red and black variation in their mid-abdominal coloration; however, a continuum from black (nearcticus) to red (bifarius s.s.) variation has led to their historic synonymization. Integrating mitochondrial and nuclear data and whole-genome sequencing, we reveal a high level of both mitochondrial and nuclear divergence delimiting two morphologically cryptic species – the red bifarius s.s. and the colour-variable (black to red) nearcticus. Population genomic analysis supports an absence of recent genomic admixture and a strong population structure between the two clades, even in sympatry. Species distribution models predict partially differentiated niches between the genetically inferred clades with annual precipitation being a leading differentiating variable. The bifarius s.s. lineage also occupies significantly higher elevations, with regions of sympatry being among the highest elevations in nearcticus. Our data also support a subspecies-level divergence between the broadly distributed nearcticus and the island population vancouverensis. In this paper, we formally recognize the two species, Bombus bifarius Cresson and Bombus vancouverensis Cresson, the latter including the subspecies B. vancouverensis vancouverensis comb.n. and B. vancouverensis nearcticus comb.n ., with vancouverensis the name bearer due to year priority.     

Variation in composition of two bumble bee species across communities affects nectar robbing but maintains pollinator visitation rate to an alpine plant,Salvia przewalskii

1. In many flowering plants, bumble bees may forage as both pollinators and nectar robbers. This mixed foraging behaviour may be influenced by community context and consequently, potentially affect pollination of the focal plant. 2. Salvia przewalskii is both pollinated and robbed exclusively by bumble bees. In the present study area, it was legitimately visited by two species of bumble bees with different tongue length, Bombus friseanus and Bombus religiosus, but it was only robbed by Bombus friseanus, the shorter‐tongued bumble bee. The intensity of nectar robbing and pollinator visitation rate to the plant were investigated across 26 communities in the Hengduan Mountains in East Himalaya during a 2‐year project. For each of these communities, the floral diversity, and the population size and floral resource of S. przewalskii were quantified. The abundances of the two bumble bee species were also recorded. 3. Both nectar robbing and pollinator visitation rate were influenced by floral diversity. However, pollinator visitation rate was not affected by nectar robbing. The results revealed that relative abundance of the two bumble bee species significantly influenced the incidence of nectar robbing but not the pollinator visitation rate. Increased abundance of B. religiosus, the legitimate visitors, exacerbated nectar robbing, possibly by causing B. friseanus to shift to robbing; however, pollinator visitation remained at a relatively high level. 4. The results may help to explain the persistence of both nectar robbing and pollination, and suggest that, in comparison to pollination, nectar robbing is a more unstable event in a community.     

Longer tongues and swifter handling: why do more bumble bees (Bombus spp.) than honey bees (Apis mellifera) forage on lavender (Lavandula spp.)?

相似文献   

The Abundance and Pollen Foraging Behaviour of Bumble Bees in Relation to Population Size of Whortleberry (Vaccinium uliginosum)

Habitat fragmentation can have severe effects on plant pollinator interactions, for example changing the foraging behaviour of pollinators. To date, the impact of plant population size on pollen collection by pollinators has not yet been investigated. From 2008 to 2010, we monitored nine bumble bee species (Bombus campestris, Bombus hortorum s.l., Bombus hypnorum, Bombus lapidarius, Bombus pascuorum, Bombus pratorum, Bombus soroensis, Bombus terrestris s.l., Bombus vestalis s.l.) on Vaccinium uliginosum (Ericaceae) in up to nine populations in Belgium ranging in size from 80 m² to over 3.1 ha. Bumble bee abundance declined with decreasing plant population size, and especially the proportion of individuals of large bumble bee species diminished in smaller populations. The most remarkable and novel observation was that bumble bees seemed to switch foraging behaviour according to population size: while they collected both pollen and nectar in large populations, they largely neglected pollen collection in small populations. This pattern was due to large bumble bee species, which seem thus to be more likely to suffer from pollen shortages in smaller habitat fragments. Comparing pollen loads of bumble bees we found that fidelity to V. uliginosum pollen did not depend on plant population size but rather on the extent shrub cover and/or openness of the site. Bumble bees collected pollen only from three plant species (V. uliginosum, Sorbus aucuparia and Cytisus scoparius). We also did not discover any pollination limitation of V. uliginosum in small populations. We conclude that habitat fragmentation might not immediately threaten the pollination of V. uliginosum, nevertheless, it provides important nectar and pollen resources for bumble bees and declining populations of this plant could have negative effects for its pollinators. The finding that large bumble bee species abandon pollen collection when plant populations become small is of interest when considering plant and bumble bee conservation.     

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

Identifying drivers of dispersal limitation and genetic differentiation is a key goal in biogeography. We examine patterns of population connectivity and genetic diversity using restriction site‐associated DNA sequencing (RADseq) in two bumble bee species, Bombus vosnesenskii and Bombus bifarius, across latitude and altitude in mountain ranges from California, Oregon and Washington, U.S.A. Bombus vosnesenskii, which occurs across a broader elevational range at most latitudes, exhibits little population structure while B. bifarius, which occupies a relatively narrow higher elevation niche across most latitudes, exhibits much stronger population differentiation, although gene flow in both species is best explained by isolation with environmental niche resistance. A relationship between elevational habitat breadth and genetic diversity is also apparent, with B. vosnesenskii exhibiting relatively consistent levels of genetic diversity across its range, while B. bifarius has reduced genetic diversity at low latitudes, where it is restricted to high‐elevation habitat. The results of this study highlight the importance of the intersect between elevational range and habitat suitability in influencing population connectivity and suggest that future climate warming will have a fragmenting effect even on populations that are presently well connected, as they track their thermal niches upward in montane systems.     

Bumble bee species exhibit divergent responses to urbanisation in a Southern California landscape

1. Bees are ecologically important pollinators that are threatened by disease, habitat degradation, and habitat loss. Bee species that share ecological traits (e.g. body size, diet breadth, and sociality) may often respond similarly to environmental disturbance; however, few studies have examined the extent to which closely related, ecologically similar species exhibit divergent responses to the same forms of disturbance. 2. In the present study, the responses of bumble bees (Apidae: genus Bombus) to urbanization are examined by combining field surveys with the quantification of local‐ and landscape‐level variables using Geographic Information Systems (GIS). During the spring and summer of 2012 and 2013, 57, bumble bees were surveyed in 1‐ha plots in San Diego County, California that spanned a gradient of urbanisation. Species‐level, occurrence (presence/absence) data and logistic regressions were then used to identify determinants of Bombus occurrence across this landscape. 3. Approximately 1000 Bombus individuals belonging to three species were observed: B. californicus Smith, B. melanopygus Cresson, and B. vosnesenskii Radoszkowski. The occurrence of all Bombus species was negatively related to impermeable surface cover within study plots, but species responded uniquely to other plot‐scale variables and exhibited distinctive patterns of scale‐dependency with respect to impermeable surface cover surrounding plots. 4. Urbanisation negatively affected all three focal bumble bee species in this study. Species‐specific responses to plot‐scale and landscape‐scale variables presumably reflect interspecific trait differences (e.g. body size, tongue length, and foraging behaviour). Unique responses to urbanisation caution against pooling species into functional groups based merely on taxonomic relationships or perceived ecological similarities.     

Cryptic species diversity in a widespread bumble bee complex revealed using mitochondrial DNA RFLPs

Cryptic species diversity is thought to be common within the class Insecta, posing problems for basic ecological and population genetic studies and conservation management. Within the temperate bumble bee (Bombus spp.) fauna, members of the subgenus Bombus sensu stricto are amongst the most abundant and widespread. However, their species diversity is controversial due to the extreme difficulty or inability morphologically to identify the majority of individuals to species. Our character-based phylogenetic analyses of partial CO1 (700 bp) from 39 individuals spread across their sympatric European ranges provided unequivocal support for five taxa (3–22 diagnostic DNA base pair sites per species). Inclusion of 20 Irish specimens to the dataset revealed ≥2.3% sequence divergence between taxa and ≤1.3% within taxa. We developed a PCR-RFLP based method for unequivocally distinguishing amongst the four cryptic European taxa of this subgenus, B. cryptarum, B. lucorum, B. magnus and B. terrestris, and used it to analyse 391 females of the former three species collected across Ireland, all of which could be unambiguously assigned to species. Bombus lucorum was the most widely distributed and abundant of the cryptarum–lucorum–magnus species complex, comprising 56% of individuals, though it was significantly less abundant at higher altitudes (>200 m) whilst B. cryptarum was relatively more abundant at higher altitudes. Bombus magnus was rarely encountered at urban sites. Both B. lucorum and B. terrestris are nowadays reared commercially for pollination and transported globally. Our RFLP approach to identify native fauna can underpin ecological studies of these important cryptic species as well as the impact of commercial bumble bees on them.     

Forests do not limit bumble bee foraging movements in a montane meadow complex

1. Understanding the roles of habitat fragmentation and resource availability in shaping animal movement are integral for promoting species persistence and conservation. For insects such as bumble bees, their movement patterns affect the survival and reproductive potential of their colonies, as well as the pollen flow of plant species. However, the understanding of their mobility or the impact of putative barriers in natural environments is limited due to the technical difficulties of studying wild populations. 2. Genetic mark–recapture was used to estimate the foraging distance, resource use, and site connectivity of two bumble bee species in a montane meadow complex composed of open meadows within a matrix of forest. 3. There was no evidence that forests or changes in landcover function as barriers to the fine-scale movement for either species. Substantially greater colony-specific foraging distances were found for Bombus vosnesenskii (maximum: 1867 m) compared to Bombus bifarius (maximum: 362 m). Despite this difference in absolute range, both species were detected across putative forest barriers at frequencies expected by uninhibited movement. Siblings separated by greater distances were more likely to be foraging on different floral species, potentially suggesting a resource-based motivation for movement. 4. These results suggest that bumble bee foraging patterns are influenced by species-specific differences in movement capacity, with little influence of matrix composition between resource patches. They also support the perspective that habitat conservation for bumble bees should prioritise providing abundant and diverse patches of resources within species-specific movement radii with less emphasis on matrix composition.     

Invasive Bombus terrestris (Hymenoptera: Apidae) parasitized by a flagellate (Euglenozoa: Kinetoplastea) and a neogregarine (Apicomplexa: Neogregarinorida)

The flagellate Crithidia bombi and the neogregarine Apicystis bombi have been found in individuals of Bombus terrestris, a Palaearctic species of bumble bee commercially reared and shipped worldwide for pollination services. B. terrestris has recently entered into the northwestern Patagonia region of Argentina from Chile, where it was introduced in 1998. Prevalence was 21.6% for C. bombi and 3.6% for A. bombi (n = 111). The pathogens were not detected in 441 bumble bees belonging to five of the eight known Argentine native species (Bombus atratus, Bombus morio, Bombus bellicosus, Bombus opifex, Bombus tucumanus) collected elsewhere in the country. Although the absence of natural occurrence of C. bombi and A. bombi in Argentine native bumble bees cannot be ascertained at present due to the limited surveys performed, it is important to report their detection in invasive B. terrestris. The invasion event is relatively recent and the accompanying pathogens are not species specific within the genus Bombus.     

Phenotypic distribution models corroborate species distribution models: A shift in the role and prevalence of a dominant prairie grass in response to climate change

Phenotypic distribution within species can vary widely across environmental gradients but forecasts of species’ responses to environmental change often assume species respond homogenously across their ranges. We compared predictions from species and phenotype distribution models under future climate scenarios for Andropogon gerardii, a widely distributed, dominant grass found throughout the central United States. Phenotype data on aboveground biomass, height, leaf width, and chlorophyll content were obtained from 33 populations spanning a ~1000 km gradient that encompassed the majority of the species’ environmental range. Species and phenotype distribution models were trained using current climate conditions and projected to future climate scenarios. We used permutation procedures to infer the most important variable for each model. The species‐level response to climate was most sensitive to maximum temperature of the hottest month, but phenotypic variables were most sensitive to mean annual precipitation. The phenotype distribution models predict that A. gerardii could be largely functionally eliminated from where this species currently dominates, with biomass and height declining by up to ~60% and leaf width by ~20%. By the 2070s, the core area of highest suitability for A. gerardii is projected to shift up to ~700 km northeastward. Further, short‐statured phenotypes found in the present‐day short grass prairies on the western periphery of the species’ range will become favored in the current core ~800 km eastward of their current location. Combined, species and phenotype models predict this currently dominant prairie grass will decline in prevalence and stature. Thus, sourcing plant material for grassland restoration and forage should consider changes in the phenotype that will be favored under future climate conditions. Phenotype distribution models account for the role of intraspecific variation in determining responses to anticipated climate change and thereby complement predictions from species distributions models in guiding climate adaptation strategies.     

Experimental viral spillover can harm Bombus terrestris workers under field conditions

1. Deformed wing virus (DWV), notorious for its virulence in the western honey bee (Apis mellifera) when vectored by the ectoparasitic mite Varroa destructor, is also widespread among wild bumble bee species, presumably through spillover from honey bees. Experimental studies on the virulence of DWV in Bombus spp. have provided equivocal results and have until now been confined to bumble bees under laboratory conditions.
2. Here, we inoculated commercially reared Bombus terrestris workers with DWV-A through feeding or injection and introduced them into experimental colonies placed in the field, thus exposing them to the environment and associated stressors. We monitored the survival of inoculated worker bumble bees and quantified their viral load at 10 days post inoculation.
3. Bombus terrestris workers injected with DWV-A supported high viral loads and exhibited significantly reduced median survival compared to controls. Bumble bees inoculated by feeding had low or zero detectable viral loads while their mortality did not differ from the control group.
4. Our results demonstrate that, although DWV-A is pathogenic for commercial B. terrestris, the risks for individual fitness from spillover of DWV-A during foraging on shared flowers appear limited.
5. The findings of this experiment also highlight the necessity to address the potential context-dependence of virulence when evaluating the impact of a pathogen in an alternative host.

     

Unraveling the habitat preferences of two closely related bumble bee species in Eastern Europe

Co‐occurrence of closely related species is often explained through resource partitioning, where key morphological or life‐history traits evolve under strong divergent selection. In bumble bees (genus Bombus), differences in tongue lengths, nest sites, and several life‐history traits are the principal factors in resource partitioning. However, the buff‐tailed and white‐tailed bumble bee (Bombus terrestris and B. lucorum respectively) are very similar in morphology and life history, but their ranges nevertheless partly overlap, raising the question how they are ecologically divergent. What little is known about the environmental factors determining their distributions stems from studies in Central and Western Europe, but even less information is available about their distributions in Eastern Europe, where different subspecies occur. Here, we aimed to disentangle the broad habitat requirements and associated distributions of these species in Romania and Bulgaria. First, we genetically identified sampled individuals from many sites across the study area. We then not only computed species distributions based on presence‐only data, but also expanded on these models using relative abundance data. We found that B. terrestris is a more generalist species than previously thought, but that B. lucorum is restricted to forested areas with colder and wetter climates, which in our study area are primarily found at higher elevations. Both vegetation parameters such as annual mean Leaf Area Index and canopy height, as well as climatic conditions, were important in explaining their distributions. Although our models based on presence‐only data suggest a large overlap in their respective distributions, results on their relative abundance suggest that the two species replace one another across an environmental gradient correlated to elevation. The inclusion of abundance enhances our understanding of the distribution of these species, supporting the emerging recognition of the importance of abundance data in species distribution modeling.     

Improvement of almond production using Bombus terrestris (Hymenoptera: Apidae) in Mediterranean conditions

Almond trees are one of the most important crops in the Balearic Islands. The pollination of almonds is limited to the activity of insects, and cross‐pollination is necessary for fruit development. Currently, honey bees and wild bee populations are declining considerably due to multiple causes, such as the use of pesticides, diseases and habitat loss. An alternative to increase the almond production is the use of commercial pollinators. In this long‐term (3 years) study, the effect of the introduction of Bombus terrestris colonies on almond production was evaluated in two orchards. Two experimental designs were carried out to study the best management of this pollinator. For 2 years, all bumble bee colonies were placed in the middle of the plot and during the last year, the bumble bee colonies were distributed homogenously in the plot. Fruit set and the foraging behaviour of bumble bees during the blossoming period was determined, and the effect of different environmental variables on the visitation rate of bumble bees was assessed by means of a generalized linear mixed model (GLMM). Moreover, for the first time, the spatial distribution of fruit set was evaluated. Our results show that fruit set was significantly higher in the fields where B. terrestris had been introduced than in the control plots. This increased production resulted in a positive economic balance for the farmer. Moreover, bumble bees showed to prefer trees in a southwest orientation that were close to their colony. The activity of bumble bees showed to be significantly influenced by wind speed (the higher the speed the more flowers are visited by B. terrestris) and time after flowering (visitation rate decreased with days after flowering). In order to improve its management and obtain the highest possible almond production, it is important to understand the activity and behaviour of this pollinator.     

Short-term effects of experimental boreal forest &;logging disturbance on bumble bees,bumble &;bee-pollinated flowers and the bee–flower match

This study examines how, over the short term, logging affects the density of bumble bees (Apidae: Bombus), the understory plants commonly visited by bumble bees, and the numerical relationship between bumble bees and flowers. In the summers before and after winter logging, bumble bees and plants were surveyed in 50 deciduous stands (each of 8–10 ha) in the boreal forest of northern Alberta, Canada. Logging was replicated at three different intensities: 0, 10–20, and 50–75% of trees remaining. There were generally more bumble bees, species of bumble bee-visited plants, and flowers in moderately (50–75%) logged sites, but this pattern depended on the time of year. Before logging, bumble bees matched resources according to an ideal free distribution (IFD). Logging affected the distribution of bumble bees across floral resources: the slope of the regression relating bumble bee and flower proportions was less than one for clearcut and control treatments (i.e., undermatching), with too many bumble bees in the flower-poor compartments and too few in the flower-rich ones. Deviations from an IFD were negative in control sites, such that fewer bumble bees occurred here than warranted by flower numbers. Controlling for flower density, bumble bee density was significantly greater in clearcuts than in the other treatments. By disproportionately visiting plants in clearcuts (relative to flower density), and by undermatching, bumble bees in clearcuts should experience higher levels of competition. Conversely, the fewer (and undermatching) bumble bees in control sites (relative to flower abundances there) may cause these plants to obtain diminished pollination service. The proximity of clearcut logging to pristine areas may therefore negatively impact plants and bumble bees in the pristine areas, at least in the season immediately following logging.     

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

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

Effective paternity in natural colonies of Japanese native bumble bees

Kin selection theory predicts a high coefficient of genetic relatedness among nestmates, explaining the frequent evolution of eusociality in a social hymenopteran colony. The bumble bee is a primitively eusocial hymenoptera whose colonies are founded by a single queen. In such a monogynous colony, mating frequency of the queen is the sole factor that affects genetic relatedness among nestmates. Although the queens of most bumble bee species are known to be monandrous, there are some species that are known to be polyandrous. Here, we estimated the effective paternity in the native colonies of Japanese bumble bees, Bombus ardens, B. diversus, and B. honshuensis, using genetic polymorphic data of microsatellites. We found no evidence for polyandry in any investigated colonies, which suggests that within-colony genetic relatedness is very high for all of these colonies.     

BUMBLE BEE POLLINATION AND FLORAL MORPHOLOGY: FACTORS INFLUENCING POLLEN DISPERSAL IN THE ALPINE SKY PILOT,POLEMONIUM VISCOSUM (POLEMONIACEAE)

Pollen dispersal success in entomophilous plants is influenced by the amount of pollen produced per flower, the fraction of pollen that is exported to other flowers during a pollinator visit, visitation frequency, and the complementarity between pollen donor and recipients. For bumble bee-pollinated Polemonium viscosum the first three determinants of male function are correlated with morphometric floral traits. Pollen production is positively related to corolla and style length, whereas pollen removal per visit by bumble bee pollinators is a positive function of corolla flare. Larger-flowered plants receive more bumble bee visits than small-flowered individuals. We found no evidence of tradeoffs between pollen export efficiency and per visit accumulation of outcross pollen; each was influenced by unique aspects of flower morphology. Individual queen bumble bees of the principal pollinator species, Bombus kirbyellus, were similar in male, female, and absolute measures of pollination effectiveness. An estimated 2.9% of the pollen that bumble bees removed from flowers during a foraging bout was, on average, deposited on stigmas of compatible recipients. Significant plant-to-plant differences in pollen production, pollen export per visit, and outcross pollen receipt were found for co-occurring individuals of P. viscosum indicating that variation in these fitness related traits can be seen by pollinator-mediated selection.     

Attractiveness of single and multiple species flower patches to beneficial insects in agroecosystems

The provision of floral resources for the enhancement of beneficial insect populations has shown promise as a strategy to enhance biological control and pollination in agroecosystems. One approach involves the provision of a single flower species while a second involves the multiple flower species, but the two have never been compared experimentally. Here we examine the influence of single and multiple species flower treatments on the abundance and foraging behaviour of key beneficial insects in two agricultural agroecosystems (broccoli and lucerne crops). The five flower treatments comprised buckwheat only, phacelia only, a simple mixture of buckwheat and phacelia, a complex mixture of buckwheat, phacelia and a commercial seed blend or the existing crop as a control. The abundance of bumble‐bees (Bombus hortorum) and honey bees (Apis mellifera) was highest in the three treatments that contained phacelia, while hoverfly (Melanostoma fasciatum) numbers were high in all four flower treatments. Bumble‐bees and honey bees probed almost exclusively phacelia flowers, even when provided with a choice of other flower species in the simple and complex mixture treatments. In contrast, hoverflies probed the flowers of all plant species in single and multiple species treatments, with no apparent difference in acceptance. However, in mixture treatments, the majority of individual bumble‐bees, honey bees and hoverflies probed the flowers from only one species, despite the presence of alternative flower species. Our results illustrate how an appreciation of insect floral attractiveness can be used to customise the species composition of floral patches to potentially maximise biological control and pollination in targeted agroecosystems.     

Allozyme variation in bumble bees (Hymenoptera: Apidae)

Allozyme variation at an average of 37.3 loci was assessed in queens of 16 Bombus and 2 Psithyrus bumble bee species from North America. The mean expected heterozygosity (H) for the Bombus species was 0.008 +/- 0.006 (95% confidence limits) and that for the Psithyrus was 0.007 +/- 0.007. These levels are significantly lower than found in other Hymenoptera but are comparable to those found in previous studies of bumble bees based on far fewer loci. Neutral mutation and random genetic drift can account for the observed variation, but this implies a very small effective population size for species of bumble bees.