A cytochrome p450 conserved in insects is involved in cuticle formation

Changes in heme oxygenase level during development affect the adult life of Drosophila melanogaster

Heme oxygenase (HO) has been shown to control various cellular processes in both mammals and Drosophila melanogaster. Here, we investigated how changes in HO levels in neurons and glial cells during development affect adult flies, by using the TARGET Drosophila system to manipulate the expression of the ho gene. The obtained data showed differences in adult survival, maximum lifespan, climbing, locomotor activity, and sleep, which depended on the level of HO (after ho up-regulation or downregulation), the timing of expression (chronic or at specific developmental stages), cell types (neurons or glia), sex (males or females), and age of flies. In addition to ho, the effects of changing the mRNA level of the Drosophila CNC factor gene (NRF2 homolog in mammals and master regulator of HO), were also examined to compare with those observed after changing ho expression. We showed that HO levels in neurons and glia must be maintained at an appropriate physiological level during development to ensure the well-being of adults. We also found that the downregulation of ho in either neurons or glia in the brain is compensated by ho expressed in the retina.

Using machine learning to predict protein-protein interactions between a zombie ant fungus and its carpenter ant host

Parasitic fungi produce proteins that modulate virulence, alter host physiology, and trigger host responses. These proteins, classified as a type of "effector," often act via protein-protein interactions (PPIs). The fungal parasite Ophiocordyceps camponoti-floridani (zombie ant fungus) manipulates Camponotus floridanus (carpenter ant) behavior to promote transmission. The most striking aspect of this behavioral change is a summit disease phenotype where infected hosts ascend and attach to an elevated position. Plausibly, interspecific PPIs drive aspects of Ophiocordyceps infection and host manipulation. Machine learning PPI predictions offer high-throughput methods to produce mechanistic hypotheses on how this behavioral manipulation occurs. Using D-SCRIPT to predict host-parasite PPIs, we found ca. 6000 interactions involving 2083 host proteins and 129 parasite proteins, which are encoded by genes upregulated during manipulated behavior. We identified multiple overrepresentations of functional annotations among these proteins. The strongest signals in the host highlighted neuromodulatory G-protein coupled receptors and oxidation-reduction processes. We also detected Camponotus structural and gene-regulatory proteins. In the parasite, we found enrichment of Ophiocordyceps proteases and frequent involvement of novel small secreted proteins with unknown functions. From these results, we provide new hypotheses on potential parasite effectors and host targets underlying zombie ant behavioral manipulation.

Genome-Wide Identification and Expression Pattern of Cytochrome P450 Genes in the Social Aphid Pseudoregma bambucicola

Cytochrome P450 monooxygenases (P450s) have a variety of functions, including involvement in the metabolism of exogenous substances and the synthesis and degradation of endogenous substances, which are important for the growth and development of insects. Pseudoregma bambucicola is a social aphid that produces genetically identical but morphologically and behaviorally distinct first-instar soldiers and normal nymphs within colonies. In this study, we identified 43 P450 genes based on P. bambucicola genome data. Phylogenetic analysis showed that these genes were classified into 4 clans, 13 families, and 23 subfamilies. The CYP3 and CYP4 clans had a somewhat decreased number of genes. In addition, differential gene expression analysis based on transcriptome data showed that several P450 genes, including CYP18A1, CYP4G332, and CYP4G333, showed higher expression levels in soldiers compared to normal nymphs and adult aphids. These genes may be candidates for causing epidermal hardening and developmental arrest in soldiers. This study provides valuable data and lays the foundation for the study of functions of P450 genes in the social aphid P. bambucicola.

Gut transcriptome analysis of P450 genes and cytochrome P450 reductase in three moth species feeding on gymnosperms or angiosperms

Cytochrome P450 enzymes (P450s, CYPs) are a superfamily of heme–thiolate proteins involved in the metabolism of endogenous and exogenous substances in insects. In this study, the identification of putative P450 proteins was done and the elimination of the repeated sequences resulted in 57 proteins from Gastropacha populifolia, 63 proteins from Dendrolimus punctatus, and 53 proteins from Dendrolimus tabulaeformis. The putative P450 proteins were aligned together with seven other insect species based on five conserved domains. A total of ten co-orthologous groups were identified. Interestingly, one co-orthologous gene, CYP4g15 in CYP4 clan, was identified and its 3D structure analysis showed that the highly conserved sites of the predicted motifs were close to the active sites of P450. Furthermore, this study revealed that insect CYP4g15 and two bacteria cytochrome P450 were monophyletic. This suggests that insects CYP4g15 are not only functionally conserved but also an ancient gene originating from different bacteria species.

Genome-wide and expression-profiling analyses of the cytochrome P450 genes in Tenebrionidea

Cytochrome P450 monooxygenases (CYPs) are present in almost all areas of the tree of life. As one of the largest and most diverse superfamilies of multifunctional enzymes, they play important roles in the metabolism of xenobiotics and biosynthesis of endogenous compounds, shaping the success of insects. In this study, the CYPome (an omics term for all the CYP genes in a genome) diversification was examined in the four Tenebrionidea species through genome-wide analysis. A total of 483 CYP genes were identified, of which 103, 157, 122, and 101 were respectively deciphered from the genomes of Tebebrio molitor, Asbolus verucosus, Hycleus cichorii and Hycleus phaleratus. These CYPs were classified into four major clans (mitochondrial, CYP2, CYP3, and CYP4), and clans CYP3 and CYP4 are most diverse. Phylogenetic analysis showed that most CYPs of these Tenebrionidea beetles from each clan had a very close 1:1 orthology to each other, suggesting that they originate closely and have evolutionally conserved function. Expression analysis at different developmental stages and in various tissues showed the life stage-, gut-, salivary gland-, fat body-, Malpighian tubule-, antennae-, ovary- and testis-specific expression patterns of T. molitor CYP genes, implying their various potential roles in development, detoxification, immune response, digestion, olfaction, and reproduction. Our studies provide a platform to understand the evolution of Tenebrionidea CYP gene superfamily, and a basis for further functional investigation of the T. molitor CYPs involved in various biological processes.

Comparative Transcriptome Analysis to Reveal Differentially Expressed Cytochrome P450 in Response to Imidacloprid in the Aphid Lion, Chrysoperla zastrowi sillemi (Esben-Petersen)

The aphid lion, Chrysoperla zastrowi sillemi (Neuroptera: Chrysopidae) is a highly effective beneficial predator of many agricultural pests and has developed resistance to several insecticides. Understanding the molecular mechanism of insecticide resistance in the predators is crucial for its effective application in IPM programs. Therefore, transcriptomes of imidacloprid-resistant and susceptible strains have been assessed using RNA-seq. Cytochrome P450 is one of the important gene families involved in xenobiotic metabolism. Hence, our study focused on the CYP gene family where mining, nomenclature, and phylogenetic analysis revealed a total of 95 unique CYP genes with considerable expansion in CYP3 and CYP4 clans. Further, differential gene expression (DGE) analysis revealed ten CYP genes from CYP3 and CYP4 clans to be differentially expressed, out of which nine genes (CYP4419A1, CYP4XK1, CYP4416A10, CYP4416A-fragment8, CYP6YL1, CYP6YH6, CYP9GK-fragment16, CYP9GN2, CYP9GK6) were downregulated and one (CYP9GK3) was upregulated in the resistant strain as compared to the susceptible strain. Expression validation by quantitative real-time PCR (qRT-PCR) is consistent with the DGE results. The expansion and differential expression of CYP genes may be an indicator of the capacity of the predator to detoxify a particular group of insecticides.

Identification and characterization of candidate detoxification genes in Pharsalia antennata Gahan (Coleoptera: Cerambycidae)

