Serial founder effects slow range expansion in an invasive social insect

Invasive populations often experience founder effects: a loss of genetic diversity relative to the source population, due to a small number of founders. Even where these founder effects do not impact colonization success, theory predicts they might affect the rate at which invasive populations expand. This is because secondary founder effects are generated at advancing population edges, further reducing local genetic diversity and elevating genetic load. We show that in an expanding invasive population of the Asian honey bee (Apis cerana), genetic diversity is indeed lowest at range edges, including at the complementary sex determiner, csd, a locus that is homozygous-lethal. Consistent with lower local csd diversity, range edge colonies had lower brood viability than colonies in the range centre. Further, simulations of a newly-founded and expanding honey bee population corroborate the spatial patterns in mean colony fitness observed in our empirical data and show that such genetic load at range edges will slow the rate of population expansion.

Post-invasion selection acts on standing genetic variation despite a severe founding bottleneck

Invasive populations often have lower genetic diversity relative to the native-range populations from which they derive.1,2 Despite this, many biological invaders succeed in their new environments, in part due to rapid adaptation.3,4,5,6 Therefore, the role of genetic bottlenecks in constraining the adaptation of invaders is debated.7,8,9,10 Here, we use whole-genome resequencing of samples from a 10-year time-series dataset, representing the natural invasion of the Asian honey bee (Apis cerana) in Australia, to investigate natural selection occurring in the aftermath of a founding event. We find that Australia's A. cerana population was founded by as few as one colony, whose arrival was followed by a period of rapid population expansion associated with an increase of rare variants.11 The bottleneck resulted in a steep loss of overall genetic diversity, yet we nevertheless detected loci with signatures of positive selection during the first years post-invasion. When we investigated the origin of alleles under selection, we found that selection acted primarily on the variation introduced by founders and not on the variants that arose post-invasion by mutation. In all, our data highlight that selection on standing genetic variation can occur in the early years post-invasion, even where founding bottlenecks are severe.

Rapid evolution, rearrangements and whole mitogenome duplication in the Australian stingless bees Tetragonula (Hymenoptera: Apidae): A steppingstone towards understanding mitochondrial function and evolution

The extreme conservation of mitochondrial genomes in metazoans poses a significant challenge to understanding mitogenome evolution. However, the presence of variation in gene order or genome structure, found in a small number of taxa, can provide unique insights into this evolution. Previous work on two stingless bees in the genus Tetragonula (T. carbonaria and T. hockingsi) revealed highly divergent CO1 regions between them and when compared to the bees from the same tribe (Meliponini), indicating rapid evolution. Using mtDNA isolation and Illumina sequencing, we elucidated the mitogenomes of both species. In both species, there has been a duplication of the whole mitogenome to give a total genome size of 30,666 bp in T. carbonaria; and 30,662 bp in T. hockingsi. These duplicated genomes present a circular structure with two identical and mirrored copies of all 13 protein coding genes and 22 tRNAs, with the exception of a few tRNAs that are present as single copies. In addition, the mitogenomes are characterized by rearrangements of two block of genes. We believe that rapid evolution is present in the whole Indo-Malay/Australasian group of Meliponini but is extraordinarily elevated in T. carbonaria and T. hockingsi, probably due to founder effect, low effective population size and the mitogenome duplication. All these features - rapid evolution, rearrangements, and duplication - deviate significantly from the vast majority of the mitogenomes described so far, making the mitogenomes of Tetragonula unique opportunities to address fundamental questions of mitogenome function and evolution.

Abundant small RNAs in the reproductive tissues and eggs of the honey bee, Apis mellifera

Background

Polyandrous social insects such as the honey bee are prime candidates for parental manipulation of gene expression in offspring. Although there is good evidence for parent-of-origin effects in honey bees the epigenetic mechanisms that underlie these effects remain a mystery. Small RNA molecules such as miRNAs, piRNAs and siRNAs play important roles in transgenerational epigenetic inheritance and in the regulation of gene expression during development.

Results

Here we present the first characterisation of small RNAs present in honey bee reproductive tissues: ovaries, spermatheca, semen, fertilised and unfertilised eggs, and testes. We show that semen contains fewer piRNAs relative to eggs and ovaries, and that piRNAs and miRNAs which map antisense to genes involved in DNA regulation and developmental processes are differentially expressed between tissues. tRNA fragments are highly abundant in semen and have a similar profile to those seen in the semen of other animals. Intriguingly we also find abundant piRNAs that target the sex determination locus, suggesting that piRNAs may play a role in honey bee sex determination.

Conclusions

We conclude that small RNAs may play a fundamental role in honey bee gametogenesis and reproduction and provide a plausible mechanism for parent-of-origin effects on gene expression and reproductive physiology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08478-9.

Males Are Capable of Long-Distance Dispersal in a Social Bee

