Cryptic Bumblebee Species: Consequences for Conservation and the Trade in Greenhouse Pollinators

Abstract

Commercial greenhouse growers in both Japan and China are increasingly using reared orange-tailed bumblebees known previously as Bombus hypocrita Pérez as pollinators. Phylogenetic analysis of the DNA (COI) barcodes with Bayesian methods shows that this “species” is a long-standing confusion of two cryptic species. We find that the orange-tailed bumblebees in North China are actually part of the widespread Russian (otherwise white-tailed) B. patagiatus Nylander (as B. patagiatus ganjsuensis Skorikov, n. comb.), whereas the orange-tailed bees in Japan are true B. hypocrita. This situation has been further complicated because two other cryptic species from North China that were previously confused with the Russian B. patagiatus are now recognised as separate: B. lantschouensis Vogt n. stat. and B. minshanensis Bischoff n. stat.. As demand for pollination services by greenhouse growers inevitably increases, these bees are more likely to be transported between countries. In order to conserve genetic resources of pollinator species for their option value for future food security, we advocate preventing trade and movement of B. patagiatus from China into Japan and of B. hypocrita from Japan into China.

About the Authors:

Paul H. Williams

* E-mail: [email protected]

Affiliation: Department of Entomology, The Natural History Museum, London, United Kingdom

Jiandong An

Affiliation: Institute of Apiculture, Chinese Academy of Agricultural Sciences, Beijing, China

Mark J. F. Brown

Affiliations School of Biological Sciences, Royal Holloway University of London, Egham, Surrey, United Kingdom, Department of Zoology, School of Natural Sciences, Trinity College Dublin, Dublin, Ireland

James C. Carolan

Affiliations Department of Zoology, School of Natural Sciences, Trinity College Dublin, Dublin, Ireland, Department of Biology, National University of Ireland Maynooth, Maynooth, County Kildare, Ireland

Dave Goulson

Affiliation: Biological & Environmental Sciences, University of Stirling, Stirling, United Kingdom

Jiaxing Huang

Affiliation: Institute of Apiculture, Chinese Academy of Agricultural Sciences, Beijing, China

Masao Ito

Affiliation: Higashi-Ku, Sapporo, Japan

Introduction

Commercial movement of species outside of their natural ranges has had significant negative impacts on biodiversity [1]–[3]. Consequently, such translocation is now regulated at national and international levels [4]. To be effective, regulation depends upon accurate taxonomy to identify the relevant species. Here we show that imprecise taxonomy poses a substantial threat to ecologically important pollinators.

Bumblebees are among the most important pollinators in wild ecosystems, but recently have been suffering worldwide declines [5], [6]. At the same time, bumblebees have become increasingly important commercially for their pollination services to agriculture, especially for tree fruits, berries, and greenhouse crops such as tomatoes [7]. This has led to the widespread movement of bumblebees between countries to provide pollination services, an industry now worth billions of dollars annually [8]–[10]. Such commercial translocation has resulted in the introduction and invasion of exotic bumblebee species (and their pathogens) into New Zealand [11], Tasmania [12], South America [13], [14], and Japan [15], [16].

The situation in Japan has been particularly well-studied. There, Bombus terrestris (Linnaeus) was introduced from Europe into greenhouses, but feral colonies were soon discovered and the species has shown not only invasive spread within Japan, but is also replacing the indigenous B. hypocrita Pérez in many areas [15], [16]. One solution is the commercial development of indigenous pollinators and in Japan, B. hypocrita is undergoing trials [17]. In China, government-funded research projects have been established to study the feasibility of what has been understood to be the same species [18]–[20]. B. hypocrita is currently recognised as an orange-tailed species, believed to be distributed in both Japan and China [18], [21]–[27]. Alongside these bees in China, some other bumblebees with white tails are also being used as greenhouse pollinators [28] and many of these have been widely referred to previously using the name B. patagiatus Nylander [18], [22], [25], [26], [29]. But just as B. hypocrita in Japan is currently seriously threatened by introduced B. terrestris (Linnaeus), so populations of all of these bumblebees in Japan and China could be threatened by introductions between countries if the bumblebees used are in fact not conspecific.

All of these commercially important species belong to the subgenus Bombus s. str. and many of these species are well known for being cryptic in Europe [30]. Not all individuals can be identified with confidence using morphological characters [31] and specialists also disagree among themselves on the precise criteria for diagnosing them [32]. Nonetheless, support for the interpretation that there are separate species continues to grow from studies of morphology [33], enzyme electrophoresis [34], [35], male labial gland secretions [36], [37], and DNA sequences [38]–[41]. However, this European work has been done against a background of very patchy knowledge of the Asian species of the group.

In this paper we use DNA barcodes to show that the orange-tailed bumblebees previously recognised as B. hypocrita in Asia are actually comprised of parts of two more geographically restricted species: B. hypocrita, and a part of B. patagiatus with an unrecognised cryptic colour pattern. Given that these bees are already being used for pollination in greenhouses in Asia, we discuss the consequences of our results for conserving the genetic resources of these commercially valuable pollinators and the need to restrict movement of bumblebees between China and Japan.

Materials and Methods

Sampling bees

