Gene expression profiles underlying alternative caste phenotypes in a highly eusocial bee, Melipona quadrifasciata

Eusociality is linked to caste-specific differences in metabolism, immune system, and somatic maintenance-related processes in an ant species

The social organization of many primate, bird and rodent species and the role of individuals within that organization are associated with specific individual physiological traits. However, this association is perhaps most pronounced in eusocial insects (e.g., termites, ants). In such species, genetically close individuals show significant differences in behavior, physiology, and life expectancy. Studies addressing the metabolic changes according to the social role are still lacking. We aimed at understanding how sociality could influence essential molecular processes in a eusocial insect, the black garden ant (Lasius niger) where queens can live up to ten times longer than workers. Using mass spectrometry-based analysis, we explored the whole metabolome of queens, nest-workers and foraging workers. A former proteomics study done in the same species allowed us to compare the findings of both approaches. Confirming the former results at the proteome level, we showed that queens had fewer metabolites related to immunity. Contrary to our predictions, we did not find any metabolite linked to reproduction in queens. Among the workers, foragers had a metabolic signature reflecting a more stressful environment and a more highly stimulated immune system. We also found that nest-workers had more digestion-related metabolites. Hence, we showed that specific metabolic signatures match specific social roles. Besides, we identified metabolites differently expressed among behavioral castes and involved in nutrient sensing and longevity pathways (e.g., sirtuins, FOXO). The links between such molecular pathways and aging being found in an increasing number of taxa, our results confirm and strengthen their potential universality.

Evaluation of reference genes for gene expression analysis by real-time quantitative PCR (qPCR) in three stingless bee species (Hymenoptera: Apidae: Meliponini)

Stingless bees are generalist pollinators distributed through the pantropical region. There is growing evidence that their wild populations are experiencing substantial decline in response to habitat degradation and pesticides. Policies for conservation of endangered species will benefit from studies focusing on genetic and molecular aspects of their development and behavior. The most common method for looking at gene expression is real-time quantitative polymerase chain reaction preceded by reverse transcription (RT-qPCR) of the mRNA of interest. This method requires the identification of reliable reference genes to correctly estimate fluctuations in transcript levels. To contribute to molecular studies on stingless bees, we used Frieseomelitta varia, Melipona quadrifasciata, and Scaptotrigona bipunctata species to test the expression stability of eight reference genes (act, ef1-α, gapdh, rpl32, rps5, rps18, tbp, and tbp-af) in RT-qPCR procedures in five physiological and experimental conditions (development, sex, tissues, bacteria injection, and pesticide exposure). In general, the rpl32, rps5 and rps18 ribosomal protein genes and tpb-af gene showed the highest stability, thus being identified as suitable reference genes for the three stingless bee species and defined conditions. Our results also emphasized the need to evaluate the stability of candidate genes for any designed experimental condition and stingless bee species.

Division of labour in the black garden ant (Lasius niger) leads to three distinct proteomes

Task specialization in social insects leads to striking intra-specific differences in behaviour, morphology, physiology and longevity, but the underlying mechanisms remain not yet fully understood. Adult colonies of black garden ants (Lasius niger) have a single queen fertilized by one or a small number of males. The inter-individual genetic variability is thus relatively low, making it easier to focus on the individual molecular differences linked to the division of labour. Mass spectrometry-based proteomics enabled us to highlight which biological functions create the difference between queens, foragers and nest-workers. The proteome of each caste reflected nicely their social role: e.g., reproduction for queens, pesticide resistance for foragers - that are the most exposed to environmental risk factors - and, interestingly, digestion for nest-workers, thus highlighting proteomic profiles differences even among workers. Furthermore, our exploratory approach suggests energy trade-off mechanisms - in connection with the theory of social immunity - that might explain the difference in longevity between queens and workers. This study brings evidence that proteomics is able to highlight the subtle mechanisms of molecular regulation induced by social organization.

Reproductive capacity and castes in eusocial stingless bees (Hymenoptera: Apidae)

Eusocial lifestyle is one of the most important transitions in the evolutionary history of some groups of organisms. In bees, there are only two eusocial groups: the honey bees (Apini) and the stingless bees (Meliponini). Despite similarities on the eusocial lifestyles of these taxa, they present profound differences related to caste determination, development, behavior, and reproductive capacity of their members. In most of them the queen has a monopoly on reproduction. However, even though workers are tipically sterile, they can contribute to producing haploid eggs that generate males, or trophic eggs, used as an additional nutrition by the queen.

Transcriptome Profile of the Asian Giant Hornet (Vespa mandarinia) Using Illumina HiSeq 4000 Sequencing: De Novo Assembly, Functional Annotation, and Discovery of SSR Markers

