Find research datasets worth reusing

Data from: Population structure of mountain pine beetle symbiont Leptographium longiclavatum and the implication on the multipartite beetle-fungi relationships

Over 18 million ha of forests have been destroyed in the past decade in Canada by the mountain pine beetle (MPB) and its fungal symbionts. Understanding their population dynamics is critical to improving modeling of beetle epidemics and providing potential clues to predict population expansion. Leptographium longiclavatum and Grosmannia clavigera are fungal symbionts of MPB that aid the beetle to colonize and kill their pine hosts. We investigated the genetic structure and demographic expansion of L. longiclavatum in populations established within the historic distribution range and in the newly colonized regions. We identified three genetic clusters/populations that coincide with independent geographic locations. The genetic profiles of the recently established populations in northern British Columbia (BC) and Alberta suggest that they originated from central and southern BC. Approximate Bayesian Computation supports the scenario that this recent expansion represents an admixture of individuals originating from BC and the Rocky Mountains. Highly significant correlations were found among genetic distance matrices of L. longiclavatum, G. clavigera, and MPB. This highlights the concordance of demographic processes in these interacting organisms sharing a highly specialized niche and supports the hypothesis of long-term multipartite beetle-fungus co-evolutionary history and mutualistic relationships.

Data from: Divergent population structure and climate associations of a chromosomal inversion polymorphism across the Mimulus guttatus species complex

Chromosomal rearrangement polymorphisms are common and increasingly found to be associated with adaptive ecological divergence and speciation. Rearrangements, such as inversions, reduce recombination in heterozygous individuals and thus can protect favorable allelic combinations at linked loci, facilitating their spread in the presence of gene flow. Recently, we identified a chromosomal inversion polymorphism that contributes to ecological adaptation and reproductive isolation between annual and perennial ecotypes of the yellow monkeyflower, Mimulus guttatus. Here we evaluate the population genetic structure of this inverted region in comparison with the collinear regions of the genome across the M. guttatus species complex. We tested whether annual and perennial M. guttatus exhibit different patterns of divergence for loci in the inverted and noninverted regions of the genome. We then evaluated whether there are contrasting climate associations with these genomic regions through redundancy analysis. We found that the inversion exhibits broadly different patterns of divergence among annual and perennial M. guttatus and is associated with environmental variation across population accessions. This study is the first widespread population genetic survey of the diversity of the M. guttatus species complex. Our findings contribute to a greater understanding of morphological, ecological, and genetic evolutionary divergence across this highly diverse group of closely related ecotypes and species. Finally, understanding species relationships among M. guttatus sp. has hitherto been stymied by accumulated evidence of substantial gene flow among populations as well as designated species. Nevertheless, our results shed light on these relationships and provide insight into adaptation in life history traits within the complex.

opencc-zeroDec 2017View details →

Data from: Structural design of the minute clypeasteroid echinoid Echinocyamus pusillus

The clypeasteroid echinoid skeleton is a multi-plated, light-weight shell construction produced by biomineralization processes. In shell constructions, joints between individual elements are considered as weak points, yet, these echinoid skeleton shows an extensive preservation potential in both Recent and fossil environments. The remarkable strength of the test is achieved by skeletal reinforcement structures and their constructional layouts. Micro-computed tomography and scanning electron microscopy are used for micro-structural and volumetric analyses of the echinoid's skeleton. It is shown, that strengthening mechanisms act on different hierarchical levels from the overall shape of the skeleton to micro-skeletal interlocking. The tight-fitting and interlocking plate joints lead to a shell considered to behave as a monolithic structure. The plate's architecture features distinct regions interpreted as a significant load transferring system. The internal support system follows the segmentation of the remaining skeleton, where sutural layout and stereom distribution are designed for effective load transfer. The structural analysis of the multi-plated, yet monolithic skeleton of Echinocyamus pusillus reveals new aspects of the micro-morphology and its structural relevance for the load-bearing behaviour. The analysed structural principles allow E. pusillus to be considered as a role model for the development of multi-element, light-weight shell constructions.

Data from: Fine-scale population genetic structure of the Bengal tiger (Panthera tigris tigris) in a human-dominated western Terai Arc Landscape, India

