Thrice out of Asia and the adaptive radiation of the western honey bee.

The origin of the western honey bee Apis mellifera has been intensely debated. Addressing this knowledge gap is essential for understanding the evolution and genetics of one of the world's most important pollinators. By analyzing 251 genomes from 18 native subspecies, we found support for an Asian origin of honey bees with at least three expansions leading to African and European lineages. The adaptive radiation of honey bees involved selection on a few genomic "hotspots." We found 145 genes with independent signatures of selection across all bee lineages, and these genes were highly associated with worker traits. Our results indicate that a core set of genes associated with worker and colony traits facilitated the adaptive radiation of honey bees across their vast distribution.

Adaptive, caste-specific changes to recombination rates in a thelytokous honeybee population.

The ability to clone oneself has clear benefits-no need for mate hunting or dilution of one's genome in offspring. It is therefore unsurprising that some populations of haplo-diploid social insects have evolved thelytokous parthenogenesis-the virgin birth of a female. But thelytokous parthenogenesis has a downside: the loss of heterozygosity (LoH) as a consequence of genetic recombination. LoH in haplo-diploid insects can be highly deleterious because female sex determination often relies on heterozygosity at sex-determining loci. The two female castes of the Cape honeybee, Apis mellifera capensis, differ in their mode of reproduction. While workers always reproduce thelytokously, queens always mate and reproduce sexually. For workers, it is important to reduce the frequency of recombination so as to not produce offspring that are homozygous. Here, we ask whether recombination rates differ between Cape workers and Cape queens that we experimentally manipulated to reproduce thelytokously. We tested our hypothesis that Cape workers have evolved mechanisms that restrain genetic recombination, whereas queens have no need for such mechanisms because they reproduce sexually. Using a combination of microsatellite genotyping and whole-genome sequencing we find that a reduction in recombination is confined to workers only.

Deformed wing virus prevalence and load in honeybees in South Africa.

Deformed wing virus (DWV) is an emerging honeybee pathogen that has appeared across the globe in the past 40 years. When transmitted by the parasitic varroa mite, it has been associated with the collapse of millions of colonies throughout the Northern Hemisphere. However, despite the presence of the mite in the Southern Hemisphere, infested colonies survive. This study investigated the prevalence of DWV genotypes A, B and C along with their viral loads in South Africa and compared the findings with recent data from Brazil, the UK and the USA. We found that DWV-B was the most prevalent genotype throughout South Africa, although the total DWV viral load was significantly lower (2.8E+07) than found in the Northern Hemisphere (2.8E+07 vs. 2.7E+10, p > 0.00001) and not significantly different to that found in Brazil (5E+06, p = 0.13). The differences in viral load can be explained by the mite resistance in Brazil and South Africa, since mite-infested cells containing high viral loads are removed by the bees, thus lowering the colony's viral burden. This behaviour is much less developed in the vast majority of honeybees in the Northern Hemisphere.

A Single Gene Causes Thelytokous Parthenogenesis, the Defining Feature of the Cape Honeybee Apis mellifera capensis.

In honeybees, the ability of workers to produce daughters asexually, i.e., thelytokous parthenogenesis, is restricted to a single subspecies inhabiting the Cape region of South Africa, Apis mellifera capensis. Thelytoky has unleashed new selective pressures and the evolution of traits such as social parasitism, invasiveness, and social cancer. Thelytoky arises from an abnormal meiosis that results in the fusion of two maternal pronuclei, restoring diploidy in newly laid eggs. The genetic basis underlying thelytoky is disputed. To resolve this controversy, we generated a backcross between thelytokous A. m. capensis and non-thelytokous A. m. scutellata from the neighboring population and looked for evidence of genetic markers that co-segregated with thelytokous reproduction in 49 backcross females. We found that markers associated with the gene GB45239 on chromosome 11, including non-synonymous variants, showed consistent co-segregation with thelytoky, whereas no other region did so. Alleles associated with thelytoky were present in all A. m. capensis genomes examined but were absent from all other honeybees worldwide including A. m. scutellata. GB45239 is derived in A. m. capensis and has a putative role in chromosome segregation. It is expressed in ovaries and is downregulated in thelytokous bees, likely because of polymorphisms in the promoter region. Our study reveals how mutations affecting the sequence and/or expression of a single gene can change the reproductive mode of a population.

