M. Alice Pinto

Standard methods for molecular research in Apis mellifera

Summary From studies of behaviour, chemical communication, genomics and developmental biology, among many others, honey bees have long been a key organism for fundamental breakthroughs in biology. With a genome sequence in hand, and much improved genetic tools, honey bees are now an even more appealing target for answering the major questions of evolutionary biology, population structure, and social organization. At the same time, agricultural incentives to understand how honey bees fall prey to disease, or evade and survive their many pests and pathogens, have pushed for a genetic understanding of individual and social immunity in this species. Below we describe and reference tools for using modern molecular-biology techniques to understand bee behaviour, health, and other aspects of their biology. We focus on DNA and RNA techniques, largely because techniques for assessing bee proteins are covered in detail in Hartfelder et al. (2013). We cover practical needs for bee sampling, transport, and storage, and then discuss a range of current techniques for genetic analysis. We then provide a roadmap for genomic resources and methods for studying bees, followed by specific statistical protocols for population genetics, quantitative genetics, and phylogenetics. Finally, we end with three important tools for predicting gene regulation and function in honey bees: Fluorescence in situ hybridization (FISH), RNA interference (RNAi), and the estimation of chromosomal methylation and its role in epigenetic gene regulation.

Temporal pattern of africanization in a feral honeybee population from Texas inferred from mitochondrial DNA

Abstract The invasion of Africanized honeybees (Apis mellifera L.) in the Americas provides a window of opportunity to study the dynamics of secondary contact of subspecies of bees that evolved in allopatry in ecologically distinctive habitats of the Old World. We report here the results of an 11‐year mitochondrial DNA survey of a feral honeybee population from southern United States (Texas). The mitochondrial haplotype (mitotype) frequencies changed radically during the 11‐year study period. Prior to immigration of Africanized honeybees, the resident population was essentially of eastern and western European maternal ancestry. Three years after detection of the first Africanized swarm there was a mitotype turnover in the population from predominantly eastern European to predominantly A. m. scutellata (ancestor of Africanized honeybees). This remarkable change in the mitotype composition coincided with arrival of the parasitic mite Varroa destructor, which was likely responsible for severe losses experienced by colonies of European ancestry. From 1997 onward the population stabilized with most colonies of A. m. scutellata maternal origin.

Genetic integrity of the Dark European honey bee (Apis mellifera mellifera) from protected populations: a genome-wide assessment using SNPs and mtDNA sequence data

Summary The recognition that the Dark European honey bee, Apis mellifera mellifera, is increasingly threatened in its native range has led to the establishment of conservation programmes and protected areas throughout western Europe. Previous molecular surveys showed that, despite management strategies to preserve the genetic integrity of A. m. mellifera, protected populations had a measurable component of their gene pool derived from commercial C-lineage honey bees. Here we used both sequence data from the tRNAleu-cox2 intergenic mtDNA region and a genome-wide scan, with over 1183 single nucleotide polymorphisms (SNPs), to assess genetic diversity and introgression levels in several protected populations of A. m. mellifera, which were then compared with samples collected from unprotected populations. MtDNA analysis of the protected populations revealed a single colony bearing a foreign haplotype, whereas SNPs showed varying levels of introgression ranging from virtually zero in Norway to about 14% in Denmark. Introgression overall was higher in unprotected (30%) than in protected populations (8%), and is reflected in larger SNP diversity levels of the former, although opposite diversity levels were observed for mtDNA. These results suggest that, despite controlled breeding, some protected populations still require adjustments to the management strategies to further purge foreign alleles, which can be identified by SNPs.

Signatures of selection in the Iberian honey bee (Apis mellifera iberiensis) revealed by a genome scan analysis of single nucleotide polymorphisms.

