Four quantitative trait loci that influence worker sterility in the honeybee (Apis mellifera)

The microRNA miR-14 Regulates Egg-Laying by Targeting EcR in Honeybees (Apis mellifera)

Honeybees (Apis mellifera) are important pollinators and are commonly used for honey production. The oviposition behavior in honeybees is complex and errors in oviposition could affect the development of the bee colony. Recent studies reported that RNA-RNA cross-talk played a critical role in several biological processes, including reproduction. Ecdysone receptor (EcR) and miR-14 were previously reported to play important roles in egg-laying. Moreover, EcR was predicted to be the target gene of miR-14 and may form miR-14-EcR cross-talk. In this study, knocking down and overexpression of miR-14 and EcR in queen model were implemented. The effect of RNA expression of miR-14 and EcR on the number of eggs laid by honeybee queens were analyzed. Further, luciferase assay was used to confirm the target relation between miR-14 and 3'UTR of EcR. The results showed that the expression of miR-14 and EcR was associated with the number of eggs laid by queens. In specific, inhibition of miR-14 expression enhanced the number of eggs laid, while overexpression of EcR enhanced the number of eggs laid. Lastly, we determined that miR-14 directly targets the mRNA of EcR. These findings suggest that the cross-talk of miR-14-EcR plays an important role in the number of eggs laid by honeybee queens.

Genetics in the Honey Bee: Achievements and Prospects toward the Functional Analysis of Molecular and Neural Mechanisms Underlying Social Behaviors

The European honey bee is a model organism for studying social behaviors. Comprehensive analyses focusing on the differential expression profiles of genes between the brains of nurse bees and foragers, or in the mushroom bodies—the brain structure related to learning and memory, and multimodal sensory integration—has identified candidate genes related to honey bee behaviors. Despite accumulating knowledge on the expression profiles of genes related to honey bee behaviors, it remains unclear whether these genes actually regulate social behaviors in the honey bee, in part because of the scarcity of genetic manipulation methods available for application to the honey bee. In this review, we describe the genetic methods applied to studies of the honey bee, ranging from classical forward genetics to recently developed gene modification methods using transposon and CRISPR/Cas9. We then discuss future functional analyses using these genetic methods targeting genes identified by the preceding research. Because no particular genes or neurons unique to social insects have been found yet, further exploration of candidate genes/neurons correlated with sociality through comprehensive analyses of mushroom bodies in the aculeate species can provide intriguing targets for functional analyses, as well as insight into the molecular and neural bases underlying social behaviors.

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.

MicroRNA signatures characterizing caste-independent ovarian activity in queen and worker honeybees (Apis mellifera L.)

Queen and worker honeybees differ profoundly in reproductive capacity. The queen of this complex society, with 200 highly active ovarioles in each ovary, is the fertile caste, whereas the workers have approximately 20 ovarioles as a result of receiving a different diet during larval development. In a regular queenright colony, the workers have inactive ovaries and do not reproduce. However, if the queen is sensed to be absent, some of the workers activate their ovaries, producing viable haploid eggs that develop into males. Here, a deep-sequenced ovary transcriptome library of reproductive workers was used as supporting data to assess the dynamic expression of the regulatory molecules and microRNAs (miRNAs) of reproductive and nonreproductive honeybee females. In this library, most of the differentially expressed miRNAs are related to ovary physiology or oogenesis. When we quantified the dynamic expression of 19 miRNAs in the active and inactive worker ovaries and compared their expression in the ovaries of virgin and mated queens, we noted that some miRNAs (miR-1, miR-31a, miR-13b, miR-125, let-7 RNA, miR-100, miR-276, miR-12, miR-263a, miR-306, miR-317, miR-92a and miR-9a) could be used to identify reproductive and nonreproductive statuses independent of caste. Furthermore, integrative gene networks suggested that some candidate miRNAs function in the process of ovary activation in worker bees.

Structure and function of gene regulatory networks associated with worker sterility in honeybees

