Genetic Structure of Asian Populations of Bombus ignitus (Hymenoptera: Apidae)

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.

Bee and floral traits affect the characteristics of the vibrations experienced by flowers during buzz pollination

During buzz pollination, bees use their indirect flight muscles to produce vibrations that are transmitted to the flowers and result in pollen release. Although buzz pollination has been known for >100 years, we are still in the early stages of understanding how bee and floral characteristics affect the production and transmission of floral vibrations. Here, we analysed floral vibrations produced by four closely related bumblebee taxa (Bombus spp.) on two buzz-pollinated plants species (Solanum spp.). We measured floral vibrations transmitted to the flower to establish the extent to which the mechanical properties of floral vibrations depend on bee and plant characteristics. By comparing four bee taxa visiting the same plant species, we found that peak acceleration, root mean-squared acceleration (RMS) and frequency vary between bee taxa, but that neither bee size (intertegular distance) nor flower biomass (dry mass) affects peak acceleration, RMS or frequency. A comparison of floral vibrations of two bee taxa visiting flowers of two plant species showed that, while bee species affects peak acceleration, RMS and frequency, plant species only affects acceleration (peak acceleration and RMS), not frequency. When accounting for differences in the transmission of vibrations across the two types of flower, using a species-specific 'coupling factor', we found that RMS acceleration and peak displacement do not differ between plant species. This suggests that bees produce the same initial acceleration in different plants but that transmission of these vibrations through the flower is affected by floral characteristics.

Reproductive Interference in an Introduced Bumblebee: Polyandry may Mitigate Negative Reproductive Impact

As a signature of reproductive interference (RI), we reviewed hybrid production in eusocial bumblebees in Japan, by comparing introduced Bombus terrestris with native B. ignitus in Honshu (main island of Japan) and with native B. hypocrita sapporoensis in Hokkaido (northern island of Japan). In this review, we present additional new data showing hybrid production between introduced B. terrestris and native B. ignitus in Honshu. Interspecific mating with introduced B. terrestris disrupts the reproduction of native B. h. sapporoensis and B. ignitus, which belong to the same subgenus of Bombus, through inviable egg production. This interference appears to facilitate species replacement on Hokkaido. Simultaneously, the mating frequencies for queens of B. terrestris have increased, suggesting that polyandry might evolve in response to the extent of RI between B. terrestris and B. h. sapporoensis. To suppress the population size of B. terrestris in Hokkaido, two methods have been proposed: the mass release of B. h. sapporoensis males to induce RI between the two species and the spraying of insecticides against foraging workers so that the workers will carry the insecticides back to their colonies, killing the immature bees within the colonies. A candidate insecticide type is insect growth regulator, which may disrupt larval development without any apparent effect on foraging workers.

Population structuring of the ubiquitous stingless bee Tetragonisca angustula in southern Brazil as revealed by microsatellite and mitochondrial markers

Tetragonisca angustula is one of the most widespread stingless bees in the Neotropics. This species swarms frequently and is extremely successful in urban environments. In addition, it is one of the most popular stingless bee species for beekeeping in Latin America, so nest transportation and trading is common. Nest transportation can change the genetic structure of the host population, reducing inbreeding and increasing homogenization. Here, we evaluate the genetic structure of 17 geographic populations of T. angustula in southern Brazil to quantify the level of genetic differentiation between populations. Analyses were conducted on partially sequenced mitochondrial genes and 11 microsatellite loci of 1002 workers from 457 sites distributed on the mainland and on 3 islands. Our results show that T. angustula populations are highly differentiated as demonstrated by mitochondrial DNA (mtDNA) and microsatellite markers. Of 73 haplotypes, 67 were population-specific. MtDNA diversity was low in 9 populations but microsatellite diversity was moderate to high in all populations. Microsatellite data suggest 10 genetic clusters and low level of gene flow throughout the studied area. However, physical barriers, such as rivers and mountain ranges, or the presence or absence of forest appear to be unrelated to population clusters. Factors such as low dispersal, different ecological conditions, and isolation by distance are most likely shaping the population structure of this species. Thus far, nest transportation has not influenced the general population structure in the studied area. However, due to the genetic structure we found, we recommend that nest transportation should only occur within and between populations that are genetically similar.

Comparative phylogeography in the Atlantic forest and Brazilian savannas: pleistocene fluctuations and dispersal shape spatial patterns in two bumblebees

Background

