Transcriptome dynamics during metamorphosis of imaginal discs into wings and thoracic dorsum in Apis mellifera castes.

BACKGROUND: Much of the complex anatomy of a holometabolous insect is built from disc-shaped epithelial structures found inside the larva, i.e., the imaginal discs, which undergo a rapid differentiation during metamorphosis. Imaginal discs-derived structures, like wings, are built through the action of genes under precise regulation. RESULTS: We analyzed 30 honeybee transcriptomes in the search for the gene expression needed for wings and thoracic dorsum construction from the larval wing discs primordia. Analyses were carried out before, during, and after the metamorphic molt and using worker and queen castes. Our RNA-seq libraries revealed 13,202 genes, representing 86.2% of the honeybee annotated genes. Gene Ontology analysis revealed functional terms that were caste-specific or shared by workers and queens. Genes expressed in wing discs and descendant structures showed differential expression profiles dynamics in premetamorphic, metamorphic and postmetamorphic developmental phases, and also between castes. At the metamorphic molt, when ecdysteroids peak, the wing buds of workers showed maximal gene upregulation comparatively to queens, thus underscoring differences in gene expression between castes at the height of the larval-pupal transition. Analysis of small RNA libraries of wing buds allowed us to build miRNA-mRNA interaction networks to predict the regulation of genes expressed during wing discs development. CONCLUSION: Together, these data reveal gene expression dynamics leading to wings and thoracic dorsum formation from the wing discs, besides highlighting caste-specific differences during wing discs metamorphosis.

PROTEINS OF THE INTEGUMENTARY SYSTEM OF THE HONEYBEE, Apis mellifera.

The integument of insects and other arthropods is composed of an inner basal lamina coated by the epidermis, which secretes the bulk of the outer integument layer, the cuticle. The genome sequencing of several insect species has allowed predicting classes of proteins integrating the cuticle. However, only a small proportion of them, as well as other proteins in the integumentary system, have been validated. Using two-dimensional gel electrophoresis coupled with mass spectrometry, we identified 45 different proteins in a total of 112 selected gel spots derived from thoracic integument samples of developing honeybee workers, including 14 cuticular proteins (AmelCPR 3, AmelCPR 12, AmelCPR 16, AmelCPR 27, apidermin 2, apidermin 3, endocuticle structural glycoprotein SgAbd-8-like, LOC100577363, LOC408365, LOC413679, LOC725454, LOC100576916, LOC725838, and peritrophin 3-C analogous). Gene ontology functional analysis revealed that the higher proportions of the identified proteins have molecular functions related to catalytic and structural molecule activities, are involved in metabolic biological processes, and pertain to the protein class of structural or cytoskeletal proteins and hydrolases. It is noteworthy that 26.7% of the identified proteins, including five cuticular proteins, were revealed as protein species resulting from allelic isoforms or derived from posttranslational modifications. Also, 66.7% of the identified cuticular proteins were expressed in more than one developmental phase, thus indicating that they are part of the larval, pupal, and adult cuticle. Our data provide experimental support for predicted honeybee gene products and new information on proteins expressed in the developing integument.

Genes involved in thoracic exoskeleton formation during the pupal-to-adult molt in a social insect model, Apis mellifera.

