Exploring integument transcriptomes, cuticle ultrastructure, and cuticular hydrocarbons profiles in eusocial and solitary bee species displaying heterochronic adult cuticle maturation

Changes in the proteome of Apis mellifera acutely exposed to sublethal dosage of glyphosate and imidacloprid

Uncontrolled use of pesticides has caused a dramatic reduction in the number of pollinators, including bees. Studies on the effects of pesticides on bees have reported effects on both metabolic and neurological levels under chronic exposure. In this study, variations in the differential expression of head and thorax-abdomen proteins in Africanized A. mellifera bees treated acutely with sublethal doses of glyphosate and imidacloprid were studied using a proteomic approach. A total of 92 proteins were detected, 49 of which were differentially expressed compared to those in the control group (47 downregulated and 2 upregulated). Protein interaction networks with differential protein expression ratios suggested that acute exposure of A. mellifera to sublethal doses of glyphosate could cause head damage, which is mainly associated with behavior and metabolism. Simultaneously, imidacloprid can cause damage associated with metabolism as well as, neuronal damage, cellular stress, and impairment of the detoxification system. Regarding the thorax-abdomen fractions, glyphosate could lead to cytoskeleton reorganization and a reduction in defense mechanisms, whereas imidacloprid could affect the coordination and impairment of the oxidative stress response.

Molecular signatures of alternative reproductive strategies in a facultatively social hover wasp

Social insect reproductives and non-reproductives represent ideal models with which to understand the expression and regulation of alternative phenotypes. Most research in this area has focused on the developmental regulation of reproductive phenotypes in obligately social taxa such as honey bees, while relatively few studies have addressed the molecular correlates of reproductive differentiation in species in which the division of reproductive labour is established only in plastic dominance hierarchies. To address this knowledge gap, we generate the first genome for any stenogastrine wasp and analyse brain transcriptomic data for non-reproductives and reproductives of the facultatively social species Liostenogaster flavolineata, a representative of one of the simplest forms of social living. By experimentally manipulating the reproductive 'queues' exhibited by social colonies of this species, we show that reproductive division of labour in this species is associated with transcriptomic signatures that are more subtle and variable than those observed in social taxa in which colony living has become obligate; that variation in gene expression among non-reproductives reflects their investment into foraging effort more than their social rank; and that genes associated with reproductive division of labour overlap to some extent with those underlying division of labour in the separate polistine origin of wasp sociality but only explain a small portion of overall variation in this trait. These results indicate that broad patterns of within-colony transcriptomic differentiation in this species are similar to those in Polistinae but offer little support for the existence of a strongly conserved 'toolkit' for sociality.

Screening and characterization of Bacillus thuringiensis isolates for high production of Vip3A and Cry proteins and high thermostability to control Spodoptera spp

Bacillus thuringiensis (Bt) is an entomopathogenic bacterium that produces crystalline (Cry and Cyt) and soluble (vegetative insecticidal proteins or Vips) proteins during the sporulation and vegetative growth phases, respectively. Combining Cry and Vip proteins could delay insect resistance development and exhibit synergistic activity against various insect pests. This study aims to screen Bt isolates collected from Thailand for high Vip3A and Cry protein production levels and high thermostability to control Spodoptera spp. Among the selected Bt isolates with high target protein synthesis, Bt isolate 506 was found to be safe for further biopesticide formulation due to the absence of non-specific metabolite, as determined by the detection of thermo-stable β-exotoxin I based on biological assays and PCR analysis. Bt isolate 506 showed the presence of Cry1A, Cry2A, and Vip3A-type proteins identified as Cry1Aa45, Cry2Aa22, and Vip3A87, respectively. The insecticidal activity of whole culture extracts containing Vip3A and Cry mixtures and culture supernatants containing secreted Vip3A protein was evaluated against the second-instar larvae of S. exigua and S. frugiperda. The Bt isolate 506 showed high toxicity against both insects, and the insecticidal proteins produced by this isolate retained their activity after heating at 50 °C. This Bt isolate is a promising candidate for further development as a biopesticide against lepidopteran pests.

Cellular heterogeneity of the developing worker honey bee (Apis mellifera) pupa: a single cell transcriptomics analysis

