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

De novo transcriptome assemblies of five major European oilseed rape insect pests

OBJECTIVE

Insect pests can cause severe losses in oilseed rape yields across Europe. Genomic and transcriptomic information is very limited for these insects. The aim of our study was to provide transcriptomic resources on several oilseed rape herbivores that will support research into their biology and help develop new methods of sustainable pest management.

DATA

Transcriptomes for larval stages of five major European pest species were de novo assembled by Trinity assembler. Total number of transcripts ranged from 112,247 for Ceutorhynchus pallidactylus to 225,110 for Ceutorhyncus napi. Intermediate numbers of 140,588, 140,998 and 144,504, were found for Psylliodes chrysocephala, Dasineura brassicae, and Brassicogethes aeneus, respectively. Bench-marking universal single-copy orthologues analyses for each dataset indicated high degree of completeness for all five species. The transcriptomes extend the list of genomic data on insect larvae that constitute major pests of oilseed rape. The data provide information on larval physiology and form a basis to develop highly specific RNA interference-based plant protection.

Transcriptomics and metagenomics of common cutworm (Spodoptera litura) and fall armyworm (Spodoptera frugiperda) demonstrate differences in detoxification and development

BACKGROUND

Spodoptera litura is an important polyphagous pest that causes significant damage to the agricultural sector. We performed RNA-seq of 15 S. litura individuals from larval (fifth and sixth instar larvae), chrysalis, and adult developmental stages. We also compared the S. litura transcriptome data with Spodoptera frugiperda across the same developmental stages, which was sequenced in our previous study.

RESULTS

A total of 101,885 differentially expressed transcripts (DETs) were identified in S. litura. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that S. litura may undergo active xenobiotic and detoxifying metabolism during its larval and adult stages, which may explain difficulties with current population control measures. We also found that DETs of single-copy orthologous genes between S. litura and S. frugiperda were involved in basic metabolism and development. However, energy and metabolic processes genes had a higher expression in S. litura, whereas nervous and olfactory function genes had a higher expression in S. frugiperda. Metagenomics analysis in larval S. litura and S. frugiperda revealed that microbiota participate in the detoxification and metabolism processes, but the relative abundance of detoxification-related microbiota was more abundant in S. frugiperda. Transcriptome results also confirmed the detoxification-related pathway of S. frugiperda was more abundant than in S. litura.

CONCLUSIONS

Significant changes at transcriptional level were identified during the different development stages of S. litura. Importantly, we also identified detoxification associated genes and gut microbiota between S. litura and S. frugiperda at different developmental stages, which will be valuable in revealing possible mechanisms of detoxification and development in these two lepidopterans.

Transcriptomic profile of the predatory mite Amblyseius swirskii (Acari: Phytoseiidae) on different host plants

Amblyseius swirskii Athias-Henriot (Acari: Phytoseiidae) is a predatory mite, effective at controlling whiteflies and thrips in protected crops. However, on tomato its efficacy as a biocontrol agent is hindered, most probably by the plant trichomes and their exudates. Our aim was to characterize the response of A. swirskii to the tomato trichome exudates and identify three major detoxification gene sets in this species: cytochromes P450 (CYPs), glutathione S-transferases (GSTs) and carboxyl/cholinesterases (CCEs). Mites were exposed separately to tomato and pepper, a favourable host plant for A. swirskii, after which their transcriptional responses were analysed and compared. The de novo transcriptome assembly resulted in 71,336 unigenes with 66.1% of them annotated. Thirty-nine A. swirskii genes were differentially expressed after transfer on tomato leaves when compared to pepper leaves; some of the expressed genes were associated with the metabolism of tomato exudates. Our results illustrate that the detoxification gene sets CYPs, GSTs and CCEs are abundant in A. swirskii, but do not play a significant role when in contact with the tomato exudates.

Bacillus thuringiensis and Chlorantraniliprole Trigger the Expression of Detoxification-Related Genes in the Larval Midgut of Plutella xylostella