Vespa mandarinia found in the forests of East Asia, including Korea, occupies the highest rank in the arthropod food web within its geographical range. It serves as a source of nutrition in the form of Vespa amino acid mixture and is listed as a threatened species, although no conservation measures have been implemented. Here, we performed de novo assembly of the V. mandarinia transcriptome by Illumina HiSeq 4000 sequencing. Over 60 million raw reads and 59,184,811 clean reads were obtained. After assembly, a total of 66,837 unigenes were clustered, 40,887, 44,455, and 22,390 of which showed homologous matches against the PANM, Unigene, and KOG databases, respectively. A total of 15,675 unigenes were assigned to Gene Ontology terms, and 5,132 unigenes were mapped to 115 KEGG pathways. The zinc finger domain (C2H2-like), serine/threonine/dual specificity protein kinase domain, and RNA recognition motif domain were among the top InterProScan domains predicted for V. mandarinia sequences. Among the unigenes, we identified 534,922 cDNA simple sequence repeats as potential markers. This is the first transcriptomic analysis of the wasp V. mandarinia using Illumina HiSeq 4000. The obtained datasets should promote the search for new genes to understand the physiological attributes of this wasp.

Insights into the dynamics of hind leg development in honey bee (Apis mellifera L.) queen and worker larvae - A morphology/differential gene expression analysis

Phenotypic plasticity is a hallmark of the caste systems of social insects, expressed in their life history and morphological traits. These are best studied in bees. In their co-evolution with angiosperm plants, the females of corbiculate bees have acquired a specialized structure on their hind legs for collecting pollen. In the highly eusocial bees (Apini and Meliponini), this structure is however only present in workers and absent in queens. By means of histological sections and cell proliferation analysis we followed the developmental dynamics of the hind legs of queens and workers in the fourth and fifth larval instars. In parallel, we generated subtractive cDNA libraries for hind leg discs of queen and worker larvae by means of a Representational Difference Analysis (RDA). From the total of 135 unique sequences we selected 19 for RT-qPCR analysis, where six of these were confirmed as differing significantly in their expression between the two castes in the larval spinning stage. The development of complex structures such as the bees' hind legs, requires diverse patterning mechanisms and signaling modules, as indicated by the set of differentially expressed genes related with cell adhesion and signaling pathways.

Identification of differentially expressed genes in female Drosophila antonietae and Drosophila meridionalis in response to host cactus odor

Background

Studies of insect-plant interactions have provided critical insights into the ecology and evolution of adaptive processes within and among species. Cactophilic Drosophila species have received much attention because larval development occurs in the necrotic tissues of cacti, and both larvae and adults feed on these tissues. Such Drosophila-cactus interactions include effects of the host plant on the physiology and behavior of the flies, especially so their nutritional status, mating condition and reproduction. The aim of this work was to compare the transcriptional responses of two species, Drosophila antonietae and Drosophila meridionalis, and identify genes potentially related to responses to odors released by their host cactus, Cereus hildmannianus. The two fly species are sympatric in most of their populations and use this same host cactus in nature.

Results

We obtained 47 unique sequences (USs) for D. antonietae in a suppression subtractive hybridization screen, 30 of these USs had matches with genes predicted for other Drosophila species. For D. meridionalis we obtained 81 USs, 46 of which were orthologous with genes from other Drosophila species. Functional information (Gene Ontology) revealed that these differentially expressed genes are related to metabolic processes, detoxification mechanisms, signaling, response to stimuli, and reproduction. The expression of 13 genes from D. meridionalis and 12 from D. antonietae were further analyzed by quantitative real time-PCR, showing that four genes were significantly overexpressed in D. antonietae and six in D. meridionalis.

Conclusions

Our results revealed the differential expression of genes related to responses to odor stimuli by a cactus, in two associated fly species. Although the majority of activated genes were similar between the two species, we also observed that certain metabolic pathways were specifically activated, especially those related to signaling pathways and detoxification mechanisms. The activation of these genes may reflect different metabolic pathways used by these flies in their interaction with this host cactus. Our findings provide insight into how the use of C. hildmannianus may have arisen independently in the two fly species, through genetic differentiation in metabolic pathways to effectively explore this cactus as a host.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0191-2) contains supplementary material, which is available to authorized users.

Transcriptome analysis of the Asian honey bee Apis cerana cerana

Background

