Twenty-five-year study of Nosema spp. in honey bees (Apis mellifera) in Serbia

A total of 7386 samples of adult honey bees from different areas of Serbia (fifteen regions and 79 municipalities) were selected for light microscopy analysis for Nosema species during 1992–2017. A selection of honey bee samples from colonies positive for microsporidian spores during 2009–2011, 2015 and 2017 were then subjected to molecular diagnosis by multiplex PCR using specific primers for a region of the 16S rRNA gene of Nosema species. The prevalence of microsporidian spore-positive bee colonies ranged between 14.4% in 2013 and 65.4% in 1992. PCR results show that Nosema ceranae is not the only Nosema species to infect honey bees in Serbia. Mixed N. apis/N. ceranae infections were detected in the two honey bee samples examined by mPCR during 2017. The beekeeping management of disease prevention, such as replacement of combs and queens and hygienic handling of colonies are useful in the prevention of Nosema infection.     

High Concentrations of the Alarm Pheromone Component,Isopentyl Acetate,Reduces Foraging and Dancing in <Emphasis Type="Italic">Apis mellifera</Emphasis> Ligustica and <Emphasis Type="Italic">Apis cerana</Emphasis> Cerana

A honeybee colony is a superorganism that has evolved precise communication systems, which allow the colony to gather information from numerous individuals and coordinate its behavior. Alarm pheromones, such as isopentyl acetate (IPA), the main component of sting alarm pheromone, play a critical role in the coordination of individual behaviors as well as colony communication in honeybee colonies. In this study, honeybees (Apis mellifera ligustica and Apis cerana cerana) were exposed to relatively high levels of IPA at a foraging site (6–8 bee equivalents) and inside their colony (28–58 bee equivalents) to investigate the influence of alarm pheromones on foraging activity and hive flight activity. IPA reduced the number of bees that flew out the hive, foraged, and waggle danced. Under both contexts in the hive and at the food source, IPA can therefore inhibit honey bee foraging and foraging communication.     

Sudden deaths and colony population decline in Greek honey bee colonies

During June and July of 2009, sudden deaths, tremulous movements and population declines of adult honey bees were reported by the beekeepers in the region of Peloponnesus (Mt. Mainalo), Greece. A preliminary study was carried out to investigate these unexplained phenomena in this region. In total, 37 bee samples, two brood frames containing honey bee brood of various ages, eight sugar samples and four sugar patties were collected from the affected colonies. The samples were tested for a range of pests, pathogens and pesticides. Symptomatic adult honey bees tested positive for Varroa destructor,Nosema ceranae, Chronic bee paralysis virus (CBPV), Acute paralysis virus (ABPV), Deformed wing virus (DWV), Sacbrood virus (SBV) and Black queen cell virus (BQCV), but negative for Acarapis woodi. American Foulbrood was absent from the brood samples. Chemical analysis revealed that amitraz, thiametoxan, clothianidin and acetamiprid were all absent from symptomatic adult bees, sugar and sugar patty samples. However, some bee samples, were contaminated with imidacloprid in concentrations between 14 ng/g and 39 ng/g tissue. We present: the infection of Greek honey bees by multiple viruses; the presence of N. ceranae in Greek honey bees and the first record of imidacloprid (neonicotonoid) residues in Greek honey bee tissues. The presence of multiple pathogens and pesticides made it difficult to associate a single specific cause to the depopulation phenomena observed in Greece, although we believe that viruses and N. ceranae synergistically played the most important role. A follow up in-depth survey across all Greek regions is required to provide context to these preliminary findings.     

Pleiotropy, epistasis and new QTL: the genetic architecture of honey bee foraging behavior