It is estimated that animals pollinate 87.5% of flowering plants worldwide and that managed honey bees (Apis mellifera) account for 30-50% of this ecosystem service to agriculture. In addition to their important role as pollinators, honey bees are well-established insect models for studying learning and memory, behavior, caste differentiation, epigenetic mechanisms, olfactory biology, sex determination, and eusociality. Despite their importance to agriculture, knowledge of honey bee biology lags behind many other livestock species. In this study, we have used scRNA-Seq to map cell types to different developmental stages of the worker honey bee (prepupa at day 11 and pupa at day 15) and sought to determine their gene expression signatures. To identify cell-type populations, we examined the cell-to-cell network based on the similarity of the single-cells transcriptomic profiles. Grouping similar cells together we identified 63 different cell clusters of which 17 clusters were identifiable at both stages. To determine genes associated with specific cell populations or with a particular biological process involved in honey bee development, we used gene coexpression analysis. We combined this analysis with literature mining, the honey bee protein atlas, and gene ontology analysis to determine cell cluster identity. Of the cell clusters identified, 17 were related to the nervous system and sensory organs, 7 to the fat body, 19 to the cuticle, 5 to muscle, 4 to compound eye, 2 to midgut, 2 to hemocytes, and 1 to malpighian tubule/pericardial nephrocyte. To our knowledge, this is the first whole single-cell atlas of honey bees at any stage of development and demonstrates the potential for further work to investigate their biology at the cellular level.

Feeding Asian honeybee queens with European honeybee royal jelly alters body color and expression of related coding and non-coding RNAs

Background and aims: The Asian honeybee (Apis cerana) and the European honeybee (Apis mellifera) are reproductively isolated. Previous studies reported that exchanging the larval food between the two species, known as nutritional crossbreeding, resulted in obvious changes in morphology, physiology and behavior. This study explored the molecular mechanisms underlying the honeybee nutritional crossbreeding. Methods: This study used full nutritional crossbreeding technology to rear A. cerana queens by feeding them with an A. mellifera royal jelly-based diet in an incubator. The body color and the expression of certain genes, microRNA, lncRNA, and circRNA among nutritional crossbred A. cerana queens (NQ), and control A. cerana queens (CQ) were compared. The biological functions of two target genes, TPH1 and KMO, were verified using RNA interference. Results: Our results showed that the NQ's body color turned yellow compared to the black control queens. Whole transcriptome sequencing results showed that a total of 1484, 311, 92, and 169 DEGs, DElncRNAs, DEmiRNAs, and DEcircRNAs, respectively, were identified in NQ and CQ, in which seven DEGs were enriched for three key pathways (tryptophan, tyrosine, and dopamine) involved in melanin synthesis. Interestingly, eight DElncRNAs and three DEmiRNAs were enriched into the key pathways regulating the above key DEGs. No circRNAs were enriched into these key pathways. Knocking down two key genes (KMO and TPH1) resulted in altered body color, suggesting that feeding NQ's an RNAi-based diet significantly downregulated the expression of TPH1 and KMO in 4-day-old larvae, which confirmed the function of key DEGs in the regulation of honeybee body color. Conclusion: These findings reveal that the larval diets from A. mellifera could change the body color of A. cerana, perhaps by altering the expression of non-coding RNAs and related key genes. This study serves as a model of epigenetic regulation in insect body color induced by environmental factors.

Insights into insecticide-resistance mechanisms in invasive species: Challenges and control strategies

Threatening the global community is a wide variety of potential threats, most notably invasive pest species. Invasive pest species are non-native organisms that humans have either accidentally or intentionally spread to new regions. One of the most effective and first lines of control strategies for controlling pests is the application of insecticides. These toxic chemicals are employed to get rid of pests, but they pose great risks to people, animals, and plants. Pesticides are heavily used in managing invasive pests in the current era. Due to the overuse of synthetic chemicals, numerous invasive species have already developed resistance. The resistance development is the main reason for the failure to manage the invasive species. Developing pesticide resistance management techniques necessitates a thorough understanding of the mechanisms through which insects acquire insecticide resistance. Insects use a variety of behavioral, biochemical, physiological, genetic, and metabolic methods to deal with toxic chemicals, which can lead to resistance through continuous overexpression of detoxifying enzymes. An overabundance of enzymes causes metabolic resistance, detoxifying pesticides and rendering them ineffective against pests. A key factor in the development of metabolic resistance is the amplification of certain metabolic enzymes, specifically esterases, Glutathione S-transferase, Cytochromes p450 monooxygenase, and hydrolyses. Additionally, insect guts offer unique habitats for microbial colonization, and gut bacteria may serve their hosts a variety of useful services. Most importantly, the detoxification of insecticides leads to resistance development. The complete knowledge of invasive pest species and their mechanisms of resistance development could be very helpful in coping with the challenges and effectively developing effective strategies for the control of invasive species. Integrated Pest Management is particularly effective at lowering the risk of chemical and environmental contaminants and the resulting health issues, and it may also offer the most effective ways to control insect pests.