We sampled bumblebees as part of a review of all of the species of the subgenus Bombus s. str. across their entire global distributions [42], which encompass most of the northern hemisphere. Progress with the taxonomy of this group using only morphological evidence has been difficult [31], but recently new insights have been gained by using DNA-sequence data [38], [40]. In insects, sequences of the mitochondrial COI (cox-1) gene have been used to discover species identical to those recognised by more traditional methods [43]. Inheritance of mitochondrial genes is reasonably well understood [44]. It avoids the problems with nuclear genes that arise from having multiple alleles, while the absence of indels makes alignment for homology straightforward. Mitochondrial genes also have a relatively high substitution rate, so that even the short COI ‘barcode’ region of 658 nucleotides can be used to distinguish the most closely-related taxa [45], [46], and COI-barcode sequences as short as 100 nucleotides can be diagnostic for 90% of the species in other animal groups [47]. Consequently, despite many potential pitfalls and the desirability of external supporting evidence, COI barcodes can be useful for recognising likely cryptic species [48]. Permissions for collection and export of bumblebee specimens from China were obtained via the CAAS Institute of Apiculture and CAS Institute of Zoology, number 2010-86 issued 15.ix.2010.

DNA-barcode data

Most specimens were extracted, amplified, and sequenced at the Biodiversity Institute of Ontario, University of Guelph, as part of the BEE-BOL campaign to barcode the bees of the world [49]. COI-barcode extraction, amplification, and sequencing used the standard protocols described by Hebert et al. [45]. Universal primers for the COI-barcode sequence for insects were used (variants LepF1 and LepR1 [50]). COI-barcode sequences (without primer sequences) from the samples were aligned using the ClustalW function within BioEdit (version 7.0.9.0; www.mbio.ncsu.edu/BioEdit/bioedit.html, accessed 2010) and trimmed to a common frame length of 658 nucleotides that is shared by most samples. Overall, we analysed 559 sequences from 279 localities in 33 countries from throughout the world-wide range of the subgenus. For the analysis reported here, sequences for species that are the focus of this paper (i.e. B. hypocrita, B. patagiatus, and cryptic relatives: 106 sequences) were reduced to 20 unique haplotypes using Collapse (version 1.2; darwin.uvigo.es/software/collapse.html, accessed 2011) after ranking sequences by their length. For the species that are not the focus of this analysis, we used sequences from among the samples available that were collected closest to the species' type localities (details of sequence data are given in Appendix S1).

Phylogenetic analysis

Phylogenetic relationships among sequences were estimated using MrBayes (version 3.1.2 [51], [52]), from 10 million generations of the metropolis-coupled Markov-chain Monte Carlo (MCMC) algorithm with four chains, chain temperature set to 0.2, and with sampling of the trees every 1000 generations. We found the nucleotide-substitution model that fitted our COI-barcode data best according to jModelTest (version 0.1.1 [53]) to be the general time-reversible model with a gamma frequency distribution of changes among sites and allowing invariant sites (GTR+G+I). Burn-in was set initially to 10%, with convergence between two separate runs of the analysis judged to have occurred when the average standard deviation of the split frequencies approached stationarity. Post burn-in stability of the log likelihood of the cold chain was confirmed using Tracer (version 1.5.0 [54]) and stability of the sample groups was confirmed using AWTY (version 0.8.1 [55]). The sample of 16,002 post burn-in trees from both replicates was combined. Trees were rooted using data for the sister groups of Bombus s. str. as sequences from B. (Pyrobombus) vagans Smith, B. (Alpinobombus) alpinus (Linnaeus), and B. (Al.) balteatus Dahlbom following the results of Cameron et al. [56] from a phylogenetic analysis of five genes across almost all bumblebee species. MEGA (version 4.0 [57]) was used to generate intra- and inter-group genetic distances and to determine the number of informative polymorphisms that support the phylogenetic grouping of the cryptic taxa and the recognition of new species. Genetic distances (sequence divergences) were calculated using the Kimura two-parameter (K2P) distance model [58]. The translated amino acid sequences were also compared to determine whether polymorphisms were non-synonymous. Treatment of taxon names follows the rules laid down by ICZN [59].

Results

Our estimate of phylogenetic relationships from COI sequences (Figure 1) shows only weak support (posterior probabilities) for several of the species groups, including especially those representing the relationships of B. hypocrita and B. magnus Vogt with other species. However, this tree shows strong support for the monophyly of the groups interpreted here as the species B. hypocrita, B. lantschouensis Vogt, B. minshanensis Bischoff, and B. patagiatus. It does not support a monophyletic group of orange-tailed bees (o) among the white tailed bees (w) within our focal group of B. hypocrita and B. patagiatus in the former broader sense. In contrast, among these bees it supports: (1) that there are two relatively distantly related groups of Asian orange-tailed bees - one group in Japan, Korea, Primorsky, and Sakhalin (B. hypocrita) and another group in North China (B. patagiatus ganjsuensis Skorikov); and (2) that the white-tailed bees considered previously to be a single species, B. patagiatus, are comprised of at least three distinct groups (B. lantschouensis, B. minshanensis, and B. patagiatus).

[Figure omitted. See PDF.]

Figure 1. Estimated Phylogenetic Tree for Asian Bumblebees Resembling B. hypocrita and B. patagiatus.

Consensus Bayesian tree for (o) orange- and (w) white-tailed examples representing the unique haplotypes from among 106 COI-barcode sample sequences of B. hypocrita, B. lantschouensis, B. minshanensis and B. patagiatus, together with examples representing all of the other species of Bombus s. str. (from a total of 559 sequences). The codes following the taxon names are the specimen identifiers from the project database and from BOLD (or from other external databases including GenBank), followed by country and province. Values next to the nodes are posterior probabilities for groups (groups with values of <0.9 are considered unreliable) and the scale bar represents 0.08 substitutions per nucleotide site.

https://doi.org/10.1371/journal.pone.0032992.g001

