A deficit of detoxification enzymes: pesticide sensitivity and environmental response in the honeybee

Comparative transcriptome profiling reveals candidate genes related to insecticide resistance of Glyphodes pyloalis

Glyphodes pyloalis Walker (Lepidoptera: Pyralididae) is a common pest in sericulture and has developed resistance to different insecticides. However, the mechanisms involved in insecticide resistance of G. pyloalis are poorly understood. Here, we present the first whole-transcriptome analysis of differential expression genes in insecticide-resistant and susceptible G. pyloalis. Clustering and enrichment analysis of DEGs revealed several biological pathways and enriched Gene Ontology terms were related to detoxification or insecticide resistance. Genes involved in insecticide metabolic processes, including cytochrome P450, glutathione S-transferases and carboxylesterase, were identified in the larval midgut of G. pyloalis. Among them, CYP324A19, CYP304F17, CYP6AW1, CYP6AB10, GSTs5, and AChE-like were significantly increased after propoxur treatment, while CYP324A19, CCE001c, and AChE-like were significantly induced by phoxim, suggesting that these genes were involved in insecticide metabolism. Furthermore, the sequence variation analysis identified 21 single nucleotide polymorphisms within CYP9A20, CYP6AB47, and CYP6AW1. Our findings reveal many candidate genes related to insecticide resistance of G. pyloalis. These results provide novel insights into insecticide resistance and facilitate the development of insecticides with greater specificity to G. pyloalis.

Impacts of dietary supplementation with p-coumaric acid and indole-3-acetic acid on survival and biochemical response of honey bees treated with tau-fluvalinate

Pollinator populations are in decline worldwide. Multiple factors have been cited as potential causes to these declines. In honey bees, a combination of stressors is known to cause colony losses. Adequate nutrition is a key factor for honey bee growth and colony development. Several studies show that the nutritional quality of the diet is directly proportional to the ability of the bee to face challenges or stressors. We explored the effect of p-coumaric (600 μM) and indole-3-acetic acid (2, 20 or 200 μM) supplementation on the survival and activity of key detoxification enzymes of honey bees exposed to tau-fluvalinate. The dietary supplementation with p-coumaric and indole-3-acetic acids (20 μM) enhanced the survival of bees exposed to tau-fluvalinate (approximately 20%). We also showed that dietary p-coumaric acid increased the levels of cytochrome P450 and glutathione reductase activity in bees treated with tau-fluvalinate, as well as in the untreated controls, while glutathione-S-transferase activity was lower in treated bees than in untreated. In bees fed with indole-3-acetic acid, cytochrome P450 showed increased levels, however, glutathione-S-transferase showed the lowest activity. Moreover, the results showed that supplementation with p-coumaric and indole-3-acetic acids did not alter acetyl cholinesterase activity, nor did treatment with tau-fluvalinate. Altogether, the enzymatic changes related to the detoxification mechanisms observed in bees that were fed with p-coumaric and indole-3-acetic acids could be responsible for the increased survival of bees treated with tau-fluvalinate compared to those that received a control diet. The results presented in this study, together with previous studies, provide evidence of the importance of dietary phytochemicals in the response of honey bees to pesticide exposure. Moreover, these results are the first report of the beneficial effect of the phytohormone indole-3-acetic acid on the survival of honey bees treated with tau-fluvalinate.

Frequently encountered pesticides can cause multiple disorders in developing worker honey bees

Pesticide exposure is regarded as a contributing factor to the high gross loss rates of managed colonies of Apis mellifera. Pesticides enter the hive through contaminated nectar and pollen carried by returning forager honey bees or placed in the hive by beekeepers when managing hive pests. We used an in vitro rearing method to characterize the effects of seven pesticides on developing brood subjected dietary exposure at worse-case environmental concentrations detected in wax and pollen. The pesticides tested included acaricides (amitraz, coumaphos, fluvalinate), insecticides (chlorpyrifos, imidacloprid), one fungicide (chlorothalonil), and one herbicide (glyphosate). The larvae were exposed chronically for six days of mimicking exposure during the entire larval feeding period, which is the worst possible scenario of larval exposure. Survival, duration of immature development, the weight of newly emerged adult, morphologies of the antenna and the hypopharyngeal gland, and gene expression were recorded. Survival of bees exposed to amitraz, coumaphos, fluvalinate, chlorpyrifos, and chlorothalonil was the most sensitive endpoint despite observed changes in many developmental and physiological parameters across the seven pesticides. Our findings suggest that pesticide exposure during larvae development may affect the survival and health of immature honey bees, thus contributing to overall colony stress or loss. Additionally, pesticide exposure altered gene expression of detoxification enzymes. However, the tested exposure scenario is unlikely to be representative of real-world conditions but emphasizes the importance of proper hive management to minimize pesticide contamination of the hive environment or simulates a future scenario of increased contamination.

Translocation of pharmaceuticals from wastewater into beehives

There has been a substantial research focus on the presence of pesticides in flowers and the subsequent exposure to honeybees. Here we demonstrate for the first time that honeybees can also be exposed to pharmaceuticals, commonly present in wastewater. Residues of carbamazepine (an anti-epileptic drug) up to 371 ng/mL and 30 µg/g were detected in nectar and pollen sampled from zucchini flowers (Cucurbita pepo) grown in carbamazepine spiked soil (0.5-20 µg/g). Under realistic exposure conditions from the use of recycled wastewater, carbamazepine concentrations were estimated to be 0.37 ng/L and 30 ng/kg in nectar and pollen, respectively. Incorporation of environmentally relevant carbamazepine residues in nectar and pollen into a modelling framework able to simulate beehive dynamics including the honeybee foraging activity at the landscape scale (BEEHAVE and BEESCOUT) enabled the simulation of carbamazepine translocation from zucchini fields into honeybee hives. Carbamazepine accumulation was modelled in 11 beehives across a 25 km2 landscape over three years chosen to represent distinct climatic conditions. During a single flowering period, carbamazepine concentrations were simulated to range between 0 and 2478 ng per beehive. The amount of carbamazepine gathered not only varied across the simulated years but there were also differences in accumulation of carbamazepine between beehives within the same year. This work illustrates a fundamental first step in assessing the risk of pharmaceuticals to bees through realistic scenarios by demonstrating a method to quantify potential exposure of honeybees at the landscape scale. Pharmaceuticals are being inadvertently but increasingly applied to agricultural lands globally via the use of wastewater for agricultural irrigation in response to water scarcity problems. We have demonstrated a route of pharmaceutical exposure to honeybees via contaminated nectar and pollen. Given the biological potency of pharmaceuticals, accumulation of these chemicals in nectar and pollen suggest potential implications for honeybee health, with unknown ecosystem consequences.

Effects of chronic dietary thiamethoxam and prothioconazole exposure on Apis mellifera worker adults and brood

BACKGROUND

Chronic exposure of honey bees (Apis mellifera Linnaeus) to the neonicotinoid thiamethoxam and the fungicide prothioconazole is common during foraging in agricultural landscapes. We evaluated the survival and hypopharyngeal gland development of adult worker honey bees, and the survival of the worker brood when chronically exposed to thiamethoxam or thiamethoxam and prothioconazole in combination.

RESULTS

We found that 30 days of exposure to 40 μg kg^-1 of thiamethoxam significantly (P < 0.001) increased the frequency of death in worker adults by four times relative to solvent control. The worker brood required 23 times higher doses of thiamethoxam (1 mg L^-1 or 909 μg kg^-1 ) before a significant (P = 0.04), 3.9 times increase in frequency of death was observed relative to solvent control. No additive effects of simultaneous exposure of worker adults or brood to thiamethoxam and prothioconazole were observed. At day 8 and day 12, the hypopharyngeal gland acinar diameter was not significantly different (P > 0.05) between controls and adult workers exposed to thiamethoxam and/or prothioconazole.

CONCLUSION

These results indicate that chronic exposure to field-realistic doses of thiamethoxam and/or prothioconazole are unlikely to affect the survival of adult workers and brood. © 2019 Society of Chemical Industry.

Host plant adaptation in the polyphagous whitefly, Trialeurodes vaporariorum, is associated with transcriptional plasticity and altered sensitivity to insecticides

BACKGROUND

The glasshouse whitefly, Trialeurodes vaporariorum, is a damaging crop pest and an invasive generalist capable of feeding on a broad range of host plants. As such this species has evolved mechanisms to circumvent the wide spectrum of anti-herbivore allelochemicals produced by its host range. T. vaporariorum has also demonstrated a remarkable ability to evolve resistance to many of the synthetic insecticides used for control.

RESULTS