Background: Diamondback moth (DBM), Plutella xylostella (L.), has developed resistance to many insecticides. The molecular mechanism of DBM resistance to Bt-G033A combined with chlorantraniliprole (CL) remains undefined. Methods: In this study, field-resistant strains of Plutella xylostella to three pesticides, namely, Bacillus thuringiensis (Bt) toxin (Bt-G033A), CL, and a mixture of Bt + CL, were selected to evaluate the resistance level. Additionally, transcriptomic profiles of a susceptible (SS-DBM), field-resistant (FOH-DBM), Bt-resistant (Bt-DBM), CL-resistant (CL-DBM), and Bt + CL-resistant (BtC-DBM) strains were performed by comparative analysis to identify genes responsible for detoxification. Results: The Bt-G033A was the most toxic chemical to all the DBM strains among the three insecticides. The comparative analysis identified 25,518 differentially expressed genes (DEGs) between pairs/combinations of strains. DEGs were enriched in pathways related to metabolic and catalytic activity and ABC transporter in resistant strains. In total, 17 metabolic resistance-related candidate genes were identified in resistance to Bt-G033A, CL, and Bt + CL by co-expression network analysis. Within candidate genes, the majority was upregulated in key genes including cytochrome P450, glutathione S-transferase (GST), carboxylesterase, and acetylcholinesterase in CL- and BtC-resistant strains. Furthermore, aminopeptidase N (APN), alkaline phosphatase (ALP), cadherin, trypsin, and ABC transporter genes were eminent as Bt-resistance-related genes. Expression patterns of key genes by the quantitative real-time PCR (qRT-PCR) proved the credibility of transcriptome data and suggest their association in the detoxification process. Conclusion: To date, this study is the most comprehensive research presenting functional transcriptome analysis of DBM using Bt-G033A and CL combined insecticidal activity.

High genetic diversity and strong genetic structure of Strongyllodes variegatus populations in oilseed rape production areas of China

Background

Strongyllodes variegatus (Fairmaire) is a major insect pest of oilseed rape in China. Despite its economic importance, the contribution of its population genetics in the development of any suitable protection control strategy for the management of oilseed rape crops is poorly studied. It is a much urgent need to prevent its spread to the rest of the world.

Results

Using the sequences of mitochondrial DNA cytochrome c oxidase subunit I (COI) and cytochrome b (Cytb) as genetic markers, we analyzed the population genetic diversity and structure of 437 individuals collected from 15 S. variegatus populations located in different oilseed rape production areas in China. In addition, we estimated the demographic history using neutrality test and mismatch distribution analysis. The high level of genetic diversity was detected among the COI and Cytb sequences of S. variegatus. The population structure analyses strongly suggested three distinct genetic and geographical regions in China with limited gene flow. The Mantel test showed that the genetic distance was greatly influenced by the geographical distance. The demographic analyses showed that S. variegatus had experienced population fluctuation during the Pleistocene Epoch, which was likely to be related to the climatic changes.

Conclusion

Overall, these results demonstrate that the strong genetic structure of S. variegatus populations in China, which is attributed by the isolation through the geographical distance among populations, their weak flight capacity and subsequent adaptation to the regional ecological conditions.

Proteogenomic insight into the basis of the insecticide tolerance/resistance of the pollen beetle Brassicogethes (Meligethes) aeneus

The pollen beetle is a major pest of oilseed rape. Although various resistance mechanisms have been identified, such as kdr (mutation in the sodium channel) and metabolic resistance (CYP overexpression), other "hidden" factors also exist. Some studies have stressed the importance of epistasis as a genetic background. The combination of kdr and metabolic resistance appears to be unfavorable under field conditions in the absence of pesticide selection. The regulation of detoxification enzymes can play an important role, but we highlight different detoxification markers compared to those emphasized in other studies. We also stress the importance of studying the role of markers identified as pathogenesis-related protein 5-like (PR5; upregulated by insecticides) and highlight the role of RNA (DEAD-box) helicases (downregulated by insecticides). Thus, we suggest the importance of epigenetic drivers of resistance/tolerance to pesticides. The key results are similar to those of our previous study, in which deltamethrin treatment of the pollen beetle was also investigated by a proteogenomic approach. Indeed, the mechanism leading to resistance of the pollen beetle may be an innate mechanism that the pollen beetle can also employ in natural habitats, but under field conditions (pesticide exposure), this mechanism is used to survive in response to insecticides. SIGNIFICANCE: Pesticide resistance is a serious problem that hampers the successful production of crops. Understanding the mechanisms of insecticide resistance is highly important for successful pest control, especially when considering integrated pest management. Here, using a proteogenomic approach, we identified novel markers for understanding pollen beetle resistance to pesticides. In addition, future studies will reveal the role of these markers in the multiresistance of pollen beetle populations. We highlight that the proteins identified as PR5, which are known to occur in beetles and are similar to those in plants, may be responsible for tolerance to multiple stresses. In addition, our results indicate that the RNA helicases that exhibited changes in expression may be the epigenetic drivers of multiresistance. The nature of these changes remains an open question, and their relevance in different situations (responses to different stresses) in natural habitats in the absence of pesticides can be proposed.

In silico identification and assessment of insecticide target sites in the genome of the small hive beetle, Aethina tumida