Despite massive global conservation strategies, tiger populations continued to decline until recently, mainly due to habitat loss, human-animal conflicts, and poaching. These factors are known to affect the genetic characteristics of tiger populations and decrease local effective population sizes. The Terai Arc Landscape (TAL) at the foothills of the Himalaya is one of the 42 source sites of tigers around the globe. Therefore, information on how landscape features and anthropogenic factors affect the fine-scale spatial genetic structure and variation of tigers in TAL is needed to develop proper management strategies for achieving long-term conservation goals. We document, for the first time, the genetic characteristics of this tiger population by genotyping 71 tiger samples using 13 microsatellite markers from the western region of TAL (WTAL) (1800 km2). Specifically, we aimed to estimate the genetic variability, population structure, and gene flow. The microsatellite markers indicated that the levels of allelic diversity (MNA = 6.6) and genetic variation (Ho =0.50, HE = 0.64) were slightly lower than those reported previously in other Bengal tiger populations. We observed moderate gene flow and significant genetic differentiation (FST= 0.060), and identified the presence of cryptic genetic structure using Bayesian and non-Bayesian approaches. There was low and significantly asymmetric migration between the two main subpopulations of the Rajaji Tiger Reserve and the Corbett Tiger Reserve in WTAL. Sibship relationships indicated that the functionality of the corridor between these subpopulations may be retained if the quality of the habitat does not deteriorate. However, we found that gene flow is not adequate in view of changing land use matrices. We discuss the need to maintain connectivity by implementing the measures that have been suggested previously to minimize the level of human disturbance, including relocation of villages and industries, prevention of encroachment, and banning sand and boulder mining in the corridors.

Data from: Past climate change drives current genetic structure of an endangered freshwater mussel species

Historical-to-recent climate change and anthropogenic disturbance affect species distributions and genetic structure. The Rio Grande watershed of the United States and Mexico encompasses ecosystems that are intensively exploited, resulting in substantial degradation of aquatic habitats. While significant anthropogenic disturbances in the Rio Grande are recent, inhospitable conditions for freshwater organisms likely existed prior to such disturbances. A combination of anthropogenic and past climate factors may contribute to current distributions of aquatic fauna in the Rio Grande basin. We used mitochondrial DNA and 18 microsatellite loci to infer evolutionary history and genetic structure of an endangered freshwater mussel, Popenaias popeii, throughout the Rio Grande drainage. We estimated spatial connectivity and gene flow across extant populations of P. popeii and used ecological niche models (ENMs) and approximate Bayesian computation (ABC) to infer its evolutionary history during the Pleistocene. structure results recovered regional and local population clusters in the Rio Grande. ENMs predicted drastic reductions in suitable habitat during the last glacial maximum. ABC analyses suggested that regional population structure likely arose in this species during the mid-to-late Pleistocene and was followed by a late Pleistocene population bottleneck in New Mexico populations. The local population structure arose relatively recently, perhaps due to anthropogenic factors. Popenaias popeii, one of the few freshwater mussel species native to the Rio Grande basin, is a case study for understanding how both geological and anthropogenic factors shape current population genetic structure. Conservation strategies for this species should account for the fragmented nature of contemporary populations.

opencc-zeroJul 2019View details →

Data from: Pattern of population structuring between Belgian and Estonian bumblebees

Several population genetic studies investigated the extent of gene flow and population connectivity in bumblebees. In general, no restriction in gene flow is considered for mainland populations of common bumblebee species. Whether this assumption holds true for all species is not known. An assessment of bumblebee genetic structure in the context of their geographic distribution is needed to prioritize conservation and management needs. Here, we conducted a genetic study on seven bumblebee species occurring in Belgium and Estonia. Using 16 microsatellite markers, we investigated genetic diversity and population structuring in each species. This is the first study investigating population structuring of both declining and stable bumblebee species on both small and large geographic scales. Our results showed no or only low population structuring between the populations of the restricted and declining bumblebee species on both scales, while significant structuring was found for populations of the common species on the larger scale. The latter result, which may be due to human or environmental changes in the landscape, implies the need for the conservation of also widespread bumblebee species. Conservation strategies to improve gene flow and connectivity of populations could avoid the isolation and future losses of populations of these important species.

Data from: When viruses don't go viral: the importance of host phylogeographic structure in the spatial spread of arenaviruses