To gain insight into the molecular mechanisms that underpin the polyphagy of T. vaporariorum and its resistance to natural and synthetic xenobiotics, we sequenced and assembled a reference genome for this species. Curation of genes putatively involved in the detoxification of natural and synthetic xenobiotics revealed a marked reduction in specific gene families between this species and another generalist whitefly, Bemisia tabaci. Transcriptome profiling of T. vaporariorum upon transfer to a range of different host plants revealed profound differences in the transcriptional response to more or less challenging hosts. Large scale changes in gene expression (> 20% of genes) were observed during adaptation to challenging hosts with a range of genes involved in gene regulation, signalling, and detoxification differentially expressed. Remarkably, these changes in gene expression were associated with significant shifts in the tolerance of host-adapted T. vaporariorum lines to natural and synthetic insecticides.

CONCLUSIONS

Our findings provide further insights into the ability of polyphagous insects to extensively reprogram gene expression during host adaptation and illustrate the potential implications of this on their sensitivity to synthetic insecticides.

Acute thiamethoxam toxicity in honeybees is not enhanced by common fungicide and herbicide and lacks stress-induced changes in mRNA splicing

Securing food supply for a growing population is a major challenge and heavily relies on the use of agrochemicals to maximize crop yield. It is increasingly recognized, that some neonicotinoid insecticides have a negative impact on non-target organisms, including important pollinators such as the European honeybee Apis mellifera. Toxicity of neonicotinoids may be enhanced through simultaneous exposure with additional pesticides, which could help explain, in part, the global decline of honeybee colonies. Here we examined whether exposure effects of the neonicotinoid thiamethoxam on bee viability are enhanced by the commonly used fungicide carbendazim and the herbicide glyphosate. We also analysed alternative splicing changes upon pesticide exposure in the honeybee. In particular, we examined transcripts of three genes: (i) the stress sensor gene X box binding protein-1 (Xbp1), (ii) the Down Syndrome Cell Adhesion Molecule (Dscam) gene and iii) the embryonic lethal/abnormal visual system (elav) gene, which are important for neuronal function. Our results showed that acute thiamethoxam exposure is not enhanced by carbendazim, nor glyphosate. Toxicity of the compounds did not trigger stress-induced, alternative splicing in the analysed mRNAs, thereby leaving dormant a cellular response pathway to these man-made environmental perturbations.

Genome-wide survey of cytochrome P450 genes in the salmon louse Lepeophtheirus salmonis (Krøyer, 1837)

Background

The salmon louse (Lepeophtheirus salmonis) infests farmed and wild salmonid fishes, causing considerable economic damage to the salmon farming industry. Infestations of farmed salmon are controlled using a combination of non-medicinal approaches and veterinary drug treatments. While L. salmonis has developed resistance to most available salmon delousing agents, relatively little is known about the molecular mechanisms involved. Members of the cytochrome P450 (CYP) superfamily are typically monooxygenases, some of which are involved in the biosynthesis and metabolism of endogenous compounds, while others have central roles in the detoxification of xenobiotics. In terrestrial arthropods, insecticide resistance can be based on the enhanced expression of CYPs. The reported research aimed to characterise the CYP superfamily in L. salmonis and assess its potential roles in drug resistance.

Methods

Lepeophtheirus salmonis CYPs were identified by homology searches of the genome and transcriptome of the parasite. CYP transcript abundance in drug susceptible and multi-resistant L. salmonis was assessed by quantitative reverse transcription PCR, taking into account both constitutive expression and expression in parasites exposed to sublethal levels of salmon delousing agents, ecdysteroids and environmental chemicals.

Results

The above strategy led to the identification of 25 CYP genes/pseudogenes in L. salmonis, making its CYP superfamily the most compact characterised for any arthropod to date. Lepeophtheirus salmonis possesses homologues of a number of arthropod CYP genes with roles in ecdysteroid metabolism, such as the fruit fly genes disembodied, shadow, shade, spook and Cyp18a1. CYP transcript expression did not differ between one drug susceptible and one multi-resistant strain of L. salmonis. Exposure of L. salmonis to emamectin benzoate or deltamethrin caused the transcriptional upregulation of certain CYPs. In contrast, neither ecdysteroid nor benzo[a]pyrene exposure affected CYP transcription significantly.

Conclusions

The parasite L. salmonis is demonstrated to possess the most compact CYP superfamily characterised for any arthropod to date. The complement of CYP genes in L. salmonis includes conserved CYP genes involved in ecdysteroid biosynthesis and metabolism, as well as drug-inducible CYP genes. The present study does not provide evidence for a role of CYP genes in the decreased susceptibility of the multiresistant parasite strain studied.

Acetylcholine and Its Receptors in Honeybees: Involvement in Development and Impairments by Neonicotinoids

Acetylcholine (ACh) is the major excitatory neurotransmitter in the insect central nervous system (CNS). However, besides the neuronal expression of ACh receptors (AChR), the existence of non-neuronal AChR in honeybees is plausible. The cholinergic system is a popular target of insecticides because the pharmacology of insect nicotinic acetylcholine receptors (nAChRs) differs substantially from their vertebrate counterparts. Neonicotinoids are agonists of the nAChR and are largely used in crop protection. In contrast to their relatively high safety for humans and livestock, neonicotinoids pose a threat to pollinating insects such as bees. In addition to its effects on behavior, it becomes increasingly evident that neonicotinoids affect developmental processes in bees that appear to be independent of neuronal AChRs. Brood food (royal jelly, worker jelly, or drone jelly) produced in the hypopharyngeal glands of nurse bees contains millimolar concentrations of ACh, which is required for proper larval development. Neonicotinoids reduce the secreted ACh-content in brood food, reduce hypopharyngeal gland size, and lead to developmental impairments within the colony. We assume that potential hazards of neonicotinoids on pollinating bees occur neuronally causing behavioral impairments on adult individuals, and non-neuronally causing developmental disturbances as well as destroying gland functioning.

Hydroxylation patterns associated with pheromone synthesis and composition in two honey bee subspecies Apis mellifera scutellata and A. m. capensis laying workers

Colony losses due to social parasitism in the form of reproductive workers of the Apis mellifera capensis clones results from the production of queen-like pheromonal signals coupled with ovarian activation in these socially parasitic honey bees. While the behavioral attributes of these social parasites have been described, their genetic attributes require more detailed exploration. Here, we investigate the production of mandibular gland pheromones in queenless workers of two sub-species of African honey bees; A. m. scutellata (low reproductive potential) and A. m. capensis clones (high reproductive potential). We used standard techniques in gas chromatography to assess the amounts of various pheromone components present, and qPCR to assess the expression of cytochrome P450 genes cyp6bd1 and cyp6as8, thought to be involved in the caste-dependent hydroxylation of acylated stearic acid in queens and workers, respectively. We found that, for both subspecies, the quality and quantity of the individual pheromone components vary with age, and that from the onset, A. m. capensis parasites make use of gene pathways typically upregulated in queens in achieving reproductive dominance. Due to the high production of 9-hydroxy-decenoic acid (9-HDA) the precursor to the queen substance 9-oxo-decenoic acid (9-ODA) in newly emerged capensis clones, we argue that clones are primed for parasitism upon emergence and develop into fully fledged parasites depending on the colony's social environment.

New approaches to old problems: Removal of phospholipase A2 results in highly active microsomal membranes from the honey bee, Apis mellifera

Over the last 50 years numerous studies were published by insect toxicologists using native microsomal membrane preparations in order to investigate in vitro cytochrome P450-(P450) mediated oxidative metabolism of xenobiotics, including insecticides. Whereas the preparation of active microsomal membranes from many pest insect species is straightforward, their isolation from honey bees, Apis mellifera (Hymenoptera: Apidae) remained difficult, if not impossible, due to the presence of a yet unidentified endogenous inhibitory factor released during abdominal gut membrane isolation. Thus hampering in vitro toxicological studies on microsomal oxidative phase 1 metabolism of xenobiotics, including compounds of ecotoxicological concern. The use of microsomal membranes rather than individually expressed P450s offers advantages and allows to develop a better understanding of phase 1 driven metabolic fate of foreign compounds. Here we biochemically investigated the problems associated with the isolation of active honey bee microsomes and developed a method resulting in highly active native microsomal preparations from adult female worker abdomens. This was achieved by removal of the abdominal venom gland sting complex prior to microsomal membrane preparation. Molecular sieve chromatography of the venom sac content leads to the identification of phospholipase A₂ as the enzyme responsible for the immediate inhibition of cytochrome P450 activity in microsomal preparations. The substrate specificity of functional honey bee microsomes was investigated with different fluorogenic substrates, and revealed a strong preference for coumarin over resorufin derivatives. Furthermore we were able to demonstrate the metabolism of insecticides by honey bee microsomes using an approach coupled to LC-MS/MS analysis of hydroxylated metabolites. Our work provides access to a new and simple in vitro tool to study honey bee phase 1 metabolism of xenobiotics utilising the entire range of microsomal cytochrome P450s.