The Eastern hive honey bee, Apis cerana cerana is a native and widely bred honey bee species in China. Molecular biology research about this honey bee species is scarce, and genomic information for A. c. cerana is not currently available. Transcriptome and expression profiling data for this species are therefore important resources needed to better understand the biological mechanisms of A. c. cerana. In this study, we obtained the transcriptome information of A. c. cerana by RNA-sequencing and compared gene expression differences between queens and workers of A. c. cerana by digital gene expression (DGE) analysis.

Results

Using high-throughput Illumina RNA sequencing we obtained 51,581,510 clean reads corresponding to 4.64 Gb total nucleotides from a single run. These reads were assembled into 46,999 unigenes with a mean length of 676 bp. Based on a sequence similarity search against the five public databases (NR, Swissport, GO, COG, KEGG) with a cut-off E-value of 10⁻⁵ using BLASTX, a total of 24,630 unigenes were annotated with gene descriptions, gene ontology terms, or metabolic pathways. Using these transcriptome data as references we analyzed the gene expression differences between the queens and workers of A. c. cerana using a tag-based digital gene expression method. We obtained 5.96 and 5.66 million clean tags from the queen and worker samples, respectively. A total of 414 genes were differentially expressed between them, with 189 up-regulated and 225 down-regulated in queens.

Conclusions

Our transcriptome data provide a comprehensive sequence resource for future A. c. cerana study, establishing an important public information platform for functional genomic studies in A. c. cerana. Furthermore, the DGE data provide comprehensive gene expression information for the queens and workers, which will facilitate our understanding of the molecular mechanisms of the different physiological aspects of the two castes.

Transcript levels of ten caste-related genes in adult diploid males of Melipona quadrifasciata (Hymenoptera, Apidae) - A comparison with haploid males, queens and workers

In Hymenoptera, homozygosity at the sex locus results in the production of diploid males. In social species, these pose a double burden by having low fitness and drawing resources normally spent for increasing the work force of a colony. Yet, diploid males are of academic interest as they can elucidate effects of ploidy (normal males are haploid, whereas the female castes, the queens and workers, are diploid) on morphology and life history. Herein we investigated expression levels of ten caste-related genes in the stingless bee Melipona quadrifasciata, comparing newly emerged and 5-day-old diploid males with haploid males, queens and workers. In diploid males, transcript levels for dunce and paramyosin were increased during the first five days of adult life, while those for diacylglycerol kinase and the transcriptional co-repressor groucho diminished. Two general trends were apparent, (i) gene expression patterns in diploid males were overall more similar to haploid ones and workers than to queens, and (ii) in queens and workers, more genes were up-regulated after emergence until day five, whereas in diploid and especially so in haploid males more genes were down-regulated. This difference between the sexes may be related to longevity, which is much longer in females than in males.

Alternative splicing of a single transcription factor drives selfish reproductive behavior in honeybee workers (Apis mellifera)

In eusocial insects the production of daughters is generally restricted to mated queens, and unmated workers are functionally sterile. The evolution of this worker sterility has been plausibly explained by kin selection theory [Hamilton W (1964) J Theor Biol 7:1-52], and many traits have evolved to prevent conflict over reproduction among the females in an insect colony. In honeybees (Apis mellifera), worker reproduction is regulated by the queen, brood pheromones, and worker policing. However, workers of the Cape honeybee, Apis mellifera capensis, can evade this control and establish themselves as social parasites by activating their ovaries, parthenogenetically producing diploid female offspring (thelytoky) and producing queen-like amounts of queen pheromones. All these traits have been shown to be strongly influenced by a single locus on chromosome 13 [Lattorff HMG, et al. (2007) Biol Lett 3:292-295]. We screened this region for candidate genes and found that alternative splicing of a gene homologous to the gemini transcription factor of Drosophila controls worker sterility. Knocking out the critical exon in a series of RNAi experiments resulted in rapid worker ovary activation-one of the traits characteristic of the social parasites. This genetic switch may be controlled by a short intronic splice enhancer motif of nine nucleotides attached to the alternative splice site. The lack of this motif in parasitic Cape honeybee clones suggests that the removal of nine nucleotides from the altruistic worker genome may be sufficient to turn a honeybee from an altruistic worker into a parasite.

Representational Difference Analysis (RDA) reveals differential expression of conserved as well as novel genes during caste-specific development of the honey bee (Apis mellifera L.) ovary