Many emerging infections are RNA virus spillovers from animal reservoirs. Reservoir identification is necessary for predicting the geographic extent of infection risk, but rarely are taxonomic levels below the animal species considered as reservoir, and only key circumstances in nature and methodology allow intrinsic virus-host associations to be distinguished from simple geographic (co-)isolation. We sampled and genetically characterized in detail a contact zone of two subtaxa of the rodent Mastomys natalensis in Tanzania. We find two distinct arenaviruses, Gairo and Morogoro virus, each spatially confined to a single M. natalensis subtaxon, only co-occurring at the contact zone's centre. Inter-subtaxon hybridization at this centre and a continuum of quality habitat for M. natalensis show that both viruses have the ecological opportunity to spread into the other substaxon's range, but do not, strongly suggesting host-intrinsic barriers. Such barriers could explain why human cases of another M. natalensis-borne arenavirus, Lassa virus, are limited to West Africa.

Data from: Fire does not strongly affect genetic diversity or structure of a common treefrog in the endangered Florida scrub

Fire regimes influence natural populations of organisms in diverse ways, via direct effects on population dynamics as well as indirect effects on habitat and ecosystem processes. Although many amphibian species have evolved to persist in fire-dependent ecosystems, the effects of fire on the genetic diversity of amphibian populations remain relatively unexplored. We examined how different aspects of fire history relate to population genetic diversity and structure of an abundant anuran, Hyla femoralis, in a large, intact area of Florida scrub containing hundreds of seasonally inundated ponds. Specifically, we assessed the overall population genetic structure and examined whether variation in time since fire, fire intensity, or historical fire frequency at breeding sites explained spatial variation in genetic diversity. Based on our sampling of 17 breeding aggregations within the 2,100-ha study area, neither recent nor frequent fire reduce genetic diversity or restrict connectivity among ponds for H. femoralis. Overall, mean effective population sizes were large (average range = 68–572). We detected a positive trend between effective population size (Ne) and average intensity of the most-recent fire, with this factor explaining 42% of the variation in Ne. Our results contrast with previous studies that consistently demonstrate strong relationships between fire history and population genetic structure of scrub-associated lizard species, suggesting that H. femoralis is resilient to a wide range of fire regimes. More generally, our study contributes to understanding the roles of life-history characteristics and environmental unpredictability in shaping organisms' responses to fire.

Data from: NetView P: a network visualization tool to unravel complex population structure using genome-wide SNPs

Network-based approaches are emerging as valuable tools for the analysis of complex genetic structure in both wild and captive populations. NetView P combines data quality control with the construction of population networks based on mutual k-nearest-neighbours thresholds applied to genome-wide SNPs. The program is cross-platform compatible, open-source and efficiently operates on data ranging from hundreds to hundreds of thousands of SNPs through multiprocessing in Python. We used the pipeline for the analysis of pedigree data from simulated (n = 750, SNPs = 1279) and captive Silver-lipped Pearl Oysters (n = 415, SNPs = 1107), wild populations of the European Hake from the Atlantic and Mediterranean (n = 834, SNPs = 380) and Gray Wolves from North America (n = 239, SNPs = 86,103). The population networks effectively visualize large- and fine-scale genetic structure within and between populations, including family-level structure and relationships. NetView P comprises a network-based addition to other population analysis tools and provides user-friendly access to a complex network analysis pipeline through implementation in Python.

opencc-zeroDec 2015View details →

Data from: Homogenous population genetic structure of the non-native raccoon dog (Nyctereutes procyonoides) in Europe as a result of rapid population expansion

The extent of gene flow during the range expansion of non-native species influences the amount of genetic diversity retained in expanding populations. Here, we analyse the population genetic structure of the raccoon dog (Nyctereutes procyonoides) in north-eastern and central Europe. This invasive species is of management concern because it is highly susceptible to fox rabies and an important secondary host of the virus. We hypothesized that the large number of introduced animals and the species' dispersal capabilities led to high population connectivity and maintenance of genetic diversity throughout the invaded range. We genotyped 332 tissue samples from seven European countries using 16 microsatellite loci. Different algorithms identified three genetic clusters corresponding to Finland, Denmark and a large 'central' population that reached from introduction areas in western Russia to northern Germany. Cluster assignments provided evidence of long-distance dispersal. The results of an Approximate Bayesian Computation analysis supported a scenario of equal effective population sizes among different pre-defined populations in the large central cluster. Our results are in line with strong gene flow and secondary admixture between neighbouring demes leading to reduced genetic structuring, probably a result of its fairly rapid population expansion after introduction. The results presented here are remarkable in the sense that we identified a homogenous genetic cluster inhabiting an area stretching over more than 1500km. They are also relevant for disease management, as in the event of a significant rabies outbreak, there is a great risk of a rapid virus spread among raccoon dog populations.

