UDP-glycosyltransferases act as key determinants of host plant range in generalist and specialist Spodoptera species.

Phytophagous insects have evolved sophisticated detoxification systems to overcome the antiherbivore chemical defenses produced by many plants. However, how these biotransformation systems differ in generalist and specialist insect species and their role in determining insect host plant range remains an open question. Here, we show that UDP-glucosyltransferases (UGTs) play a key role in determining the host range of insect species within the Spodoptera genus. Comparative genomic analyses of Spodoptera species that differ in host plant breadth identified a relatively conserved number of UGT genes in generalist species but high levels of UGT gene pseudogenization in the specialist Spodoptera picta. CRISPR-Cas9 knockouts of the three main UGT gene clusters of Spodoptera frugiperda revealed that UGT33 genes play an important role in allowing this species to utilize the poaceous plants maize, wheat, and rice, while UGT40 genes facilitate utilization of cotton. Further functional analyses in vivo and in vitro identified the UGT SfUGT33F32 as the key mechanism that allows generalist S. frugiperda to detoxify the benzoxazinoid DIMBOA (2,4-dihydroxy-7-methoxy-2H-1,4-benzoxazin-3(4H)-one), a potent insecticidal phytotoxin produced by poaceous plants. However, while this detoxification capacity is conserved in several generalist Spodoptera species, Spodoptera picta, which specializes on Crinum plants, is unable to detoxify DIMBOA due to a nonfunctionalizing mutation in SpUGT33F34. Collectively, these findings provide insight into the role of insect UGTs in host plant adaptation, the mechanistic basis of evolutionary transitions between generalism and specialism and offer molecular targets for controlling a group of notorious insect pests.

Dual mutations in the whitefly nicotinic acetylcholine receptor ß1 subunit confer target-site resistance to multiple neonicotinoid insecticides.

Neonicotinoid insecticides, which target insect nicotinic acetylcholine receptors (nAChRs), have been widely and intensively used to control the whitefly, Bemisia tabaci, a highly damaging, globally distributed, crop pest. This has inevitably led to the emergence of populations with resistance to neonicotinoids. However, to date, there have been no reports of target-site resistance involving mutation of B. tabaci nAChR genes. Here we characterize the nAChR subunit gene family of B. tabaci and identify dual mutations (A58T&R79E) in one of these genes (BTß1) that confer resistance to multiple neonicotinoids. Transgenic D. melanogaster, where the native nAChR Dß1 was replaced with BTß1A58T&R79E, were significantly more resistant to neonicotinoids than flies where Dß1 were replaced with the wildtype BTß1 sequence, demonstrating the causal role of the mutations in resistance. The two mutations identified in this study replace two amino acids that are highly conserved in >200 insect species. Three-dimensional modelling suggests a molecular mechanism for this resistance, whereby A58T forms a hydrogen bond with the R79E side chain, which positions its negatively-charged carboxylate group to electrostatically repulse a neonicotinoid at the orthosteric site. Together these findings describe the first case of target-site resistance to neonicotinoids in B. tabaci and provide insight into the molecular determinants of neonicotinoid binding and selectivity.

Transposon accumulation at xenobiotic gene family loci in aphids.

The evolution of resistance is a major challenge for the sustainable control of pests and pathogens. Thus, a deeper understanding of the evolutionary and genomic mechanisms underpinning resistance evolution is required to safeguard health and food production. Several studies have implicated transposable elements (TEs) in xenobiotic-resistance evolution in insects. However, analyses are generally restricted to one insect species and/or one or a few xenobiotic gene families (XGFs). We examine evidence for TE accumulation at XGFs by performing a comparative genomic analysis across 20 aphid genomes, considering major subsets of XGFs involved in metabolic resistance to insecticides: cytochrome P450s, glutathione S-transferases, esterases, UDP-glucuronosyltransferases, and ABC transporters. We find that TEs are significantly enriched at XGFs compared with other genes. XGFs show similar levels of TE enrichment to those of housekeeping genes. But unlike housekeeping genes, XGFs are not constitutively expressed in germline cells, supporting the selective enrichment of TEs at XGFs rather than enrichment owing to chromatin availability. Hotspots of extreme TE enrichment occur around certain XGFs. We find, in aphids of agricultural importance, particular enrichment of TEs around cytochrome P450 genes with known functions in the detoxification of synthetic insecticides. Our results provide evidence supporting a general role for TEs as a source of genomic variation at host XGFs and highlight the existence of considerable variability in TE content across XGFs and host species. These findings show the need for detailed functional verification analyses to clarify the significance of individual TE insertions and elucidate underlying mechanisms at TE-XGF hotspots.