Pollinator conservation is aided by knowledge of dispersal behavior, which shapes gene flow and population structure. In many bees, dispersal is thought to be male-biased, and males’ movements may be critical to maintaining gene flow in disturbed and fragmented habitats. Yet male bee movements are challenging to track directly and male dispersal ability remains poorly understood in most species. Here, we combine field manipulations and models to assess male dispersal ability in a stingless bee (Tetragonula carbonaria). We placed colonies with virgin queens at varying distances apart (1–48 km), genotyped the males that gathered at mating aggregations outside each colony, and used pairwise sibship assignment to determine the distribution of likely brothers across aggregations. We then compared simulations of male dispersal to our observed distributions and found best-fit models when males dispersed an average of 2–3 km (>2-fold female flight ranges), and maximum of 20 km (30-fold female flight ranges). Our data supports the view that male bee dispersal can facilitate gene flow over long-distances, and thus play a key role in bee populations’ resilience to habitat loss and fragmentation. In addition, we show that the number of families contributing to male aggregations can be used to estimate local stingless bee colony densities, allowing population monitoring of these important tropical pollinators.

Reproductive plasticity and oogenesis in the queen honey bee (Apis mellifera)

In the honey bee (Apis mellifera), queen and worker castes originate from identical genetic templates but develop into different phenotypes. Queens lay up to 2000 eggs daily whereas workers are sterile in the queen's presence. Periodically queens stop laying: during swarming, when resources are scarce in winter, and when they are confined to a cage by beekeepers. We used confocal microscopy and gene expression assays to investigate the control of oogenesis in the ovaries of honey bee queens that were caged inside and outside the colony. We find evidence that queens use a different combination of 'checkpoints' to regulate oogenesis compared to honey bee workers and other insect species. However, both queen and worker castes likely use the same programmed cell death pathways to terminate oocyte development at their caste-specific checkpoints. Our results also suggest that a key factor driving the termination of oogenesis in queens is nutritional stress. Thus, queens may regulate oogenesis via the same regulatory pathways that were utilised by ancestral solitary species but likely have adjusted physiological checkpoints to suit their highly-derived life history.

DNA methylation is not a driver of gene expression reprogramming in young honey bee workers

The presence of DNA methylation marks within genic intervals, also called gene body methylation, is an evolutionarily-conserved epigenetic hallmark of animal and plant methylomes. In social insects, gene body methylation is thought to contribute to behavioural plasticity, for example between foragers and nurse workers, by modulating gene expression. However, recent studies have suggested that the majority of DNA methylation is sequence-specific, and therefore cannot act as a flexible mediator between environmental cues and gene expression. To address this paradox, we examined whole-genome methylation patterns in the brains and ovaries of young honey bee workers that had been subjected to divergent social contexts: the presence or absence of the queen. Although these social contexts are known to bring about extreme changes in behavioral and reproductive traits through differential gene expression, we found no significant differences between the methylomes of workers from queenright and queenless colonies. In contrast, thousands of regions were differentially methylated between colonies, and these differences were not associated with differential gene expression in the subset of genes examined. Methylation patterns were highly similar between brain and ovary tissues and only differed in nine regions. These results strongly indicate that DNA methylation is not a driver of differential gene expression between tissues or behavioral morphs. Finally, despite the lack of difference in methylation patterns, queen presence affected the expression of all four DNA methyltransferase genes, suggesting that these enzymes have roles beyond DNA methylation. Therefore, the functional role of DNA methylation in social insect genomes remains an open question.

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.

How does epigenetics influence the course of evolution?

Epigenetics is the study of changes in gene activity that can be transmitted through cell divisions but cannot be explained by changes in the DNA sequence. Epigenetic mechanisms are central to gene regulation, phenotypic plasticity, development and the preservation of genome integrity. Epigenetic mechanisms are often held to make a minor contribution to evolutionary change because epigenetic states are typically erased and reset at every generation, and are therefore, not heritable. Nonetheless, there is growing appreciation that epigenetic variation makes direct and indirect contributions to evolutionary processes. First, some epigenetic states are transmitted intergenerationally and affect the phenotype of offspring. Moreover, bona fide heritable 'epialleles' exist and are quite common in plants. Such epialleles could, therefore, be subject to natural selection in the same way as conventional DNA sequence-based alleles. Second, epigenetic variation enhances phenotypic plasticity and phenotypic variance and thus can modulate the effect of natural selection on sequence-based genetic variation. Third, given that phenotypic plasticity is central to the adaptability of organisms, epigenetic mechanisms that generate plasticity and acclimation are important to consider in evolutionary theory. Fourth, some genes are under selection to be 'imprinted' identifying the sex of the parent from which they were derived, leading to parent-of-origin-dependent gene expression and effects. These effects can generate hybrid disfunction and contribute to speciation. Finally, epigenetic processes, particularly DNA methylation, contribute directly to DNA sequence evolution, because they act as mutagens on the one hand and modulate genome stability on the other by keeping transposable elements in check. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

The role of epigenetics, particularly DNA methylation, in the evolution of caste in insect societies

Eusocial insects can be defined as those that live in colonies and have distinct queens and workers. For most species, queens and workers arise from a common genome, and so caste-specific developmental trajectories must arise from epigenetic processes. In this review, we examine the epigenetic mechanisms that may be involved in the regulation of caste dimorphism. Early work on honeybees suggested that DNA methylation plays a causal role in the divergent development of queen and worker castes. This view has now been challenged by studies that did not find consistent associations between methylation and caste in honeybees and other species. Evidence for the involvement of methylation in modulating behaviour of adult workers is also inconsistent. Thus, the functional significance of DNA methylation in social insects remains equivocal. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Vitellogenin expression in the ovaries of adult honeybee workers provides insights into the evolution of reproductive and social traits