The leafcutter bee, Megachile rotundata, is more sensitive to N-cyanoamidine neonicotinoid and butenolide insecticides than other managed bees

Recent research has shown that several managed bee species have specific P450 enzymes that are preadapted to confer intrinsic tolerance to some insecticides including certain neonicotinoids. However, the universality of this finding across managed bee pollinators is unclear. Here we show that the alfalfa leafcutter bee, Megachile rotundata, lacks such P450 enzymes and is >2,500-fold more sensitive to the neonicotinoid thiacloprid and 170-fold more sensitive to the butenolide insecticide flupyradifurone than other managed bee pollinators. These findings have important implications for the safe use of insecticides in crops where M. rotundata is used for pollination, and ensuring that regulatory pesticide risk assessment frameworks are protective of this species.

A systemic study of indoxacarb resistance in Spodoptera litura revealed complex expression profiles and regulatory mechanism

The tobacco cutworm, Spodoptera litura, is an important pest of crop and vegetable plants worldwide, and its resistance to insecticides have quickly developed. However, the resistance mechanisms of this pest are still unclear. In this study, the change in mRNA and miRNA profiles in the susceptible, indoxacarb-resistant and field indoxacarb-resistant strains of S. litura were characterized. Nine hundred and ten co-up-regulated and 737 co-down-regulated genes were identified in the resistant strains. Further analysis showed that 126 co-differentially expressed genes (co-DEGs) (cytochrome P450, carboxy/cholinesterase, glutathione S-transferase, ATP-binding cassette transporter, UDP-glucuronosyl transferase, aminopeptidase N, sialin, serine protease and cuticle protein) may play important roles in indoxacarb resistance in S. litura. In addition, a total of 91 known and 52 novel miRNAs were identified, and 10 miRNAs were co-differentially expressed in the resistant strains of S. litura. Furthermore, 10 co-differentially expressed miRNAs (co-DEmiRNAs) had predicted co-DEGs according to the expected miRNA-mRNA negative regulation pattern and 37 indoxacarb resistance-related co-DEGs were predicted to be the target genes. These results not only broadened our understanding of molecular mechanisms of insecticide resistance by revealing complicated profiles, but also provide important clues for further study on the mechanisms of miRNAs involved in indoxacarb resistance in S. litura.

Biopesticide spinosad induces transcriptional alterations in genes associated with energy production in honey bees (Apis mellifera) at sublethal concentrations

Bees experience substantial colony losses, which are often associated with pesticides. Besides synthetic insecticides biological compounds such as spinosad are used in agriculture and organic farming against insect pests. However, potential adverse effect at sublethal concentrations to pollinators are poorly known. Here we aim to determine potential adverse outcome pathways of spinosad and to identify molecular effects by investigating transcriptional alterations in the brain of honey bees. We experimentally exposed bees to three sublethal concentrations of 0.05, 0.5 and 5 ng spinosad/bee, and assessed transcriptional alterations of target genes. Additionally, we evaluated whether spinosad-induced transcriptional alterations were influenced by the time of the year. In April, alterations were most pronounced after 24 h exposure, while in June alterations occurred mostly after 48 h. In July, expressional alterations were often lower but the pattern was more similar to that in June than that in April. Down-regulation of genes encoding acetylcholine receptors, enzymes involved in oxidative phosphorylation (cox5a, ndufb7 and cox17), cytochrome P450 dependent monooxygenases (cyp9q1, cyp9q2 and cyp9q3) and insulin-like peptide-1 were among the most significant transcriptional alterations. This suggests adverse effects of spinosad to energy production and metabolism and thus negative consequences on foraging. Together, our study indicates that spinosad causes adverse effects at environmentally realistic concentrations, which may pose a risk to bee populations.

Genome-Wide Analysis of Carboxylesterases (COEs) in the Whitefly, Bemisia tabaci (Gennadius)

The whitefly (Bemisia tabaci), an important invasive pest that causes severe damage to crops worldwide, has developed resistance to a variety of insecticides. Carboxylesterases (COEs) are important multifunctional enzymes involved in the growth, development, and xenobiotic metabolism of insects. However, systematic studies on the COEs of B. tabaci are scarce. Here, 42 putative COEs in different functional categories were identified in the Mediterranean species of B. tabaci (B. tabaci MED) based on a genome database and neighbor-joining phylogeny. The expression patterns of the COEs were affected by the development of B. tabaci. The expression levels of six COEs were positively correlated with the concentration of imidacloprid to which B. tabaci adults were exposed. The mortality of B. tabaci MED adults fed dsBTbe5 (67.5%) and dsBTjhe2 (58.4%) was significantly higher than the adults fed dsEGFP (41.1%) when treated with imidacloprid. Our results provide a basis for functional research on COEs in B. tabaci and provide new insight into the imidacloprid resistance of B. tabaci.

Targeting the potassium ion channel genes SK and SH as a novel approach for control of insect pests: efficacy and biosafety

BACKGROUND

Potassium ion channels play a critical role in the generation of electrical signals and thus provide potential targets for control of insect pests by RNA interference.

RESULTS

Genes encoding the small conductance calcium-activated potassium channel (SK) and the voltage-gated potassium channel (SH) were knocked down in Tribolium castaneum by injection and oral delivery of dsRNA (dsTcSK and dsTcSH, respectively). Irrespective of the delivery mechanism a dose-dependent effect was observed for knockdown (KD) of gene expression and insect mortality for both genes. Larvae fed a 400 ng dsRNA mg^-1 diet showed significant gene (P < 0.05) knockdown (98% and 83%) for SK and SH, respectively, with corresponding mortalities of 100% and 98% after 7 days. When injected (248.4 ng larva^-1 ), gene KD was 99% and 98% for SK and SH, causing 100% and 73.4% mortality, respectively. All developmental stages tested (larvae, early- and late-stage pupae and adults) showed an RNAi-sensitive response for both genes. LC50 values were lower for SK than SH, irrespective of delivery method, demonstrating that the knockdown of SK had a greater effect on larval mortality. Biosafety studies using adult honeybee Apis mellifera showed that there were no significant differences either in expression levels or mortality of honeybees orally dosed with dsTcSK and dsTcSH compared to control-fed bees. Similarly, there was no significant difference in the titre of deformed wing virus, used as a measure of immune suppression, between experimental and control bees.

CONCLUSION

This study demonstrates the potential of using RNAi targeting neural receptors as a technology for the control of T. castaneum. © 2019 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Identification of glutathione-S-transferase genes by transcriptome analysis in Meteorus pulchricornis (Hymenoptera: Braconidae) and their expression patterns under stress of phoxim and cypermethrin

Meteorus pulchricornis (Wesmael) (Hymenoptera: Braconidae) is a preponderant endoparasitoid wasp, attacking the larvae of many lepidopteran pests. We present the first body transcriptome dataset for M. pulchricornis. In total, 50,781,796 clean reads were obtained and 33,144 unigenes were assembled; 15,458 unigenes showed a significant similarity (E value < 10^-5) to known proteins in the NCBI non-redundant protein database. Gene ontology and cluster of orthologous group analyses were performed to classify the functions of genes. To better understand the role of glutathione-S-transferases (GSTs) in detoxification mechanism in M. pulchricornis, we identified seventeen GST genes (MpulGSTs) from the body transcriptome. Among these, fifteen MpulGSTs belonged to cytosolic GSTs and the other two belonged to microsomal classes. The cytosolic GSTs were classified into four different clades: four in delta, three in omega, seven in sigma, and one in zeta. The expression levels of these MpulGSTs after exposure to sub-lethal concentrations of phoxim and cypermethrin were determined using real-time quantitative polymerase chain reaction: seven MpulGSTs (MpulGSTD3, MpulGSTS1, MpulGSTS2, MpulGSTS4, MpulGSTS6 MpulGSTO3, and MpulGSTmic1) and 11 MpulGSTs (MpulGSTD1, MpulGSTD2, MpulGSTD3, MpulGSTO2, MpulGSTS1, MpulGSTS2, MpulGSTS3, MpulGSTS4, MpulGSTS5, MpulGSTS7, and MpulGSTmic1) were highly expressed, respectively. These results suggested that GST genes may play a pivotal role in detoxification process in M. pulchricornis. Our findings would provide a theoretical base for elucidating insecticide susceptibility and should promote functional research on specific GST genes in parasitoid wasps.

Survey of the Bradysia odoriphaga Transcriptome Using PacBio Single-Molecule Long-Read Sequencing

