Sequencing, De Novo assembly and annotation of the Colorado Potato Beetle, Leptinotarsa decemlineata, Transcriptome

Identification of chitinase genes and roles in the larval-pupal transition of Leptinotarsa decemlineata

BACKGROUND

Insect chitinases play crucial roles in degrading chitin in the extracellular matrix, affecting insect development and molting. However, our understanding of the specific functions of various chitinases in Leptinotarsa decemlineata is limited, hindering the deployment of novel gene-targeting technologies as pest management strategies.

RESULTS

We identified and characterized 19 full-length complementary DNA (cDNA) sequences of chitinase genes (LdChts) in Leptinotarsa decemlineata. Despite having varying domain architectures, all these chitinases contained at least one chitinase catalytic domain. Phylogenetic analysis classified the chitinase proteins into ten distinct clusters (groups I-X). Expression profiles showed the highest expression in chitin-rich tissues or during specific developmental stages from the larva-to-pupa transition. Gene-specific RNA interference (RNAi) experiments provided valuable insight into chitinase gene function. Silencing of group II LdCht10 prevented larval-larval molting, larval-prepupal, and prepupal-pupal processes. Moreover, our study revealed that LdCht5, LdCht2, LdCht11, LdCht1, and LdCht3 from groups I and VII-X were specifically essential for the transition from prepupal to pupal stage, whereas LdIDGF2 from group V was necessary for the larval-prepupal metamorphic process. The chitinase gene LdCht7 from group III and LdIDGF4 from group V were involved in both the larva-to-prepupa and the prepupa-to-pupa shift. Additionally, our findings also shed light on the exclusive expression of nine chitinase genes within group IV in the digestive system, suggesting their potential role in regulating larval body weight and larva-to-pupa transition.

CONCLUSION

Our results provide a comprehensive understanding of the functional specialization of chitinase genes during the molting process of various stages and identify potential targets for RNAi-based management of Leptinotarsa decemlineata. © 2023 Society of Chemical Industry.

The effects of short-term glyphosate-based herbicide exposure on insect gene expression profiles

Glyphosate-based herbicides (GBHs) are the most frequently used herbicides worldwide. The use of GBHs is intended to tackle weeds, but GBHs have been shown to affect the life-history traits and antioxidant defense system of invertebrates found in agroecosystems. Thus far, the effects of GBHs on detoxification pathways among invertebrates have not been sufficiently investigated. We performed two different experiments-1) the direct pure glyphosate and GBH treatment, and 2) the indirect GBH experiment via food-to examine the possible effects of environmentally relevant GBH levels on the survival of the Colorado potato beetle (Leptinotarsa decemlineata) and the expression profiles of their detoxification genes. As candidate genes, we selected four cytochrome P450 (CYP), three glutathione-S-transferase (GST), and two acetylcholinesterase (AChE) genes that are known to be related to metabolic or target-site resistances in insects. We showed that environmentally relevant levels of pure glyphosate and GBH increased the probability for higher mortality in the Colorado potato beetle larvae in the direct experiment, but not in the indirect experiment. The GBHs or glyphosate did not affect the expression profiles of the studied CYP, GST, or AChE genes; however, we found a large family-level variation in expression profiles in both the direct and indirect treatment experiments. These results suggest that the genes selected for this study may not be the ones expressed in response to glyphosate or GBHs. It is also possible that the relatively short exposure time did not affect gene expression profiles, or the response may have already occurred at a shorter exposure time. Our results show that glyphosate products may affect the survival of the herbivorous insect already at lower levels, depending on their sensitivity to pesticides.

The Metagenomic Analysis of Viral Diversity in Colorado Potato Beetle Public NGS Data

The Colorado potato beetle (CPB) is one of the most serious insect pests due to its high ecological plasticity and ability to rapidly develop resistance to insecticides. The use of biological insecticides based on viruses is a promising approach to control insect pests, but the information on viruses which infect leaf feeding beetles is scarce. We performed a metagenomic analysis of 297 CPB genomic and transcriptomic samples from the public National Center for Biotechnology Information Sequence Read Archive (NCBI SRA) database. The reads that were not aligned to the reference genome were assembled with metaSPAdes, and 13314 selected contigs were analyzed with BLAST tools. The contigs and non-aligned reads were also analyzed with Kraken2 software. A total of 3137 virus-positive contigs were attributed to different viruses belonging to 6 types, 17 orders, and 32 families, matching over 97 viral species. The annotated sequences can be divided into several groups: those that are homologous to genetic sequences of insect viruses (Adintoviridae, Ascoviridae, Baculoviridae, Dicistroviridae, Chuviridae, Hytrosaviridae, Iflaviridae, Iridoviridae, Nimaviridae, Nudiviridae, Phasmaviridae, Picornaviridae, Polydnaviriformidae, Xinmoviridae etc.), plant viruses (Betaflexiviridae, Bromoviridae, Kitaviridae, Potyviridae), and endogenous retroviral elements (Retroviridae, Metaviridae). Additionally, the full-length genomes and near-full length genome sequences of several viruses were assembled. We also found sequences belonging to Bracoviriform viruses and, for the first time, experimentally validated the presence of bracoviral genetic fragments in the CPB genome. Our work represents the first attempt to discover the viral genetic material in CPB samples, and we hope that further studies will help to identify new viruses to extend the arsenal of biopesticides against CPB.

The N-glycosylation-related genes as potential targets for RNAi-mediated pest control of the Colorado potato beetle (Leptinotarsa decemlineata)

BACKGROUND

N-glycosylation is one of the most common and important post-translational modifications in the eukaryotic cell. The study of protein N-glycosylation in several model insects confirmed the importance of this process in insect development, immunity, survival and fertility. The Colorado potato beetle (Leptinotarsa decemlineata) (CPB) is a common pest of Solanaceae crops. With the infamous title of champion of insecticide resistance, novel pest control strategies for this insect are needed. Luckily this pest insect is reported as very sensitive for the post-genomic technology of RNA interference (RNAi).

RESULTS

In this project, we investigated the importance of N-glycosylation in the survival and development of CPB using RNAi-mediated gene silencing of N-glycosylation-related genes (NGRGs) during the different transition steps from the larva, through the pupa to the adult stage. High mortality was observed in the larval stage with the silencing of early NGRGs, as STT3a, DAD1 and GCS1. With dsRNA against middle NGRGs, abnormal phenotypes at the ecdysis process and adult formation were observed, while the silencing of late NGRGs did not cause mortality.

CONCLUSION

The lethal phenotypes observed on silencing of the genes involved in the early processing steps of the N-glycosylation pathway suggest these genes are good candidates for RNAi-mediated control of CPB. Next to the gene-specific mechanism of RNAi for biosafety and possible implementation in integrated pest management, we believe these early NGRGs provide a possible alternative to the well-known target genes Snf7 and vacuolar ATPases that are now used in the first commercial RNAi-based products and thus they may be useful in the context of proactive resistance management. © 2021 Society of Chemical Industry.

Sequence variation and regulatory variation in acetylcholinesterase genes contribute to insecticide resistance in different populations of Leptinotarsa decemlineata