A characteristic of eusocial bees is a reproductive division of labor in which one or a few queens monopolize reproduction, while her worker daughters take on reproductively altruistic roles within the colony. The evolution of worker reproductive altruism involves indirect selection for the coordinated expression of genes that regulate personal reproduction, but evidence for this type of selection remains elusive. In this study, we tested whether genes coexpressed under queen‐induced worker sterility show evidence of adaptive organization within a model brain transcriptional regulatory network (TRN). If so, this structured pattern would imply that indirect selection on nonreproductive workers has influenced the functional organization of genes within the network, specifically to regulate the expression of sterility. We found that literature‐curated sets of candidate genes for sterility, ranging in size from 18 to 267, show strong evidence of clustering within the three‐dimensional space of the TRN. This finding suggests that our candidate sets of genes for sterility form functional modules within the living bee brain's TRN. Moreover, these same gene sets colocate to a single, albeit large, region of the TRN's topology. This spatially organized and convergent pattern contrasts with a null expectation for functionally unrelated genes to be haphazardly distributed throughout the network. Our meta‐genomic analysis therefore provides first evidence for a truly “social transcriptome” that may regulate the conditional expression of honeybee worker sterility.

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.

Transcriptome comparison between inactivated and activated ovaries of the honey bee Apis mellifera L

Ovarian activity not only influences fertility, but is also involved with the regulation of division of labour between reproductive and behavioural castes of female honey bees. In order to identify candidate genes associated with ovarian activity, we compared the gene expression patterns between inactivated and activated ovaries of queens and workers by means of high-throughput RNA-sequencing technology. A total of 1615 differentially expressed genes (DEGs) was detected between ovaries of virgin and mated queens, and more than 5300 DEGs were detected between inactivated and activated worker ovaries. Intersection analysis of DEGs amongst five libraries revealed that a similar set of genes (824) participated in the ovary activation of both queens and workers. A large number of these DEGs were predominantly related to cellular, cell and cell part, binding, biological regulation and metabolic processes. In addition, over 1000 DEGs were linked to more than 230 components of Kyoto Encyclopedia of Genes and Genomes pathways, including 25 signalling pathways. The reliability of the RNA-sequencing results was confirmed by means of quantitative real-time PCR. Our results provide new insights into the molecular mechanisms involved in ovary activation and reproductive division of labour.

Reproduction of honeybee workers is regulated by epidermal growth factor receptor signaling

Eusocial insect societies display a remarkable reproductive division of labor between a single fertile queen and thousands of largely sterile workers. In most species, however, the workers retain the capacity to reproduce, particularly in queenless colonies where typically many workers lay eggs. As yet, the molecular determinants that initiate this shift in worker fertility are still poorly documented. By using RNA interference we here demonstrate that the knockdown of epidermal growth factor receptor, a gene which was previously shown to be involved in queen-worker caste differentiation, also induces reproduction in worker honeybees (Apis mellifera). These data show that worker fertility and queen-worker caste determination partly rely on the same gene regulatory networks, thereby providing a major breakthrough in our understanding of the molecular determinants of the social insects' spectacular reproductive division of labor.

Genes underlying altruism

William D. Hamilton postulated the existence of 'genes underlying altruism', under the rubric of inclusive fitness theory, a half-century ago. Such genes are now poised for discovery. In this article, we develop a set of intuitive criteria for the recognition and analysis of genes for altruism and describe the first candidate genes affecting altruism from social insects and humans. We also provide evidence from a human population for genetically based trade-offs, underlain by oxytocin-system polymorphisms, between alleles for altruism and alleles for non-social cognition. Such trade-offs between self-oriented and altruistic behaviour may influence the evolution of phenotypic diversity across all social animals.

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

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

Kin conflict in insect societies: a new epigenetic perspective

The social hymenopterans (ants, wasps and bees) have all the enzymatic and genetic mechanisms necessary for the functional modification of DNA by methylation. Methylation appears to play a central role in shaping the developmental processes that give rise to the different castes. However, could DNA methylation have other roles in social insects? Theoretical arguments predict that male and female hymenopterans can be in conflict over the reproductive potential of their female offspring. An exciting prospect for future research is to examine the possibility that queens and males imprint the genomes of their gametes using DNA methylation to manipulate the reproductive potential of their progeny in ways that favour the inclusive fitness of the parent.

Genome-wide analysis of alternative reproductive phenotypes in honeybee workers