Manipulating multi-level selection in a fungal entomopathogen reveals social conflicts and a method for improving biocontrol traits.

Changes in parasite virulence are commonly expected to lead to trade-offs in other life history traits that can affect fitness. Understanding these trade-offs is particularly important if we want to manipulate the virulence of microbial biological control agents. Theoretically, selection across different spatial scales, i.e. between- and within-hosts, shapes these trade-offs. However, trade-offs are also dependent on parasite biology. Despite their applied importance the evolution of virulence in fungal parasites is poorly understood: virulence can be unstable in culture and commonly fails to increase in simple passage experiments. We hypothesized that manipulating selection intensity at different scales would reveal virulence trade-offs in a fungal pathogen of aphids, Akanthomyces muscarius. Starting with a genetically diverse stock we selected for speed of kill, parasite yield or infectivity by manipulating competition within and between hosts and between-populations of hosts over 7 rounds of infection. We characterized ancestral and evolved lineages by whole genome sequencing and by measuring virulence, growth rate, sporulation and fitness. While several lineages showed increases in virulence, we saw none of the trade-offs commonly found in obligately-killing parasites. Phenotypically similar lineages within treatments often shared multiple single-nucleotide variants, indicating strong convergent evolution. The most dramatic phenotypic changes were in timing of sporulation and spore production in vitro. We found that early sporulation led to reduced competitive fitness but could increase yield of spores on media, a trade-off characteristic of social conflict. Notably, the selection regime with strongest between-population competition and lowest genetic diversity produced the most consistent shift to early sporulation, as predicted by social evolution theory. Multi-level selection therefore revealed social interactions novel to fungi and showed that these biocontrol agents have the genomic flexibility to improve multiple traits-virulence and spore production-that are often in conflict in other parasites.

Symbiotic bacteria confer insecticide resistance by metabolizing buprofezin in the brown planthopper, Nilaparvata lugens (Stål).

Buprofezin, a chitin synthesis inhibitor, is widely used to control several economically important insect crop pests. However, the overuse of buprofezin has led to the evolution of resistance and exposed off-target organisms present in agri-environments to this compound. As many as six different strains of bacteria isolated from these environments have been shown to degrade buprofezin. However, whether insects can acquire these buprofezin-degrading bacteria from soil and enhance their own resistance to buprofezin remains unknown. Here we show that field strains of the brown planthopper, Nilaparvata lugens, have acquired a symbiotic bacteria, occurring naturally in soil and water, that provides them with resistance to buprofezin. We isolated a symbiotic bacterium, Serratia marcescens (Bup_Serratia), from buprofezin-resistant N. lugens and showed it has the capacity to degrade buprofezin. Buprofezin-susceptible N. lugens inoculated with Bup_Serratia became resistant to buprofezin, while antibiotic-treated N. lugens became susceptible to this insecticide, confirming the important role of Bup_Serratia in resistance. Sequencing of the Bup_Serratia genome identified a suite of candidate genes involved in the degradation of buprofezin, that were upregulated upon exposure to buprofezin. Our findings demonstrate that S. marcescens, an opportunistic pathogen of humans, can metabolize the insecticide buprofezin and form a mutualistic relationship with N. lugens to enhance host resistance to buprofezin. These results provide new insight into the mechanisms underlying insecticide resistance and the interactions between bacteria, insects and insecticides in the environment. From an applied perspective they also have implications for the control of highly damaging crop pests.

Phylogenomic and functional characterization of an evolutionary conserved cytochrome P450-based insecticide detoxification mechanism in bees.