The recognition of B. hypocrita, B. lantschouensis, B. minshanensis, and B. patagiatus as separate species is supported by a number of unique and diagnostic polymorphisms (Table 1). Fifteen, four, and five unique polymorphisms are present in COI barcodes for each of B. hypocrita (8 haplotypes), B. lantschouensis (3 haplotypes), and B. minshanensis (3 haplotypes), including non-synonymous polymorphisms. Within-group K2P distances for these species are low, ranging from 0.002 to 0.004 (Table 2). B. patagiatus is characterised by seven unique polymorphisms, although a clear divergence exists within this group, as indicated by the high intra-group distance value of 0.014 (synonymous polymorphisms). Two subgroups are recognised: Group 1, of haplotypes from white-tailed bumblebees from Russia and the far north east of China; and Group 2, of haplotypes from both white-tailed and orange-tailed bumblebees from North China. Analyses of these groups separately showed intra-group distances of 0.002 (Group 1) and <0.001 (Group2). The latter value equates to just one polymorphism among the 42 specimens and three haplotypes of B. patagiatus found in this study. Eleven unique polymorphisms are shared by the COI barcodes for B. patagiatus Group 2 (Table 1), providing strong support for the cryptic status of the orange-tailed B. patagiatus ganjsuensis and the placement of these orange-tailed bumblebees as part of B. patagiatus and not B. hypocrita.

[Figure omitted. See PDF.]

Table 1. COI-Barcode Polymorphisms among Haplotypes of B. patagiatus and B. hypocrita.

https://doi.org/10.1371/journal.pone.0032992.t001

[Figure omitted. See PDF.]

Table 2. COI-Barcode Characteristics for B. hypocrita, B. lantschouensis, B. minshanensis, and B. patagiatus.

https://doi.org/10.1371/journal.pone.0032992.t002

Discussion

Taxonomic revisions of bumblebees and other commercially important insect groups are needed in part to prevent people from inadvertently transporting cryptic non-native species. Such revisions crucially need to be global in scope, both to include all of the known species as well as to include broadly representative samples from across the entire breadths of all of the species' geographic ranges. Attempts to revise groups from more restricted geographic sampling have led to problems. For example, the name B. minshanicola Bischoff has recently been used [60] for a ‘neglected taxon’ of Bombus s. str. from China. This name does not appear in our Figure 1 because we agree from a global analysis [42] with earlier studies [61] (see also [33]) that recognised B. minshanicola from Gansu as part of the broader species B. longipennis Friese, which extends further west within China and the Himalaya. While we also note that there are many potential pitfalls when using COI barcodes for exploring relationships among taxa [62], there is as yet no evidence that these have affected studies of bumblebees [60], [63].

Disambiguation of species

Figure 2 uses our sequenced samples to map the former concept of B. hypocrita, as revised by Tkalcu [21] (map on his page 89) as a taxon of orange-tailed bumblebees, distributed from Gansu to Japan. From our COI-based estimate of relationships (Figure 1), we recognise two species within Tkalcu's [21] broader concept in Fig. 2: B. hypocrita in Japan and adjacent Pacific provinces of Asia (syntypes examined, Muséum National d'Histoire Naturelle, Paris), and a second species in North China. For the latter Chinese bumblebees, the oldest available name is Bombus ikonnikovi ganjsuensis Skorikov [64] (holotype examined, Zoological Institute of the Russian Academy of Sciences, St Petersburg), revised here to Bombus patagiatus ganjsuensis Skorikov (new combination).

[Figure omitted. See PDF.]

Figure 2. Orange-Tailed B. hypocrita Sensu Tkalcu Recognised as B. patagiatus ganjsuensis and B. hypocrita.

Distribution of the orange-tailed taxon B. hypocrita in the sense of Tkalcu [21], revised from the results in Figure 1 to two species: black spots, B. patagiatus ganjsuensis; and black triangles, B. hypocrita. Records shown include only specimens identified in our data from COI barcodes. Map showing shaded relief and international boundaries, Cartesian orthonormal projection.

https://doi.org/10.1371/journal.pone.0032992.g002

Similarly, Figure 3 maps the former concept of B. patagiatus, as revised by Tkalcu [29] (figure on his page 52) as a taxon of white-tailed bumblebees, distributed from near the Finnish border of Russia in the west to Sachalin in the east, and then south westwards through North China to Gansu and Sichuan. From our COI-based estimate of relationships (Figure 1), we recognise three species within Tkalcu's [29] concept in Fig. 3: B. patagiatus patagiatus s. str. (type presumed lost [29], although the taxon concept is not in doubt), B. lantschouensis (new status, syntypes examined, Zoölogisch Museum der Universiteit van Amsterdam), and B. minshanensis (new status, syntypes examined, Museum für Naturkunde an der Humboldt-Universität, Berlin). B. lantschouensis and B. minshanensis are reciprocally monophyletic with non-synonymous polymorphisms and no individuals with intermediate colour patterns are known to us, so we recognise them as separate species for the first time. Assessment of evidence from DNA barcodes has uncovered multiple cryptic species in other insects [45], [48], suggesting that this situation could be widespread among insects.

[Figure omitted. See PDF.]

Figure 3. White-Tailed B. patagiatus Sensu Tkalcu Recognised as B. patagiatus patagiatus, B. lantschouensis, and B. minshanensis.

Asian distribution of the white-tailed taxon B. patagiatus in the sense of Tkalcu [29], revised from the results in Figure 1 to three species: white spots, B. patagiatus patagiatus; grey spots, B. lantschouensis; and white triangles, B. minshanensis. Records shown include only specimens identified in our data from COI barcodes. Map showing shaded relief and international boundaries, Cartesian orthonormal projection.

https://doi.org/10.1371/journal.pone.0032992.g003

