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

Phenotypic and genetic characterization of Africanized Apis mellifera colonies with natural tolerance to Varroa destructor and contrasting defensive behavior

Africanized Apis mellifera colonies with promising characteristics for beekeeping have been detected in northern Argentina (subtropical climate) and are considered of interest for breeding programs. Integral evaluation of this feral material revealed high colony strength and resistance/tolerance to brood diseases. However, these Africanized honeybees (AHB) also showed variable negative behavioral traits for beekeeping, such as defensiveness, tendency to swarm and avoidance behavior. We developed a protocol for the selection of AHB stocks based on defensive behavior and characterized contrasting colonies for this trait using NGS technologies. For this purpose, population and behavioral parameters were surveyed throughout a beekeeping season in nine daughter colonies obtained from a mother colony (A1 mitochondrial haplotype) with valuable characteristics (tolerance to the mite Varroa destructor, high colony strength and low defensiveness). A Defensive Behavior Index was developed and tested in the colonies under study. Mother and two daughter colonies displaying contrasting defensive behavior were analyzed by ddRADseq. High-quality DNA samples were obtained from 16 workers of each colony. Six pooled samples, including two replicates of each of the three colonies, were processed. A total of 12,971 SNPs were detected against the reference genome of A. mellifera, 142 of which showed significant differences between colonies. We detected SNPs in coding regions, lncRNA, miRNA, rRNA, tRNA, among others. From the original data set, we also identified 647 SNPs located in protein-coding regions, 128 of which are related to 21 genes previously associated with defensive behavior, such as dop3 and dopR2, CaMKII and ADAR, obp9 and obp10, and members of the 5-HT family. We discuss the obtained results by considering the influence of polyandry and paternal lineages on the defensive behavior in AHB and provide baseline information to use this innovative molecular approach, ddRADseq, to assist in the selection and evaluation of honey bee stocks showing low defensive behavior for commercial uses.

Examining parent-of-origin effects on transcription and RNA methylation in mediating aggressive behavior in honey bees (Apis mellifera)

Conflict between genes inherited from the mother (matrigenes) and the father (patrigenes) is predicted to arise during social interactions among offspring if these genes are not evenly distributed among offspring genotypes. This intragenomic conflict drives parent-specific transcription patterns in offspring resulting from parent-specific epigenetic modifications. Previous tests of the kinship theory of intragenomic conflict in honey bees (Apis mellifera) provided evidence in support of theoretical predictions for variation in worker reproduction, which is associated with extreme variation in morphology and behavior. However, more subtle behaviors - such as aggression - have not been extensively studied. Additionally, the canonical epigenetic mark (DNA methylation) associated with parent-specific transcription in plant and mammalian model species does not appear to play the same role as in honey bees, and thus the molecular mechanisms underlying intragenomic conflict in this species is an open area of investigation. Here, we examined the role of intragenomic conflict in shaping aggression in honey bee workers through a reciprocal cross design and Oxford Nanopore direct RNA sequencing. We attempted to probe the underlying regulatory basis of this conflict through analyses of parent-specific RNA m6A and alternative splicing patterns. We report evidence that intragenomic conflict occurs in the context of honey bee aggression, with increased paternal and maternal allele-biased transcription in aggressive compared to non-aggressive bees, and higher paternal allele-biased transcription overall. However, we found no evidence to suggest that RNA m6A or alternative splicing mediate intragenomic conflict in this species.

Lineage and Parent-of-Origin Effects in DNA Methylation of Honey Bees (Apis mellifera) Revealed by Reciprocal Crosses and Whole-Genome Bisulfite Sequencing

Parent-of-origin methylation arises when the methylation patterns of a particular allele are dependent on the parent it was inherited from. Previous work in honey bees has shown evidence of parent-of-origin-specific expression, yet the mechanisms regulating such pattern remain unknown in honey bees. In mammals and plants, DNA methylation is known to regulate parent-of-origin effects such as genomic imprinting. Here, we utilize genotyping of reciprocal European and Africanized honey bee crosses to study genome-wide allele-specific methylation patterns in sterile and reproductive individuals. Our data confirm the presence of allele-specific methylation in honey bees in lineage-specific contexts but also importantly, though to a lesser degree, parent-of-origin contexts. We show that the majority of allele-specific methylation occurs due to lineage rather than parent-of-origin factors, regardless of the reproductive state. Interestingly, genes affected by allele-specific DNA methylation often exhibit both lineage and parent-of-origin effects, indicating that they are particularly labile in terms of DNA methylation patterns. Additionally, we re-analyzed our previous study on parent-of-origin-specific expression in honey bees and found little association with parent-of-origin-specific methylation. These results indicate strong genetic background effects on allelic DNA methylation and suggest that although parent-of-origin effects are manifested in both DNA methylation and gene expression, they are not directly associated with each other.

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?'