Social insects are notable for having two female castes that exhibit extreme differences in their reproductive capacity. The molecular basis of these differences is largely unknown. Vitellogenin (Vg) is a powerful antioxidant and insulin-signalling regulator used in oocyte development. Here we investigate how Royal Jelly (the major food of honeybee queens) and queen mandibular pheromone (a major regulator of worker fertility), affect the longevity and reproductive status of honey bee workers, the expression of Vg, its receptor VgR and associated regulatory proteins. We find that Vg is expressed in the ovaries of workers and that workers fed a queen diet of Royal Jelly have increased Vg expression in the ovaries. Surprisingly, we find that expression of Vg is not associated with ovary activation in workers, suggesting that this gene has potentially acquired non-reproductive functions. Therefore, Vg expression in the ovaries of honeybee workers provides further support for the Ovarian Ground Plan Hypothesis, which argues that genes implicated in the regulation of reproduction have been co-opted to regulate behavioural differences between queens and workers.

Parent-of-origin effects, allele-specific expression, genomic imprinting and paternal manipulation in social insects

Haplo-diploidy and the relatedness asymmetries it generates mean that social insects are prime candidates for the evolution of genomic imprinting. In single-mating social insect species, some genes may be selected to evolve genomic mechanisms that enhance reproduction by workers when they are inherited from a female. This situation reverses in multiple mating species, where genes inherited from fathers can be under selection to enhance the reproductive success of daughters. Reciprocal crosses between subspecies of honeybees have shown strong parent-of-origin effects on worker reproductive phenotypes, and this could be evidence of such genomic imprinting affecting genes related to worker reproduction. It is also possible that social insect fathers directly affect gene expression in their daughters, for example, by placing small interfering RNA molecules in semen. Gene expression studies have repeatedly found evidence of parent-specific gene expression in social insects, but it is unclear at this time whether this arises from genomic imprinting, paternal manipulation, an artefact of cyto-nuclear interactions, or all of these. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Intergenerational transfer of DNA methylation marks in the honey bee

The evolutionary significance of epigenetic inheritance is controversial. While epigenetic marks such as DNA methylation can affect gene function and change in response to environmental conditions, their role as carriers of heritable information is often considered anecdotal. Indeed, near-complete DNA methylation reprogramming, as occurs during mammalian embryogenesis, is a major hindrance for the transmission of nongenetic information between generations. Yet it remains unclear how general DNA methylation reprogramming is across the tree of life. Here we investigate the existence of epigenetic inheritance in the honey bee. We studied whether fathers can transfer epigenetic information to their daughters through DNA methylation. We performed instrumental inseminations of queens, each with four different males, retaining half of each male's semen for whole genome bisulfite sequencing. We then compared the methylation profile of each father's somatic tissue and semen with the methylation profile of his daughters. We found that DNA methylation patterns were highly conserved between tissues and generations. There was a much greater similarity of methylomes within patrilines (i.e., father-daughter subfamilies) than between patrilines in each colony. Indeed, the samples' methylomes consistently clustered by patriline within colony. Samples from the same patriline had twice as many shared methylated sites and four times fewer differentially methylated regions compared to samples from different patrilines. Our findings indicate that there is no DNA methylation reprogramming in bees and, consequently, that DNA methylation marks are stably transferred between generations. This points to a greater evolutionary potential of the epigenome in invertebrates than there is in mammals.

Queen pheromone modulates the expression of epigenetic modifier genes in the brain of honeybee workers

Pheromones are used by many insects to mediate social interactions. In the highly eusocial honeybee (Apis mellifera), queen mandibular pheromone (QMP) is involved in the regulation of the reproductive and other behaviour of workers. The molecular mechanisms by which QMP acts are largely unknown. Here, we investigate how genes responsible for epigenetic modifications to DNA, RNA and histones respond to the presence of QMP in the environment. We show that several of these genes are upregulated in the honeybee brain when workers are exposed to artificial QMP. We propose that pheromonal communication systems, such as those used by social insects, evolved to respond to environmental signals by making use of existing epigenomic machineries.

Thermal adaptation in the honeybee (Apis mellifera) via changes to the structure of malate dehydrogenase

In honeybees there are three alleles of cytosolic malate dehydrogenase gene: F, M and S. Allele frequencies are correlated with environmental temperature, suggesting that the alleles have temperature-dependent fitness benefits. We determined the enzyme activity of each allele across a range of temperatures in vitro The F and S alleles have higher activity and are less sensitive to high temperatures than the M allele, which loses activity after incubation at temperatures found in the thorax of foraging bees in hot climates. Next, we predicted the protein structure of each allele and used molecular dynamics simulations to investigate their molecular flexibility. The M allozyme is more flexible than the S and F allozymes at 50°C, suggesting a plausible explanation for its loss of activity at high temperatures, and has the greatest structural flexibility at 15°C, suggesting that it can retain some enzyme activity at cooler temperatures. MM bees recovered from 2 h of cold narcosis significantly better than all other genotypes. Combined, these results explain clinal variation in malate dehydrogenase allele frequencies in the honeybee at the molecular level.

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.