In consequence, our revised concept of B. patagiatus (Figure 4) reunites white-tailed bees (open circles) from the Taiga forests of Russia to the north of the Mongolian desert, with orange-tailed bees (black spots) previously misidentified as B. hypocrita, from the hills of North China on the southern side of the Mongolian desert. During specially-targeted fieldwork in 2010, specimens were also collected from an isolated patch of suitable habitat on a mountain top (near Huanggangliang) in between these two regions, within the barrier zone formed by the Mongolian arid belt (in Neimenggu or Inner Mongolia). These individuals (our Group 2) all have the COI haplotype of B. patagiatus ganjsuensis, but have either varying degrees of intermediate colour patterns with narrow dorsal pale bands (like B. patagiatus ganjsuensis) and white tails (like B. patagiatus patagiatus), or even have broad dorsal pale bands and white tails resembling B. patagiatus patagiatus (our Group 1). Therefore, while the orange-tailed colour pattern appears to be restricted to the south of the Mongolian arid zone, closely related haplotypes (Group 2) and individuals with intermediate colour patterns extend into suitable habitats within the arid zone. We recognise them all here as parts of one species for the first time.

[Figure omitted. See PDF.]

Figure 4. Asian Distribution of the Revised Concept of B. patagiatus.

These bees are revised from the results in Figure 1 to two subspecies: white spots, white-tailed B. patagiatus patagiatus; and black spots, orange-tailed B. patagiatus ganjsuensis. Records shown include only specimens identified in our data from COI barcodes. Map showing shaded relief and international boundaries, Cartesian orthonormal projection.

https://doi.org/10.1371/journal.pone.0032992.g004

Movement of cryptic commercial pollinators between countries

Our results are significant because conserving the genetic resources of commercially important pollinator species for their option value for future food security in a changing world means that we need to consider the potential consequences of moving bumblebees and their pathogens outside their indigenous ranges and into new contact with other species. It is already known that introductions of B. terrestris into Japan have resulted in the spread of feral populations, the growth of which has coincided on a local scale with declines in populations of indigenous B. hypocrita [15], [16]. Therefore with due regard to the precautionary principle [5], [65], we advocate preventing any movement of B. patagiatus from China into Japan (where there are no indigenous populations of B. patagiatus) and of B. hypocrita from Japan into China (where B. hypocrita is, as far as we know, absent) until it can be proven safe. We also recommend research to assess the potential of the indigenous species for use in commercial pollination in greenhouses within each region. Deliberate movement of bumblebee species between countries will otherwise become more likely as increasing demand for food puts pressure on commercial growers to deliver reliable pollination in ever larger areas of greenhouse crops. Critically, minimising the risk from commercial movement of species and their pathogens (to which different bumblebee populations may have differing resistance, even within a species) depends upon a full understanding of their taxonomy.

Supporting Information

[Figure omitted. See PDF.]

Appendix S1.

Accession numbers for sequence data for the samples used in Fig. 1, including IDs from the BOLD database.

https://doi.org/10.1371/journal.pone.0032992.s001

(DOC)

Acknowledgments

We thank A. Byvaltsev, B. Cederberg, R. Dawson, F. Gusenleitner, G. Hancock, R. Hart, W. Hogenes, F. Koch, M. Proshchalykin, H. Schmidt, C. Sheffield, D. Sheppard, Z. Xie for specimens; the Biodiversity Institute of Ontario, Guelph, for extraction and sequencing; L. Bailey and anonymous referees for comments.

Author Contributions

Conceived and designed the experiments: PW. Performed the experiments: PW JC JH. Analyzed the data: PW JC. Contributed reagents/materials/analysis tools: JA MB JC DG JH MI. Wrote the paper: PW MB. Co-wrote grant application: DG. Obtained information for manuscript on Japanese species: MI.

Citation: Williams PH, An J, Brown MJF, Carolan JC, Goulson D, Huang J, et al. (2012) Cryptic Bumblebee Species: Consequences for Conservation and the Trade in Greenhouse Pollinators. PLoS ONE7(3): e32992. https://doi.org/10.1371/journal.pone.0032992

References

1. Mooney HA, Cleland EE (2001) The evolutionary impact of invasive species. Proceedings of the National Academy of Sciences USA 98: 5446–5451.HA MooneyEE Cleland2001The evolutionary impact of invasive species.Proceedings of the National Academy of Sciences USA9854465451

2. Strayer DL, Eviner VT, Jeschke JM, Pace ML (2006) Understanding the long-term effects of species invasions. Trends in Ecology and Evolution 21: 645–651.DL StrayerVT EvinerJM JeschkeML Pace2006Understanding the long-term effects of species invasions.Trends in Ecology and Evolution21645651

3. Molnar JL, Gamboa RL, Revenga C, Spalding MD (2008) Assessing the global threat of invasive species to marine biodiversity. Frontiers in Ecology and the Environment 6: 485–492.JL MolnarRL GamboaC. RevengaMD Spalding2008Assessing the global threat of invasive species to marine biodiversity.Frontiers in Ecology and the Environment6485492

4. McNeely JA, Mooney HA, Neville LE, Schei PJ, Waage JK, editors. (in press) JA McNeelyHA MooneyLE NevillePJ ScheiJK Waagein press(in press) A Global Strategy on Non-Native Species: IUCN, Switzerland and Cambridge, UK. (in press) A Global Strategy on Non-Native Species: IUCN, Switzerland and Cambridge, UK.

5. Williams PH, Osborne JL (2009) Bumblebee vulnerability and conservation world-wide. Apidologie 40: 367–387.PH WilliamsJL Osborne2009Bumblebee vulnerability and conservation world-wide.Apidologie40367387

6. Cameron SA, Lozier JD, Strange JP, Koch JB, Cordes N, et al. (2011) Patterns of widespread decline in North American bumble bees. Proceedings of the National Academy of Sciences of the United States of America 108: 662–667.SA CameronJD LozierJP StrangeJB KochN. Cordes2011Patterns of widespread decline in North American bumble bees.Proceedings of the National Academy of Sciences of the United States of America108662667