Paternally-biased gene expression follows kin-selected predictions in female honey bee embryos.

The Kinship Theory of Genomic Imprinting (KTGI) posits that, in species where females mate with multiple males, there is selection for a male to enhance the reproductive success of his offspring at the expense of other males and his mating partner. Reciprocal crosses between honey bee subspecies show parent-of-origin effects for reproductive traits, suggesting that males modify the expression of genes related to female function in their female offspring. This effect is likely to be greater in the Cape honey bee (Apis mellifera capensis), because a male's daughters have the unique ability to produce female offspring that can develop into reproductive workers or the next queen without mating. We generated reciprocal crosses between Capensis and another subspecies and used RNA-seq to identify transcripts that are over- or underexpressed in the embryos, depending on the parental origin of the gene. As predicted, 21 genes showed expression bias towards the Capensis father's allele in colonies with a Capensis father, with no such bias in the reciprocal cross. A further six genes showed a consistent bias towards expression of the father's allele across all eight colonies examined, regardless of the direction of the cross. Consistent with predictions of the KTGI, six of the 21 genes are associated with female reproduction. No gene consistently showed overexpression of the maternal allele.

Strikingly high levels of heterozygosity despite 20 years of inbreeding in a clonal honey bee.

Inbreeding (the mating between closely related individuals) often has detrimental effects that are associated with loss of heterozygosity at overdominant loci, and the expression of deleterious recessive alleles. However, determining which loci are detrimental when homozygous, and the extent of their phenotypic effects, remains poorly understood. Here, we utilize a unique inbred population of clonal (thelytokous) honey bees, Apis mellifera capensis, to determine which loci reduce individual fitness when homozygous. This asexual population arose from a single worker ancestor approximately 20 years ago and has persisted for at least 100 generations. Thelytokous parthenogenesis results in a 1/3 of loss of heterozygosity with each generation. Yet, this population retains heterozygosity throughout its genome due to selection against homozygotes. Deep sequencing of one bee from each of the three known sub-lineages of the population revealed that 3,766 of 10,884 genes (34%) have retained heterozygosity across all sub-lineages, suggesting that these genes have heterozygote advantage. The maintenance of heterozygosity in the same genes and genomic regions in all three sub-lineages suggests that nearly every chromosome carries genes that show sufficient heterozygote advantage to be selectively detrimental when homozygous.

Viable Triploid Honey Bees (Apis mellifera capensis) Are Reliably Produced in the Progeny of CO2 Narcotised Queens.

The haplodiploid system of sex determination of Hymenoptera acts as an exaptation for species to evolve novel forms of asexual reproduction including thelytoky (clonal offspring of the mother). During normal reproduction in Hymenoptera, three of the four products of meiosis that are present in newly-laid eggs are lost as polar bodies, while the remaining pronucleus either develops as a haploid male or fuses with a sperm nucleus to produce a diploid zygote. In contrast, in thelytokous reproduction, which is uncommon but taxonomically widespread, two of the four products of meiosis fuse, as if one acted as a sperm. Queenless workers of Apis mellifera capensis, a subspecies of honey bee from South Africa, routinely reproduce thelytokously. Unmated A. m. capensis queens can also be induced to lay thelytokously by narcosis with carbon dioxide, but mated queens are never thelytokous. We artificially inseminated A. m. capensis queens using CO2 narcosis. Up to 1/3 of offspring workers carried two maternal alleles and an allele of one father whereas no three-allele progeny were seen in control queens of the arrhenotokous (unfertilized eggs result in males) subspecies A. m. scutellata Flow cytometry of three-allele individuals revealed that they were triploid and arose from the fertilization of a thelytokous fusion nucleus. We then reared six queens from a narcotized A. m. capensis queen and determined the ploidy of the offspring queens based on microsatellites. One of the five daughters was triploid. Following artificial insemination, this queen produced unfertilized thelytokous diploid eggs at high frequency, and unfertilized triploid eggs at much lower frequency. If fertilized, thelytokous diploid eggs were non-viable, even though triploidy in itself does not impede normal development. In contrast, when the rarer triploid eggs were fertilized, a proportion developed into viable tetraploids. Our study highlights the extraordinary developmental flexibility of haplo-diploid systems.