A defining feature of social insects is the reproductive division of labour, in which workers usually forego all reproduction to help their mother queen to reproduce. However, little is known about the molecular basis of this spectacular form of altruism. Here, we compared gene expression patterns between nonreproductive, altruistic workers and reproductive, non-altruistic workers in queenless honeybee colonies using a whole-genome microarray analysis. Our results demonstrate massive differences in gene expression patterns between these two sets of workers, with a total of 1292 genes being differentially expressed. In nonreproductive workers, genes associated with energy metabolism and respiration, flight and foraging behaviour, detection of visible light, flight and heart muscle contraction and synaptic transmission were overexpressed relative to reproductive workers. This implies they probably had a higher whole-body energy metabolism and activity rate and were most likely actively foraging, whereas same-aged reproductive workers were not. This pattern is predicted from evolutionary theory, given that reproductive workers should be less willing to compromise their reproductive futures by carrying out high-risk tasks such as foraging or other energetically expensive tasks. By contrast, reproductive workers mainly overexpressed oogenesis-related genes compared to nonreproductive ones. With respect to key switches for ovary activation, several genes involved in steroid biosynthesis were upregulated in reproductive workers, as well as genes known to respond to queen and brood pheromones, genes involved in TOR and insulin signalling pathways and genes located within quantitative trait loci associated with reproductive capacity in honeybees. Overall, our results provide unique insight into the molecular mechanisms underlying alternative reproductive phenotypes in honeybee workers.

Support for the reproductive ground plan hypothesis of social evolution and major QTL for ovary traits of Africanized worker honey bees (Apis mellifera L.)

BACKGROUND

The reproductive ground plan hypothesis of social evolution suggests that reproductive controls of a solitary ancestor have been co-opted during social evolution, facilitating the division of labor among social insect workers. Despite substantial empirical support, the generality of this hypothesis is not universally accepted. Thus, we investigated the prediction of particular genes with pleiotropic effects on ovarian traits and social behavior in worker honey bees as a stringent test of the reproductive ground plan hypothesis. We complemented these tests with a comprehensive genome scan for additional quantitative trait loci (QTL) to gain a better understanding of the genetic architecture of the ovary size of honey bee workers, a morphological trait that is significant for understanding social insect caste evolution and general insect biology.

RESULTS

Back-crossing hybrid European x Africanized honey bee queens to the Africanized parent colony generated two study populations with extraordinarily large worker ovaries. Despite the transgressive ovary phenotypes, several previously mapped QTL for social foraging behavior demonstrated ovary size effects, confirming the prediction of pleiotropic genetic effects on reproductive traits and social behavior. One major QTL for ovary size was detected in each backcross, along with several smaller effects and two QTL for ovary asymmetry. One of the main ovary size QTL coincided with a major QTL for ovary activation, explaining 3/4 of the phenotypic variance, although no simple positive correlation between ovary size and activation was observed.

CONCLUSIONS

Our results provide strong support for the reproductive ground plan hypothesis of evolution in study populations that are independent of the genetic stocks that originally led to the formulation of this hypothesis. As predicted, worker ovary size is genetically linked to multiple correlated traits of the complex division of labor in worker honey bees, known as the pollen hoarding syndrome. The genetic architecture of worker ovary size presumably consists of a combination of trait-specific loci and general regulators that affect the whole behavioral syndrome and may even play a role in caste determination. Several promising candidate genes in the QTL intervals await further study to clarify their potential role in social insect evolution and the regulation of insect fertility in general.

Hormone response to bidirectional selection on social behavior

Behavior is a quantitative trait determined by multiple genes. Some of these genes may have effects from early development and onward by influencing hormonal systems that are active during different life-stages leading to complex associations, or suites, of traits. Honey bees (Apis mellifera) have been used extensively in experiments on the genetic and hormonal control of complex social behavior, but the relationships between their early developmental processes and adult behavioral variation are not well understood. Bidirectional selective breeding on social food-storage behavior produced two honey bee strains, each with several sublines, that differ in an associated suite of anatomical, physiological, and behavioral traits found in unselected wild type bees. Using these genotypes, we document strain-specific changes during larval, pupal, and early adult life-stages for the central insect hormones juvenile hormone (JH) and ecdysteroids. Strain differences correlate with variation in female reproductive anatomy (ovary size), which can be influenced by JH during development, and with secretion rates of ecdysteroid from the ovaries of adults. Ovary size was previously assigned to the suite of traits of honey bee food-storage behavior. Our findings support that bidirectional selection on honey bee social behavior acted on pleiotropic gene networks. These networks may bias a bee's adult phenotype by endocrine effects on early developmental processes that regulate variation in reproductive traits.