Bombus morio and B. pauloensis are sympatric widespread bumblebee species that occupy two major Brazilian biomes, the Atlantic forest and the savannas of the Cerrado. Differences in dispersion capacity, which is greater in B. morio, likely influence their phylogeographic patterns. This study asks which processes best explain the patterns of genetic variation observed in B. morio and B. pauloensis, shedding light on the phenomena that shaped the range of local populations and the spatial distribution of intra-specific lineages.

Results

Results suggest that Pleistocene climatic oscillations directly influenced the population structure of both species. Correlative species distribution models predict that the warmer conditions of the Last Interglacial contributed to population contraction, while demographic expansion happened during the Last Glacial Maximum. These results are consistent with physiological data suggesting that bumblebees are well adapted to colder conditions. Intra-specific mitochondrial genealogies are not congruent between the two species, which may be explained by their documented differences in dispersal ability.

Conclusions

While populations of the high-dispersal B. morio are morphologically and genetically homogeneous across the species range, B. pauloensis encompasses multiple (three) mitochondrial lineages, and show clear genetic, geographic, and morphological differences. Because the lineages of B. pauloensis are currently exposed to distinct climatic conditions (and elevations), parapatric diversification may occur within this taxon. The eastern portion of the state of São Paulo, the most urbanized area in Brazil, represents the center of genetic diversity for B. pauloensis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0803-0) contains supplementary material, which is available to authorized users.

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.

Hitting an Unintended Target: Phylogeography of Bombus brasiliensis Lepeletier, 1836 and the First New Brazilian Bumblebee Species in a Century (Hymenoptera: Apidae)

This work tested whether or not populations of Bombus brasiliensis isolated on mountain tops of southeastern Brazil belonged to the same species as populations widespread in lowland areas in the Atlantic coast and westward along the Paraná-river valley. Phylogeographic and population genetic analyses showed that those populations were all conspecific. However, they revealed a previously unrecognized, apparently rare, and potentially endangered species in one of the most threatened biodiversity hotspots of the World, the Brazilian Atlantic Forest. This species is described here as Bombus bahiensis sp. n., and included in a revised key for the identification of the bumblebee species known to occur in Brazil. Phylogenetic analyses based on two mtDNA markers suggest this new species to be sister to B. brasiliensis, from which its workers and queens can be easily distinguished by the lack of a yellow hair-band on the first metasomal tergum. The results presented here are consistent with the hypothesis that B. bahiensis sp. n. may have originated from an ancestral population isolated in an evergreen-forest refuge (the so-called Bahia refuge) during cold, dry periods of the Pleistocene. This refuge is also known as an important area of endemism for several animal taxa, including other bees. Secondary contact between B. bahiensis and B. brasiliensis may be presently prevented by a strip of semi-deciduous forest in a climate zone characterized by relatively long dry seasons. Considering the relatively limited range of this new species and the current anthropic pressure on its environment, attention should be given to its conservation status.

Patterns of genetic and reproductive traits differentiation in Mainland vs. Corsican populations of bumblebees

Populations on islands often exhibit lower levels of genetic variation and ecomorphological divergence compared to their mainland relatives. While phenotypic differentiation in characters, such as size or shape among insular organisms, has been well studied, insular differentiation in quantitative reproductive traits involved in chemical communication has received very little attention to date. Here, we investigated the impact of insularity on two syntopic bumblebee species pairs: one including species that are phylogenetically related (Bombus terrestris and B. lucorum), and the other including species that interact ecologically (B. terrestris and its specific nest inquiline B. vestalis). For each bumblebee species, we characterized the patterns of variation and differentiation of insular (Corsican) vs. mainland (European) populations (i) with four genes (nuclear and mitochondrial, 3781 bp) and (ii) in the chemical composition of male marking secretions (MMS), a key trait for mate attraction in bumblebees, by gas chromatography-mass spectrometry (GC-MS). Our results provide evidence for genetic differentiation in Corsican bumblebees and show that, contrary to theoretical expectations, island populations of bumblebees exhibit levels of genetic variation similar to the mainland populations. Likewise, our comparative chemical analyses of MMS indicate that Corsican populations of bumblebees are significantly differentiated from the mainland yet they hold comparative levels of within-population MMS variability compared to the mainland. Therefore, insularity has led Corsican populations to diverge both genetically and chemically from their mainland relatives, presumably through genetic drift, but without a decrease of genetic diversity in island populations. We hypothesize that MMS divergence in Corsican bumblebees was driven by a persistent lack of gene flow with mainland populations and reinforced by the preference of Corsican females for sympatric (Corsican) MMS. The impoverished Corsican bumblebee fauna has not led to relaxation of stabilizing selection on MMS but to consistent differentiation chemical reproductive traits on the island.

