Population genetic structure of Bombus terrestris in Europe: Isolation and genetic differentiation of Irish and British populations

Does coevolution in refugia drive mimicry in bumble bees? Insights from a South Asian mimicry group

Müllerian mimicry was proposed to be an example of a coevolved mutualism promoted by population isolation in glacial refugia. This, however, has not been well supported in butterfly models. Here, we use genomic data to test this theory while examining the population genetics behind mimetic diversification in a pair of co-mimetic bumble bees, Bombus breviceps Smith and Bombus trifasciatus Smith. In both lineages, populations were structured by geography but not as much by color pattern, suggesting sharing of color alleles across regions of restricted gene flow and formation of mimicry complexes in the absence of genetic differentiation. Demographic analyses showed mismatches between historical effective population size changes and glacial cycles, and niche modeling revealed only mild habitat retraction during glaciation. Moreover, mimetic subpopulations of the same color form in the two lineages only in some cases exhibit similar population history and genetic divergence. Therefore, the current study supports a more complex history in this comimicry than a simple refugium-coevolution model.

Limited introgression from non-native commercial strains and signatures of adaptation in the key pollinator Bombus terrestris

Insect pollination is fundamental for natural ecosystems and agricultural crops. The bumblebee species Bombus terrestris has become a popular choice for commercial crop pollination worldwide due to its effectiveness and ease of mass rearing. Bumblebee colonies are mass produced for the pollination of more than 20 crops and imported into over 50 countries including countries outside their native ranges, and the risk of invasion by commercial non-native bumblebees is considered an emerging issue for global conservation and biological diversity. Here, we use genome-wide data from seven wild populations close to and far from farms using commercial colonies, as well as commercial populations, to investigate the implications of utilizing commercial bumblebee subspecies in the UK. We find evidence for generally low levels of introgression between commercial and wild bees, with higher admixture proportions in the bees occurring close to farms. We identify genomic regions putatively involved in local and global adaptation, and genes in locally adaptive regions were found to be enriched for functions related to taste receptor activity, oxidoreductase activity, fatty acid and lipid biosynthetic processes. Despite more than 30 years of bumblebee colony importation into the UK, we observe low impact on the genetic integrity of local B. terrestris populations, but we highlight that even limited introgression might negatively affect locally adapted populations.

Genomic divergence and a lack of recent introgression between commercial and wild bumblebees (Bombus terrestris)

The global movement of bees for agricultural pollination services can affect local pollinator populations via hybridization. When commercial bumblebees are of the same species but of different geographic origin, intraspecific hybridization may result in beneficial integration of new genetic variation, or alternatively may disrupt locally adapted gene complexes. However, neither the existence nor the extent of genomic introgression and evolutionary divergence between wild and commercial bumblebees is fully understood. We obtained whole-genome sequencing data from wild and commercial Bombus terrestris collected from sites in Southern Sweden with and without long-term use of commercially imported B. terrestris. We search for evidence of introgression, dispersal and genome-wide differentiation in a comparative genomic analysis of wild and commercial bumblebees. Commercial B. terrestris were found in natural environments near sites where commercial bumblebees were used, as well as drifting wild B. terrestris in commercial bumblebee colonies. However, we found no evidence for widespread, recent genomic introgression of commercial B. terrestris into local wild conspecific populations. We found that wild B. terrestris had significantly higher nucleotide diversity (Nei's pi, π), while the number of segregating sites (Watterson's theta, θw) was higher in commercial B. terrestris. A highly divergent region on chromosome 11 was identified in commercial B. terrestris and found to be enriched with structural variants. The genes present in this region are involved in flight muscle contraction and structure and pathogen immune response, providing evidence for differing evolutionary processes operating in wild and commercial B. terrestris. We did not find evidence for recent introgression, suggesting that co-occurring commercial B. terrestris have not disrupted evolutionary processes in wild B. terrestris populations.

Genomic Signatures of Recent Adaptation in a Wild Bumblebee