BACKGROUND

The small hive beetle, Aethina tumida, is a rapidly emerging global pest of honey bee colonies. Small hive beetle infestation can be extremely destructive, which may cause honey bees to abscond and render colony infrastructure unusable. Due to the impacts small hive beetles have on honey bees, a wide variety of physical, cultural, and chemical control measures have been implemented to manage small hive beetle infestations. The use of insecticides to control small hive beetle populations is an emerging management tactic. Currently, very little genomic information exists on insecticide target sites in the small hive beetle. Therefore, the objective of this study is to utilize focused in silico comparative genomics approaches to identify and assess the potential insecticide sensitivity of the major insecticide target sites in the small hive beetle genome.

RESULTS

No previously described resistance mutations were identified in any orthologs of insecticide target sites. Alternative exon use and A-to-I RNA editing were absent in AtumSC1. The ryanodine receptor in small hive beetle (Atum_Ryr) was highly conserved and no previously described resistance mutations were identified. A total of 12 nAChR subunits were identified with similar alternative exon use in other insects. Alternative exon use and critical structural features of the GABA-gated chloride channel subunits (Atum_RDL, Atum_GRD, and Atum_LCCH3) were conserved. Five splice variants were found for the glutamate-gated chloride channel subunit. Exon 3c of Atum_GluCl may be a beetle-specific alternative exon. The co-occurrence of exons 9a and 9b in the pH-sensitive chloride channel (Atum_pHCl) is a unique combination that introduces sites of post-translational modification. The repertoire and alternative exon use for histamine-gated chloride channels (Atum-HisCl), octopamine (Atum_OctR) and tyramine receptors (Atum_TAR) were conserved.

CONCLUSIONS

The recently published small hive beetle genome likely serves as a reference for insecticide-susceptible versions of insecticide target sites. These comparative in silico studies are the first step in discovering targets that can be exploited for small hive beetle-specific control as well as tracking changes in the frequency of resistance alleles as part of a resistance monitoring program. Comparative toxicity alongside honey bees is required to verify these in silico predictions.

Identification of transcriptome and fluralaner responsive genes in the common cutworm Spodoptera litura Fabricius, based on RNA-seq

Background

Fluralaner is a novel isoxazoline insecticide with a unique action site on the γ-aminobutyric acid receptor (GABAR), shows excellent activity on agricultural pests including the common cutworm Spodoptera litura, and significantly influences the development and fecundity of S. litura at either lethal or sublethal doses. Herein, Illumina HiSeq Xten (IHX) platform was used to explore the transcriptome of S. litura and to identify genes responding to fluralaner exposure.

Results

A total of 16,572 genes, including 451 newly identified genes, were observed in the S. litura transcriptome and annotated according to the COG, GO, KEGG and NR databases. These genes included 156 detoxification enzyme genes [107 cytochrome P450 enzymes (P450s), 30 glutathione S-transferases (GSTs) and 19 carboxylesterases (CarEs)] and 24 insecticide-targeted genes [5 ionotropic GABARs, 1 glutamate-gated chloride channel (GluCl), 2 voltage-gated sodium channels (VGSCs), 13 nicotinic acetylcholine receptors (nAChRs), 2 acetylcholinesterases (AChEs) and 1 ryanodine receptor (RyR)]. There were 3275 and 2491 differentially expressed genes (DEGs) in S. litura treated with LC₃₀ or LC₅₀ concentrations of fluralaner, respectively. Among the DEGs, 20 related to detoxification [16 P450s, 1 GST and 3 CarEs] and 5 were growth-related genes (1 chitin and 4 juvenile hormone synthesis genes). For 26 randomly selected DEGs, real-time quantitative PCR (RT-qPCR) results showed that the relative expression levels of genes encoding several P450s, GSTs, heat shock protein (HSP) 68, vacuolar protein sorting-associated protein 13 (VPSAP13), sodium-coupled monocarboxylate transporter 1 (SCMT1), pupal cuticle protein (PCP), protein takeout (PT) and low density lipoprotein receptor adapter protein 1-B (LDLRAP1-B) were significantly up-regulated. Conversely, genes encoding esterase, sulfotransferase 1C4, proton-coupled folate transporter, chitinase 10, gelsolin-related protein of 125 kDa (GRP), fibroin heavy chain (FHC), fatty acid synthase and some P450s were significantly down-regulated in response to fluralaner.

Conclusions

The transcriptome in this study provides more effective resources for the further study of S. litura whilst the DEGs identified sheds further light on the molecular response to fluralaner.

Limited genetic structure and demographic expansion of the Brassicogethes aeneus populations in France and in Europe