The damage caused by Bradysia odoriphaga is the main factor threatening the production of vegetables in the Liliaceae family. However, few genetic studies of B. odoriphaga have been conducted because of a lack of genomic resources. Many long-read sequencing technologies have been developed in the last decade; therefore, in this study, the transcriptome including all development stages of B. odoriphaga was sequenced for the first time by Pacific single-molecule long-read sequencing. Here, 39,129 isoforms were generated, and 35,645 were found to have annotation results when checked against sequences available in different databases. Overall, 18,473 isoforms were distributed in 25 various Clusters of Orthologous Groups, and 11,880 isoforms were categorized into 60 functional groups that belonged to the three main Gene Ontology classifications. Moreover, 30,610 isoforms were assigned into 44 functional categories belonging to six main Kyoto Encyclopedia of Genes and Genomes functional categories. Coding DNA sequence (CDS) prediction showed that 36,419 out of 39,129 isoforms were predicted to have CDS, and 4319 simple sequence repeats were detected in total. Finally, 266 insecticide resistance and metabolism-related isoforms were identified as candidate genes for further investigation of insecticide resistance and metabolism in B. odoriphaga.

The impact of four widely used neonicotinoid insecticides on Tetragonisca angustula (Latreille) (Hymenoptera: Apidae)

Application of neonicotinoid insecticides on crops can reduce the pollination services and population levels of the stingless bee Tetragonisca angustula (Latreille) (Hymenoptera: Apidae) in Neotropical agroecosystems. However, the impact of these insecticides on this bee has not been fully investigated. This study assessed the susceptibility levels of T. angustula to four neonicotinoid insecticides (acetamiprid, imidacloprid, thiacloprid and thiamethoxam), widely used to manage pests on different crops, and their effects on locomotion of the bee. Neonicotinoids with the cyano radical caused lower bee mortality (assessed by mean lethal concentration, LC₅₀), while those compounds with the nitro radical were highly harmful to T. angustula. Locomotion activity was strongly impaired in bees treated with thiacloprid and imidacloprid, while acetamiprid did not affect the locomotion activity, which was similar to the control. Application of thiamethoxam caused hyperactivity, as observed by increases of ∼4.5 and 5.0-fold in mean speed and distance traveled, respectively. These results suggest that applications of neonicotinoid insecticides can negatively affect the pollination activity of T. angustula. These results can also help to develop strategies to conserve these pollinators in agroecosystems.

Differential Expression of Cytochrome P450 CYP6 Genes in the Brown Marmorated Stink Bug, Halyomorpha halys (Hemiptera: Pentatomidae)

Cytochrome (CYP) P450s are a superfamily of enzymes that detoxify xenobiotics and regulate numerous physiological processes in insects. The genomes of phytophagous insects usually contain large numbers of P450s, especially within the CYP3 clan. Within this clan, CYP6 subfamily members help detoxify plant host secondary metabolites. In this study, we analyzed three CYP6 genes in the highly polyphagous invasive pest, Halyomorpha halys (Stål), commonly known as brown marmorated stink bug. We characterized and validated the expression of HhCYP6BQ27, HhCYP6BK13, and HhCYP6BK24 among sexes, tissues (gut, fat body, and Malpighian tubules) and hosts (apple, corn, soybean). Sequence characterization by amino acid alignments confirmed the presence of conserved motifs typical of the P450 superfamily. No significant differences existed in gene expression among sexes or when fed different hosts, suggesting that these transcripts might have broad substrate specificities. However, significant differences in gene expression were observed among the tissues studied and were gene-dependent. Collectively, the results show that H. halys differentially expressed CYP6 genes among tissues, which may be related to important and specific physiological functions. This study has increased our understanding of H. halys biology that can be useful for functional studies and can potentially be exploited in developing sustainable pest management strategies.

MALDI-imaging analyses of honeybee brains exposed to a neonicotinoid insecticide

BACKGROUND

Toxicological studies evaluating the possible harmful effects of pesticides on bees are important and allow the emergence of protection and pollinator conservation strategies. This study aimed to evaluate the effects of exposure to a sublethal concentration of imidacloprid (LC_50/100 : 0.014651 ng imidacloprid µL^-1 diet) on the distribution of certain proteins identified in the brain of Apis mellifera worker bees using a MALDI-imaging approach. This technique enables proteomic analysis of tissues in situ by monitoring the spatiotemporal dynamics of the biochemical processes occurring at a specific time in specific brain neuropils. For this purpose, foraging bees were exposed to an 8-day diet containing a sublethal concentration of imidacloprid corresponding to the LC_50/100 . Bees were collected on day 8 of exposure, and their brains analyzed using protein density maps.

RESULTS

The results showed that exposure to imidacloprid led to a series of biochemical changes, including alterations in synapse regulation, apoptosis regulation and oxidative stress, which may adversely impair the physiology of these colony bees.

CONCLUSION

Worker bee contact with even tiny amounts of imidacloprid had potent effects leading to the overexpression of a series of proteins related to important cellular processes that were possibly damaged by the insecticide. © 2018 Society of Chemical Industry.

Laboratory bioassays on the response of honey bee (Apis mellifera L.) glutathione S-transferase and acetylcholinesterase to the oral exposure to copper, cadmium, and lead

In the present study, the influence of cadmium, copper, and lead on two enzymes often used as biomarkers in toxicological analysis was investigated. Bees were fed with 1 M sucrose solution containing 10-fold serial dilutions of CuCl₂ (1000 mg L^-1, 100 mg L^-1, and 10 mg L^-1), CdCl₂ (0.1 mg L^-1, 0.01 mg L^-1, and 0.001 mg L^-1), or PbCl₂ (10 mg L^-1, 1 mg L^-1, and 0.1 mg L^-1) during 48 h. Our results showed that the total glutathione S-transferase activity was not changed under the influence of cadmium and lead, and it was decreased with the highest concentration of copper. The level of gene expression of the three analyzed classes of glutathione S-transferase was significantly increased with increasing concentrations of copper and cadmium. Lead did not cause significant changes in glutathione S-transferase activity and gene expression, while it showed biphasic effect on acetylcholinesterase activity: lower concentration of lead, 0.1 mg L^-1 inhibited and higher dose, 10 mg L^-1 induced acetylcholinesterase activity in honey bees. Furthermore, our results showed a significant decrease of the acetylcholinesterase activity in honey bees treated with 0.001 and 0.01 mg L^-1 CdCl₂. Our results indicate the influence of cadmium, copper, and lead on GST and AChE in the honey bees. These results form the basis for future research on the impact of metallic trace element pollution on honey bees.

Detoxifying enzyme complements and host use phenotypes in 160 insect species

We use the genomes of 160 insect species to test the hypothesis that the size of detoxifying enzyme families is greater in species using more chemically diverse food resources. Phylogenetically appropriate contrasts in subsamples of the data generally support the hypothesis. We find relatively high numbers of cytochrome P450, glutathione S-transferase and carboxyl/choline esterase genes in omnivores and herbivores feeding on chemically complex tissues and relatively low numbers of these genes in specialists on relatively simple diets, including plant sap, nectar and pollen, and blood. Among Lepidoptera feeding on green plant tissue and Condylognatha feeding on sap we also find more of these genes in highly polyphagous species, many of which are major agricultural pests. These genomic signatures of food resource use are consistent with the hypothesis that some taxa are preadapted for insecticide resistance evolution.

Genomic insights into neonicotinoid sensitivity in the solitary bee Osmia bicornis

The impact of pesticides on the health of bee pollinators is determined in part by the capacity of bee detoxification systems to convert these compounds to less toxic forms. For example, recent work has shown that cytochrome P450s of the CYP9Q subfamily are critically important in defining the sensitivity of honey bees and bumblebees to pesticides, including neonicotinoid insecticides. However, it is currently unclear if solitary bees have functional equivalents of these enzymes with potentially serious implications in relation to their capacity to metabolise certain insecticides. To address this question, we sequenced the genome of the red mason bee, Osmia bicornis, the most abundant and economically important solitary bee species in Central Europe. We show that O. bicornis lacks the CYP9Q subfamily of P450s but, despite this, exhibits low acute toxicity to the N-cyanoamidine neonicotinoid thiacloprid. Functional studies revealed that variation in the sensitivity of O. bicornis to N-cyanoamidine and N-nitroguanidine neonicotinoids does not reside in differences in their affinity for the nicotinic acetylcholine receptor or speed of cuticular penetration. Rather, a P450 within the CYP9BU subfamily, with recent shared ancestry to the Apidae CYP9Q subfamily, metabolises thiacloprid in vitro and confers tolerance in vivo. Our data reveal conserved detoxification pathways in model solitary and eusocial bees despite key differences in the evolution of specific pesticide-metabolising enzymes in the two species groups. The discovery that P450 enzymes of solitary bees can act as metabolic defence systems against certain pesticides can be leveraged to avoid negative pesticide impacts on these important pollinators.