Environmental changes threaten insect pollinators, creating risks for agriculture and ecosystem stability. Despite their importance, we know little about how wild insects respond to environmental pressures. To understand the genomic bases of adaptation in an ecologically important pollinator, we analyzed genomes of Bombus terrestris bumblebees collected across Great Britain. We reveal extensive genetic diversity within this population, and strong signatures of recent adaptation throughout the genome affecting key processes including neurobiology and wing development. We also discover unusual features of the genome, including a region containing 53 genes that lacks genetic diversity in many bee species, and a horizontal gene transfer from a Wolbachia bacteria. Overall, the genetic diversity we observe and how it is distributed throughout the genome and the population should support the resilience of this important pollinator species to ongoing and future selective pressures. Applying our approach to more species should help understand how they can differ in their adaptive potential, and to develop conservation strategies for those most at risk.

Environmental differences explain subtle yet detectable genetic structure in a widespread pollinator

BACKGROUND

The environment is a strong driver of genetic structure in many natural populations, yet often neglected in population genetic studies. This may be a particular problem in vagile species, where subtle structure cannot be explained by limitations to dispersal. Consequently, these species might falsely be considered quasi-panmictic and hence potentially mismanaged. A species this might apply to, is the buff-tailed bumble bee (Bombus terrestris), an economically important and widespread pollinator, which is considered to be quasi-panmictic at mainland continental scales. Here we aimed to (i) quantify genetic structure in 21+ populations of the buff-tailed bumble bee, sampled throughout two Eastern European countries, and (ii) analyse the degree to which structure is explained by environmental differences, habitat permeability and geographic distance. Using 12 microsatellite loci, we characterised populations of this species with Fst analyses, complemented by discriminant analysis of principal components and Bayesian clustering approaches. We then applied generalized dissimilarity modelling to simultaneously assess the informativeness of geographic distance, habitat permeability and environmental differences among populations in explaining divergence.

RESULTS

Genetic structure of the buff-tailed bumble bee quantified by means of Fst was subtle and not detected by Bayesian clustering. Discriminant analysis of principal components suggested insignificant but still noticeable structure that slightly exceeded estimates obtained through Fst analyses. As expected, geographic distance and habitat permeability were not informative in explaining the spatial pattern of genetic divergence. Yet, environmental variables related to temperature, vegetation and topography were highly informative, explaining between 33 and 39% of the genetic variation observed.

CONCLUSIONS

In contrast to previous studies reporting quasi-panmixia in continental populations of this species, we demonstrated the presence of subtle population structure related to environmental heterogeneity. Environmental data proved to be highly useful in unravelling the drivers of genetic structure in this vagile and opportunistic species. We highlight the potential of including these data to obtain a better understanding of population structure and the processes driving it in species considered to be quasi-panmictic.

Spatial and temporal patterns of genetic diversity in Bombus terrestris populations of the Iberian Peninsula and their conservation implications

The bumblebee Bombus terrestris is used worldwide for crop pollination. Despite its positive impact on crop yield, it has become a widespread threat to biodiversity due to its interactions with local bumblebee populations. Commercial subspecies introduced to the Iberian Peninsula since the 1990s without any regulation have colonized the environment, with evidence of naturalization and introgression with the endemic subspecies Bombus terrestris lusitanicus. We have used mitochondrial and nuclear genetic data to describe the current genetic diversity of the Iberian population and to estimate the expansion of commercial bumblebees. Samples from the natural distribution range of the commercial subspecies, the natural intergradation area between the two subspecies and from a period prior to the use of commercial colonies (i.e., before the 1990s) have been used for comparison. Our results show that the mitochondrial haplotype of the commercial breeds has spread throughout the territory, which, together with subtle changes observed in the nuclear genetic diversity of the populations, indicates that hybridization and consequent introgression are occurring in most of the peninsula. It is, therefore, necessary to improve the existing legislation concerning the management and exportation of commercial bumblebees to conserve locally adapted populations.

On the Trail of Spatial Patterns of Genetic Variation

