Open access journal: http://periodicos.uefs.br/ojs/index.php/sociobiology ISSN: 0361-6525 DOI: 10.13102/sociobiology.v68i3.5910Sociobiology 68(3): e5910 (September, 2021) Bombus terrestris is the most abundant and widespread bumblebee in the West Palearctic where it has differentiated into eight subspecies that differ in morphology and genetic characters, behaviour and phenological traits (Rasmont et al., 2008). This species has a great economic and ecological importance because several B. terrestris subspecies have been massively reared and exported worldwide for crop pollination (Velthuis & van Doorn, 2004). This commercialization has resulted in reproductive interference with native species (Tsuchida et al., 2019) or hybridization with native subspecies (Cejas et al., 2018, 2020; Seabra et al., 2019; Bartomeus et al., 2020), multiple invasions with competitive displacement of native taxa (Matsumura et al., 2004; Ings et al., 2006; Inoue et al., 2008; Nagamitsu et al., 2010; Morales et al., 2013) and possibly pathogen spillover to local bumblebees (Goka et al., 2001; Arbetman et al., 2013; Schmid-Hempel et al., 2014). Abstract The taxonomic status of Bombus terrestris subspecies is complex and has deep implications in the management of commercial bumblebees for crop pollination as well as in the establishment of appropriate conservation plans. Herein, the complete mitogenome of the endemic Canary Islands subspecies Bombus terrestris canariensis is newly sequenced and compared with available mitochondrial sequences in order to shed light into its taxonomic status. The obtained sequence of the mitochondrial genome was 17,300 bp in length and contained 37 genes, including 13 protein-coding genes (PCGs), two rRNAs, and 22 tRNAs and a partial sequence of the AT rich control region. The phylogenetic analysis of PCGs of the mitogenome was congruent with its subspecies status and a close relationship with the North African subspecies africanus as previously suggested. The sequencing of the mitogenome of B. t. canariensis provides useful genetic information to study the conservation genetics and genetic diversity of these island bumblebee populations. Sociobiology An international journal on social insects Carlos Ruiz1, Diego Cejas2, Irene Muñoz3, Pilar De la Rua4 Article History Edited by Marco Antonio Costa, UESC, Brazil Received 08 October 2020 Initial acceptance 07 February 2021 Final acceptance 16 May 2021 Publication date 13 August 2021 Keywords Pollinators, Mitochondrial DNA, Conservation, Phylogeny, Bombus terrestris canariensis, Insularity. Corresponding author Pilar De la Rua University of Murcia Calle Campus Universitario, 11 30100 Murcia, Espanha. E-Mail: pdelarua@um.es Given these problems, the trade of the different commercial B. terrestris subspecies has been restricted in few places such as the Israel, Norway, Turkey, the United Kingdom and the Canary Islands, where only the local subspecies B. terrestris canariensis is allowed to be commercialized (Velthuis & van Doorn, 2006; Lecocq et al., 2016). Molecular, morphological and pheromone markers have been used to clarify the taxonomy of the terrestris subspecies complex. In particular, the endemic taxon of the Canary Island was firstly described as a B. terrestris subspecies (Perez, 1895) and later elevated to species status by Erlandsson (1979) based on its distinct coloration pattern and geographic isolation. Several mtDNA and nuclear markers have supported both species (Widmer, 1998) and subspecies status (Estoup et al., 1996) or remain unclear (Bertsch, 2010). Integrative taxonomic approaches combining independent 1 - Dpto. Biología Animal, Edafología y Geología, Facultad de Ciencias, Universidad de La Laguna. La Laguna, 38206 Tenerife, Espanha 2 - Université de Mons, Laboratoire de Zoologie. Place du Parc 23, 7000 Mons, Belgique 3 - Dpto. de Zoología y Antropología Física, Facultad de Veterinaria, Universidad de Murcia. 30100 Murcia, Espanha 4 - University of Murcia. 30100 Murcia, Espanha SHORT NOTE Characterizing the Mitogenome of the Endemic Bumblebee Subspecies from the Canary Islands for Conservation Purposes Carlos Ruiz, Diego Cejas, Irene Muñoz, Pilar De la Rua – Mitogenome sequence of endemic bumblebee subspecies2 markers such as pheromones, morphological and molecular data generated contrasting results supporting both species (De Meulemeester, 2012) and subspecies status (Lecocq et al., 2016). Furthermore, crossing experiments showed that canariensis interbreed with other B. terrestris subspecies leading to fertile offspring (van den Eijnde & de Ruijter, 2000), what has been used as an evidence of its subspecies status (Velthuis & van Doorn, 2006), although interspecific mating cannot be rule out (Kondo et al., 2009; Yoon et al., 2009). The taxonomy of B. terrestris subspecies has deep implications also in the establishment of appropriate conservation plans (Lecocq et al., 2016). Therefore, a solid taxonomic status must be achieved using novel molecular approaches. Mitogenome sequence combined with other molecular markers has successfully resolved phylogenetic and taxonomic issues in several taxa, resulting in an effective tool for bee phylogeny and conservation genetics (Nishimoto et al., 2018; Du et al., 2016; Lin et al., 2019). Here we have sequenced the complete mitogenome of the canariensis taxon and compared it with available mitochondrial sequences to shed light into its taxonomic status. To ensure an adequate amount of mitochondrial DNA for sequencing (Cejas et al., 2020), muscle tissue was dissected from each of ten B. t. canariensis individuals collected from flowers at a single locality on La Gomera (Canary Islands, Spain) and pooled. DNA extraction and mtDNA enrichment was done using the minipred kit of Qiagen (Hilden, Germany). Sequencing was performed on a HiSeq2000 (Illumina) using an Illumina TruSeq Nano DNA Library Prep Kit for 350 base pairs (bp) (Macrogen, South Korea). An agilent Technologies 2100 Bioanalyzer with a DNA 1000 chip was used to measure the size of the raw library reads. Reads were filtered with the software FastQC (Babraham Bioinformatics 2012) before assembling them into contigs and scaffolds as in Cejas et al. (2020). Annotation of PCGs, transfer RNA (tRNA) genes and ribosomal RNA (rRNA) genes for the consensus sequence was obtained by comparison with previous annotated genomes in the web server MITOS (Bernt et al., 2013). In order to resolve the phylogenetic position of canariensis, ten mitogenomes of related species and subspecies were obtained from GenBank (Supplementary Material 1). Two mitochondrial markers (cox1 and cytb) comprising 1,307 bp were used for all the subspecies with no mitogenome available (except for africanus with only sequence data of cox1). The phylogenetic position of canariensis was estimated from a concatenated dataset including the 13 protein coding genes (PCGs) using the software MrBayes 3.2.6 (Ronquist & Huelsenbeck, 2003). The best-fitting substitution model was assessed with IQ-TREE. A phylogenetic tree was constructed using B. consobrinus (MF995069) as outgroup. Fig 1. Physical map of the obtained sequence mitochondrial genome of Bombus terrestris canariensis in absence of the complete A+T control region. Protein coding genes (PCG) are indicated in light blue, tRNA genes in green, and rRNA genes in dark blue. The GC content is indicated in black. Sociobiology 68(3): e5910 (September, 2021) 3 The A+T control region was not completely sequenced, owing to its extreme variability and the pooling approach followed. Thus, the mitogenome of B. t. canariensis, was sequenced to a length of 17,300 bp (GenBank accession number MW959771) (Fig 1). The sequenced part varies in length in comparison to other B. terrestris subspecies (B. t. terrestris: 17,232 bp; B. t. lusitanicus 17,049 pb) mainly due to the presence of indels in intergenic regions. Gene order was consistent with published data (Cejas et al., 2020). It contained 13 protein-coding genes (PCGs), two rRNAs, and 22 tRNAs and a partial sequence of the AT rich control region. The average A+T content was 86.4%, slightly higher than that in other B. terrestris subspecies (B. t. lusitanicus and B. t. terrestris 86%). This A+T content was higher in non PCGs (87.5%) than in the PCGs region (83.7%). Most of the variation observed in the available mitogenomes of terrestris subspecies occurred in non PCGs. PCGs of B. t. canariensis showed 193 SNPs, especially in NAD4 (34), NAD5 (33) and cox1 (27) genes. Transitions (83%) where more frequent than transversions in the SNPs. No indels were found within the PCGs whereas a 6 bp deletion was observed in the large RNA gene sequence. Fig 2. Bayesian phylogenetic tree showing the relationship between Bombus terrestris canariensis (in red) and other eight B. terrestris subspecies. An arrow indicates the B. t. canariensis mitogenome sequenced in this study and bold names the mitogenomes used. Empty and filled circles indicates posterior probabilities of 0.90-0.95 and >0.95 respectively. Carlos Ruiz, Diego Cejas, Irene Muñoz, Pilar De la Rua – Mitogenome sequence of endemic bumblebee subspecies4 The combined analysis of concatenated PCGs of mitogenomes and available markers of the different subspecies showed high support for all the basal nodes and low support for the nodes within terrestris clade (Fig 2). B. terrestris appeared within the subgenus Bombus, related with species such as B. lucorum, B. hypocrita and B. cryptarum. All the studied terrestris subspecies including canariensis appeared as monophyletic in a supported clade, thus reinforcing the subspecies status of canariensis. However, other taxon such as xanthopus that has been previously elevated to species (Lecocq et al., 2015) appeared within the B. terrestris clade. These results should be taken with caution as only three mitogenomes were available for the terrestris subspecies and only partial mitochondrial data (1,307 bp) have been analysed for the remaining subspecies. Given its importance as a commercial species, further mitogenomes of the remaining subspecies will establish a solid taxonomy of the group. The sequencing of the B. t. canariensis genome itself provides additional genetic information useful for studying the conservation genetics of these island bumblebee populations, creating a framework for establishing conservation programs for pollination networks of B. terrestris subspecies with locally endemic flora. Acknowledgments This work was supported by the projects E-RTA2014- 00003-C03 (INIA, Spain; European Regional Development Fund) and 19908/GERM/2015 of Regional Excellence (Fundación Séneca, CARM). IM is supported by a MINECO Spanish postdoctoral grant “Juan de la Cierva- Incorporación” (JCI2018-036614-I). Authors’ Contribution All the authors have given their consent to participate in the redaction of the manuscript. CR, IM and PDlR conceived the ideas and designed the experiments. DC did the DNA extractions and performed the assembly of the mitogenome. CR performed most of the analyses with the assistance of DC. CR and PDlR wrote the manuscript with input from IM. All the authors critically reviewed the manuscript. References Arbetman, M.P., Meeus, I., Morales, C.L., Aizen, M.A. & Smagghe, G. (2013). Alien parasite hitchhikes to Patagonia on invasive bumblebee. Biological Invasions, 15: 489 - 494. doi: 10.1007/s10530-012-0311-0 Bernt, M., Donath, A., Jühling, F., Externbrink, F., Florentz, C., Fritzsch, G., ... & Stadler, P. F. (2013). MITOS: improved de novo metazoan mitochondrial genome annotation. Molecular Phylogenetics and Evolution, 69: 313-319. doi: 10.1016/j. ympev.2012.08.023 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. (Hymenoptera: Apidae: Bombus). Beiträge Zur Entomologie, 60: 471-487. doi: 10.21248/contrib.entomol.60.2.471-487 Bartomeus, I., Molina, F.P., Hidalgo-Galiana, A. & Ortego, J. (2020). Safeguarding the genetic integrity of native pollinators requires stronger regulations on commercial lines. Ecological Solutions and Evidence, 1(1), e12012 Cejas, D., Ornosa, C., Muñoz, I. & De la Rúa, P. (2018). Searching for molecular markers to differentiate Bombus terrestris (Linnaeus) subspecies in the Iberian Peninsula. Sociobiology, 65: 558-565. doi: 10.13102/sociobiology.v65i4. 3442 Cejas, D., López-López, A., Muñoz, I., Ornosa, C. & De la Rúa, P. (2020). Unveiling introgression in bumblebee (Bombus terrestris) populations through mitogenome-based markers. Animal Genetics, 51: 70 -77. doi: 10.1111/age.12874 De Meulemeester, T. (2012). Approche intégrative dans la systématique de taxons complexes: bourdons et abeilles fossiles. Disertation, Université de Mons Du, Q., Bi, G,. Zhao, E., Yang, J., Zhang, Z. & Liu, G. (2016). Complete mitochondrial genome of Bombus terrestris (Hymenoptera: Apidae). Mitochondrial DNA Part A, 27: 4455-4456. doi: 10.3109/19401736.2015.1089568 Erlandsson, S. (1979). Bombus canariensis Pérez, 1895 n. stat and B. maderensis n. sp. from the Macaronesian Islands. Insect Systematics and Evolution, 10: 187-192 Estoup, A., Solignac, M., Cornuet, J.M., Goudet, J. & Scholl, A. (1996). Genetic differentiation of continental and island populations of Bombus terrestris (Hymenoptera: Apidae) in Europe. Molecular Ecology, 5: 19-31 Goka, K., Okabe, K., Yoneda, M. & Niwa, S. (2001). Bumblebee commercialization will cause worldwide migration of parasitic mites. Molecular Ecology, 10: 2095-2099. doi: 10.10 46/j.0962-1083.2001.01323.x Ings, T.C., Ward, N.L.L. & Chittka, L. (2006). Can commercially imported bumble bees out-compete their native conspecifics? Journal of Applied Ecology, 43: 940-948. doi: 10.1111/j.1365- 2664.2006.01199.x Inoue, M.N., 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. doi: 10.1007/ s10841-007-9071-z Kondo, N.I., Yamanaka, D., Kanbe, Y., Kunitake, Y.K., Yoneda, M., Tsuchida, K. & Goka, K. (2009). Reproductive disturbance of Japanese bumblebees by the introduced https://doi.org/10.1007/s10530-012-0311-0. https://doi.org/10.1016/j.ympev.2012.08.023 https://doi.org/10.1016/j.ympev.2012.08.023 Sociobiology 68(3): e5910 (September, 2021) 5 European bumblebee Bombus terrestris. Naturwissenschaften, 96: 467-475. doi: 10.1007/s00114-008-0495-4 Lecocq, T., Brasero, N., De Meulemeester, T., Michez, D., Dellicour, S., Lhomme, P., et al (2015). An integrative taxonomic approach to assess the status of Corsican bumblebees: implications for conservation. Animal Conservation, 18: 236-248. doi: 10.1111/acv.12164 Lecocq, T., Coppée, A., Michez, D., Brasero, N., Rasplus, J.Y., Valterová, I. & Rasmont, P. (2016). The alien’s identity: consequences of taxonomic status for the international bumblebee trade regulations. Biological Conservation, 195: 169-176. doi: 10.1016/j.biocon.2016.01.004 Lin, G., Jiang, K., Su, T., He, B., Huang, Z. & Zhao, F. (2019). Complete mitochondrial genome of Bombus waltoni (Hymenoptera: Apidae). Mitochondrial DNA Part B, 4: 1484- 1485. doi: 10.1080/23802359.2019.1601528 Matsumura, C., Yokoyama, J. & Washitani I. (2004). Invasion status and potential ecological impacts of an invasive alien bumblebee, Bombus terrestris L. (Hymenoptera: Apidae) naturalized in Southern Hokkaido, Japan Global Environmental Research, 8: 51-66 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: 529-534. doi: 10.1890/120321 Nagamitsu, T., Yamagishi, H., Kenta, T., Inari, N. & Kato, E. (2010). Competitive effects of the exotic Bombus terrestris on native bumble bees revealed by a field removal experiment. Population Ecology, 52: 123-136. doi: 10.1007/s10144-009- 0151-7 Nishimoto, M., Umezawa, M., Okuyama, H., Kumano, N., Nomura, T. & Takahashi, J.I. (2018). The complete mitochondrial genome of the bumblebee, Bombus hypocrita sapporensis (Insecta: Hymenoptera: Apidae) from Hokkaido Island, Japan. Mitochondrial DNA Part B, 3: 354-356. doi: 10.1080/23802359.2018.1450673 Rasmont, P., Coppée, A., Michez, D. & De Meulemeester, T. (2008). An overview of the Bombus terrestris (L. 1758) subspecies (Hymenoptera: Apidae). Annales de la Société Entomologique de France, 44: 243-250. doi: 10.1080/0037 9271.2008.10697559 Ronquist, F. & Huelsenbeck, J.P. (2003). MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics, 19: 1572-1574. doi: 10.1093/bioinformatics/btg180 Schmid-Hempel, R., Eckhardt, M., Goulson, D., Heinzmann, D., Lange, C., Plischuk, S., et al (2014) The invasion of southern South America by imported bumblebees and associated parasites. Journal of Animal Ecology, 83: 823-837. doi: 10.1111/1365-2656.12185 Seabra, S.G., Silva, S.E., Nunes, V.L., Sousa, V.C., Martins, J., Marabuto, E. et al (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. doi: 10.1111/eva.12732 Tsuchida, K., Yamaguchi, A., Kanbe, Y., & Goka, K. (2019). Reproductive interference in an introduced bumblebee: Polyandry may mitigate negative reproductive impact. Insects, 10: 59. doi: 10.3390/insects10020059 van den Eijnde, J. & de Ruijter, A. (2000). Bumblebees from the Canary Islands: mating experiments with Bombus terrestris L. from the Netherlands. Proceedings of the Section Experimental and Applied Entomology, 11: 159-161 Velthuis, H.H. & van Doorn, A. (2004). The breeding, commercialization and economic value of bumblebees. In: Freitas BM, Pereira JOP (Eds) Solitary bees: Conservation, Rearing and Management in Pollination. Universitaria, Fortaleza, pp 135-149 Velthuis, H.H. & Van Doorn, A. (2006). A century of advances in bumblebee domestication and the economic and environmental aspects of its commercialization for pollination. Apidologie, 37: 421-451. doi: 10.1051/apido:2006019 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: 563-572. Yoon, H.J., Kim, S.Y., Lee, K.Y., Lee, S.B., Park, I.G. & Kim, I.S. (2009). Interspecific hybridization of the bumblebees Bombus ignitus and B. terrestris. International Journal of Industrial Entomology, 18: 41-48 Electronic Supplementary Material Supplementary file 1: Sampling information of the individuals used for the phylogenetic analysis.