The wood-boring beetles, including the majority of Cerambycidae, have developed the ability to metabolize a variety of toxic compounds derived from host plants and the surrounding environment. However, detoxification mechanisms underlying the evolutionary adaptation of a cerambycid beetle Pharsalia antennata to hosts and habitats are largely unexplored. Here, we characterized three key gene families in relation to detoxification (cytochrome P450 monooxygenases: P450s, carboxylesterases: COEs and glutathione-S-transferases: GSTs), by combinations of transcriptomics, gene identification, phylogenetics and expression profiles. Illumina sequencing generated 668,701,566 filtered reads in 12 tissues of P. antennata, summing to 100.28 gigabases data. From the transcriptome, 215 genes encoding 106 P450s, 77 COEs and 32 GSTs were identified, of which 107 relatives were differentially expressed genes. Of the identified 215 genes, a number of relatives showed the orthology to those in Anoplophora glabripennis, revealing 1:1 relationships in 94 phylogenetic clades. In the trees, P. antennata detoxification genes mainly clustered into one or two subfamilies, including 64 P450s in the CYP3 clan, 33 COEs in clade A, and 20 GSTs in Delta and Epsilon subclasses. Combining transcriptomic data and PCR approaches, the numbers of detoxification genes expressed in abdomens, antennae and legs were 188, 148 and 141, respectively. Notably, some genes exhibited significantly sex-biased levels in antennae or legs of both sexes. The findings provide valuable reference resources for further exploring xenobiotics metabolism and odorant detection in P. antennata.

The insecticidal capacity of ethanol extract from Cascabela peruviana (L.) Lippold against fruit fly

Cascabela peruviana (L.) Lippold (C. peruviana) has been extensively used for its antifungal and antibacterial properties. However, its role in anti-insect is still under investigation. To investigate the ability of the ethanol extract of C. peruviana against insects, we used the fruit fly (Drosophila melanogaster) as a model to gain more insight into the toxic effects of this extract. We found that the ethanol extract from the stem and leaves of C. peruviana was effective against insects and contained polyphenol and flavonoid compounds. C. peruviana could induce mortality of 2^nd-instar larvae and reduce growth and reproduction of fruit flies. Interestingly, the toxicity of C. peruviana extract has been remained to affect the development of the next generation of fruit flies. The locomotor activity and feeding ability of the F1 generation of this insect were significantly reduced by C. peruviana. In addition, flavonoids and polyphenols, as well as saponins and tannins were detected in the ethanol extract of C. peruviana. We assume that the ability of the extract of C. peruviana to control insects may be related to the presence of high levels of these compounds. The findings highlighted that the extract from the leaves of Cascabela peruviana has the potential to be used as an insecticide.

A near-chromosome level genome assembly of the European hoverfly, Sphaerophoria rueppellii (Diptera: Syrphidae), provides comparative insights into insecticide resistance-related gene family evolution

BACKGROUND

Sphaerophoria rueppellii, a European species of hoverfly, is a highly effective beneficial predator of hemipteran crop pests including aphids, thrips and coleopteran/lepidopteran larvae in integrated pest management (IPM) programmes. It is also a key pollinator of a wide variety of important agricultural crops. No genomic information is currently available for S. rueppellii. Without genomic information for such beneficial predator species, we are unable to perform comparative analyses of insecticide target-sites and genes encoding metabolic enzymes potentially responsible for insecticide resistance, between crop pests and their predators. These metabolic mechanisms include several gene families - cytochrome P450 monooxygenases (P450s), ATP binding cassette transporters (ABCs), glutathione-S-transferases (GSTs), UDP-glycosyltransferases (UGTs) and carboxyl/choline esterases (CCEs).

METHODS AND FINDINGS

In this study, a high-quality near-chromosome level de novo genome assembly (as well as a mitochondrial genome assembly) for S. rueppellii has been generated using a hybrid approach with PacBio long-read and Illumina short-read data, followed by super scaffolding using Hi-C data. The final assembly achieved a scaffold N50 of 87Mb, a total genome size of 537.6Mb and a level of completeness of 96% using a set of 1,658 core insect genes present as full-length genes. The assembly was annotated with 14,249 protein-coding genes. Comparative analysis revealed gene expansions of CYP6Zx P450s, epsilon-class GSTs, dietary CCEs and multiple UGT families (UGT37/302/308/430/431). Conversely, ABCs, delta-class GSTs and non-CYP6Zx P450s showed limited expansion. Differences were seen in the distributions of resistance-associated gene families across subfamilies between S. rueppellii and some hemipteran crop pests. Additionally, S. rueppellii had larger numbers of detoxification genes than other pollinator species.

CONCLUSION AND SIGNIFICANCE

This assembly is the first published genome for a predatory member of the Syrphidae family and will serve as a useful resource for further research into selectivity and potential tolerance of insecticides by beneficial predators. Furthermore, the expansion of some gene families often linked to insecticide resistance and selectivity may be an indicator of the capacity of this predator to detoxify IPM selective insecticides. These findings could be exploited by targeted insecticide screens and functional studies to increase effectiveness of IPM strategies, which aim to increase crop yields by sustainably and effectively controlling pests without impacting beneficial predator populations.

The flying insect thoracic cuticle is heterogenous in structure and in thickness-dependent modulus gradation

The thorax is a specialized structure central to insect flight. In the thorax, flight muscles are surrounded by a thin layer of cuticle. The structure, composition, and material properties of this chitinous structure may influence the efficiency of the thorax in flight. However, these properties, as well as their variation throughout the thorax and between insect taxa, are not known. We provide a multi-faceted assessment of thorax cuticle for fliers with asynchronous (honey bee; Apis mellifera) and synchronous (hawkmoth; Manduca sexta) muscles. These muscle types are defined by the relationship between their activation frequency and the insect's wingbeat frequency. We investigated cuticle structure using histology, resilin distribution through confocal laser scanning microscopy, and modulus gradation with nanoindentation. Our results suggest that thorax cuticle properties are highly dependent on anatomical region and species. Modulus gradation, but not mean modulus, differed between the two types of fliers. In some regions, A. mellifera had a positive linear modulus gradient from cuticle interior to exterior of about 2 GPa. In M. sexta, modulus values through cuticle thickness were not well represented by linear fits. We utilized finite element modeling to assess how measured modulus gradients influenced maximum stress in cuticle. Stress was reduced when cuticle with a linear gradient was compressed from the high modulus side. These results support the protective role of the A. mellifera thorax cuticle. Our multi-faceted assessment advances our understanding of thorax cuticle structural and material heterogeneity and the potential benefits of material gradation to flying insects. STATEMENT OF SIGNIFICANCE: The insect thorax is essential for efficient flight but questions remain about the contribution of the exoskeletal cuticle. We investigated the microscale properties of the thorax cuticle, a crucial step to determine its role in flight. Techniques including histology, nanoindentation, and confocal laser scanning microscopy revealed that cuticle properties vary through cuticle thickness, by thorax region, and between species with asynchronous (honey bee; Apis mellifera) and synchronous (hawkmoth; Manduca sexta) muscles. This variation highlights the importance of high resolution cuticle assessment for flying insect lineages and points to factors that may (modulus gradation) and may not (mean modulus) contribute to different flight forms. Understanding material variation in the thorax may inform design of technologies inspired by insects, such as mobile micro robots.

Gut transcriptome of two bark beetle species stimulated with the same kairomones reveals molecular differences in detoxification pathways

Dendroctonus bark beetles are the most destructive agents in coniferous forests. These beetles come into contact with the toxic compounds of their host's chemical defenses throughout their life cycle, some of which are also used by the insects as kairomones to select their host trees during the colonization process. However, little is known about the molecular mechanisms by which the insects counteract the toxicity of these compounds. Here, two sibling species of bark beetles, D. valens and D. rhizophagus, were stimulated with vapors of a blend of their main kairomones (α-pinene, β-pinene and 3-carene), in order to compare the transcriptional response of their gut. A total of 48 180 unigenes were identified in D. valens and 43 704 in D. rhizophagus, in response to kairomones blend. The analysis of differential gene expression showed a transcriptional response in D. valens (739 unigenes, 0.58–10.36 Log2FC) related to digestive process and in D. rhizophagus (322 unigenes 0.87–13.08 Log2FC) related to xenobiotics metabolism. The expression profiles of detoxification genes mainly evidenced the up-regulation of COEs and GSTs in D. valens, and the up-regulation of P450s in D. rhizophagus. Results suggest that terpenes metabolism comes accompanied by an integral hormetic response, result of compensatory mechanisms, including the activation of other metabolic pathways, to ensure the supply of energy and the survival of organisms which is specific for each species, according to its life history and ecological strategy.