In highly eusocial insects, such as the honey bee, Apis mellifera, the reproductive bias has become embedded in morphological caste differences. These are most expressively denoted in ovary size, with adult queens having large ovaries consisting of 150-200 ovarioles each, while workers typically have only 1-20 ovarioles per ovary. This morphological differentiation is a result of hormonal signals triggered by the diet change in the third larval instar, which eventually generate caste-specific gene expression patterns. To reveal these we produced differential gene expression libraries by Representational Difference Analysis (RDA) for queen and worker ovaries in a developmental stage when cell death is a prominent feature in the ovarioles of workers, whereas all ovarioles are maintained and extend in length in queens. In the queen library, 48% of the gene set represented homologs of known Drosophila genes, whereas in the worker ovary, the largest set (59%) were ESTs evidencing novel genes, not even computationally predicted in the honey bee genome. Differential expression was confirmed by quantitative RT-PCR for a selected gene set, denoting major differences for two queen and two worker library genes. These included two unpredicted genes located in chromosome 11 (Group11.35 and Group11.31, respectively) possibly representing long non-coding RNAs. Being candidates as modulators of ovary development, their expression and functional analysis should be a focal point for future studies.

Developmental decoupling of alternative phenotypes: insights from the transcriptomes of horn-polyphenic beetles

Developmental mechanisms play an important role in determining the costs, limits, and evolutionary consequences of phenotypic plasticity. One issue central to these claims is the hypothesis of developmental decoupling, where alternate morphs result from evolutionarily independent developmental pathways. We address this assumption through a microarray study that tests whether differences in gene expression between alternate morphs are as divergent as those between sexes, a classic example of developmental decoupling. We then examine whether genes with morph-biased expression are less conserved than genes with shared expression between morphs, as predicted if developmental decoupling relaxes pleiotropic constraints on divergence. We focus on the developing horns and brains of two species of horned beetles with impressive sexual- and morph-dimorphism in the expression of horns and fighting behavior. We find that patterns of gene expression were as divergent between morphs as they were between sexes. However, overall patterns of gene expression were also highly correlated across morphs and sexes. Morph-biased genes were more evolutionarily divergent, suggesting a role of relaxed pleiotropic constraints or relaxed selection. Together these results suggest that alternate morphs are to some extent developmentally decoupled, and that this decoupling has significant evolutionary consequences. However, alternative morphs may not be as developmentally decoupled as sometimes assumed and such hypotheses of development should be revisited and refined.

Next-generation small RNA sequencing for microRNAs profiling in the honey bee Apis mellifera

MicroRNAs (miRNAs) are key regulators in various physiological and pathological processes via post-transcriptional regulation of gene expression. The honey bee (Apis mellifera) is a key model for highly social species, and its complex social behaviour can be interpreted theoretically as changes in gene regulation, in which miRNAs are thought to be involved. We used the SOLiD sequencing system to identify the repertoire of miRNAs in the honey bee by sequencing a mixed small RNA library from different developmental stages. We obtained a total of 36,796,459 raw sequences; of which 5,491,100 short sequences were fragments of mRNA and other noncoding RNAs (ncRNA), and 1,759,346 reads mapped to the known miRNAs. We predicted 267 novel honey bee miRNAs representing 380,182 short reads, including eight miRNAs of other insects in 14,107,583 genome-mapped sequences. We verified 50 of them using stem-loop reverse-transcription PCR (RT-PCR), in which 35 yielded PCR products. Cross-species analyses showed 81 novel miRNAs with homologues in other insects, suggesting that they were authentic miRNAs and have similar functions. The results of this study provide a basis for studies of the miRNA-modulating networks in development and some intriguing phenomena such as caste differentiation in A. mellifera.

Suppression subtractive hybridization analysis reveals expression of conserved and novel genes in male accessory glands of the ant Leptothorax gredleri

Background

During mating, insect males eject accessory gland proteins (Acps) into the female genital tract. These substances are known to affect female post-mating behavior and physiology. In addition, they may harm the female, e.g., in reducing its lifespan. This is interpreted as a consequence of sexual antagonistic co-evolution. Whereas sexual conflict abounds in non-social species, the peculiar life history of social insects (ants, bees, wasps) with lifelong pair-bonding and no re-mating aligns the reproductive interests of the sexes. Harming the female during mating would negatively affect male fitness and sexual antagonism is therefore not expected. Indeed, mating appears to increase female longevity in at least one ant species. Acps are presumed to play a role in this phenomenon, but the underlying mechanisms are unknown. In this study, we investigated genes, which are preferentially expressed in male accessory glands of the ant Leptothorax gredleri, to determine which proteins might be transferred in the seminal fluid.

Results

