Genomic signatures of introgression between commercial and native bumblebees, Bombus terrestris, in western Iberian Peninsula-Implications for conservation and trade regulation

Human-mediated introductions of species may have profound impacts on native ecosystems. One potential impact with largely unforeseen consequences is the potential admixture of introduced with autochthonous species through hybridization. Throughout the world, bumblebees have been deliberately introduced for crop pollination with known negative impacts on native pollinators. Given the likely allochthonous origin of commercial bumblebees used in Portugal (subspecies Bombus terrestris terrestris and B. t. dalmatinus), our aim was to assess their putative introgression with the native Iberian subspecies B. terrestris lusitanicus. We analysed one mitochondrial gene, cytochrome c oxidase subunit I (COX1) and genomic data involving thousands of genome-wide restriction-site-associated DNA markers (RAD-seq). In the mitochondrial COX1 analyses, we detected one relatively common haplotype in commercial bumblebees, also present in wild samples collected nearby the greenhouses where the commercial hives are used. In the RAD-seq analysis, we found a clear genetic differentiation between native and commercial lineages. Furthermore, we detected candidate hybrids in the wild, as well as putatively escaped commercial bumblebees, some of which being potentially fertile males. Although we cannot assess directly the fitness effects of introgressed alleles, there is a risk of maladaptive allele introgression to the local bumblebee subspecies, which can negatively impact autochthon populations. One immediate recommendation to farmers is for the proper disposal of hive boxes, after their use in greenhouses, so as to minimize the risk of escapees contaminating native populations. On the other hand, the feasibility of using local subspecies B. t. lusitanicus, preferably with local production, should be evaluated.