The regulatory process for assessing the risks of pesticides to bees relies heavily on the use of the honeybee, Apis mellifera, as a model for other bee species. However, the validity of using A. mellifera as a surrogate for other Apis and non-Apis bees in pesticide risk assessment has been questioned. Related to this line of research, recent work on A. mellifera has shown that specific P450 enzymes belonging to the CYP9Q subfamily act as critically important determinants of insecticide sensitivity in this species by efficiently detoxifying certain insecticide chemotypes. However, the extent to which the presence of functional orthologs of these enzymes is conserved across the diversity of bees is unclear. Here we used a phylogenomic approach to identify > 100 putative CYP9Q functional orthologs across 75 bee species encompassing all major bee families. Functional analysis of 26 P450s from 20 representative bee species revealed that P450-mediated detoxification of certain systemic insecticides, including the neonicotinoid thiacloprid and the butenolide flupyradifurone, is conserved across all major bee pollinator families. However, our analyses also reveal that CYP9Q-related genes are not universal to all bee species, with some Megachilidae species lacking such genes. Thus, our results reveal an evolutionary conserved capacity to metabolize certain insecticides across all major bee families while identifying a small number of bee species where this function may have been lost. Furthermore, they illustrate the potential of a toxicogenomic approach to inform pesticide risk assessment for nonmanaged bee species by predicting the capability of bee pollinator species to break down synthetic insecticides.

P450 gene duplication and divergence led to the evolution of dual novel functions and insecticide cross-resistance in the brown planthopper Nilaparvata lugens.

The sustainable control of many highly damaging insect crop pests and disease vectors is threatened by the evolution of insecticide resistance. As a consequence, strategies have been developed that aim to prevent or delay resistance development by rotating or mixing insecticides with different modes of action (MoA). However, these approaches can be compromised by the emergence of mechanisms that confer cross-resistance to insecticides with different MoA. Despite the applied importance of cross-resistance, its evolutionary underpinnings remain poorly understood. Here we reveal how a single gene evolved the capacity to detoxify two structurally unrelated insecticides with different MoA. Using transgenic approaches we demonstrate that a specific variant of the cytochrome P450 CYP6ER1, previously shown to confer resistance to the neonicotinoid imidacloprid in the brown planthopper, N. lugens, also confers cross-resistance to the phenylpyrazole ethiprole. CYP6ER1 is duplicated in resistant strains, and we show that while the acquisition of mutations in two encoded substrate recognition sites (SRS) of one of the parologs led to resistance to imidacloprid, a different set of mutations, outside of known SRS, are primarily responsible for resistance to ethiprole. Epistatic interactions between these mutations and their genetic background suggest that the evolution of dual resistance from the same gene copy involved functional trade-offs in respect to CYP6ER1 catalytic activity for ethiprole versus imidacloprid. Surprisingly, the mutations leading to ethiprole and imidacloprid resistance do not confer the ability to detoxify the insecticide fipronil, another phenylpyrazole with close structural similarity to ethiprole. Taken together, these findings reveal how gene duplication and divergence can lead to the evolution of multiple novel functions from a single gene. From an applied perspective they also demonstrate how cross-resistance to structurally unrelated insecticides can evolve, and illustrate the difficulty in predicting cross-resistance profiles mediated by metabolic mechanisms.

Multiple toxins and a protease contribute to the aphid-killing ability of Pseudomonas fluorescens PpR24.

Aphids are globally important pests causing damage to a broad range of crops. Due to insecticide resistance, there is an urgent need to develop alternative control strategies. In our previous work, we found Pseudomonas fluorescens PpR24 can orally infect and kill the insecticide-resistant green-peach aphid (Myzus persicae). However, the genetic basis of the insecticidal capability of PpR24 remains unclear. Genome sequencing of PpR24 confirmed the presence of various insecticidal toxins such as Tc (toxin complexes), Rhs (rearrangement hotspot) elements, and other insect-killing proteases. Upon aphids infection with PpR24, RNA-Seq analysis revealed 193 aphid genes were differentially expressed with down-regulation of 16 detoxification genes. In addition, 1325 PpR24 genes (542 were upregulated and 783 downregulated) were subject to differential expression, including genes responsible for secondary metabolite biosynthesis, the iron-restriction response, oxidative stress resistance, and virulence factors. Single and double deletion of candidate virulence genes encoding a secreted protease (AprX) and four toxin components (two TcA-like; one TcB-like; one TcC-like insecticidal toxins) showed that all five genes contribute significantly to aphid killing, particularly AprX. This comprehensive host-pathogen transcriptomic analysis provides novel insight into the molecular basis of bacteria-mediated aphid mortality and the potential of PpR24 as an effective biocontrol agent.