Microsatellite variation in two subspecies of cynomolgus monkeys (Macaca fascicularis)

To estimate the genetic variability of two subspecies of cynomolgus monkeys (Macaca fascicularis fascicularis and M. f. aurea) using microsatellite markers, 26 microsatellite markers were selected from previous reports. Seventeen markers showed high polymorphism in a subset of monkeys and were used for the assessment of genetic diversity in the larger sample. The effective number of alleles, the polymorphism information content (PIC) and the expected heterozygosity of M. f. aurea monkeys were all statistically significantly higher than those of M. f. fascicularis monkeys (P < 0.05), suggesting the M. f. aurea monkeys had a higher degree of genetic variation than the M. f. fascicularis monkeys. Substantial differences in allele distribution were also detected between the two subspecies of cynomolgus monkeys. Private alleles restricted to the M. f. fascicularis or the M. f. aurea monkeys were found throughout the selected 17 loci. These private alleles may allow the discrimination of the two subspecies of cynomolgus monkeys. The selected markers could also be used to estimate the genetic variation for other subspecies of cynomolgus monkeys. Further work using additional animals obtained from native or independent sources will be important for a more complete understanding of the genetic differences between these two subgroups.

Investigation of the genetic diversity of an invasive whitefly (Bemisia tabaci) in China using both mitochondrial and nuclear DNA markers

It is often considered that reduced genetic variation due to bottlenecks and founder effects limits the capacity for species to establish in new environments and subsequently spread. The recent invasion (during the past five years) of an alien whitefly, one member of Bemisia tabaci cryptic species complex, referred to as Mediterranean (herein referred to as Q-type) in Shandong Province, China, provides an ideal opportunity to study the changes in genetic variation between its home range in the Mediterranean region and its invasion range. Using both the mitochondrial cytochrome oxidase I (mtCOI) and nuclear (microsatellite) DNA, we show that Q in Shandong likely originated in the western Mediterranean. We also found that the haplotype diversity was low compared with its presumed geographic origin, whereas microsatellite allele diversity showed no such decline. A key factor in invasions is the establishment of females and so bottleneck and founder events can lead to a very rapid and considerable loss of mitochondrial diversity. The lack of haplotype diversity in Shandong supports the interpretation that, at one or more points between the western Mediterranean and China, the invading Q lost haplotype diversity, most probably through the serial process of establishment and redistribution through trade in ornamental plants. However, the loss in haplotype diversity does not necessarily mean that nuclear allelic diversity should also decline. Provided females can mate freely with whichever males are available, allelic diversity can be maintained or even increased relative to the origin of the invader. Our findings may offer some explanation to the apparent paradox between the concept of reduced genetic variation limiting adaptation to new environments and the observed low diversity in successful invaders.

Sequence determinants of human microsatellite variability

BACKGROUND

Microsatellite loci are frequently used in genomic studies of DNA sequence repeats and in population studies of genetic variability. To investigate the effect of sequence properties of microsatellites on their level of variability we have analyzed genotypes at 627 microsatellite loci in 1,048 worldwide individuals from the HGDP-CEPH cell line panel together with the DNA sequences of these microsatellites in the human RefSeq database.

RESULTS

Calibrating PCR fragment lengths in individual genotypes by using the RefSeq sequence enabled us to infer repeat number in the HGDP-CEPH dataset and to calculate the mean number of repeats (as opposed to the mean PCR fragment length), under the assumption that differences in PCR fragment length reflect differences in the numbers of repeats in the embedded repeat sequences. We find the mean and maximum numbers of repeats across individuals to be positively correlated with heterozygosity. The size and composition of the repeat unit of a microsatellite are also important factors in predicting heterozygosity, with tetra-nucleotide repeat units high in G/C content leading to higher heterozygosity. Finally, we find that microsatellites containing more separate sets of repeated motifs generally have higher heterozygosity.

CONCLUSIONS

These results suggest that sequence properties of microsatellites have a significant impact in determining the features of human microsatellite variability.

Cryptic differences in dispersal lead to differential sensitivity to habitat fragmentation in two bumblebee species