Understanding the genetic mechanisms of adaptive population divergence is one of the most fundamental endeavours in evolutionary biology and is becoming increasingly important as it will allow predictions about how organisms will respond to global environmental crisis. This is particularly important for the honey bee, a species of unquestionable ecological and economical importance that has been exposed to increasing human‐mediated selection pressures. Here, we conducted a single nucleotide polymorphism (SNP)‐based genome scan in honey bees collected across an environmental gradient in Iberia and used four FST‐based outlier tests to identify genomic regions exhibiting signatures of selection. Additionally, we analysed associations between genetic and environmental data for the identification of factors that might be correlated or act as selective pressures. With these approaches, 4.4% (17 of 383) of outlier loci were cross‐validated by four FST‐based methods, and 8.9% (34 of 383) were cross‐validated by at least three methods. Of the 34 outliers, 15 were found to be strongly associated with one or more environmental variables. Further support for selection, provided by functional genomic information, was particularly compelling for SNP outliers mapped to different genes putatively involved in the same function such as vision, xenobiotic detoxification and innate immune response. This study enabled a more rigorous consideration of selection as the underlying cause of diversity patterns in Iberian honey bees, representing an important first step towards the identification of polymorphisms implicated in local adaptation and possibly in response to recent human‐mediated environmental changes.

Spatial and Temporal Distribution and Nest Site Characteristics of Feral Honey Bee (Hymenoptera: Apidae) Colonies in a Coastal Prairie Landscape

Abstract We evaluated the distribution and abundance of feral honey bee, Apis mellifera L., colonies in a coastal prairie landscape by examining nest site characteristics, population trends, and spatial and temporal patterns in cavity use. The colony densities of up to 12.5 colonies per km2 were the highest reported in the literature for an area including both suitable and unsuitable patches of nesting habitat. The measured cavity attributes were similar to those reported from other areas. The time occupied and turnover indices provided useful information about cavity quality, although none of the measured cavity attributes were correlated with these indices. Unmeasurable cavity characteristics, such as cavity volume, may provide a better estimate of cavity quality. Spatial patterns existed in cavity use by the feral colonies, with the colonies showing an aggregated pattern of distribution throughout the study. Colony aggregations probably resulted from the distribution of resources, especially cavities. Two years after the arrival of Africanized honey bees, cavities used by Africanized and European colonies were aggregated in distribution. During what seemed to be a transition period, both Africanized and European colonies were randomly distributed. After that time, European colonies remained randomly distributed, whereas Africanized colonies were aggregated. Therefore, the invasion of Africanized honey bees seemed to fragment the existing European population, corresponding to a decrease in the overall number of European colonies in the study area.

Identification of Africanized Honey Bee (Hymenoptera: Apidae) Mitochondrial DNA: Validation of a Rapid Polymerase Chain Reaction-Based Assay

Abstract Polymerase chain reaction (PCR)-amplified mitochondrial DNA (mtDNA) assays have been used in studies of the Africanization process in neotropical feral and managed honey bee populations. The approach has been adopted, in conjunction with morphometric analysis, to identify Africanized bees for regulatory purposes in the United States such as in California. In this study, 211 Old World colonies, representing all known introduced subspecies in the United States, and 451 colonies from non-Africanized areas of the southern United States were screened to validate a rapid PCR-based assay for identification of Africanized honey bee mtDNA. This PCR-based assay requires a single enzyme digestion (BglII) of a single PCR-amplified segment of the cytochrome b gene. The BglII polymorphism discriminates the mitochondrial haplotype (mitotype) of Apis mellifera scutellata L. (ancestor of Africanized bees) from that of A. m. mellifera, A. m. caucasia, A. m. ligustica, A. m. carnica, A. m. lamarcki, A. m. cypria, A. m. syriaca, and some A. m. iberiensis, but not from that of A. m. intermissa and some A. m. iberiensis. Nonetheless, given the very low frequency (<1%) of African non-A. m. scutellata mitotype present before arrival of Africanized bees in the United States, cytochrome b/BglII assay can be used to identify maternally Africanized bees with a high degree of reliability.

Reduced SNP Panels for Genetic Identification and Introgression Analysis in the Dark Honey Bee (Apis mellifera mellifera)

Beekeeping activities, especially queen trading, have shaped the distribution of honey bee (Apis mellifera) subspecies in Europe, and have resulted in extensive introductions of two eastern European C-lineage subspecies (A. m. ligustica and A. m. carnica) into the native range of the M-lineage A. m. mellifera subspecies in Western Europe. As a consequence, replacement and gene flow between native and commercial populations have occurred at varying levels across western European populations. Genetic identification and introgression analysis using molecular markers is an important tool for management and conservation of honey bee subspecies. Previous studies have monitored introgression by using microsatellite, PCR-RFLP markers and most recently, high density assays using single nucleotide polymorphism (SNP) markers. While the latter are almost prohibitively expensive, the information gained to date can be exploited to create a reduced panel containing the most ancestry-informative markers (AIMs) for those purposes with very little loss of information. The objective of this study was to design reduced panels of AIMs to verify the origin of A. m. mellifera individuals and to provide accurate estimates of the level of C-lineage introgression into their genome. The discriminant power of the SNPs using a variety of metrics and approaches including the Weir & Cockerham’s FST, an FST-based outlier test, Delta, informativeness (In), and PCA was evaluated. This study shows that reduced AIMs panels assign individuals to the correct origin and calculates the admixture level with a high degree of accuracy. These panels provide an essential tool in Europe for genetic stock identification and estimation of admixture levels which can assist management strategies and monitor honey bee conservation programs.