Transposon-mediated insertional mutagenesis unmasks recessive insecticide resistance in the aphid Myzus persicae.

The evolution of resistance to insecticides threatens the sustainable control of many of the world's most damaging insect crop pests and disease vectors. To effectively combat resistance, it is important to understand its underlying genetic architecture, including the type and number of genetic variants affecting resistance and their interactions with each other and the environment. While significant progress has been made in characterizing the individual genes or mutations leading to resistance, our understanding of how genetic variants interact to influence its phenotypic expression remains poor. Here, we uncover a mechanism of insecticide resistance resulting from transposon-mediated insertional mutagenesis of a genetically dominant but insecticide-susceptible allele that enables the adaptive potential of a previously unavailable recessive resistance allele to be unlocked. Specifically, we identify clones of the aphid pest Myzus persicae that carry a resistant allele of the essential voltage-gated sodium channel (VGSC) gene with the recessive M918T and L1014F resistance mutations, in combination with an allele lacking these mutations but carrying a Mutator-like element transposon insertion that disrupts the coding sequence of the VGSC. This results in the down-regulation of the dominant susceptible allele and monoallelic expression of the recessive resistant allele, rendering the clones resistant to the insecticide bifenthrin. These findings are a powerful example of how transposable elements can provide a source of evolutionary potential that can be revealed by environmental and genetic perturbation, with applied implications for the control of highly damaging insect pests.

Changes in both trans- and cis-regulatory elements mediate insecticide resistance in a lepidopteron pest, Spodoptera exigua.

The evolution of insect resistance to insecticides is frequently associated with overexpression of one or more cytochrome P450 enzyme genes. Although overexpression of CYP450 genes is a well-known mechanism of insecticide resistance, the underlying regulatory mechanisms are poorly understood. Here we uncovered the mechanisms of overexpression of the P450 gene, CYP321A8 in a major pest insect, Spodoptera exigua that is resistant to multiple insecticides. CYP321A8 confers resistance to organophosphate (chlorpyrifos) and pyrethroid (cypermethrin and deltamethrin) insecticides in this insect. Constitutive upregulation of transcription factors CncC/Maf are partially responsible for upregulated expression of CYP321A8 in the resistant strain. Reporter gene assays and site-directed mutagenesis analyses demonstrated that CncC/Maf enhanced the expression of CYP321A8 by binding to specific sites in the promoter. Additional cis-regulatory elements resulting from a mutation in the CYP321A8 promoter in the resistant strain facilitates the binding of the orphan nuclear receptor, Knirps, and enhances the promoter activity. These results demonstrate that two independent mechanisms; overexpression of transcription factors and mutations in the promoter region resulting in a new cis-regulatory element that facilitates binding of the orphan nuclear receptor are involved in overexpression of CYP321A8 in insecticide-resistant S. exigua.

Aberrant splicing of a nicotinic acetylcholine receptor alpha 6 subunit is associated with spinosad tolerance in the thrips predator Orius laevigatus.

Susceptibility to insecticides is one of the limiting factors preventing wider adoption of natural enemies to control insect pest populations. Identification and selective breeding of insecticide tolerant strains of commercially used biological control agents (BCAs) is one of the approaches to overcome this constraint. Although a number of beneficial insects have been selected for increased tolerance to insecticides the molecular mechanisms underpinning these shifts in tolerance are not well characterised. Here we investigated the molecular mechanisms of enhanced tolerance of a lab selected strain of Orius laevigatus (Fieber) to the commonly used biopesticide spinosad. Transcriptomic analysis showed that spinosad tolerance is not a result of overexpressed detoxification genes. Molecular analysis of the target site for spinosyns, the nicotinic acetylcholine receptor (nAChR), revealed increased expression of truncated transcripts of the nAChR α6 subunit in the spinosad selected strain, a mechanism of resistance which was described previously in insect pest species. Collectively, our results demonstrate the mechanisms by which some beneficial biological control agents can evolve insecticide tolerance and will inform the development and deployment of insecticide-tolerant natural enemies in integrated pest management strategies.