Data from: Influence of drift and admixture on population structure of American black bears (Ursus americanus) in the Central Interior Highlands, U.S.A. 50 years after translocation

Bottlenecks, founder events, and genetic drift often result in decreased genetic diversity and increased population differentiation. These events may follow abundance declines due to natural or anthropogenic perturbations, where translocations may be an effective conservation strategy to increase population size. American black bears (Ursus americanus) were nearly extirpated from the Central Interior Highlands, USA by 1920. In an effort to restore bears, 254 individuals were translocated from Minnesota, USA and Manitoba, Canada, into the Ouachita and Ozark Mountains from 1958 to 1968. Using 15 microsatellites and mitochondrial haplotypes, we observed contemporary genetic diversity and differentiation between the source and supplemented populations. We inferred four genetic clusters: Source, Ouachitas, Ozarks, and a cluster in Missouri where no individuals were translocated. Coalescent models using approximate Bayesian computation identified an admixture model as having the highest posterior probability (0.942) over models where the translocation was unsuccessful or acted as a founder event. Nuclear genetic diversity was highest in the source (A = 9.11) and significantly lower in the translocated populations (A = 7.07 - 7.34; P = 0.004). The Missouri cluster had the lowest genetic diversity (A = 5.48) and served as a natural experiment showing the utility of translocations to increase genetic diversity following demographic bottlenecks. Differentiation was greater between the two admixed populations than between the source, suggesting that genetic drift acted strongly over the eight generations since the translocation. The Ouachitas and Missouri were previously hypothesized to be remnant lineages. We observed a pre-translocation remnant signature in Missouri but not in the Ouachitas.

opencc-zeroDec 2015View details →

Data from: Predator-guided sampling reveals biotic structure in the bathypelagic

We targeted habitat used differentially by deep-diving, air-breathing predators to empirically sample their prey's distributions off southern California. Fine-scale measurements of the spatial variability of potential prey animals from the surface to 1200 m were obtained using conventional fisheries echosounders aboard a surface ship and uniquely integrated into a deep-diving autonomous vehicle. Significant spatial variability in the size, composition, total biomass, and spatial organization of biota was evident over all spatial scales examined and was consistent with the general distribution patterns of foraging Cuvier's beaked whales (Ziphius cavirostris) observed in separate studies. Striking differences were found in prey characteristics between regions at depth, however, did not reflect differences observed in surface layers. These differences in deep pelagic structure horizontally and relative to surface structure, absent clear physical differences, change our long-held views of this habitat as uniform. The revelation that animals deep in the water column are so spatially heterogeneous at scales from 10 m to 50 km critically affects our understanding of the processes driving predator-prey interactions, energy transfer, biogeochemical cycling and other ecological processes in the deep sea, and the connections between the productive surface mixed layer and the deep water column.

opencc-zeroJun 2019View details →

Data from: Landscape structure mediates zoochorous-dispersed seed rain under isolated pasture trees across distinct tropical regions

Context: Isolated pasture trees play an important role in forest recovery within fragmented tropical landscapes by attracting seed dispersers and facilitating seedling growth. However, studies with conflicting results have led to confusion about what drives variation in zoochorous-dispersed seed rain patterns under isolated tree canopies. Objectives: To assess the role of landscape and biological factors impacting zoochorous-dispersed seed rain under isolated pasture trees across three tropical regions of the world. Methods: We measured seed dispersal under 144 isolated pasture trees found in 12 fragmented tropical and sub-tropical landscapes in Australia, Colombia and Nigeria. Using linear mixed effect models, we modeled seed diversity, abundance, richness and evenness as functions of the biological features and landscape context of isolated trees. Results: Throughout all regions, the amount of woody vegetation surrounding trees in pastures was negatively related to rainforest seed diversity, evenness and abundance. Seed diversity and evenness increased significantly with the distance of isolated trees to forest fragments in the Australian sub-tropics, but elsewhere, seed diversity and evenness tended to decline with distance to forest, though not significantly. Conclusions: Our results suggest that the tree composition of landscapes surrounding isolated pasture trees is important for influencing zoochorous-dispersed seed rain, regardless of the region studied. Our study highlights the prominent role of landscape-scale, rather than local-scale factors on seed dispersal to isolated pasture trees, while providing strong evidence that early stage successional processes involving isolated pasture trees are similar throughout global tropical regions.