Although insect herbivores are known to evolve resistance to insecticides through multiple genetic mechanisms, resistance in individual species has been assumed to follow the same mechanism. While both mutations in the target site insensitivity and increased amplification are known to contribute to insecticide resistance, little is known about the degree to which geographic populations of the same species differ at the target site in a response to insecticides. We tested structural (e.g., mutation profiles) and regulatory (e.g., the gene expression of Ldace1 and Ldace2, AChE activity) differences between two populations (Vermont, USA and Belchow, Poland) of the Colorado potato beetle, Leptinotarsa decemlineata in their resistance to two commonly used groups of insecticides, organophosphates, and carbamates. We established that Vermont beetles were more resistant to azinphos-methyl and carbaryl insecticides than Belchow beetles, despite a similar frequency of resistance-associated alleles (i.e., S291G) in the Ldace2 gene. However, the Vermont population had two additional amino acid replacements (G192S and F402Y) in the Ldace1 gene, which were absent in the Belchow population. Moreover, the Vermont population showed higher expression of Ldace1 and was less sensitive to AChE inhibition by azinphos-methyl oxon than the Belchow population. Therefore, the two populations have evolved different genetic mechanisms to adapt to organophosphate and carbamate insecticides.

Effect of food source availability in the salivary gland transcriptome of the unique burying beetle Nicrophorus pustulatus (Coleoptera: Silphidae)

Nicrophorus is a genus of beetles that bury and transform small vertebrate carcasses into a brood ball coated with their oral and anal secretions to prevent decay and that will serve as a food source for their young. Nicrophorus pustulatus is an unusual species with the ability to overtake brood of other burying beetles and whose secretions, unlike other Nicrophorus species, has been reported not to exhibit antimicrobial properties. This work aims to better understand how the presence or absence of a food source influences the expression of genes involved in the feeding process of N. pustulatus. To achieve that, total RNA was extracted from pooled samples of salivary gland tissue from N. pustulatus and sequenced using an Illumina platform. The resulting reads were used to assemble a de novo transcriptome using Trinity. Duplicates with more than 95% similarity were removed to obtain a "unigene" set. Annotation of the unigene set was done using the Trinotate pipeline. Transcript abundance was determined using Kallisto and differential gene expression analysis was performed using edgeR. A total of 651 genes were found to be differentially expressed, including 390 upregulated and 261 downregulated genes in fed insects compared to starved. Several genes upregulated in fed beetles are associated with the insect immune response and detoxification processes with only one transcript encoding for the antimicrobial peptide (AMP) defensin. These results confirm that N. pustulatus does not upregulate the production of genes encoding AMPs during feeding. This study provides a snapshot of the changes in gene expression in the salivary glands of N. pustulatus following feeding while providing a well described transcriptome for the further analysis of this unique burying beetle.

Transcriptome analysis unravels RNAi pathways genes and putative expansion of CYP450 gene family in cotton leafhopper Amrasca biguttula (Ishida)

Cotton Leafhopper, Amrasca biguttula is an important pest of cotton and okra in the Indian subcontinent. Presently limited genomic/transcriptomic information is available for this insect in any of open source databases. The present study reports the first assembled and annotated de novo transcriptome of cotton leafhopper. Out of 75,551 transcripts, 39,613 CDS (Coding Sequence) were predicted with 35,282 showing positive blast hits with NCBI nr database. The Gene ontology (GO) analysis annotated 7431 CDS  with KEGG pathway categorizing these CDS into 22 different functional groups. The majority of CDS were annotated in signal transduction and transport catabolism pathways. The sequence data was screened for RNAi pathway genes and presence of 37 transcripts associated with this process confirmed the existence of robust RNAi machinery. The role of core RNAi machinery genes (Dicer-2, Ago-2, Piwi and Staufen) has been validated through dsRNA feeding studies. The data resource has also been used to identify potential RNAi targets and genes associated with insecticide detoxification specifically CYP 450 family. The current study provides a useful sequence resource which can be used to initiate molecular studies in this insect with emphasis on insecticide resistance, RNAi and functional genomics.

Transcriptomics-Based Approach Identifies Spinosad-Associated Targets in the Colorado Potato Beetle, Leptinotarsa decemlineata

Simple Summary

The Colorado potato beetle Leptinotarsa decemlineata is a potato pest that can cause substantial damages to potato crops worldwide. Multiple approaches have been leveraged to control this pest including the use of a variety of insecticides. Resistance to different insecticides aimed at controlling this insect has been reported and much work has been conducted in recent years to elucidate the underlying molecular changes associated with insecticide resistance in L. decemlineata. However, information is sparse regarding the molecular impact associated with spinosad treatment in this insect pest. The current study thus explores transcriptional changes associated with spinosad response in L. decemlineata exposed to this compound using high-throughput sequencing. Results presented show multiple transcripts of interest that exhibit differential expression in spinosad-treated L. decemlineata and provide a preliminary footprint of transcripts affected by this insecticide in this potato pest. Select targets identified in this signature should be further explored in follow-up studies to better characterize their contribution, if any, in the process of spinosad resistance.

Abstract

The Colorado potato beetle Leptinotarsa decemlineata is an insect pest that threatens potato crops globally. The primary method to control its damage on potato plants is the use of insecticides, including imidacloprid, chlorantraniliprole and spinosad. However, insecticide resistance has been frequently observed in Colorado potato beetles. The molecular targets and the basis of resistance to imidacloprid and chlorantraniliprole have both been previously quantified. This work was undertaken with the overarching goal of better characterizing the molecular changes associated with spinosad exposure in this insect pest. Next-generation sequencing was conducted to identify transcripts that were differentially expressed between Colorado potato beetles exposed to spinosad versus control insects. Results showed several transcripts that exhibit different expression levels between the two conditions, including ones coding for venom carboxylesterase-6, chitinase 10, juvenile hormone esterase and multidrug resistance-associated protein 4. In addition, several microRNAs, such as miR-12-3p and miR-750-3p, were also modulated in the investigated conditions. Overall, this work reveals a molecular footprint underlying spinosad response in Colorado potato beetles and provides novel leads that could be targeted as part of RNAi-based approaches to control this insect pest.

Gene Sequences of Potential Targets of Insecticidal PF2 Lectin Identified from the Larval De Novo Transcriptome of the Mexican Bean Weevil (Zabrotes Subfasciatus; Boheman 1833)

Simple Summary

The Mexican bean weevil Zabrotes subfasciatus is a major insect pest of stored beans. We have previously reported that the PF2 lectin, which is a protein found in the desert wild legume Olneya tesota (Palo Fierro), is toxic to Z. subfasciatus by inhibiting its early larval development. The use of proteomic means allowed us to identify PF2 targets in the midgut of Z. subfasciatus larvae. However, efforts to completely elucidate the insecticidal mechanism of PF2, as well as novel potential targets for insecticidal compounds, have been hindered by the lack of available genomic and proteomic information of non-model insects. Therefore, in this work we massively sequenced and analyzed the transcripts expressed in the larval stage of Z. subfasciatus, which is the first transcriptome reported for this insect. A total of 29,029 transcript sequences were identified, of which 30 sequences encode putative targets of PF2. The functional characteristics and biochemical, biological, or molecular roles for 15,124 sequences were established by means of bioinformatics tools. This study significantly increased the available genetic resources for Zabrotes and related insect species and will be helpful for any kind of future study that requires information on genes or protein sequences.

Abstract

The available genomic and proteomic information of non-model organisms is often underrepresented in public databases hindering their study at molecular, cellular, and physiological levels. Information on Zabrotes subfasciatus (Mexican bean weevil) is poorly represented in databases, yet it is a major pest of common beans. We report the transcriptome of Z. subfasciatus larvae; transcripts were sequenced using an Illumina RNA-Seq technology and assembled de novo identifying 29,029 unigenes with an average size of 1168 bp and an N50 value of 2196 bp. About 15,124 unigenes (52%) were functionally annotated and categorized. Further analysis revealed 30 unigene sequences encoding putative targets of the insecticidal PF2 lectin. The complete deduced amino acid sequences of eight selected proteins potentially related to insecticidal mechanism of Palo Fierro 2 (PF2) were used for predicting probable N-glycosylation sites and analyzing phylogenetic relationships with insect sequences. This work provides a dramatic increase in the genetic resources available for Coleopterans and set the basis for developing future studies on biological aspects and potential control strategies for Z. subfasciatus.