BACKGROUND

The pollen beetle, Brassicogethes aeneus (Fabricius, 1775), is one of the most significant pests of oilseed rape. To shed light on past and current pollen beetle demography (dispersal, population size), 12 microsatellite markers were developed, and population genetic diversity and structure were analysed at different spatial scales in France and in Europe from 433 individuals collected in 18 winter oilseed rape fields.

RESULTS

Genetic differentiation among the population samples was low but was significant between the Estonian sample and the rest of Europe. Isolation by distance was significant only at the European scale. Genetic variability was similar among the 18 population samples. Demographic inferences suggested a recent expansion of B. aeneus population size over Europe, possibly corresponding to an increase in oilseed rape crop area during past decades.

CONCLUSION

Current population size and dispersal are not straightforward to estimate from the distribution of genetic variability in B. aeneus over Europe because of the complexity of the demographic history of this pest. Nevertheless, because gene flow was important enough to prevent strong genetic differentiation at large geographical scales, the management of pollen beetle populations should likely be thought of at a continental Europe level. © 2018 Society of Chemical Industry.

Detection of tau-fluvalinate resistance in the mite Varroa destructor based on the comparison of vial test and PCR-RFLP of kdr mutation in sodium channel gene

Varroa destructor is the major cause of honey bee (Apis mellifera) colony losses. Mite control is limited to several miticides. The overuse of tau-fluvalinate has resulted in resistance via a knockdown resistance (kdr) mutation in the sodium channel gene Na_VChs (L925V/I/M). In this study, we used the discriminating concentration of tau-fluvalinate (0.25 µg/mL) to detect the resistance of mites in a bioassay. Further, we verified the presence of the kdr mutation in mites from the bioassay via PCR amplification of a fragment of the voltage-gated sodium channel gene (Na_VCh), restriction fragment length polymorphisms (RFLPs), and densitometry analyses in pools of surviving or dead mites. Resistance values corresponding to the densitometry of the resistant allele were related to mite survival. In the vial test, the survival of the control group was significantly higher (70.4%) than that of the tau-fluvalinate-treated group (34.3%). Mite survival in the vial test was significantly correlated with the mean proportion of resistance values. Individuals that died after tau-fluvalinate application exhibited an average resistance value of 0.0783, whereas individuals that survived exhibited an average resistance of 0.400. The concentration of tau-fluvalinate in the vials was checked using high performance liquid chromatography under different temperatures and exposure times, and indicates that the stability of tau-fluvalinate stored in the refrigerator (4 ± 1 °C) is at least 14 days. PCR-RFLP of the Na_VCh gene fragment verified that the vial test is a suitable, rapid, and cost-effective method for the identification of tau-fluvalinate resistance based on kdr mutation in V. destructor in apiaries.

Seven-year monitoring of pyrethroid resistance in the pollen beetle (Brassicogethes aeneus F.) during implementation of insect resistance management

BACKGROUND

An increase in the spread of pyrethroid resistance in the pollen beetle has been documented in many European countries. Pyrethroid resistance in the pollen beetle in the Czech Republic has been detected using a topical application bioassay. Resistance monitoring has been carried out during the implementation of insect resistance management.

RESULTS

The susceptibilities of nine pollen beetle populations to four pyrethroids in 2009 - 2015 in the Czech Republic are presented in this paper. The highest resistance ratio (RR) values [based on the lethal dose for 50% of the population (LD₅₀ )] obtained for deltamethrin, lambda-cyhalothrin, tau-fluvalinate and etofenprox were 500, 299, 108 and 66.9, respectively. Pollen beetle mortality after application of deltamethrin or lambda-cyhalothrin gradually decreased from 2009 to 2013. High cross-resistance between lambda-cyhalothrin and deltamethrin and low cross-resistance between lambda-cyhalothrin and etofenprox were demonstrated. A kdr mutation known to cause resistance in the pollen beetle was not detected.

CONCLUSION

A high level of resistance was recorded in most of the nine populations of the pollen beetle in 2013. In the following 2 years, a decrease in resistance associated with a population density decrease of the pollen beetle was observed in all nine populations as a result of insect resistance management based on pyrethroid reduction in winter oilseed rape. © 2017 Society of Chemical Industry.

Transcriptome analysis of the threatened snail Ellobium chinense reveals candidate genes for adaptation and identifies SSRs for conservation genetics

