Search

Made available in DSpace on 2022-04-29T08:35:40Z (GMT). No. of bitstreams: 0 Previous issue date: 2021-12-01 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. Departamento de Genética Faculdade de Medicina de Ribeirão Preto Universidade de São Paulo, Av. Bandeirantes 3900 Departamento de Tecnologia Faculdade de Ciências Agrárias e Veterinárias Universidade Estadual Paulista Júlio de Mesquita Filho Instituto Carlos Chagas FIOCRUZ Departamento de Biologia Faculdade de Filosofia Ciências e Letras de Ribeirão Preto Universidade de São Paulo, Av. Bandeirantes 3900 Departamento de Tecnologia Faculdade de Ciências Agrárias e Veterinárias Universidade Estadual Paulista Júlio de Mesquita Filho

Epitranscriptomics is an emerging field of investigation dedicated to the study of post-transcriptional RNA modifications. RNA methylations regulate RNA metabolism and processing, including changes in response to environmental cues. Although RNA modifications are conserved from bacteria to eukaryotes, there is little evidence of a epitranscriptomic pathway in insects. Here we identified genes related to RNA m6 A (N6-methyladenine) and m5 C (5-methylcytosine) methylation machinery in seven bee genomes (Apis mellifera, Melipona quadrifasciata, Frieseomelitta varia, Eufriesea mexicana, Bombus terrestris, Megachile rotundata and Dufourea novaeangliae). In A. mellifera, we validated the expression of methyltransferase genes and found that the global levels of m6 A and m5 C measured in the fat body and brain of adult workers differ significantly. Also, m6 A levels were differed significantly mainly between the fourth larval instar of queens and workers. Moreover, we found a conserved m5 C site in the honeybee 28S rRNA. Taken together, we confirm the existence of epitranscriptomic machinery acting in bees and open avenues for future investigations on RNA epigenetics in a wide spectrum of hymenopteran species. This article is protected by copyright. All rights reserved.

During the honeybee larval stage, queens develop larger brains than workers, with morphological differentiation appearing at the fourth larval phase (L4), just after a boost in nutritional difference both prospective females experience. The molecular promoters of this caste-specific brain development are already ongoing in previous larval phases. Transcriptomic analyses revealed a set of differentially expressed genes in the L3 brains of queens and workers, which represents the early molecular response to differential feeding females receive during larval development. Three genes of this set, hex70b, hex70c and hex110, are more highly transcribed in the brain of workers than in queens. The microRNAs miR-34, miR-210 and miR-317 are in higher levels in the queens' brain at the same phase of larval development. Here, we tested the hypothesis that the brain of workers expresses higher levels of hexamerins than that of queens during key phases of larval development and that this differential hexamerin genes expression is further enhanced by the repressing activity of miR-34, miR-210 and miR-317. Our transcriptional analyses showed that hex70b, hex70c and hex110 genes are differentially expressed in the brain of L3 and L4 larval phases of honeybee queens and workers. In silico reconstructed miRNA-mRNA interaction networks were validated using luciferase assays, which showed miR-34 and miR-210 negatively regulate hex70b and hex110 genes by directly and redundantly binding their 3'UTR (untranslated region) sequences. Taken together, our results suggest that miR-34 and miR-210 act together promoting differential brain development in honeybee castes by downregulating the expression of the putative antineurogenic hexamerin genes hex70b and hex110.

Background Meiotic recombination is a fundamental genetic process that shuffles allele combinations and promotes accurate segregation of chromosomes. Analyses of the ubiquitous variation of recombination rates within and across species suggest that recombination is evolving adaptively. All studied insects with advanced eusociality have shown exceptionally high recombination rates, which may represent a prominent case of adaptive evolution of recombination. However, our understanding of the relationship between social evolution and recombination rates is incomplete, partly due to lacking empirical data. Here, we present a linkage map of the monandrous, advanced eusocial Brazilian stingless bee, Frieseomelitta varia, providing the first recombination analysis in the diverse Meliponini (Hymenoptera, Apidae). Results Our linkage map includes 1417 markers in 19 linkage groups. This map spans approximately 2580 centimorgans, and comparisons to the physical genome assembly indicate that it covers more than 75 % of the 275 Megabasepairs (Mbp) F. varia genome. Thus, our study results in a genome-wide recombination rate estimate of 9.3–12.5 centimorgan per Mbp. This value is higher than estimates from nonsocial insects and comparable to other highly social species, although it does not support our prediction that monandry and strong queen-worker caste divergence of F. varia lead to even higher recombination rates than other advanced eusocial species. Conclusions Our study expands the association between elevated recombination and sociality in the order Hymenoptera and strengthens the support for the hypothesis that advanced social evolution in hymenopteran insects invariably selects for high genomic recombination rates.

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.

Nosema ceranae is a microsporidium that infects Apis mellifera, causing diverse physiological and behavioral alterations. Given the existence of individual and social mechanisms to reduce infection and fungal spread in the colony, bees may respond differently to infection depending on their rearing conditions. In this study, we investigated the effect of N. ceranae in honey bee foragers naturally infected with different fungal loads in a tropical region. In addition, we explored the effects of N. ceranae artificially infected young bees placed in a healthy colony under field conditions. Honey bees naturally infected with higher loads of N. ceranae showed downregulation of genes from Toll and IMD immune pathways and antimicrobial peptide (AMP) genes, but hemolymph total protein amount and Vitellogenin (Vg) titers were not affected. Artificially infected bees spread N. ceranae to the controls in the colony, but fungal loads were generally lower than those observed in cages, probably because of social immunity. Although no significant changes in mRNA levels of AMP-encoding were observed, N. ceranae artificially infected bees showed downregulation of miR-989 (an immune-related microRNA), lower vitellogenin gene expression, and decreased hemolymph Vg titers. Our results demonstrate for the first time that natural infection by N. ceranae suppresses the immune system of honey bee foragers in the field. This parasite is detrimental to the immune system of young and old bees, and disease spread, mitigation and containment will depend on the colony environment.

