Transcriptional signature of host shift in the seed beetle Zabrotes subfasciatus.

In phytophagous insects, adaptation to a new host is a dynamic process, in which early and later steps may be underpinned by different features of the insect genome. Here, we tested the hypothesis that early steps of this process are underpinned by a shift in gene expression patterns. We set up a short-term artificial selection experiment (10 generations) for the use of an alternative host (Cicer arietinum) on populations of the bean beetle Zabrotes subfasciatus. Using Illumina sequencing on young adult females, we show the selected populations differ in the expression of genes associated to stimuli, signalling, and developmental processes. Particularly, the "C. arietinum" population shows upregulation of histone methylation genes, which may constitute a strategy for fine-tuning the insect global gene expression network. Using qPCR on body regions, we demonstrated that the "Phaseolus vulgaris" population upregulates the genes polygalacturonase and egalitarian and that the expression of an odorant receptor transcript variant changes over generations. Moreover, in this population we detected the existence of vitellogenin (Vg) variants in both males and females, possibly harbouring canonical reproductive function in females and extracellular unknown functions in males. This study provides the basis for future genomic investigations seeking to shed light on the nature of the proximate mechanisms involved in promoting differential gene expression associated to insect development and adaptation to new hosts.

Transcriptional expression of m6A and m5C RNA methyltransferase genes in the brain and fat body of honey bee adult workers.

Honey bee (Apis mellifera) adult workers change behaviors and nutrition according to age progression. Young workers, such as nurses, perform in-hive tasks and consume protein-rich pollen, while older workers (foragers) leave the colony to search for food, and consume carbohydrate-rich nectar. These environmentally stimulated events involve transcriptional and DNA epigenetic marks alterations in worker tissues. However, post-transcriptional RNA modifications (epitranscriptomics) are still poorly explored in bees. We investigated the transcriptional profiles of m6A and m5C RNA methyltransferases in the brain and fat body of adult workers of 1) different ages and performing different tasks [nurses of 8 days-old (N-8D) and foragers of 29 days-old (F-29D), sampled from wild-type colonies], and 2) same-aged young workers caged in an incubator and treated with a pollen-rich [PR] or a pollen-deprived [PD] diet for 8 days. In the brain, METTL3, DNMT2, NOP2, NSUN2, NSUN5, and NSUN7 genes increased expression during adulthood (from N-8D to F-29D), while the opposite pattern was observed in the fat body for METTL3, DNMT2, and NSUN2 genes. Regarding diet treatments, high expression levels were observed in the brains of the pollen-deprived group (DNMT2, NOP2, and NSUN2 genes) and the fat bodies of the pollen-rich group (NOP2, NSUN4, and NSUN5 genes) compared to the brains of the PR group and the fat bodies of the PD group, respectively. Our data indicate that RNA epigenetics may be an important regulatory layer in the development of adult workers, presenting tissue-specific signatures of RNA methyltransferases expression in response to age, behavior, and diet content.

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.

Recombination mapping of the Brazilian stingless bee Frieseomelitta varia confirms high recombination rates in social hymenoptera.

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.

SSD - a free software for designing multimeric mono-, bi- and trivalent shRNAs.

RNA interference (RNAi) is a powerful gene silencing technology, widely used in analyses of reverse genetics, development of therapeutic strategies and generation of biotechnological products. Here we present a free software tool for the rational design of RNAi effectors, named siRNA and shRNA designer (SSD). SSD incorporates our previously developed software Strand Analysis to construct template DNAs amenable for the large scale production of mono-, bi- and trivalent multimeric shRNAs, via in vitro rolling circle transcription. We tested SSD by creating a trivalent multimeric shRNA against the vitellogenin gene of Apis mellifera. RT-qPCR analysis revealed that our molecule promoted a decrease in more than 50% of the target mRNA, in a dose-dependent manner, when compared to the control group. Thus, SSD software allows the easy design of multimeric shRNAs, for single or multiple simultaneous knockdowns, which is especially interesting for studies involving large amounts of double-stranded molecules.

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.

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.

Genomewide analysis of admixture and adaptation in the Africanized honeybee.

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.

Transcriptome Analysis of Honeybee (Apis Mellifera) Haploid and Diploid Embryos Reveals Early Zygotic Transcription during Cleavage.

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.

Social evolution. Genomic signatures of evolutionary transitions from solitary to group living.

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.

A worldwide survey of genome sequence variation provides insight into the evolutionary history of the honeybee Apis mellifera.

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.

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.

Evidence of at least two evolutionary lineages in Melipona subnitida (Apidae, Meliponini) suggested by mtDNA variability and geometric morphometrics of forewings.

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.