7. Velthuis HHW, Doorn Av (2006) A century of advances in bumblebee domestication and the economic and environmental aspects of its commercialization for pollination. Apidologie 37: 421–451.HHW VelthuisAv Doorn2006A century of advances in bumblebee domestication and the economic and environmental aspects of its commercialization for pollination.Apidologie37421451

8. Dias BSF, Raw A, Imperatri-Fonseca VL (1999) International Pollinators Initiative: the São Paulo declaration on pollinators. Report on the recommendations of the workshop on the conservation and sustainable use of pollinators in agriculture with emphasis on bees. Brasília: Brazilian Ministry of the Environment. 79 p.BSF DiasA. RawVL Imperatri-Fonseca1999International Pollinators Initiative: the São Paulo declaration on pollinators. Report on the recommendations of the workshop on the conservation and sustainable use of pollinators in agriculture with emphasis on beesBrasíliaBrazilian Ministry of the Environment79

9. D. Goulson 2010 Bumblebees, behaviour, ecology, and conservation Oxford Oxford University Press x+317

10. Winter K, Adams L, Thorp RW, Inouye DW, Day L, et al. (2006) Importation of non-native bumble bees into North America: potential consequences of using Bombus terrestris and other non-native bumble bees for greenhouse crop pollination in Canada, Mexico, and the United States. San Francisco: 33 p.K. WinterL. AdamsRW ThorpDW InouyeL. Day2006Importation of non-native bumble bees into North America: potential consequences of using Bombus terrestris and other non-native bumble bees for greenhouse crop pollination in Canada, Mexico, and the United StatesSan Francisco33

11. Macfarlane RP, Gurr L (1995) Distribution of bumble bees in New Zealand. New Zealand Entomologist 18: 29–36.RP MacfarlaneL. Gurr1995Distribution of bumble bees in New Zealand.New Zealand Entomologist182936

12. Stout J-C, Goulson D (2000) Bumble bees in Tasmania: their distribution and potential impact on Australian flora and fauna. Bee-World 81: 80–86.J-C StoutD. Goulson2000Bumble bees in Tasmania: their distribution and potential impact on Australian flora and fauna.Bee-World818086

13. Arretz PV, Macfarlane RP (1982) The introduction of Bombus ruderatus to Chile for red clover pollination. Bee World 67: 15–22.PV ArretzRP Macfarlane1982The introduction of Bombus ruderatus to Chile for red clover pollination.Bee World671522

14. Torretta JP, Medan D, Abrahamovich AH (2006) First record of the invasive bumblebee Bombus terrestris (L.) (Hymenoptera, Apidae) in Argentina. Transactions of the American Entomological Society 132: 285–289.JP TorrettaD. MedanAH Abrahamovich2006First record of the invasive bumblebee Bombus terrestris (L.) (Hymenoptera, Apidae) in Argentina.Transactions of the American Entomological Society132285289

15. Inoue MN, Yokoyama J, Washitani I (2008) Displacement of Japanese native bumblebees by the recently introduced Bombus terrestris (L.) (Hymenoptera: Apidae). Journal of Insect Conservation 12: 135–146.MN InoueJ. YokoyamaI. Washitani2008Displacement of Japanese native bumblebees by the recently introduced Bombus terrestris (L.) (Hymenoptera: Apidae).Journal of Insect Conservation12135146

16. Matsumura C, Yokoyama J, Washitani I (2004) Invasion status and potential impacts of an invasive alien bumblebee, Bombus terrestris L. (Hymenoptera: Apidae) naturalized in southern Hokkaido, Japan. Global Environment Research 8: 51–66.C. MatsumuraJ. YokoyamaI. Washitani2004Invasion status and potential impacts of an invasive alien bumblebee, Bombus terrestris L. (Hymenoptera: Apidae) naturalized in southern Hokkaido, Japan.Global Environment Research85166

17. Asada S, Ono M (1996) Crop pollination by Japanese bumblebees, Bombus spp. (Hymenoptera: Apidae): tomato foraging behaviour and pollination efficiency Applied Entomology and Zoology 31: 581–586.S. AsadaM. Ono1996Crop pollination by Japanese bumblebees, Bombus spp.(Hymenoptera: Apidae): tomato foraging behaviour and pollination efficiency Applied Entomology and Zoology31581586

18. An J-D, Huang J-X, Williams PH, Wu J, Zhou B-F (2010) Species diversity and colony characteristics of bumblebees in the Hebei region of North China. Chinese Journal of Applied Ecology 21: 1542–1550.J-D AnJ-X HuangPH WilliamsJ. WuB-F Zhou2010Species diversity and colony characteristics of bumblebees in the Hebei region of North China.Chinese Journal of Applied Ecology2115421550

19. Zhao Y-Z, An J-D, Zhou Z-Y, Dong J, Xing Y-H, et al. (2011) Pollination behavior of Apis mellifera ligustica and Bombus hypocrita (Hymenoptera, Apidae) and the influencing factors in peach greenhouse. Acta Entomologica Sinica 54: 89–96.Y-Z ZhaoJ-D AnZ-Y ZhouJ. DongY-H Xing2011Pollination behavior of Apis mellifera ligustica and Bombus hypocrita (Hymenoptera, Apidae) and the influencing factors in peach greenhouse.Acta Entomologica Sinica548996

20. Chen W-F, An J-D, Dong J, Ding K-F, Gao S (2011) Flower-visiting behavior and pollination ecology of different bee species on greenhouse strawberry. Chinese Journal of Ecology 30: 290–296.W-F ChenJ-D AnJ. DongK-F DingS. Gao2011Flower-visiting behavior and pollination ecology of different bee species on greenhouse strawberry.Chinese Journal of Ecology30290296