RNA-sequencing of a citrus bud-feeder, Podagricomela weisei (Coleoptera: Chrysomelidae), reveals xenobiotic metabolism/core RNAi machinery-associated genes and conserved miRNAs

The citrus leaf-mining beetle, Podagricomela weisei Heikertinger, is an important citrus pest that ingests the mesophyll and new shoots. The mechanism underlying the xenobiotic metabolism of P. weisei is not well understood, in part because of a lack of available genomic and transcriptomic data, which has hampered the development of novel pest management approaches [e.g., RNA interference (RNAi)]. In this study, we completed the deep sequencing of the P. weisei transcriptome to identify factors potentially involved in xenobiotic metabolism and the core RNAi machinery. The sequencing of the P. weisei transcriptome generated >27 million clean reads, ultimately yielding 90,410 unigenes with an N50 of 1065 bp. The unigenes were used as queries to search the Nr database, which revealed that 21,847 unigenes were homologous to known genes in various species. Transcripts encoding genes involved in xenobiotic metabolism were identified, including genes encoding cytochrome P450 monooxygenase (P450, 47 unigenes), glutathione S-transferase (GST, 12 unigenes), esterase (EST, 25 unigenes), and the ATP-binding cassette transporter (ABC transporter, 32 unigenes). A parallel sequencing of small RNAs detected 30 conserved miRNAs, with the most abundant being Pwe-miR-1-3p, with an expression level reaching 517,996 reads in the prepared library, followed by Pwe-miR-8-3p (149,402 reads). Genes encoding components of the miRNA, siRNA, and piRNA pathways were also identified, and the results indicated that P. weisei possesses only one of each gene in all three pathways. In summary, this is the first detailed analysis of the transcriptome and small RNAs of P. weisei. The datasets presented herein may form the basis for future molecular characterizations of P. weisei as well as the development of enhanced pest control strategies.

Combined Toxicity of Insecticides and Fungicides Applied to California Almond Orchards to Honey Bee Larvae and Adults

Beekeepers providing pollination services for California almond orchards have reported observing dead or malformed brood during and immediately after almond bloom—effects that they attribute to pesticide exposure. The objective of this study was to test commonly used insecticides and fungicides during almond bloom on honey bee larval development in a laboratory bioassay. In vitro rearing of worker honey bee larvae was performed to test the effect of three insecticides (chlorantraniliprole, diflubenzuron, and methoxyfenozide) and three fungicides (propiconazole, iprodione, and a mixture of boscalid-pyraclostrobin), applied alone or in insecticide-fungicide combinations, on larval development. Young worker larvae were fed diets contaminated with active ingredients at concentration ratios simulating a tank-mix at the maximum label rate. Overall, larvae receiving insecticide and insecticide-fungicide combinations were less likely to survive to adulthood when compared to the control or fungicide-only treatments. The insecticide chlorantraniliprole increased larval mortality when combined with the fungicides propiconazole or iprodione, but not alone; the chlorantraniliprole-propiconazole combination was also found to be highly toxic to adult workers treated topically. Diflubenzuron generally increased larval mortality, but no synergistic effect was observed when combined with fungicides. Neither methoxyfenozide nor any methoxyfenozide-fungicide combination increased mortality. Exposure to insecticides applied during almond bloom has the potential to harm honey bees and this effect may, in certain instances, be more damaging when insecticides are applied in combination with fungicides.

Influence of microbiota in the susceptibility of parasitic wasps to abamectin insecticide: deep sequencing, esterase and toxicity tests

BACKGROUND

The parasitic wasps Eretmocerus mundus, Eretmocerus eremicus and Encarsia formosa are important natural enemies of whiteflies. A broad understanding of their biology, ecology and behavior has been achieved, but the composition and role of their microbiota is not fully determined. The knowledge of the bacteria present in insects might be useful to manage species of human concern such as natural enemies or pests. Here, we performed a residual contact test to study a possible change in the susceptibility of E. mundus adults to abamectin insecticide after antibiotic treatment. Moreover, we assessed the microbiota present in adults of E. eremicus, E. formosa and two strains of E. mundus by MiSeq 16S rRNA amplicon sequencing. Finally, enzymatic tests were done to determine the influence of Arthrobacter species in the susceptibility of E. mundus to pesticides.

RESULTS

The assays showed that when E. mundus adults were pretreated with antibiotic, the toxicity of abamectin was significantly higher. Among the different bacteria associated with parasitic wasps, Arthrobacter has been shown to be involved in the degradation of several kinds of pesticides. Four Arthrobacter species were detected in all the studied insects and the presence of esterases in this bacterial species was confirmed.

CONCLUSIONS

The results suggest that the microbiota can modify the susceptibility of E. mundus to pesticides, which in turn supports the importance of the microbial community in natural enemies that it should be considered as a factor in risk assessment tests of pesticides. © 2018 Society of Chemical Industry.

Analysis of Dendrobium huoshanense transcriptome unveils putative genes associated with active ingredients synthesis

BACKGROUND

Dendrobium huoshanense C.Z. Tang et S.J. Cheng is a traditional Chinese herbal medicine with high medicinal value in China. Polysaccharides and alkaloids are its main active ingredients. To understand the difference of main active ingredients in different tissues, we determined the contents of polysaccharides and alkaloids in the roots, stems and leaves of D. huoshanense. In order to explore the reasons for the differences of active ingredients at the level of transcription, we selected roots, stems and leaves of D. huoshanenese for transcriptome sequencing and pathway mining.

RESULTS

The contents of polysaccharides and alkaloids of D. huoshanense were determined and it was found that there were significant differences in different tissues. A total of 716,634,006 clean reads were obtained and 478,361 unigenes were assembled by the Illumina platform sequencing. We identified 1407 carbohydrate-active related unigenes against CAZy database including 447 glycosyltransferase genes (GTs), 818 glycoside hydrolases (GHs), 60 carbohydrate esterases (CEs), 62 carbohydrate-binding modules (CBMs), and 20 polysaccharide lyases (PLs). In the glycosyltransferases (GTs) family, 315 differential expression genes (DEGs) were identified. In total, 124 and 58 DEGs were associated with the biosynthesis of alkaloids in Dh_L vs. Dh_S and Dh_R vs. Dh_L, respectively. A total of 62 DEGs associated with the terpenoid pathway were identified between Dh_R and Dh_S. Five key enzyme genes involved in the terpenoids pathway were identified, and their expression patterns in different tissues was validated using quantitative real-time PCR.

CONCLUSIONS

In summary, our study presents a transcriptome profile of D. huoshanense. These data contribute to our deeper relevant researches on active ingredients and provide useful insights into the molecular mechanisms regulating polysaccharides and alkaloids in Dendrobium.

Analysis of Dendrobium huoshanense transcriptome unveils putative genes associated with active ingredients synthesis

Background

Dendrobium huoshanense C.Z. Tang et S.J. Cheng is a traditional Chinese herbal medicine with high medicinal value in China. Polysaccharides and alkaloids are its main active ingredients. To understand the difference of main active ingredients in different tissues, we determined the contents of polysaccharides and alkaloids in the roots, stems and leaves of D. huoshanense. In order to explore the reasons for the differences of active ingredients at the level of transcription, we selected roots, stems and leaves of D. huoshanenese for transcriptome sequencing and pathway mining.

Results

The contents of polysaccharides and alkaloids of D. huoshanense were determined and it was found that there were significant differences in different tissues. A total of 716,634,006 clean reads were obtained and 478,361 unigenes were assembled by the Illumina platform sequencing. We identified 1407 carbohydrate-active related unigenes against CAZy database including 447 glycosyltransferase genes (GTs), 818 glycoside hydrolases (GHs), 60 carbohydrate esterases (CEs), 62 carbohydrate-binding modules (CBMs), and 20 polysaccharide lyases (PLs). In the glycosyltransferases (GTs) family, 315 differential expression genes (DEGs) were identified. In total, 124 and 58 DEGs were associated with the biosynthesis of alkaloids in Dh_L vs. Dh_S and Dh_R vs. Dh_L, respectively. A total of 62 DEGs associated with the terpenoid pathway were identified between Dh_R and Dh_S. Five key enzyme genes involved in the terpenoids pathway were identified, and their expression patterns in different tissues was validated using quantitative real-time PCR.

Conclusions

In summary, our study presents a transcriptome profile of D. huoshanense. These data contribute to our deeper relevant researches on active ingredients and provide useful insights into the molecular mechanisms regulating polysaccharides and alkaloids in Dendrobium.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5305-6) contains supplementary material, which is available to authorized users.