Draft genome of the blister beetle, Epicauta chinensis

Existence of cantharidin (CTD) in blister beetles is a significant ecological adaptive mechanism that defends against predators and regulates courtship and mating behaviors. To better understand CTD biosynthetic information as well as its biology and pharmacology, we assembled a genome of 151.88 Mb for Epicauta chinensis using PacBio sequencing technology. Gene annotation yielded 249,238 repeats, 527 non-coding RNAs and 12,520 protein-coding genes. Compared to other 11 insects, expansions of gene families in E. chinensis for most core gene families likely associated with environmental adaptation, such as chemoreception, immunity, and detoxification. We further annotated P450s and immune-related genes, a total of 117 putative P450s comprising 7 CYP2, 67 CYP3, 36 CYP4, and 7 mitochondrial P450s and 281 immune-related genes were identified. Comparative analysis of the insect immune repertoires indicated presence of immune genes detected only from Coleopteran insects such as MD2-like. This suggested a lineage-specific gene evolution for Coleopteran insects. Based on the gene family evolution analysis, we identified two probable candidate genes including CYP4TT1 and phytanoyl-CoA dioxygenase for CTD biosynthesis. The high-quality reference genome of E. chinensis provides the genetic basis for further investigation of CTD biosynthesis and in-depth studies of the development and evolution of blister beetles.

First evidence for thermal tolerance benefits of the bacterial symbiont Cardinium in an invasive whitefly, Bemisia tabaci

BACKGROUD

Cardinium symbiont is a maternally inherited bacterial endosymbiont and widely spreads in arthropods including Bemisia tabaci (Hemiptera: Aleyrodidae). However, the potential role of Cardinium played in the biology of their hosts is largely unknown. In two genetic lines (i.e. LS and SG lines) of B. tabaci MED, collected from different locations in China, we tested the effects of Cardinium on the performance of the host whitefly under a constant high temperature (31 °C) using the age-stage two-sex life table method, and explored the genes influenced by Cardinium-infection by RNA-sequencing.

RESULTS

We found that Cardinium did provide protection of B. tabaci against heat stress under 31 °C. However, there was a significant connection between Cardinium-infection and whitefly genetic backgrounds. Performance revealed that Cardinium infection can increase the longevity of both female and male adults and oviposition periods in both lines, but it also conferred benefits of fecundity and pre-adult period to LS line. Additionally, the population parameters such as intrinsic rate of increase (r), finite rate of increase (λ) and mean generation time (T) demonstrated that Cardinium infection conferred fitness benefits to LS line but not to SG line. Transcriptome analysis indicated that several genes related to homeostasis and metamorphosis such as ubiquitin-related genes were highly expressed in Cardinium-infected B. tabaci.

CONCLUSION

The research provided the first evidence that Cardinium can increase the thermal tolerance of whitefly, which may be associated with host genetic background.

Up-regulation of apoptotic- and cell survival-related gene pathways following exposures of western corn rootworm to B. thuringiensis crystalline pesticidal proteins in transgenic maize roots

BACKGROUND

Resistance of pest insect species to insecticides, including B. thuringiensis (Bt) pesticidal proteins expressed by transgenic plants, is a threat to global food security. Despite the western corn rootworm, Diabrotica virgifera virgifera, being a major pest of maize and having populations showing increasing levels of resistance to hybrids expressing Bt pesticidal proteins, the cell mechanisms leading to mortality are not fully understood.

RESULTS

Twenty unique RNA-seq libraries from the Bt susceptible D. v. virgifera inbred line Ped12, representing all growth stages and a range of different adult and larval exposures, were assembled into a reference transcriptome. Ten-day exposures of Ped12 larvae to transgenic Bt Cry3Bb1 and Gpp34/Tpp35Ab1 maize roots showed significant differential expression of 1055 and 1374 transcripts, respectively, compared to cohorts on non-Bt maize. Among these, 696 were differentially expressed in both Cry3Bb1 and Gpp34/Tpp35Ab1 maize exposures. Differentially-expressed transcripts encoded protein domains putatively involved in detoxification, metabolism, binding, and transport, were, in part, shared among transcripts that changed significantly following exposures to the entomopathogens Heterorhabditis bacteriophora and Metarhizium anisopliae. Differentially expressed transcripts in common between Bt and entomopathogen treatments encode proteins in general stress response pathways, including putative Bt binding receptors from the ATP binding cassette transporter superfamily. Putative caspases, pro- and anti-apoptotic factors, as well as endoplasmic reticulum (ER) stress-response factors were identified among transcripts uniquely up-regulated following exposure to either Bt protein.

CONCLUSIONS

Our study suggests that the up-regulation of genes involved in ER stress management and apoptotic progression may be important in determining cell fate following exposure of susceptible D. v. virgifera larvae to Bt maize roots. This study provides novel insights into insect response to Bt intoxication, and a possible framework for future investigations of resistance mechanisms.

Crowdsourced Identification of Potential Target Genes for CTV Induced Gene Silencing for Controlling the Citrus Greening Vector Diaphorina citri

Citrus Greening or Huanglongbing (HLB) is a disease of citrus, causing high reduction in citrus production and is transmitted by the Asian citrus psyllid Diaphorina citri Kuwayama vectoring a phloem-limited bacterium Candidatus Liberibacter sp. We report research results using crowdsourcing challenge strategy identifying potential gene targets in D. citri to control the insect using RNA interference (RNAi). From 63 submitted sequences, 43 were selected and tested by feeding them to D. citri using artificial diet assays. After feeding on artificial diet, the three most effective dsRNAs causing 30% mortality above control silenced genes expressing iron-sulfur cluster subunit of the mitochondrial electron transport chain complex (Rieske), heme iron-binding terminal oxidase enzyme (Cytochrome P450) and tetrahydrobiopterin (BH₄) pathway enzyme (Pterin 4α-Carbinolamine Dehydratase). These sequences were cloned into a citrus phloem-limited virus (Citrus tristeza virus, CTV T36) expressing dsRNA against these target genes in citrus. The use of a viral mediated “para-transgenic” citrus plant system caused higher mortality to adult D. citri than what was observed using artificial diet, reaching 100% when detached citrus leaves with the engineered CTV expressing dsRNA were fed to adult D. citri. Using this approach, a virus-induced gene silencing (VIGS) can be used to test future transgenic cultivars before genetically engineering citrus. RNA Seq analysis after feeding D. citri CTV-RIE on infected leaves identified transcriptionally modified genes located upstream and downstream of the targeted RIE gene. These genes were annotated showing that many are associated with the primary function of the Rieske gene that was targeted by VIGS.

Matrix metalloproteinases are involved in eclosion and wing expansion in the American cockroach, Periplaneta americana

Matrix metalloproteinases (MMPs) are the major proteinases that process or degrade numerous extracellular matrix (ECM) components and are evolutionarily conserved from nematodes to humans. During molting in insects, the old cuticle is removed and replaced by a new counterpart. Although the regulatory mechanisms of hormones and nutrients in molting have been well studied, very little is known about the roles of ECM-modifying enzymes in this process. Here, we found that MMPs are necessary for imaginal molting of the American cockroach, Periplaneta americana. Inhibition of Mmp activity via inhibitor treatment led to the failure of eclosion and wing expansion. Five Mmps genes were identified from the P. americana genome, and PaMmp2 played the dominant roles during molting. Further microscopic investigations showed that newly formed adult cuticles were attenuated and that then chitin content was reduced upon Mmp inhibition. Transcriptomic analysis of the integument demonstrated that multiple signaling and metabolic pathways were changed. Microscopic investigation of the wings showed that epithelial cells were restrained together because they were incapable of degrading the ECM upon Mmp inhibition. Transcriptomic analysis of the wing identified dozens of possible genes functioned in wing expansion. This is the first study to show the essential roles of Mmps in the nymph-adult transition of hemimetabolous insects.