The regulation of division of labor in social insects, particularly in the honey bee (Apis mellifera L.), has received considerable attention from a number of biological subdisciplines, including quantitative and behavioral genetics, because of the high complexity of the behavioral traits involved. The foraging choices of honey bee workers can be accurately quantified, and previous studies have made the foraging behavior of honey bees one of the best studied naturally occurring behavioral phenotypes. Three quantitative trait loci (QTL) have been identified that influence a set of foraging variables, including the concentration of nectar collected and the amount of pollen and nectar brought back to the hive. This study extends previous genetic investigations and represents the most comprehensive investigation of the genetic architecture of these foraging variables. We examined the effects of markers for the three established QTL and for one further candidate gene (Amfor), in two reciprocal backcross populations. These populations were also used to carry out two new QTL mapping studies, with over 400 Amplified Fragment Length Polymorphism (AFLP) markers in each. We detected a variety of effects of the genetic markers for the established QTL and the candidate gene, which were mostly epistatic in nature. A few new QTL could be detected with a variety of mapping techniques. Our results add complexity to the genetic architecture of the foraging behavior of the honey bee. Specifically, we support the hypotheses that pln1, pln2, pln3, and Amfor are involved in the regulation of foraging behavior in the honey bee and add some new factors that deserve further study in the future.     

Mitochondrial DNA Diversity of Honey Bees (Apis mellifera) from Unmanaged Colonies and Swarms in the United States

To study the genetic diversity of honey bees (Apis mellifera L.) from unmanaged colonies in the United States, we sequenced a portion of the mitochondrial DNA COI–COII region. From the 530 to 1,230 bp amplicon, we observed 23 haplotypes from 247 samples collected from 12 states, representing three of the four A. mellifera lineages known to have been imported into the United States (C, M, and O). Six of the 13 C lineage haplotypes were not found in previous queen breeder studies in the United States. The O lineage accounted for 9% of unmanaged colonies which have not yet been reported in queen breeder studies. The M lineage accounted for a larger portion of unmanaged samples (7%) than queen breeder samples (3%). Based on our mitochondrial DNA data, the genetic diversity of unmanaged honey bees in the United States differs significantly from that of queen breeder populations (p < 0.00001). The detection of genetically distinct maternal lineages of unmanaged honey bees suggests that these haplotypes may have existed outside the managed honey bee population for a long period.     

Widespread Occurrence of Chemical Residues in Beehive Matrices from Apiaries Located in Different Landscapes of Western France

Background

The honey bee, Apis mellifera, is frequently used as a sentinel to monitor environmental pollution. In parallel, general weakening and unprecedented colony losses have been reported in Europe and the USA, and many factors are suspected to play a central role in these problems, including infection by pathogens, nutritional stress and pesticide poisoning. Honey bee, honey and pollen samples collected from eighteen apiaries of western France from four different landscape contexts during four different periods in 2008 and in 2009 were analyzed to evaluate the presence of pesticides and veterinary drug residues.

Methodology/Findings

A multi-residue analysis of 80 compounds was performed using a modified QuEChERS method, followed by GC-ToF and LC−MS/MS. The analysis revealed that 95.7%, 72.3% and 58.6% of the honey, honey bee and pollen samples, respectively, were contaminated by at least one compound. The frequency of detection was higher in the honey samples (n = 28) than in the pollen (n = 23) or honey bee (n = 20) samples, but the highest concentrations were found in pollen. Although most compounds were rarely found, some of the contaminants reached high concentrations that might lead to adverse effects on bee health. The three most frequent residues were the widely used fungicide carbendazim and two acaricides, amitraz and coumaphos, that are used by beekeepers to control Varroa destructor. Apiaries in rural-cultivated landscapes were more contaminated than those in other landscape contexts, but the differences were not significant. The contamination of the different matrices was shown to be higher in early spring than in all other periods.

Conclusions/Significance

Honey bees, honeys and pollens are appropriate sentinels for monitoring pesticide and veterinary drug environmental pollution. This study revealed the widespread occurrence of multiple residues in beehive matrices and suggests a potential issue with the effects of these residues alone or in combination on honey bee health.     

Three QTL in the honey bee Apis mellifera L. suppress reproduction of the parasitic mite Varroa destructor

Varroa destructor is a highly virulent ectoparasitic mite of the honey bee Apis mellifera and a major cause of colony losses for global apiculture. Typically, chemical treatment is essential to control the parasite population in the honey bee colony. Nevertheless a few honey bee populations survive mite infestation without any treatment. We used one such Varroa mite tolerant honey bee lineage from the island of Gotland, Sweden, to identify quantitative trait loci (QTL) controlling reduced mite reproduction. We crossed a queen from this tolerant population with drones from susceptible colonies to rear hybrid queens. Two hybrid queens were used to produce a mapping population of haploid drones. We discriminated drone pupae with and without mite reproduction, and screened the genome for potential QTL using a total of 216 heterozygous microsatellite markers in a bulk segregant analysis. Subsequently, we fine mapped three candidate target regions on chromosomes 4, 7, and 9. Although the individual effect of these three QTL was found to be relatively small, the set of all three had significant impact on suppression of V. destructor reproduction by epistasis. Although it is in principle possible to use these loci for marker-assisted selection, the strong epistatic effects between the three loci complicate selective breeding programs with the Gotland Varroa tolerant honey bee stock.     