Made available in DSpace on 2020-12-10T20:02:24Z (GMT). No. of bitstreams: 0 Previous issue date: 2020-06-03 Brazilian National Science Fund (Conselho Nacional de Desenvolvimento Cientifico and Tecnologico) Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) BackgroundMost of our understanding on the social behavior and genomics of bees and other social insects is centered on the Western honey bee, Apis mellifera. The genus Apis, however, is a highly derived branch comprising less than a dozen species, four of which genomically characterized. In contrast, for the equally highly eusocial, yet taxonomically and biologically more diverse Meliponini, a full genome sequence was so far available for a single Melipona species only. We present here the genome sequence of Frieseomelitta varia, a stingless bee that has, as a peculiarity, a completely sterile worker caste.ResultsThe assembly of 243,974,526 high quality Illumina reads resulted in a predicted assembled genome size of 275Mb composed of 2173 scaffolds. A BUSCO analysis for the 10,526 predicted genes showed that these represent 96.6% of the expected hymenopteran orthologs. We also predicted 169,371 repetitive genomic components, 2083 putative transposable elements, and 1946 genes for non-coding RNAs, largely long non-coding RNAs. The mitochondrial genome comprises 15,144bp, encoding 13 proteins, 22 tRNAs and 2 rRNAs. We observed considerable rearrangement in the mitochondrial gene order compared to other bees. For an in-depth analysis of genes related to social biology, we manually checked the annotations for 533 automatically predicted gene models, including 127 genes related to reproductive processes, 104 to development, and 174 immunity-related genes. We also performed specific searches for genes containing transcription factor domains and genes related to neurogenesis and chemosensory communication.ConclusionsThe total genome size for F. varia is similar to the sequenced genomes of other bees. Using specific prediction methods, we identified a large number of repetitive genome components and long non-coding RNAs, which could provide the molecular basis for gene regulatory plasticity, including worker reproduction. The remarkable reshuffling in gene order in the mitochondrial genome suggests that stingless bees may be a hotspot for mtDNA evolution. Hence, while being just the second stingless bee genome sequenced, we expect that subsequent targeting of a selected set of species from this diverse clade of highly eusocial bees will reveal relevant evolutionary signals and trends related to eusociality in these important pollinators. Univ Sao Paulo, Fac Med Ribeirao Preto, Dept Genet, Ribeirao Preto, SP, Brazil Univ Fed Alfenas, Inst Ciencias Biomed, Dept Biol Celular & Desenvolvimento, Alfenas, MG, Brazil Univ Fed Vales Jequitinhonha & Mucuri, Fac Ciencias Biol & Saude, Dept Ciencias Biol, Diamantina, MG, Brazil Univ Fed Sao Carlos, Ctr Ciencias Biol & Saude, Dept Genet & Evolucao, Sao Carlos, SP, Brazil Univ Tecnol Fed Parana, Cornelio Procopio, PR, Brazil Univ Sao Paulo, Fac Med Ribeirao Preto, Dept Biol Celular & Mol & Bioagentes Patogen, Av Bandeirantes 3900, BR-14049900 Ribeirao Preto, SP, Brazil Univ, Estadual Paulista Julio de Mesquita Filho, Fac Ciencias Agr & Vet, Dept Tecnol, Jaboticabal, SP, Brazil Univ Sao Paulo, Fac Filosofia Ciencias & Letras Ribeirao Preto, Dept Biol, Ribeirao Preto, SP, Brazil Univ, Estadual Paulista Julio de Mesquita Filho, Fac Ciencias Agr & Vet, Dept Tecnol, Jaboticabal, SP, Brazil Brazilian National Science Fund (Conselho Nacional de Desenvolvimento Cientifico and Tecnologico): 454103/2014-0 FAPESP: 2015/06657-0 CAPES: 001

Made available in DSpace on 2021-06-25T10:13:53Z (GMT). No. of bitstreams: 0 Previous issue date: 2021-02-01 Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) Fundação de Amparo à Pesquisa do Estado de Minas Gerais (FAPEMIG) Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) Financiadora de Estudos e Projetos Apis mellifera adult workers feature more developed key brain regions than queens, which allows them to cope with the broad range of duties they need to perform in a colony. However, at the end of larval development, the brain of queens is largely more developed than that of workers. Major morphogenetic changes take place after metamorphosis that shift caste-specific brain development. Here, we tested the hypothesis that this phenomenon is hormonally governed and involves differential gene expression. Our molecular screening approach revealed a set of differentially expressed genes in Pp (first pharate-adult phase) brains between castes mainly coding for tissue remodelling and energy-converting proteins (e.g. hex 70a and ATPsynβ). An in-depth qPCR analysis of the transcriptional behaviour during pupal and pharate-adult developmental stage in both castes and in response to artificially augmented hormone titres of 18 genes/variants revealed that: i. subtle differences in hormone titres between castes might be responsible for the differential expression of the EcR and insulin/insulin-like signalling (IIS) pathway genes; ii. the morphogenetic activity of the IIS in brain development must be mediated by ILP-2, iii. which together with the tum, mnb and caspase system, can constitute the molecular effectors of the caste-specific opposing brain developmental trajectories. Departamento de Biologia Celular e do Desenvolvimento Instituto de Ciências Biomédicas Universidade Federal de Alfenas UNIFAL-MG Faculdade de Ciências Agrárias e Veterinárias UNESP - Universidade Estadual Paulista Departamento de Genética Faculdade de Medicina de Ribeirão Preto Universidade de São Paulo Departamento de Biologia Faculdade de Filosofia Ciências e Letras de Ribeirão Preto Universidade de São Paulo Faculdade de Ciências Agrárias e Veterinárias UNESP - Universidade Estadual Paulista FAPESP: 2005/03926-5 FAPESP: 2011/03171-5 CNPq: 406734/2018-6 FAPEMIG: APQ-02134-14 CAPES: Finance Code 001 Financiadora de Estudos e Projetos: MCTI/FINEP/CT-INFRA-PROINFRA 01/2008 Financiadora de Estudos e Projetos: MCTI/FINEP/CT-INFRA-PROINFRA 01/2011

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.

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.