Genetic Characterization of Exotic Commercial Honey Bee (Hymenoptera: Apidae) Populations in Thailand Reveals High Genetic Diversity and Low Population Substructure

Domestication of animal species is often associated with a reduction in genetic diversity. The honey bee, Apis mellifera Linnaeus, 1758, has been managed by beekeepers for millennia for both honey and wax production and for crop pollination. Here we use both microsatellite markers and sequence data from the mitochondrial COI gene to evaluate genetic variation of managed A. mellifera in Thailand, where the species is introduced. Microsatellite analysis revealed high average genetic diversity with expected heterozygosities ranging from 0.620 ± 0.184 to 0.734 ± 0.071 per locus per province. Observed heterozygosities were generally lower than those expected under Hardy-Weinberg equilibrium, both locally and across the population as a whole. Mitochondrial sequencing revealed that the frequency of two evolutionary linages (C-Eastern European and O-Middle Eastern) are similar to those observed in a previous survey 10 yr ago. Our results suggest that Thai beekeepers are managing their A. mellifera in ways that retain overall genetic diversity, but reduce genetic diversity between apiaries.

Conserved numts mask a highly divergent mitochondrial-COI gene in a species complex of Australian stingless bees Tetragonula (Hymenoptera: Apidae)

Tetragonula carbonaria, Tetragonula davenporti, Tetragonula hockingsi and Tetragonula mellipes comprise a species complex of Australian stingless bee species known as the 'Carbonaria' group. The species are difficult to distinguish morphologically and the major species-defining characters relate to comb architecture and nest entrance ornamentation. The taxonomy of the group is further complicated by likely nuclear mitochondrial pseudogenes (numts) and inter-specific hybrids. Here we demonstrate the existence of COI numts and isolate and characterize the 'true' mt-COI gene in T. carbonaria and T. hockingsi. Numts were isolated from enriched-nuclear DNA extraction followed by PCR amplification and Sanger sequencing, and were recognized by the presence of deletions and/or premature stop codons in the translated sequences. The mt-COI sequences were obtained from NGS sequencing using purified mtDNA. In T. carbonaria, two numts (numt1 and numt2) were identified and a third (numt3) was identified in T. hockingsi. Numt2 and numt3 are similar (1.2% sequence divergence), indicating a recent common origin. The genetic distance between the mt-COI of the two Tetragonula species was higher than might be expected for closely related species, 16.5%, corroborating previous studies in which T. carbonaria and T. hockingsi were regarded as separate species. The three numts are more similar to the COI of other stingless bee species, including Australian Austroplebia australis and South American Melipona bicolor (81.7-83.9%) than to the mt-COI of their own species (70-71.4%). This is because the mt-COI of T. carbonaria and T. hockingsi differ greatly from other Meliponinae. Our findings explain some formerly puzzling aspects of Carbonaria biogeography, and misinterpreted amplifications.

Unique DNA Methylation Profiles Are Associated with cis-Variation in Honey Bees

DNA methylation is an important epigenetic modification that mediates diverse processes such as cellular differentiation, phenotypic plasticity, and genomic imprinting. Mounting evidence suggests that local DNA sequence variation can be associated with particular DNA methylation states, indicating that the interplay between genetic and epigenetic factors may contribute synergistically to the phenotypic complexity of organisms. Social insects such as ants, bees, and wasps have extensive phenotypic plasticity manifested in their different castes, and this plasticity has been associated with variation in DNA methylation. Yet, the influence of genetic variation on DNA methylation state remains mostly unknown. Here we examine the importance of sequence-specific methylation at the genome-wide level, using whole-genome bisulfite sequencing of the semen of individual honey bee males. We find that individual males harbor unique DNA methylation patterns in their semen, and that genes that are more variable at the epigenetic level are also more likely to be variable at the genetic level. DNA sequence variation can affect DNA methylation by modifying CG sites directly, but can also be associated with local variation in cis that is not CG-site specific. We show that covariation in sequence polymorphism and DNA methylation state contributes to the individual-specificity of epigenetic marks in social insects, which likely promotes their retention across generations, and their capacity to influence evolutionary adaptation.

Conflict and major transitions - why we need true queens

In contrast to human societies, where kings and queens can be sources of conflict, we argue that the morphologically distinct queens of insect colonies are central to the minimization of conflict within their societies. Thus, the evolution of irreversible queen and worker castes represents a major transition in social evolution. Queens are selected to become better reproducers, and workers are selected to become better workers. The reproductive success of queens and workers are, therefore, inextricably linked. Workers achieve reproductive success by assisting the queen, whereas the queen needs her workers to provide her with the wherewithal to raise her brood. The tighter the mutual dependence, the lower conflict, and the larger insect societies can become. As the queen becomes a better breeder, workers are selected to become better at raising their siblings. Yet, nothing in nature is ever free of conflict and with the evolution of a true worker caste a new set of conflicts arises. Multiple mating by queens in particular opens the door to a new set of conflicts. Ironically, multiple mating can only evolve once within-colony conflict is reduced by evolving a true worker caste.

Social Parasitism in the Honeybee (Apis mellifera) Is Not Controlled by a Single SNP