Modern Techniques in Colorado Potato Beetle (Leptinotarsa decemlineata Say) Control and Resistance Management: History Review and Future Perspectives

Simple Summary

The Colorado potato beetle (CPB) is one of the most important potato pest worldwide. It is native to U.S. but during the 20th century it has dispersed through Europe, Asia and western China. It continues to expand in an east and southeast direction. Damages are caused by larvae and adults. Their feeding on potato plant leaves can cause complete defoliation and lead to a large yield loss. After the long period of using only chemical control measures, the emergence of resistance increased and some new and different methods come to the fore. The main focus of this review is on new approaches to the old CPB control problem. We describe the use of Bacillus thuringiensis and RNA interference (RNAi) as possible solutions for the future in CPB management. RNAi has proven successful in controlling many pests and shows great potential for CPB control. Better understanding of the mechanisms that affect efficiency will enable the development of this technology and boost potential of RNAi to become part of integrated plant protection in the future. We described also the possibility of using single nucleotide polymorphisms (SNPs) as a way to go deeper into our understanding of resistance and how it influences genotypes.

Abstract

Colorado potato beetle, CPB (Leptinotarsa decemlineata Say), is one of the most important pests of the potato globally. Larvae and adults can cause complete defoliation of potato plant leaves and can lead to a large yield loss. The insect has been successfully suppressed by insecticides; however, over time, has developed resistance to insecticides from various chemical groups, and its once successful control has diminished. The number of available active chemical control substances is decreasing with the process of testing, and registering new products on the market are time-consuming and expensive, with the possibility of resistance ever present. All of these concerns have led to the search for new methods to control CPB and efficient tools to assist with the detection of resistant variants and monitoring of resistant populations. Current strategies that may aid in slowing resistance include gene silencing by RNA interference (RNAi). RNAi, besides providing an efficient tool for gene functional studies, represents a safe, efficient, and eco-friendly strategy for CPB control. Genetically modified (GM) crops that produce the toxins of Bacillus thuringiensis (Bt) have many advantages over agro-technical, mechanical, biological, and chemical measures. However, pest resistance that may occur and public acceptance of GM modified food crops are the main problems associated with Bt crops. Recent developments in the speed, cost, and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs) and field of population genomics. There is a need for effective resistance monitoring programs that are capable of the early detection of resistance and successful implementation of integrated resistance management (IRM). The main focus of this review is on new technologies for CPB control (RNAi) and tools (SNPs) for detection of resistant CPB populations.

The ABCB Multidrug Resistance Proteins Do Not Contribute to Ivermectin Detoxification in the Colorado Potato Beetle, Leptinotarsa decemlineata (Say)

The Colorado potato beetle, Leptinotarsa decemlineata (Say), is a significant agricultural pest that has developed resistance to many insecticides that are used to control it. Investigating the mechanisms of insecticide detoxification in this pest is important for ensuring its continued control, since they may be contributors to such resistance. Multidrug resistance (MDR) genes that code for the ABCB transmembrane efflux transporters are one potential source of insecticide detoxification activity that have not been thoroughly examined in L. decemlineata. In this study, we annotated the ABCB genes found in the L. decemlineata genome and then characterized the expression profiles across midgut, nerve, and Malpighian tubule tissues of the three full transporters identified. To investigate if these genes are involved in defense against the macrocyclic lactone insecticide ivermectin in this insect, each gene was silenced using RNA interference or MDR protein activity was inhibited using a chemical inhibitor, verapamil, before challenging the insects with a dose of ivermectin. Survival of the insects did not significantly change due to gene silencing or protein inhibition, suggesting that MDR transporters do not significantly contribute to defense against ivermectin in L. decemlineata.

Interactive Responses of Solanum Dulcamara to Drought and Insect Feeding are Herbivore Species-Specific

In nature, plants are frequently subjected to multiple biotic and abiotic stresses, resulting in a convergence of adaptive responses. We hypothesised that hormonal signalling regulating defences to different herbivores may interact with drought responses, causing distinct resistance phenotypes. To test this, we studied the hormonal and transcriptomic responses of Solanum dulcamara subjected to drought and herbivory by the generalist Spodoptera exigua (beet armyworm; BAW) or the specialist Leptinotarsa decemlineata (Colorado potato beetle; CPB). Bioassays showed that the performance of BAW, but not CPB, decreased on plants under drought compared to controls. While drought did not alter BAW-induced hormonal responses, it enhanced the CPB-induced accumulation of jasmonic acid and salicylic acid (SA), and suppressed ethylene (ET) emission. Microarray analyses showed that under drought, BAW herbivory enhanced several herbivore-induced responses, including cell-wall remodelling and the metabolism of carbohydrates, lipids, and secondary metabolites. In contrast, CPB herbivory enhanced several photosynthesis-related and pathogen responses in drought-stressed plants. This may divert resources away from defence production and increase leaf nutritive value. In conclusion, while BAW suffers from the drought-enhanced defences, CPB may benefit from the effects of enhanced SA and reduced ET signalling. This suggests that the fine-tuned interaction between the plant and its specialist herbivore is sustained under drought.

Coordinated RNA-Seq and peptidomics identify neuropeptides and G-protein coupled receptors (GPCRs) in the large pine weevil Hylobius abietis, a major forestry pest

Hylobius abietis (Linnaeus), or large pine weevil (Coleoptera, Curculionidae), is a pest of European coniferous forests. In order to gain understanding of the functional physiology of this species, we have assembled a de novo transcriptome of H. abietis, from sequence data obtained by Next Generation Sequencing. In particular, we have identified genes encoding neuropeptides, peptide hormones and their putative G-protein coupled receptors (GPCRs) to gain insights into neuropeptide-modulated processes. The transcriptome was assembled de novo from pooled paired-end, sequence reads obtained from RNA from whole adults, gut and central nervous system tissue samples. Data analysis was performed on the transcripts obtained from the assembly including, annotation, gene ontology and functional assignment as well as transcriptome completeness assessment and KEGG pathway analysis. Pipelines were created using Bioinformatics tools and techniques for prediction and identification of neuropeptides and neuropeptide receptors. Peptidomic analysis was also carried out using a combination of MALDI-TOF as well as Q-Exactive Orbitrap mass spectrometry to confirm the identified neuropeptide. 41 putative neuropeptide families were identified in H. abietis, including Adipokinetic hormone (AKH), CAPA and DH31. Neuropeptide F, which has not been yet identified in the model beetle T. castaneum, was identified. Additionally, 24 putative neuropeptide and 9 leucine-rich repeat containing G protein coupled receptor-encoding transcripts were determined using both alignment as well as non-alignment methods. This information, submitted to the NCBI sequence read archive repository (SRA accession: SRP133355), can now be used to inform understanding of neuropeptide-modulated physiology and behaviour in H. abietis; and to develop specific neuropeptide-based tools for H. abietis control.

Cap 'n' collar C regulates genes responsible for imidacloprid resistance in the Colorado potato beetle, Leptinotarsa decemlineata