Made available in DSpace on 2019-10-04T12:36:19Z (GMT). No. of bitstreams: 0 Previous issue date: 2019-03-14 Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) Differences in the timing of exoskeleton melanization and sclerotization are evident when comparing eusocial and solitary bees. This cuticular maturation heterochrony may be associated with life style, considering that eusocial bees remain protected inside the nest for many days after emergence, while the solitary bees immediately start outside activities. To address this issue, we characterized gene expression using large-scale RNA sequencing (RNA-seq), and quantified cuticular hydrocarbon (CHC) through gas chromatography-mass spectrometry in comparative studies of the integument (cuticle plus its underlying epidermis) of two eusocial and a solitary bee species. In addition, we used transmission electron microscopy (TEM) for studying the developing cuticle of these and other three bee species also differing in life style. We found 13,200, 55,209 and 30,161 transcript types in the integument of the eusocial Apis mellifera and Frieseomelitta varia, and the solitary Centris analis, respectively. In general, structural cuticle proteins and chitin-related genes were upregulated in pharate-adults and newly-emerged bees whereas transcripts for odorant binding proteins, cytochrome P450 and antioxidant proteins were overrepresented in foragers. Consistent with our hypothesis, a distance correlation analysis based on the differentially expressed genes suggested delayed cuticle maturation in A. mellifera in comparison to the solitary bee. However, this was not confirmed in the comparison with F. varia. The expression profiles of 27 of 119 genes displaying functional attributes related to cuticle formation/differentiation were positively correlated between A. mellifera and F. varia, and negatively or non-correlated with C. analis, suggesting roles in cuticular maturation heterochrony. However, we also found transcript profiles positively correlated between each one of the eusocial species and C. analis. Gene co-expression networks greatly differed between the bee species, but we identified common gene interactions exclusively between the eusocial species. Except for F. varia, the TEM analysis is consistent with cuticle development timing adapted to the social or solitary life style. In support to our hypothesis, the absolute quantities of n-alkanes and unsaturated CHCs were significantly higher in foragers than in the earlier developmental phases of the eusocial bees, but did not discriminate newly-emerged from foragers in C. analis. By highlighting differences in integument gene expression, cuticle ultrastructure, and CHC profiles between eusocial and solitary bees, our data provided insights into the process of heterochronic cuticle maturation associated to the way of life. Univ Sao Paulo, Dept Genet, Fac Med Ribeirao Preto, Ribeirao Preto, Brazil Hosp Clin Porto Alegre, Nucleo Bioinformat, Porto Alegre, RS, Brazil Univ Estadual Paulista, Fac Ciencias Agr & Vet, Dept Tecnol, Jaboticabal, Brazil Univ Sao Paulo, Dept Biol, Fac Filosofia Ciencias & Letras Ribeirao Preto, Ribeirao Preto, Brazil Univ Sao Paulo, Dept Fis & Quim, Fac Ciencias Farmaceut Ribeirao Preto, Ribeirao Preto, Brazil Univ Estadual Paulista, Fac Ciencias Agr & Vet, Dept Tecnol, Jaboticabal, Brazil FAPESP: 2014/13136-0 FAPESP: 2011/03171-5 FAPESP: 2014/50265-3 CAPES: 88887.160608/2017-00 CAPES: 001 : FAPESP 2012/24284-5

Like all other insects, two key signalling pathways [Toll and immune deficiency (Imd)] regulate the induction of honey bee immune effectors that target microbial pathogens. Amongst these effectors are antimicrobial peptides (AMPs) that are presumed to be produced by the nuclear factors kappa B (NF-κB) Dorsal and Relish from the Toll and Imd pathways, respectively. Using in silico analysis, we previously proposed that the honey bee AMP defensin-1 was regulated by the Toll pathway, whereas hymenoptaecin was regulated by Imd and abaecin by both the Toll and Imd pathways. Here we use an RNA interference (RNAi) assay to determine the role of Dorsal in regulating abaecin and defensin-1. Honey bees have two dorsal genes (dorsal-1 and dorsal-2) and two splicing isoforms of dorsal-1 (dorsal-1A and dorsal-1B). Accordingly, we used both single and multiple (double or triple) isoform knockdown strategies to clarify the roles of dorsal proteins and their isoforms. Down-regulation of defensin-1 was observed for dorsal-1A and dorsal-2 knockdowns, but abaecin expression was not affected by dorsal RNAi. We conclude that defensin-1 is regulated by Dorsal (Toll pathway).

Genetic exchange by hybridization or admixture can make an important contribution to evolution, and introgression of favourable alleles can facilitate adaptation to new environments. A small number of honeybees (Apis mellifera) with African ancestry were introduced to Brazil ~60 years ago, which dispersed and hybridized with existing managed populations of European origin, quickly spreading across much of the Americas in an example of a massive biological invasion. Here, we analyse whole-genome sequences of 32 Africanized honeybees sampled from throughout Brazil to study the effect of this process on genome diversity. By comparison with ancestral populations from Europe and Africa, we infer that these samples have 84% African ancestry, with the remainder from western European populations. However, this proportion varies across the genome and we identify signals of positive selection in regions with high European ancestry proportions. These observations are largely driven by one large gene-rich 1.4-Mbp segment on chromosome 11 where European haplotypes are present at a significantly elevated frequency and likely confer an adaptive advantage in the Africanized honeybee population. This region has previously been implicated in reproductive traits and foraging behaviour in worker bees. Finally, by analysing the distribution of ancestry tract lengths in the context of the known time of the admixture event, we are able to infer an average generation time of 2.0 years. Our analysis highlights the processes by which populations of mixed genetic ancestry form and adapt to new environments.

MicroRNAs (miRNAs) are key regulators of developmental processes, such as cell fate determination and differentiation. Previous studies showed Dicer knockdown in honeybee embryos disrupt the processing of functional mature miRNAs and impairs embryo patterning. Here we investigated the expression profiles of miRNAs in honeybee embryogenesis and the role of the highly conserved miR-34-5p in the regulation of genes involved in insect segmentation. A total of 221 miRNAs were expressed in honey bee embryogenesis among which 97 mature miRNA sequences have not been observed before. Interestingly, we observed a switch in dominance between the 5-prime and 3-prime arm of some miRNAs in different embryonic stages; however, most miRNAs present one dominant arm across all stages of embryogenesis. Our genome-wide analysis of putative miRNA-target networks and functional pathways indicates miR-34-5p is one of the most conserved and connected miRNAs associated with the regulation of genes involved in embryonic patterning and development. In addition, we experimentally validated that miR-34-5p directly interacts to regulatory elements in the 3′-untranslated regions of pair-rule (even-skipped, hairy, fushi-tarazu transcription factor 1) and cytoskeleton (actin5C) genes. Our study suggests that miR-34-5p may regulate the expression of pair-rule and cytoskeleton genes during early development and control insect segmentation.

The honeybee Apis mellifera has major ecological and economic importance. We analyze patterns of genetic variation at 8.3 million SNPs, identified by sequencing 140 honeybee genomes from a worldwide sample of 14 populations at a combined total depth of 634×. These data provide insight into the evolutionary history and genetic basis of local adaptation in this species. We find evidence that population sizes have fluctuated greatly, mirroring historical fluctuations in climate, although contemporary populations have high genetic diversity, indicating the absence of domestication bottlenecks. Levels of genetic variation are strongly shaped by natural selection and are highly correlated with patterns of gene expression and DNA methylation. We identify genomic signatures of local adaptation, which are enriched in genes expressed in workers and in immune system- and sperm motility-related genes that might underlie geographic variation in reproduction, dispersal and disease resistance. This study provides a framework for future investigations into responses to pathogens and climate change in honeybees.