Genome-enabled insights into the biology of thrips as crop pests

BACKGROUND

The western flower thrips, Frankliniella occidentalis (Pergande), is a globally invasive pest and plant virus vector on a wide array of food, fiber, and ornamental crops. The underlying genetic mechanisms of the processes governing thrips pest and vector biology, feeding behaviors, ecology, and insecticide resistance are largely unknown. To address this gap, we present the F. occidentalis draft genome assembly and official gene set.

RESULTS

We report on the first genome sequence for any member of the insect order Thysanoptera. Benchmarking Universal Single-Copy Ortholog (BUSCO) assessments of the genome assembly (size = 415.8 Mb, scaffold N50 = 948.9 kb) revealed a relatively complete and well-annotated assembly in comparison to other insect genomes. The genome is unusually GC-rich (50%) compared to other insect genomes to date. The official gene set (OGS v1.0) contains 16,859 genes, of which ~ 10% were manually verified and corrected by our consortium. We focused on manual annotation, phylogenetic, and expression evidence analyses for gene sets centered on primary themes in the life histories and activities of plant-colonizing insects. Highlights include the following: (1) divergent clades and large expansions in genes associated with environmental sensing (chemosensory receptors) and detoxification (CYP4, CYP6, and CCE enzymes) of substances encountered in agricultural environments; (2) a comprehensive set of salivary gland genes supported by enriched expression; (3) apparent absence of members of the IMD innate immune defense pathway; and (4) developmental- and sex-specific expression analyses of genes associated with progression from larvae to adulthood through neometaboly, a distinct form of maturation differing from either incomplete or complete metamorphosis in the Insecta.

CONCLUSIONS

Analysis of the F. occidentalis genome offers insights into the polyphagous behavior of this insect pest that finds, colonizes, and survives on a widely diverse array of plants. The genomic resources presented here enable a more complete analysis of insect evolution and biology, providing a missing taxon for contemporary insect genomics-based analyses. Our study also offers a genomic benchmark for molecular and evolutionary investigations of other Thysanoptera species.

Genomic and transcriptomic analyses in Drosophila suggest that the ecdysteroid kinase-like (EcKL) gene family encodes the 'detoxification-by-phosphorylation' enzymes of insects

Phosphorylation is a phase II detoxification reaction that, among animals, occurs near exclusively in insects, but the enzymes responsible have never been cloned or otherwise identified. We propose the hypothesis that members of the arthropod-specific ecdysteroid kinase-like (EcKL) gene family encode detoxicative kinases. To test this hypothesis, we annotated the EcKL gene family in 12 species of Drosophila and explored their evolution within the genus. Many ancestral EcKL clades are evolutionarily unstable and have experienced repeated gene gain and loss events, while others are conserved as single-copy orthologs. Leveraging multiple published gene expression datasets from D. melanogaster, and using the cytochrome P450s-a classical detoxification family-as a test case, we demonstrate relationships between xenobiotic induction, detoxification tissue-enriched expression and evolutionary instability in the EcKLs and the P450s. We devised a systematic method for identifying candidate detoxification genes in large gene families that is concordant with experimentally determined functions of P450 genes in D. melanogaster. Applying this method to the EcKLs suggested a significant proportion of these genes play roles in detoxification, and that the EcKLs may constitute a detoxification gene family in insects. Additionally, we estimate that between 11 and 16 uncharacterised D. melanogaster P450s are strong detoxification candidates. Lastly, we also found previously unreported genomic and transcriptomic variation in a number of EcKLs and P450s associated with toxic stress phenotypes using a targeted phenome-wide association study (PheWAS) approach in D. melanogaster, presenting multiple future avenues of research for detoxification genetics in this species.

Brown marmorated stink bug, Halyomorpha halys (Stål), genome: putative underpinnings of polyphagy, insecticide resistance potential and biology of a top worldwide pest

BACKGROUND

Halyomorpha halys (Stål), the brown marmorated stink bug, is a highly invasive insect species due in part to its exceptionally high levels of polyphagy. This species is also a nuisance due to overwintering in human-made structures. It has caused significant agricultural losses in recent years along the Atlantic seaboard of North America and in continental Europe. Genomic resources will assist with determining the molecular basis for this species' feeding and habitat traits, defining potential targets for pest management strategies.

RESULTS

Analysis of the 1.15-Gb draft genome assembly has identified a wide variety of genetic elements underpinning the biological characteristics of this formidable pest species, encompassing the roles of sensory functions, digestion, immunity, detoxification and development, all of which likely support H. halys' capacity for invasiveness. Many of the genes identified herein have potential for biomolecular pesticide applications.

CONCLUSIONS

Availability of the H. halys genome sequence will be useful for the development of environmentally friendly biomolecular pesticides to be applied in concert with more traditional, synthetic chemical-based controls.

Identification and characterization of detoxification genes in two cerambycid beetles, Rhaphuma horsfieldi and Xylotrechus quadripes (Coleoptera: Cerambycidae: Clytini)

The longhorned beetles, Rhaphuma horsfieldi and Xylotrechus quadripes, are two polyphagous insects with larvae feeding on different host plants. In this study, we identified and characterized three gene superfamilies of cytochrome P450s (CYPs), carboxylesterases (COEs) and glutathione-S-transferases (GSTs) involved in the detoxification of endobiotics (e.g., hormones and steroids) and xenobiotics (e.g., insecticides, sex pheromones and plant allelochemicals) through a combination approach of bioinformatics, phylogenetics, expression profiles and genomics. Transcriptome analyses led to the identification of 281 transcripts encoding 135 P450s, 108 COEs and 38 GSTs from the two beetles, coupled with comparative studies of detoxification genes among coleopteran species, suggesting a correlation between host range and the sizes of P450 or COE gene repertoires. The P450s of two beetles were phylogenetically classified into four clades, representing the majority of genes in the CYP3 clan. The COEs from R. horsfieldi and X. quadripes were separately grouped into 11 and 10 clades, and the GST superfamily was assigned into six clades. Expression profiles revealed that the detoxification genes were broadly expressed in various tissues as an implication of functional diversities. Ultimately and more importantly, five alternative splicing events in the Epsilon GSTs, including RhorGSTe7.1/GSTe7.2 and XquaGSTe3.1/GST3.2, were acquired in Coleoptera, in which these genes and their orthologs shared highly conserved gene structure. Our current study has complemented the resources for the detoxification genes in the family Cerambycidae, and allows for functional experiments to identify candidate molecular targets involved in pest resistance to insecticides like organophosphates, organochlorines and pyrethroids.

CYP303A1 has a conserved function in adult eclosion in Locusta migratoria and Drosophila melanogaster

Insect cytochrome P450 monooxygenases (CYPs) play essential roles in both xenobiotic metabolism and developmental processes. However, the exact physiological function of many CYP genes remains largely unknown. Screening the expression of the CYP genes from the CYP2 and mitochondrial CYP clans of Drosophila melanogaster revealed that Cyp303a1 is highly expressed in the pupal stage. Knockdown of CYP303A1 transcripts by RNAi using the Gal4/UAS system with a ubiquitous driver (tubulin-Gal4) in Drosophila or by dsRNA injection in the last nymph stage of Locusta migratoria resulted in severe defects in eclosion and lethality during and after adult emergence. In Drosophila, tissue-specific RNAi of Cyp303a1 with a wing-specific driver (MS1096-Gal4) revealed that Cyp303a1 was essential for wing extension. Stage-specific RNAi of Cyp303a1 using Gal80^ts for thermal-dependent-suppression found that the expression of Cyp303a1 at the middle pupal stage was absolutely required. Meanwhile, Cyp303a1 mutants exhibited more than 80% lethality at the late embryonic development stages. Embryonic lethality of the Cyp303a1 mutants was fully rescued by the ubiquitous overexpression of exogenous Cyp303a1. Taken together, we conclude that Cyp303a1 is indispensable for embryonic development and adult eclosion in D. melanogaster, the latter role being conserved over 400 million years of insect evolution.