后基因组时代的蜜蜂QTL研究

继果蝇、按蚊和家蚕之后,意大利蜜蜂Apis mellifera（膜翅目： 蜜蜂科）成为又一种被完整测得基因组序列的昆虫。从此,蜜蜂研究进入后基因组时代。作为一种典型的社会性昆虫,许多和蜜蜂社会生活紧密相关的性状都是数量性状。这些性状研究中广泛涉及到了数量性状位点（quantitative traits loci,QTL）定位研究。本文综述了应用QTL对蜜蜂取食行为、自卫行为、体长、逆转学习等的研究现状,同时结合国内外最新研究进展,总结并展望了后基因组时代蜜蜂QTL的研究方向。     

A revision of subspecies structure of western honey bee Apis mellifera

The taxonomy of honey bee A. mellifera contains a lot of issues due to the specificity of population structure, features of biology and resolutions of honey bee subspecies discrimination methods. There are a lot of transition zones between ranges of subspecies which led to the gradual changes of characteristics among neighbor subspecies. The modern taxonomic pattern of honey bee Apis mellifera is given in this paper. Thirty-three distinct honey bee subspecies are distributed across all Africa (11 subspecies), Western Asia and the Middle East (9 subspecies), and Europe (13 subspecies). All honey bee subspecies are subdivided into 5 evolutionary lineages: lineage A (10 subspecies) and its sublineage Z (3 subspecies), lineage M (3 subspecies), lineage C (10 subspecies), lineage O (3 subspecies), lineage Y (1 subspecies), lineage C or O (3 subspecies).     

A Cell Line Resource Derived from Honey Bee (Apis mellifera) Embryonic Tissues

A major hindrance to the study of honey bee pathogens or the effects of pesticides and nutritional deficiencies is the lack of controlled in vitro culture systems comprised of honey bee cells. Such systems are important to determine the impact of these stress factors on the developmental and cell biology of honey bees. We have developed a method incorporating established insect cell culture techniques that supports sustained growth of honey bee cells in vitro. We used honey bee eggs mid to late in their embryogenesis to establish primary cultures, as these eggs contain cells that are progressively dividing. Primary cultures were initiated in modified Leibovitz’s L15 medium and incubated at 32^°C. Serial transfer of material from several primary cultures was maintained and has led to the isolation of young cell lines. A cell line (AmE-711) has been established that is composed mainly of fibroblast-type cells that form an adherent monolayer. Most cells in the line are diploid (2n = 32) and have the Apis mellifera karyotype as revealed by Giemsa stain. The partial sequence for the mitochondrial-encoded cytochrome c oxidase subunit I (Cox 1) gene in the cell line is identical to those from honey bee tissues and a consensus sequence for A. mellifera. The population doubling time is approximately 4 days. Importantly, the cell line is continuously subcultured every 10–14 days when split at a 1:3 ratio and is cryopreserved in liquid nitrogen. The cell culture system we have developed has potential application for studies aimed at honey bee development, genetics, pathogenesis, transgenesis, and toxicology.     

A SNP Based High-Density Linkage Map of Apis cerana Reveals a High Recombination Rate Similar to Apis mellifera

Background

The Eastern honey bee, Apis cerana Fabricius, is distributed in southern and eastern Asia, from India and China to Korea and Japan and southeast to the Moluccas. This species is also widely kept for honey production besides Apis mellifera. Apis cerana is also a model organism for studying social behavior, caste determination, mating biology, sexual selection, and host-parasite interactions. Few resources are available for molecular research in this species, and a linkage map was never constructed. A linkage map is a prerequisite for quantitative trait loci mapping and for analyzing genome structure. We used the Chinese honey bee, Apis cerana cerana to construct the first linkage map in the Eastern honey bee.

Results

