Search

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.

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 m 6 A and m 5 C 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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Bataglia, Simões and Nunes.)

Pollinator populations, including bees, are in rapid decline in many parts of the world, raising concerns over the future of ecosystems and food production. Among the factors involved in these declines, poor nutrition deserves attention. The diet consumed by adult worker honeybees (Apis mellifera) is crucial for their behavioral maturation, i.e., the progressive division of labor they perform, such as nurse bees initially and later in life as foragers. Poor pollen nutrition is known to reduce the workers' lifespan, but the underlying physiological and genetic mechanisms are not fully understood. Here we investigate how the lack of pollen in the diet of workers during their first week of adult life can affect age-related phenotypes. During the first seven days of adult life, newly emerged workers were fed either a pollen-deprived (PD) diet mimicking that of an older bee, or a control pollen-rich (PR) diet, as typically consumed by young bees. The PD-fed bees showed alterations in their fat body transcriptome, such as a switch from a protein-lipid based metabolism to a carbohydrate-based metabolism, and a reduced expression of genes involved with immune response. The absence of pollen in the diet also led to an accumulation of oxidative stress markers in fat body tissue and alterations in the cuticular hydrocarbon profiles, which became similar to those of chronologically older bees. Together, our data indicate that the absence of pollen during first week of adulthood triggers the premature onset of an aging-related worker phenotype., (Copyright © 2022 Elsevier Ltd. 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., (© 2021 Royal Entomological Society.)

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 an epitranscriptomic pathway in insects. Here we identified genes related to RNA m 6 A (N6-methyladenine) and m 5 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 m 6 A and m 5 C measured in the fat body and brain of adult workers differ significantly. Also, m 6 A levels were differed significantly mainly between the fourth larval instar of queens and workers. Moreover, we found a conserved m 5 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., (© 2021 Royal Entomological Society.)

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., (© 2021. The Author(s).)

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., (© 2021. The Author(s).)

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., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

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., (© 2020 Royal Entomological Society.)

Background: Most 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., Results: The assembly of 243,974,526 high quality Illumina reads resulted in a predicted assembled genome size of 275 Mb 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,144 bp, 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., Conclusions: The 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.

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.

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.

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., (Copyright © 2019 Elsevier Inc. All rights reserved.)

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., Competing Interests: The authors have declared that no competing interests exist.

Ftz-f1 is an orphan member of the nuclear hormone receptor superfamily. A 20-hydroxyecdysone pulse allows ftz-f1 gene expression, which then regulates the activity of downstream genes involved in major developmental progression events. In honeybees, the expression of genes like vitellogenin (vg), prophenoloxidase and juvenile hormone-esterase during late pharate-adult development is known to be hormonally controlled in both queens and workers by increasing juvenile hormone (JH) titres in the presence of declining levels of ecdysteroids. Since Ftz-f1 is known for mediating intracellular JH signalling, we hypothesized that ftz-f1 could mediate JH action during the pharate-adult development of honeybees, thus controlling the expression of these genes. Here, we show that ftz-f1 has caste-specific transcription profiles during this developmental period, with a peak coinciding with the increase in JH titre, and that its expression is upregulated by JH and downregulated by ecdysteroids. RNAi-mediated knock down of ftz-f1 showed that the expression of genes essential for adult development (e.g. vg and cuticular genes) depends on ftz-f1 expression. Finally, a double-repressor hypothesis-inspired vg gene knock-down experiment suggests the existence of a positive molecular loop between JH, ftz-f1 and vg., (© 2018 The Royal Entomological Society.)

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)., (© 2018 The Royal Entomological Society.)

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., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

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., (© 2017 John Wiley & Sons Ltd.)