The Cape bee (Apis mellifera capensis) is a subspecies of the honeybee, in which workers commonly lay diploid unfertilized eggs via a process known as thelytoky. A recent study aimed to map the genetic basis of this trait in the progeny of a single capensis queen where workers laid either diploid (thelytokous) or haploid (arrhenotokous) eggs. A nonsynonymous single nucleotide polymorphism (SNP) in a gene of unknown function was reported to be strongly associated with thelytoky in this colony. Here, we analyze genome sequences from a global sample of A. mellifera and identify populations where the proposed thelytoky allele at this SNP is common but thelytoky is absent. We also analyze genome sequences of three capensis queens produced by thelytoky and find that, contrary to predictions, they do not carry the proposed thelytoky allele. The proposed SNP is therefore neither sufficient nor required to produce thelytoky in A. mellifera.

Different bees, different needs: how nest-site requirements have shaped the decision-making processes in homeless honeybees (Apis spp.)

During reproductive swarming, a honeybee swarm needs to decide on a new nest site and then move to the chosen site collectively. Most studies of swarming and nest-site selection are based on one species, Apis mellifera Natural colonies of A. mellifera live in tree cavities. The quality of the cavity is critical to the survival of a swarm. Other honeybee species nest in the open, and have less strict nest-site requirements, such as the open-nesting dwarf honeybee Apis floreaApis florea builds a nest comprised of a single comb suspended from a twig. For a cavity-nesting species, there is only a limited number of potential nest sites that can be located by a swarm, because suitable sites are scarce. By contrast, for an open-nesting species, there is an abundance of equally suitable twigs. While the decision-making process of cavity-nesting bees is geared towards selecting the best site possible, open-nesting species need to coordinate collective movement towards areas with potential nest sites. Here, we argue that the nest-site selection processes of A. florea and A. mellifera have been shaped by each species' specific nest-site requirements. Both species use the same behavioural algorithm, tuned to allow each species to solve their species-specific problem.This article is part of the theme issue 'Collective movement ecology'.

Genetic Diversity in the Progeny of Commercial Australian Queen Honey Bees (Hymenoptera: Apidae) Produced in Autumn and Early Spring

Honey bee [Apis mellifera L. (Hymenoptera: Apidae)] queens are polyandrous, mating with an average 12 males (drones). Polyandry has been shown to confer benefits to queens and the colonies they head, including avoidance of inviable brood that can arise via sex locus homozygosity, increased resilience to pests and pathogens, and increased survival and productivity, leading to improved colony-level fitness. Queens with an effective mating frequency (ke) greater than 7 are considered adequately mated, whereas queens that fall below this threshold head colonies that have increased likelihood of failure and may be less productive for beekeepers. We determined ke in queens produced in early Spring and Autumn by five Australian commercial queen producers to determine whether the queens they produced were suitably mated. Drone populations are low at these times of year, and therefore, there is an increased risk that queens would fall below the ke > 7 threshold. We found that 33.8% of Autumn-produced queens did not meet the threshold, whereas 93.8% of Spring queens were adequately mated. The number of colonies contributing drones to the mating pool was similarly high in both seasons, suggesting that although many colonies have drones, their numbers may be decreased in Autumn and management strategies may be required to boost drone numbers at this time. Finally, queens had similar levels of homozygosity to workers, and inbreeding coefficients were very low, suggesting that inbreeding is not a problem.

Sex mosaics in the honeybee: how haplodiploidy makes possible the evolution of novel forms of reproduction in social Hymenoptera

Hymenoptera are haplodiploid: females arise from fertilized, diploid eggs, while males arise from unfertilized, haploid eggs. The cytogenetic mechanisms underlying haplodiploidy enable remarkable phenomena including female cloning, male cloning and gynandromorphy (sex mosaics). We collected 11 newly emerged putative gynandromorph honeybees from a single colony, assessed the sex of various tissues morphologically and determined the genetic origin (maternal or paternal) of each tissue by genotyping. Ten bees were gynandromorphs with one to three distinct paternal origins. Remarkably, one bee carried no maternal alleles. This bee had female organs throughout, and arose from the fusion of two sperm nuclei. This is the first reported case in the Hymenoptera of sperm fusion resulting in a female, emphasizing the flexibility for social insect reproduction and potentially novel colony-level social structures.

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.

Androgenesis: where males hijack eggs to clone themselves

Androgenesis is a form of quasi-sexual reproduction in which a male is the sole source of the nuclear genetic material in the embryo. Two types of androgenesis occur in nature. Under the first type, females produce eggs without a nucleus and the embryo develops from the male gamete following fertilization. Evolution of this type of androgenesis is poorly understood as the parent responsible for androgenesis (the mother) gains no benefit from it. Ultimate factors driving the evolution of the second type of androgenesis are better understood. In this case, a zygote is formed between a male and a female gamete, but the female genome is eliminated. When rare, androgenesis with genome elimination is favoured because an androgenesis-determining allele has twice the reproductive success of an allele that determines sexual reproduction. Paradoxically, except in hermaphrodites, a successful androgenetic strain can drive such a male-biased sex ratio that the population goes extinct. This likely explains why androgenesis with genome elimination appears to be rarer than androgenesis via non-nucleate eggs, although both forms are either very rare or remain largely undetected in nature. Nonetheless, some highly invasive species including ants and freshwater clams are androgenetic, for reasons that are largely unexplained.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