BACKGROUND: The insect exoskeleton provides shape, waterproofing, and locomotion via attached somatic muscles. The exoskeleton is renewed during molting, a process regulated by ecdysteroid hormones. The holometabolous pupa transforms into an adult during the imaginal molt, when the epidermis synthe3sizes the definitive exoskeleton that then differentiates progressively. An important issue in insect development concerns how the exoskeletal regions are constructed to provide their morphological, physiological and mechanical functions. We used whole-genome oligonucleotide microarrays to screen for genes involved in exoskeletal formation in the honeybee thoracic dorsum. Our analysis included three sampling times during the pupal-to-adult molt, i.e., before, during and after the ecdysteroid-induced apolysis that triggers synthesis of the adult exoskeleton. RESULTS: Gene ontology annotation based on orthologous relationships with Drosophila melanogaster genes placed the honeybee differentially expressed genes (DEGs) into distinct categories of Biological Process and Molecular Function, depending on developmental time, revealing the functional elements required for adult exoskeleton formation. Of the 1,253 unique DEGs, 547 were upregulated in the thoracic dorsum after apolysis, suggesting induction by the ecdysteroid pulse. The upregulated gene set included 20 of the 47 cuticular protein (CP) genes that were previously identified in the honeybee genome, and three novel putative CP genes that do not belong to a known CP family. In situ hybridization showed that two of the novel genes were abundantly expressed in the epidermis during adult exoskeleton formation, strongly implicating them as genuine CP genes. Conserved sequence motifs identified the CP genes as members of the CPR, Tweedle, Apidermin, CPF, CPLCP1 and Analogous-to-Peritrophins families. Furthermore, 28 of the 36 muscle-related DEGs were upregulated during the de novo formation of striated fibers attached to the exoskeleton. A search for cis-regulatory motifs in the 5'-untranslated region of the DEGs revealed potential binding sites for known transcription factors. Construction of a regulatory network showed that various upregulated CP- and muscle-related genes (15 and 21 genes, respectively) share common elements, suggesting co-regulation during thoracic exoskeleton formation. CONCLUSIONS: These findings help reveal molecular aspects of rigid thoracic exoskeleton formation during the ecdysteroid-coordinated pupal-to-adult molt in the honeybee.

Immunosenescence in honey bees (Apis mellifera L.) is caused by intrinsic senescence and behavioral physiology.

Young honey bee workers (0 to 2-3 weeks old) perform tasks inside the colony, including brood care (nursing), whereas older workers undergo foraging tasks during the next 3-4 weeks, when an intrinsic senescence program culminates in worker death. We hypothesized that foragers are less able to react to immune system stimulation than nurse bees and that this difference is due to an inefficient immune response in foragers. To test this hypothesis, we used an experimental design that allowed us to uncouple chronological age and behavior status (nursing/foraging). Worker bees from a normal age demography colony (where workers naturally transit from nursing to foraging tasks as they age) and of a single-cohort colony setup (composed of same-aged workers performing nursing or foraging tasks) were tested for survival and capability of activation of the immune system after bacterial injection. Expression of an antimicrobial peptide gene, defensin-1 (def-1), was used to assess immune system activation. We then checked whether the immune response includes changes in the expression of aging- and behavior-related genes, specifically vitellogenin (vg), juvenile hormone esterase (jhe), and insulin-like peptide-1 (ilp-1). We found a significant difference in survival rate between bees of different ages but carrying out the same tasks. Our results thus indicate that the bees' immune response is negatively affected by intrinsic senescence. Additionally, independent of age, foragers had a shorter lifespan than nurses after bacterial infection, although both were able to induce def-1 transcription. In the normal age demography colony, the immune system activation resulted in a reduction in the expression of vg, jhe and ilp-1 genes in foragers, but not in the nurse bees, demonstrating that age and behavior are both important influences on the bees' immune response. By disentangling the effects of age and behavior in the single-cohort colony, we found that vg, jhe and ilp-1 response to immune system stimulation was independent of behavior. Younger bees were able to mount a stronger immune response than older bees, thus highlighting age as an important factor for immunity. Taken together, our results provide new insights into how age and behavior affect the honey bee's immune response.

A honeybee storage protein gene, hex 70a, expressed in developing gonads and nutritionally regulated in adult fat body.

In preparing for metamorphosis, insect larvae store a huge amount of proteins in hemolymph, mainly hexamerins. Out of the four hexamerins present in the honeybee larvae, one, HEX 70a, exhibited a distinct developmental pattern, especially since it is also present in adults. Here, we report sequence data and experimental evidence suggesting alternative functions for HEX 70a, besides its well-known role as an amino acid resource during metamorphosis. The hex 70a gene consists of 6 exons and encodes a 684 amino acid chain containing the conserved hemocyanin N, M, and C domains. HEX 70a classifies as an arylphorin since it contains more than 15% of aromatic amino acids. In the fat body of adult workers, hex 70a expression turned out to be a nutrient-limited process. However, the fat body is not the only site for hex 70a expression. Both, transcript and protein subunits were also detected in developing gonads from workers, queens and drones, suggesting a role in ovary differentiation and testes maturation and functioning. In its putative reproductive role, HEX 70a however differs from the yolk protein, vitellogenin, since it was not detected in eggs or embryos.