21. Tkalcu B (1962) Contribution à l'étude des bourdons du Japon (I) (Hymenoptera, Apoidea). Bulletin de la Société entomologique de Mulhouse 1962: 81–100.B. Tkalcu1962Contribution à l'étude des bourdons du Japon (I) (Hymenoptera, Apoidea).Bulletin de la Société entomologique de Mulhouse196281100

22. Yao J, Wang S-F (2005) Hymenoptera, Apidae, Bombini. In: Yang X, editor. Insect fauna of middle-west Qiling range and south mountains of Gansu province. Beijing: Science Press. pp. 890–896.J. YaoS-F Wang2005Hymenoptera, Apidae, Bombini.X. YangInsect fauna of middle-west Qiling range and south mountains of Gansu provinceBeijingScience Press890896

23. Krüger E (1958) Phaenoanalytische Studien an einigen Arten der Untergattung Terrestribombus O. Vogt (Hymenoptera, Bombiidae). III. Teil. Tijdschrift voor Entomologie 101: 283–344.E. Krüger1958Phaenoanalytische Studien an einigen Arten der Untergattung Terrestribombus O. Vogt (Hymenoptera, Bombiidae). III. Teil.Tijdschrift voor Entomologie101283344

24. Sakagami SF, Ishikawa R (1969) Note préliminaire sur la répartition géographique des bourdons japonais, avec descriptions et remarques sur quelques formes nouvelles ou peu connues. Journal of the Faculty of Science, Hokkaido University (Zoology) 17: 152–196.SF SakagamiR. Ishikawa1969Note préliminaire sur la répartition géographique des bourdons japonais, avec descriptions et remarques sur quelques formes nouvelles ou peu connues.Journal of the Faculty of Science, Hokkaido University (Zoology)17152196

25. Bischoff H (1936) Schwedisch-chinesische wissenschaftliche Expedition nach den nordwestlichen Provinzen Chinas, unter Leitung von Dr. Sven Hedin und Prof. Sü Ping-chang. Insekten gesammelt vom schwedischen Arzt der Expedition Dr. David Hummel 1927–1930. 56. Hymenoptera. 10. Bombinae. Arkiv för zoologi 27A: 1–27.H. Bischoff1936Schwedisch-chinesische wissenschaftliche Expedition nach den nordwestlichen Provinzen Chinas, unter Leitung von Dr. Sven Hedin und Prof. Sü Ping-chang. Insekten gesammelt vom schwedischen Arzt der Expedition Dr. David Hummel 1927–1930. 56. Hymenoptera. 10. Bombinae.Arkiv för zoologi27A127

26. Williams PH (1991) The bumble bees of the Kashmir Himalaya (Hymenoptera: Apidae, Bombini). Bulletin of the British Museum (Natural History) (Entomology) 60: 1–204.PH Williams1991The bumble bees of the Kashmir Himalaya (Hymenoptera: Apidae, Bombini).Bulletin of the British Museum (Natural History) (Entomology)601204

27. Ito M, Sakagami SF (1980) The bumblebee fauna of the Kurile Islands (Hymenoptera: Apidae). Low Temperature Science (B) 38: 23–51.M. ItoSF Sakagami1980The bumblebee fauna of the Kurile Islands (Hymenoptera: Apidae).Low Temperature Science (B)382351

28. Peng W-J, Huang J-X, Wu J, An J-D (2009) Geographic distribution and bionomics of six bumblebee species in North China. Chinese Bulletin of Entomology 46: 115–120.W-J PengJ-X HuangJ. WuJ-D An2009Geographic distribution and bionomics of six bumblebee species in North China.Chinese Bulletin of Entomology46115120

29. Tkalcu B (1967) Sur deux espèces de bourdons décrites par William Nylander (Hymenoptera, Apoidea: Bombus). Bulletin de la Société entomologique de Mulhouse 1967: 41–58.B. Tkalcu1967Sur deux espèces de bourdons décrites par William Nylander (Hymenoptera, Apoidea: Bombus).Bulletin de la Société entomologique de Mulhouse19674158

30. Carolan JC, Murray TE, Fitzpatrick U, Crossley J, Schmidt H, et al. (2012) Colour patterns do not diagnose species: quantitative evaluation of a DNA barcoded cryptic bumblebee complex. PLoS ONE 7: 662–667.JC CarolanTE MurrayU. FitzpatrickJ. CrossleyH. Schmidt2012Colour patterns do not diagnose species: quantitative evaluation of a DNA barcoded cryptic bumblebee complex.PLoS ONE7662667

31. Rasmont P (1984) Les bourdons du genre Bombus Latreille sensu stricto en Europe occidentale et centrale (Hymenoptera, Apidae). Spixiana 7: 135–160.P. Rasmont1984Les bourdons du genre Bombus Latreille sensu stricto en Europe occidentale et centrale (Hymenoptera, Apidae).Spixiana7135160

32. Williams PH (2000) Are Bombus lucorum and magnus separate species? BWARS Newsletter 2000: 15–17.PH Williams2000Are Bombus lucorum and magnus separate species?BWARS Newsletter20001517

33. Rasmont P, Scholl A, De Jonghe R, Obrecht E, Adamski A (1986) Identité et variabilité des mâles de bourdons du genre Bombus Latreille sensu stricto en Europe occidentale et centrale (Hymenoptera, Apidae, Bombinae). Revue suisse Zoologique 93: 661–682.P. RasmontA. SchollR. De JongheE. ObrechtA. Adamski1986Identité et variabilité des mâles de bourdons du genre Bombus Latreille sensu stricto en Europe occidentale et centrale (Hymenoptera, Apidae, Bombinae).Revue suisse Zoologique93661682