The Colorado potato beetle (CPB), Leptinotarsa decemlineata developed resistance to imidacloprid after exposure to this insecticide for multiple generations. Our previous studies showed that xenobiotic transcription factor, cap 'n' collar isoform C (CncC) regulates the expression of multiple cytochrome P450 genes, which play essential roles in resistance to plant allelochemicals and insecticides. In this study, we sought to obtain a comprehensive picture of the genes regulated by CncC in imidacloprid-resistant CPB. We performed sequencing of RNA isolated from imidacloprid-resistant CPB treated with dsRNA targeting CncC or gene coding for green fluorescent protein (control). Comparative transcriptome analysis showed that CncC regulated the expression of 1798 genes, out of which 1499 genes were downregulated in CncC knockdown beetles. Interestingly, expression of 79% of imidacloprid induced P450 genes requires CncC. We performed quantitative real-time PCR to verify the reduction in the expression of 20 genes including those coding for detoxification enzymes (P450s, glutathione S-transferases, and esterases) and ABC transporters. The genes coding for ABC transporters are induced in insecticide resistant CPB and require CncC for their expression. Knockdown of genes coding for ABC transporters simultaneously or individually caused an increase in imidacloprid-induced mortality in resistant beetles confirming their contribution to insecticide resistance. These studies identified CncC as a transcription factor involved in regulation of genes responsible for imidacloprid resistance. Small molecule inhibitors of CncC or suppression of CncC by RNAi could provide effective synergists for pest control or management of insecticide resistance.

Revealing the complex genetic structure of cultivated amaryllis (Hippeastrum hybridum) using transcriptome-derived microsatellite markers

Although amaryllis (Hippeastrum hybridum) plants are commonly used in physiological and ecological research, the extent of their genomic and genetic resources remains limited. The development of molecular markers is therefore of great importance to accelerate genetic improvements in Hippeastrum species. In this study, a total of 269 unique genes were defined that might regulate the flower spathe development of amaryllis. In addition, 2000 simple sequence repeats (SSRs) were detected based on 171,462 de novo assembled unigenes from transcriptome data, and 66,4091 single nucleotide polymorphisms (SNPs) were also detected as putative molecular markers. Twenty-one SSR markers were screened to evaluate the genetic diversity and population structure of 104 amaryllis accessions. A total of 98 SSR loci were amplified for all accessions. The results reveal that Nei’s gene diversity (H) values of these markers ranged between 0.055 and 0.394, whereas the average values of Shannon’s Information index (I) ranged between 0.172 and 0.567. Genetic tree analysis further demonstrates that all accessions can be grouped into three main clusters, which can be further divided into two subgroups. STRUCTURE-based analysis revealed that the highest ΔK values were observed when K = 5, K = 6, K = 7 and K = 8. The results of this study enable large-scale transcriptomics and classification of Hippeastrum genetic polymorphisms and will be useful in the future for resource conservation and production.

Characterization of ladybird Henosepilachna vigintioctopunctata transcriptomes across various life stages

Henosepilachna vigintioctopunctata is a vegetable pest that has spread worldwide. It belongs to the Coccinellidae family, whose members exhibit remarkable diversity, both in terms of their diets and the colored spots that appear on the elytra in the adult stage. Transcriptomic data from H. vigintioctopunctata at different life stages would be useful for further investigating the genetic basis of this dietary diversity and the formation of the colored spots in ladybird beetles, as well as revealing the population dynamics of H. vigintioctopunctata, which could be useful in pest control. Here, we generated a comprehensive RNA-seq data set (a total of ~24 Gb of clean data) for H. vigintioctopunctata by sequencing samples collected at different life stages. We characterized the transcriptomes of each of the four life stages (egg, larva, pupa, adult) and generated a high-coverage pool by combining all the RNA-seq reads. Furthermore, we identified a catalog of simple sequence repeat (SSR) markers. This represents the first study to collect transcriptome data from all life stages of a ladybird beetle.

A model species for agricultural pest genomics: the genome of the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae)

The Colorado potato beetle is one of the most challenging agricultural pests to manage. It has shown a spectacular ability to adapt to a variety of solanaceaeous plants and variable climates during its global invasion, and, notably, to rapidly evolve insecticide resistance. To examine evidence of rapid evolutionary change, and to understand the genetic basis of herbivory and insecticide resistance, we tested for structural and functional genomic changes relative to other arthropod species using genome sequencing, transcriptomics, and community annotation. Two factors that might facilitate rapid evolutionary change include transposable elements, which comprise at least 17% of the genome and are rapidly evolving compared to other Coleoptera, and high levels of nucleotide diversity in rapidly growing pest populations. Adaptations to plant feeding are evident in gene expansions and differential expression of digestive enzymes in gut tissues, as well as expansions of gustatory receptors for bitter tasting. Surprisingly, the suite of genes involved in insecticide resistance is similar to other beetles. Finally, duplications in the RNAi pathway might explain why Leptinotarsa decemlineata has high sensitivity to dsRNA. The L. decemlineata genome provides opportunities to investigate a broad range of phenotypes and to develop sustainable methods to control this widely successful pest.

Ancestry & molecular evolutionary analyses of heat shock protein 47 kDa (HSP47/SERPINH1)

HSP47/SERPINH1 is key-regulator for collagen biosynthesis and its structural assembly. To date, there is no comprehensive study on the phylogenetic history of HSP47. Herein we illustrate the evolutionary history of HSP47/SERPINH1 along with sequence, structural and syntenic traits for HSP47/SERPINH1. We have identified ancestral HSP47/SERPINH1 locus in Japanese lamprey (Lethenteron japonicum). This gene remains on the same or similar locus for ~500 million years (MY), but chromosomal duplication was observed in ray-finned fishes, leading into three sets of three sets (I-III) of HSP47/SERPINH1. Two novel introns were inserted at the positions 36b and 102b in the first exon of only HSP47_1 gene from the selected ray-finned fishes. On the evolutionary time scale, the events of HSP47 duplications took placed between 416–360 MY ago (MYA) while intron insertion dates back to 231–190 MYA after early divergence of ray-finned fishes.

Overexpression of a cytochrome P450 and a UDP-glycosyltransferase is associated with imidacloprid resistance in the Colorado potato beetle, Leptinotarsa decemlineata

Current control of insect pests relies on chemical insecticides, however, insecticide resistance development by pests is a growing concern in pest management. The main mechanisms for insecticide resistance typically involve elevated activity of detoxifying enzymes and xenobiotic transporters that break-down and excrete insecticide molecules. In this study, we investigated the molecular mechanisms of imidacloprid resistance in the Colorado potato beetle, Leptinotarsa decemlineata (Say) (Coleoptera: Chrysomelidae), an insect pest notorious for its capacity to develop insecticide resistance rapidly. We compared the transcriptome profiles of imidacloprid-resistant and sensitive beetle strains and identified 102 differentially expressed transcripts encoding detoxifying enzymes and xenobiotic transporters. Of these, 74 were up-regulated and 28 were down-regulated in the resistant strain. We then used RNA interference to knock down the transcript levels of seven up-regulated genes in the resistant beetles. Ingestion of double-stranded RNA successfully knocked down the expression of the genes for three cytochrome P450s (CYP6BQ15, CYP4Q3 and CYP4Q7), one ATP binding cassette (ABC) transporter (ABC-G), one esterase (EST1), and two UDP-glycosyltransferases (UGT1 and UGT2). Further, we demonstrated that silencing of CYP4Q3 and UGT2 significantly increased susceptibility of resistant beetles to imidacloprid, indicating that overexpression of these two genes contributes to imidacloprid resistance in this resistant strain.

Temporal patterns of imidacloprid resistance throughout a growing season in Leptinotarsa decemlineata populations

BACKGROUND

The Colorado potato beetle, Leptinotarsa decemlineata (Say), is a major agricultural pest of commercial potatoes. Pest managers use a combination of control tactics to limit populations, including multiple insecticides. Finding a window of insecticide susceptibility and understanding genetic responses to insecticide exposure during a growing season may provide novel management recommendations for L. decemlineata.

RESULTS