Worker bees (Apis mellifera) deprived of pollen in the first week of adulthood exhibit signs of premature aging

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.

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.

Identification of Bacillus thuringiensis Strains for the Management of Lepidopteran Pests

Bacillus thuringiensis (Bt)-based bioinsecticides and transgenic plants expressing proteins with insecticidal activity (Cry and Vip) have been successfully used in several integrated pest management programs worldwide. Lepidoptera comprise some of the most economically important insect pests of the major agricultural crops. In this study, the toxicity of 150 Bt strains was evaluated against Helicoverpa armigera (Hübner) larvae. Eight strains (426, 520B, 1636, 1641, 1644, 1648, 1657 and 1658) showed high insecticide activity against H. armigera and were therefore tested against Anticarsia gemmatalis (Hübner), Spodoptera cosmioides (Walker), Chrysodeixis includens (Walker), and Diatraea saccharalis (Fabricius) larvae. Our results showed that most of the Bt strains were also toxic to these lepidopteran species. The biochemical and molecular analyses of these strains revealed that they had a similar protein profile; however, their cry and vip gene contents were variable. In addition, the median lethal concentration (LC₅₀) of the selected strains indicated that the strains 1636, 1641, and 1658 were the most effective against H. armigera, showing LC₅₀ values of 185.02, 159.44, and 192.98 ng/cm², respectively. Our results suggest that the selected Bt strains have great potential to control the lepidopteran pests H. armigera, A. gemmatalis, D. saccharalis, S. cosmioides, and C. includes.

Disentangling multiple PM emission sources in the Po Valley (Italy) using honey bees

Particulate matter (PM) is a complex mixture of airborne chemical compounds commonly classified by their aerodynamic diameter. Although PM toxicity strongly depends on the morphology, chemical composition, and dimensions of particles, exposure limits set by environmental organisations only refer to the mean mass concentration of PM sampled daily or annually by monitoring stations.

In this study, we used honey bees as sensors of airborne PM₁₀ and PM_2.5 in a highly polluted area of the Po Valley, northern Italy. Honey bees are an efficient sampler of airborne PM because, during flight and foraging activities, their pubescence promotes the accumulation of electrical charge on the body surface owing to air resistance, thus enhancing airborne PM attraction. Particles attached to the body of bees are readily accessible for physico-chemical characterisation using a scanning electron microscope coupled with X-ray spectroscopy (SEM/EDX). Our results demonstrate that residents in the study area are intermittently but chronically exposed to a well-defined spectrum of metal-bearing particles and mineral phases known to induce specific health outcomes.

The morphology, size, and chemical composition of PM₁₀ and PM_2.5 detected on bees in the monitoring area were indicative of traffic, agricultural operations, and high-temperature combustion processes. The contribution of the A1 Milano-Bologna highway, local wheat and alfalfa cultivation, and the Parma incineration plant were clearly distinguishable. Our data also demonstrated that PM exposure levels may vary sharply throughout the year based on recurrent local activities.

Molecular identification, expression and function analysis of peroxidasin in Chilo suppressalis

Peroxidasin plays a unique role in the formation and stability of extracellular matrix (ECM) in the animal kingdom; however, it was only characterized in Diptera, not in other insect orders. In this study peroxidasin (CsPxd) was first identified and characterized from Chilo suppressalis, a lepidopteran pest. CsPxd complementary DNA with a 4080 bp open reading frame encodes a peptide of 1359 amino acids; the derived amino acid sequence of CsPxd harbors the typical structural characteristics of peroxidasin family in heme-peroxidase superfamily, including the signal peptide at N-terminal, leucine-rich repeat domain, Ig-loop motifs and peroxidase domain, signifying the extracellular location of protein and the involvement in ECM formation. Eukaryotic expression reveals CsPxd protein displays peroxidase activity on H2 O2 , justifying the membership of peroxidase. Phyletic analysis shows the monophyletic evolution pattern of peroxidasin in insect phyle, and moreover only one peroxidasin is present in each species of insects, suggesting its evolutionary conservation on function. Peroxidasin messenger RNA is mainly expressed in egg and the final instar larvae stage. Injection of peroxidasin double-stranded RNA into the final instar larvae impacts the cuticle sclerotization during the metamorphosis from larvae to pupa, and eventually lead to lethality of larvae and pupa. These results suggest the presence of collagen crosslink in chorion and cuticle of insects, and indicate peroxidasin plays a role in the development of chorion and cuticle; furthermore peroxidasin might be the one of potential target genes for pest control using RNA interference.