34. Scholl A, Obrecht E (1983) Enzymelektrophoretische Untersuchungen zur Artabgrenzung im Bombus lucorum - Komplex (Apidae, Bombini). Apidologie 14: 65–78.A. SchollE. Obrecht1983Enzymelektrophoretische Untersuchungen zur Artabgrenzung im Bombus lucorum - Komplex (Apidae, Bombini).Apidologie146578

35. Pamilo P, Varvio-Aho S-L, Pekkarinen A (1984) Genetic variation in bumblebees (Bombus, Psithyrus) and putative sibling species of Bombus lucorum. Hereditas 101: 245–251.P. PamiloS-L Varvio-AhoA. Pekkarinen1984Genetic variation in bumblebees (Bombus, Psithyrus) and putative sibling species of Bombus lucorum.Hereditas101245251

36. Pamilo P, Tengö J, Rasmont P, Pirhonen K, Pekkarinen A, et al. (1997) Pheromonal and enzyme genetic characteristics of the Bombus lucorum species complex in northern Europe. Entomologica fennica 7: 187–194.P. PamiloJ. TengöP. RasmontK. PirhonenA. Pekkarinen1997Pheromonal and enzyme genetic characteristics of the Bombus lucorum species complex in northern Europe.Entomologica fennica7187194

37. Bertsch A, Schweer H, Titze A (2004) Discrimination of the bumblebee species Bombus lucorum, B. cryptarum and B. magnus by morphological characters and male labial gland secretions. Beiträge zur Entomologie 54: 365–386.A. BertschH. SchweerA. Titze2004Discrimination of the bumblebee species Bombus lucorum, B. cryptarum and B. magnus by morphological characters and male labial gland secretions.Beiträge zur Entomologie54365386

38. Bertsch A, Schweer H, Titze A, Tanaka H (2005) Male labial gland secretions and mitochondrial DNA markers support species status of Bombus cryptarum and B. magnus (Hymenoptera, Apidae). Insectes Sociaux 52: 45–54.A. BertschH. SchweerA. TitzeH. Tanaka2005Male labial gland secretions and mitochondrial DNA markers support species status of Bombus cryptarum and B. magnus (Hymenoptera, Apidae).Insectes Sociaux524554

39. Pedersen BV (2002) European bumblebees (Hymenoptera: Bombini) - phylogenetic relationships inferred from DNA sequences. Insect Systematics & Evolution 33: 361–386.BV Pedersen2002European bumblebees (Hymenoptera: Bombini) - phylogenetic relationships inferred from DNA sequences.Insect Systematics & Evolution33361386

40. Murray TE, Fitzpatrick U, Brown MJF, Paxton RJ (2008) Cryptic species diversity in a widespread bumble bee complex revealed using mitochondrial DNA RFLPs. Conservation Genetics 9: 653–666.TE MurrayU. FitzpatrickMJF BrownRJ Paxton2008Cryptic species diversity in a widespread bumble bee complex revealed using mitochondrial DNA RFLPs.Conservation Genetics9653666

41. Waters J, Darvill B, Lye GC, Goulson D (2010) Niche differentiation of a cryptic bumblebee complex in the Western Isles of Scotland. Insect Conservation and Diversity 2010: 1–7.J. WatersB. DarvillGC LyeD. Goulson2010Niche differentiation of a cryptic bumblebee complex in the Western Isles of Scotland.Insect Conservation and Diversity201017

42. Williams PH, Brown MJF, Carolan JC, An J-D, Goulson D, et al. (2012) Unveiling cryptic species of the bumblebee subgenus Bombus s. str. world-wide with COI barcodes (Hymenoptera: Apidae). Systematics and Biodiversity 10: PH WilliamsMJF BrownJC CarolanJ-D AnD. Goulson2012Unveiling cryptic species of the bumblebee subgenus Bombus s. str. world-wide with COI barcodes (Hymenoptera: Apidae).Systematics and Biodiversity10in press. in press.

43. Monaghan MT, Wild R, Elliot M, Fujisawa T, Balke M, et al. (2009) Accelerated species inventory on Madagascar using coalescent-based models of species delineation. Systematic Biology 58: 298–311.MT MonaghanR. WildM. ElliotT. FujisawaM. Balke2009Accelerated species inventory on Madagascar using coalescent-based models of species delineation.Systematic Biology58298311

44. Ballard JWO, Whitlock MC (2004) The incomplete natural history of mitochondria. Molecular Ecology 13: 729–744.JWO BallardMC Whitlock2004The incomplete natural history of mitochondria.Molecular Ecology13729744

45. Hebert PDN, Ratnasingham S, deWaard JR (2003) Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proceedings of the Royal Society of London (B) 270: S96–S99.PDN HebertS. RatnasinghamJR deWaard2003Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species.Proceedings of the Royal Society of London (B)270S96S99

46. Baker AJ, Tavares ES, Elbourne RF (2009) Countering criticisms of single mitochondrial DNA gene barcoding in birds. Molecular Ecology Resources 9(supplement 1): 257–268.AJ BakerES TavaresRF Elbourne2009Countering criticisms of single mitochondrial DNA gene barcoding in birds.Molecular Ecology Resources9supplement 1257268

47. Meusnier I, Singer GAC, Landry J-F, Hickey DA, Hebert PDN, et al. (2008) A universal DNA mini-barcode for biodiversity analysis. Bio Med Central Genomics 9: 214[214 pp.].I. MeusnierGAC SingerJ-F LandryDA HickeyPDN Hebert2008A universal DNA mini-barcode for biodiversity analysis.Bio Med Central Genomics9214[214 pp.]