Tissue-specific transcription patterns support the kinship theory of intragenomic conflict in honey bees (Apis mellifera)

Kin selection may act differently on genes inherited from parents (matrigenes and patrigenes), resulting in intragenomic conflict. This conflict can be observed as differential expression of matrigenes and patrigenes, or parent-specific gene expression (PSGE). In honey bees (Apis mellifera), intragenomic conflict is hypothesized to occur in multiple social contexts. Previously, we found that patrigene-biased expression in reproductive tissues was associated with increased reproductive potential in worker honey bees, consistent with the prediction that patrigenes are selected to promote selfish behaviour in this context. Here, we examined brain gene expression patterns to determine if PSGE is also found in other tissues. As before, the number of transcripts showing patrigene expression bias was significantly greater in the brains of reproductive vs. sterile workers, while the number of matrigene-biased transcripts was not significantly different. Twelve transcripts out of the 374 showing PSGE in either tissue showed PSGE in both brain and reproductive tissues; this overlap was significantly greater than expected by chance. However, the majority of transcripts show PSGE only in one tissue, suggesting the epigenetic mechanisms mediating PSGE exhibit plasticity between tissues. There was no significant overlap between transcripts that showed PSGE and transcripts that were significantly differentially expressed. Weighted gene correlation network analysis identified modules which were significantly enriched in both types of transcripts, suggesting that these genes may influence each other through gene networks. Our results provide further support for the kin selection theory of intragenomic conflict, and provide valuable insights into the mechanisms which may mediate this process.

Parent of origin gene expression in the bumblebee, Bombus terrestris, supports Haig's kinship theory for the evolution of genomic imprinting

Genomic imprinting is the differential expression alleles in diploid individuals, with the expression being dependent on the sex of the parent from which it was inherited. Haig's kinship theory hypothesizes that genomic imprinting is due to an evolutionary conflict of interest between alleles from the mother and father. In social insects, it has been suggested that genomic imprinting should be widespread. One recent study identified parent-of-origin expression in honey bees and found evidence supporting the kinship theory. However, little is known about genomic imprinting in insects and multiple theoretical predictions must be tested to avoid single-study confirmation bias. We, therefore, tested for parent-of-origin expression in a primitively eusocial bee. We found equal numbers of maternally and paternally biased expressed genes. The most highly biased genes were maternally expressed, offering support for the kinship theory. We also found low conservation of potentially imprinted genes with the honey bee, suggesting rapid evolution of genomic imprinting in Hymenoptera.

A novel method for using RNA-seq data to identify imprinted genes in social Hymenoptera with multiply mated queens

Genomic imprinting results in parent-of-origin-dependent gene expression biased towards either the maternally or paternally derived allele at the imprinted locus. The kinship theory of genomic imprinting argues that this unusual expression pattern can be a manifestation of intra-genomic conflict between the maternally and paternally derived halves of the genome that arises because they are not equally related to the genomes of social partners. The theory thus predicts that imprinting may evolve wherever there are close interactions among asymmetrically related kin. The social Hymenoptera with permanent caste differentiation are suitable candidates for testing the kinship theory because haplodiploid sex determination creates strong relatedness asymmetries and nursing workers interact closely with kin. However, progress in the search for imprinted genes in the social Hymenoptera has been slow, in part because tests for imprinting rely on reciprocal crosses that are impossible in most species. Here, we develop a method to systematically search for imprinting in haplodiploid social insects without crosses, using instead samples of pooled individuals collected from natural colonies. We tested this protocol using data available for the leaf-cutting ant Acromyrmex echinatior, providing the first genome-wide search for imprinting in any ant. Although we identified several genes as potentially imprinted, none of the four genes tested could be verified as imprinted using digital droplet PCR, highlighting the need for higher quality genomic assemblies that accurately map duplicated genes.

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.

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.

A Single SNP Turns a Social Honey Bee (Apis mellifera) Worker into a Selfish Parasite