Melipona subnitida, a tropical stingless bee, is an endemic species of the Brazilian northeast and exhibits great potential for honey and pollen production in addition to its role as one of the main pollinators of the Caatinga biome. To understand the genetic structure and better assist in the conservation of this species, we characterized the population variability of M. subnitida using geometric morphometrics of the forewing and cytochrome c oxidase I gene fragment sequencing. We collected workers from six localities in the northernmost distribution. Both methodologies indicated that the variability among the sampled populations is related both to the environment in which samples were collected and the geographical distance between the sampling sites, indicating that differentiation among the populations is due to the existence of at least evolutionary lineages. Molecular clock data suggest that this differentiation may have begun in the middle Pleistocene, approximately 396 kya. The conservation of all evolutionary lineages is important since they can present differential resistance to environmental changes, as resistance to drought and diseases.

Vitellogenin (Vg) and lipophorin (Lp) are lipoproteins which play important roles in female reproductive physiology of insects. Both are actively taken up by growing oocytes and especially Vg and its receptor are considered as female-specifically expressed. The finding that the fat body of in honey bee (Apis mellifera) drones synthesizes Vg and is present in hemolymph has long been viewed as a curiosity. The recent paradigm change concerning the role played by Vg in honey bee life history, especially social division of labor, has now led us to investigate whether a physiological constellation similar to that seen in female reproduction may also be represented in the male sex. By means of Western blot analysis we could show that both Vg and Lp are present in the reproductive tract of adult drones, including the accessory (mucus) glands, but apparently are not secreted. Furthermore, we analyzed the transcript levels of the genes encoding these proteins (vg and lp), as well as their putative receptors (Amvgr and Amlpr) in fat body and accessory glands. Whereas lp, vg and Amlpr transcript levels decreased with age in both tissues, Amvgr mRNA levels increased with age in fat body. To our knowledge this is the first report that vitellogenin and its receptor are co-expressed in the reproductive system of a male insect. We interpret these findings as a cross-sexual transfer of a social physiological trait, associated with the rewiring of the juvenile hormone/vitellogenin circuitry that occurred in the female sex of honey bees.

In insects, exoskeleton (cuticle) formation at each molt cycle includes complex biochemical pathways wherein the laccase enzymes (EC 1.10.3.2) may have a key role. We identified an Amlac2 gene that encodes a laccase2 in the honey bee, Apis mellifera, and investigated its function in exoskeleton differentiation. The Amlac2 gene consists of nine exons resulting in an ORF of 2193 nucleotides. The deduced translation product is a 731 amino acid protein of 81.5 kDa and a pI of 6.05. Amlac2 is highly expressed in the integument of pharate adults, and the expression precedes the onset of cuticle pigmentation and the intensification of sclerotization. In accordance with the temporal sequence of exoskeleton differentiation from anterior to posterior direction, the levels of Amlac2 transcript increase earlier in the thoracic than in the abdominal integument. The gene expression lasts even after the bees emerge from brood cells and begin activities in the nest, but declines after the transition to foraging stage, suggesting that maturation of the exoskeleton is completed at this stage. Post-transcriptional knockdown of Amlac2 gene expression resulted in structural abnormalities in the exoskeleton and drastically affected adult eclosion. By setting a ligature between the thorax and abdomen of early pupae we could delay the increase in hemolymph ecdysteroid levels in the abdomen. This severely impaired the increase in Amlac2 transcript levels and also the differentiation of the abdominal exoskeleton. Taken together, these results indicate that Amlac2 expression is controlled by ecdysteroids and has a critical role in the differentiation of the adult exoskeleton of honey bees.

In the Apis mellifera post-genomic era, RNAi protocols have been used in functional approaches. However, sample manipulation and invasive methods such as injection of double-stranded RNA (dsRNA) can compromise physiology and survival. To circumvent these problems, we developed a non-invasive method for honeybee gene knockdown, using a well-established vitellogenin RNAi system as a model. Second instar larvae received dsRNA for vitellogenin (dsVg-RNA) in their natural diet. For exogenous control, larvae received dsRNA for GFP (dsGFP-RNA). Untreated larvae formed another control group. Around 60% of the treated larvae naturally developed until adult emergence when 0.5 microg of dsVg-RNA or dsGFP-RNA was offered while no larvae that received 3.0 microg of dsRNA reached pupal stages. Diet dilution did not affect the removal rates. Viability depends not only on the delivered doses but also on the internal conditions of colonies. The weight of treated and untreated groups showed no statistical differences. This showed that RNAi ingestion did not elicit drastic collateral effects. Approximately 90% of vitellogenin transcripts from 7-day-old workers were silenced compared to controls. A large number of samples are handled in a relatively short time and smaller quantities of RNAi molecules are used compared to invasive methods. These advantages culminate in a versatile and a cost-effective approach.

Though the replacement of European bees by Africanized honey bees in tropical America has attracted considerable attention, little is known about the temporal changes in morphological and genetic characteristics in these bee populations. We examined the changes in the morphometric and genetic profiles of an Africanized honey bee population collected near where the original African swarms escaped, after 34 years of Africanization. Workers from colonies sampled in 1968 and in 2002 were morphometrically analyzed using relative warps analysis and an Automatic Bee Identification System (ABIS). All the colonies had their mitochondrial DNA identified. The subspecies that mixed to form the Africanized honey bees were used as a comparison for the morphometric analysis. The two morphometric approaches showed great similarity of Africanized bees with the African subspecies, Apis mellifera scutellata, corroborating with other markers. We also found the population of 1968 to have the pattern of wing venation to be more similar to A. m. scutellata than the current population. The mitochondrial DNA of European origin, which was very common in the 1968 population, was not found in the current population, indicating selective pressure replacing the European with the African genome in this tropical region. Both morphometric methodologies were very effective in discriminating the A. mellifera groups; the non-linear analysis of ABIS was the most successful in identifying the bees, with more than 94% correct classifications.