Data from: Structure of the epiphyte community in a tropical montane forest in SW China

Vascular epiphytes are an understudied and particularly important component of tropical forest ecosystems. However, owing to the difficulties of access, little is known about the properties of epiphyte-host tree communities and the factors structuring them, especially in Asia. We investigated factors structuring the vascular epiphyte-host community and its network properties in a tropical montane forest in Xishuangbanna, SW China. Vascular epiphytes were surveyed in six plots located in mature forests. Six host and four micro-site environmental factors were investigated. Epiphyte diversity was strongly correlated with host size (DBH, diameter at breast height), while within hosts the highest epiphyte diversity was in the middle canopy and epiphyte diversity was significantly higher in sites with canopy soil or a moss mat than on bare bark. DBH, elevation and stem height explained 22% of the total variation in the epiphyte species assemblage among hosts, and DBH was the most important factor which alone explained 6% of the variation. Within hosts, 51% of the variation in epiphyte assemblage composition was explained by canopy position and substrate, and the most important single factor was substrate which accounted for 16% of the variation. Analysis of network properties indicated that the epiphyte host community was highly nested, with a low level of epiphyte specialization, and an almost even interaction strength between epiphytes and host trees. Together, these results indicate that large trees harbor a substantial proportion of the epiphyte community in this forest.

opencc-zeroDec 2018View details →

Data from: Accounting for interspecific competition and age structure in demographic analyses of density dependence improves predictions of fluctuations in population size

Understanding species coexistence has long been a major goal of ecology. Coexistence theory for two competing species posits that intraspecific density dependence should be stronger than interspecific density dependence. Great tits and blue tits are two bird species that compete for food resources and nesting cavities. Based on long-term monitoring of these two competing species at sites across Europe, combining observational and manipulative approaches, we show that the strength of density regulation is similar for both species, and that individuals have contrasting abilities to compete depending on their age. For great tits, density regulation is driven mainly by intraspecific competition. In contrast, for blue tits, interspecific competition contributes as much as intraspecific competition, consistent with asymmetric competition between the two species. In addition, including age-specific effects of intra- and interspecific competition in density-dependence models improves predictions of fluctuations in population size by up to three times.

Data from: Range-wide analysis of genetic structure in a widespread, highly mobile species (Odocoileus hemionus) reveals the importance of historical biogeography

Highly mobile species that thrive in a wide range of habitats are expected to show little genetic differentiation across their range. A limited but growing number of studies have revealed that patterns of broad-scale genetic differentiation can and do emerge in vagile, continuously distributed species. However, these patterns are complex and often shaped by both historical and ecological factors. Comprehensive surveys of genetic variation at a broad scale and at high resolution are useful for detecting cryptic spatial genetic structure, and for investigating the relative roles of historical and ecological processes in structuring widespread, highly mobile species. In this study, we analyzed 10 microsatellite loci from over 1,900 samples collected across the full range of mule deer (Odocoileus hemionus), one of the most widely distributed and abundant of all large mammal species in North America. Through both individual- and population-based analyses we found evidence for three main genetic lineages, one corresponding to the 'mule deer' morphological type and two to the 'black-tailed deer' type. Historical biogeographic events likely are the primary drivers of genetic divergence in this species; boundaries of the three lineages correspond well with predictions based on Pleistocene glacial cycles and substructure within each lineage demonstrates island vicariance. However, across large geographic areas, including the entire mule deer lineage, we found that genetic variation fit an isolation-by-distance pattern rather than discrete clusters. A lack of genetic structure across wide geographic areas of the continental west indicates that ecological processes have not resulted in restrictions to gene flow sufficient for spatial genetic structure to emerge. Our results have important implications for our understanding of evolutionary mechanisms of divergence, as well as for taxonomy, conservation, and management.

Data from: Comparative landscape genetics of pond-breeding amphibians in Mediterranean temporal wetlands: the positive role of structural heterogeneity in promoting gene flow