The evolution of altruism in complex insect societies is arguably one of the major transitions in evolution and inclusive fitness theory plausibly explains why this is an evolutionary stable strategy. Yet, workers of the South African Cape honey bee (Apis mellifera capensis) can reverse to selfish behavior by becoming social parasites and parthenogenetically producing female offspring (thelytoky). Using a joint mapping and population genomics approach, in combination with a time-course transcript abundance dynamics analysis, we show that a single nucleotide polymorphism at the mapped thelytoky locus (Th) is associated with the iconic thelytokous phenotype. Th forms a linkage group with the ecdysis-triggering hormone receptor (Ethr) within a nonrecombining region under strong selection in the genome. A balanced detrimental allele system plausibly explains why the trait is specific to A. m. capensis and cannot easily establish itself into genomes of other honey bee subspecies.

Do social insects support Haig's kin theory for the evolution of genomic imprinting?

Although numerous imprinted genes have been described in several lineages, the phenomenon of genomic imprinting presents a peculiar evolutionary problem. Several hypotheses have been proposed to explain gene imprinting, the most supported being Haig's kinship theory. This theory explains the observed pattern of imprinting and the resulting phenotypes as a competition for resources between related individuals, but despite its relevance it has not been independently tested. Haig's theory predicts that gene imprinting should be present in eusocial insects in many social scenarios. These lineages are therefore ideal for testing both the theory's predictions and the mechanism of gene imprinting. Here we review the behavioral evidence of genomic imprinting in eusocial insects, the evidence of a mechanism for genomic imprinting and finally we evaluate recent results showing parent of origin allele specific expression in honeybees in the light of Haig's theory.

A Genomic Imprinting Model of Termite Caste Determination: Not Genetic but Epigenetic Inheritance Influences Offspring Caste Fate

Eusocial insects exhibit the most striking example of phenotypic plasticity. There has been a long controversy over the factors determining caste development of individuals in social insects. Here we demonstrate that parental phenotypes influence the social status of offspring not through genetic inheritance but through genomic imprinting in termites. Our extensive field survey and genetic analysis of the termite Reticulitermes speratus show that its breeding system is inconsistent with a genetic caste determination model. We therefore developed a genomic imprinting model, in which queen- and king-specific epigenetic marks antagonistically influence sexual development of offspring. The model accounts for all known empirical data on caste differentiation of R. speratus and other related species. By conducting colony-founding experiments and additively incorporating relevant socio-environmental factors into our genomic imprinting model, we show the relative importance of genomic imprinting and environmental factors in caste determination. The idea of epigenetic inheritance of sexual phenotypes solves the puzzle of why parthenogenetically produced daughters carrying only maternal chromosomes exclusively develop into queens and why parental phenotypes (nymph- or worker-derived reproductives) strongly influence caste differentiation of offspring. According to our model, the worker caste is seen as a "neuter" caste whose sexual development is suppressed due to counterbalanced maternal and paternal imprinting and opens new avenues for understanding the evolution of caste systems in social insects.

Asymmetric paternal effect on offspring size linked to parent-of-origin expression of an insulin-like growth factor

Sexual reproduction brings together reproductive partners whose long‐term interests often differ, raising the possibility of conflict over their reproductive investment. Males that enhance maternal investment in their offspring gain fitness benefits, even if this compromises future reproductive investment by iteroparous females. When the conflict occurs at a genomic level, it may be uncovered by crossing divergent populations, as a mismatch in the coevolved patterns of paternal manipulation and maternal resistance may generate asymmetric embryonic growth. We report such an asymmetry in reciprocal crosses between populations of the fish Girardinichthys multiradiatus. We also show that a fragment of a gene which can influence embryonic growth (Insulin‐Like Growth Factor 2; igf2) exhibits a parent‐of‐origin methylation pattern, where the maternally inherited igf2 allele has much more 5′ cytosine methylation than the paternally inherited allele. Our findings suggest that male manipulation of maternal investment may have evolved in fish, while the parent‐of‐origin methylation pattern appears to be a potential candidate mechanism modulating this antagonistic coevolution process. However, disruption of other coadaptive processes cannot be ruled out, as these can lead to similar effects as conflict.

Pathogen-induced maternal effects result in enhanced immune responsiveness across generations