We examined temporal changes (during one growing season) in phenotypic response between a susceptible population and an imidacloprid-resistant population. Beetles remained more susceptible to imidacloprid in the susceptible population throughout the growing season. Estimated mean LC₅₀ values varied throughout the growing season in the resistant population, with increased susceptibility among overwintered and recently emerged adult beetles compared with a heightened level of resistance in the second generation. RNA transcript abundance was compared among multiple time points through the growing season, showing that cuticular proteins and cytochrome p450s were highly upregulated during peaks of measured resistance.

CONCLUSION

Temporal variation in imidacloprid susceptibility of L. decemlineata was observed, which included early time points of susceptibility and later peaks in resistance. Heightened resistance occurred during the second generation and correlated to increased transcript abundance of multiple mechanisms of resistance, including multiple cuticular protein and cytochrome p450 transcripts. © 2016 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

RNA interference of three up-regulated transcripts associated with insecticide resistance in an imidacloprid resistant population of Leptinotarsa decemlineata

The Colorado potato beetle, Leptinotarsa decemlineata (Say), is a major agricultural pest of potatoes in the Central Sands production region of Wisconsin. Previous studies have shown that populations of L. decemlineata have become resistant to many classes of insecticides, including the neonicotinoid insecticide, imidacloprid. Furthermore, L. decemlineata has multiple mechanisms of resistance to deal with a pesticide insult, including enhanced metabolic detoxification by cytochrome p450s and glutathione S-transferases. With recent advances in the transcriptomic analysis of imidacloprid susceptible and resistant L. decemlineata populations, it is possible to investigate the role of candidate genes involved in imidacloprid resistance. A recently annotated transcriptome analysis of L. decemlineata was obtained from select populations of L. decemlineata collected in the Central Sands potato production region, which revealed a subset of mRNA transcripts constitutively up-regulated in resistant populations. We hypothesize that a portion of the up-regulated transcripts encoding for genes within the resistant populations also encode for pesticide resistance and can be suppressed to re-establish a susceptible phenotype. In this study, a discrete set of three up-regulated targets were selected for RNA interference experiments using a resistant L. decemlineata population. Following the successful suppression of transcripts encoding for a cytochrome p450, a cuticular protein, and a glutathione synthetase protein in a select L. decemlineata population, we observed reductions in measured resistance to imidacloprid that strongly suggest these genes control essential steps in imidacloprid metabolism in these field populations.

Deep, Staged Transcriptomic Resources for the Novel Coleopteran Models Atrachya menetriesi and Callosobruchus maculatus

Despite recent efforts to sample broadly across metazoan and insect diversity, current sequence resources in the Coleoptera do not adequately describe the diversity of the clade. Here we present deep, staged transcriptomic data for two coleopteran species, Atrachya menetriesi (Faldermann 1835) and Callosobruchus maculatus (Fabricius 1775). Our sampling covered key stages in ovary and early embryonic development in each species. We utilized this data to build combined assemblies for each species which were then analysed in detail. The combined A. menetriesi assembly consists of 228,096 contigs with an N50 of 1,598 bp, while the combined C. maculatus assembly consists of 128,837 contigs with an N50 of 2,263 bp. For these assemblies, 34.6% and 32.4% of contigs were identified using Blast2GO, and 97% and 98.3% of the BUSCO set of metazoan orthologs were present, respectively. We also carried out manual annotation of developmental signalling pathways and found that nearly all expected genes were present in each transcriptome. Our analyses show that both transcriptomes are of high quality. Lastly, we performed read mapping utilising our timed, stage specific RNA samples to identify differentially expressed contigs. The resources presented here will provide a firm basis for a variety of experimentation, both in developmental biology and in comparative genomic studies.

Digestive proteolysis in the Colorado potato beetle, Leptinotarsa decemlineata: Activity-based profiling and imaging of a multipeptidase network

The Colorado potato beetle (CPB), Leptinotarsa decemlineata, is a major pest of potato plants, and its digestive system is a promising target for development of pest control strategies. This work focuses on functional proteomic analysis of the digestive proteolytic enzymes expressed in the CPB gut. We identified a set of peptidases using imaging with specific activity-based probes and activity profiling with selective substrates and inhibitors. The secreted luminal peptidases were classified as: (i) endopeptidases of cathepsin D, cathepsin L, and trypsin types and (ii) exopeptidases with aminopeptidase (cathepsin H), carboxypeptidase (serine carboxypeptidase, prolyl carboxypeptidase), and carboxydipeptidase (cathepsin B) activities. The proteolytic arsenal also includes non-luminal peptidases with prolyl oligopeptidase and metalloaminopeptidase activities. Our results indicate that the CPB gut employs a multienzyme network of peptidases with complementary specificities to efficiently degrade ingested proteins. This proteolytic system functions in both CPB larvae and adults and is controlled mainly by cysteine and aspartic peptidases and supported by serine and metallopeptidases. The component enzymes identified here are potential targets for inhibitors with tailored specificities that could be engineered into potato plants to confer resistance to CPB.

Sequencing and de novo assembly of the Asian gypsy moth transcriptome using the Illumina platform

The Asian gypsy moth (Lymantria dispar) is a serious pest of forest and shade trees in many Asian and some European countries. However, there have been few studies of L. dispar genetic information and comprehensive genetic analyses of this species are needed in order to understand its genetic and metabolic sensitivities, such as the molting mechanism during larval development. In this study, high-throughput sequencing technology was used to sequence the transcriptome of the Asian subspecies of the gyspy moth, after which a comprehensive analysis of chitin metabolism was undertaken. We generated 37,750,380 high-quality reads and assembled them into contigs. A total of 37,098 unigenes were identified, of which 15,901 were annotated in the NCBI non-redundant protein database and 9,613 were annotated in the Swiss-Prot database. We mapped 4,329 unigenes onto 317 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. Chitin metabolism unigenes were found in the transcriptome and the data indicated that a variety of enzymes was involved in chitin catabolic and biosynthetic pathways.

Extraordinary Adaptive Plasticity of Colorado Potato Beetle: "Ten-Striped Spearman" in the Era of Biotechnological Warfare

Expanding from remote areas of Mexico to a worldwide scale, the ten-striped insect, the Colorado potato beetle (CPB, Leptinotarsa decemlineata Say), has risen from being an innocuous beetle to a prominent global pest. A diverse life cycle, phenotypic plasticity, adaptation to adverse conditions, and capability to detoxify or tolerate toxins make this insect appear to be virtually “indestructible”. With increasing advances in molecular biology, tools of biotechnological warfare were deployed to combat CPB. In the last three decades, genetically modified potato has created a new challenge for the beetle. After reviewing hundreds of scientific papers dealing with CPB control, it became clear that even biotechnological means of control, if used alone, would not defeat the Colorado potato beetle. This control measure once again appears to be provoking the potato beetle to exhibit its remarkable adaptability. Nonetheless, the potential for adaptation to these techniques has increased our knowledge of this pest and thus opened possibilities for devising more sustainable CPB management programs.

The De Novo Transcriptome and Its Functional Annotation in the Seed Beetle Callosobruchus maculatus

Despite their unparalleled biodiversity, the genomic resources available for beetles (Coleoptera) remain relatively scarce. We present an integrative and high quality annotated transcriptome of the beetle Callosobruchus maculatus, an important and cosmopolitan agricultural pest as well as an emerging model species in ecology and evolutionary biology. Using Illumina sequencing technology, we sequenced 492 million read pairs generated from 51 samples of different developmental stages (larvae, pupae and adults) of C. maculatus. Reads were de novo assembled using the Trinity software, into a single combined assembly as well as into three separate assemblies based on data from the different developmental stages. The combined assembly generated 218,192 transcripts and 145,883 putative genes. Putative genes were annotated with the Blast2GO software and the Trinotate pipeline. In total, 33,216 putative genes were successfully annotated using Blastx against the Nr (non-redundant) database and 13,382 were assigned to 34,100 Gene Ontology (GO) terms. We classified 5,475 putative genes into Clusters of Orthologous Groups (COG) and 116 metabolic pathways maps were predicted based on the annotation. Our analyses suggested that the transcriptional specificity increases with ontogeny. For example, out of 33,216 annotated putative genes, 51 were only expressed in larvae, 63 only in pupae and 171 only in adults. Our study illustrates the importance of including samples from several developmental stages when the aim is to provide an integrative and high quality annotated transcriptome. Our results will represent an invaluable resource for those working with the ecology, evolution and pest control of C. maculatus, as well for comparative studies of the transcriptomics and genomics of beetles more generally.