Subfamily-dependent alternative reproductive strategies in worker honeybees

Functional worker sterility is the defining feature of insect societies. Yet, workers are sometimes found reproducing in their own or foreign colonies. The proximate mechanisms underlying these alternative reproductive phenotypes are keys to understanding how reproductive altruism and selfishness are balanced in eusocial insects. In this study, we show that in honeybee (Apis mellifera) colonies, the social environment of a worker, that is, the presence and relatedness of the queens in a worker's natal colony and in surrounding colonies, significantly influences her fertility and drifting behaviour. Furthermore, subfamilies vary in the frequency of worker ovarian activation, propensity to drift and the kind of host colony that is targeted for reproductive parasitism. Our results show that there is an interplay between a worker's subfamily, reproductive state and social environment that substantially affects her reproductive phenotype. Our study further indicates that honeybee populations show substantial genetic variance for worker reproductive strategies, suggesting that no one strategy is optimal under all the circumstances that a typical worker may encounter.

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.

Extreme polyandry aids the establishment of invasive populations of a social insect

Although monandry is believed to have facilitated the evolution of eusociality, many highly eusocial insects have since evolved extreme polyandry. The transition to extreme polyandry was likely driven by the benefits of within-colony genetic variance to task specialization and/or disease resistance, but the extent to which it confers secondary benefits, once evolved, is unclear. Here we investigate the consequences of extreme polyandry on the invasive potential of the Asian honey bee, Apis cerana. In honey bees and other Hymenoptera, small newly founded invasive populations must overcome the genetic constraint of their sex determination system that requires heterozygosity at a sex-determining locus to produce viable females. We find A. cerana queens in an invasive population mate with an average of 27 males (range 16-42) that would result in the founding queen/s carrying 75% of their source population's sex alleles in stored sperm. This mating frequency is similar to native-range Chinese A. cerana (mean 29 males, range 19-46). Simulations reveal that extreme polyandry reduces the risk, relative to monandry or moderate polyandry, that colonies produce a high incidence of inviable brood in populations that have experienced a founder event, that is, when sex allele diversity is low and/or allele frequencies are unequal. Thus, extreme polyandry aids the invasiveness of A. cerana in two ways: (1) by increasing the sex locus allelic richness carried to new populations with each founder, thereby increasing sex locus heterozygosity; and (2) by reducing the population variance in colony fitness following a founder event.

THE EVOLUTION OF WORKER STERILITY IN HONEY BEES: AN INVESTIGATION INTO A BEHAVIORAL MUTANT CAUSING FAILURE OF WORKER POLICING

Normally, worker honey bees (Apis mellifera) only lay eggs when their colony is queenless. When a queen is present, worker egg-laying is controlled by mutual "policing" behavior in which any rare worker-laid eggs are eaten by other workers. However, an extremely rare behavioral phenotype arises in which workers develop functional ovaries and lay large numbers of eggs despite the presence of the queen. In this study, microsatellite analysis was used to determine the maternity of drones produced in such a colony under various conditions. One subfamily was found to account for about 90% of drone progeny, with the remainder being laid by other subfamilies or the queen. No evidence of queen policing was found. After a one-month period of extreme worker oviposition in spring, the colony studied reverted to normal behavior and showed no signs of worker oviposition. However, upon removal of the queen, workers commenced oviposition very quickly. Significantly, the subfamily that laid eggs when the queen was present did not contribute to the drone production when the colony was queenless. However, another subfamily contributed a disproportionately large number of drones. The frequency of worker oviposition appears to be determined by opposing selective forces. Individual bees benefit from personal reproduction, whereas other bees and the colony are disadvantaged by it. Thus a behavioral polymorphism can be maintained in the population in which some workers can escape worker policing, with balancing selection at the colony level to detect and eliminate these mutations.

Queen pheromone regulates programmed cell death in the honey bee worker ovary

In social insect colonies the presence of a queen, secreting her pheromones, is a key environmental cue for regulating the reproductive state of workers. However, until recently the proximate molecular mechanisms underlying facultative worker sterility were unidentified. Studies into worker oogenesis in the honey bee (Apis mellifera) have indicated that programmed cell death is central to the regulation of oogenesis. Here we investigate how queen pheromone, age of the worker and ovary state affect both programmed cell death and cell number in worker ovaries. We describe a novel method to simultaneously measure programmed cell death (caspase activity) and live cell number (estimated from the amount of adenosine triphosphate) in an insect tissue. Workers exposed to queen pheromone have higher levels of caspase activity in the ovary than those not exposed. Our results suggest that queen pheromone triggers programmed cell death at the mid-oogenesis checkpoint causing the abortion of worker oocytes and reproductive inhibition of the worker caste. Nonetheless, high caspase activity is present in activated ovaries from workers not exposed to queen pheromone. This caspase activity is most likely to be from the nurse cells undergoing programmed cell death, in late oogenesis, for normal oocyte development. Our study shows that the social environment of an organism can influence programmed cell death within a tissue.

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.

Entomology: A Bee Farming a Fungus