Parental investment theory postulates that adults can accurately perceive cues from their surroundings, anticipate the needs of future offspring based on those cues, and selectively allocate nongenetic resources to their progeny. Such context‐dependent parental contributions can result in phenotypically variable offspring. Consistent with these predictions, we show that bacterially exposed Manduca sexta mothers oviposited significantly more variable embryos (as measured by mass, volume, hatching time, and hatching success) relative to naïve and control mothers. By using an in vivo “clearance of infection” assay, we also show that challenged larvae born to heat‐killed‐ or live‐Serratia‐injected mothers, supported lower microbial loads and cleared the infection faster than progeny of control mothers. Our data support the notion that mothers can anticipate the future pathogenic risks and immunological needs of their unborn offspring, providing progeny with enhanced immune protection likely through transgenerational immune priming. Although the inclusion of live Serratia into oocytes does not appear to be the mechanism by which mothers confer protection to their young, other mechanisms, including epigenetic modifications in the progeny due to maternal pathogenic stress, may be at play. The adaptive nature of maternal effects in the face of pathogenic stress provides insights into parental investment, resource allocation, and life‐history theories and highlights the significant role that pathogen‐induced maternal effects play as generators and modulators of evolutionary change.

Sexual selection modulates genetic conflicts and patterns of genomic imprinting

Recent years have seen a surge of interest in linking the theories of kin selection and sexual selection. In particular, there is a growing appreciation that kin selection, arising through demographic factors such as sex-biased dispersal, may modulate sexual conflicts, including in the context of male-female arms races characterized by coevolutionary cycles. However, evolutionary conflicts of interest need not only occur between individuals, but may also occur within individuals, and sex-specific demography is known to foment such intragenomic conflict in relation to social behavior. Whether and how this logic holds in the context of sexual conflict-and, in particular, in relation to coevolutionary cycles-remains obscure. We develop a kin-selection model to investigate the interests of different genes involved in sexual and intragenomic conflict, and we show that consideration of these conflicting interests yields novel predictions concerning parent-of-origin specific patterns of gene expression and the detrimental effects of different classes of mutation and epimutation at loci underpinning sexually selected phenotypes.

Molecular Evolution of Insect Sociality: An Eco-Evo-Devo Perspective

The evolution of eusociality is a perennial issue in evolutionary biology, and genomic advances have fueled steadily growing interest in the genetic changes underlying social evolution. Along with a recent flurry of research on comparative and evolutionary genomics in different eusocial insect groups (bees, ants, wasps, and termites), several mechanistic explanations have emerged to describe the molecular evolution of eusociality from solitary behavior. These include solitary physiological ground plans, genetic toolkits of deeply conserved genes, evolutionary changes in protein-coding genes, cis regulation, and the structure of gene networks, epigenetics, and novel genes. Despite this proliferation of ideas, there has been little synthesis, even though these ideas are not mutually exclusive and may in fact be complementary. We review available data on molecular evolution of insect sociality and highlight key biotic and abiotic factors influencing social insect genomes. We then suggest both phylogenetic and ecological evolutionary developmental biology (eco-evo-devo) perspectives for a more synthetic view of molecular evolution in insect societies.

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.

Intragenomic Conflict over Soldier Allocation in Polyembryonic Parasitoid Wasps

Understanding the selection pressures that have driven the evolution of sterile insect castes has been the focus of decades of intense scientific debate. An amenable empirical test bed for theory on this topic is provided by the sterile-soldier caste of polyembryonic parasitoid wasps. The function of these soldiers has been a source of controversy, with two basic hypotheses emerging: the "brood-benefit" hypothesis that they provide an overall benefit for their siblings and the "sex-ratio-conflict" hypothesis that the soldiers mediate a conflict between brothers and sisters by killing their opposite-sex siblings. Here, we investigate the divergent sex-ratio optima of a female embryo's maternal-origin and paternal-origin genes, to determine the potential for, and direction of, intragenomic conflict over soldiering. We then derive contrasting empirically testable predictions concerning the patterns of genomic imprinting that are expected to arise out of this intragenomic conflict, for the brood-benefit versus the sex-ratio-conflict hypothesis of soldier function.

Testing the kinship theory of intragenomic conflict in honey bees (Apis mellifera)

Sexual reproduction brings genes from two parents (matrigenes and patrigenes) together into one individual. These genes, despite being unrelated, should show nearly perfect cooperation because each gains equally through the production of offspring. However, an individual's matrigenes and patrigenes can have different probabilities of being present in other relatives, so kin selection could act on them differently. Such intragenomic conflict could be implemented by partial or complete silencing (imprinting) of an allele by one of the parents. Evidence supporting this theory is seen in offspring-mother interactions, with patrigenes favoring acquisition of more of the mother's resources if some of the costs fall on half-siblings who do not share the patrigene. The kinship theory of intragenomic conflict is little tested in other contexts, but it predicts that matrigene-patrigene conflict may be rife in social insects. We tested the hypothesis that honey bee worker reproduction is promoted more by patrigenes than matrigenes by comparing across nine reciprocal crosses of two distinct genetic stocks. As predicted, hybrid workers show reproductive trait characteristics of their paternal stock, (indicating enhanced activity of the patrigenes on these traits), greater patrigenic than matrigenic expression, and significantly increased patrigenic-biased expression in reproductive workers. These results support both the general prediction that matrigene-patrigene conflict occurs in social insects and the specific prediction that honey bee worker reproduction is driven more by patrigenes. The success of these predictions suggests that intragenomic conflict may occur in many contexts where matrigenes and patrigenes have different relatednesses to affected kin.