The accurate determination of the spatial trends on the variability of a species’ gene pool is essential to elucidate the underlying demographic-evolutionary events, thus helping to unravel the microevolutionary history of the population under study. Herein we present a new software called GenoCline, mainly addressed to detect genetic clines from allele, haplotype, and genome-wide data. This program package allows identifying the geographic orientation of clinal genetic variation through a system of iterative rotation of a virtual coordinate axis. Besides, GenoCline can perform complementary analyses to explore the potential origin of the genetic clines observed, including spatial autocorrelation, isolation by distance, centroid method, multidimensional scaling and Sammon projection. Among the advantages of this software is the ease in data entry and potential interconnection with other programs. Genetic and geographic data can be entered in spreadsheet table formatting (.xls), whereas genome-wide data can be imported in Eigensoft format. Genetic frequencies can also be exported in a format compatible with other programs dealing with population genetic and evolutionary biology analyses. All illustrations of results are saved in.svg format so that there will be high quality and easily editable vectorial graphs available for the researcher. Being implemented in Java, GenoCline is highly portable, thus working in different operating systems.

Population genomics of Bombus terrestris reveals high but unstructured genetic diversity in a potential glacial refugium

Ongoing climate change is expected to cause an increase in temperature and a reduction of precipitation levels in the Mediterranean region, which might cause changes in many species distributions. These effects negatively influence species gene pools, decreasing genetic variability and adaptive potential. Here, we use mitochondrial DNA and RADseq to analyse population genetic structure and genetic diversity of the bumblebee species Bombus terrestris (subspecies Bombus terrestris lusitanicus), in the Iberian Peninsula. Although this subspecies shows a panmictic pattern of population structure across Iberia and beyond, we found differentiation between subspecies B. t. lusitanicus and B. t. africanus, probably caused by the existence of barriers to gene flow between Iberia and North Africa. Furthermore, the results revealed that the Iberian Peninsula harbours a large fraction of B. terrestris intraspecific genetic variation, with the highest number of mitochondrial haplotypes found when compared with any other region in Europe studied so far, suggesting a potential role for the Iberian Peninsula as a glacial refugium. Our findings strengthen the idea that Iberia is a very important source of diversity for the global genetic pool of this species, because rare alleles might play a role in population resilience against human- or climate-mediated changes.

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.

Quaternary climate instability is correlated with patterns of population genetic variability in Bombus huntii

Climate oscillations have left a significant impact on the patterns of genetic diversity observed in numerous taxa. In this study, we examine the effect of Quaternary climate instability on population genetic variability of a bumble bee pollinator species, Bombus huntii in western North America. Pleistocene and contemporary B. huntii habitat suitability (HS) was estimated with an environmental niche model (ENM) by associating 1,035 locality records with 10 bioclimatic variables. To estimate genetic variability, we genotyped 380 individuals from 33 localities at 13 microsatellite loci. Bayesian inference was used to examine population structure with and without a priori specification of geographic locality. We compared isolation by distance (IBD) and isolation by resistance (IBR) models to examine population differentiation within and among the Bayesian inferred genetic clusters. Furthermore, we tested for the effect of environmental niche stability (ENS) on population genetic diversity with linear regression. As predicted, high-latitude B. huntii habitats exhibit low ENS when compared to low-latitude habitats. Two major genetic clusters of B. huntii inhabit western North America: (a) a north genetic cluster predominantly distributed north of 28°N and (b) a south genetic cluster distributed south of 28°N. In the south genetic cluser, both IBD and IBR models are significant. However, in the north genetic cluster, IBD is significant but not IBR. Furthermore, the IBR models suggest that low-latitude montane populations are surrounded by habitat with low HS, possibly limiting dispersal, and ultimately gene flow between populations. Finally, we detected high genetic diversity across populations in regions that have been climatically unstable since the last glacial maximum (LGM), and low genetic diversity across populations in regions that have been climatically stable since the LGM. Understanding how species have responded to climate change has the potential to inform management and conservation decisions of both ecological and economic concerns.

Genome-wide single nucleotide polymorphism scan suggests adaptation to urbanization in an important pollinator, the red-tailed bumblebee (Bombus lapidarius L.)