Latrophilin participates in insecticide susceptibility through positively regulating CSP10 and partially compensated by OBPC01 in Tribolium castaneum

Latrophilin (LPH) is an adhesion G protein-coupled receptor (aGPCR) that participates in multiple essential physiological processes. Our previous studies have shown that lph is not only indispensable for the development and reproduction of red flour beetles (Tribolium castaneum), but also for their resistance against dichlorvos or carbofuran insecticides. However, the regulatory mechanism of lph-mediated insecticide susceptibility remains unclear. Here, we revealed that knockdown of lph in beetles resulted in opposing changes in two chemoreception genes, chemosensory protein 10 (CSP10) and odorant-binding protein C01 (OBPC01), in which the expression of TcCSP10 was downregulated, whereas the expression of TcOBPC01 was upregulated. TcCSP10 and TcOBPC01 were expressed at the highest levels in early pupal and late larval stages, respectively. High levels of expression of both these genes were observed in the heads (without antennae) of adults. TcCSP10 and TcOBPC01 were significantly induced by dichlorvos or carbofuran between 12 and 72 h (hrs) after exposure, suggesting that they are likely associated with increasing the binding affinity of insecticides, leading to a decrease in sensitivity to the insecticides. Moreover, once these two genes were knocked down, the susceptibility of the beetles to dichlorvos or carbofuran was enhanced. Additionally, RNA interference (RNAi) targeting of lph followed by exposure to dichlorvos or carbofuran also caused the opposing expression levels of TcCSP10 and TcOBPC01 compared to the expression levels of wild-type larvae treated with insecticides alone. All these results indicate that lph is involved in insecticide susceptibility through positively regulating TcCSP10; and the susceptibility could also further partially compensated for through the negative regulation of TcOBPC01 when lph was knockdown in the red flour beetle. Our studies shed new light on the molecular regulatory mechanisms of lph related to insecticide susceptibility.

T929I and K1774N mutation pair and M918L single mutation identified in the voltage-gated sodium channel gene of pyrethroid-resistant Thrips tabaci (Thysanoptera: Thripidae) in Japan

Pyrethroid-resistance in onion thrips, Thrips tabaci, has been reported in many countries including Japan. Identifying factors of the resistance is important to correctly monitoring the resistance in field populations. To identify pyrethroid-resistance related genes in T. tabaci in Japan, we performed RNA-Seq analysis of seven T. tabaci strains including two pyrethroid-resistant and five pyrethroid-susceptible strains. We identified a pair of single point mutations, T929I and K1774N, introducing two amino acid mutations, in the voltage-gated sodium channel gene, a pyrethroid target gene, in the two resistant strains. The K1774N is a newly identified mutation located in the fourth repeat domain of the sodium channel. Genotyping analysis of field-collected populations showed that most of the T. tabaci individuals in resistant populations carried the mutation pair, indicating that the mutation pair is closely associated with pyrethroid-resistance in Japan. Another resistance-related mutation, M918L, was also identified in part of the resistant populations. Most of the individuals with the mutation pair were arrhenotokous while all individuals with the M918L single mutation were thelytokous. The result of differentially expressed gene analysis revealed a small number of up-regulated detoxification genes in each resistant strain which might be involved in resistance to pyrethroid. However, no up-regulated detoxification genes common to the two resistant strains were detected. Our results indicate that the mutation pair in the sodium channel gene is the most important target for monitoring pyrethroid-resistance in T. tabaci, and that pyrethroid-resistant arrhenotokous individuals with the mutation pair are likely to be widely distributed in Japan.

Global transcriptome profiling and functional analysis reveal that tissue-specific constitutive overexpression of cytochrome P450s confers tolerance to imidacloprid in palm weevils in date palm fields

Background

Cytochrome P450-dependent monooxygenases (P450s), constituting one of the largest and oldest gene superfamilies found in many organisms from bacteria to humans, play a vital role in the detoxification and inactivation of endogenous toxic compounds. The use of various insecticides has increased over the last two decades, and insects have developed resistance to most of these compounds through the detoxifying function of P450s. In this study, we focused on the red palm weevil (RPW), Rhynchophorus ferrugineus, the most devastating pest of palm trees worldwide, and demonstrated through functional analysis that upregulation of P450 gene expression has evolved as an adaptation to insecticide stress arising from exposure to the neonicotinoid-class systematic insecticide imidacloprid.

Results

Based on the RPW global transcriptome analysis, we identified 101 putative P450 genes, including 77 likely encoding protein coding genes with ubiquitous expression. A phylogenetic analysis revealed extensive functional and species-specific diversification of RPW P450s, indicating that multiple CYPs actively participated in the detoxification process. We identified highly conserved paralogs of insect P450s that likely play a role in the development of resistance to imidacloprid: Drosophila Cyp6g1 (CYP6345J1) and Bemisia tabaci CYP4C64 (CYP4LE1). We performed a toxicity bioassay and evaluated the induction of P450s, followed by the identification of overexpressed P450s, including CYP9Z82, CYP6fra5, CYP6NR1, CYP6345J1 and CYP4BD4, which confer cross-resistance to imidacloprid. In addition, under imidacloprid insecticide stress in a date palm field, we observed increased expression of various P450 genes, with CYP9Z82, CYP4BD4, CYP6NR1 and CYP6345J1 being the most upregulated detoxification genes in RPWs. Expression profiling and cluster analysis revealed P450 genes with multiple patterns of induction and differential expression. Furthermore, we used RNA interference to knock down the overexpressed P450s, after which a toxicity bioassay and quantitative expression analysis revealed likely candidates involved in metabolic resistance against imidacloprid in RPW. Ingestion of double-stranded RNA (dsRNA) successfully knocked down the expression of CYP9Z82, CYP6NR1 and CYP345J1 and demonstrated that silencing of CYP345J1 and CYP6NR1 significantly decreased the survival rate of adult RPWs treated with imidacloprid, indicating that overexpression of these two P450s may play an important role in developing tolerance to imidacloprid in a date palm field.

Conclusion

Our study provides useful background information on imidacloprid-specific induction and overexpression of P450s, which may enable the development of diagnostic tools/markers for monitoring the spread of insecticide resistant RPWs. The observed trend of increasing tolerance to imidacloprid in the date palm field therefore indicated that strategies for resistance management are urgently needed.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5837-4) contains supplementary material, which is available to authorized users.

Expression and functional analysis of cytochrome P450 genes in the wolf spider Pardosa pseudoannulata under cadmium stress

Cytochrome P450 enzymes (CYPs), encoded by Halloween genes, mediate the biosynthesis of molting hormone, ecdysteroids, in arthropods. In this report, the effect of heavy metal cadmium (Cd) stress on the expression of cytochrome P450 genes in the wolf spider Pardosa pseudoannulata was analyzed. The results showed the expression levels of genes encoding for Cd transporters including ABC transporters, zinc transporters, calcium channel proteins and calcium binding proteins were inhibited or induced by Cd stress. In addition, the increase in metallothionein (MT) content and glutathione peroxidase (GPX) activity and decrease in total acetylcholine esterase (AChE) activity were also detected. Apparently, these detoxification methods did not completely protect the spider from the cytotoxicity of Cd stress. Increased mortality of P. pseudoannulata was observed when they were under Cd tress. In total 569 CYP genes belonging to 62 CYP subfamilies were obtained from P. pseudoannulata RNA-seq databases. BlaxtX analysis showed that 150, 161, 11, and 40 CYP genes were similar to the genes dib, phm, sad and shd, respectively, which are thought to catalyze the biosynthesis of ecdysteroids. Gene expression analysis suggested that 25 dib encoding genes, 27 phm encoding genes, 2 sad encoding genes, and 6 shd encoding genes were differentially expressed in TS2 vs. S2 comparison (Cd-treated 2nd instar spider vs. 2nd instar spider), respectively. There were 70 dib, 70 phm and 19 shd encoding genes either upregulated or downregulated, while 3 sad encoding genes were upregulated in TS5 vs. S5 (Cd-treated 5nd instar spider vs. 5nd instar spider). Genes related to heme binding and essential for activating the CYPs were also differentially expressed. Expression levels of cuticle related genes were significant differentially expressed, implying the changes in activities of chitin synthases and chitinase. Therefore we assume that unsuccessful molting process may occur on P. pseudoannulata due to influenced ecdysteroids levels, thus increasing mortality of spider.

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.

