References
Belkhir, K., Borsa, P., Chikhi, L., Raufaste, N., & Bonhomme, F. (2004). GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université de Montpellier II, Montpellier, France. http://kimura.univ-montp2.fr/genetix/.
Buttermore, R.E. (1997). Observations of successful Bombus terrestris (L.) (Hymenoptera: Apidae) colonies in Southern Tasmania.Australian Journal of Entomology, 36 (3), 251 – 254. https://doi.org/10.1111/j.1440-6055.1997.tb01463.x
Cejas, D., Ornosa, C., Munoz, I., & De la Rua, P. (2018). Searching for molecular markers to differentiate Bombus terrestris (Linnaeus) subspecies in the Iberian Peninsula. Sociobiology , 65 (4), 558–565. https://doi.org/ 10.13102/sociobiology.v65i4.3442
Dafni, A., Kevan, P., Gross, C.L., & Goka, K. (2010). Bombus terrestris , pollinator, invasive and pest: An assessment of problems associated with its widespread introductions for commercial purposes.Applied Entomology and Zoology, 45 (1), 101–113. https://doi.org/10.1303/aez.2010.101
Earl, D.A., & vonHoldt, B.M. (2012). STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method.Conservation Genetics Resources,  4, 359–361.
https://doi.org/10.1007/s12686-011-9548-7
Estoup, A., Scholl, A., Pouvreau, A.,& Solignac, M. (1995). Monoandry and polyandry in
bumble bees (Hymenoptera; Bombinae) as evidenced by highly variable microsatellites. Molecular Ecology, 4, 89-93. https://doi.org/10.1111/j.1365-294X.1995.tb00195.x
Estoup, A., Solignac, M., Cournet, J.M., Goudet, J., & Scholl, A.(1996). Genetic differentiation of continental and island populations of Bombus terrestris (Hymenoptera: Apidae) in Europe.Molecular Ecology, 5 (1), 19-31.
https://doi.org/10.1111/j.1365-294X.1996.tb00288.x
Excoffier, L., Laval, G., & Schneider, S. (2005). Arlequin ver. 3.0: An integrated software package for population genetics data analysis.Evolutionary Bioinformatics Online 1, 47-50.
Evanno, G., Regnaut, S., & Goudet, J. (2005).Detecting the number of clusters of individuals using the software structure: a simulation study. Molecular Ecology, 14 , 2611 –2620,
https://doi.org/10.1111/j.1365-294X.2005.02553.x
Evans, E. (2017). From humble bee to greenhouse pollination workhorse: Can we mitigate risks for bumble bees? Bee World, 94 (2), 34-41.
https://doi.org/ 10.1080/0005772X.2017.1290892
Folmer,O., Black, M., Hoeh, W., Lutz, R.,& Vrijenhoek, R. (1994). DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Molecular Marine Biology and Biotechnology, 3 (5), 294-299.
Glaubitz, J.C. (2004). convert: A user‐friendly program to reformat diploid genotypic data for commonly used population genetic software packages, Molecular Ecology Notes, 4, 309-310.https://doi.org/10.1111/j.1471-8286.2004.00597.x
Goka, K. (2010). Introduction to the special feature for ecological risk assessment of
introduced bumblebees: Status of the European bumblebee, Bombus terrestris, in Japan as a beneficial pollinator and an invasive alien species. Journal of Applied Entomology and Zoology, 45(1), 1–6. https://doi.org/10.1303/aez.2010.1
Gosterit, A., & Gurel, F. (2005) Comparison of development patterns of imported and native Bombus terrestris L. (Hymenoptera: Apidae) colonies in the Mediterranean coastal region. Turkish Journal of Veterinary and Animal Sciences , 29 , 393–398.
Gosterit, A.,& Baskar, V C. (2016). Impacts of commercialization on the developmental characteristics of native Bombus terrestris (L.) colonies. Insectes Sociaux , 63, 609–614.
https://doi.org/10.1007/s00040-016-0507-x
Goudet, J. (1995). FSTAT (Version 1.2): A computer program to calculate F-Statistics. Journal of the Heredity , 86 (6), 485-486.
https://doi.org/10.1093/oxfordjournals.jhered.a111627
Goulson, D. (2010). Impacts of non‐native bumblebees in Western Europe and North America. Applied Entomology and Zoology , 45 (1), 7–12. https://doi.org/10.1303/aez.2010.7
Gurel, F., Gosterit, A., & Eren, Ö. (2008). Life-cycle and foraging patterns of native Bombus terrestris (L.) (Hymenoptera, Apidae) in the Mediterranean region. Insectes Sociaux ,55 , 123–128. https://doi.org/10.1007/s00040-008-0984-7
Huson, D.H., & Bryant, D. (2006). Application of phylogenetic networks in evolutionary studies. Molecular Biology Evolution,23 (2), 254-267. https://doi.org/10.1093/molbev/msj030
Ings, T.C.,Ward, N.L.,& Chittka, L. (2006) Can commercially imported bumble bees out-compete their native conspecifics? Journal of Applied Ecology, 43 , 940–948.
https://doi.org/10.1111/j.1365-2664.2006.01199.x
Jermiin, L.S., & Crozier, R.H. (1994).The cytochrome b region in the mitochondrial DNA of the ant Tetraponera rufoniger: Sequence divergence in hymenoptera may be associated with nucleotide content.Journal Molecular Evolution 38 , 282–294. https://doi.org/10.1007/BF00176090
Koide, T., Yamada, Y., Yabe, K.,& Yamashita, F. (2008). Methods of netting greenhouses to prevent the escape of bumblebees. Japanese Journal of Applied Entomology and Zoology, 52 , 19–26. https://doi.org/10.1303/jjaez.2008.19
Kraus, F.B., Szentgyorgyi, H., Rozej, E., Rhode, M., Moron, D., Woyciechowski, M.,& Moritz, R.F.A. (2011).Greenhouse bumblebees (Bombus terrestris ) spread their genes into the wild.Conservation Genetics, 12 (1), 187-192. https://doi.org/10.1007/s10592-010-0131-7
Kumar, S., Stecher, G., Li, M., Knyaz. C., & Tamura, K. (2018). MEGA X: Molecular
evolutionary genetics analysis across computing platforms. Molecular Biology and Evolution, 35(6), 1547-1549. https://doi.org 10.1093/molbev/msy096
Larkin, M.A., Blackshields, G., Brown, N.P,, Chenna, R., McGettigan, P.A., McWilliam, H,…. Higgins, D.G. (2007). Clustal W and Clustal version 2.0. Bioinformatics, 23 (21), 2947-2948. https://doi.org/10.1093/bioinformatics/btm404
Lecocq, T., Vereecken, N.J., Michez, D., Dellicour, S., Lhomme, P., Valterova, I.,.. Rasmont, P.(2013). Patterns of genetic and reproductive traits differentiation in mainland vs. Corsican populations of bumblebees. Plos One, 8 (6), e65642
https://doi.org/10.1371/journal.pone.0065642
Lecocq, T., Rasmont, P., Harpke, A., & Schweiger, O. (2016). Improving international trade regulation by considering intraspecific variation for invasion risk assessment of commercially traded species: TheBombus terrestris case. Conservation Letters, 9 (4), 281–289. https://doi.org 10.1111/conl.12215
Lecocq, T., Coppee, A., Michez, D., Brasero N., Rasplus, J.Y., Valterova, I., & Rasmont, P. (2016). The alien’s identity: consequences of taxonomic status for the international bumble bee trade regulations.Biological Conservation, 195 , 169-176.
https ://doi. org/10.1016/j.biocon.2016.01.004
Librado, P.J.R., & Rozas, J. (2009). DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25, 1451-1452.
https ://doi. org/10.1093/bioinformatics/btp187
Michener, C.D. (2000) . The Bees of the World . The Johns Hopkins University Press, Baltimore, Maryland.
Montalva, J., Dudley, L., Arroyo, M.K., Retamales, H., & Abrahamovich, A.H. (2011). Geographic distribution and associated flora of native and introduced bumble bees (Bombus spp.) in Chile. Journal of Apicultural Research, 50 (1), 11-21. https://doi.org/10.3896/IBRA.1.50.1.02
Morales, C.L., Arbetman, M.P., Cameron, S.A.,& Aizen, M.A. (2013). Rapid ecological replacement of a native bumble bee by invasive species.Frontiers in Ecology and the Environment, 11 (10), 529–534. https://doi.org/10.1890/120321
Moreira, A.S., Horgan, F.G., Murray, T.,& Kakouli-Duarte, T.(2015). Population genetic
structure of Bombus terrestris in Europe: Isolation and genetic differentiation of Irish and British populations. Molecular Ecology, 24(13), 3257–3268. https://doi.org/10.1111/mec.13235
Ozbek, H. 1997. Bumblebee fauna of Turkey with distribution maps (Hymenoptera: Apidae,
Bombinae) Part 1: Alpigenobombus Skorikov, Bombias Robertson andBombus latreille. Turkish Journal of Entomology, 21(1), 37-56.
Park, S.D.E. (2001). Microsatellite toolkit for Microsoft Excel, Version 3.1. Animal
Genomics Laboratory, University College Dublin, Ireland. http://animalgenomics.ucd.ie/sdepark/ms-toolkit
Pritchard, J.K., Matthew S., & Donnelly, P. (2000). Inference of population structure using
multilocus genotype data. Genetics, 155(2), 945-959.
Rasmont, P., Coppee, A., Michez, D., & Meullemeester, T. (2008). An overview of the Bombus terrestris (L. 1758) subspecies (Hymenoptera: Apidae.) Annales de la Societe Entomologique de France, 44 (2), 243-25. https://doi.org/10.1080/00379271.2008.10697559
Reber-Funk, C., Schmid-Hempel, R.,& Schmid-Hempel, P. (2006). Microsatellite loci for Bombus spp. Molecular Ecology Notes , 6 (1), 83-86.
https://doi.org/10.1111/j.1471-8286.2005.01147.x
Schneider, .S., Roessli, D.,& Excoffier, L. 2000. ARLEQUIN ver. 2.00: A software for population genetics data analysis. Genetics and Biometry Laboratory, University of Geneva, Switzerland. https://anthro.unige.ch/arlequin
Seabra , S.G., Silva, S.E., Nunes V.L., Sousa, V.C., Martins, J., Marabuto,.E….Paulo, O.S. (2019). Genomic signatures of introgression between commercial and native bumblebees, Bombus terrestris , in western Iberian Peninsula—Implications for conservation and trade regulation. Evolutionary Applications, 12,679-691. https://doi.org10.1111/eva.12732
Stolle, E., Rohde, M., Vautrin, D., Solignac, M., Schmid-Hempel, P., Schmid-Hempel, R.,&
Moritz, R.F.A. (2009). Novel microsatellite DNA loci for Bombus terrestris (Linnaeus, 1758). Molecular Ecology Resources,9(5), 1345-1352. https://doi.org/10.1111/j.1755-0998.2009.02610.x
Tamura, K., & Nei, M. (1993). Estimation of the number of nucleotide substitutions in the
control region of mitochondrial DNA in humans and chimpanzees. Molecular Biology and Evolution, 10, 512-526. https://doi.org/10.1093/oxfordjournals.molbev.a040023
Tamura, K. Stecher, G., Peterson, D., Filipski, A.,& Kumar, S. (2013). MEGA6: Molecular evolutionary genetics analysis, Version 6.0.Molecular Biology and Evolution , 30 (12), 2725-2729. https://doi.org/10.1093/molbev/mst197
Velthuis,H.H.W., & van Doorn, A. (2006). A century of advances in bumblebee domestication and the economic and environmental aspects of its commercialization for pollination. Apidologie, 37 (4), 421-451. http://dx.doi.org/10.1051/apido:2006019
Walsh, P.S., Metzger, D.A.,& Higuchi, R. (1991). Chelex 100 as a medium for simple
extraction of DNA for PCR-based typing from forensic material.BioTechniques, 10(4), 506-513.
Widmer, A., Schmid-Hempel, P., Estoup, A., & Scholl, A. (1998). Population genetic structure and colonization history of Bombus terrestris s.l. (Hymenoptera: Apidae) from the Canary Islands and Madeira. Heredity, 81 (5) 563-572.
https://doi.org/10.1046/j.1365-2540.1998.00407
Williams, P.H. (1998). An annotated checklist of bumble bees with an analysis of patterns of description (Hymenoptera: Apidae, Bombini).Bulletin of the Natural History Museum Entomology, 67 (1), 79-152.
Yeh, F.C., Yang, R.C., Boyle, T.B.J., Ye, Z.H.,& Mao, J.X. (1997).POPGENE, the user-friendly shareware for population genetic analysis . Molecular Biology and Biotechnology Centre, University of Alberta, Canada. http://www.ualberta.ca/ fyeh/.
Yeninar, H., Duchateau, M.J., Kaftanoglu, O., & Velthuis, H. (2000). Colony developmental patterns in different local populations of the Turkish bumble bee, Bombus terrestris dalmatinus. Journal of Apicultural Research, 39 (3-4), 107-116. https://doi.org/10.1080/00218839.2000.11101030
Zhang, Z., Schwartz, S., Wagner, L., & Miller, W. (2000). A greedy algorithm for aligning
DNA sequences. Journal of Computational Biology, 7(1-2), 203-214. https://doi.org/10.1089/10665270050081478
Table 1 The names, codes and GPS positions of the sampling locations