The role of the Bemisia tabaci and Trialeurodes vaporariorum cytochrome-P450 clade CYP6DPx in resistance to nicotine and neonicotinoids.

The alkaloid, nicotine, produced by tobacco and other Solanaceae as an anti-herbivore defence chemical is one of the most toxic natural insecticides in nature. However, some insects, such as the whitefly species, Trialeurodes vaporariorum and Bemisia tabaci show strong tolerance to this allelochemical and can utilise tobacco as a host. Here, we used biological, molecular and functional approaches to investigate the role of cytochrome P450 enzymes in nicotine tolerance in T. vaporariorum and B. tabaci. Insecticide bioassays revealed that feeding on tobacco resulted in strong induced tolerance to nicotine in both species. Transcriptome profiling of both species reared on tobacco and bean hosts revealed profound differences in the transcriptional response these host plants. Interrogation of the expression of P450 genes in the host-adapted lines revealed that P450 genes belonging to the CYP6DP subfamily are strongly upregulated in lines reared on tobacco. Functional characterisation of these P450s revealed that CYP6DP1 and CYP6DP2 of T. vaporariorum and CYP6DP3 of B. tabaci confer resistance to nicotine in vivo. These three genes, in addition to the B. tabaci P450 CYP6DP5, were also found to confer resistance to the neonicotinoid imidacloprid. Our data provide new insight into the molecular basis of nicotine resistance in insects and illustrates how divergence in the evolution of P450 genes in this subfamily in whiteflies may have impacted the extent to which different species can tolerate a potent natural insecticide.

Bulk segregant mapping and transcriptome analyses reveal the molecular mechanisms of spinetoram resistance in Spodoptera frugiperda.

The evolution of resistance to insecticides poses a significant threat to pest management programs. Understanding the molecular mechanisms underlying insecticide resistance is essential to design sustainable pest control and resistance management programs. The fall armyworm, Spodoptera frugiperda, is an important insect pest of many crops and has a remarkable ability to evolve resistance to insecticides. In this study, we employed bulk segregant analysis (BSA) combined with DNA and RNA sequencing to characterize the molecular basis of spinetoram resistance in S. frugiperda. Analysis of genomic data derived from spinetoram selected and unselected bulks and the spinetoram-resistant and susceptible parental strains led to the identification of a three-nucleotide deletion in the gene encoding the nicotinic acetylcholine receptor α6 subunit (nAChR α6). Transcriptome profiling identified the upregulation of few genes encoding detoxification enzymes associated with spinetoram resistance. Thus, spinetoram resistance in S. frugiperda appears to be mediated mainly by target site insensitivity with a minor role of detoxification enzymes. Our findings provide insight into the mechanisms underpinning resistance to spinetoram in S. frugiperda and will inform the development of strategies to control this highly damaging, globally distributed crop pest.

Whole-chromosome hitchhiking driven by a male-killing endosymbiont.

Neo-sex chromosomes are found in many taxa, but the forces driving their emergence and spread are poorly understood. The female-specific neo-W chromosome of the African monarch (or queen) butterfly Danaus chrysippus presents an intriguing case study because it is restricted to a single 'contact zone' population, involves a putative colour patterning supergene, and co-occurs with infection by the male-killing endosymbiont Spiroplasma. We investigated the origin and evolution of this system using whole genome sequencing. We first identify the 'BC supergene', a broad region of suppressed recombination across nearly half a chromosome, which links two colour patterning loci. Association analysis suggests that the genes yellow and arrow in this region control the forewing colour pattern differences between D. chrysippus subspecies. We then show that the same chromosome has recently formed a neo-W that has spread through the contact zone within approximately 2,200 years. We also assembled the genome of the male-killing Spiroplasma, and find that it shows perfect genealogical congruence with the neo-W, suggesting that the neo-W has hitchhiked to high frequency as the male-killer has spread through the population. The complete absence of female crossing-over in the Lepidoptera causes whole-chromosome hitchhiking of a single neo-W haplotype, carrying a single allele of the BC supergene and dragging multiple non-synonymous mutations to high frequency. This has created a population of infected females that all carry the same recessive colour patterning allele, making the phenotypes of each successive generation highly dependent on uninfected male immigrants. Our findings show how hitchhiking can occur between the physically unlinked genomes of host and endosymbiont, with dramatic consequences.