Variation in Mitochondria-Derived Transcript Levels Associated With DDT Resistance in the 91-R Strain of Drosophila melanogaster (Diptera: Drosophilidae)

The organochloride insecticide dichlorodiphenyltrichloroethane (DDT) and its metabolites can increase cellular levels of reactive oxygen species (ROS), cause mitochondrial dysfunction, and induce apoptosis. The highly DDT-resistant Drosophila melanogaster Meigen 1830 (Drosophila) strain, 91-R, and its susceptible control, 91-C, were used to investigate functional and structural changes among mitochondrial-derived pathways. Resequencing of mitochondrial genomes (mitogenomes) detected no structural differences between 91-R and 91-C, whereas RNA-seq suggested the differential expression of 221 mitochondrial-associated genes. Reverse transcriptase-quantitative PCR validation of 33 candidates confirmed that transcripts for six genes (Cyp12d1-p, Cyp12a4, cyt-c-d, COX5BL, COX7AL, CG17140) were significantly upregulated and two genes (Dif, Rel) were significantly downregulated in 91-R. Among the upregulated genes, four genes are duplicated within the reference genome (cyt-c-d, CG17140, COX5BL, and COX7AL). The predicted functions of the differentially expressed genes, or known functions of closely related genes, suggest that 91-R utilizes existing ROS regulation pathways of the mitochondria to combat increased ROS levels from exposure to DDT. This study represents, to our knowledge, the initial investigation of mitochondrial genome sequence variants and functional adaptations in responses to intense DDT selection and provides insights into potential adaptations of ROS management associated with DDT selection in Drosophila.

Transcriptional and Translational Relationship in Environmental Stress: RNAseq and ITRAQ Proteomic Analysis Between Sexually Reproducing and Parthenogenetic Females in Moina micrura

Moina micrura is a kind of small-bodied water flea within the family Moinidae. Similar to Daphnia, M. micrura could also switch its reproduction mode from parthenogenetic female (PF) to sexual female (SF) to adapt to the external environment. To uncover the mechanisms of reproductive switching in M. micrura, we used both RNA-Seq and iTRAQ analyses to investigate the differentially expressed genes (DEGs) and their protein products between SF and PF in M. micrura. A total of 1665 DEGs (702 up-regulated, 963 down-regulated) and 600 differentially expressed proteins (DEPs) (102 up-regulated, 498 down-regulated) were detected in SF. Correlation analyses indicated that 31 genes were expressed significantly differentially at both transcriptomic and proteomic levels, including 15 up-regulated genes and 16 down-regulated genes in SF. Meanwhile, our data also showed that 528 DEPs have discordant expression at transcript level, implying post-transcriptional (including translational) regulation. These top up-regulated genes and their protein products in SF were mainly grouped into the globin-related family, vitellogenin-related family, cuticle-related family, Hsp-related family and methyltransferases-related family, which were all involved in the reproductive switching in Daphnia. In contrast, a cluster of orthologous groups revealed that up-regulated genes and their protein products in PF were strongly associated with the metabolic process, which may be responsible for rapid population proliferation in M. micrura.

Transcriptome-based identification and characterization of genes commonly responding to five different insecticides in the diamondback moth, Plutella xylostella

When the 3rd instar larvae of the diamondback moth (DBM), Plutella xylostella, were pretreated with sublethal doses (LC₁₀) and then subsequently exposed to lethal doses (LC₅₀) of chlorantraniliprole, cypermethrin, dinotefuran, indoxacarb and spinosad via leaf dipping, their tolerance to insecticides was significantly enhanced. To identify genes that commonly respond to the treatment of different insecticides and are responsible for the tolerance enhancement, transcriptomic profiles of larvae treated with sublethal doses of the five insecticides were compared with that of untreated control. A total of 117,181 transcripts with a mean length of 662bp were generated by de novo assembly, of which 35,329 transcripts were annotated. Among them, 125, 143, 182, 215 and 149 transcripts were determined to be up-regulated whereas 67, 45, 60, 60 and 38 genes were down-regulated following treatments with chlorantraniliprole, cypermethrin, dinotefuran, indoxacarb and spinosad, respectively. Gene ontology (GO) analysis of differentially expressed genes (DEGs) revealed little differences in their GO profiles between treatments with different insecticides except for spinosad. Finally, the DEGs commonly responding to all insecticides were selected for further characterization, and some of their over-transcription levels were confirmed by quantitative PCR. The most notable examples of commonly responding over-transcribed genes were two cytochrome P450 genes (Cyp301a1 and Cyp9e2) and nine cuticular protein genes. In contrast, several genes composing the mitochondrial energy generation system were significantly down-regulated in all treated larvae. Considering the distinct structure and mode of action of the five insecticides tested, the differentially expressed genes identified in this study appear to be involved in general chemical defense at the initial stage of intoxication. Their possible roles in the tolerance/resistance development were discussed.

Induced thiacloprid insensitivity in honeybees (Apis mellifera L.) is associated with up-regulation of detoxification genes

Honey bees, Apis mellifera, are markedly less sensitive to neonicotinoid insecticides containing a cyanoimino pharmacophore than to those with a nitroimino group. Although previous work has suggested that this results from enhanced metabolism of the former by detoxification enzymes, the specific enzyme(s) involved remain to be characterized. In this work, a pretreatment of honey bees with a sublethal dose of thiacloprid resulted in induced insensitivity to the same compound immediately following thiacloprid feeding. A longer pretreatment time resulted in no, or increased, sensitivity. Transcriptome profiling, using microarrays, identified a number of genes encoding detoxification enzymes that were over-expressed significantly in insecticide-treated bees compared with untreated controls. These included five P450s, CYP6BE1, CYP305D1, CYP6AS5, CYP315A1, CYP301A1, and a carboxyl/cholinesterase (CCE) CCE8. Four of these P450s were functionally expressed in Escherichia coli and their ability to metabolize thiacloprid examined by liquid chromatography-mass spectrometry (LC-MS) analysis.

Isolation of the monooxygenase complex from Rhipicephalus (Boophilus) microplus - clues to understanding acaricide resistance

The monooxygenase complex is composed of three key proteins, a cytochrome P450 (CYP), the cytochrome P450 oxidoreductase (CPR) and cytochrome b5 and plays a key role in the metabolism and detoxification of xenobiotic substances, including pesticides. In addition, overexpression of these components has been linked to pesticide resistance in several important vectors of disease. Despite this, the monooxygenase complex has not been isolated from the Southern cattle tick Rhipicephalus (Boophilus) microplus, a major disease vector in livestock. Using bioinformatics 115 transcriptomic sequences were analyzed to identify putative pesticide metabolizing CYPs. RACE-PCR was used to amplify the full length sequence of one CYP; CYP3006G8 which displays a high degree of homology to members of the CYP6 and 9 subfamilies, known to metabolize pyrethroids. mRNA expression levels of CYP3006G8 were investigated in 11 strains of R. microplus with differing resistance profiles by qPCR, the results of which indicated a correlation with pyrethroid metabolic resistance. In addition to this gene, the sequences for CPR and cytochrome b5 were also identified and subsequently isolated from R. microplus using PCR. CYP3006G8 is only the third CYP gene isolated from R. microplus and the first to putatively metabolize pesticides. The initial results of expression analysis suggest that CYP3006G8 metabolizes pyrethroids but further biochemical characterization is required to confirm this. Differences in the kinetic parameters of human and mosquito CPR in terms of NADPH binding have been demonstrated and could potentially be used to design species specific pesticides. Similar differences in the tick CPR would confirm that this is a characteristic of heamatophagous arthropods.

Identifying Loci Contributing to Natural Variation in Xenobiotic Resistance in Drosophila