The dynamic association between ovariole loss and sterility in adult honeybee workers.

In the social insects, ovary state (the presence or absence of mature oocytes) and ovary size (the number of ovarioles) are often used as proxies for the reproductive capacity of an individual worker. Ovary size is assumed to be fixed post-eclosion whereas ovary state is demonstrably plastic post-eclosion. Here, we show that in fact ovary size declines as honeybee workers age. This finding is robust across two honeybee species: Apis mellifera and A. cerana The ovariole loss is likely to be due to the regression of particular ovarioles via programmed cell death. We also provide further support for the observation that honeybee workers with activated ovaries (mature oocytes present) most commonly have five ovarioles rather than a greater or smaller number. This result suggests that workers with more than five ovarioles are unable to physiologically support more than five activated ovarioles and that workers with fewer than five ovarioles are below a threshold necessary for ovary activation. As a worker's ovariole number declines with age, studies on worker ovariole number need to take this plasticity into account.

Parent-of-origin effects on genome-wide DNA methylation in the Cape honey bee (Apis mellifera capensis) may be confounded by allele-specific methylation.

BACKGROUND: Intersexual genomic conflict sometimes leads to unequal expression of paternal and maternal alleles in offspring, resulting in parent-of-origin effects. In honey bees reciprocal crosses can show strong parent-of-origin effects, supporting theoretical predictions that genomic imprinting occurs in this species. Mechanisms behind imprinting in honey bees are unclear but differential DNA methylation in eggs and sperm suggests that DNA methylation could be involved. Nonetheless, because DNA methylation is multifunctional, it is difficult to separate imprinting from other roles of methylation. Here we use a novel approach to investigate parent-of-origin DNA methylation in honey bees. In the subspecies Apis mellifera capensis, reproduction of females occurs either sexually by fertilization of eggs with sperm, or via thelytokous parthenogenesis, producing female embryos derived from two maternal genomes. RESULTS: We compared genome-wide methylation patterns of sexually-produced, diploid embryos laid by a queen, with parthenogenetically-produced diploid embryos laid by her daughters. Thelytokous embryos inheriting two maternal genomes had fewer hypermethylated genes compared to fertilized embryos, supporting the prediction that fertilized embryos have increased methylation due to inheritance of a paternal genome. However, bisulfite PCR and sequencing of a differentially methylated gene, Stan (GB18207) showed strong allele-specific methylation that was maintained in both fertilized and thelytokous embryos. For this gene, methylation was associated with haplotype, not parent of origin. CONCLUSIONS: The results of our study are consistent with predictions from the kin theory of genomic imprinting. However, our demonstration of allele-specific methylation based on sequence shows that genome-wide differential methylation studies can potentially confound imprinting and allele-specific methylation. It further suggests that methylation patterns are heritable or that specific sequence motifs are targets for methylation in some genes.

A SNP test to identify Africanized honeybees via proportion of 'African' ancestry.