Farming is done not only by humans, but also by some ant, beetle and termite species. With the discovery of a stingless bee farming a fungus that provides benefits to its larvae, bees can be added to this list.

Associations between reproduction and work in workers of the Asian hive bee Apis cerana

If a honey bee (Apis spp.) colony becomes queenless, about 1/3 of young workers activate their ovaries and produce haploid male-producing eggs. In doing so queenless workers maximize their inclusive fitness because the normal option of vicarious production of relatives via their queen's eggs is no longer available. But if many workers are engaged in reproduction, how does a queenless colony continue to feed its brood and forage? Here we show that in the Asian hive bee Apis cerana hypopharyngeal gland (HPG) size is larger in queenless workers than in queenright workers and that bees undertaking brood-rearing tasks have larger HPG than same-aged bees that are foraging. In queenless colonies, workers with a smaller number of ovarioles are more likely to have activated ovaries. This reinforces the puzzling observation that a large number of ovarioles reduces reproductive success in queenless A. cerana. It further suggests that reproductive workers either avoid foraging or transition to foraging later in life than non-reproductive workers. Finally, our study also showed that ovary activation and larger-than-average numbers of ovarioles had no statistically detectable influence on foraging specialization for pollen or nectar.

Regulation of oogenesis in honey bee workers via programed cell death

Reproductive division of labour characterises eusociality. Currently little is known about the mechanisms that underlie the 'sterility' of the worker caste, but queen pheromone plays a major role in regulating the reproductive state. Here we investigate oogenesis in the young adult honey bee worker ovary in the presence of queen pheromone and in its absence. When queen pheromone is absent, workers can activate their ovaries and have well-developed follicles. When queen pheromone is present, even though workers have non-activated ovaries, they continually produce oocytes which are aborted at an early stage. Therefore, irrespective of the presence of the queen, the young adult worker ovary contains oocytes. By this means young workers retain reproductive plasticity. The degeneration of the germ cells in the ovarioles of workers in the presence of queen pheromone has the morphological hallmarks of programmed cell death. Therefore the mechanistic basis of 'worker sterility' relies in part on the regulation of oogenesis via programmed cell death. Our results suggest that honey bees have co-opted a highly conserved checkpoint at mid-oogenesis to regulate the fertility of the worker caste.

RNA-sequencing elucidates the regulation of behavioural transitions associated with the mating process in honey bee queens

Background

Mating is a complex process, which is frequently associated with behavioural and physiological changes. However, understanding of the genetic underpinnings of these changes is limited. Honey bees are both a model system in behavioural genomics, and the dominant managed pollinator of human crops; consequently understanding the mating process has both pure and applied value. We used next-generation transcriptomics to probe changes in gene expression in the brains of honey bee queens, as they transition from virgin to mated reproductive status. In addition, we used CO₂-narcosis, which induces oviposition without mating, to isolate the process of reproductive maturation.

Results

The mating process produced significant changes in the expression of vision, chemo-reception, metabolic, and immune-related genes. Differential expression of these genes maps clearly onto known behavioural and physiological changes that occur during the transition from being a virgin queen to a newly-mated queen. A subset of these changes in gene expression were also detected in CO₂-treated queens, as predicted from previous physiological studies. In addition, we compared our results to previous studies that used microarray techniques across a range of experimental time-points. Changes in expression of immune- and vision-related genes were common to all studies, supporting an involvement of these groups of genes in the mating process.

Conclusions

Our study is an important step in understanding the molecular mechanisms regulating post-mating behavioural transitions in a natural system. The weak overlap in patterns of gene expression with previous studies demonstrates the high sensitivity of genome-wide approaches. Thus, while we build on previous microarray studies that explored post-mating changes in honey bees, the broader experimental design, use of RNA-sequencing, and focus on Australian honey bees, which remain free from the devastating parasite Varroa destructor, in the current study, provide unique insights into the biology of the mating process in honey bees.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1750-7) contains supplementary material, which is available to authorized users.

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.

Pheromones affecting ovary activation and ovariole loss in the Asian honey bee Apis cerana

The Asian hive bee Apis cerana has similar queen mandibular pheromones (QMP) to the Western honey bee Apismellifera. However the effects of individual QMP components have never been tested to determine their effects on the reproductive physiology of A. cerana workers. We fed one queen equivalent of each of the major components of A. cerana QMP to groups of c.a. 500 day-old, caged, workers twice a day until the workers were 10 days old. Half of the cages were also provided with 10% royal jelly in the food. Workers were sampled each day and dissected to determine the number of ovarioles and the degree of ovary activation (egg development). In cages treated with 9-carbon fatty acids ovary activation was minimal, whereas the 10-carbon acids suppressed ovary activation very little. Royal jelly enhanced ovary activation, especially in cages treated with 10-carbon acids. The number of ovarioles declined with bee age, but the rate of decline was slowed by the 9-carbon acids in particular. The results show conservation of the composition and function of QMP between A. cerana and A. mellifera and support the hypothesis that QMP is an honest signal of queen fecundity rather than a chemical castrator of workers.

Inheritance of thelytoky in the honey bee Apis mellifera capensis

