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.

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.

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.

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.

Global allele polymorphism indicates a high rate of allele genesis at a locus under balancing selection.

When selection favours rare alleles over common ones (balancing selection in the form of negative frequency-dependent selection), a locus may maintain a large number of alleles, each at similar frequency. To better understand how allelic richness is generated and maintained at such loci, we assessed 201 sequences of the complementary sex determiner (csd) of the Asian honeybee (Apis cerana), sampled from across its range. Honeybees are haplodiploid; hemizygotes at csd develop as males and heterozygotes as females, while homozygosity is lethal. Thus, csd is under strong negative frequency-dependent selection because rare alleles are less likely to end up in the lethal homozygous form. We find that in A. cerana, as in other Apis, just a few amino acid differences between csd alleles in the hypervariable region are sufficient to trigger female development. We then show that while allelic lineages are spread across geographical regions, allelic differentiation is high between populations, with most csd alleles (86.3%) detected in only one sample location. Furthermore, nucleotide diversity in the hypervariable region indicates an excess of recently arisen alleles, possibly associated with population expansion across Asia since the last glacial maximum. Only the newly invasive populations of the Austral-Pacific share most of their csd alleles. In all, the geographic patterns of csd diversity in A. cerana indicate that high mutation rates and balancing selection act together to produce high rates of allele genesis and turnover at the honeybee sex locus, which in turn leads to its exceptionally high local and global polymorphism.

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.

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.

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.

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.

Workers' sons rescue genetic diversity at the sex locus in an invasive honey bee population.

The hallmark of eusociality is the division of labour between reproductive (queen) and nonreproductive (worker) females. Yet in many eusocial insects, workers retain the ability to produce haploid male offspring from unfertilized eggs. The reproductive potential of workers has well-documented consequences for the structure and function of insect colonies, but its implications at the population level are less often considered. We show that worker reproduction in honey bees can have an important role in maintaining genetic diversity at the sex locus in invasive populations. The honey bee sex locus is homozygous-lethal, and, all else being equal, a higher allele number in the population lead to higher mean brood survival. In an invasive population of the honey bee Apis cerana in Australia, workers contribute significantly to male production: 38% of male-producing colonies are queenless, and these contribute one-third of all males at mating congregations. Using a model, we show that such male production by queenless workers will increase the number of sex alleles retained in nascent invasive populations following founder events, relative to a scenario in which only queens reproduce. We conclude that by rescuing sex locus diversity that would otherwise be lost, workers' sons help honey bee populations to minimize the negative effects of inbreeding after founder events and so contribute to their success as invaders.

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.

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.

Anarchy Is a Molecular Signature of Worker Sterility in the Honey Bee.

Worker sterility is a defining characteristic of eusociality. The existence of the sterile worker caste remains a fundamental question for evolutionary biology as it requires the existence of genes that reduce personal reproduction. Currently, little is known about the proximate mechanisms underpinning worker sterility. Studies into a mutant "anarchistic" strain (in which workers can activate their ovaries) of honey bee, Apis mellifera, identified a list of candidate genes that regulate ovary activation. We quantified the expression of the four most promising candidate genes (Anarchy, Pdk1, S6k, and Ulk3) in nonactivated and activated ovaries of wild-type workers. Ovarian expression of Anarchy, a peroxisomal membrane protein, predicts the ovary state of workers with 88.2% accuracy. Increased expression of Anarchy in the ovary is strongly associated with suppression of oogenesis and its expression is sensitive to the presence of the queen. Therefore, Anarchy satisfies key criteria for a "gene underlying altruism". When we knocked down expression of Anarchy in the ovary using RNA interference (RNAi) we altered the expression of Buffy, a gene that regulates programmed cell death. Whole-mount multiplex fluorescent in situ hybridization (mFISH) shows Anarchy transcripts localize to degenerating oocytes within the ovary. Our results suggest that Anarchy is involved in the regulation of oogenesis through programmed cell death. The evolution of facultative worker sterility most likely occurred when the conserved mechanism of programmed cell death was co-opted to regulate ovary activation. Anarchy may therefore be the first example of a gene that has evolved through kin selection to regulate worker sterility.

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.

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.

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.

Only full-sibling families evolved eusociality.

Arising from M. A. Nowak, C. E. Tarnita & E. O. Wilson 466, 1057-1062 (2010); Nowak et al. reply. The paper by Nowak et al. has the evolution of eusociality as its title, but it is mostly about something else. It argues against inclusive fitness theory and offers an alternative modelling approach that is claimed to be more fundamental and general, but which, we believe, has no practical biological meaning for the evolution of eusociality. Nowak et al. overlook the robust empirical observation that eusociality has only arisen in clades where mothers are associated with their full-sibling offspring; that is, in families where the average relatedness of offspring to siblings is as high as to their own offspring, independent of population structure or ploidy. We believe that this omission makes the paper largely irrelevant for understanding the evolution of eusociality.

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.

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 CO2-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 CO2-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.