Physiological and molecular mechanisms underlying photoperiodism in the spider mite: comparisons with insects

Photoperiodism is an adaptive, seasonal timing system that enables organisms to coordinate their development and physiology to annual changes in the environment using day length (photoperiod) as a cue. This review summarizes our knowledge of the physiological mechanisms underlying photoperiodism in spider mites. In particular, the two-spotted spider mite Tetranychus urticae is focussed, which has long been used as a model species for studying photoperiodism. Photoperiodism is established by several physiological modules, such as the photoreceptor, photoperiodic time measurement system, counter system, and endocrine effector. It is now clear that retinal photoreception through the ocelli is indispensable for the function of photoperiodism, at least in T. urticae. Visual pigment, which comprised opsin protein and a vitamin A-based pigment, is involved in photoreception. The physiological basis of the photoperiodic time measurement system is still under debate, and we have controversial evidence for the hourglass-based time measurement and the oscillator-based time measurement. Less attention has been centred on the counter system in insects and mites. Mite reproduction is possibly regulated by the ecdysteroid, ponasterone A. Prior physiological knowledge has laid the foundation for the next steps essential for the elucidation of the molecular mechanisms driving photoperiodism.

Comparison of De Novo Transcriptome Assemblers and k-mer Strategies Using the Killifish, Fundulus heteroclitus

Background

De novo assembly of non-model organism’s transcriptomes has recently been on the rise in concert with the number of de novo transcriptome assembly software programs. There is a knowledge gap as to what assembler software or k-mer strategy is best for construction of an optimal de novo assembly. Additionally, there is a lack of consensus on which evaluation metrics should be used to assess the quality of de novo transcriptome assemblies.

Result

Six different assembly strategies were evaluated from four different assemblers. The Trinity assembly was used in its default 25 single k-mer value while Bridger, Oases, and SOAPdenovo-Trans were performed with multiple k-mer strategies. Bridger, Oases, and SOAPdenovo-Trans used a small multiple k-mer (SMK) strategy consisting of the k-mer lengths of 21, 25, 27, 29, 31, and 33. Additionally, Oases and SOAPdenovo-Trans were performed using a large multiple k-mer (LMK) strategy consisting of k-mer lengths of 25, 35, 45, 55, 65, 75, and 85. Eleven metrics were used to evaluate each assembly strategy including three genome related evaluation metrics (contig number, N50 length, Contigs >1 kb, reads) and eight transcriptome evaluation metrics (mapped back to transcripts (RMBT), number of full length transcripts, number of open reading frames, Detonate RSEM-EVAL score, and percent alignment to the southern platyfish, Amazon molly, BUSCO and CEGMA databases). The assembly strategy that performed the best, that is it was within the top three of each evaluation metric, was the Bridger assembly (10 of 11) followed by the Oases SMK assembly (8 of 11), the Oases LMK assembly (6 of 11), the Trinity assembly (4 of 11), the SOAP LMK assembly (4 of 11), and the SOAP SMK assembly (3 of 11).

Conclusion

This study provides an in-depth multi k-mer strategy investigation concluding that the assembler itself had a greater impact than k-mer size regardless of the strategy employed. Additionally, the comprehensive performance transcriptome evaluation metrics utilized in this study identified the need for choosing metrics centered on user defined research goals. Based on the evaluation metrics performed, the Bridger assembly was able to construct the best assembly of the testis transcriptome in Fundulus heteroclitus.

A specialist herbivore pest adaptation to xenobiotics through up-regulation of multiple Cytochrome P450s

The adaptation of herbivorous insects to their host plants is hypothesized to be intimately associated with their ubiquitous development of resistance to synthetic pesticides. However, not much is known about the mechanisms underlying the relationship between detoxification of plant toxins and synthetic pesticides. To address this knowledge gap, we used specialist pest Colorado potato beetle (CPB) and its host plant, potato, as a model system. Next-generation sequencing (454 pyrosequencing) was performed to reveal the CPB transcriptome. Differential expression patterns of cytochrome P450 complement (CYPome) were analyzed between the susceptible (S) and imidacloprid resistant (R) beetles. We also evaluated the global transcriptome repertoire of CPB CYPome in response to the challenge by potato leaf allelochemicals and imidacloprid. The results showed that more than half (51.2%) of the CBP cytochrome P450 monooxygenases (P450s) that are up-regulated in the R strain are also induced by both host plant toxins and pesticide in a tissue-specific manner. These data suggest that xenobiotic adaptation in this specialist herbivore is through up-regulation of multiple P450s that are potentially involved in detoxifying both pesticide and plant allelochemicals.

Are feeding preferences and insecticide resistance associated with the size of detoxifying enzyme families in insect herbivores?

The size of gene families associated with xenobiotic detoxification in insects may be associated with the complexity of their diets and their propensities to develop insecticide resistance. We test these hypotheses by collating the annotations of cytochrome P450, carboxyl/cholinesterase and glutathione S-transferase genes in 65 insect species with data on their host use and history of insecticide resistance. We find 2-4 fold variation across the species in the numbers of these genes and, in some orders, especially the Hymenoptera, there is a clear relationship between the numbers of genes and feeding preferences. However in other orders, in particular the Lepidoptera, no such relationship is apparent. The size of these three gene families also tend to correlate with insecticide resistance propensity but this may not be an independent effect because species with broader host ranges are more likely to be pests that are heavily sprayed with insecticides.

Characterizing Molecular Mechanisms of Imidacloprid Resistance in Select Populations of Leptinotarsa decemlineata in the Central Sands Region of Wisconsin

The Colorado potato beetle, Leptinotarsa decemlineata (Say), is a major agricultural pest in the Central Sands region of Wisconsin. Imidacloprid, a neonicotinoid insecticide, has commonly been used for control of L. decemlineata since its registration in 1995. In the last 10 years, many field populations of L. decemlineata have begun to show increasing imidacloprid resistance. We studied resistance phenotype as a phenomenon that reduces neonicotinoid efficacy and has practical consequences for potato pest management. Although we have not observed complete field failure following the use of these products, multiple studies have demonstrated that the lethal concentration to kill 50% of the test organisms (LC50) in different field populations of L. decemlineata varies greatly which may suggest that resistance of L. decemlineata is heritable and involves genetic changes. An important challenge in understanding resistance is assessing the genetic mechanisms associated with resistance and classifying up-regulated genes that may be involved in combating an insecticide insult. In this study we uncovered trends in imidacloprid phenotypic response that have developed in the region by estimating the LC50 values among different field populations against a range of imidacloprid doses. The LC50 values collected in 2008-2011, and more recently in 2013 and 2014, show that some field locations remain susceptible to imidacloprid, while nearby fields (<100km) have developed high levels of resistance. We also sought to uncover potential mechanisms of resistance at each field location. We compiled a transcriptome for populations, characterized as phenotypically 'susceptible' and 'resistant', by isolating mRNA from adult beetles and analyzing gene expression level differences. Strong differences were observed in constituently up and down-regulated genes among different field populations. Most significantly, the up-regulation of 3 cytochrome p450s and a glutathione synthetase related protein in multiple resistant populations provide a mechanistic explanation of resistance evolution in L. decemlineata.