By a suppression subtractive hybridization protocol we obtained 20 unique sequences (USs). Twelve had mutual best matches with genes predicted for Apis mellifera and Nasonia vitripennis. Functional information (Gene Ontology) was available only for seven of these, including intracellular signaling, energy-dependent transport and metabolic enzyme activities. The remaining eight USs did not match sequences from other species. Six genes were further analyzed by quantitative RT-PCR in three life cycle stages of male ants. A gene with carboxy-lyase activity and one of unpredicted function were significantly overexpressed in accessory glands of sexually mature males.

Conclusions

Our study is the first one to investigate differential gene expression in ants in a context related to mating. Our findings indicate that male accessory glands of L. gredleri express a series of genes that are unique to this species, possibly representing novel genes, in addition to conserved ones for which functions can be predicted. Identifying differentially expressed genes might help to better understand molecular mechanisms involved in reproductive processes in eusocial Hymenoptera. While the novel genes could account for rapidly evolving ones driven by intra-sexual conflict between males, conserved genes imply that rather beneficial traits might get fixed by a process described as inter-sexual cooperation between males and females.

Differential gene expression and protein abundance evince ontogenetic bias toward castes in a primitively eusocial wasp

Polistes paper wasps are models for understanding conditions that may have characterized the origin of worker and queen castes and, therefore, the origin of paper wasp sociality. Polistes is “primitively eusocial” by virtue of having context-dependent caste determination and no morphological differences between castes. Even so, Polistes colonies have a temporal pattern in which most female larvae reared by the foundress become workers, and most reared by workers become future-reproductive gynes. This pattern is hypothesized to reflect development onto two pathways, which may utilize mechanisms that regulate diapause in other insects. Using expressed sequence tags (ESTs) for Polistes metricus we selected candidate genes differentially expressed in other insects in three categories: 1) diapause vs. non-diapause phenotypes and/or worker vs. queen differentiation, 2) behavioral subcastes of worker honey bees, and 3) no a priori expectation of a role in worker/gyne development. We also used a non-targeted proteomics screen to test for peptide/protein abundance differences that could reflect larval developmental divergence. We found that foundress-reared larvae (putative worker-destined) and worker-reared larvae (putative gyne-destined) differed in quantitative expression of sixteen genes, twelve of which were associated with caste and/or diapause in other insects, and they also differed in abundance of nine peptides/proteins. Some differentially-expressed genes are involved in diapause regulation in other insects, and other differentially-expressed genes and proteins are involved in the insulin signaling pathway, nutrient metabolism, and caste determination in highly social bees. Differential expression of a gene and a peptide encoding hexameric storage proteins is especially noteworthy. Although not conclusive, our results support hypotheses of 1) larval developmental pathway divergence that can lead to caste bias in adults and 2) nutritional differences as the foundation of the pathway divergence. Finally, the differential expression in Polistes larvae of genes and proteins also differentially expressed during queen vs. worker caste development in honey bees may indicate that regulatory mechanisms of caste outcomes share similarities between primitively eusocial and advanced eusocial Hymenoptera.

Phenotypic plasticity and diversity in insects

Phenotypic plasticity in general and polyphenic development in particular are thought to play important roles in organismal diversification and evolutionary innovation. Focusing on the evolutionary developmental biology of insects, and specifically that of horned beetles, I explore the avenues by which phenotypic plasticity and polyphenic development have mediated the origins of novelty and diversity. Specifically, I argue that phenotypic plasticity generates novel targets for evolutionary processes to act on, as well as brings about trade-offs during development and evolution, thereby diversifying evolutionary trajectories available to natural populations. Lastly, I examine the notion that in those cases in which phenotypic plasticity is underlain by modularity in gene expression, it results in a fundamental trade-off between degree of plasticity and mutation accumulation. On one hand, this trade-off limits the extent of plasticity that can be accommodated by modularity of gene expression. On the other hand, it causes genes whose expression is specific to rare environments to accumulate greater variation within species, providing the opportunity for faster divergence and diversification between species, compared with genes expressed across environments. Phenotypic plasticity therefore contributes to organismal diversification on a variety of levels of biological organization, thereby facilitating the evolution of novel traits, new species and complex life cycles.

Toward a population genetic framework of developmental evolution: the costs, limits, and consequences of phenotypic plasticity

Adaptive phenotypic plasticity allows organisms to cope with environmental variability, and yet, despite its adaptive significance, phenotypic plasticity is neither ubiquitous nor infinite. In this review, we merge developmental and population genetic perspectives to explore costs and limits on the evolution of plasticity. Specifically, we focus on the role of modularity in developmental genetic networks as a mechanism underlying phenotypic plasticity, and apply to it lessons learned from population genetic theory on the interplay between relaxed selection and mutation accumulation. We argue that the environmental specificity of gene expression and the associated reduction in pleiotropic constraints drive a fundamental tradeoff between the range of plasticity that can be accommodated and mutation accumulation in alternative developmental networks. This tradeoff has broad implications for understanding the origin and maintenance of plasticity and may contribute to a better understanding of the role of plasticity in the origin, diversification, and loss of phenotypic diversity.