Natural populations exhibit a great deal of interindividual genetic variation in the response to toxins, exemplified by the variable clinical efficacy of pharmaceutical drugs in humans, and the evolution of pesticide resistant insects. Such variation can result from several phenomena, including variable metabolic detoxification of the xenobiotic, and differential sensitivity of the molecular target of the toxin. Our goal is to genetically dissect variation in the response to xenobiotics, and characterize naturally-segregating polymorphisms that modulate toxicity. Here, we use the Drosophila Synthetic Population Resource (DSPR), a multiparent advanced intercross panel of recombinant inbred lines, to identify QTL (Quantitative Trait Loci) underlying xenobiotic resistance, and employ caffeine as a model toxic compound. Phenotyping over 1,700 genotypes led to the identification of ten QTL, each explaining 4.5-14.4% of the broad-sense heritability for caffeine resistance. Four QTL harbor members of the cytochrome P450 family of detoxification enzymes, which represent strong a priori candidate genes. The case is especially strong for Cyp12d1, with multiple lines of evidence indicating the gene causally impacts caffeine resistance. Cyp12d1 is implicated by QTL mapped in both panels of DSPR RILs, is significantly upregulated in the presence of caffeine, and RNAi knockdown robustly decreases caffeine tolerance. Furthermore, copy number variation at Cyp12d1 is strongly associated with phenotype in the DSPR, with a trend in the same direction observed in the DGRP (Drosophila Genetic Reference Panel). No additional plausible causative polymorphisms were observed in a full genomewide association study in the DGRP, or in analyses restricted to QTL regions mapped in the DSPR. Just as in human populations, replicating modest-effect, naturally-segregating causative variants in an association study framework in flies will likely require very large sample sizes.

Genome of the house fly, Musca domestica L., a global vector of diseases with adaptations to a septic environment

Background

Adult house flies, Musca domestica L., are mechanical vectors of more than 100 devastating diseases that have severe consequences for human and animal health. House fly larvae play a vital role as decomposers of animal wastes, and thus live in intimate association with many animal pathogens.

Results

We have sequenced and analyzed the genome of the house fly using DNA from female flies. The sequenced genome is 691 Mb. Compared with Drosophila melanogaster, the genome contains a rich resource of shared and novel protein coding genes, a significantly higher amount of repetitive elements, and substantial increases in copy number and diversity of both the recognition and effector components of the immune system, consistent with life in a pathogen-rich environment. There are 146 P450 genes, plus 11 pseudogenes, in M. domestica, representing a significant increase relative to D. melanogaster and suggesting the presence of enhanced detoxification in house flies. Relative to D. melanogaster, M. domestica has also evolved an expanded repertoire of chemoreceptors and odorant binding proteins, many associated with gustation.

Conclusions

This represents the first genome sequence of an insect that lives in intimate association with abundant animal pathogens. The house fly genome provides a rich resource for enabling work on innovative methods of insect control, for understanding the mechanisms of insecticide resistance, genetic adaptation to high pathogen loads, and for exploring the basic biology of this important pest. The genome of this species will also serve as a close out-group to Drosophila in comparative genomic studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0466-3) contains supplementary material, which is available to authorized users.

Draft genome of the most devastating insect pest of coffee worldwide: the coffee berry borer, Hypothenemus hampei

The coffee berry borer, Hypothenemus hampei, is the most economically important insect pest of coffee worldwide. We present an analysis of the draft genome of the coffee berry borer, the third genome for a Coleopteran species. The genome size is ca. 163 Mb with 19,222 predicted protein-coding genes. Analysis was focused on genes involved in primary digestion as well as gene families involved in detoxification of plant defense molecules and insecticides, such as carboxylesterases, cytochrome P450, gluthathione S-transferases, ATP-binding cassette transporters, and a gene that confers resistance to the insecticide dieldrin. A broad range of enzymes capable of degrading complex polysaccharides were identified. We also evaluated the pathogen defense system and found homologs to antimicrobial genes reported in the Drosophila genome. Ten cases of horizontal gene transfer were identified with evidence for expression, integration into the H. hampei genome, and phylogenetic evidence that the sequences are more closely related to bacterial rather than eukaryotic genes. The draft genome analysis broadly expands our knowledge on the biology of a devastating tropical insect pest and suggests new pest management strategies.

Increased expression of CSP and CYP genes in adult silkworm females exposed to avermectins

We analyzed 20 chemosensory protein (CSP) genes of the silkworm Bombyx mori. We found a high number of retrotransposons inserted in introns. We then analyzed expression of the 20 BmorCSP genes across tissues using quantitative real-time polymerase chain reaction (PCR). Relatively low expression levels of BmorCSPs were found in the gut and fat body tissues. We thus tested the effects of endectocyte insecticide abamectin (B1a and B1b avermectins) on BmorCSP gene expression. Quantitative real-time PCR experiments showed that a single brief exposure to insecticide abamectin increased dramatically CSP expression not only in the antennae but in most tissues, including gut and fat body. Furthermore, our study showed coordinate expression of CSPs and metabolic cytochrome P450 enzymes in a tissue-dependent manner in response to the insecticide. The function of CSPs remains unknown. Based on our results, we suggest a role in detecting xenobiotics that are then detoxified by cytochrome P450 anti-xenobiotic enzymes.

Characterization and expression of the cytochrome P450 gene family in diamondback moth, Plutella xylostella (L.)

Cytochrome P450 monooxygenases are present in almost all organisms and can play vital roles in hormone regulation, metabolism of xenobiotics and in biosynthesis or inactivation of endogenous compounds. In the present study, a genome-wide approach was used to identify and analyze the P450 gene family of diamondback moth, Plutella xylostella, a destructive worldwide pest of cruciferous crops. We identified 85 putative cytochrome P450 genes from the P. xylostella genome, including 84 functional genes and 1 pseudogene. These genes were classified into 26 families and 52 subfamilies. A phylogenetic tree constructed with three additional insect species shows extensive gene expansions of P. xylostella P450 genes from clans 3 and 4. Gene expression of cytochrome P450s was quantified across multiple developmental stages (egg, larva, pupa and adult) and tissues (head and midgut) using P. xylostella strains susceptible or resistant to insecticides chlorpyrifos and fiprinol. Expression of the lepidopteran specific CYP367s predominantly occurred in head tissue suggesting a role in either olfaction or detoxification. CYP340s with abundant transposable elements and relatively high expression in the midgut probably contribute to the detoxification of insecticides or plant toxins in P. xylostella. This study will facilitate future functional studies of the P. xylostella P450s in detoxification.

A de novo transcriptome of European pollen beetle populations and its analysis, with special reference to insecticide action and resistance

The pollen beetle Meligethes aeneus is the most important coleopteran pest in European oilseed rape cultivation, annually infesting millions of hectares and responsible for substantial yield losses if not kept under economic damage thresholds. This species is primarily controlled with insecticides but has recently developed high levels of resistance to the pyrethroid class. The aim of the present study was to provide a transcriptomic resource to investigate mechanisms of resistance. cDNA was sequenced on both Roche (Indianapolis, IN, USA) and Illumina (LGC Genomics, Berlin, Germany) platforms, resulting in a total of ∼53 m reads which assembled into 43 396 expressed sequence tags (ESTs). Manual annotation revealed good coverage of genes encoding insecticide target sites and detoxification enzymes. A total of 77 nonredundant cytochrome P450 genes were identified. Mapping of Illumina RNAseq sequences (from susceptible and pyrethroid-resistant strains) against the reference transcriptome identified a cytochrome P450 (CYP6BQ23) as highly overexpressed in pyrethroid resistance strains. Single-nucleotide polymorphism analysis confirmed the presence of a target-site resistance mutation (L1014F) in the voltage-gated sodium channel of one resistant strain. Our results provide new insights into the important genes associated with pyrethroid resistance in M. aeneus. Furthermore, a comprehensive EST resource is provided for future studies on insecticide modes of action and resistance mechanisms in pollen beetle.

The molecular evolution of cytochrome P450 genes within and between drosophila species