CHARACTERIZATION AND FUNCTIONAL STUDY OF A PUTATIVE JUVENILE HORMONE DIOL KINASE IN THE COLORADO POTATO BEETLE Leptinotarsa decemlineata (Say)

Juvenile hormone diol kinase (JHDK) is an enzyme involved in JH degradation. In the present article, a putative JHDK cDNA (LdJHDK) was cloned from the Colorado potato beetle Leptinotarsa decemlineata. The cDNA consists of 814 bp, containing a 555 bp open reading frame encoding a 184 amino acid protein. LdJHDK reveals a high degree of identity to the previously reported insect JHDKs. It possesses three conserved purine nucleotide-binding elements, and contains three EF-hand motifs (helix-loop-helix structural domains). LdJHDK mRNA was mainly detected in hindgut and Malpighian tubules. Besides, a trace amount of LdJHDK mRNA was also found in thoracic muscles, brain-corpora cardiaca-corpora allata complex, foregut, midgut, ventral ganglia, fat body, epidermis, and hemocytes. Moreover, LdJHDK was expressed throughout all developmental stages. Within the first, second, and third larval instar, the expression levels of LdJHDK were higher just before and right after the molt, and were lower in the intermediate instar. In the fourth larval instar, the highest peak of LdJHDK occurred 56 h after ecdysis. Ingestion of double-stranded RNA (dsRNA) against LdJHDK successfully knocked down the target gene, increased JH titer, and significantly upregulated LdKr-h1 mRNA level. Knockdown of LdJHDK significantly impaired adult emergence. Thus, we provide a line of experimental evidence in L. decemlineata to support that LdJHDK encodes function protein involved in JH degradation.

De Novo Assembly and Genome Analyses of the Marine-Derived Scopulariopsis brevicaulis Strain LF580 Unravels Life-Style Traits and Anticancerous Scopularide Biosynthetic Gene Cluster

The marine-derived Scopulariopsis brevicaulis strain LF580 produces scopularides A and B, which have anticancerous properties. We carried out genome sequencing using three next-generation DNA sequencing methods. De novo hybrid assembly yielded 621 scaffolds with a total size of 32.2 Mb and 16298 putative gene models. We identified a large non-ribosomal peptide synthetase gene (nrps1) and supporting pks2 gene in the same biosynthetic gene cluster. This cluster and the genes within the cluster are functionally active as confirmed by RNA-Seq. Characterization of carbohydrate-active enzymes and major facilitator superfamily (MFS)-type transporters lead to postulate S. brevicaulis originated from a soil fungus, which came into contact with the marine sponge Tethya aurantium. This marine sponge seems to provide shelter to this fungus and micro-environment suitable for its survival in the ocean. This study also builds the platform for further investigations of the role of life-style and secondary metabolites from S. brevicaulis.

Knocking down a putative Δ(1) -pyrroline-5-carboxylate dehydrogenase gene by RNA interference inhibits flight and causes adult lethality in the Colorado potato beetle Leptinotarsa decemlineata (Say)

BACKGROUND

Leptinotarsa decemlineata is an able disperser by flight. Novel control strategies must be explored to control the damage and inhibit the dispersal efficiently. Proline is a major energy substrate during flight. Δ-Pyrroline-5-carboxylate dehydrogenase (P5CDh) catalyses the second step of proline degradation for the production of ATP.

RESULTS

A full-length Ldp5cdh cDNA was cloned. Ldp5cdh was ubiquitously expressed in the eggs, the first through fourth larval instars, wandering larvae, pupae and adults. In the adults, Ldp5cdh mRNA was widely distributed in thorax muscles, midgut, foregut, hindgut, Malpighian tubules, ventral ganglion, fat body and epidermis, with the expression levels from the highest to the lowest. Two double-stranded RNAs (dsRNAs) (dsLdp5cdh1 and dsLdp5cdh2) targeting Ldp5cdh were constructed and bacterially expressed. Ingestion of dsLdp5cdh1 and dsLdp5cdh2 successfully silenced Ldp5cdh, significantly increased the contents of proline, arginine and alanine, but strongly decreased the contents of asparate, asparagine, glutamate and glutamine in the haemolymph. Moreover, knocking down Ldp5cdh significantly reduced ATP content, decreased flight speed, shortened flight distance and increased adult mortality.

CONCLUSIONS

It seems that identified Ldp5cdh encodes a functional P5CDh enzyme, and Ldp5cdh may serve as a potential target for dsRNA-based pesticide for controlling the damage and dispersal of L. decemlineata adults. © 2014 Society of Chemical Industry.

Instar-dependent systemic RNA interference response in Leptinotarsa decemlineata larvae

RNA interference (RNAi) is a promising approach to control Leptinotarsa decemlineata. In this study, RNAi efficiency by double-stranded RNA (dsRNA) targeting S-adenosyl-L-homocysteine hydrolase (LdSAHase) was compared among L. decemlineata first- to fourth-instar larvae. Ingesting dsLdSAHase successfully decreased the target gene expression, caused lethality, inhibited growth and impaired pupation in an instar- and concentration-dependent manner. To study the role of Dicer2 and Argonaute2 genes in RNAi efficiency, we identified LdDcr2a, LdDcr2b, LdAgo2a and LdAgo2b. Their expression levels were higher in young larvae than those in old ones. Exposure to dsegfp for 6 h significantly elevated LdDcr2a, LdDcr2b, LdAgo2a and LdAgo2b mRNA levels in the first-, second-, third- and fourth-instar larvae. When the exposure periods were extended, however, the expression levels were gradually reduced. Continuous exposure for 72 h significantly repressed the expression of LdAgo2a and LdAgo2b in the first, second and third larval instars, and the four genes in final instars. Moreover, we found that dsLdSAHase-caused LdSAHase suppressions and larval mortalities were influenced by previous dsegfp exposure: 12 h of previous exposure increased LdSAHase silencing and mortality of the final instar larvae, whereas 72 h of exposure reduced LdSAHase silencing and mortality. Thus, it seems the activities of core RNAi-machinery proteins affect RNAi efficiency in L. decemlineata.

Functions of nuclear receptor HR3 during larval-pupal molting in Leptinotarsa decemlineata (Say) revealed by in vivo RNA interference

Our previous results revealed that RNA interference-aided knockdown of Leptinotarsa decemlineata FTZ-F1 (LdFTZ-F1) reduced 20E titer, and impaired pupation. In this study, we characterized a putative LdHR3 gene, an early-late 20E-response gene upstream of LdFTZ-F1. Within the first, second and third larval instars, three expression peaks of LdHR3 occurred just before the molt. In the fourth (final) larval instar 80 h after ecdysis and prepupal stage 3 days after burying into soil, two LdHR3 peaks occurred. The LdHR3 expression peaks coincide with the peaks of circulating 20E level. In vitro midgut culture and in vivo bioassay revealed that 20E and an ecdysteroid agonist halofenozide (Hal) enhanced LdHR3 expression in the final larval instars. Conversely, a decrease in 20E by feeding a double-stranded RNA (dsRNA) against an ecdysteroidogenesis gene Ldshd repressed the expression. Moreover, Hal rescued the transcript levels in the Ldshd-silenced larvae. Thus, 20E peaks activate the expression of LdHR3. Furthermore, ingesting dsRNA against LdHR3 successfully knocked down the target gene, and impaired pupation. Finally, knockdown of LdHR3 upregulated the transcription of three ecdysteroidogenesis genes (Ldphm, Lddib and Ldshd), increased 20E titer, and activated the expression of two 20E-response genes (LdEcR and LdFTZ-F1). Thus, LdHR3 functions in regulation of pupation in the Colorado potato beetle.