Ellobium chinense (Pfeiffer, 1854) is a brackish pulmonate species that inhabits the bases of mangrove trees and is most commonly found in salt grass meadows. Threats to mangrove ecosystems due to habitat degradation and overexploitation have threatened the species with extinction. In South Korea, E. chinense has been assessed as vulnerable, but there are limited data on its population structure and distribution. The nucleotide and protein sequences for this species are not available in databases, which limits the understanding of adaptation-related traits. We sequenced an E. chinense cDNA library using the Illumina platform, and the subsequent bioinformatics analysis yielded 227,032 unigenes. Of these unigenes, 69,088 were annotated to matched protein and nucleotide sequences in databases, for an annotation rate of 30.42%. Among the predominant gene ontology terms, cellular and metabolic processes (under the biological process category), membrane and cell (under the cellular component category), and binding and catalytic activity (under the molecular function category) were noteworthy. In addition, 4850 unigenes were distributed to 15 Kyoto Encyclopaedia of Genes and Genomes based enrichment categories. Among the candidate genes related to adaptation, angiotensin I converting enzyme, adenylate cyclase activating polypeptide, and AMP-activated protein kinase were the most prominent. A total of 15,952 simple sequence repeats (SSRs) were identified in sequences of > 1 kb in length. The di- and trinucleotide repeat motifs were the most common. Among the repeat motif types, AG/CT, AC/GT, and AAC/GTT dominated. Our study provides the first comprehensive genomics dataset for E. chinense, which favors conservation programs for the restoration of the species and provides sufficient evidence for genetic variability among the wild populations.

Comparative transcriptome analysis reveals host-associated differentiation in Chilo suppressalis (Lepidoptera: Crambidae)

The striped stem borer, Chilo suppressalis Walker (Lepidoptera: Crambidae), is one of the most serious rice pests. Besides attacking rice, it also feeds on an economically important vegetable crop, water-oat Zizania latifolia. The species feeding on water-oat has higher growth and survival rate than those on rice, suggesting their success in adaptation to the new host plant. However, little is known about the molecular mechanisms of host plant adaptation. Here we investigated the midgut transcriptome responses of C. suppressalis larvae reared on rice and water-oat. A total of 1,633 differentially expressed genes were identified, with a greater number up-regulated on the more delicious new host. The up-regulation of most digestive and detoxification-related genes may be the result of adaptation to the changes in nutritional requirements and toxic chemicals during host shift. In contrast, down-regulation of ribosomal genes may be related to their better development performance when feeding on the new host. In conclusion, our results suggest that transcriptional regulation of genes related to digestion, detoxification and ribosome may play an important role in adaptation of C. suppressalis to a new host plant.

Characterization of transcriptome in the Indian meal moth Plodia interpunctella (Lepidoptera: Pyralidae) and gene expression analysis during developmental stages

The Indian meal moth Plodia interpunctella (Lepidoptera: Pyralidae) is a worldwide pest that causes serious damage to stored foods. Although many efforts have been conducted on this species due to its economic importance, the study of genetic basis of development, behavior and insecticide resistance has been greatly hampered due to lack of genomic information. In this study, we used high throughput sequencing platform to perform a de novo transcriptome assembly and tag-based digital gene expression profiling (DGE) analyses across four different developmental stages of P. interpunctella (egg, third-instar larvae, pupae and adult). We obtained approximate 9gigabyte (GB) of clean data and recovered 84,938 unigenes, including 37,602 clusters and 47,336 singletons. These unigenes were annotated using BLAST against the non-redundant protein databases and then functionally classified based on Gene Ontology (GO), Clusters of Orthologous Groups (COG), and Kyoto Encyclopedia of Genes and Genomes databases (KEGG). A large number of differentially expressed genes were identified by pairwise comparisons among different developmental stages. Gene expression profiles dramatically changed between developmental stage transitions. Some of these differentially expressed genes were related to digestion and cuticularization. Quantitative real-time PCR results of six randomly selected genes conformed the findings in the DGEs. Furthermore, we identified over 8000 microsatellite markers and 97,648 single nucleotide polymorphisms which will be useful for population genetics studies of P. interpunctella. This transcriptomic information provided insight into the developmental basis of P. interpunctella and will be helpful for establishing integrated management strategies and developing new targets of insecticides for this serious pest.

De novo transcriptome and expression profile analyses of the Asian corn borer (Ostrinia furnacalis) reveals relevant flubendiamide response genes

Background

The Asian corn borer (ACB), Ostrinia furnacalis (Guenée), has become the most damaging insect pest of corn in Asia. However, the lack of genome or transcriptome information heavily hinders our further understanding of ACB in every aspect at a molecular level and on a genome-wide scale. Here, we used the Ion Torrent Personal Genome Machine (PGM) Sequencer to explore the ACB transcriptome and to identify relevant genes in response to flubendiamide, showing high selective activity against ACB.

Results