Identification of a juvenile hormone esterase-like gene in the honey bee, Apis mellifera L.--expression analysis and functional assays.

Tight control over circulating juvenile hormone (JH) levels is of prime importance in an insect's life cycle. Consequently, enzymes involved in JH metabolism, especially juvenile hormone esterases (JHEs), play major roles during metamorphosis and reproduction. In the highly eusocial Hymenoptera, JH has been co-opted into additional functions, primarily in the development of the queen and worker castes and in age-related behavioral development of workers. Within a set of 21 carboxylesterases predicted in the honey bee genome we identified one gene (Amjhe-like) that contained the main functional motifs of insect JHEs. Its transcript levels during larval development showed a maximum at the switch from feeding to spinning behavior, coinciding with a JH titer minimum. In adult workers, the highest levels were observed in nurse bees, where a low JH titer is required to prevent the switch to foraging. Functional assays showed that Amjhe-like expression is induced by JH-III and suppressed by 20-hydroxyecdysone. RNAi-mediated silencing of Amjhe-like gene function resulted in a six-fold increase in the JH titer in adult worker bees. The temporal profile of Amjhe-like expression in larval and adult workers, the pattern of hormonal regulation and the knockdown phenotype are consistent with the function of this gene as an authentic JHE.

Methyl farnesoate epoxidase (mfe) gene expression and juvenile hormone titers in the life cycle of a highly eusocial stingless bee, Melipona scutellaris.

In social insects, juvenile hormone (JH) has acquired novel functions related to caste determination and division of labor among workers, and this is best evidenced in the honey bee. In contrast to honey bees, stingless bees are a much more diverse group of highly eusocial bees, and the genus Melipona has long called special attention due to a proposed genetic mechanism of caste determination. Here, we examined methyl farnesoate epoxidase (mfe) gene expression, encoding an enzyme relevant for the final step in JH biosynthesis, and measured the hemolymph JH titers for all life cycle stages of Melipona scutellaris queens and workers. We confirmed that mfe is exclusively expressed in the corpora allata. The JH titer is high in the second larval instar, drops in the third, and rises again as the larvae enter metamorphosis. During the pupal stage, mfe expression is initialy elevated, but then gradually drops to low levels before adult emergence. No variation was, however, seen in the JH titer. In adult virgin queens, mfe expression and the JH titer are significantly elevated, possibly associated with their reproductive potential. For workers we found that JH titers are lower in foragers than in nurse bees, while mfe expression did not differ. Stingless bees are, thus, distinct from honey bee workers, suggesting that they have maintained the ancestral gonadotropic function for JH. Hence, the physiological circuitries underlying a highly eusocial life style may be variable, even within a monophyletic clade such as the corbiculate bees.

The HEX 110 Hexamerin Is a Cytoplasmic and Nucleolar Protein in the Ovaries of Apis mellifera.

Hexamerins are insect storage proteins abundantly secreted by the larval fat body into the haemolymph. The canonical role of hexamerins consists of serving as an amino acid reserve for development toward the adult stage. However, in Apis mellifera, immunofluorescence assays coupled to confocal laser-scanning microscopy, and high-throughput sequencing, have recently shown the presence of hexamerins in other organs than the fat body. These findings have led us to study these proteins with the expectation of uncovering additional functions in insect development. We show here that a honeybee hexamerin, HEX 110, localizes in the cytoplasm and nucleus of ovarian cells. In the nucleus of somatic and germline cells, HEX 110 colocalized with a nucleolar protein, fibrillarin, suggesting a structural or even regulatory function in the nucleolus. RNase A provoked the loss of HEX 110 signals in the ovarioles, indicating that the subcellular localization depends on RNA. This was reinforced by incubating ovaries with pyronin Y, a RNA-specific dye. Together, the colocalization with fibrillarin and pyronin Y, and the sensitivity to RNase, highlight unprecedented roles for HEX110 in the nucleolus, the nuclear structure harbouring the gene cluster involved in ribosomal RNA production. However, the similar patterns of HEX 110 foci distribution in the active and inactive ovaries of queens and workers preclude its association with the functional status of these organs.