Identification of cytochrome P450 monooxygenase genes and their expression in response to high temperature in the alligatorweed flea beetle Agasicles hygrophila (Coleoptera: Chrysomelidae)

Cytochrome P450 monooxygenases (P450s) are a large class of enzymes that play essential roles in metabolic processes such as hormone synthesis and the catabolism of toxins and other chemicals in insects. In the present study, we identified 82 P450 genes using comprehensive RNA sequencing in the flea beetle Agasicles hygrophila, and all of the sequences were validated by cloning and sequencing. Phylogenetic analysis showed that the P450 genes in A. hygrophila fell into the mitochondrial clan, CYP2 clan, CYP3 clan and CYP4 clan and were classified into 20 families and 48 subfamilies. Most A. hygrophila P450 genes had high sequence homology with those from other coleopteran insects. To understand the effects of high temperatures on the metabolic processes of female and male adults, we studied the effects of two temperature regimes (constant temperature of 28 °C for 20 h with a 4-h period of high temperatures of 30 °C and 39 °C) on the expression levels of P450 genes in A. hygrophila using RT-PCR and qRT-PCR. The results showed that there were no differences in expression in 30 P450 genes between the control and high-temperature-treated A. hygrophila adults, while 22 P450 genes showed up-regulated expression and 19 P450 genes were down-regulated in A. hygrophila female adults after high-temperature treatment. For A. hygrophila male adults exposed to high temperatures, we found that 8 P450 genes had higher expression levels and 12 P450 genes had lower expression levels under the same conditions. The P450 genes are candidates that showed significantly different expression levels after high-temperature treatments in A. hygrophila adults, and further studies are needed to determine their possible roles in metabolic processes during the response to elevated temperatures.

Molecular evolution of juvenile hormone esterase-like proteins in a socially exchanged fluid

Socially exchanged fluids are a direct means by which an organism can influence conspecifics. It was recently shown that when workers of the carpenter ant Camponotus floridanus feed larval offspring via trophallaxis, they transfer Juvenile Hormone III (JH), a key developmental regulator, as well as paralogs of JH esterase (JHE), an enzyme that catalyzes the hydrolysis of JH. Here we combine proteomic, phylogenetic and selection analyses to investigate the evolution of this esterase subfamily. We show that Camponotus JHE-like proteins have undergone multiple duplications, experienced positive selection, and changed tissue localization to become abundantly and selectively present in trophallactic fluid. The Camponotus trophallactic esterases have maintained their catalytic triads and contain a number of positively-selected amino acid changes distributed throughout the protein, which possibly reflect an adaptation to the highly acidic trophallactic fluid of formicine ants. To determine whether these esterases might regulate larval development, we fed workers with a JHE-specific pharmacological inhibitor to introduce it into the trophallactic network. This inhibitor increased the likelihood of pupation of the larvae reared by these workers, similar to the influence of food supplementation with JH. Together, these findings suggest that JHE-like proteins have evolved a new role in the inter-individual regulation of larval development in the Camponotus genus.

Structure and expression of sulfatase and sulfatase modifying factor genes in the diamondback moth, Plutella xylostella

The diamondback moth, Plutella xylostella (L.), uses sulfatases (SULF) to counteract the glucosinolate-myrosinase defensive system that cruciferous plants have evolved to deter insect feeding. Sulfatase activity is regulated by post-translational modification of a cysteine residue by sulfatase modifying factor 1 (SUMF1). We identified 12 SULF genes (PxylSulfs) and two SUMF1 genes (PxylSumf1s) in the P. xylostella genome. Phylogenetic analysis of SULFs and SUMFs from P. xylostella, Bombyx mori, Manduca sexta, Heliconius melpomene, Danaus plexippus, Drosophila melanogaster, Tetranychus urticae and Homo sapiens showed that the SULFs were clustered into five groups, and the SUMFs could be divided into two groups. Profiling of the expression of PxylSulfs and PxylSumfs by RNA-seq and by quantitative real-time polymerase chain reaction showed that two glucosinolate sulfatase genes (GSS), PxylSulf2 and PxylSulf3, were primarily expressed in the midgut of 3rd- and 4th-instar larvae. Moreover, expression of sulfatases PxylSulf2, PxylSulf3 and PxylSulf4 were correlated with expression of the sulfatases modifying factor PxylSumf1a. The findings from this study provide new insights into the structure and expression of SUMF1 and PxylSulf genes that are considered to be key factors for the evolutionary success of P. xylostella as a specialist herbivore of cruciferous plants.

Genome of the small hive beetle (Aethina tumida, Coleoptera: Nitidulidae), a worldwide parasite of social bee colonies, provides insights into detoxification and herbivory

Background

The small hive beetle (Aethina tumida; ATUMI) is an invasive parasite of bee colonies. ATUMI feeds on both fruits and bee nest products, facilitating its spread and increasing its impact on honey bees and other pollinators. We have sequenced and annotated the ATUMI genome, providing the first genomic resources for this species and for the Nitidulidae, a beetle family that is closely related to the extraordinarily species-rich clade of beetles known as the Phytophaga. ATUMI thus provides a contrasting view as a neighbor for one of the most successful known animal groups.

Results

We present a robust genome assembly and a gene set possessing 97.5% of the core proteins known from the holometabolous insects. The ATUMI genome encodes fewer enzymes for plant digestion than the genomes of wood-feeding beetles but nonetheless shows signs of broad metabolic plasticity. Gustatory receptors are few in number compared to other beetles, especially receptors with known sensitivity (in other beetles) to bitter substances. In contrast, several gene families implicated in detoxification of insecticides and adaptation to diverse dietary resources show increased copy numbers. The presence and diversity of homologs involved in detoxification differ substantially from the bee hosts of ATUMI.

Conclusions

Our results provide new insights into the genomic basis for local adaption and invasiveness in ATUMI and a blueprint for control strategies that target this pest without harming their honey bee hosts. A minimal set of gustatory receptors is consistent with the observation that, once a host colony is invaded, food resources are predictable. Unique detoxification pathways and pathway members can help identify which treatments might control this species even in the presence of honey bees, which are notoriously sensitive to pesticides.

Lethality of synthetic and natural acaricides to worker honey bees (Apis mellifera) and their impact on the expression of health and detoxification-related genes

In this study, honey bees (Apis mellifera L.) were exposed to LD₀₅ and LD₅₀ doses of five commonly used acaricides for controlling the parasitic mite, Varroa destructor. LD₅₀ values at 48 h post-treatment showed that tau-fluvalinate was the most toxic, followed by amitraz, coumaphos, thymol, and formic acid. However, the hazard ratios, which estimate the hive risk level based on a ratio of a standard dose of acaricide per hive to the LD₅₀ of the acaricide, revealed that tau-fluvalinate was the most hazardous followed by formic acid, coumaphos, amitraz, and thymol. The expression of the honey bee acetylcholinesterase gene increased after treatment with the LD₀₅ and LD₅₀ acaricide doses and could distinguish three patterns in the timing and level of increased expression between acaricides: one for amitraz, one for tau-fluvalinate and formic acid, and one for coumaphos and thymol. Conversely, changes in cytochrome P450 gene expression could also be detected in response to all five acaricides, but there were no significant differences between them. Changes in vitellogenin gene expression could only detect the effects of tau-fluvalinate, amitraz, or coumaphos treatment, which were not significantly different from each other. Among the acaricides tested, coumaphos, amitraz, and thymol appear to be the safest acaricides based on their hazard ratios, and a good marker to detect differences between the effects of sub-lethal doses of acaricides is monitoring changes in acetylcholinesterase gene expression.

Modified neonicotinoid insecticide with bi-directional selective toxicity and drug resistance

A three-dimensional quantitative structure-activity relationship (3D-QSAR) model was established based on the molecular structures and the negative logarithm of experimental lethal concentration 50 values (pLC₅₀) of neonicotinoid insecticides. Then, the mechanisms of bi-directional selective toxic effects and drug resistance were determined using homology modeling and molecular docking analyses. The results of the model showed that the 1-, 2-, 4-, and 12- positions of neonicotinoid insecticides strongly affected their toxicity, and that the introduction of bulky or electropositive groups at these positions could increase the pLC₅₀ values. Using Compound 19 as a template, we designed 37 derivatives with greater toxicity (increased by 0.04-11.45%). Among them, 20 derivatives had bioconcentrations lower than that of Compound 19 (reduced by 0.38-147.88%). Further screening of Compound 19 and the 20 derivatives mentioned above by homology modeling and acetylcholine receptors (AChRs) molecular docking analyses showed that 10 derivatives had bi-directional selective toxic effects against pests and bees. Further docking analyses of Compound 19 and these 10 derivatives identified that Derivative-33 showed decreased docking with superoxide dismutase (SOD) and glutathione S transferase (GST) in pests and enhanced docking with these enzymes in bees, indicating bi-directional selective resistance for pests and bees. Accordingly, Derivative-33 was selected as a new insecticide with high toxicity to pests and low toxicity to bees (bi-directional selective toxicity), low resistance in pest populations, and high resistance in bee populations. This study provides valuable reference data and will be useful for the development of strategies to produce new environmentally friendly pesticides.