Biased Allele Expression and Aggression in Hybrid Honeybees may be Influenced by Inappropriate Nuclear-Cytoplasmic Signaling

Hybrid effects are often exhibited asymmetrically between reciprocal families. One way this could happen is if silencing of one parent’s allele occurs in one lineage but not the other, which could affect the phenotypes of the hybrids asymmetrically by silencing that allele in only one of the hybrid families. We have previously tested for allele-specific expression biases in hybrids of European and Africanized honeybees and we found that there was an asymmetric overabundance of genes showing a maternal bias in the family with a European mother. Here, we further analyze allelic bias in these hybrids to ascertain whether they may underlie previously described asymmetries in metabolism and aggression in similar hybrid families and we speculate on what mechanisms may produce this biased allele usage. We find that there are over 500 genes that have some form of biased allele usage and over 200 of these are biased toward the maternal allele but only in the family with European maternity, mirroring the pattern observed for aggression and metabolic rate. This asymmetrically biased set is enriched for genes in loci associated with aggressive behavior and also for mitochondrial-localizing proteins. It contains many genes that play important roles in metabolic regulation. Moreover we find genes relating to the piwi-interacting RNA (piRNA) pathway, which is involved in chromatin modifications and epigenetic regulation and may help explain the mechanism underlying this asymmetric allele use. Based on these findings and previous work investigating aggression and metabolism in bees, we propose a novel hypothesis; that the asymmetric pattern of biased allele usage in these hybrids is a result of inappropriate use of piRNA-mediated nuclear-cytoplasmic signaling that is normally used to modulate aggression in honeybees. This is the first report of widespread asymmetric effects on allelic expression in hybrids and may represent a novel mechanism for gene regulation.

Differentially methylated obligatory epialleles modulate context-dependent LAM gene expression in the honeybee Apis mellifera

Differential intragenic methylation in social insects has been hailed as a prime mover of environmentally driven organismal plasticity and even as evidence for genomic imprinting. However, very little experimental work has been done to test these ideas and to prove the validity of such claims. Here we analyze in detail differentially methylated obligatory epialleles of a conserved gene encoding lysosomal α-mannosidase (AmLAM) in the honeybee. We combined genotyping of progenies derived from colonies founded by single drone inseminated queens, ultra-deep allele-specific bisulfite DNA sequencing, and gene expression to reveal how sequence variants, DNA methylation, and transcription interrelate. We show that both methylated and non-methylated states of AmLAM follow Mendelian inheritance patterns and are strongly influenced by polymorphic changes in DNA. Increased methylation of a given allele correlates with higher levels of context-dependent AmLAM expression and appears to affect the transcription of an antisense long noncoding RNA. No evidence of allelic imbalance or imprinting involved in this process has been found. Our data suggest that by generating alternate methylation states that affect gene expression, sequence variants provide organisms with a high level of epigenetic flexibility that can be used to select appropriate responses in various contexts. This study represents the first effort to integrate DNA sequence variants, gene expression, and methylation in a social insect to advance our understanding of their relationships in the context of causality.

A Search for Parent-of-Origin Effects on Honey Bee Gene Expression

Parent-specific gene expression (PSGE) is little known outside of mammals and plants. PSGE occurs when the expression level of a gene depends on whether an allele was inherited from the mother or the father. Kin selection theory predicts that there should be extensive PSGE in social insects because social insect parents can gain inclusive fitness benefits by silencing parental alleles in female offspring. We searched for evidence of PSGE in honey bees using transcriptomes from reciprocal crosses between European and Africanized strains. We found 46 transcripts with significant parent-of-origin effects on gene expression, many of which overexpressed the maternal allele. Interestingly, we also found a large proportion of genes showing a bias toward maternal alleles in only one of the reciprocal crosses. These results indicate that PSGE may occur in social insects. The nonreciprocal effects could be largely driven by hybrid incompatibility between these strains. Future work will help to determine if these are indeed parent-of-origin effects that can modulate inclusive fitness benefits.

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.

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.

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.