RNAi-Mediated Functional Analysis of Bursicon Genes Related to Adult Cuticle Formation and Tanning in the Honeybee, Apis mellifera.

Bursicon is a heterodimeric neurohormone that acts through a G protein-coupled receptor named rickets (rk), thus inducing an increase in cAMP and the activation of tyrosine hydroxylase, the rate-limiting enzyme in the cuticular tanning pathway. In insects, the role of bursicon in the post-ecdysial tanning of the adult cuticle and wing expansion is well characterized. Here we investigated the roles of the genes encoding the bursicon subunits during the adult cuticle development in the honeybee, Apis mellifera. RNAi-mediated knockdown of AmBurs α and AmBurs ß bursicon genes prevented the complete formation and tanning (melanization/sclerotization) of the adult cuticle. A thinner, much less tanned cuticle was produced, and ecdysis toward adult stage was impaired. Consistent with these results, the knockdown of bursicon transcripts also interfered in the expression of genes encoding its receptor, AmRk, structural cuticular proteins, and enzymes in the melanization/sclerotization pathway, thus evidencing roles for bursicon in adult cuticle formation and tanning. Moreover, the expression of AmBurs α, AmBurs ß and AmRk is contingent on the declining ecdysteroid titer that triggers the onset of adult cuticle synthesis and deposition. The search for transcripts of AmBurs α, AmBurs ß and candidate targets in RNA-seq libraries prepared with brains and integuments strengthened our data on transcript quantification through RT-qPCR. Together, our results support our premise that bursicon has roles in adult cuticle formation and tanning, and are in agreement with other recent studies pointing for roles during the pharate-adult stage, in addition to the classical post-ecdysial ones.

Molecular characterization of a cDNA encoding prophenoloxidase and its expression in Apis mellifera.

Phenoloxidase (PO), a melanin-synthesizing enzyme known to play an important role in insect defense, is found as a zymogen (ProPO) in hemolymph and cuticle, where it is activated by proteolysis. We characterized the first proPO cDNA in an eusocial insect, the Apis mellifera honey bee. The AmproPO cDNA contains an ORF of 2079 bp encoding 693 amino acids, and is composed of 9 exons and 8 introns. Southern blot of digested genomic DNA suggested that only one copy of the proPO gene is present in A. mellifera. The molecular mass of the deduced ProPO and the active enzyme was predicted to be 80.1 and 74.4 kDa, respectively. The calculated pI was 6.28. BLASTp search of the deduced amino acid sequence, and neighbor-joining analysis, showed similarity with ProPOs from other insects, ranging from 47% to 63%. Protein signature analyses revealed four conserved regions, including the two copper binding sites characteristic of arthropod ProPOs. RT-PCR and Southern blot showed the highest amount of AmproPO transcripts in workers whole body, followed by queens and drones. Expression was also detected in hemocytes and integument. Real time RT-PCR showed higher amounts of AmproPO transcripts in adults and older pupae than in younger pupae and larvae, suggesting a function of AmproPO in adult exoskeleton melanization and differentiation.

Effects of juvenile hormone and ecdysone on the timing of vitellogenin appearance in hemolymph of queen and worker pupae of Apis mellifera.