Revisiting the Iberian honey bee (Apis mellifera iberiensis) contact zone: maternal and genome-wide nuclear variations provide support for secondary contact from historical refugia.

Dissecting diversity patterns of organisms endemic to Iberia has been truly challenging for a variety of taxa, and the Iberian honey bee is no exception. Surveys of genetic variation in the Iberian honey bee are among the most extensive for any honey bee subspecies. From these, differential and complex patterns of diversity have emerged, which have yet to be fully resolved. Here, we used a genome‐wide data set of 309 neutrally tested single nucleotide polymorphisms (SNPs), scattered across the 16 honey bee chromosomes, which were genotyped in 711 haploid males. These SNPs were analysed along with an intergenic locus of the mtDNA, to reveal historical patterns of population structure across the entire range of the Iberian honey bee. Overall, patterns of population structure inferred from nuclear loci by multiple clustering approaches and geographic cline analysis were consistent with two major clusters forming a well‐defined cline that bisects Iberia along a northeastern–southwestern axis, a pattern that remarkably parallels that of the mtDNA. While a mechanism of primary intergradation or isolation by distance could explain the observed clinal variation, our results are more consistent with an alternative model of secondary contact between divergent populations previously isolated in glacial refugia, as proposed for a growing list of other Iberian taxa. Despite current intense honey bee management, human‐mediated processes have seemingly played a minor role in shaping Iberian honey bee genetic structure. This study highlights the complexity of the Iberian honey bee patterns and reinforces the importance of Iberia as a reservoir of Apis mellifera diversity.

Mitochondrial DNA variation of Apis mellifera iberiensis: further insights from a large-scale study using sequence data of the tRNAleu-cox2 intergenic region

A large-scale survey of the Iberian honey bee (Apis mellifera iberiensis) diversity patterns, using sequence data of the tRNAleu-cox2 mitochondrial DNA (mtDNA) region, demonstrates that earlier studies based on the DraI test missed significant components of genetic variation. Based on results from this survey, existing haplotype names were revised and updated following a nomenclature system established earlier and extended herein for the intergenic region. A more complete picture of the complex diversity patterns of IHBs is revealed that includes 164 novel haplotypes, 113 belonging to lineage A and 51 to lineage M and within lineage A and 69 novel haplotypes that belong to sub-lineage AI, 13 to AII, and 31 to AIII. Within lineage M, two novel haplotypes show a striking architecture with features of lineages A and M, which based on sequence comparisons and relationships among haplotypes are seemingly ancestral. These data expand our knowledge of the complex architecture of the tRNAleu-cox2 intergenic region in Apis mellifera and re-emphasizes the importance of Iberia as a source of honey bee mtDNA diversity.

Novel diagnostic tools for Asian (Apis cerana) and European (Apis mellifera) honey authentication

Honey can be produced by different species of honeybees, with two being of economic importance due to their use in apiculture, namely Apis mellifera (known as European honeybee) and Apis cerana (known as Asian honeybee). Due to the decline of the wild populations of the Asian honeybee, this honey generally attains much higher market value, being prone to adulteration. This work aims at proposing new tools, based on the use of molecular markers, for the entomological authentication of honey. To this end, new species-specific primers were designed targeting the tRNAleu-cox2 intergenic region and allowing the detection of A. cerana DNA by qualitative polymerase chain reaction (PCR). Additionally, a novel real-time PCR method with high resolution melting analysis was developed to target the 16S rRNA gene of both bee species, allowing their discrimination in different clusters. The proposed methodologies were further applied with success in the authentication of Asian and European honey samples by the identification of honeybee DNA, demonstrating the usefulness of these simple and cost-effective new approaches.