CYP6AE gene cluster knockout in Helicoverpa armigera reveals role in detoxification of phytochemicals and insecticides

The cotton bollworm Helicoverpa armigera, is one of the world's major pest of agriculture, feeding on over 300 hosts in 68 plant families. Resistance cases to most insecticide classes have been reported for this insect. Management of this pest in agroecosystems relies on a better understanding of how it copes with phytochemical or synthetic toxins. We have used genome editing to knock out a cluster of nine P450 genes and show that this significantly reduces the survival rate of the insect when exposed to two classes of host plant chemicals and two classes of insecticides. Functional expression of all members of this gene cluster identified the P450 enzymes capable of metabolism of these xenobiotics. The CRISPR-Cas9-based reverse genetics approach in conjunction with in vitro metabolism can rapidly identify the contributions of insect P450s in xenobiotic detoxification and serve to identify candidate genes for insecticide resistance.

Identification and Characterization of Three New Cytochrome P450 Genes and the Use of RNA Interference to Evaluate Their Roles in Antioxidant Defense in Apis cerana cerana Fabricius

Cytochrome P450s play critical roles in maintaining redox homeostasis and protecting organisms from the accumulation of toxic reactive oxygen species (ROS). The biochemical functions of the P450 family have essentially been associated with the metabolism of xenobiotics. Here, we sequenced and characterized three P450 genes, AccCYP314A1, AccCYP4AZ1, and AccCYP6AS5, from Apis cerana cerana Fabricius; these genes play a critical role in maintaining biodiversity. Quantitative PCR (qPCR) analysis indicated that the three genes were all predominantly expressed in the epidermis (EP), followed by the brain (BR) and midgut (MG). In addition, the highest expression levels were detected in the dark-eyed pupae and adult stages. The three genes were induced by temperature (4°C and 44°C), heavy metals (CdCl₂ and HgCl₂), pesticides (DDV, deltamethrin, and paraquat) and UV treatments. Furthermore, Western blot analysis indicated that the protein expression levels could be induced by some abiotic stressors, a result that complements the qPCR results. We analyzed the silencing of these three genes and found that silencing these genes enhanced the enzymatic activities of peroxidase (POD) and catalase (CAT). Additionally, we investigated the expression of other antioxidant genes and found that some were upregulated, while others were downregulated, suggesting that the upregulated genes may be involved in compensating for the silencing of AccCYP314A1, AccCYP4AZ1, and AccCYP6AS5. Our findings suggest that AccCYP314A1, AccCYP4AZ1, and AccCYP6AS5 may play very significant roles in the antioxidant defense against damage caused by ROS.

Glyphosate affects the larval development of honey bees depending on the susceptibility of colonies

As the main agricultural insect pollinator, the honey bee (Apis mellifera) is exposed to a number of agrochemicals, including glyphosate (GLY), the most widely used herbicide. Actually, GLY has been detected in honey and bee pollen baskets. However, its impact on the honey bee brood is poorly explored. Therefore, we assessed the effects of GLY on larval development under chronic exposure during in vitro rearing. Even though this procedure does not account for social compensatory mechanisms such as brood care by adult workers, it allows us to control the herbicide dose, homogenize nutrition and minimize environmental stress. Our results show that brood fed with food containing GLY traces (1.25-5.0 mg per litre of food) had a higher proportion of larvae with delayed moulting and reduced weight. Our assessment also indicates a non-monotonic dose-response and variability in the effects among colonies. Differences in genetic diversity could explain the variation in susceptibility to GLY. Accordingly, the transcription of immune/detoxifying genes in the guts of larvae exposed to GLY was variably regulated among the colonies studied. Consequently, under laboratory conditions, the response of honey bees to GLY indicates that it is a stressor that affects larval development depending on individual and colony susceptibility.

Functional characterization of CYP4G11-a highly conserved enzyme in the western honey bee Apis mellifera

Determining the functionality of CYP4G11, the only CYP4G in the genome of the western honey bee Apis mellifera, can provide insight into its reduced CYP4 inventory. Toward this objective, CYP4G11 transcripts were quantified, and CYP4G11 was expressed as a fusion protein with housefly CPR in Sf9 cells. Transcript levels varied with age, task, and tissue type in a manner consistent with the need for cuticular hydrocarbon production to prevent desiccation or with comb wax production. Young larvae, with minimal need for desiccation protection, expressed CYP4G11 at very low levels. Higher levels were observed in nurses, and even higher levels in wax producers and foragers, the latter of which risk desiccation upon leaving the hive. Recombinant CYP4G11 readily converted octadecanal to n-heptadecane in a time-dependent manner, demonstrating its functions as an oxidative decarbonylase. CYP4G11 expression levels are high in antennae; heterologously expressed CYP4G11 converted tetradecanal to n-tridecane, demonstrating that it metabolizes shorter-chain aldehydes. Together, these findings confirm the involvement of CYP4G11 in cuticular hydrocarbon production and suggest a possible role in clearing pheromonal and phytochemical compounds from antennae. This possible dual functionality of CYP4G11, i.e., cuticular hydrocarbon and comb wax production and antennal odorant clearance, may explain how honey bees function with a reduced CYP4G inventory.

Pesticide residues in honey bees, pollen and beeswax: Assessing beehive exposure

In order to study the distribution of pesticide residues in beekeeping matrices, samples of live in-hive worker honey bees (Apis mellifera), fresh stored pollen and beeswax were collected during 2016-2017 from 45 apiaries located in different landscape contexts in Spain. A total of 133 samples were screened for 63 pesticides or their degradation products to estimate the pesticide exposure to honey bee health through the calculation of the hazard quotient (HQ). The influence of the surrounding environment on the content of pesticides in pollen was assessed by comparing the concentrations of pesticide residues found in apiaries from intensive farming landscapes to those found in apiaries located in mountainous, grassland and urban contexts. Beeswax revealed high levels of miticides used in beekeeping such as coumaphos, chlorfenvinphos, fluvalinate and acrinathrin, which were detected in more than 75% of samples. Pollen was predominantly contaminated by miticides but also by insecticides used in agriculture such as chlorpyrifos and acetamiprid, which showed concentrations significantly higher in apiaries located in intensive farming contexts. Pesticides residues were less frequent and at lower concentrations in live honey bees. Beeswax showed the highest average hazard scores (HQ > 5000) to honey bees. Pollen samples contained the largest number of pesticide residues and relevant hazard (HQ > 50) to bees. Acrinathrin was the most important contributor to the hazard quotient scores in wax and pollen samples. The contributions of the pesticides dimethoate and chlorpyrifos to HQ were considered relevant in samples.

Identification and expression profiles of twenty-six glutathione S-transferase genes from rice weevil, Sitophilus oryzae (Coleoptera: Curculionidae)

The rice weevil, Sitophilus oryzae, is one of the most destructive pests in stored cereals products. In this study, 26 cDNAs encoding glutathione S-transferases (GSTs) were sequenced and characterized in S. oryzae. Phylogenetic analysis displayed the categorization of 26 GSTs into six different cytosolic classes, including two in the delta, twelve in epsilon, three in omega, six in sigma, two in theta, and one in zeta class. RT-qPCR assay illustrated that the relative expression of ten GST genes was significantly higher in adult stages than in larval and pupal developmental stages. Tissue-specific expression analysis revealed that the SoGSTe5, SoGSTe7, SoGSTe12, and SoGSTz1 were up-regulated in the midgut, SoGSTe2, SoGSTe6, and SoGSTs2 were up-regulated in the fat body, and three GSTs (SoGSTd1, SoGSTd2 and SoGSTe4) were up-regulated in Malpighian tubules. RT-qPCR indicated that five GST genes were over expressed after exposure to phosphine at various times and concentrations. The increase in GST gene expressions after phosphine exposure in S. oryzae may lead to an improved tolerance for fumigations and xenobiotics.

Responses of Honey Bees to Lethal and Sublethal Doses of Formulated Clothianidin Alone and Mixtures