Comparative landscape genetics studies can provide key information to implement cost-effective conservation measures favoring a broad set of taxa. These studies are scarce, particularly in Mediterranean areas, which include diverse but threatened biological communities. Here we focus on Mediterranean wetlands in central Iberia and perform a multi-level, comparative study of two endemic pond-breeding amphibians, a salamander (Pleurodeles waltl) and a toad (Pelobates cultripes). We genotyped 411 salamanders from 20 populations and 306 toads from 16 populations at 18 and 16 microsatellite loci, respectively, and identified major factors associated with population connectivity through the analysis of three sets of variables potentially affecting gene flow at increasingly finer levels of spatial resolution. Topographic, land use/cover, and remotely sensed vegetation/moisture indices were used to derive optimized resistance surfaces for the two species. We found contrasting patterns of genetic structure, with stronger, finer-scale genetic differentiation in Pleurodeles waltl, and notable differences in the role of fine-scale patterns of heterogeneity in vegetation cover and water content in shaping patterns of regional genetic structure in the two species. Overall, our results suggest a positive role of structural heterogeneity in population connectivity in pond-breeding amphibians, with habitat patches of Mediterranean scrubland and open oak woodlands ("dehesas") facilitating gene flow. Our study highlights the usefulness of remotely sensed continuous variables of land cover, vegetation and water content (e.g., NDVI, NDMI) in conservation-oriented studies aimed at identifying major drivers of population connectivity.

opencc-zeroDec 2012View details →

Data from: Geographic population structure of the African malaria vector Anopheles gambiae suggests a role for the forest-savannah biome transition as a barrier to gene flow

The primary Afrotropical malaria mosquito vector Anopheles gambiae sensu stricto has a complex population structure. In western Africa, this species is split into two molecular forms and displays local and regional variation in chromosomal arrangements and behaviours. To investigate patterns of macro-geographic population substructure, 25 An. gambiae samples from 12 African countries were genotyped at 13 microsatellite loci. This analysis detected the presence of additional population structuring, with the M-form being subdivided into distinct west, central and southern African genetic clusters. These clusters are coincident with the central African rainforest belt and northern and southern savannah biomes, which suggests restrictions to gene flow associated with the transition between these biomes. By contrast geographically patterned population substructure appears much weaker within the S-form.

opencc-zeroDec 2015View details →

Data from: Fine-scale spatial genetic structure in the frankincense tree Boswellia papyrifera (Del.) Hochst. and implications for conservation

The fine-scale genetic structure and how it varies between generations depends on the spatial scale of gene dispersal and other fundamental aspects of species' biology, such as the mating system. Such knowledge is crucial for the design of genetic conservation strategies. This is particularly relevant for species that are increasingly fragmented such as Boswellia papyrifera. This species occurs in dry tropical forests from Ethiopia, Eritrea and Sudan and is an important source of frankincense, a highly valued aromatic resin obtained from the bark of the tree. This study assessed the genetic diversity and fine-scale spatial genetic structure (FSGS) of two cohorts (adults and seedlings) from two populations (Guba-Arenja and Kurmuk) in Western Ethiopia and inferred intra-population gene dispersal in the species, using microsatellite markers. The expected heterozygosity (HE) was 0.664–0.724. The spatial analyses based on kinship coefficient (Fij) revealed a significant positive genetic correlation up to a distance of 130 m. Spatial genetic structure was relatively weak (Sp = 0.002–0.014) indicating that gene dispersal is extensive within the populations. Based on the FSGS patterns found, we estimate indirectly gene dispersal distances of 103 and 124 m for the two populations studied. The high heterozygosity, the low fixation index and the low Sp values found in this study are consistent with outcrossing as the (predominant) mating system in B. papyrifera. We suggest that seed collection for ex situ conservation and reforestation programmes of B. papyrifera should use trees separated by distances of at least 100 m but preferably 150 m to limit genetic relatedness among seeds from different trees.

Data from: Population genetic structure of Bombus terrestris in Europe: isolation and genetic differentiation of Irish and British populations

The genetic structure of the earth bumblebee (Bombus terrestris L.) was examined across 22 wild populations and two commercially reared populations using eight microsatellite loci and two mitochondrial genes. Our study included wild bumblebee samples from six populations in Ireland, one from the Isle of Man, four from Britain and 11 from mainland Europe. A further sample was acquired from New Zealand. Observed levels of genetic variability and heterozygosity were low in Ireland and the Isle of Man, but relatively high in continental Europe and among commercial populations. Estimates of Fst revealed significant genetic differentiation among populations. Bayesian cluster analysis indicated that Irish populations were highly differentiated from British and continental populations, the latter two showing higher levels of admixture. The data suggest that the Irish Sea and prevailing south westerly winds act as a considerable geographical barrier to gene flow between populations in Ireland and Britain; however, some immigration from the Isle of Man to Ireland was detected. The results are discussed in the context of the recent commercialization of bumblebees for the European horticultural industry.