MAPK-directed activation of the whitefly transcription factor CREB leads to P450-mediated imidacloprid resistance.

The evolution of insect resistance to pesticides poses a continuing threat to agriculture and human health. While much is known about the proximate molecular and biochemical mechanisms that confer resistance, far less is known about the regulation of the specific genes/gene families involved, particularly by trans-acting factors such as signal-regulated transcription factors. Here we resolve in fine detail the trans-regulation of CYP6CM1, a cytochrome P450 that confers resistance to neonicotinoid insecticides in the whitefly Bemisia tabaci, by the mitogen-activated protein kinase (MAPK)-directed activation of the transcription factor cAMP-response element binding protein (CREB). Reporter gene assays were used to identify the putative promoter of CYP6CM1, but no consistent polymorphisms were observed in the promoter of a resistant strain of B. tabaci (imidacloprid-resistant, IMR), which overexpresses this gene, compared to a susceptible strain (imidacloprid-susceptible, IMS). Investigation of potential trans-acting factors using in vitro and in vivo assays demonstrated that the bZIP transcription factor CREB directly regulates CYP6CM1 expression by binding to a cAMP-response element (CRE)-like site in the promoter of this gene. CREB is overexpressed in the IMR strain, and inhibitor, luciferase, and RNA interference assays revealed that a signaling pathway of MAPKs mediates the activation of CREB, and thus the increased expression of CYP6CM1, by phosphorylation-mediated signal transduction. Collectively, these results provide mechanistic insights into the regulation of xenobiotic responses in insects and implicate both the MAPK-signaling pathway and a transcription factor in the development of pesticide resistance.

The Role of Cytochrome P450s in Insect Toxicology and Resistance.

Insect cytochrome P450 monooxygenases (P450s) perform a variety of important physiological functions, but it is their role in the detoxification of xenobiotics, such as natural and synthetic insecticides, that is the topic of this review. Recent advances in insect genomics and postgenomic functional approaches have provided an unprecedented opportunity to understand the evolution of insect P450s and their role in insect toxicology. These approaches have also been harnessed to provide new insights into the genomic alterations that lead to insecticide resistance, the mechanisms by which P450s are regulated, and the functional determinants of P450-mediated insecticide resistance. In parallel, an emerging body of work on the role of P450s in defining the sensitivity of beneficial insects to insecticides has been developed. The knowledge gained from these studies has applications for the management of P450-mediated resistance in insect pests and can be leveraged to safeguard the health of important beneficial insects.

Genome-wide transcriptomic changes reveal the genetic pathways involved in insect migration.

Insects are capable of extraordinary feats of long-distance movement that have profound impacts on the function of terrestrial ecosystems. The ability to undertake these movements arose multiple times through the evolution of a suite of traits that make up the migratory syndrome, however the underlying genetic pathways involved remain poorly understood. Migratory hoverflies (Diptera: Syrphidae) are an emerging model group for studies of migration. They undertake seasonal movements in huge numbers across large parts of the globe and are important pollinators, biological control agents and decomposers. Here, we assembled a high-quality draft genome of the marmalade hoverfly (Episyrphus balteatus). We leveraged this genomic resource to undertake a genome-wide transcriptomic comparison of actively migrating Episyrphus, captured from a high mountain pass as they flew south to overwinter, with the transcriptomes of summer forms which were non-migratory. We identified 1543 genes with very strong evidence for differential expression. Interrogation of this gene set reveals a remarkable range of roles in metabolism, muscle structure and function, hormonal regulation, immunity, stress resistance, flight and feeding behaviour, longevity, reproductive diapause and sensory perception. These features of the migrant phenotype have arisen by the integration and modification of pathways such as insulin signalling for diapause and longevity, JAK/SAT for immunity, and those leading to octopamine production and fuelling to boost flight capabilities. Our results provide a powerful genomic resource for future research, and paint a comprehensive picture of global expression changes in an actively migrating insect, identifying key genomic components involved in this important life-history strategy.

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.