F2 workers (N = 103) were genotyped for 126,990 single nucleotide polymorphisms (SNPs). After filtering low quality and those not passing the Mendel test, we obtained 3,000 SNPs, 1,535 of these were informative and used to construct a linkage map. The preliminary map contains 19 linkage groups, we then mapped the 19 linkage groups to 16 chromosomes by comparing the markers to the genome of A. mellfiera. The final map contains 16 linkage groups with a total of 1,535 markers. The total genetic distance is 3,942.7 centimorgans (cM) with the largest linkage group (180 loci) measuring 574.5 cM. Average marker interval for all markers across the 16 linkage groups is 2.6 cM.

Conclusion

We constructed a high density linkage map for A. c. cerana with 1,535 markers. Because the map is based on SNP markers, it will enable easier and faster genotyping assays than randomly amplified polymorphic DNA or microsatellite based maps used in A. mellifera.     

Flight and fight: a comparative view of the neurophysiology and genetics of honey bee defensive behavior

Honey bee nest defense involves guard bees that specialize in olfaction-based nestmate recognition and alarm-pheromone-mediated recruitment of nestmates to sting. Stinging is influenced by visual, tactile and olfactory stimuli. Both quantitative trait locus (QTL) mapping and behavioral studies point to guarding behavior as a key factor in colony stinging response. Results of reciprocal F1 crosses show that paternally inherited genes have a greater influence on colony stinging response than maternally inherited genes. The most active alarm pheromone component, isoamyl acetate (IAA) causes increased respiration and may induce stress analgesia in bees. IAA primes worker bees for 'fight or flight', possibly through actions of neuropeptides and/or biogenic amines. Studies of aggression in other species lead to an expectation that octopamine or 5-HT might play a role in honey bee defensive response. Genome sequence and QTL mapping identified 128 candidate genes for three regions known to influence defensive behavior. Comparative bioinformatics suggest possible roles of genes involved in neurogenesis and central nervous system (CNS) activity, and genes involved in sensory tuning through G-protein coupled receptors (GPCRs), such as an arrestin (AmArr4) and the metabotropic GABA(B) receptor (GABA-B-R1).     

Evaluation of the Distribution and Impacts of Parasites,Pathogens, and Pesticides on Honey Bee (Apis mellifera) Populations in East Africa

In East Africa, honey bees (Apis mellifera) provide critical pollination services and income for small-holder farmers and rural families. While honey bee populations in North America and Europe are in decline, little is known about the status of honey bee populations in Africa. We initiated a nationwide survey encompassing 24 locations across Kenya in 2010 to evaluate the numbers and sizes of honey bee colonies, assess the presence of parasites (Varroa mites and Nosema microsporidia) and viruses, identify and quantify pesticide contaminants in hives, and assay for levels of hygienic behavior. Varroa mites were present throughout Kenya, except in the remote north. Levels of Varroa were positively correlated with elevation, suggesting that environmental factors may play a role in honey bee host-parasite interactions. Levels of Varroa were negatively correlated with levels of hygienic behavior: however, while Varroa infestation dramatically reduces honey bee colony survival in the US and Europe, in Kenya Varroa presence alone does not appear to impact colony size. Nosema apis was found at three sites along the coast and one interior site. Only a small number of pesticides at low concentrations were found. Of the seven common US/European honey bee viruses, only three were identified but, like Varroa, were absent from northern Kenya. The number of viruses present was positively correlated with Varroa levels, but was not correlated with colony size or hygienic behavior. Our results suggest that Varroa, the three viruses, and Nosema have been relatively recently introduced into Kenya, but these factors do not yet appear to be impacting Kenyan bee populations. Thus chemical control for Varroa and Nosema are not necessary for Kenyan bees at this time. This study provides baseline data for future analyses of the possible mechanisms underlying resistance to and the long-term impacts of these factors on African bee populations.     

Simulating a base population in honey bee for molecular genetic studies

Background