Urbanization is considered a global threat to biodiversity; the growth of cities results in an increase in impervious surfaces, soil and air pollution, fragmentation of natural vegetation and invasion of non-native species, along with numerous environmental changes, including the heat island phenomenon. The combination of these effects constitutes a challenge for both the survival and persistence of many native species, while also imposing altered selective regimes. Here, using 110 314 single nucleotide polymorphisms generated by restriction-site-associated DNA sequencing, we investigated the genome-wide effects of urbanization on putative neutral and adaptive genomic diversity in a major insect pollinator, Bombus lapidarius, collected from nine German cities and nine paired rural sites. Overall, genetic differentiation among sites was low and there was no obvious genome-wide genetic structuring, suggesting the absence of strong effects of urbanization on gene flow. We nevertheless identified several loci under directional selection, a subset of which was associated with urban land use, including the percentage of impervious surface surrounding each sampling site. Overall, our results provide evidence of local adaptation to urbanization in the face of gene flow in a highly mobile insect pollinator.

Development of Multiple Polymorphic Microsatellite Markers for Ceratina calcarata (Hymenoptera: Apidae) Using Genome-Wide Analysis

The small carpenter bee, Ceratina calcarata (Robertson), is a widespread native pollinator across eastern North America. The behavioral ecology and nesting biology of C. calcarata has been relatively well-studied and the species is emerging as a model organism for both native pollinator and social evolution research. C. calcarata is subsocial: reproductively mature females provide extended maternal care to their brood. As such, studies of C. calcarata may also reveal patterns of relatedness and demography unique to primitively social Hymenoptera. Here, we present 21 microsatellite loci, isolated from the recently completed C. calcarata genome. Screening in 39 individuals across their distribution revealed that no loci were in linkage disequilibrium, nor did any deviate significantly from Hardy-Weinberg following sequential Bonferroni correction. Allele count ranged from 2 to 14, and observed and expected heterozygosities ranged from 0.08 to 0.82 (mean 0.47) and 0.26 to 0.88 (mean 0.56), respectively. These markers will enable studies of population-wide genetic structuring across C. calcarata's distribution. Such tools will also allow for exploration of between and within-colony relatedness in this subsocial native pollinator.

Population structure of honey bees in the Carpathian Basin (Hungary) confirms introgression from surrounding subspecies

Carniolan honey bees (Apis mellifera carnica) are considered as an indigenous subspecies in Hungary adapted to most of the ecological and climatic conditions in this area. However, during the last decades Hungarian beekeepers have recognized morphological signs of the Italian honey bee (Apis mellifera ligustica). As the natural distribution of the honey bee subspecies can be affected by the importation of honey bee queens or by natural gene flow, we aimed at determining the genetic structure and characteristics of the local honey bee population using molecular markers. All together, 48 Hungarian and 84 foreign (Italian, Polish, Spanish, Liberian) pupae and/or workers were used for mitochondrial DNA analysis. Additionally, 53 sequences corresponding to 10 subspecies and the Buckfast hybrid were downloaded from GenBank. For the nuclear analysis, 236 Hungarian and 106 foreign honey bees were genotyped using nine microsatellites. Heterozygosity values, population-specific alleles, FST values, principal coordinate analysis, assignment tests, structure analysis, and dendrograms were calculated. Haplotype and nucleotide diversity values showed moderate values. We found that one haplotype (H9) was dominant in Hungary. The presence of the black honey bee (Apis mellifera mellifera) was negligible, but a few individuals resembling other subspecies were identified. We proved that the Hungarian honey bee population is nearly homogeneous but also demonstrated introgression from the foreign subspecies. Both mitochondrial DNA and microsatellite analyses corroborated the observations of the beekeepers. Molecular analyses suggested that Carniolan honey bee in Hungary is slightly affected by Italian and black honey bee introgression. Genetic differences were detected between Polish and Hungarian Carniolan honey bee populations, suggesting the existence of at least two different gene pools within A. m. carnica.