Nutritionally driven differential gene expression leads to heterochronic brain development in honeybee castes.

The differential feeding regimes experienced by the queen and worker larvae of the honeybee Apis mellifera shape a complex endocrine response cascade that ultimately gives rise to differences in brain morphologies. Brain development analyzed at the morphological level from the third (L3) through fifth (L5) larval instars revealed an asynchrony between queens and workers. In the feeding phase of the last larval instar (L5F), two well-formed structures, pedunculi and calyces, are identifiable in the mushroom bodies of queens, both of which are not present in workers until a later phase (spinning phase, L5S). Genome-wide expression analyses and normalized transcript expression experiments monitoring specific genes revealed that this differential brain development starts earlier, during L3. Analyzing brains from L3 through L5S1 larvae, we identified 21 genes with caste-specific transcription patterns (e.g., APC-4, GlcAT-P, fax, kr-h1 and shot), which encode proteins that are potentially involved in the development of brain tissues through controlling the cell proliferation rate (APC4, kr-h1) and fasciculation (GlcAT-P, fax, and shot). Shot, whose expression is known to be required for axon extension and cell proliferation, was found to be transcribed at significantly higher levels in L4 queens compared with worker larvae. Moreover, the protein encoded by this gene was immunolocalized to the cytoplasm of cells near the antennal lobe neuropiles and proximal to the Kenyon cells in the brains of L4 queens. In conclusion, during the larval period, the brains of queens are larger and develop more rapidly than workers' brains, which represents a developmental heterochrony reflecting the effect of the differential feeding regime of the two castes on nervous system development. Furthermore, this differential development is characterized by caste-specific transcriptional profiles of a set of genes, thus pointing to a link between differential nutrition and differential neurogenesis via genes that control cell proliferation and fasciculation.

Hox gene expression leads to differential hind leg development between honeybee castes.

Beyond the physiological and behavioural, differences in appendage morphology between the workers and queens of Apis mellifera are pre-eminent. The hind legs of workers, which are highly specialized pollinators, deserve special attention. The hind tibia of worker has an expanded bristle-free region used for carrying pollen and propolis, the corbicula. In queens this structure is absent. Although the morphological differences are well characterized, the genetic inputs driving the development of this alternative morphology remain unknown. Leg phenotype determination takes place between the fourth and fifth larval instar and herein we show that the morphogenesis is completed at brown-eyed pupa. Using results from the hybridization of whole genome-based oligonucleotide arrays with RNA samples from hind leg imaginal discs of pre-pupal honeybees of both castes we present a list of 200 differentially expressed genes. Notably, there are castes preferentially expressed cuticular protein genes and members of the P450 family. We also provide results of qPCR analyses determining the developmental transcription profiles of eight selected genes, including abdominal-A, distal-less and ultrabithorax (Ubx), whose roles in leg development have been previously demonstrated in other insect models. Ubx expression in workers hind leg is approximately 25 times higher than in queens. Finally, immunohistochemistry assays show that Ubx localization during hind leg development resembles the bristles localization in the tibia/basitarsus of the adult legs in both castes. Our data strongly indicate that the development of the hind legs diphenism characteristic of this corbiculate species is driven by a set of caste-preferentially expressed genes, such as those encoding cuticular protein genes, P450 and Hox proteins, in response to the naturally different diets offered to honeybees during the larval period.

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.

Non-coding transcription characterization and annotation: a guide and web resource for non-coding RNA databases.

Large-scale transcriptome projects have shown that the number of RNA transcripts not coding for proteins (non-coding RNAs) is much larger than previously recognized. High-throughput technologies, coupled with bioinformatics approaches, have produced increasing amounts of data, highlighting the role of non-coding RNAs (ncRNAs) in biological processes. Data generated by these studies include diverse non-coding RNA classes from organisms of different kingdoms, which were obtained using different experimental and computational assays. This has led to a rapid increase of specialized RNA databases. The fast growth in the number of available databases makes integration of stored information a difficult task. We present here NRDR, a Non-coding RNA Databases Resource for information retrieval on ncRNA databases (www.ncrnadatabases.org). We performed a survey of 102 public databases on ncRNAs and we have introduced four categorizations to classify these databases and to help researchers quickly search and find the information they need: RNA family, information source, information content and available search mechanisms. NRDR is a useful databases searching tool that will facilitate research on ncRNAs.

Mars is close to venus--female reproductive proteins are expressed in the fat body and reproductive tract of honey bee (Apis mellifera L.) drones.

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.

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.

A non-invasive method for silencing gene transcription in honeybees maintained under natural conditions.

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.