Sub-lethal effects of the consumption of Eupatorium buniifolium essential oil in honeybees

When developing new products to be used in honeybee colonies, further than acute toxicity, it is imperative to perform an assessment of risks, including various sublethal effects. The long-term sublethal effects of xenobiotics on honeybees, more specifically of acaricides used in honeybee hives, have been scarcely studied, particularly so in the case of essential oils and their components. In this work, chronic effects of the ingestion of Eupatorium buniifolium (Asteraceae) essential oil were studied on nurse honeybees using laboratory assays. Survival, food consumption, and the effect on the composition of cuticular hydrocarbons (CHC) were assessed. CHC were chosen due to their key role as pheromones involved in honeybee social recognition. While food consumption and survival were not affected by the consumption of the essential oil, CHC amounts and profiles showed dose-dependent changes. All groups of CHC (linear and branched alkanes, alkenes and alkadienes) were altered when honeybees were fed with the highest essential oil dose tested (6000 ppm). The compounds that significantly varied include n-docosane, n-tricosane, n-tetracosane, n-triacontane, n-tritriacontane, 9-tricosene, 7-pentacosene, 9-pentacosene, 9-heptacosene, tritriacontene, pentacosadiene, hentriacontadiene, tritriacontadiene and all methyl alkanes. All of them but pentacosadiene were up-regulated. On the other hand, CHC profiles were similar in healthy and Nosema-infected honeybees when diets included the essential oil at 300 and 3000 ppm. Our results show that the ingestion of an essential oil can impact CHC and that the effect is dose-dependent. Changes in CHC could affect the signaling process mediated by these pheromonal compounds. To our knowledge this is the first report of changes in honeybee cuticular hydrocarbons as a result of essential oil ingestion.

Synthesis, structural revision, and tyrosinase inhibitory activity of proposed phloretin-4-O-β-D-glucopyranoside from Homalium stenophyllum

Phloretin-4-O-β-D-glucopyranoside (1), isolated from Homalium stenophyllum, was synthesized for the first time through aldol condensation and Schmidt glycosylation reactions aiming to develop a novel hydrophilic tyrosinase inhibitor. However, the specific rotation of synthetic 1 was found to be negative and different from that reported for natural product 1. Thus, L-glucoside 2 was chemically synthesized using the established synthetic route of 1, suggesting that the configuration of the natural product 1 was the same as that of 2, as their specific rotation and spectroscopic data were also the same. In addition, the evaluation of the inhibitory activity of 1 and 2 against tyrosinase indicated that 2 was 1.4 times more potent than 1, but they were both relatively weak. Therefore, the enantiomeric analogues 1 and 2 were proved to be unique tyrosinase inhibitors due to the chiral recognition from the tyrosinase active site.

The nuclear and mitochondrial genomes of Frieseomelitta varia - a highly eusocial stingless bee (Meliponini) with a permanently sterile worker caste

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.

Meta-analysis of honey bee neurogenomic response links Deformed wing virus type A to precocious behavioral maturation

Crop pollination by the western honey bee Apis mellifera is vital to agriculture but threatened by alarmingly high levels of colony mortality, especially in Europe and North America. Colony loss is due, in part, to the high viral loads of Deformed wing virus (DWV), transmitted by the ectoparasitic mite Varroa destructor, especially throughout the overwintering period of a honey bee colony. Covert DWV infection is commonplace and has been causally linked to precocious foraging, which itself has been linked to colony loss. Taking advantage of four brain transcriptome studies that unexpectedly revealed evidence of covert DWV-A infection, we set out to explore whether this effect is due to DWV-A mimicking naturally occurring changes in brain gene expression that are associated with behavioral maturation. Consistent with this hypothesis, we found that brain gene expression profiles of DWV-A infected bees resembled those of foragers, even in individuals that were much younger than typical foragers. In addition, brain transcriptional regulatory network analysis revealed a positive association between DWV-A infection and transcription factors previously associated with honey bee foraging behavior. Surprisingly, single-cell RNA-Sequencing implicated glia, not neurons, in this effect; there are relatively few glial cells in the insect brain and they are rarely associated with behavioral plasticity. Covert DWV-A infection also has been linked to impaired learning, which together with precocious foraging can lead to increased occurrence of infected bees from one colony mistakenly entering another colony, especially under crowded modern apiary conditions. These findings provide new insights into the mechanisms by which DWV-A affects honey bee health and colony survival.

Perspectives on plasma-assisted synthesis of N-doped nanoparticles as nanopesticides for pest control in crops

Green plasma-based technology production of N-doped NPs for a new agri-tech revolution in pest control.

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

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