We map 114 gene gains and 74 gene losses in the P450 gene family across the phylogeny of 12 Drosophila species by examining the congruence of gene trees and species trees. Although the number of P450 genes varies from 74 to 94 in the species examined, we infer that there were at least 77 P450 genes in the ancestral Drosophila genome. One of the most striking observations in the data set is the elevated loss of P450 genes in the Drosophila sechellia lineage. The gain and loss events are not evenly distributed among the P450 genes-with 30 genes showing no gene gains or losses whereas others show as many as 20 copy number changes among the species examined. The P450 gene clades showing the fewest number of gene gain and loss events tend to be those evolving with the most purifying selection acting on the protein sequences, although there are exceptions, such as the rapid rate of amino acid replacement observed in the single copy phantom (Cyp306a1) gene. Within D. melanogaster, we observe gene copy number polymorphism in ten P450 genes including multiple cases of interparalog chimeras. Nonallelic homologous recombination (NAHR) has been associated with deleterious mutations in humans, but here we provide a second possible example of an NAHR event in insect P450s being adaptive. Specifically, we find that a polymorphic Cyp12a4/Cyp12a5 chimera correlates with resistance to an insecticide. Although we observe such interparalog exchange in our within-species data sets, we have little evidence of it between species, raising the possibility that such events may occur more frequently than appreciated but are masked by subsequent sequence change.

Crustacean oxi-reductases protein sequences derived from a functional genomic project potentially involved in ecdysteroid hormones metabolism - a starting point for function examination

A transcriptomic assembly originated from hypodermis and Y organ of the crustacean Pontastacus leptodactylus is used here for in silico characterization of oxi-reductase enzymes potentially involved in the metabolism of ecdysteroid molting hormones. RNA samples were extracted from male Y organ and its neighboring hypodermis in all stages of the molt cycle. An equimolar RNA mix from all stages was sequenced using next generation sequencing technologies and de novo assembled, resulting with 74,877 unique contigs. These transcript sequences were annotated by examining their resemblance to all GenBank translated transcripts, determining their Gene Ontology terms and their characterizing domains. Based on the present knowledge of arthropod ecdysteroid metabolism and more generally on steroid metabolism in other taxa, transcripts potentially related to ecdysteroid metabolism were identified and their longest possible conceptual protein sequences were constructed in two stages, correct reading frame was deduced from BLASTX resemblances, followed by elongation of the protein sequence by identifying the correct translation frame of the original transcript. The analyzed genes belonged to several oxi-reductase superfamilies including the Rieske non heme iron oxygenases, cytochrome P450s, short-chained hydroxysteroid oxi-reductases, aldo/keto oxireductases, lamin B receptor/sterol reductases and glucose-methanol-cholin oxi-reductatses. A total of 68 proteins were characterized and the most probable participants in the ecdysteroid metabolism where indicated. The study provides transcript and protein structural information, a starting point for further functional studies, using a variety of gene-specific methods to demonstrate or disprove the roles of these proteins in relation to ecdysteroid metabolism in P. leptodactylus.

De novo characterization of transcriptome and gene expression dynamics in epidermis during the larval-pupal metamorphosis of common cutworm

Larval cuticle is degraded and replaced by the pupal counterpart during larval-pupal metamorphosis in the holometabolous insects. In addition to the extrinsic transformation, the epidermis goes through significant changes at molecular levels. To elucidate the intrinsic mechanism of epidermal metamorphosis, the dynamics of chitin content in the cuticle was examined in an important agricultural lepidopteran, the common cutworm, and the transcriptome was analyzed using Illumina sequencing technology. Gene expression profiles during the metamorphosis were further studied by both the digital gene expression (DGE) system and real-time quantitative PCR. The results showed that the chitin content decreased in prepupae and then increased in pupae. A total of 58 million sequencing reads were obtained and assembled into 70,346 unigenes. Over 9000 unigenes were identified to express differentially during the transformation process. As compared with the 6th instar feeding larvae, the most significant changes took place in the proteasome and metabolic pathways in prepupae and pupae, respectively. The cytochrome P450s, VHDLs, chitinase, serine protease and genes involved in sex pheromone biosynthesis changed their mRNA levels remarkably. Three chitinolytic enzymes (chitinase, β-N-acetylglucosaminidase and chitin deacetylase) showed distinct mRNA expression patterns, the former two enzymes revealed the highest expression in prepupae, however the latter one showed its climax mRNA level in pupae. The gene expression patterns suggest that chitinase and β-N-acetylglucosaminidase may be responsible for the degradation of larval cuticles, whereas chitin deacetylase may help to degrade the pupal counterparts. Gene expression dynamics also implied that the chitin of pupal cuticle might be formed by recycling of the degraded chitin of larval cuticle rather than through de novo synthesis. The 20E-induced nuclear receptors seem to be important factors regulating chitin metabolic enzymes during the cuticle remodeling. Our data provide a comprehensive resource for exploring the molecular mechanism of epidermal metamorphosis in insects.

dusky-like is required to maintain the integrity and planar cell polarity of hairs during the development of the Drosophila wing

The cuticular hairs and sensory bristles that decorate the adult Drosophila epidermis and the denticles found on the embryo have been used in studies on planar cell polarity and as models for the cytoskeletal mediated morphogenesis of cellular extensions. ZP domain proteins have recently been found to be important for the morphogenesis of both denticles and bristles. Here we show that the ZP domain protein Dusky-like is a key player in hair morphogenesis. As is the case in bristles, in hairs dyl mutants display a dramatic phenotype that is the consequence of a failure to maintain the integrity of the extension after outgrowth. Hairs lacking dyl function are split, thinned, multipled and often very short. dyl is required for normal chitin deposition in hairs, but chitin is not required for the normal accumulation of Dyl, hence dyl acts upstream of chitin. A lack of chitin however, does not mimic the dyl hair phenotype, thus Dyl must have other targets in hair morphogenesis. One of these appears to be the actin cytoskeleton. Interestingly, dyl mutants also display a unique planar cell polarity phenotype that is distinct from that seen with mutations in the frizzled/starry night or dachsous/fat pathway genes. Rab11 was previously found to be essential for Dyl plasma membrane localization in bristles. Here we found that the expression of a dominant negative Rab11 can mimic the dyl hair morphology phenotype consistent with Rab11 also being required for Dyl function in hairs. We carried out a small directed screen to identify genes that might function with dyl and identified Chitinase 6 (Cht6) as a strong candidate, as knocking down Cht6 function led to weak versions of all of the dyl hair phenotypes.

Draft genome of the mountain pine beetle, Dendroctonus ponderosae Hopkins, a major forest pest

BACKGROUND

The mountain pine beetle, Dendroctonus ponderosae Hopkins, is the most serious insect pest of western North American pine forests. A recent outbreak destroyed more than 15 million hectares of pine forests, with major environmental effects on forest health, and economic effects on the forest industry. The outbreak has in part been driven by climate change, and will contribute to increased carbon emissions through decaying forests.

RESULTS

We developed a genome sequence resource for the mountain pine beetle to better understand the unique aspects of this insect's biology. A draft de novo genome sequence was assembled from paired-end, short-read sequences from an individual field-collected male pupa, and scaffolded using mate-paired, short-read genomic sequences from pooled field-collected pupae, paired-end short-insert whole-transcriptome shotgun sequencing reads of mRNA from adult beetle tissues, and paired-end Sanger EST sequences from various life stages. We describe the cytochrome P450, glutathione S-transferase, and plant cell wall-degrading enzyme gene families important to the survival of the mountain pine beetle in its harsh and nutrient-poor host environment, and examine genome-wide single-nucleotide polymorphism variation. A horizontally transferred bacterial sucrose-6-phosphate hydrolase was evident in the genome, and its tissue-specific transcription suggests a functional role for this beetle.

CONCLUSIONS

Despite Coleoptera being the largest insect order with over 400,000 described species, including many agricultural and forest pest species, this is only the second genome sequence reported in Coleoptera, and will provide an important resource for the Curculionoidea and other insects.