Habitat loss has led to fragmentation of populations of many invertebrates, but social hymenopterans may be particularly sensitive to habitat fragmentation due to their low effective population sizes. The impacts of fragmentation depend strongly on dispersal abilities, but these are difficult to quantify. Here, we quantify and compare dispersal abilities of two bumblebee species, Bombus muscorum and Bombus jonellus, in a model island system. We use microsatellites to investigate population genetic structuring, dispersal and spatial patterns in genetic diversity. Populations of both species showed significant structuring, and isolation by distance, but this was markedly greater in B. muscorum (theta = 0.13) than in B. jonellus (theta = 0.034). This difference could reflect a higher effective population size in B. jonellus compared to B. muscorum, but this is not consistent with the observed abundance of the two species. We argue that it is more likely that B. jonellus has a higher propensity to disperse than B. muscorum. This will influence their relative susceptibility to habitat fragmentation and may in part explain differential declines of mainland populations of these and other bumblebee species.

Complex patterns of genetic variability in populations of the New World screwworm fly revealed by mitochondrial DNA markers

Cochliomyia hominivorax (Coquerel) (Diptera: Calliphoridae), the New World screwworm fly, is an important agent of traumatic myiasis, which is endemic in the Neotropical region and which has great economic impact on the development of the livestock industry. International efforts have been aimed at designing programmes to control and eradicate this species from endemic areas. Thorough knowledge of the population genetics of an insect pest is a fundamental component to ensuring the success of a pest management strategy because it enables the determination of an appropriate geographic scale for carrying out effective treatments. This study undertook an analysis of mtDNA polymerase chain reaction-restricted fragment length polymorphism (PCR-RFLP) in 34 populations of C. hominivorax from 10 countries, encompassing almost all the current distribution of the species. Results showed high levels of mitochondrial DNA variability (pi= 2.9%) and a complex pattern of population genetic structure for this species. Significant population structure (Phi st= 0.5234) and low variability were found in Caribbean populations, suggesting that, in general, islands constitute independent evolutionary entities connected by restricted gene flow. By contrast, high variability and low, but significant, differentiation was found among mainland populations (Phi st= 0.0483), which could not be attributed to geographic distance. Several processes may be acting to maintain the observed patterns, with different implications for establishing control programmes.

Interspecific mating of the introduced bumblebee Bombus terrestris and the native Japanese bumblebee Bombus hypocrita sapporoensis results in inviable hybrids

The bumblebee Bombus terrestris is not only an effective pollinator, but also a potential invasive alien species outside its native range. Recently, nearly 30% of queens of the Japanese native species Bombus hypocrita sapporoensis and B. hypocrita hypocrita were estimated to copulate with B. terrestris males in the field, suggesting that indigenous bumblebees could be genetically deteriorated through hybrid production with the introduced species. In this study, we evaluated hybrid production between the introduced B. terrestris and the indigenous B. hypocrita sapporoensis under laboratory conditions. The hatching rate of eggs derived from interspecific matings was 0% and 8.6% depending on the direction of the cross, which was significantly lower than that from intraspecific matings of B. terrestris (76.9%) and B. hypocrita sapporoensis (78.9%). Genetic studies using microsatellite markers revealed that both haploid and diploid individuals were present in the egg stage, whereas all hatched larvae were haploid. In addition, histological studies revealed that eggs derived from interspecific matings terminated development 2 days after oviposition. These results strongly suggested that eggs derived from interspecific matings are inviable due to post-mating isolation mechanisms. Mass release of exotic pollinators could cause serious population declines of native bumblebee species.

Decline and conservation of bumble bees

Declines in bumble bee species in the past 60 years are well documented in Europe, where they are driven primarily by habitat loss and declines in floral abundance and diversity resulting from agricultural intensification. Impacts of habitat degradation and fragmentation are likely to be compounded by the social nature of bumble bees and their largely monogamous breeding system, which renders their effective population size low. Hence, populations are susceptible to stochastic extinction events and inbreeding. In North America, catastrophic declines of some bumble bee species since the 1990s are probably attributable to the accidental introduction of a nonnative parasite from Europe, a result of global trade in domesticated bumble bee colonies used for pollination of greenhouse crops. Given the importance of bumble bees as pollinators of crops and wildflowers, steps must be taken to prevent further declines. Suggested measures include tight regulation of commercial bumble bee use and targeted use of environmentally comparable schemes to enhance floristic diversity in agricultural landscapes.

Molecular evolution study in China: progress and future promise