Over the past years, reports have indicated that honey bee populations are declining and that infestation by an ecto-parasitic mite (Varroa destructor) is one of the main causes. Selective breeding of resistant bees can help to prevent losses due to the parasite, but it requires that a robust breeding program and genetic evaluation are implemented. Genomic selection has emerged as an important tool in animal breeding programs and simulation studies have shown that it yields more accurate breeding value estimates, higher genetic gain and low rates of inbreeding. Since genomic selection relies on marker data, simulations conducted on a genomic dataset are a pre-requisite before selection can be implemented. Although genomic datasets have been simulated in other species undergoing genetic evaluation, simulation of a genomic dataset specific to the honey bee is required since this species has a distinct genetic and reproductive biology. Our software program was aimed at constructing a base population by simulating a random mating honey bee population. A forward-time population simulation approach was applied since it allows modeling of genetic characteristics and reproductive behavior specific to the honey bee.

Results

Our software program yielded a genomic dataset for a base population in linkage disequilibrium. In addition, information was obtained on (1) the position of markers on each chromosome, (2) allele frequency, (3) χ² statistics for Hardy-Weinberg equilibrium, (4) a sorted list of markers with a minor allele frequency less than or equal to the input value, (5) average r² values of linkage disequilibrium between all simulated marker loci pair for all generations and (6) average r² value of linkage disequilibrium in the last generation for selected markers with the highest minor allele frequency.

Conclusion

We developed a software program that takes into account the genetic and reproductive biology specific to the honey bee and that can be used to constitute a genomic dataset compatible with the simulation studies necessary to optimize breeding programs. The source code together with an instruction file is freely accessible at http://msproteomics.org/Research/Misc/honeybeepopulationsimulator.html     

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.     

The influence of nest site selection on the population dynamics of Africanized honey bees in an urban landscape

Urban landscapes provide habitat for many species, including domesticated and feral honey bees, Apis mellifera L. (Hymenoptera: Apidae). With recent losses of managed honey bee colonies, there is increasing interest in feral honey bee colonies and their potential contribution to pollination services in agricultural, natural, and urban settings. However, in some regions the feral honey bee population consists primarily of Africanized honey bees. Africanized honey bees (AHB) are hybrids between European honey bees and the African honey bee, Apis mellifera scutellataLepeletier, and have generated economic, ecological, and human health concerns because of their aggressive behavior. In this study, we used two long‐term datasets (7–10 years) detailing the spatial and temporal distribution of AHB colonies in Tucson, AZ, USA, where feral colonies occupy a variety of cavities including water meter boxes. A stage‐structured matrix model was used to elucidate the implications of nest site selection and the effects of colony terminations on the structure and dynamics of the AHB population. Our results suggest that Tucson's AHB population is driven by a relatively small number of ‘source’ colonies that escape termination (ca. 0.165 colonies per km² or 125 colonies in total), although immigrating swarms and absconding colonies from the surrounding area may have also contributed to the stability of the Tucson AHB population. Furthermore, the structure of the population has likely been impacted by the number and spatial distribution of water meter boxes across the city. The study provides an example of how urban wildlife populations are driven by interactions among landscape structure, human management, and behavioral traits conferred by an invasive genotype.     

Development and Evaluation of PCR Assays for the Detection of Paenibacillus larvae in Honey Samples: Comparison with Isolation and Biochemical Characterization

PCR assays were developed for the direct detection of Paenibacillus larvae in honey samples and compared with isolation and biochemical characterization procedures. Different primer pairs, designed from the 16S rRNA and the metalloproteinase precursor gene regions, and different DNA extraction methods were tested and compared. The sensitivity of the reactions was evaluated by serial dilutions of DNA extracts obtained from P. larvae cultures. The specificity of the primers was assessed by analyzing related Paenibacillus and Bacillus strains isolated from honey. The PCR assays also amplified these related bacteria, but at lower sensitivity. In the next step, the PCR assays were applied to contaminated honey and other bee products originating from 15 countries. Lysozyme treatment followed by proteinase K digestion was determined to be the best DNA extraction method for P. larvae spores. The most sensitive primer pair detected P. larvae in 18 of 23 contaminated honey samples, as well as in pollen, wax, and brood. Honey specimens containing saprophyte bacilli and paenibacilli, but not P. larvae, were PCR negative. Although the isolation and biochemical identification method (BioLog) showed higher sensitivity and specificity, PCR proved to be a valuable technique for large-scale screening of honey samples for American foulbrood, especially considering its rapidity and moderate costs.     

Population structure of North African honey bees is influenced by both biological and anthropogenic factors