Aphid wing dimorphisms: linking environmental and genetic control of trait variation

Both genetic and environmental factors underlie phenotypic variation. While research at the interface of evolutionary and developmental biology has made excellent advances in understanding the contribution of genes to morphology, less well understood is the manner in which environmental cues are incorporated during development to influence the phenotype. Also virtually unexplored is how evolutionary transitions between environmental and genetic control of trait variation are achieved. Here, I review investigations into molecular mechanisms underlying phenotypic plasticity in the aphid wing dimorphism system. Among aphids, some species alternate between environmentally sensitive (polyphenic) and genetic (polymorphic) control of wing morph determination in their life cycle. Therefore, a traditional molecular genetic approach into understanding the genetically controlled polymorphism may provide a unique avenue into not only understanding the molecular basis of polyphenic variation in this group, but also the opportunity to compare and contrast the mechanistic basis of environmental and genetic control of similar dimorphisms.

Genetic and genomic analyses of the division of labour in insect societies

Division of labour--individuals specializing in different activities--features prominently in the spectacular success of the social insects. Until recently, genetic and genomic analyses of division of labour were limited to just a few species. However, research on an ever-increasing number of species has provided new insight, from which we highlight two results. First, heritable influences on division of labour are more pervasive than previously imagined. Second, different forms of division of labour, in lineages in which eusociality has arisen independently, have evolved through changes in the regulation of highly conserved molecular pathways associated with several basic life-history traits, including nutrition, metabolism and reproduction.

Reference Gene Selection for Insect Expression Studies Using Quantitative Real-Time PCR: The Head of the Honeybee, Apis mellifera, After a Bacterial Challenge

In this study an important and often neglected aspect of gene expression studies in insects, the validation of appropriate reference genes with stable expression levels between sample groups, is addressed. Although in this paper the reference gene selection for the honeybee, Apis mellifera L. (Hymenoptera: Apidae) head was tested in the context of bacterial challenge with Escherichia coli, this work can serve as a resource to help select and screen insect reference genes for gene expression studies in any tissue and under any experimental manipulation. Since it is recommended to use multiple reference genes for accurate normalization, we analyzed the expression of eleven candidate reference genes in the honeybee head, for their potential use in the analysis of differential gene expression following bacterial challenge. Three software programs, BestKeeper, Normfinder and geNorm, were used to assess candidate reference genes. GeNorm recommended the use of four reference genes. Both geNorm and Normfinder identified the genes GAPDH, RPS18, actin and RPL13a as the most stable ones, only differing in their ranking order. BestKeeper identified RPS18 as being the reference gene with the least overall variation, but also actin and GAPDH were found to be the second and third most stable expressed gene. By a combination of three software programs the genes actin, RPS18 and GAPDH were found suitable reference genes in the honeybee head in the context of bacterial infection.

Control of reproductive dominance by the thelytoky gene in honeybees

Differentiation into castes and reproductive division of labour are a characteristics of eusocial insects. Caste determination occurs at an early stage of larval development in social bees and is achieved via differential nutrition irrespective of the genotype. Workers are usually subordinate to the queen and altruistically refrain from reproduction. Workers of the Cape honeybee (Apis mellifera capensis) do not necessarily refrain from reproduction. They have the unique ability to produce female offspring parthenogenetically (thelytoky) and can develop into 'pseudoqueens'. Although these are morphologically workers, they develop a queen-like phenotype with respect to physiology and behaviour. Thelytoky is determined by a single gene (th) and we show that this gene also influences other traits related to the queen phenotype, including egg production and queen pheromone synthesis. Using 566 microsatellite markers, we mapped this gene to chromosome 13 and identified a candidate locus thelytoky, similar to grainy head (a transcription factor), which has been shown to be highly expressed in queens of eusocial insects. We therefore suggest that this gene is not only important for determining the pseudoqueen phenotype in A. m. capensis workers, but is also of general importance in regulating the gene cascades controlling reproduction and sterility in female social bees.

Genome-wide analysis reveals differences in brain gene expression patterns associated with caste and reproductive status in honey bees (Apis mellifera)