The honeybee, Apis mellifera, is the world's most important pollinator and is ubiquitous in most agricultural ecosystems. Four major evolutionary lineages and at least 24 subspecies are recognized. Commercial populations are mainly derived from subspecies originating in Europe (75-95%). The Africanized honeybee is a New World hybrid of A. m. scutellata from Africa and European subspecies, with the African component making up 50-90% of the genome. Africanized honeybees are considered undesirable for bee-keeping in most countries, due to their extreme defensiveness and poor honey production. The international trade in honeybees is restricted, due in part to bans on the importation of queens (and semen) from countries where Africanized honeybees are extant. Some desirable strains from the United States of America that have been bred for traits such as resistance to the mite Varroa destructor are unfortunately excluded from export to countries such as Australia due to the presence of Africanized honeybees in the USA. This study shows that a panel of 95 single nucleotide polymorphisms, chosen to differentiate between the African, Eastern European and Western European lineages, can detect Africanized honeybees with a high degree of confidence via ancestry assignment. Our panel therefore offers a valuable tool to mitigate the risks of spreading Africanized honeybees across the globe and may enable the resumption of queen and bee semen imports from the Americas.

Selection on overdominant genes maintains heterozygosity along multiple chromosomes in a clonal lineage of honey bee.

Correlations between fitness and genome-wide heterozygosity (heterozygosity-fitness correlations, HFCs) have been reported across a wide range of taxa. The genetic basis of these correlations is controversial: do they arise from genome-wide inbreeding ("general effects") or the "local effects" of overdominant loci acting in linkage disequilibrium with neutral loci? In an asexual thelytokous lineage of the Cape honey bee (Apis mellifera capensis), the effects of inbreeding have been homogenized across the population, making this an ideal system in which to detect overdominant loci, and to make inferences about the importance of overdominance on HFCs in general. Here we investigate the pattern of zygosity along two chromosomes in 42 workers from the clonal Cape honey bee population. On chromosome III (which contains the sex-locus, a gene that is homozygous-lethal) and chromosome IV we show that the pattern of zygosity is characterized by loss of heterozygosity in short regions followed by the telomeric restoration of heterozygosity. We infer that at least four selectively overdominant genes maintain heterozygosity on chromosome III and three on chromosome IV via local effects acting on neutral markers in linkage disequilibrium. We conclude that heterozygote advantage and local effects may be more common and evolutionarily significant than is generally appreciated.

A parent-of-origin effect on honeybee worker ovary size.

Apis mellifera capensis is unique among honeybees in that unmated workers can produce pseudo-clonal female offspring via thelytokous parthenogenesis. Workers use this ability to compete among themselves and with their queen to be the mother of new queens. Males could therefore enhance their reproductive success by imprinting genes that enhance fertility in their daughter workers. This possibility sets the scene for intragenomic conflict between queens and drones over worker reproductive traits. Here, we show a strong parent-of-origin effect for ovary size (number of ovarioles) in reciprocal crosses between two honeybee subspecies, A. m. capensis and Apis mellifera scutellata. In this cross, workers with an A. m. capensis father had 30% more ovarioles than genotypically matched workers with an A. m. scutellata father. Other traits we measured (worker weight at emergence and the presence/absence of a spermatheca) are influenced more by rearing conditions than by parent-of-origin effects. Our study is the first to show a strong epigenetic (or, less likely, cytoplasmic maternal) effect for a reproductive trait in the honeybee and suggests that a search for parent-of-origin effects in other social insects may be fruitful.

Cheaters sometimes prosper: targeted worker reproduction in honeybee (Apis mellifera) colonies during swarming.

Kin selection theory predicts that honeybee (Apis mellifera) workers should largely refrain from producing their own offspring, as the workers collectively have higher inclusive fitness if they rear the sons of their mother, the queen. Studies that have quantified levels of ovary activation and reproduction among workers have largely supported this prediction. We sampled pre-emergent male pupae and adult workers from seven colonies at regular intervals throughout the reproductive part of the season. We show that the overall contribution of workers to male (drone) production is 4.2%, nearly 40 times higher than is generally reported, and is highest during reproductive swarming, when an average of 6.2% of the males genotyped are worker-produced. Similarly, workers in our samples were 100 times more likely to have active ovaries than previously assumed. Worker reproduction is seasonally influenced and peaks when colonies are rearing new queens. Not all worker subfamilies contribute equally to reproduction. Instead, certain subfamilies are massively over-represented in drone brood. By laying eggs within the period in which many colonies produce virgin queens, these rare worker subfamilies increase their direct fitness via their well-timed sons.