Honey bee diversity is under threat due to anthropogenic factors as the use of pesticides and the replacement of local colonies to recover from colony losses. To assess the effect of these activities on the genetic diversity and structure of North African honey bee colonies, we studied colonies from the north (Tellian) and the south (Saharan) regions in Algeria, by determining their mitochondrial haplotype and the variation of ten microsatellite loci. Particular haplotypes have been found with a high frequency in each region that may constitute subspecies-specific markers for Apis mellifera intermissa (haplotype A8 at the north) and A. m. sahariensis (haplotype A9 at the south). Moreover, the presence of the haplotype A8 in some Saharan colonies may reflect recent introductions of A. m. intermissa. Structure analysis suggests that a natural differentiation between honey bee populations from Saharan and Tellian regions still exists despite increased colony movements (migration, queen purchases, etc.) during the last decades. One apiary established for the conservation of A. m. sahariensis showed no indication of maternal introgression since all the colonies bear the same haplotype A9. Furthermore, Hardy–Weinberg equilibrium observed in this population indicates that this apiary is appropriate for conservation programs of A. m. sahariensis.     

Iridovirus and microsporidian linked to honey bee colony decline

Background

In 2010 Colony Collapse Disorder (CCD), again devastated honey bee colonies in the USA, indicating that the problem is neither diminishing nor has it been resolved. Many CCD investigations, using sensitive genome-based methods, have found small RNA bee viruses and the microsporidia, Nosema apis and N. ceranae in healthy and collapsing colonies alike with no single pathogen firmly linked to honey bee losses.

Methodology/Principal Findings

We used Mass spectrometry-based proteomics (MSP) to identify and quantify thousands of proteins from healthy and collapsing bee colonies. MSP revealed two unreported RNA viruses in North American honey bees, Varroa destructor-1 virus and Kakugo virus, and identified an invertebrate iridescent virus (IIV) (Iridoviridae) associated with CCD colonies. Prevalence of IIV significantly discriminated among strong, failing, and collapsed colonies. In addition, bees in failing colonies contained not only IIV, but also Nosema. Co-occurrence of these microbes consistently marked CCD in (1) bees from commercial apiaries sampled across the U.S. in 2006–2007, (2) bees sequentially sampled as the disorder progressed in an observation hive colony in 2008, and (3) bees from a recurrence of CCD in Florida in 2009. The pathogen pairing was not observed in samples from colonies with no history of CCD, namely bees from Australia and a large, non-migratory beekeeping business in Montana. Laboratory cage trials with a strain of IIV type 6 and Nosema ceranae confirmed that co-infection with these two pathogens was more lethal to bees than either pathogen alone.

Conclusions/Significance

These findings implicate co-infection by IIV and Nosema with honey bee colony decline, giving credence to older research pointing to IIV, interacting with Nosema and mites, as probable cause of bee losses in the USA, Europe, and Asia. We next need to characterize the IIV and Nosema that we detected and develop management practices to reduce honey bee losses.     

Development of single nucleotide polymorphism markers specific to Apis mellifera (Hymenoptera: Apidae) line displaying high hygienic behavior against Varroa destructor,an ectoparasitic mite

To control Varroa destructor, an ectoparasitic mite, a honey bee line possessing high hygienic behavior (HHB) against this mite has been bred in South Korea. However, a method that can diagnose and assess the HHB line from control (the low hygienic behavior, LHB) line has not been reported yet. Thus, the objective of this study was to develop single nucleotide polymorphism (SNP) markers through whole-genome sequencing of worker bees from HHB line of A. mellifera caucasica and LHB line of A. m. carnica (Hymenoptera: Apidae). A total of 319,445,977 sequence reads were mapped to the known A. mellifera reference genome (average coverage of 87.46%). In 2,316,128 and 3,266,756 SNPs from HHB and LHB line, respectively, 20 SNPs that showed homozygosity in each line were selected and eight SNPs were used to diagnose the HHB line either by typical PCR-restriction fragment length polymorphism or allele-specific PCR. Six of remaining SNPs were of different sizes, enabling relatively easy differentiation of these two honey bee lines on typical agarose gel. Another remaining six SNPs had different sequences, including SNP sites. These SNP markers can be used to diagnose and assess V. destructor-specific HHB line of honey bees.