A key characteristic of eusocial species is reproductive division of labour. Honey bee colonies typically have a single reproductive queen and thousands of sterile workers. Adult queens differ dramatically from workers in anatomy, physiology, behaviour and lifespan. Young female workers can activate their ovaries and initiate egg laying; these 'reproductive' workers differ from sterile workers in anatomy, physiology, and behaviour. These differences, however, are on a much smaller scale than those observed between the queen and worker castes. Here, we use microarrays to monitor expression patterns of several thousand genes in the brains of same-aged virgin queens, sterile workers, and reproductive workers. We found large differences in expression between queens and both worker groups (~2000 genes), and much smaller differences between sterile and reproductive workers (221 genes). The expression patterns of these 221 genes in reproductive workers are more queen-like, and may represent a core group of genes associated with reproductive physiology. Furthermore, queens and reproductive workers preferentially up-regulate genes associated with the nurse bee behavioural state, which supports the hypothesis of an evolutionary link between worker division of labour and molecular pathways related to reproduction. Finally, several functional groups of genes associated with longevity in other species are significantly up-regulated in queens. Identifying the genes that underlie the differences between queens, sterile workers, and reproductive workers will allow us to begin to characterize the molecular mechanisms underlying the evolution of social behaviour and large-scale remodelling of gene networks associated with polyphenisms.

Endocrine analysis in evolutionary-developmental studies of insect polymorphism: hormone manipulation versus direct measurement of hormonal regulators

"Hormone manipulation" is being used increasingly in evo-devo studies as the sole or primary technique to investigate the regulation of insect polymorphism by hormones, most notably juvenile hormone (JH). This manuscript critically evaluates the limitations and strengths of this indirect method for inferring aspects of endocrine regulation, and conclusions derived from recent endocrine studies of evolution and development in which data have been obtained primarily or exclusively by this method. The main conclusions of this critique are as follows: first, when used alone, or as the primary empirical technique, hormone manipulation is a superficial method that is fraught with problems with respect to identifying a hormone that regulates developmental-morphological variation, let alone identifying its mode of action. Second, conclusions reported in studies using this technique as the exclusive, or nearly exclusive experimental approach, most notably recent studies of JH regulation of horn polymorphism in dung beetles, and some studies of wing polymorphism should be considered, at best, weakly supported until substantiated by well-validated, direct methods. Finally, there are many reliable and well-validated techniques that can be used to directly and accurately quantify JH levels, and activities of JH regulators, in many insects, even in small, nonmodel species. Some of the most important of these assays will be briefly described and their strengths and weaknesses will be discussed.

Molecular basis for the reproductive division of labour in a lower termite

Background

Polyphenism, the expression of different phenotypes with the same genetic background, is well known for social insects. The substantial physiological and morphological differences among the castes generally are the result of differential gene expression. In lower termites, workers are developmentally flexible to become neotenic replacement reproductives via a single moult after the death of the founding reproductives. Thus, both castes (neotenics and workers) are expected to differ mainly in the expression of genes linked to reproductive division of labour, which constitutes the fundamental basis of insect societies.

Results

Representational difference analysis of cDNAs was used to study differential gene expression between neotenics and workers in the drywood termite Cryptotermes secundus (Kalotermitidae). We identified and, at least partially cloned five novel genes that were highly expressed in female neotenics. Quantitative real-time PCR analysis of all five genes in different castes (neotenics, founding reproductives, winged sexuals and workers of both sexes) confirmed the differential expression patterns. In addition, the relative expression of these genes was determined in three body parts of female neotenics (head, thorax, and abdomen) using quantitative real-time PCR.

Conclusion

The identified genes could be involved in the control and regulation of reproductive division of labour. Interestingly, this study revealed an expression pattern partly similar to social Hymenoptera indicating both common and species-specific regulatory mechanisms in hemimetabolous and holometabolous social insects.

Molecular determinants of caste differentiation in the highly eusocial honeybee Apis mellifera

Background

In honeybees, differential feeding of female larvae promotes the occurrence of two different phenotypes, a queen and a worker, from identical genotypes, through incremental alterations, which affect general growth, and character state alterations that result in the presence or absence of specific structures. Although previous studies revealed a link between incremental alterations and differential expression of physiometabolic genes, the molecular changes accompanying character state alterations remain unknown.

Results