Maintenance and loss of heterozygosity in a thelytokous lineage of honey bees (Apis mellifera capensis).

An asexual lineage that reproduces by automictic thelytokous parthenogenesis has a problem: rapid loss of heterozygosity resulting in effective inbreeding. Thus, the circumstances under which rare asexual lineages thrive provide insights into the trade-offs that shape the evolution of alternative reproductive strategies across taxa. A socially parasitic lineage of the Cape honey bee, Apis mellifera capensis, provides an example of a thelytokous lineage that has endured for over two decades. It has been proposed that cytological adaptations slow the loss of heterozygosity in this lineage. However, we show that heterozygosity at the complementary sex determining (csd) locus is maintained via selection against homozygous diploid males that arise from recombination. Further, because zygosity is correlated across the genome, it appears that selection against diploid males reduces loss of homozygosity at other loci. Selection against homozygotes at csd results in substantial genetic load, so that if a thelytokous lineage is to endure, unusual ecological circumstances must exist in which asexuality permits such a high degree of fecundity that the genetic load can be tolerated. Without these ecological circumstances, sex will triumph over asexuality. In A. m. capensis, these conditions are provided by the parasitic interaction with its conspecific host, Apis mellifera scutellata.

Asexually produced Cape honeybee queens (Apis mellifera capensis) reproduce sexually.

Unmated workers of the Cape honeybee Apis mellifera capensis can produce female offspring including daughter queens. As worker-laid queens are produced asexually, we wondered whether these asexually produced individuals reproduce asexually or sexually. We sampled 11 colonies headed by queens known to be the clonal offspring of workers and genotyped 23 worker offspring from each queen at 5 microsatellite loci. Without exception, asexually produced queens produced female worker offspring sexually. In addition, we report the replacement of a queen by her asexually produced granddaughter, with this asexually produced queen also producing offspring sexually. Hence, once a female larva is raised as a queen, mating and sexual reproduction appears to be obligatory in this subspecies, despite the fact that worker-laid queens are derived from asexual lineages.

A thelytokous lineage of socially parasitic honey bees has retained heterozygosity despite at least 10 years of inbreeding.

The honey bee population of South Africa is divided into two subspecies: a northern population in which queenless workers reproduce arrhenotokously and a southern one in which workers reproduce thelytokously. A hybrid zone separates the two, but on at least three occasions the northern population has become infested by reproductive workers derived from the southern population. These parasitic workers lay in host colonies parthenogenetically, resulting in yet more parasites. The current infestation is 20-year old--surprising because an asexual lineage is expected to show a decline in vigor over time due to increasing homozygosity. The decline is expected to be acute in honey bees, where homozygosity at the sex locus is lethal. We surveyed colonies from the zone of infestation and genotyped putative parasites at two sets of linked microsatellite loci. We confirm that there is a single clonal lineage of parasites that shows minor variations arising from recombination events. The lineage shows high levels of heterozygosity, which may be maintained by selection against homozygotes, or by a reduction in recombination frequency within the lineage. We suggest that the clonal lineage can endure the costs of asexual reproduction because of the fitness benefits of its parasitic life history.

Maternity of emergency queens in the Cape honey bee, Apis mellifera capensis.

During reproductive swarming, some workers of the Cape honey bee, Apis mellifera capensis, lay eggs in queen cells, many of which are reared to maturity. However, it is unknown if workers are able to lay in queen cells immediately after queen loss during an episode of emergency queen rearing. In this study we experimentally de-queened colonies and determined the maternity of larvae and pupae that were reared as queens. This allowed us to determine how soon after queen loss workers contribute to the production of new queens. We were further interested to see if workers would preferentially raise new queens from queen-laid brood if this was introduced later. We performed our manipulations in two different settings: an apiary setting where colonies were situated close together and a more natural situation in which the colonies were well separated. This allowed us to determine how the vicinity of other colonies affects the presence of parasites. We found that workers do indeed contribute to queen cell production immediately after the loss of their queen, thus demonstrating that some workers either have activated ovaries even when their colony has a queen or are able to activate their ovaries extremely rapidly. Queen-laid brood introduced days after queen loss was ignored, showing that workers do not prefer to raise new queens from queen brood when given a choice. We also detected non-natal parasitism of queen cells in both settings. We therefore conclude that some A. m. capensis genotypes specialize in parasitizing queen cells.