China has a large land area with highly diverse topography, climate and vegetation, and animal resources and is ranked eighth in the world and first in the Northern Hemisphere on richness of biodiversity. Even though little work on molecular evolution had been reported a decade ago, studies on both the evolution of macromolecules and the molecular phylogeny have become active in China in recent years. This review highlights some of the interesting and important developments in molecular evolution study in China. Chinese scientists have made significant contribution on the methods inferring phylogeny and biogeography of animals and plants in East Asia using molecular data. Studies on population and conservation genetics of animals and plants, such as Golden monkey and Chinese sturgeon, provided useful information for conserving the endangered species. East and South Asia has been demonstrated to be one of the centres of domestication. Origin and evolution of genes and gene families have been explored, which shed new insight on the genetic mechanism of adaptation. In the genomic era, Chinese researchers also made a transition from single-gene to a genomic investigation approach. Considering the fact that amazing progress has been made in the past few years, and more and more talented young scientists are entering field, the future of molecular evolution study in China holds much promise.

Extremely low effective population sizes, genetic structuring and reduced genetic diversity in a threatened bumblebee species, Bombus sylvarum (Hymenoptera: Apidae)

Habitat fragmentation may severely affect survival of social insect populations as the number of nests per population, not the number of individuals, represents population size, hence they may be particularly prone to loss of genetic diversity. Erosion of genetic diversity may be particularly significant among social Hymenoptera such as bumblebees (Bombus spp.), as this group may be susceptible to diploid male production, a suggested direct cost of inbreeding. Here, for the first time, we assess genetic diversity and population structuring of a threatened bumblebee species (Bombus sylvarum) which exists in highly fragmented habitat (rather than oceanic) islands. Effective population sizes, estimated from identified sisterhoods, were very low (range 21-72) suggesting that isolated populations will be vulnerable to loss of genetic variation through drift. Evidence of significant genetic structuring between populations (theta = 0.084) was found, but evidence of a bottleneck was detected in only one population. Comparison across highly fragmented UK populations and a continental population (where this species is more widespread) revealed significant differences in allelic richness attributable to a high degree of genetic diversity in the continental population. While not directly related to population size, this is perhaps explained by the high degree of isolation between UK populations relative to continental populations. We suggest that populations now existing on isolated habitat islands were probably linked by stepping-stone populations prior to recent habitat loss.

Panmixia: an example from Dawson's burrowing bee (Amegilla dawsoni) (Hymenoptera: Anthophorini)

Dawson's burrowing bee is a large, fast-flying solitary nesting bee endemic to the arid zone of Western Australia. In this study the population structure of the species was examined with molecular markers. Using eight microsatellite loci, we genotyped 531 adult female bees collected from 13 populations of Dawson's burrowing bee, Amegilla dawsoni, across the species range. The mean number of alleles per locus ranged from 4 to 38 and expected heterozygosity was uniformly high with a mean of 0.602. Pairwise comparisons of F(ST) among all 13 populations ranged from 0.0071 to 0.0122 with only one significant estimate and an overall F(ST) of 0.001. The entire sample collection was in Hardy-Weinberg equilibrium and there was no evidence of inbreeding with a mean F(IS) of 0.010. The mating and nesting behaviour of this bee suggests that gene flow would be limited by monandry and the fact that almost 90% of females mate immediately on emergence. Nevertheless there is obviously sufficient gene flow to maintain panmixia, and we suggest that this results from infrequent and unreliable rainfall in the species range, which causes the bees to congregate at limited food resources, allowing a small number of unmated females from one emergence site to come into contact with males from another population. In addition, when drought eliminates food resources near an emergence site, the whole population may move elsewhere, increasing gene flow across the species range.

Population structure and inbreeding in a rare and declining bumblebee, Bombus muscorum (Hymenoptera: Apidae)

Owing to habitat loss populations of many organisms have declined and become fragmented. Vertebrate conservation strategies routinely consider genetic factors, but their importance in invertebrate populations is poorly understood. Bumblebees are important pollinators, and many species have undergone dramatic declines. As monoandrous social hymenopterans they may be particularly susceptible to inbreeding due to low effective population sizes. We study fragmented populations of a bumblebee species, on a model island system, and on mainland Great Britain where it is rare and declining. We use microsatellites to study: population genetic structuring and gene flow; the relationships between genetic diversity, population size and isolation; and frequencies of (sterile) diploid males - an indicator of inbreeding. We find significant genetic structuring (theta = 0.12) and isolation by distance. Populations > 10 km apart are all significantly differentiated, both on oceanic islands and on the mainland. Genetic diversity is reduced relative to closely related common species, and isolated populations exhibit further reductions. Of 16 populations, 10 show recent bottlenecking, and 3 show diploid male production. These results suggest that surviving populations of this rare insect suffer from inbreeding as a result of geographical isolation. Implications for the conservation of social hymenopterans are discussed.