The widespread use of neonicotinoid insecticides has sparked concern over the toxicity risk to honey bees (Apis mellifera L. (Hymenoptera: Apidae)). In this study, feeding treatments with the clothianidin formulation at 2.6 ppb (residue concentration) or its binary mixtures with five representative pesticides (classes) did not influence on adult survivorship, but all treatments caused significantly lower body weight than controls. Most binary mixtures at residue levels showed minor or no interaction on body weight loss, and synergistic interaction was detected only from the mixture of clothianidin + λ-cyhalothrin. Chlorpyrifos alone and the mixture of clothianidin + chlorpyrifos significantly suppressed esterase (EST) activity, while most treatments of individual pesticides and mixtures had no effect on EST and glutathione S-transferase (GST) activities. However, ingestion of clothianidin at 2.6 ppb significantly enhanced P450 oxidase activity by 19%. The LC50 of formulated clothianidin was estimated at 0.53 ppm active ingredient, which is equivalent to 25.4 ng clothianidin per bee (LD50) based on the average sugar consumption of 24 µl per bee per day. In addition to mortality, ingestion of clothianidin at LC50 significantly reduced bee body weight by 12%. P450 activities were also significantly induced at 24 and 48 h in clothianidin-treated bees, while no significant difference was found in GST and EST activities. Further examinations revealed that the expression of an important CYP9q1 detoxification gene was significantly induced by clothianidin. Thus, data consistently indicated that P450s were involved in clothianidin detoxification in honey bees. Although the honey bee population in Stoneville (MS, United States) had sixfold lower susceptibility than other reported populations, clothianidin had very high oral toxicity to bees.

Genome-wide and expression-profiling analyses suggest the main cytochrome P450 genes related to pyrethroid resistance in the malaria vector, Anopheles sinensis (Diptera Culicidae)

BACKGROUND

Anopheles sinensis is one of the major malaria vectors. However, pyrethroid resistance in An. sinensis is threatening malaria control. Cytochrome P450-mediated detoxification is an important pyrethroid resistance mechanism that has been unexplored in An. sinensis. In this study, we performed a comprehensive analysis of the An. sinensis P450 gene superfamily with special attention to their role in pyrethroid resistance using bioinformatics and molecular approaches.

RESULTS

Our data revealed the presence of 112 individual P450 genes in An. sinensis, which were classified into four major clans (mitochondrial, CYP2, CYP3 and CYP4), 18 families and 50 subfamilies. Sixty-seven genes formed nine gene clusters, and genes within the same cluster and the same gene family had a similar gene structure. Phylogenetic analysis showed that most of An. sinensis P450s (82/112) had very close 1: 1 orthology with Anopheles gambiae P450s. Five genes (AsCYP6Z2, AsCYP6P3v1, AsCYP6P3v2, AsCYP9J5 and AsCYP306A1) were significantly upregulated in three pyrethroid-resistant populations in both RNA-seq and RT-qPCR analyses, suggesting that they could be the most important P450 genes involved in pyrethroid resistance in An. sinensis.

CONCLUSION

Our study provides insight on the diversity of An. sinensis P450 superfamily and basis for further elucidating pyrethroid resistance mechanism in this mosquito species. © 2018 Society of Chemical Industry.

Enhancement of chronic bee paralysis virus levels in honeybees acute exposed to imidacloprid: A Chinese case study

Though honeybee populations have not yet been reported to be largely lost in China, many stressors that affect the health of honeybees have been confirmed. Honeybees inevitably come into contact with environmental stressors that are not intended to target honeybees, such as pesticides. Although large-scale losses of honeybee colonies are thought to be associated with viruses, these viruses usually lead to covert infections and to not cause acute damage if the bees do not encounter outside stressors. To reveal the potential relationship between acute pesticides and viruses, we applied different doses of imidacloprid to adult bees that were primarily infected with low levels (4.3×10⁵ genome copies) of chronic bee paralysis virus (CBPV) to observe whether the acute oral toxicity of imidacloprid was able to elevate the level of CBPV. Here, we found that the titer of CBPV was significantly elevated in adult bees after 96h of acute treatment with imidacloprid at the highest dose 66.9ng/bee compared with other treatments and controls. Our study provides clear evidence that exposure to acute high doses of imidacloprid in honeybees persistently infected by CBPV can exert a remarkably negative effect on honeybee survival. These results imply that acute environmental stressors might be one of the major accelerators causing rapid viral replication, which may progress to cause mass proliferation and dissemination and lead to colony decline. The present study will be useful for better understanding the harm caused by this pesticide, especially regarding how honeybee tolerance to the viral infection might be altered by acute pesticide exposure.

Transfer of the Active Ingredients of Some Plant Protection Products from Raspberry Plants to Beehives

Plant protection products (PPPs) have been found increasingly in the environment. They pose a huge threat to bees, contributing to honeybee colony losses and consequently to enormous economic losses. Therefore, this field investigation was designed to determine whether their active ingredients (AIs) were transferred from raspberry plants to beehives located in the immediate neighbourhood of the crop and to what extent they were transferred. Every week for 2 months, samples of soil, raspberry leaves, flowers and fruits, worker bees, honeybee brood, and honey were collected and analysed for the presence of propyzamide, chlorpyrifos, iprodione, pyraclostrobin, boscalid, cypermethrin, difenoconazole, azoxystrobin, and pyrimethanil residues. Five of these substances were found in the worker bee bodies. Chlorpyrifos, applied to only the soil through the irrigation system, also was detected in the brood. A small amount of boscalid was noted in the honey, but its residues did not exceed the maximum residue level. For chlorpyrifos, boscalid, and pyrimethanil, a positive correlation between the occurrence of PPPs in the crops and the beehives was found. Statistical methods confirmed that the application of PPPs on a raspberry plantation, as an example of nectar-secreting plants, was linked to the transfer of their AIs to beehives.

Acute Toxicity of Permethrin, Deltamethrin, and Etofenprox to the Alfalfa Leafcutting Bee

Current regulatory requirements for insecticide toxicity to nontarget insects focus on the honey bee, Apis mellifera (L.; Hymenoptera: Apidae), but this species cannot represent all insect pollinator species in terms of response to insecticides. Therefore, we characterized the toxicity of pyrethroid insecticides used for adult mosquito management (permethrin, deltamethrin, and etofenprox) on a nontarget insect, the adult alfalfa leafcutting bee, Megachile rotundata (F.; Hymenoptera: Megachilidae) in two separate studies. In the first study, the doses causing 50 and 90% mortality (LD50 and LD90, respectively) were used as endpoints and 2-d-old adult females were exposed to eight concentrations ranging from 0.0075 to 0.076 μg/bee for permethrin and etofenprox, and 0.0013-0.0075 μg/bee for deltamethrin. For the second study, respiration rates of female M. rotundata were also recorded for 2 h after bees were dosed at the LD50 values to give an indication of stress response. Results indicated a relatively similar LD50 for permethrin and etofenprox, 0.057 and 0.051 μg/bee, respectively, and a more toxic response, 0.0016 μg/bee for deltamethrin. Comparatively, female A. mellifera workers have a LD50 value of 0.024 μg/bee for permethrin and 0.015 μg/bee for etofenprox indicating that female M. rotundata are less susceptible to topical doses of these insecticides, except for deltamethrin, where both A. mellifera and M. rotundata have an identical LD50 of 0.0016 μg/bee. Respiration rates comparing each active ingredient to control groups, as well as rates between each active ingredient, were statistically different (P < 0.0001). The addition of these results to existing information on A. mellifera may provide more insights on how other economically beneficial and nontarget bees respond to pyrethroids.

Time-to-death approach to reveal chronic and cumulative toxicity of a fungicide for honeybees not revealed with the standard ten-day test

Synthetic fungicides are pesticides widely used in agriculture to control phytopathogenic fungi. The systemicity, persistency and intense application of some of these fungicides, such as boscalid, leads to long periods of exposure for honeybees via contaminated water, pollen and nectar. We exposed adult honeybees in the lab to food contaminated with boscalid for 33 days instead of the standard 10-day test. Most of the toxic effects were observed after 10 days. The median time to death (LT₅₀) ranged from 24.9 days (lowest concentration) to 7.1 days (highest concentration) and was significantly shorter in all cases than with the control (32.0 days). The concentration and dietary doses of boscalid inducing 50% mortality (LC₅₀ and LDD₅₀, respectively) decreased strongly with the time of exposure: LC₅₀ = 14,729 and 1,174 mg/l and LDD₅₀ = 0.318 and 0.0301 mg bee⁻¹ day⁻¹ at days 8 and 25, respectively. We found evidence of reinforced toxicity when exposure is prolonged, but with an unusual pattern: no cumulative toxicity is observed until 17–18 days, when a point of inflexion appears that suggests a reduced capacity of bees to deal with the toxicant. Our results show the importance of time-to-death experiments rather than fixed-duration studies for evaluating chronic toxicity.