Identification of carboxylesterase genes and their expression profiles in the Colorado potato beetle Leptinotarsa decemlineata treated with fipronil and cyhalothrin

Based on the Leptinotarsa decemlineata transcriptome dataset and the GenBank sequences, 70 novel carboxylesterases and 2 acetylcholinesterases were found. The 72 members belong to a multifunctional carboxylesterase/cholinesterase superfamily (CCE). A phylogenetic tree including the 72 LdCCEs and the CCEs from Tribolium castaneum, Drosophila melanogaster and Apis mellifera revealed that all CCEs fell into three main phylogenetic groups: dietary/detoxification, hormone/semiochemical processing, and neurodevelopmental classes. Numbers of L. decemlineata CCEs in the three classes were 52, 12 and 8, respectively. The dietary/detoxification class includes two clades: coleopteran xenobiotic metabolizing and α-esterase type CCEs. CCEs in the two clades have independently expanded in L. decemlineata. The hormone/semiochemical processing class has three clades: integument CCEs, β- and pheromone CCEs and juvenile hormone CCEs. Integument CCEs in L. decemlineata have also expanded. The neurodevelopmental CCEs are implicated the most ancient class, containing acetylcholinesterase, neuroligin, neurotactin, glutactin, gliotactin and others. Among the 70 novel CCE genes, KM220566, KM220530, KM220576, KM220527 and KM220541 were fipronil-inducible, and KM220578, KM220566, KM220542, KM220564, KM220561, KM220554, KM220527, KM220538 and KM220541 were cyhalothrin-inducible. They were the candidates involving in insecticide detoxification. Moreover, our results also provided a platform to understand the functions and evolution of L. decemlineata CCE genes.

Bayesian phylogeny analysis of vertebrate serpins illustrates evolutionary conservation of the intron and indels based six groups classification system from lampreys for ∼500 MY

The serpin superfamily is characterized by proteins that fold into a conserved tertiary structure and exploits a sophisticated and irreversible suicide-mechanism of inhibition. Vertebrate serpins are classified into six groups (V1-V6), based on three independent biological features-genomic organization, diagnostic amino acid sites and rare indels. However, this classification system was based on the limited number of mammalian genomes available. In this study, several non-mammalian genomes are used to validate this classification system using the powerful Bayesian phylogenetic method. This method supports the intron and indel based vertebrate classification and proves that serpins have been maintained from lampreys to humans for about 500 MY. Lampreys have fewer than 10 serpins, which expand into 36 serpins in humans. The two expanding groups V1 and V2 have SERPINB1/SERPINB6 and SERPINA8/SERPIND1 as the ancestral serpins, respectively. Large clusters of serpins are formed by local duplications of these serpins in tetrapod genomes. Interestingly, the ancestral HCII/SERPIND1 locus (nested within PIK4CA) possesses group V4 serpin (A2APL1, homolog of α 2-AP/SERPINF2) of lampreys; hence, pointing to the fact that group V4 might have originated from group V2. Additionally in this study, details of the phylogenetic history and genomic characteristics of vertebrate serpins are revisited.

A de novo transcriptome and valid reference genes for quantitative real-time PCR in Colaphellus bowringi

Background

The cabbage beetle Colaphellus bowringi Baly is a serious insect pest of crucifers and undergoes reproductive diapause in soil. An understanding of the molecular mechanisms of diapause regulation, insecticide resistance, and other physiological processes is helpful for developing new management strategies for this beetle. However, the lack of genomic information and valid reference genes limits knowledge on the molecular bases of these physiological processes in this species.

Results

Using Illumina sequencing, we obtained more than 57 million sequence reads derived from C. bowringi, which were assembled into 39,390 unique sequences. A Clusters of Orthologous Groups classification was obtained for 9,048 of these sequences, covering 25 categories, and 16,951 were assigned to 255 Kyoto Encyclopedia of Genes and Genomes pathways. Eleven candidate reference gene sequences from the transcriptome were then identified through reverse transcriptase polymerase chain reaction. Among these candidate genes, EF1α, ACT1, and RPL19 proved to be the most stable reference genes for different reverse transcriptase quantitative polymerase chain reaction experiments in C. bowringi. Conversely, aTUB and GAPDH were the least stable reference genes.

Conclusion

The abundant putative C. bowringi transcript sequences reported enrich the genomic resources of this beetle. Importantly, the larger number of gene sequences and valid reference genes provide a valuable platform for future gene expression studies, especially with regard to exploring the molecular mechanisms of different physiological processes in this species.

A putative Δ1-pyrroline-5-carboxylate synthetase involved in the biosynthesis of proline and arginine in Leptinotarsa decemlineata

Delta 1-pyrroline-5-carboxylate synthetase (P5CS) catalyzes the conversion of glutamate (Glu) to Glu semialdehyde (GSA). GSA spontaneously cyclizes to form P5C. P5C is then reduced to proline (Pro) or is converted to ornithine, the intermediate for arginine (Arg) biosynthesis. In the present study, a full-length Ldp5cs complementary DNA was cloned from the Colorado potato beetle Leptinotarsa decemlineata, a notorious insect defoliator of potato in most potato-growing regions of the world. Ldp5cs encodes a 792-amino-acid protein which shares high identity to homologues from other insect species. Quantitative reverse transcription polymerase chain reaction revealed that Ldp5cs was ubiquitously expressed in the eggs, first to fourth-instar larvae, wandering larvae, pupae and sexually mature adults. In the adults, Ldp5cs mRNA levels were higher in the fat body, foregut, midgut and hindgut, moderate in the ventral ganglion, lower in the thorax muscles, epidermis and Malpighian tubules. Two double-stranded RNAs (dsRNAs) (dsLdp5cs1 and dsLdp5cs2) targeting Ldp5cs were constructed and bacterially expressed. Ingestion during 3 consecutive days of dsLdp5cs1 or dsLdp5cs2 successfully silenced Ldp5cs, significantly reduced the contents of Pro and Arg in the hemolymph, decreased flight speed and shortened flight distance of the resulting adults. Furthermore, knocking down Ldp5cs significantly increased adult mortality. Thus, our results suggest that identified Ldp5cs encodes a functional P5CS enzyme that is involved in the biosynthesis of Pro and Arg in L. decemlineata.

Photoperiodic effects on diapause-associated gene expression trajectories in European Leptinotarsa decemlineata populations

Behavioural and physiological changes during diapause, an important strategy of insects for surviving harsh seasonal conditions, have been intensively studied. The genetic and molecular mechanisms underpinning diapause development are less well known. We took a candidate gene approach to study prediapause gene expression patterns in the Colorado potato beetle (Leptinotarsa decemlineata), an invasive insect that has rapidly spread northwards to high seasonality environments. Newly eclosed beetles originating from southern (Italy) and northern (Russia) Europe were reared under short- [12 h light (L):12 h dark (D)] and long-day (18L:6D) photoperiods for 10 days. This time period includes the sensitive period for the photoperiodic induction and initiation of diapause. Gene expression trajectories of 12 diapause-related genes (regulatory, metabolic and stress-resistance) were analysed from 0-, 5- and 10-day-old beetles. Gene expression differences increased with age, deviating significantly between populations and photoperiods in 10-day-old beetles. The gene expression profiles, particularly those related to energy metabolism and stress-resistance, indicate that beetles originating from Russia also prepare for diapause under the long-day photoperiod and show qualitative differences in the diapausing phenotype. Our study shows that population-dependent differences seen in behavioural and physiological traits connected with diapause in L. decemlineata are also evident in the expression trajectories of diapause-related genes.