A quantitative study of worker reproduction in queenright colonies of the Cape honey bee, Apis mellifera capensis.

Reproduction by workers is rare in honey bee colonies that have an active queen. By not producing their own offspring and preventing other workers from producing theirs, workers are thought to increase their inclusive fitness due to their higher average relatedness towards queen-produced male offspring compared with worker-produced male offspring. But there is one exception. Workers of the Cape honey bee, Apis mellifera capensis, are able to produce diploid female offspring via thelytokous parthenogenesis and thus produce clones of themselves. As a result, worker reproduction and tolerance towards worker-produced offspring is expected to be more permissive than in arrhenotokous (sub)species where worker offspring are male. Here we quantify the extent to which A. m. capensis workers contribute to reproduction in queenright colonies using microsatellite analyses of pre-emergent brood. We show that workers produced 10.5% of workers and 0.48% of drones. Most of the workers' contribution towards the production of new workers coincided with the colonies producing new queens during reproductive swarming.

Thelytokous parthenogenesis in unmated queen honeybees (Apis mellifera capensis): central fusion and high recombination rates.

The subspecies of honeybee indigenous to the Cape region of South Africa, Apis mellifera capensis, is unique because a high proportion of unmated workers can lay eggs that develop into females via thelytokous parthenogenesis involving central fusion of meiotic products. This ability allows pseudoclonal lineages of workers to establish, which are presently widespread as reproductive parasites within the honeybee populations of South Africa. Successful long-term propagation of a parthenogen requires the maintenance of heterozygosity at the sex locus, which in honeybees must be heterozygous for the expression of female traits. Thus, in successful lineages of parasitic workers, recombination events are reduced by an order of magnitude relative to meiosis in queens of other honeybee subspecies. Here we show that in unmated A. m. capensis queens treated to induce oviposition, no such reduction in recombination occurs, indicating that thelytoky and reduced recombination are not controlled by the same gene. Our virgin queens were able to lay both arrhenotokous male-producing haploid eggs and thelytokous female-producing diploid eggs at the same time, with evidence that they have some voluntary control over which kind of egg was laid. If so, they are able to influence the kind of second-division meiosis that occurs in their eggs post partum.

Inheritance of traits associated with reproductive potential in Apis mellifera capensis and Apis mellifera scutellata workers.

When workers of the thelytokous Cape honeybee, Apis mellifera capensis, come into contact with colonies of the neighboring arrhenotokous subspecies Apis mellifera scutellata, they can become lethal social parasites. We examined the inheritance of 3 traits (number of ovarioles, number of basitarsal hairs, and size of spermatheca) that are thought to be associated with reproductive potential in A. m. capensis workers. To do so, we produced hybrid A. m. scutellata/A. m. capensis queens and backcrossed them to either A. m. capensis or A. m. scutellata drones. We then measured the 3 traits in parental, hybrid, and backcross offspring. We show that the 3 traits are phenotypically correlated. We also show that the expression of ovariole number, basitarsal hairs, and size of spermatheca is influenced by the genotype of the individual and the rearing environment but that the influence of the rearing environment is less important to the number of ovarioles. We hypothesize a single recessive allele (l), present at high frequency in natural A. m. capensis populations, which when homozygous causes larvae to elicit more food. This increased feeding as larvae causes resulting adult workers to develop more queen-like morphology and increased reproductive potential. The number of ovarioles, in contrast, appears to be under independent genetic control.