We obtained 35,430 unigenes, with an average length of 716 bp, representing a dramatic expansion of existing cDNA sequences available for ACB. These sequences were annotated with Non-redundant Protein (Nr), Gene Ontology (GO), Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) to better understand their functions. A total of 31 cytochrome P450 monooxygenases (P450s), 27 carboxyl/cholinesterases (CCEs) and 19 glutathione S-transferases (GSTs) were manually curated to construct phylogenetic trees, and 25 unigenes encoding target proteins (acetylcholinesterase, nicotinic acetylcholine receptor, gamma-aminobutyric acid receptor, glutamate-gated chloride channel, voltage-gated sodium channel and ryanodine receptor) were identified. In addition, we compared and validated the differentially expressed unigenes upon flubendiamide treatment, revealing that the genes for detoxification enzymes (P450s and esterase), calcium signaling pathways and muscle control pathways (twitchin and tropomyosin), immunoglobulin (hemolin), chemosensory protein and heat shock protein 70 were significantly overexpressed in response to flubendiamide, while the genes for cuticular protein, protease and oxidoreductase showed much lower expression levels.

Conclusion

The obtained transcriptome information provides large genomic resources available for further studies of ACB. The differentially expressed gene data will elucidate the molecular mechanisms of ACB in response to the novel diamide insecticide, flubendiamide. In particular, these findings will facilitate the identification of the genes involved in insecticide resistance and the development of new compounds to control the ACB.

Electronic supplementary material

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

The search for novel insecticide targets in the post-genomics era, with a specific focus on G-protein coupled receptors

Insects are considered pests globally, implicated in the destruction of agricultural fields and transmission of pathogens that cause deadly human diseases, such as dengue, Zika and malaria. The diversity of the insecticide arsenal has remained stagnant for decades, but the recent rise of insecticide resistance fueled the discovery of novel modes of action, and the power of genomics has reinvigorated this search. This review discusses the importance of comparative and functional insect genomics in the identification of potential gene targets for an insecticidal mode of action with low off-target toxicity. Due to the global participation in the sequencing and annotation of insect genomes, the targeting of specific genes with molecular tools like RNAi and CRISPR/Cas9 for genome engineering and consequent functional identification and validation has become more efficient. While there are multiple avenues to explore for insecticidal candidates, this review identifies G-protein coupled receptors as attractive targets, and hones in on the octopamine and dopamine receptors due to their potential.

Transcriptome Analysis of the Tadpole Shrimp (Triops longicaudatus) by Illumina Paired-End Sequencing: Assembly, Annotation, and Marker Discovery

The tadpole shrimp (Triops longicaudatus) is an aquatic crustacean that helps control pest populations. It inhabits freshwater ponds and pools and has been described as a living fossil. T. longicaudatus was officially declared an endangered species South Korea in 2005; however, through subsequent protection and conservation management, it was removed from the endangered species list in 2012. The limited number of available genetic resources on T. longicaudatus makes it difficult to obtain valuable genetic information for marker-aided selection programs. In this study, whole-transcriptome sequencing of T. longicaudatus generated 39.74 GB of clean data and a total of 269,822 contigs using the Illumina HiSeq 2500 platform. After clustering, a total of 208,813 unigenes with an N₅₀ length of 1089 bp were generated. A total of 95,105 unigenes were successfully annotated against Protostome (PANM), Unigene, Eukaryotic Orthologous Groups (KOG), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases using BLASTX with a cut-off of 1E−5. A total of 57,731 unigenes were assigned to GO terms, and 7247 unigenes were mapped to 129 KEGG pathways. Furthermore, 1595 simple sequence repeats (SSRs) were detected from the unigenes with 1387 potential SSR markers. This is the first report of high-throughput transcriptome analysis of T. longicaudatus, and it provides valuable insights for genetic research and molecular-assisted breeding of this important species.

Beyond the survival and death of the deltamethrin-threatened pollen beetle Meligethes aeneus: An in-depth proteomic study employing a transcriptome database