The switch from within-hive activities to foraging behavior is a major transition in the life cycle of a honeybee (Apis mellifera) worker. A prominent regulatory role in this switch has long been attributed to juvenile hormone (JH), but recent evidence also points to the yolk precursor protein vitellogenin as a major player in behavioral development. In the present study, we injected vitellogenin double-stranded RNA (dsVg) into newly emerged worker bees of Africanized genetic origin and introduced them together with controls into observation hives to record flight behavior. RNA interference-mediated silencing of vitellogenin gene function shifted the onset of long-duration flights (>10 min) to earlier in life (by 3–4 days) when compared with sham and untreated control bees. In fact, dsVg bees were observed conducting such flights extremely precociously, when only 3 days old. Short-duration flights (

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.

A cDNA encoding a cuticle protein containing the R&R Consensus was characterized in the honeybee integument. AmelCPR14 developmental expression is distinguished by an on-off-on pattern, the transition from a low to a high level of transcripts occurring as the ecdysteroid titer is declining after the peak that triggers the onset of pharate (pupal and adult) development. The transcript is abundant during cuticle tanning and sclerotization, and persists even in the adult integument, suggesting that the corresponding protein is required for differentiation and maintenance of the adult cuticle. Such developmental pattern suggested that AmelCPR14 gene might be regulated by the titer of ecdysteroids. We confirmed this hypothesis using different experimental strategies. By tying a ligature in early pupae to prevent exposure of abdominal integument to a high ecdysteroid titer, we delayed the accumulation of AmelCPR14 transcripts in the abdominal integument. This is consistent with ecdysteroid priming being required in pupae for the increase in AmelCPR14 expression in pharate adults. By injecting 20-hydroxyecdysone (20E) in early pupae we demonstrated that hormone titer decay after the peak is critical for AmelCPR14 expression induction. Exposure of pupal integument in vitro to a 20E concentration mimicking the pupal ecdysteroid peak repressed AmelCPR14 expression, which was recovered by hormone removal. Taken together, these data are consistent with an ecdysteroid pulse (increase in hormone titer followed by its decline) being critical for a high AmelCPR14 gene expression in pharate adults.

In honeybees, the haplodiploid sex determination system promotes a unique embryogenesis process wherein females develop from fertilized eggs and males develop from unfertilized eggs. However, the developmental strategies of honeybees during early embryogenesis are virtually unknown. Similar to most animals, the honeybee oocytes are supplied with proteins and regulatory elements that support early embryogenesis. As the embryo develops, the zygotic genome is activated and zygotic products gradually replace the preloaded maternal material. The analysis of small RNA and mRNA libraries of mature oocytes and embryos originated from fertilized and unfertilized eggs has allowed us to explore the gene expression dynamics in the first steps of development and during the maternal-to-zygotic transition (MZT). We localized a short sequence motif identified as TAGteam motif and hypothesized to play a similar role in honeybees as in fruit flies, which includes the timing of early zygotic expression (MZT), a function sustained by the presence of the zelda ortholog, which is the main regulator of genome activation. Predicted microRNA (miRNA)-target interactions indicated that there were specific regulators of haploid and diploid embryonic development and an overlap of maternal and zygotic gene expression during the early steps of embryogenesis. Although a number of functions are highly conserved during the early steps of honeybee embryogenesis, the results showed that zygotic genome activation occurs earlier in honeybees than in Drosophila based on the presence of three primary miRNAs (pri-miRNAs) (ame-mir-375, ame-mir-34 and ame-mir-263b) during the cleavage stage in haploid and diploid embryonic development.

An N-terminal amino acid sequence of a previously reported honey bee hexamerin, HEX 110 [Danty et al., Insect Biochem Mol Biol 28:387-397 (1998)], was used as reference to identify the predicted genomic sequence in a public GenBank database. In silico analysis revealed an ORF of 3,033 nucleotides that encompasses eight exons. The conceptual translation product is a glutamine-rich polypeptide with a predicted molecular mass of 112.2 kDa and pI of 6.43, which contains the conserved M and C hemocyanin domains. Semiquantitative and quantitative RT-PCR with specific primers allowed for an analysis of mRNA levels during worker bee development and under different physiological conditions. Concomitantly, the abundance of the respective polypeptide in the hemolymph was examined by SDS-PAGE. Hex 110 transcripts were found in high levels during the larval stages, then decreased gradually during the pupal stage, and increased again in adults. HEX 110 subunits were highly abundant in larval hemolymph, decreased at the spinning-stage, and remained at low levels in pupae and adults. In 5th instar larvae, neither starvation nor supplementation of larval food with royal jelly changed the Hex 110 transcript levels or the amounts of HEX 110 subunit in hemolymph. In adult workers, high levels of Hex 110 mRNA, but not of the respective subunit, were related to ovary activation, and also to the consumption of a pollen-rich diet.

The evolution of eusociality is one of the major transitions in evolution, but the underlying genomic changes are unknown. We compared the genomes of 10 bee species that vary in social complexity, representing multiple independent transitions in social evolution, and report three major findings. First, many important genes show evidence of neutral evolution as a consequence of relaxed selection with increasing social complexity. Second, there is no single road map to eusociality; independent evolutionary transitions in sociality have independent genetic underpinnings. Third, though clearly independent in detail, these transitions do have similar general features, including an increase in constrained protein evolution accompanied by increases in the potential for gene regulation and decreases in diversity and abundance of transposable elements. Eusociality may arise through different mechanisms each time, but would likely always involve an increase in the complexity of gene networks.

In the honeybee, Apis mellifera, vitellogenin (Vg) expression has been detected in the ovary of queens, but not in that of workers. In addition, larvae of both sexes produce Vg in significant amounts, which suggest that Vg serves for functions additional to oocyte growth and energy supply to the embryo. In vivo hormone treatment experiments suggest that the decrease of 20-hydroxyecdysone concentration occurring in previtellogenic phases allows Vg production. Southern analysis indicates that the Vg gene is present as a single copy in the honeybee genome. Arch. Insect Biochem. Physiol. 59:211–218, 2005. © 2005 Wiley-Liss, Inc.

Phenoloxidase (monophenol, l-dopa: oxygen oxidoreductase, EC 1.14.18.1) is a multicopper oxidase, which plays an important role in melanin synthesis, necessary for defense against intruding microorganisms and parasites, wound healing and cuticle pigmentation. A phenoloxidase from the hemolymph of honey bee pupae exhibited an apparent molecular mass of 70 kDa, as estimated by gel filtration and SDS-PAGE. Optimal pH and temperature were 6.5 and 20 degrees C, respectively. Activity was fully stable for 30 min at 50 degrees C. Like phenoloxidases from the hemolymph of other insects, the honey bee enzyme was activated by trypsin and inhibited by protease inhibitors and phenylthiourea. Only high concentrations of sodium azide effectively inhibited the detected activity. A low concentration (5 microM) of Ca2+, Mg2+, and Mn2+ had a stimulatory effect on the activity. Single Michaelis-Menten curves were observed for l-dopa and dopamine oxidation, but the affinity of the enzyme for dopamine was greater than for L-dopa. Semiquantitative RT-PCR and Southern blot analysis using a 359 bp labeled probe, and quantification of the prophenoloxidase mRNA levels by real-time PCR showed increased amounts of transcripts in hemocytes and integument from young pupae injected with 20-hydroxyecdysone.

Social life is prone to invasion by microorganisms, and binding of ferric ions by transferrin is an efficient strategy to restrict their access to iron. In this study, we isolated cDNA and genomic clones encoding an Apis mellifera transferrin (AmTRF) gene. It has an open reading frame (ORF) of 2136 bp spread over nine exons. The deduced protein sequence comprises 686 amino acid residues plus a 26 residues signal sequence, giving a predicted molecular mass of 76 kDa. Comparison of the deduced AmTRF amino acid sequence with known insect transferrins revealed significant similarity extending over the entire sequence. It clusters with monoferric transferrins, with which it shares putative iron-binding residues in the N-terminal lobe. In a functional analysis of AmTRF expression in honey bee development, we monitored its expression profile in the larval and pupal stages. The negative regulation of AmTRF by ecdysteroids deduced from the developmental expression profile was confirmed by experimental treatment of spinning-stage honey bee larvae with 20-hydroxyecdysone, and of fourth instar-larvae with juvenile hormone. A juvenile hormone application to spinning-stage larvae, in contrast, had only a minor effect on AmTRF transcript levels. This is the first study implicating ecdysteroids in the developmental regulation of transferrin expression in an insect species.

A striking example of plasticity in life span is seen in social insects such as ants and bees, where different castes may display distinct ageing patterns. In particular, the honeybee offers an intriguing illustration of environmental control on ageing rate. Honeybee workers display a temporal division of labour where young bees (or 'hive bees') perform tasks within the brood nest, and older bees forage for nectar, pollen propolis and water. When bees switch from the hive bee to the forager stage, their cellular defence machinery is down-regulated by a dramatic reduction in the number of functioning haemocytes (immunocytes). This study documents that the yolk precursor vitellogenin is likely to be involved in a regulatory pathway that controls the observed decline in somatic maintenance function of honeybee foragers. An association between the glyco-lipoprotein vitellogenin and immune function has not previously been reported for any organism. Honeybee workers are functionally sterile, and via the expression of juvenile hormone, a key gonotrophic hormone in adult insects, their vitellogenin levels are influenced by social interactions with other bees. Our results therefore suggest that in terms of maintenance of the cellular immune system, senescence of the honeybee worker is under social control.

The cDNA of Apis mellifera vitellogenin was cloned and sequenced. It is 5440 bp long and contains an ORF of 1770 amino acids (including a putative signal peptide of 16 residues). The deduced amino acid sequence shows significant similarity with other hymenopteran vitellogenins (58% with Pimpla nipponica and 54% with Athalia rosae). The alignment with 19 insect vitellogenins shows a high number of conserved motifs; for example, close to the C-terminus there is a GL/ICG motif followed by nine cysteines, as occurs in all hymenopteran species, and, as in other insect vitellogenins, a DGXR motif is located 18 residues upstream the GL/ICG motif. Phylogenetic analysis of vitellogenin sequences available in insects gave a tree that is congruent with the currently accepted insect phylogenetic schemes. Using two fragments of the vitellogenin cDNA as probes, we analyzed by Northern blot the sex- and caste-specific patterns of vitellogenin expression in pupae and adults of A. mellifera. In queens, vitellogenin mRNA was first detected in mid-late pupal stage, whereas in workers it was first detected in late pupal stage. Vitellogenin mRNA was also observed in drones, although it was first detected not in pupae but in freshly molted adults.

Evidence from field wasps and bumblebees appoints the endocrine system as a mediator between dominance status and ovarian activity in primitively social Hymenoptera. In this comparative study on ecdysteroid titers in the highly social honey bee, Apis mellifera, and a stingless bee, Melipona quadrifasciata, we focussed on the relationship between the ecdysteroid titer, social conditions (presence or absence of the queen), and ovary activity. In contrast to bumblebees, ecdysteroid titers in honey bee and stingless bee workers were either not altered, or dropped to even lower levels after the queen was removed. We also did not detect differences between virgin queens and mated, egg laying queens. These results suggest that ecdysteroids may have lost most of their reproductive functions - yet gained functions in larval caste differentiation - as higher levels of social organization were attained in the evolution of social insects. The observation that ecdysteroid titers are transiently elevated in young workers adds a new, yet functionally still speculative facet to hormonal regulation in insect societies.

Increasing evidence suggests small non-coding RNAs (ncRNAs) such as microRNAs (miRNAs) control levels of mRNA expression during experience-related remodelling of the brain. Here we use an associative olfactory learning paradigm in the honeybee Apis mellifera to examine gene expression changes in the brain during memory formation. Brain transcriptome analysis reveals a general downregulation of protein-coding genes, including asparagine synthetase and actin, and upregulation of ncRNAs. miRNA-mRNA network predictions together with PCR validation suggest miRNAs including miR-210 and miR-932 target the downregulated protein-coding genes. Feeding cholesterol-conjugated antisense RNA to bees results in the inhibition of miR-210 and of miR-932. Loss of miR-932 impairs long-term memory formation, but not memory acquisition. Functional analyses show that miR-932 interacts with Act5C, providing evidence for direct regulation of actin expression by an miRNA. An activity-dependent increase in miR-932 expression may therefore control actin-related plasticity mechanisms and affect memory formation in the brain. published

We developed a rapid method for extraction of DNA from honey bees, Apis mellifera, and from the parasitic bee mite, Varroa destructor. The advantages include fast processing and low toxicity of the substances that are utilized. We used lysis buffer with nonionic detergents to lyse cell walls and proteinase K to digest proteins. We tested whole thorax, thoracic muscle mass, legs, and antennae from individual bees; the mites were processed whole (1 mite/sample). Each thorax was incubated whole, without cutting, because exocuticle color pigment darkened the extraction solution, interfering with PCR results. The procedure was performed with autoclaved equipment and laboratory gloves. For each sample, we used 100 µL lysis buffer (2 mL stock solution of 0.5 M Tris/HCl, pH 8.5, 10 mL stock solution of 2 M KCl, 500 µL solution of 1 M MgCl2, 2 mL NP40, and 27.6 g sucrose, completed to 200 mL with bidistilled water and autoclaved) and 2 µL proteinase K (10 mg/mL in bidistilled water previously autoclaved, as proteinase K cannot be autoclaved). Tissues were incubated in the solutions for 1-2 h in a water bath (62°-68 °C) or overnight at 37 °C. After incubation, the tissues were removed from the extraction solution (lysis buffer + proteinase K) and the solution heated to 92 °C for 10 min, for proteinase K inactivation. Then, the solution with the extracted DNA was stored in a refrigerator (4°-8 °C) or a freezer (-20 °C). This method does not require centrifugation or phenol/chloroform extraction. The reduced number of steps allowed us to sample many individuals/day. Whole mites and bee antennae were the most rapidly processed. All bee tissues gave the same quality DNA. This method, even using a single bee antenna or a single mite, was adequate for extraction and analysis of bee genomic and mitochondrial DNA and mite genomic DNA.

In honey bees, vitellogenin (Vg) is hypothesized to be a major factor affecting hormone signaling, food-related behavior, immunity, stress resistance and lifespan. MicroRNAs, which play important roles in post-transcriptional gene regulation, likewise affect many biological processes. The actions of microRNAs and Vg are known to intersect in the context of reproduction; however, the role of these associations on social behavior is unknown. The phenotypic effects of Vg knockdown are best established and studied in the forager stage of workers. Thus, we exploited the well-established RNA interference (RNAi) protocol for Vg knockdown to investigate its downstream effects on microRNA population in honey bee foragers' brain and fat body tissue. To identify microRNAs that are differentially expressed between tissues in control and knockdown foragers, we used μParaflo microfluidic oligonucleotide microRNA microarrays. Our results showed that 76 and 74 microRNAs were expressed in the brain of control and knockdown foragers whereas 66 and 69 microRNAs were expressed in the fat body of control and knockdown foragers, respectively. Target prediction identified potential seed matches for a differentially expressed subset of microRNAs affected by Vg knockdown. These candidate genes are involved in a broad range of biological processes including insulin signaling, juvenile hormone (JH) and ecdysteroid signaling previously shown to affect foraging behavior. Thus, here we demonstrate a causal link between the Vg knockdown forager phenotype and variation in the abundance of microRNAs in different tissues, with possible consequences for the regulation of foraging behavior. * dsRNA : double-stranded RNA GFBr : control forager brain GFFb : control forager fat body JH : juvenile hormone RNAi : RNA interference VFBr : knockdown forager brain VFFb : knockdown forager fat body Vg : vitellogenin

In insects, a rapid and massive synthesis of antimicrobial peptides (AMPs) is activated through signaling pathways (Toll and Imd) to combat invading microbial pathogens. However, it is still unclear whether different types of bacteria provoke specific responses. Immune response mechanisms and the activation of specific genes were investigated by challenging Apis mellifera workers with the Gram-negative bacterium Serratia marcescens or the Gram-positive bacterium Micrococcus luteus. The immune system responded by activating most genes of the Toll and Imd pathways, particularly AMP genes. However, genes specifically regulated by M. luteus or S. marcescens were not detected, suggesting an interaction between the signaling pathways that lead to immune effectors synthesis. Despite this finding, kappaB motifs in the 5'-UTRs of selected genes suggest a pathway-specific control of AMP and transferrin-1 gene expression. Regulation by miRNAs was also investigated and revealed a number of candidates for the post-transcriptional regulation of immune genes in bees.

The fat body (FB) consists of two types of cells: throphocytes and oenocytes. Throphocytes are related to intermediary metabolism storing lipids, carbohydrates, and proteins while oenocytes play role in the lipids and lipoproteins production. The vitellogenin is the precursor of egg yolk (vitelline) and is synthesized on FB. The aim of this work was to analyze the effects of hormones acting in bee reproduction, as juvenile hormone (JH) and ecdisteroids (20 HE) on FB cells, where vitellogenin is synthesized. For the study were chose nurse workers that in Melipona quadrifasciata anthidioides present activated ovaries and produce eggs, and virgin queens whose ovaries are not yet activated, presenting only previtellogenic follicles. FB trophocytes from these classes of bees were cultivated in media containing different amounts of JH and 20-HE. The effects on trophocytes cytoplasm reserves of lipids, proteins, and activity of acid phosphatase were compared by observing preparations from cultured FB, treated and control, by transmission electron microscopy (TEM). The results showed that the hormones effects are related to the bee's caste and functional ovary stage. The role of acid phosphatase on mobilization of the trophocyte reserves was also determined.

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.

Cuticle renewal is a complex biological process that depends on the cross talk between hormone levels and gene expression. This study characterized the expression of two genes encoding cuticle proteins sharing the four conserved amino acid blocks of the Tweedle family, AmelTwdl1 and AmelTwdl2, and a gene encoding a cuticle peroxidase containing the Animal haem peroxidase domain, Ampxd, in the honey bee. Gene sequencing and annotation validated the formerly predicted tweedle genes, and revealed a novel gene, Ampxd, in the honey bee genome. Expression of these genes was studied in the context of the ecdysteroid-coordinated pupal-to-adult molt, and in different tissues. Higher transcript levels were detected in the integument after the ecdysteroid peak that induces apolysis, coinciding with the synthesis and deposition of the adult exoskeleton and its early differentiation. The effect of this hormone was confirmed in vivo by tying a ligature between the thorax and abdomen of early pupae to prevent the abdominal integument from coming in contact with ecdysteroids released from the prothoracic gland. This procedure impaired the natural increase in transcript levels in the abdominal integument. Both tweedle genes were expressed at higher levels in the empty gut than in the thoracic integument and trachea of pharate adults. In contrast, Ampxd transcripts were found in higher levels in the thoracic integument and trachea than in the gut. Together, the data strongly suggest that these three genes play roles in ecdysteroid-dependent exoskeleton construction and differentiation and also point to a possible role for the two tweedle genes in the formation of the cuticle (peritrophic membrane) that internally lines the gut.

Insect hexamerins have long been known as storage proteins that are massively synthesized by the larval fat body and secreted into hemolymph. Following the larval-to-pupal molt, hexamerins are sequestered by the fat body via receptor-mediated endocytosis, broken up, and used as amino acid resources for metamorphosis. In the honey bee, the transcript and protein subunit of a hexamerin, HEX 70a, were also detected in ovaries and testes. Aiming to identify the subcellular localization of HEX 70a in the female and male gonads, we used a specific antibody in whole mount preparations of ovaries and testes for analysis by confocal laser-scanning microscopy. Intranuclear HEX 70a foci were evidenced in germ and somatic cells of ovarioles and testioles of pharate-adult workers and drones, suggesting a regulatory or structural role. Following injection of the thymidine analog EdU we observed co-labeling with HEX 70a in ovariole cell nuclei, inferring possible HEX 70a involvement in cell proliferation. Further support to this hypothesis came from an injection of anti-HEX 70a into newly ecdysed queen pupae where it had a negative effect on ovariole thickening. HEX 70a foci were also detected in ovarioles of egg laying queens, particularly in the nuclei of the highly polyploid nurse cells and in proliferating follicle cells. Additional roles for this storage protein are indicated by the detection of nuclear HEX 70a foci in post-meiotic spermatids and spermatozoa. Taken together, these results imply undescribed roles for HEX 70a in the developing gonads of the honey bee and raise the possibility that other hexamerins may also have tissue specific functions.

Background Hexamerins are hemocyanin-derived proteins that have lost the ability to bind copper ions and transport oxygen; instead, they became storage proteins. The current study aimed to broaden our knowledge on the hexamerin genes found in the honey bee genome by exploring their structural characteristics, expression profiles, evolution, and functions in the life cycle of workers, drones and queens. Results The hexamerin genes of the honey bee (hex 70a, hex 70b, hex 70c and hex 110) diverge considerably in structure, so that the overall amino acid identity shared among their deduced protein subunits varies from 30 to 42%. Bioinformatics search for motifs in the respective upstream control regions (UCRs) revealed six overrepresented motifs including a potential binding site for Ultraspiracle (Usp), a target of juvenile hormone (JH). The expression of these genes was induced by topical application of JH on worker larvae. The four genes are highly transcribed by the larval fat body, although with significant differences in transcript levels, but only hex 110 and hex 70a are re-induced in the adult fat body in a caste- and sex-specific fashion, workers showing the highest expression. Transcripts for hex 110, hex 70a and hex70b were detected in developing ovaries and testes, and hex 110 was highly transcribed in the ovaries of egg-laying queens. A phylogenetic analysis revealed that HEX 110 is located at the most basal position among the holometabola hexamerins, and like HEX 70a and HEX 70c, it shares potential orthology relationship with hexamerins from other hymenopteran species. Conclusions Striking differences were found in the structure and developmental expression of the four hexamerin genes in the honey bee. The presence of a potential binding site for Usp in the respective 5' UCRs, and the results of experiments on JH level manipulation in vivo support the hypothesis of regulation by JH. Transcript levels and patterns in the fat body and gonads suggest that, in addition to their primary role in supplying amino acids for metamorphosis, hexamerins serve as storage proteins for gonad development, egg production, and to support foraging activity. A phylogenetic analysis including the four deduced hexamerins and related proteins revealed a complex pattern of evolution, with independent radiation in insect orders.

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.

Proteins stored in insect hemolymph may serve as a source of amino acids and energy for metabolism and development. The expression of the main storage proteins was assessed in bacterial-challenged honey bees using real-time (RT)-PCR and Western blot. After ensuring that the immune system had been activated by measuring the ensuing expression of the innate immune response genes, defensin-1 (def-1) and prophenoloxidase (proPO), we verified the expression of four genes encoding storage proteins. The levels of vitellogenin (vg) mRNA and of the respective protein were significantly lowered in bees injected with bacteria or water only (injury). An equivalent response was observed in orally-infected bees. The levels of apolipophorin II/I (apoLp-II/I) and hexamerin (hex 70a) mRNAs did not significantly change, but levels of Hex 70a protein subunit showed a substantial decay after bacterial challenge or injury. Infection also caused a strong reduction in the levels of apoLp-III transcripts. Our findings are consistent with a down-regulation of the expression and accumulation of storage proteins as a consequence of activation of the immune system, suggesting that this phenomenon represents a strategy to redirect resources to combat injury or infection. © 2009 Wiley Periodicals, Inc.

Two members of the low density lipoprotein receptor (LDLR) family were identified as putative orthologs for a vitellogenin receptor (Amvgr) and a lipophorin receptor (Amlpr) in the Apis mellifera genome. Both receptor sequences have the structural motifs characteristic of LDLR family members and show a high degree of similarity with sequences of other insects. RT-PCR analysis of Amvgr and Amlpr expression detected the presence of both transcripts in different tissues of adult female (ovary, fat body, midgut, head and specifically hypopharyngeal gland), as well as in embryos. In the head RNA samples we found two variant forms of AmLpR: a full length one and a shorter one lacking 29 amino acids in the O-linked sugar domain. In ovaries the expression levels of the two honey bee LDLR members showed opposing trends: whereas Amvgr expression was upregulated as the ovaries became activated, Amlpr transcript levels gradually declined. In situ hybridization analysis performed on ovaries detected Amvgr mRNA exclusively in germ line cells and corroborated the qPCR results showing an increase in Amvgr gene expression concomitant with follicle growth.

Ecdysteroids regulate many aspects of insect physiology after binding to a heterodimer composed of the nuclear hormone receptor proteins ecdysone receptor (EcR) and ultraspiracle (Usp). Several lines of evidence have suggested that the latter also plays important roles in mediating the action of juvenile hormone (JH) and, thus, integrates signaling by the two morphogenetic hormones. By using an RNAi approach, we show here that Usp participates in the mechanism that regulates the progression of pupal development in Apis mellifera, as indicated by the observed pupal developmental delay in usp knocked-down bees. Knock-down experiments also suggest that the expression of regulatory genes such as ftz transcription factor 1 (ftz-f1) and juvenile hormone esterase (jhe) depend on Usp. Vitellogenin (vg), the gene coding the main yolk protein in honeybees, does not seem to be under Usp regulation, thus suggesting that the previously observed induction of vg expression by JH during the last stages of pupal development is mediated by yet unknown transcription factor complexes.

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.

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.

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.

A comparison of the most conserved sex-determining genes between the fruit fly, Drosophila melanogaster, and the honey bee, Apis mellifera, was performed with bioinformatics tools developed for computational molecular biology. An initial set of protein sequences already described in the fruit fly as participants of the sex-determining cascade was retrieved from the Gene Ontology database (http://www.geneontology.org/) and aligned against a database of protein sequences predicted from the honey bee genome. The doublesex (dsx) gene is considered one of the most conserved sex-determining genes among metazoans, and a male-specific partial cDNA of putative A. mellifera dsx gene (Amdsx) was identified experimentally. The theoretical predictions were developed in the context of sequence similarity. Experimental evidence indicates that dsx is present in embryos and larvae, and that it encodes a transcription factor widely conserved in metazoans, containing a DM DNA-binding domain implicated in the regulation of the expression of genes involved in sexual phenotype formation.

Ce travail résulte de la collaboration de très nombreux chercheurs. Seuls les auteurs de la rubrique Physical and Genetic Mapping sont cités explicitement.; Here we report the genome sequence of the honeybee Apis mellifera, a key model for social behaviour and essential to global ecology through pollination. Compared with other sequenced insect genomes, the A. mellifera genome has high A1T and CpG contents, lacks major transposon families, evolves more slowly, and is more similar to vertebrates for circadian rhythm, RNA interference and DNA methylation genes, among others. Furthermore, A.mellifera has fewer genes for innate immunity, detoxification enzymes, cuticle-forming proteins and gustatory receptors, more genes for odorant receptors, and novel genes for nectar and pollen utilization, consistent with its ecology and social organization. Compared to Drosophila, genes in early developmental pathways differ in Apis, whereas similarities exist for functions that differ markedly, such as sex determination, brain function and behaviour. Population genetics suggests a novel African origin for the species A.mellifera and insights into whether Africanized bees spread throughout the New World via hybridization or displacement