Adaptation of Helicoverpa armigera to Soybean Peptidase Inhibitors Is Associated with the Transgenerational Upregulation of Serine Peptidases.

Molecular phenotypes induced by environmental stimuli can be transmitted to offspring through epigenetic inheritance. Using transcriptome profiling, we show that the adaptation of Helicoverpa armigera larvae to soybean peptidase inhibitors (SPIs) is associated with large-scale gene expression changes including the upregulation of genes encoding serine peptidases in the digestive system. Furthermore, approximately 60% of the gene expression changes induced by SPIs persisted in the next generation of larvae fed on SPI-free diets including genes encoding regulatory, oxidoreductase, and protease functions. To investigate the role of epigenetic mechanisms in regulating SPI adaptation, the methylome of the digestive system of first-generation larvae (fed on a diet with and without SPIs) and of the progeny of larvae exposed to SPIs were characterized. A comparative analysis between RNA-seq and Methyl-seq data did not show a direct relationship between differentially methylated and differentially expressed genes, while trypsin and chymotrypsin genes were unmethylated in all treatments. Rather, DNA methylation potential epialleles were associated with transcriptional and translational controls; these may play a regulatory role in the adaptation of H. armigera to SPIs. Altogether, our findings provided insight into the mechanisms of insect adaptation to plant antiherbivore defense proteins and illustrated how large-scale transcriptional reprograming of insect genes can be transmitted across generations.

A toxicogenomics approach reveals characteristics supporting the honey bee (Apis mellifera L.) safety profile of the butenolide insecticide flupyradifurone.

Flupyradifurone, a novel butenolide insecticide, selectively targets insect nicotinic acetylcholine receptors (nAChRs), comparable to structurally different insecticidal chemotypes such as neonicotinoids and sulfoximines. However, flupyradifurone was shown in acute toxicity tests to be several orders of magnitude less toxic to western honey bee (Apis mellifera L.) than many other insecticides targeting insect nAChRs. The underlying reasons for this difference in toxicity remains unknown and were investigated here. Pharmacokinetic studies after contact application of [14C]flupyradifurone to honey bees revealed slow uptake, with internalized compound degraded into a few metabolites that are all practically non-toxic to honey bees in both oral and contact bioassays. Furthermore, receptor binding studies revealed a lack of high-affinity binding of these metabolites to honey bee nAChRs. Screening of a library of 27 heterologously expressed honey bee cytochrome P450 enzymes (P450s) identified three P450s involved in the detoxification of flupyradifurone: CYP6AQ1, CYP9Q2 and CYP9Q3. Transgenic Drosophila lines ectopically expressing CYP9Q2 and CYP9Q3 were significantly less susceptible to flupyradifurone when compared to control flies, confirming the importance of these P450s for flupyradifurone metabolism in honey bees. Biochemical assays using the fluorescent probe substrate 7-benzyloxymethoxy-4-(trifluoromethyl)-coumarin (BOMFC) indicated a weak, non-competitive inhibition of BOMFC metabolism by flupyradifurone. In contrast, the azole fungicides prochloraz and propiconazole were strong nanomolar inhibitors of these flupyradifurone metabolizing P450s, explaining their highly synergistic effects in combination with flupyradifurone as demonstrated in acute laboratory contact toxicity tests of adult bees. Interestingly, the azole fungicide prothioconazole is only slightly synergistic in combination with flupyradifurone - an observation supported by molecular P450 inhibition assays. Such molecular assays have value in the prediction of potential risks posed to bees by flupyradifurone mixture partners under applied conditions. Quantitative PCR confirmed the expression of the identified P450 genes in all honey bee life-stages, with highest expression levels observed in late larvae and adults, suggesting honey bees have the capacity to metabolize flupyradifurone across all life-stages. These findings provide a biochemical explanation for the low intrinsic toxicity of flupyradifurone to honey bees and offer a new, more holistic approach to support bee pollinator risk assessment by molecular means.