Insecticide resistance is an increasingly global problem that hampers pest control. We sought the mechanism responsible for survival following pyrethroid treatment and the factors connected to paralysis/death of the pollen beetle Meligethes aeneus through a proteome-level analysis using nanoLC coupled with Orbitrap Fusion™ Tribrid™ mass spectrometry. A tolerant field population of beetles was treated with deltamethrin, and the ensuing proteome changes were observed in the survivors (resistant), dead (paralyzed) and control-treated beetles. The protein database consisted of the translated transcriptome, and the resulting changes were manually annotated via BLASTP. We identified a number of high-abundance changes in which there were several dominant proteins, e.g., the electron carrier cytochrome b5, ribosomal proteins 60S RPL28, 40S RPS23 and RPS26, eIF4E-transporter, anoxia up-regulated protein, 2 isoforms of vitellogenin and pathogenesis-related protein 5. Deltamethrin detoxification was influenced by different cytochromes P450, which were likely boosted by increased cytochrome b5, but glutathione-S-transferase ε and UDP-glucuronosyltransferases also contributed. Moreover, we observed changes in proteins related to RNA interference, RNA binding and epigenetic modifications. The high changes in ribosomal proteins and associated factors suggest specific control of translation. Overall, we showed modulation of expression processes by epigenetic markers, alternative splicing and translation. Future functional studies will benefit.

BIOLOGICAL SIGNIFICANCE

Insects develop pesticide resistance, which has become one of the key issues in plant protection. This growing resistance increases the demand for pesticide applications and the development of new substances. Knowledge in the field regarding the resistance mechanism and its responses to pesticide treatment provides us the opportunity to propose a solution for this issue. Although the pollen beetle Meligethes aeneus was effectively controlled with pyrethroids for many years, there have been reports of increasing resistance. We show protein changes including production of isoforms in response to deltamethrin at the protein level. These results illustrate the insect's survival state as a resistant beetle and in its paralyzed state (evaluated as dead) relative to resistant individuals.

Transcriptomic Analysis of Eucryptorrhynchus chinensis (Coleoptera: Curculionidae) Using 454 Pyrosequencing Technology

Eucryptorrhynchus chinensis Olivier (Coleoptera: Curculionidae) is one of the most important pests of Ailanthus altissima; however, so far, no studies on the genome or transcriptome of E. chinensis have been reported. Using the Roche 454 FLX Titanium platform, an RNA pool obtained from E. chinensis eggs, larva, pupae, and adults was sequenced and assembled de novo to achieve maximum diversity of sampled transcripts. We obtained 1,441,137 (∼518 Mb) raw reads with an average length of 360 bp. After trimming, 89% qualified reads were produced and assembled into 35,509 isotigs with an average length of 440 bp, N50 of 1,048 bp, and 111,643 singletons. We generated 87,894 unigenes following a cluster analysis of the isotigs and singletons, and then functionally annotated the unigenes with gene descriptions. We obtained 23,363 GO assignments, and 12,724 unigenes were assigned to KOG. Based on these annotations, 294 biochemical pathways involved in growth, reproduction, and stress or immune responses were predicted. A total of 659,026 single nucleotide variants and 6,112 simple sequence repeats were detected. Our data provide comprehensive information on the sequence and possible functions of E. chinensis transcripts.

Molecular taxonomic analysis of the plant associations of adult pollen beetles (Nitidulidae: Meligethinae), and the population structure of Brassicogethes aeneus

Pollen beetles (Nitidulidae: Meligethinae) are among the most abundant flower-visiting insects in Europe. While some species damage millions of hectares of crops annually, the biology of many species is little known. We assessed the utility of a 797 base pair fragment of the cytochrome oxidase 1 gene to resolve molecular operational taxonomic units (MOTUs) in 750 adult pollen beetles sampled from flowers of 63 plant species sampled across the UK and continental Europe. We used the same locus to analyse region-scale patterns in population structure and demography in an economically important pest, Brassicogethes aeneus. We identified 44 Meligethinae at ∼2% divergence, 35 of which contained published sequences. A few specimens could not be identified because the MOTUs containing them included published sequences for multiple Linnaean species, suggesting either retention of ancestral haplotype polymorphism or identification errors in published sequences. Over 90% of UK specimens were identifiable as B. aeneus. Plant associations of adult B. aeneus were found to be far wider taxonomically than for their larvae. UK B. aeneus populations showed contrasting affiliations between the north (most similar to Scandinavia and the Baltic) and south (most similar to western continental Europe), with strong signatures of population growth in the south.

Transcriptome and Difference Analysis of Fenpropathrin Resistant Predatory Mite, Neoseiulus barkeri (Hughes)

Several fenpropathrin-resistant predatory mites have been reported. However, the molecular mechanism of the resistance remains unknown. In the present study, the Neoseiulus barkeri (N. barkeri) transcriptome was generated using the Illumina sequencing platform, 34,211 unigenes were obtained, and 15,987 were manually annotated. After manual annotation, attentions were attracted to resistance-related genes, such as voltage-gated sodium channel (VGSC), cytochrome P450s (P450s), and glutathione S-transferases (GSTs). A polymorphism analysis detected two point mutations (E1233G and S1282G) in the linker region between VGSC domain II and III. In addition, 43 putative P450 genes and 10 putative GST genes were identified from the transcriptome. Among them, two P450 genes, NbCYP4EV2 and NbCYP4EZ1, and four GST genes, NbGSTd01, NbGSTd02, NbGSTd03 and NbGSTm03, were remarkably overexpressed 3.64-46.69-fold in the fenpropathrin resistant strain compared to that in the susceptible strain. These results suggest that fenpropathrin resistance in N. barkeri is a complex biological process involving many genetic changes and provide new insight into the N. barkeri resistance mechanism.