By using cDNA microarray analyses of >6,000 Apis mellifera ESTs, we found 240 differentially expressed genes (DEGs) between developing queens and workers. Many genes recorded as up-regulated in prospective workers appear to be unique to A. mellifera, suggesting that the workers' developmental pathway involves the participation of novel genes. Workers up-regulate more developmental genes than queens, whereas queens up-regulate a greater proportion of physiometabolic genes, including genes coding for metabolic enzymes and genes whose products are known to regulate the rate of mass-transforming processes and the general growth of the organism (e.g., tor). Many DEGs are likely to be involved in processes favoring the development of caste-biased structures, like brain, legs and ovaries, as well as genes that code for cytoskeleton constituents. Treatment of developing worker larvae with juvenile hormone (JH) revealed 52 JH responsive genes, specifically during the critical period of caste development. Using Gibbs sampling and Expectation Maximization algorithms, we discovered eight overrepresented cis-elements from four gene groups. Graph theory and complex networks concepts were adopted to attain powerful graphical representations of the interrelation between cis-elements and genes and objectively quantify the degree of relationship between these entities.

Conclusion

We suggest that clusters of functionally related DEGs are co-regulated during caste development in honeybees. This network of interactions is activated by nutrition-driven stimuli in early larval stages. Our data are consistent with the hypothesis that JH is a key component of the developmental determination of queen-like characters. Finally, we propose a conceptual model of caste differentiation in A. mellifera based on gene-regulatory networks.

Differential gene expression between adult queens and workers in the ant Lasius niger

Ants and other social insects forming large societies are generally characterized by marked reproductive division of labour. Queens largely monopolize reproduction whereas workers have little reproductive potential. In addition, some social insect species show tremendous lifespan differences between the queen and worker caste. Remarkably, queens and workers are usually genotypically identical, meaning that any phenotypic differences between the two castes arise from caste-specific gene expression. Using a combination of differential display, microarrays and reverse Northern blots, we found 16 genes that were differentially expressed between adult queens and workers in the ant Lasius niger, a species with highly pronounced reproductive division of labour and a several-fold lifespan difference between queens and workers. RNA ligase mediated rapid amplification of cDNA ends (RLM-RACE) and gene walking were used to further characterize these genes. On the basis of the molecular function of their nearest homologues, three genes appear to be involved in reproductive division of labour. Another three genes, which were exclusively overexpressed in queens, are possibly involved in the maintenance and repair of the soma, a candidate mechanism for lifespan determination. In-depth functional analyses of these genes are now needed to reveal their exact role.

Gene expression and the evolution of phenotypic diversity in social wasps

Background

Organisms are capable of developing different phenotypes by altering the genes they express. This phenotypic plasticity provides a means for species to respond effectively to environmental conditions. One of the most dramatic examples of phenotypic plasticity occurs in the highly social hymenopteran insects (ants, social bees, and social wasps), where distinct castes and sexes all arise from the same genes. To elucidate how variation in patterns of gene expression affects phenotypic variation, we conducted a study to simultaneously address the influence of developmental stage, sex, and caste on patterns of gene expression in Vespula wasps. Furthermore, we compared the patterns found in this species to those found in other taxa in order to investigate how variation in gene expression leads to phenotypic evolution.

Results

We constructed 11 different cDNA libraries derived from various developmental stages and castes of Vespula squamosa. Comparisons of overall expression patterns indicated that gene-expression differences distinguishing developmental stages were greater than expression differences differentiating sex or caste. Furthermore, we determined that certain sets of genes showed similar patterns of expression in the same phenotypic forms of different species. Specifically, larvae upregulated genes related to metabolism and genes possessing structural activity. Surprisingly, our data indicated that at least a few specific gene functions and at least one specific gene family are important components of caste differentiation across social insect taxa.

Conclusion

Despite research on various aspects of development originating from model systems, growth in understanding how development is related to phenotypic diversity relies on a growing literature of contrasting studies in non-model systems. In this study, we found that comparisons of patterns of gene expression with model systems highlighted areas of conserved and convergent developmental evolution across diverse taxa. Indeed, conserved biological functions across species implicated key functions related to how phenotypes are built. Finally, overall differences between social insect taxa suggest that the independent evolution of caste arose via distinct developmental trajectories.

Caste development and reproduction: a genome-wide analysis of hallmarks of insect eusociality

The honey bee queen and worker castes are a model system for developmental plasticity. We used established expressed sequence tag information for a Gene Ontology based annotation of genes that are differentially expressed during caste development. Metabolic regulation emerged as a major theme, with a caste-specific difference in the expression of oxidoreductases vs. hydrolases. Motif searches in upstream regions revealed group-specific motifs, providing an entry point to cis-regulatory network studies on caste genes. For genes putatively involved in reproduction, meiosis-associated factors came out as highly conserved, whereas some determinants of embryonic axes either do not have clear orthologs (bag of marbles, gurken, torso), or appear to be lacking (trunk) in the bee genome. Our results are the outcome of a first genome-based initiative to provide an annotated framework for trends in gene regulation during female caste differentiation (representing developmental plasticity) and reproduction.