The caste-specific regulation of vitellogenin synthesis in the honeybee represents a problem with many yet unresolved details. We carried out experiments to determine when levels of vitellogenin are first detected in hemolymph of female castes of Apis mellifera, and whether juvenile hormone and ecdysteroids modulate this process. Vitellogenin levels were measured in hemolymph using immunological techniques. We show that in both castes the appearance of vitellogenin in the hemolymph occurs during the pupal period, but the timing was different in the queen and worker. Vitellogenin appears in queens during an early phase of cuticle pigmentation approximately 60h before eclosion, while in workers the appearance of vitellogenin is more delayed, initiating in the pharate adult stage, approximately 10h before eclosion. The timing of vitellogenin appearance in both castes coincides with a slight increase in endogenous levels of juvenile hormone that occurs at the end of pupal development. The correlation between these events was corroborated by topical application of juvenile hormone. Exogenous juvenile hormone advanced the timing of vitellogenin appearance in both castes, but caste-specific differences in timing were maintained. Injection of actinomycin D prevented the response to juvenile hormone. In contrast, queen and worker pupae that were treated with ecdysone showed a delay in the appearance of vitellogenin. These data suggest that queens and workers share a common control mechanism for the timing of vitellogenin synthesis, involving an increase in juvenile hormone titers in the presence of low levels of ecdysteroids.

Dimorphic ovary differentiation in honeybee (Apis mellifera) larvae involves caste-specific expression of homologs of ark and buffy cell death genes.

The establishment of the number of repeated structural units, the ovarioles, in the ovaries is one of the critical events that shape caste polyphenism in social insects. In early postembryonic development, honeybee (Apis mellifera) larvae have a pair of ovaries, each one consisting of almost two hundred ovariole primordia. While practically all these ovarioles continue developing in queen-destined larvae, they undergo massive programmed cell death (PCD) in worker-destined larvae. So as to gain insight into the molecular basis of this fundamental process in caste differentiation we used quantitative PCR (qPCR) and fluorescent in situ hybridization (FISH) to investigate the expression of the Amark and Ambuffy genes in the ovaries of the two honeybee castes throughout the fifth larval instar. These are the homologs of ark and buffy Drosophila melanogaster genes, respectively, involved in activating and inhibiting PCD. Caste-specific expression patterns were found during this time-window defining ovariole number. Amark transcript levels were increased when ovariole resorption was intensified in workers, but remained at low levels in queen ovaries. The transcripts were mainly localized at the apical end of all the worker ovarioles, but appeared in only a few queen ovarioles, thus strongly suggesting a function in mediating massive ovariolar cell death in worker larvae. Ambuffy was mainly expressed in the peritoneal sheath cells covering each ovariole. The levels of Ambuffy transcripts increased earlier in the developing ovaries of queens than in workers. Consistent with a protective role against cell death, Ambuffy transcripts were localized in practically all queen ovarioles, but only in few worker ovarioles. The results are indicative of a functional relationship between the expression of evolutionary conserved cell death genes and the morphological events leading to caste-specific ovary differentiation in a social insect.

Juvenile hormone biosynthesis gene expression in the corpora allata of honey bee (Apis mellifera L.) female castes.

Juvenile hormone (JH) controls key events in the honey bee life cycle, viz. caste development and age polyethism. We quantified transcript abundance of 24 genes involved in the JH biosynthetic pathway in the corpora allata-corpora cardiaca (CA-CC) complex. The expression of six of these genes showing relatively high transcript abundance was contrasted with CA size, hemolymph JH titer, as well as JH degradation rates and JH esterase (jhe) transcript levels. Gene expression did not match the contrasting JH titers in queen and worker fourth instar larvae, but jhe transcript abundance and JH degradation rates were significantly lower in queen larvae. Consequently, transcriptional control of JHE is of importance in regulating larval JH titers and caste development. In contrast, the same analyses applied to adult worker bees allowed us inferring that the high JH levels in foragers are due to increased JH synthesis. Upon RNAi-mediated silencing of the methyl farnesoate epoxidase gene (mfe) encoding the enzyme that catalyzes methyl farnesoate-to-JH conversion, the JH titer was decreased, thus corroborating that JH titer regulation in adult honey bees depends on this final JH biosynthesis step. The molecular pathway differences underlying JH titer regulation in larval caste development versus adult age polyethism lead us to propose that mfe and jhe genes be assayed when addressing questions on the role(s) of JH in social evolution.

Nuclear Immunolocalization of Hexamerins in the Fat Body of Metamorphosing Honey Bees.