Insecticide resistance in pollen beetles over 7 years - a landscape approach

BACKGROUND

In spite of considerable interest in the impact of pesticides on pest populations, few attempts have been made to link resistance patterns of insect pests to land-use features across spatial and temporal scales. We hypothesise that pollen beetle pesticide resistance increases in areas with a high proportion of oilseed rape and with an even mixture of winter and spring oilseed rape owing to high pesticide selection pressure in such areas.

RESULTS

Here, we investigated 7 years of lambda-cyhalothrin (Karate(®) ) resistance in field-collected pollen beetle adults from a total of 180 sampling points across ten regions in Sweden. We found a positive effect on pollen beetle pesticide resistance of proportion of oilseed rape and even spring-winter oilseed rape mixture. However, this was true only for the regional spatial scale. Significant land-use effects in the long-term models, with oilseed rape data averaged over a longer (4 years) period of time, suggested an effect of regional landscape history on current pest resistance.

CONCLUSION

For successful control of pollen beetle pesticide resistance, we suggest a long-term regional strategy for oilseed rape management. This land-use approach provides a framework for further investigations that integrate resistance management into landscape research.

Insight into the Meligethes aeneus voltage-sensitive sodium channel structure and an attempt to select the best pyrethroid ligands

BACKGROUND

The voltage-sensitive sodium channel (VSSC) is a target for the pharmacological action of pyrethroids which are used in controlling pests, including those of agricultural importance. Among these is the pollen beetle (Meligethes aeneus F.) - the most serious pest of Brassica napus. Owing to the heavy use of pyrethroids, a widespread build-up of resistance has occurred. The main cause of pyrethroid insensitivity in M. aeneus is considered to be an increased oxidative metabolism; however, the additional mechanism of resistance associated with mutations in the VSSC might contribute to this phenomenon.

RESULTS

We generated a VSSC 3D model to study the docking affinities of pyrethroids to their target site within the channel. Our goal was to identify the pyrethroids for which docking affinity scores were high and not affected by potential mutations in the VSSC. We found that the docking scores of cypermethrin are hardly influenced by the appearance of point mutations. Additionally, tau-fluvalinate, deltamethrin and bifenthrin are VSSC ligands with high affinity scores.

CONCLUSIONS

Our docking models suggest that point mutations in the VSSC binding pocket might affect the stability of ligand interactions and change the pattern of ligand docking locations, which might have a potential effect on VSSC gating properties.

De novo assembly of the Carcinus maenas transcriptome and characterization of innate immune system pathways

Background

The European shore crab, Carcinus maenas, is used widely in biomonitoring, ecotoxicology and for studies into host-pathogen interactions. It is also an important invasive species in numerous global locations. However, the genomic resources for this organism are still sparse, limiting research progress in these fields. To address this resource shortfall we produced a C. maenas transcriptome, enabled by the progress in next-generation sequencing technologies, and applied this to assemble information on the innate immune system in this species.

Results

We isolated and pooled RNA for twelve different tissues and organs from C. maenas individuals and sequenced the RNA using next generation sequencing on an Illumina HiSeq 2500 platform. After de novo assembly a transcriptome was generated encompassing 212,427 transcripts (153,699 loci). The transcripts were filtered, annotated and characterised using a variety of tools (including BLAST, MEGAN and RSEM) and databases (including NCBI, Gene Ontology and KEGG). There were differential patterns of expression for between 1,223 and 2,741 transcripts across tissues and organs with over-represented Gene Ontology terms relating to their specific function. Based on sequence homology to immune system components in other organisms, we show both the presence of transcripts for a series of known pathogen recognition receptors and response proteins that form part of the innate immune system, and transcripts representing the RNAi, Toll-like receptor signalling, IMD and JAK/STAT pathways.

Conclusions

We have produced an assembled transcriptome for C. maenas that provides a significant molecular resource for wide ranging studies in this species. Analysis of the transcriptome has revealed the presence of a series of known targets and functional pathways that form part of their innate immune system and illustrate tissue specific differences in their expression patterns.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1667-1) contains supplementary material, which is available to authorized users.