48. Brower AVZ (2006) Problems with DNA barcodes for species delimitation: ‘ten species’ of Astraptes fulgerator reassessed (Lepidoptera: Hesperiidae). Systematics and Biodiversity 4: 127–132.AVZ Brower2006Problems with DNA barcodes for species delimitation: ‘ten species’ of Astraptes fulgerator reassessed (Lepidoptera: Hesperiidae).Systematics and Biodiversity4127132

49. Packer L, Sheffield CS, Gibbs J, de Silva N, Best LR, et al. The campaign to barcode the bees of the world: progress, problems, prognosis. In: Yurrita CL, editor. 2009; Centro de Estudios Conservacionistas. Universidad de San Carlos de Guatemala. pp. 178–180.L. PackerCS SheffieldJ. GibbsN. de SilvaLR BestThe campaign to barcode the bees of the world: progress, problems, prognosis.CL Yurrita2009; Centro de Estudios ConservacionistasUniversidad de San Carlos de Guatemala178180

50. Hebert PDN, Penton EH, Burns JM, Janzen DH, Hallwachs W (2004) Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proceedings of the National Academy of Sciences 101: 14812–14817.PDN HebertEH PentonJM BurnsDH JanzenW. Hallwachs2004Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator.Proceedings of the National Academy of Sciences1011481214817

51. Huelsenbeck JP, Ronquist F (2001) MRBAYES: Bayesian inference of phylogeny. Bioinformatics 17: 754–755.JP HuelsenbeckF. Ronquist2001MRBAYES: Bayesian inference of phylogeny.Bioinformatics17754755

52. Ronquist F, Huelsenbeck JP (2003) MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.F. RonquistJP Huelsenbeck2003MRBAYES 3: Bayesian phylogenetic inference under mixed models.Bioinformatics1915721574

53. Posada D (2008) jModelTest: phylogenetic model averaging. Molecular Biology and Evolution 25: 1253–1256.D. Posada2008jModelTest: phylogenetic model averaging.Molecular Biology and Evolution2512531256

54. Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BioMed Central Evolutionary Biology 7: 214.AJ DrummondA. Rambaut2007BEAST: Bayesian evolutionary analysis by sampling trees.BioMed Central Evolutionary Biology7214

55. Nylander JAA, Wilgenbusch JC, Warren DL, Swofford DL (2008) AWTY (are we there yet?): a system for graphical exploration of MCMC convergence in Bayesian phylogenetics. Bioinformatics 24: 581–583.JAA NylanderJC WilgenbuschDL WarrenDL Swofford2008AWTY (are we there yet?): a system for graphical exploration of MCMC convergence in Bayesian phylogenetics.Bioinformatics24581583

56. Cameron SA, Hines HM, Williams PH (2007) A comprehensive phylogeny of the bumble bees (Bombus). Biological Journal of the Linnean Society 91: 161–188.SA CameronHM HinesPH Williams2007A comprehensive phylogeny of the bumble bees (Bombus).Biological Journal of the Linnean Society91161188

57. Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Molecular Biology and Evolution 24: 1596–1599.K. TamuraJ. DudleyM. NeiS. Kumar2007MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0.Molecular Biology and Evolution2415961599

58. Kimura M (1980) A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. Journal of Molecular Evolution 16: 111–120.M. Kimura1980A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences.Journal of Molecular Evolution16111120

59. ICZN 1999 International code of zoological nomenclature London International Commission on Zoological Nomenclature xxix+306

60. Bertsch A (2010) A phylogenetic framework for the bumblebee species of the subgenus Bombus sensu stricto based on mitochondrial DNA markers, with a short description of the neglected taxon B. minshanicola Bischoff, 1936 n. status. Beiträge zur Entomologie 60: 471–487.A. Bertsch2010A phylogenetic framework for the bumblebee species of the subgenus Bombus sensu stricto based on mitochondrial DNA markers, with a short description of the neglected taxon B. minshanicola Bischoff, 1936 n. status.Beiträge zur Entomologie60471487

61. Tkalcu B (1974) Eine Hummel-Ausbeute aus dem Nepal-Himalaya (Insecta, Hymenoptera, Apoidea, Bombinae). Senckenbergiana biologica 55: 311–349.B. Tkalcu1974Eine Hummel-Ausbeute aus dem Nepal-Himalaya (Insecta, Hymenoptera, Apoidea, Bombinae).Senckenbergiana biologica55311349

62. Rubinoff D, Cameron SA, Will K (2006) A genomic perspective on the shortcomings of mitochondrial DNA for “barcoding” identification. Journal of Heredity 97: 581–594.D. RubinoffSA CameronK. Will2006A genomic perspective on the shortcomings of mitochondrial DNA for “barcoding” identification.Journal of Heredity97581594

63. Williams PH, An J-D, Huang J-X (2011) The bumblebees of the subgenus Subterraneobombus: integrating evidence from morphology and DNA barcodes (Hymenoptera, Apidae, Bombus). Zoological Journal of the Linnean Society 163: 813–862.PH WilliamsJ-D AnJ-X Huang2011The bumblebees of the subgenus Subterraneobombus: integrating evidence from morphology and DNA barcodes (Hymenoptera, Apidae, Bombus).Zoological Journal of the Linnean Society163813862

64. Skorikov AS (1913) Neue Hummelformen (Hymenoptera, Bombidae). V Russkoe éntomologicheskoe Obozrênie 13: 171–175.AS Skorikov1913Neue Hummelformen (Hymenoptera, Bombidae).V Russkoe éntomologicheskoe Obozrênie13171175

65. Goulson D (2010) Impacts of non-native bumblebees in Western Europe and North America. Applied Entomology and Zoology 45: 7–12.D. Goulson2010Impacts of non-native bumblebees in Western Europe and North America.Applied Entomology and Zoology45712