Genetic architecture of ovary size and asymmetry in European honeybee workers

The molecular basis of complex traits is increasingly understood but a remaining challenge is to identify their co-regulation and inter-dependence. Pollen hoarding (pln) in honeybees is a complex trait associated with a well-characterized suite of linked behavioral and physiological traits. In European honeybee stocks bidirectionally selected for pln, worker (sterile helper) ovary size is pleiotropically affected by quantitative trait loci that were initially identified for their effect on foraging behavior. To gain a better understanding of the genetic architecture of worker ovary size in this model system, we analyzed a series of crosses between the selected strains. The crossing results were heterogeneous and suggested non-additive effects. Three significant and three suggestive quantitative trait loci of relatively large effect sizes were found in two reciprocal backcrosses. These loci are not located in genome regions of known effects on foraging behavior but contain several interesting candidate genes that may specifically affect worker-ovary size. Thus, the genetic architecture of this life history syndrome may be comprised of pleiotropic, central regulators that influence several linked traits and other genetic factors that may be downstream and trait specific.

Expression of genes related to reproduction and pollen foraging in honey bees (Apis mellifera) narcotized with carbon dioxide

It has been proposed that a honey bee (Apis mellifera) worker's preference for foraging for pollen or nectar is modulated by a gene network that was originally involved in regulating the reproductive cycles of an ancestral solitary species. We used carbon dioxide to induce narcosis in queens and workers. This treatment is known to initiate oogenesis in queens, reduce oogenesis in queenless workers and to change worker foraging preference. We then assessed changes in gene expression of genes suspected to be involved in either foraging behaviour or reproduction. We show that some genes change expression in the opposite direction between castes in response to treatment. Our results therefore support the hypothesis that reproductive and foraging traits are causally related in the honey bee.

Characterization of quantitative trait loci for the age of first foraging in honey bee workers

Identifying the basis of quantitative trait loci (QTL) remains challenging for the study of complex traits, such as behavior. The honey bee is a good model combining interesting social behavior with a high recombination rate that facilitates this identification. Several studies have focused on the pollen hoarding syndrome, identifying multiple QTL as the genetic basis of its behavioral components. One component, the age of first foraging, is central for colony organization and four QTL were previously described without identification of their genomic location. Enabled by the honey bee genome project, this study provides data from multiple experiments to scrutinize these QTL, including individual and pooled SNP mapping, sequencing of AFLP markers, and microsatellite genotyping. The combined evidence confirms and localizes two of the previous QTL on chromosome four and five, dismisses the other two, and suggests one novel genomic region on chromosome eleven to influence the age of first foraging. Among the positional candidates the Ank2, PKC, Erk7, and amontillado genes stand out due to corroborating functional evidence. This study thus demonstrates the power of combined, genome-based approaches to enable targeted studies of a manageable set of candidate genes for natural behavioral variation in the important, complex social trait "age of first foraging".

Genome-wide analysis of genes related to ovary activation in worker honey bees

A defining characteristic of eusocial animals is their division of labour into reproductive and nonreproductive specialists. Here, we used a microarray study to identify genes associated with functional sterility in the worker honey bee Apis mellifera. We contrasted gene expression in workers from a functionally sterile wild-type strain with that in a mutant (anarchist) strain selected for high rates of ovary activation. We identified a small set of genes from the brain (n = 7) and from the abdomen (n = 5) that are correlated in their expression with early stages of ovary activation. Sterile wild-type workers up-regulated two unknown genes and a homologue of Drosophila CG6004. By contrast, reproductive anarchist workers up-regulated genes for the yolk protein vitellogenin, venom peptides and a member of the AdoHycase superfamily, among others. The differentially expressed genes identified are likely to be involved in early differentiation into sterile and reproductive worker phenotypes and may therefore form part of the gene networks associated with the regulation of honey bee worker sterility. Our study may have lacked sufficient power to detect all but a minority of biologically relevant changes taking place; however, the differential expression of vitellogenin and a putative AdoHycase suggests that our screen has captured core reproductive genes and that ovary activation may involve an epigenetic mechanism.