Asexual reproduction via thelytokous parthenogenesis is widespread in the Hymenoptera, but its genetic underpinnings have been described only twice. In the wasp Lysiphlebus fabarum and the Cape honey bee Apis mellifera capensis the origin of thelytoky have each been traced to a single recessive locus. In the Cape honey bee it has been argued that thelytoky (th) controls the thelytoky phenotype and that a deletion of 9 bp in the flanking intron downstream of exon 5 (tae) of the gemini gene switches parthenogenesis from arrhenotoky to thelytoky. To further explore the mode of inheritance of thelytoky, we generated reciprocal backcrosses between thelytokous A. m. capensis and the arrhenotokous A. m. scutellata. Ten genetic markers were used to identify 108 thelytokously produced offspring and 225 arrhenotokously produced offspring from 14 colonies. Patterns of appearance of thelytokous parthenogenesis were inconsistent with a single locus, either th or tae, controlling thelytoky. We further show that the 9 bp deletion is present in the arrhenotokous A. m. scutellata population in South Africa, in A. m. intermissa in Morocco and in Africanized bees from Brazil and Texas, USA, where thelytoky has not been reported. Thus the 9  p deletion cannot be the cause of thelytoky. Further, we found two novel tae alleles. One contains the previously described 9 bp deletion and an additional deletion of 7 bp nearby. The second carries a single base insertion with respect to the wild type. Our data are consistent with the putative th locus increasing reproductive capacity.

Queen signals in a stingless bee: suppression of worker ovary activation and spatial distribution of active compounds

In most species of social insect the queen signals her presence to her workers via pheromones. Worker responses to queen pheromones include retinue formation around the queen, inhibition of queen cell production and suppression of worker ovary activation. Here we show that the queen signal of the Brazilian stingless bee Friesella schrottkyi is a mixture of cuticular hydrocarbons. Stingless bees are therefore similar to ants, wasps and bumble bees, but differ from honey bees in which the queen's signal mostly comprises volatile compounds originating from the mandibular glands. This shows that cuticular hydrocarbons have independently evolved as the queen's signal across multiple taxa, and that the honey bees are exceptional. We also report the distribution of four active queen-signal compounds by Matrix-assisted laser desorption/ionization (MALDI) imaging. The results indicate a relationship between the behavior of workers towards the queen and the likely site of secretion of the queen's pheromones.

Phantom alternatives influence food preferences in the eastern honeybee Apis cerana

Most models of animal choice behaviour assume that desirable but unavailable options, such as a high quality, but inhabited nest sites, do not influence an individual's preferences for the remaining options. However, experiments suggest that in mammals, the mere presence of such 'phantom' alternatives can alter, and even reverse, an individual's preferences for other items in a choice set. Phantom alternatives may be widespread in nature, as they occur whenever a resource is visible, but unavailable at the time of choice. They are particularly relevant for nectar-foraging animals, where previously rewarding flowers may sometimes be empty. Here, we investigate the effect of phantom alternatives on feeder preferences in the eastern honeybee, Apis cerana. First, we tested the effects of unattractive and attractive phantom alternatives by presenting individual bees with either a binary choice set containing two feeders that differed strongly in two qualities, but were equally preferred overall ('option 1' and 'option 2'), or a ternary choice set containing option 1, option 2 and one of two phantom types (unattractive and attractive). Secondly, we determined whether phantoms increase (similarity effect) or decrease (dissimilarity effect) preference for phantom-similar choices. In binary trials, bees had no significant preference for option 1 or option 2. However, after encountering an attractive phantom alternative, individual bees preferred option 2. The unattractive phantom did not influence bee preferences. Phantoms consistently changed individual bee preferences in favour of the phantom-similar choice. This means that the presence of an attractive food source, even if it is unavailable, can influence preference relationships between remaining items in the choice set. Our findings highlight the importance of considering the potential for phantom effects when studying the foraging behaviour of animals. Our results are particularly relevant for nectarivores, where empty but previously rewarding flowers are a common occurrence. Since an increase in pollinator visits can result in higher seed set, our results open up the possibility that by shifting pollinator preferences, empty flowers could have otherwise-unpredicted influences on community composition, plant-pollinator interactions and pollinator behaviour.

The dynamic DNA methylation cycle from egg to sperm in the honey bee Apis mellifera

In honey bees (Apis mellifera), the epigenetic mark of DNA methylation is central to the developmental regulation of caste differentiation, but may also be involved in additional biological functions. In this study, we examine the whole genome methylation profiles of three stages of the haploid honey bee genome: unfertilised eggs, the adult drones that develop from these eggs and the sperm produced by these drones. These methylomes reveal distinct patterns of methylation. Eggs and sperm show 381 genes with significantly different CpG methylation patterns, with the vast majority being more methylated in eggs. Adult drones show greatly reduced levels of methylation across the genome when compared with both gamete samples. This suggests a dynamic cycle of methylation loss and gain through the development of the drone and during spermatogenesis. Although fluxes in methylation during embryogenesis may account for some of the differentially methylated sites, the distinct methylation patterns at some genes suggest parent-specific epigenetic marking in the gametes. Extensive germ line methylation of some genes possibly explains the lower-than-expected frequency of CpG sites in these genes. We discuss the potential developmental and evolutionary implications of methylation in eggs and sperm in this eusocial insect species.