Hexamerins are storage proteins with primordial functions in insect metamorphosis. They are actively secreted by the larval fat body and stored in the hemolymph. During metamorphosis, they return to the fat body to be processed. For decades, these proteins were thought to exclusively function as an amino acid source for tissue reconstruction during the non-feeding pupal and pharate adult stages and, in some species, for egg production. Recently, new findings have linked the hexamerins to caste polyphenism and gonad development in social insects. To explore the roles of hexamerins during the honey bee metamorphosis, we used specific antibodies in expression analysis by western blot, in situ immunolocalization by confocal laser-scanning microscopy and in vivo injections to lower their endogenous levels. Our expression analysis highlighted the changing expression patterns in the fat body and hemolymph during development, which is consistent with the temporal dynamics of hexamerin secretion, storage and depletion. Confocal microscopy showed hexamerin expression in the cytoplasm of both types of fat body cells, trophocytes and oenocytes. Notably, hexamerin foci were also found in the nuclei of these cells, thus confirming our western blot analysis of fat body nuclear-enriched fractions. We also observed that the decrease in soluble hexamerins in antibody-treated pharate adults led to a precocious adult ecdysis, perhaps in response to the lack (or decrease) in hexamerin-derived amino acids. Taken together, these findings indicate that hexamerins have other functions in addition to their well-established role as amino acid sources for development.

Potential costs of bacterial infection on storage protein gene expression and reproduction in queenless Apis mellifera worker bees on distinct dietary regimes.

Insects are able to combat infection by initiating an efficient immune response that involves synthesizing antimicrobial peptides and a range of other defense molecules. These responses may be costly to the organism, resulting in it exploiting endogenous resources to maintain homeostasis or support defense to the detriment of other physiological needs. We used queenless worker bees on distinct dietary regimes that may alter hemolymph protein storage and ovary activation to investigate the physiological costs of infection with Serratia marcescens. The expression of the genes encoding the storage proteins vitellogenin and hexamerin 70a, the vitellogenin receptor, and vasa (which has a putative role in reproduction), was impaired in the infected bees. This impairment was mainly evident in the bees fed beebread, which caused significantly higher expression of these genes than did royal jelly or syrup, and this was confirmed at the vitellogenin and hexamerin 70a protein levels. Beebread was also the only diet that promoted ovary activation in the queenless bees, but this activation was significantly impaired by the infection. The expression of the genes encoding the storage proteins apolipophorins-I and -III and the lipophorin receptor was not altered by infection regardless the diet provided to the bees. Similarly, the storage of apolipophorin-I in the hemolymph was only slightly impaired by the infection, independently of the supplied diet. Taken together these results indicate that, infection demands a physiological cost from the transcription of specific protein storage-related genes and from the reproductive capacity.

The juvenile hormone (JH) epoxide hydrolase gene in the honey bee (Apis mellifera) genome encodes a protein which has negligible participation in JH degradation.

Epoxide hydrolases are multifunctional enzymes that are best known in insects for their role in juvenile hormone (JH) degradation. Enzymes involved in JH catabolism can play major roles during metamorphosis and reproduction, such as the JH epoxide hydrolase (JHEH), which degrades JH through hydration of the epoxide moiety to form JH diol, and JH esterase (JHE), which hydrolyzes the methyl ester to produce JH acid. In the honey bee, JH has been co-opted for additional functions, mainly in caste differentiation and in age-related behavioral development of workers, where the activity of both enzymes could be important for JH titer regulation. Similarity searches for jheh candidate genes in the honey bee genome revealed a single Amjheh gene. Sequence analysis, quantification of Amjheh transcript levels and Western blot assays using an AmJHEH-specific antibody generated during this study revealed that the AmJHEH found in the fat body shares features with the microsomal JHEHs from several insect species. Using a partition assay we demonstrated that AmJHEH has a negligible role in JH degradation, which, in the honey bee, is thus performed primarily by JHE. High AmJHEH levels in larvae and adults were related to the ingestion of high loads of lipids, suggesting that AmJHEH has a role in dietary lipid catabolism.