RNA-sequencing analysis reveals abundant developmental stage-specific and immunity-related genes in the pollen beetle Meligethes aeneus

De novo transcriptome assembly and differential gene expression analysis in different developmental stages of Agriotes sputator (click beetle)

Wireworms, the larva of click beetle (Agriotes species), are one of the most destructive pests of horticultural crops in North America, responsible for considerable economic losses in Canada. Agriotes sputator (A. sputator) species is a predominant wireworm pest attacking potato fields in Eastern Canada. However, no information about its genome-wide gene expression profile, specifically for the genes involved with development is available to date. Therefore, we generated the transcriptome profile of A. sputator during five developmental stages, including the three larval stages and adult male and female click beetle. Out of 714.7 million raw reads, de novo assembly generated 564,561 transcripts. The data were subjected to differential expression analysis using DESeq2, gene ontology, annotation, and pathway analyses. A total of 34,709 differentially expressed genes (DEGs) were significant (log2 fold change > 2, padj < 0.05) across the developmental stages. Functional analysis of DEGs identified development signaling, metabolism, transport, cellular mechanisms, and drug metabolism (cytochrome p450) pathways. This study provides comprehensive sequence resources and potential gene differences at different developmental stages of A. sputator. These findings will represent a major step towards developing sustainable methods to control this widely distributed pest in agricultural fields.

Divergent venom effectors correlate with ecological niche in neuropteran predators

Neuropteran larvae are fierce predators that use venom to attack and feed on arthropod prey. Neuropterans have adapted to diverse and sometimes extreme habitats, suggesting their venom may have evolved accordingly, but the ecology and evolution of venom deployment in different families is poorly understood. We applied spatial transcriptomics, proteomics, morphological analysis, and bioassays to investigate the venom systems in the antlion Euroleon nostras and the lacewing Chrysoperla carnea, which occupy distinct niches. Although the venom system morphology was similar in both species, we observed remarkable differences at the molecular level. E. nostras produces particularly complex venom secreted from three different glands, indicating functional compartmentalization. Furthermore, E. nostras venom and digestive tissues were devoid of bacteria, strongly suggesting that all venom proteins are of insect origin rather than the products of bacterial symbionts. We identified several toxins exclusive to E. nostras venom, including phospholipase A2 and several undescribed proteins with no homologs in the C. carnea genome. The compositional differences have significant ecological implications because only antlion venom conferred insecticidal activity, indicating its use for the immobilization of large prey. Our results indicate that molecular venom evolution plays a role in the adaptation of antlions to their unique ecological niche.

Cereal weevils' antimicrobial peptides: at the crosstalk between development, endosymbiosis and immune response

Interactions between animals and microbes are ubiquitous in nature and strongly impact animal physiology. These interactions are shaped by the host immune system, which responds to infections and contributes to tailor the associations with beneficial microorganisms. In many insects, beneficial symbiotic associations not only include gut commensals, but also intracellular bacteria, or endosymbionts. Endosymbionts are housed within specialized host cells, the bacteriocytes, and are transmitted vertically across host generations. Host-endosymbiont co-evolution shapes the endosymbiont genome and host immune system, which not only fights against microbial intruders, but also ensures the preservation of endosymbionts and the control of their load and location. The cereal weevil Sitophilus spp. is a remarkable model in which to study the evolutionary adaptation of the immune system to endosymbiosis owing to its binary association with a unique, relatively recently acquired nutritional endosymbiont, Sodalis pierantonius. This Gram-negative bacterium has not experienced the genome size shrinkage observed in long-term endosymbioses and has retained immunogenicity. We focus here on the sixteen antimicrobial peptides (AMPs) identified in the Sitophilus oryzae genome and their expression patterns in different tissues, along host development or upon immune challenges, to address their potential functions in the defensive response and endosymbiosis homeostasis along the insect life cycle. This article is part of the theme issue 'Sculpting the microbiome: how host factors determine and respond to microbial colonization'.

An Overview of Ovarian Calyx Fluid Proteins of Toxoneuron nigriceps (Viereck) (Hymenoptera: Braconidae): An Integrated Transcriptomic and Proteomic Approach

The larval stages of the tobacco budworm, Heliothis virescens (Fabricius) (Lepidoptera: Noctuidae), are parasitized by the endophagous parasitoid wasp, Toxoneuron nigriceps (Viereck) (Hymenoptera: Braconidae). During the injections of eggs, this parasitoid wasp also injects into the host body the secretion of the venom gland and the calyx fluid, which contains a polydnavirus (T. nigriceps BracoVirus: TnBV) and the Ovarian calyx fluid Proteins (OPs). The effects of the OPs on the host immune system have recently been described. In particular, it has been demonstrated that the OPs cause hemocytes to undergo a number of changes, such as cellular oxidative stress, actin cytoskeleton modifications, vacuolization, and the inhibition of hemocyte encapsulation capacity, which results in both a loss of hemocyte functionality and cell death. In this study, by using a combined transcriptomic and proteomic analysis, the main components of T. nigriceps ovarian calyx fluid proteins were identified and their possible role in the parasitic syndrome was discussed. This study provides useful information to support the analysis of the function of ovarian calyx fluid proteins, to better understand T. nigriceps parasitization success and for a more thorough understanding of the components of ovarian calyx fluid proteins and their potential function in combination with other parasitoid factors.

An Assassin's Secret: Multifunctional Cytotoxic Compounds in the Predation Venom of the Assassin Bug Psytalla horrida (Reduviidae, Hemiptera)

Predatory assassin bugs produce venomous saliva that enables them to overwhelm, kill, and pre-digest large prey animals. Venom from the posterior main gland (PMG) of the African assassin bug Psytalla horrida has strong cytotoxic effects, but the responsible compounds are yet unknown. Using cation-exchange chromatography, we fractionated PMG extracts from P. horrida and screened the fractions for toxicity. Two venom fractions strongly affected insect cell viability, bacterial growth, erythrocyte integrity, and intracellular calcium levels in Drosophila melanogaster olfactory sensory neurons. LC-MS/MS analysis revealed that both fractions contained gelsolin, redulysins, S1 family peptidases, and proteins from the uncharacterized venom protein family 2. Synthetic peptides representing the putative lytic domain of redulysins had strong antimicrobial activity against Escherichia coli and/or Bacillus subtilis but only weak toxicity towards insect or mammalian cells, indicating a primary role in preventing the intake of microbial pathogens. In contrast, a recombinant venom protein family 2 protein significantly reduced insect cell viability but exhibited no antibacterial or hemolytic activity, suggesting that it plays a role in prey overwhelming and killing. The results of our study show that P. horrida secretes multiple cytotoxic compounds targeting different organisms to facilitate predation and antimicrobial defense.

You are what you eat-ecological niche and microhabitat influence venom activity and composition in aquatic bugs

True water bugs (Nepomorpha) are mostly predacious insects that live in aquatic habitats. They use their piercing-sucking mouthparts to inject venomous saliva that facilitates the capture and extra-oral digestion of prey animals, but their venom can also be deployed for defence. In Central Europe, nepomorph species representing different families coexist in the same habitat. However, their feeding ecology, including venom composition and deployment, has not been investigated in detail. We used an integrated proteotranscriptomic and bioactivity-based approach to test whether venom composition and activity differ between four water bug species sharing the same habitat but occupying different ecological niches. We found considerable species-dependent differences in the composition of digestive enzymes and venom components that probably evolved as adaptations to particular food sources, foraging strategies and/or microhabitats. The venom of Corixa punctata differed substantially from that of the three strictly predatory species (Ilyocoris cimicoides, Notonecta glauca and Nepa cinerea), and the abundance of herbivory-associated proteins confirms a mostly plant-based diet. Our findings reveal independent adaptations of the digestive and defensive enzyme repertoires accompanied by the evolution of distinct feeding strategies in aquatic bugs.

Transcriptome analysis of the endangered dung beetle Copris tripartitus (Coleoptera: Scarabaeidae) and characterization of genes associated to immunity, growth, and reproduction

BACKGROUND

Dung beetles recycle organic matter through the decomposition of feces and support ecological balance. However, these insects are threatened by the indiscriminate use of agrochemicals and habitat destruction. Copris tripartitus Waterhouse (Coleoptera: Scarabaeidae), a dung beetle, is listed as a class-II Korean endangered species. Although the genetic diversity of C. tripartitus populations has been investigated through analysis of mitochondrial genes, genomic resources for this species remain limited. In this study, we analyzed the transcriptome of C. tripartitus to elucidate functions related to growth, immunity and reproduction for the purpose of informed conservation planning.

RESULTS

The transcriptome of C. tripartitus was generated using next-generation Illumina sequencing and assembled de novo using a Trinity-based platform. In total, 98.59% of the raw sequence reads were processed as clean reads. These reads were assembled into 151,177 contigs, 101,352 transcripts, and 25,106 unigenes. A total of 23,450 unigenes (93.40%) were annotated to at least one database. The largest proportion of unigenes (92.76%) were annotated to the locally curated PANM-DB. A maximum of 5,512 unigenes had homologous sequences in Tribolium castaneum. Gene Ontology (GO) analysis revealed a maximum of 5,174 unigenes in the Molecular function category. Further, in Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, a total of 462 enzymes were associated with established biological pathways. Based on sequence homology to known proteins in PANM-DB, representative immunity, growth, and reproduction-related genes were screened. Potential immunity-related genes were categorized into pattern recognition receptors (PRRs), the Toll-like receptor signaling pathway, the MyD88- dependent pathway, endogenous ligands, immune effectors, antimicrobial peptides, apoptosis, and adaptation-related transcripts. Among PRRs, we conducted detailed in silico characterization of TLR-2, CTL, and PGRP_SC2-like. Repetitive elements such as long terminal repeats, short interspersed nuclear elements, long interspersed nuclear elements and DNA elements were enriched in the unigene sequences. A total of 1,493 SSRs were identified among all unigenes of C. tripartitus.

CONCLUSIONS

This study provides a comprehensive resource for analysis of the genomic topography of the beetle C. tripartitus. The data presented here clarify the fitness phenotypes of this species in the wild and provide insight to support informed conservation planning.

Cloning and Functional Characterization of a Novel β-GRP Gene From Melanotus cribricollis

In this study, a novel β-1,3-glucan recognition protein gene (β-GRP) was identified from Melanotus cribricollis, and its potential role in the immune response was investigated. The full length of β-GRP cDNA (Accession Number: MT941530) was 1644 bp, encoding a protein composed of 428 amino acids. The theoretical molecular weight and the isoelectric point were 51.53 kDa and 6.17, respectively. The amino acid sequence of β-GRP from M. cribricollis was closely related to that of. β-GRP-like from Photinus pyralis, and was predicted to contain conserved GH16 domain without glucanase active site. The results of real-time quantitative PCR showed that fungal infection of Metarhizium pingshaense may significantly upregulated the expression level of β-GRP gene. The RNAi suppression of β-GRP gene expression significantly increased the corrected cumulative mortality. Meanwhile, antimicrobial peptide genes defensin and lysozyme were significantly downregulated after interference of β-GRP. Taken together, these results suggest that β-GRP of M. cribricollis probably participates in the host immune system by mediating the expression of antimicrobial peptides. This study provides comprehensive insights into the innate immune system of insect larvae.

Storage of carbon reserves in spruce trees is prioritized over growth in the face of carbon limitation

Climate change is expected to pose a global threat to forest health by intensifying extreme events like drought and insect attacks. Carbon allocation is a fundamental process that determines the adaptive responses of long-lived late-maturing organisms like trees to such stresses. However, our mechanistic understanding of how trees coordinate and set allocation priorities among different sinks (e.g., growth and storage) under severe source limitation remains limited. Using flux measurements, isotopic tracing, targeted metabolomics, and transcriptomics, we investigated how limitation of source supply influences sink activity, particularly growth and carbon storage, and their relative regulation in Norway spruce (Picea abies) clones. During photosynthetic deprivation, absolute rates of respiration, growth, and allocation to storage all decline. When trees approach neutral carbon balance, i.e., daytime net carbon gain equals nighttime carbon loss, genes encoding major enzymes of metabolic pathways remain relatively unaffected. However, under negative carbon balance, photosynthesis and growth are down-regulated while sucrose and starch biosynthesis pathways are up-regulated, indicating that trees prioritize carbon allocation to storage over growth. Moreover, trees under negative carbon balance actively increase the turnover rate of starch, lipids, and amino acids, most likely to support respiration and mitigate stress. Our study provides molecular evidence that trees faced with severe photosynthetic limitation strategically regulate storage allocation and consumption at the expense of growth. Understanding such allocation strategies is crucial for predicting how trees may respond to extreme events involving steep declines in photosynthesis, like severe drought, or defoliation by heat waves, late frost, or insect attack.

Identification of a Sulfatase that Detoxifies Glucosinolates in the Phloem-Feeding Insect Bemisia tabaci and Prefers Indolic Glucosinolates

Cruciferous plants in the order Brassicales defend themselves from herbivory using glucosinolates: sulfur-containing pro-toxic metabolites that are activated by hydrolysis to form compounds, such as isothiocyanates, which are toxic to insects and other organisms. Some herbivores are known to circumvent glucosinolate activation with glucosinolate sulfatases (GSSs), enzymes that convert glucosinolates into inactive desulfoglucosinolates. This strategy is a major glucosinolate detoxification pathway in a phloem-feeding insect, the silverleaf whitefly Bemisia tabaci, a serious agricultural pest of cruciferous vegetables. In this study, we identified and characterized an enzyme responsible for glucosinolate desulfation in the globally distributed B. tabaci species MEAM1. In in vitro assays, this sulfatase showed a clear preference for indolic glucosinolates compared with aliphatic glucosinolates, consistent with the greater representation of desulfated indolic glucosinolates in honeydew. B. tabaci might use this detoxification strategy specifically against indolic glucosinolates since plants may preferentially deploy indolic glucosinolates against phloem-feeding insects. In vivo silencing of the expression of the B. tabaci GSS gene via RNA interference led to lower levels of desulfoglucosinolates in honeydew. Our findings expand the knowledge on the biochemistry of glucosinolate detoxification in phloem-feeding insects and suggest how detoxification pathways might facilitate plant colonization in a generalist herbivore.

Genomic analysis of novel Yarrowia-like yeast symbionts associated with the carrion-feeding burying beetle Nicrophorus vespilloides

Background

Mutualistic interactions with microbes can help insects adapt to extreme environments and unusual diets. An intriguing example is the burying beetle Nicrophorus vespilloides, which feeds and reproduces on small vertebrate carcasses. Its fungal microbiome is dominated by yeasts that potentially facilitate carcass utilization by producing digestive enzymes, eliminating cadaver-associated toxic volatiles (that would otherwise attract competitors), and releasing antimicrobials to sanitize the microenvironment. Some of these yeasts are closely related to the biotechnologically important species Yarrowia lipolytica.

Results

To investigate the roles of these Yarrowia-like yeast (YLY) strains in more detail, we selected five strains from two different phylogenetic clades for third-generation sequencing and genome analysis. The first clade, represented by strain B02, has a 20-Mb genome containing ~ 6400 predicted protein-coding genes. The second clade, represented by strain C11, has a 25-Mb genome containing ~ 6300 predicted protein-coding genes, and extensive intraspecific variability within the ITS–D1/D2 rDNA region commonly used for species assignments. Phenotypic microarray analysis revealed that both YLY strains were able to utilize a diverse range of carbon and nitrogen sources (including microbial metabolites associated with putrefaction), and can grow in environments with extreme pH and salt concentrations.

Conclusions

The genomic characterization of five yeast strains isolated from N. vespilloides resulted in the identification of strains potentially representing new YLY species. Given their abundance in the beetle hindgut, and dominant growth on beetle-prepared carcasses, the analysis of these strains has revealed the genetic basis of a potential symbiotic relationship between yeasts and burying beetles that facilitates carcass digestion and preservation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07597-z.

An Integrated-Omics/Chemistry Approach Unravels Enzymatic and Spontaneous Steps to Form Flavoalkaloidal Nudicaulin Pigments in Flowers of Papaver nudicaule L

Flower colour is an important trait for plants to attract pollinators and ensure their reproductive success. Among yellow flower pigments, the nudicaulins in Papaver nudicaule L. (Iceland poppy) are unique due to their rarity and unparalleled flavoalkaloid structure. Nudicaulins are derived from pelargonidin glycoside and indole, products of the flavonoid and indole/tryptophan biosynthetic pathway, respectively. To gain insight into the molecular and chemical basis of nudicaulin biosynthesis, we combined transcriptome, differential gel electrophoresis (DIGE)-based proteome, and ultra-performance liquid chromatography-high resolution mass spectrometry (UPLC-HRMS)-based metabolome data of P. nudicaule petals with chemical investigations. We identified candidate genes and proteins for all biosynthetic steps as well as some key metabolites across five stages of petal development. Candidate genes of amino acid biosynthesis showed a relatively stable expression throughout petal development, whereas most candidate genes of flavonoid biosynthesis showed increasing expression during development followed by downregulation in the final stage. Notably, gene candidates of indole-3-glycerol-phosphate lyase (IGL), sharing characteristic sequence motifs with known plant IGL genes, were co-expressed with flavonoid biosynthesis genes, and are probably providing free indole. The fusion of indole with pelargonidin glycosides was retraced synthetically and promoted by high precursor concentrations, an excess of indole, and a specific glycosylation pattern of pelargonidin. Thus, nudicaulin biosynthesis combines the enzymatic steps of two different pathways with a spontaneous fusion of indole and pelargonidin glycoside under precisely tuned reaction conditions.

Immune-related genes of the larval Holotrichia parallela in response to entomopathogenic nematodes Heterorhabditis beicherriana LF

Background

Entomopathogenic nematodes (EPNs) emerge as compatible alternatives to conventional insecticides in controlling Holotrichia parallela larvae (Coleoptera: Scarabaeidae). However, the immune responses of H. parallela against EPNs infection remain unclear.

Results

In present research, RNA-Seq was firstly performed. A total of 89,427 and 85,741 unigenes were achieved from the midgut of H. parallela larvae treated with Heterorhabditis beicherriana LF for 24 and 72 h, respectively; 2545 and 3156 unigenes were differentially regulated, respectively. Among those differentially expressed genes (DEGs), 74 were identified potentially related to the immune response. Notably, some immune-related genes, such as peptidoglycan recognition protein SC1 (PGRP-SC1), pro-phenoloxidase activating enzyme-I (PPAE-I) and glutathione s-transferase (GST), were induced at both treatment points. Bioinformatics analysis showed that PGRP-SC1, PPAE-I and GST were all involved in anti-parasitic immune process. Quantitative real-time PCR (qRT-PCR) results showed that the three immune-related genes were expressed in all developmental stages; PGRP-SC1 and PPAE-I had higher expressions in midgut and fat body, respectively, while GST exhibited high expression in both of them. Moreover, in vivo silencing of them resulted in increased susceptibility of H. parallela larvae to H. beicherriana LF.

Conclusion

These results suggest that H. parallela PGRP-SC1, PPAE-I and GST are involved in the immune responses to resist H. beicherriana LF infection. This study provides the first comprehensive transcriptome resource of H. parallela exposure to nematode challenge that will help to support further comparative studies on host-EPN interactions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07506-4.

An integrated transcriptomic and proteomic approach to identify the main Torymus sinensis venom components

During oviposition, ectoparasitoid wasps not only inject their eggs but also a complex mixture of proteins and peptides (venom) in order to regulate the host physiology to benefit their progeny. Although several endoparasitoid venom proteins have been identified, little is known about the components of ectoparasitoid venom. To characterize the protein composition of Torymus sinensis Kamijo (Hymenoptera: Torymidae) venom, we used an integrated transcriptomic and proteomic approach and identified 143 venom proteins. Moreover, focusing on venom gland transcriptome, we selected additional 52 transcripts encoding putative venom proteins. As in other parasitoid venoms, hydrolases, including proteases, phosphatases, esterases, and nucleases, constitute the most abundant families in T. sinensis venom, followed by protease inhibitors. These proteins are potentially involved in the complex parasitic syndrome, with different effects on the immune system, physiological processes and development of the host, and contribute to provide nutrients to the parasitoid progeny. Although additional in vivo studies are needed, initial findings offer important information about venom factors and their putative host effects, which are essential to ensure the success of parasitism.

The molecular mechanisms that determine different degrees of polyphagy in the Bemisia tabaci species complex

The whitefly Bemisia tabaci is a closely related group of >35 cryptic species that feed on the phloem sap of a broad range of host plants. Species in the complex differ in their host-range breadth, but the mechanisms involved remain poorly understood. We investigated, therefore, how six different B. tabaci species cope with the environmental unpredictability presented by a set of four common and novel host plants. Behavioral studies indicated large differences in performances on the four hosts and putative specialization of one of the species to cassava plants. Transcriptomic analyses revealed two main insights. First, a large set of genes involved in metabolism (>85%) showed differences in expression between the six species, and each species could be characterized by its own unique expression pattern of metabolic genes. However, within species, these genes were constitutively expressed, with a low level of environmental responsiveness (i.e., to host change). Second, within each species, sets of genes mainly associated with the super-pathways "environmental information processing" and "organismal systems" responded to the host switching events. These included genes encoding for proteins involved in sugar homeostasis, signal transduction, membrane transport, and immune, endocrine, sensory and digestive responses. Our findings suggested that the six B. tabaci species can be divided into four performance/transcriptomic "Types" and that polyphagy can be achieved in multiple ways. However, polyphagy level is determined by the specific identity of the metabolic genes/pathways that are enriched and overexpressed in each species (the species' individual metabolic "tool kit").

Identify immune-related genes of adzuki bean weevil (Callosobruchus chinensis) in response to bacteria challenge by transcriptome analysis

BACKGROUND

Callosobruchus chinensis is one of the important postharvest pests in legume growing areas. Bacterial pesticide is a potential alternative method to control storage pests. However, the effect of these pathogen bacteria on storage pests, and the molecular mechanisms of insect response remain to be to investigated.

RESULTS

Using the next generation sequencing technology, we established a transcriptomic library for C. chinensis larvae in response to Escherichia coli. Total of 355 differential expressed genes (DEGs) were identified, which 178 DEGs were upregulated, and 177 DEGs were downregulated compared to control group. To validate the RNA-seq analysis, 20 DEGs and 14 immune-related genes were selected to perform quantitative polymerase chain reaction (RT-qPCR). These immune-related genes were involved in recognition (peptidoglycan recognition proteins), signal transduction (fibrinogen-related proteins, serine proteinases and NF-κB), and execution effectors (phenoloxidase, defensin, attacin, and antimicrobial peptide). In addition, genes that encode digestive and respiratory enzymes were altered in C. chinensis larvae in response to infection. Some genes that involved in juvenile hormone and insulin pathway appeared to express differentially, suggesting that pathogen infection might lead to developmental arrest. Furthermore, iron homeostasis and chitin metabolism appeared significantly altered after infection.

CONCLUSION

In this study, we characterized the immune response of C. chinensis larvae in response to E. coli using RNA-seq, from pathogen recognition, signal transduction, to execution. Some other identified genes were involved in iron homeostasis, respiration, and digestion. A better understanding of molecular response of beetle to pathogen will facilitate us to develop an available strategy to control storage pests.

Isolation, Identification, and Bioinformatic Analysis of Antibacterial Proteins and Peptides from Immunized Hemolymph of Red Palm Weevil Rhynchophorus ferrugineus

Red palm weevil (Rhynchophorus ferrugineus Olivier, 1791, Coleoptera: Curculionidae) is a destructive pest of palms, rapidly extending its native geographical range and causing large economic losses worldwide. The present work describes isolation, identification, and bioinformatic analysis of antibacterial proteins and peptides from the immunized hemolymph of this beetle. In total, 17 different bactericidal or bacteriostatic compounds were isolated via a series of high-pressure liquid chromatography steps, and their partial amino acid sequences were determined by N-terminal sequencing or by mass spectrometry. The bioinformatic analysis of the results facilitated identification and description of corresponding nucleotide coding sequences for each peptide and protein, based on the recently published R. ferrugineus transcriptome database. The identified compounds are represented by several well-known bactericidal factors: two peptides similar to defensins, one cecropin-A1-like peptide, and one attacin-B-like protein. Interestingly, we have also identified some unexpected compounds comprising five isoforms of pheromone-binding proteins as well as seven isoforms of odorant-binding proteins. The particular role of these factors in insect response to bacterial infection needs further investigation.

Transcriptome and gene expression analysis of Rhynchophorus ferrugineus (Coleoptera: Curculionidae) during developmental stages

Background

Red palm weevil, Rhynchophorus ferrugineus Olivier, is one of the most destructive pests harming palm trees. However, genomic resources for R. ferrugineus are still lacking, limiting the ability to discover molecular and genetic means of pest control.

Methods

In this study, PacBio Iso-Seq and Illumina RNA-seq were used to generate transcriptome from three developmental stages of R. ferrugineus (pupa, 7th-instar larva, adult) to increase the understanding of the life cycle and molecular characteristics of the pest.

Results

Sequencing generated 625,983,256 clean reads, from which 63,801 full-length transcripts were assembled with N50 of 3,547 bp. Expression analyses revealed 8,583 differentially expressed genes (DEGs). Moreover, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that these DEGs were mainly related to the peroxisome pathway which associated with metabolic pathways, material transportation and organ tissue formation. In summary, this work provides a valuable basis for further research on the growth and development, gene expression and gene prediction, and pest control of R. ferrugineus.

A bioinformatic study of antimicrobial peptides identified in the Black Soldier Fly (BSF) Hermetia illucens (Diptera: Stratiomyidae)

Antimicrobial peptides (AMPs) play a key role in the innate immunity, the first line of defense against bacteria, fungi, and viruses. AMPs are small molecules, ranging from 10 to 100 amino acid residues produced by all living organisms. Because of their wide biodiversity, insects are among the richest and most innovative sources for AMPs. In particular, the insect Hermetia illucens (Diptera: Stratiomyidae) shows an extraordinary ability to live in hostile environments, as it feeds on decaying substrates, which are rich in microbial colonies, and is one of the most promising sources for AMPs. The larvae and the combined adult male and female H. illucens transcriptomes were examined, and all the sequences, putatively encoding AMPs, were analysed with different machine learning-algorithms, such as the Support Vector Machine, the Discriminant Analysis, the Artificial Neural Network, and the Random Forest available on the CAMP database, in order to predict their antimicrobial activity. Moreover, the iACP tool, the AVPpred, and the Antifp servers were used to predict the anticancer, the antiviral, and the antifungal activities, respectively. The related physicochemical properties were evaluated with the Antimicrobial Peptide Database Calculator and Predictor. These analyses allowed to identify 57 putatively active peptides suitable for subsequent experimental validation studies.

Context-dependent venom deployment and protein composition in two assassin bugs

The Heteroptera are a diverse suborder of phytophagous, hematophagous, and zoophagous insects. The shift to zoophagy can be traced back to the transformation of salivary glands into venom glands, but the venom is used not only to kill and digest invertebrate prey but also as a defense strategy, mainly against vertebrates. In this study, we used an integrated transcriptomics and proteomics approach to compare the composition of venoms from the anterior main gland (AMG) and posterior main gland (PMG) of the reduviid bugs Platymeris biguttatus L. and Psytalla horrida Stål. In both species, the AMG and PMG secreted distinct protein mixtures with few interspecific differences. PMG venom consisted mostly of S1 proteases, redulysins, Ptu1-like peptides, and uncharacterized proteins, whereas AMG venom contained hemolysins and cystatins. There was a remarkable difference in biological activity between the AMG and PMG venoms, with only PMG venom conferring digestive, neurotoxic, hemolytic, antibacterial, and cytotoxic effects. Proteomic analysis of venom samples revealed the context-dependent use of AMG and PMG venom. Although both species secreted PMG venom alone to overwhelm their prey and facilitate digestion, the deployment of defensive venom was species-dependent. P. biguttatus almost exclusively used PMG venom for defense, whereas P. horrida secreted PMG venom in response to mild harassment but AMG venom in response to more intense harassment. This intriguing context-dependent use of defensive venom indicates that future research should focus on species-dependent differences in venom composition and defense strategies among predatory Heteroptera.

Mechanisms of Novel Host Use by Bactrocera tau (Tephritid: Diptera) Revealed by RNA Transcriptomes

Use of novel plant hosts can facilitate the establishment and range expansion of herbivorous invasive species. However, the inherent mechanisms of novel host use are still unclear in many herbivorous species. Here, we examine mechanisms of novel host use in the invasive tephritid fruit fly Bactrocera tau (Walker)(Diptera: Tephritidae) by documenting changes in the RNA transcriptomes associated with a novel host. RNA transcripts of B. tau were obtained with high-throughput sequencing from samples continuously reared on two traditional Cucurbitaceae hosts and a novel host (banana). We found transcriptome variation was strongly associated with feeding on banana. Moreover, B. tau feeding on banana contained more differentially expressed genes (DEGs) and more annotated categories of DEGs in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database with 1,595 DEGs and 21 major annotated pathways. The annotated categories of DEGs in individuals reared on banana differed with from those individuals feeding on other hosts and were enriched in oxidative phosphorylation, citrate cycle pathway, and four other carbohydrate pathways. For B. tau feeding on banana, the predominant numbers of upregulated genes in the mitochondrial NADH (56 on average) and a relatively higher numbers of upregulated genes (13 on average) were found in oxidative phosphorylation and the TCA pathway, respectively. Changes in RNA transcriptomes associated with novel host use, especially for genes related to energy and carbohydrate metabolism, help to explain how B. tau can be successful in use of novel hosts and may be useful in developing novel strategies for control of tephritid flies.

Developmental and sexual divergence in the olfactory system of the marine insect Clunio marinus

An animal's fitness strongly depends on successful feeding, avoidance of predators and reproduction. All of these behaviours commonly involve chemosensation. As a consequence, when species' ecological niches and life histories differ, their chemosensory abilities need to be adapted accordingly. The intertidal insect Clunio marinus (Diptera: Chironomidae) has tuned its olfactory system to two highly divergent niches. The long-lived larvae forage in a marine environment. During the few hours of terrestrial adult life, males have to find the female pupae floating on the water surface, free the cryptic females from their pupal skin, copulate and carry the females to the oviposition sites. In order to explore the possibility for divergent olfactory adaptations within the same species, we investigated the chemosensory system of C. marinus larvae, adult males and adult females at the morphological and molecular level. The larvae have a well-developed olfactory system, but olfactory gene expression only partially overlaps with that of adults, likely reflecting their marine vs. terrestrial lifestyles. The olfactory system of the short-lived adults is simple, displaying no glomeruli in the antennal lobes. There is strong sexual dimorphism, the female olfactory system being particularly reduced in terms of number of antennal annuli and sensilla, olfactory brain centre size and gene expression. We found hints for a pheromone detection system in males, including large trichoid sensilla and expression of specific olfactory receptors and odorant binding proteins. Taken together, this makes C. marinus an excellent model to study within-species evolution and adaptation of chemosensory systems.

Transcriptomic analysis of the interactions between the Spodoptera exigua midgut and nucleopolyhedrovirus

Spodoptera exigua nucleopolyhedrovirus (SeNPV) has been successfully applied as a bioinsecticide against S. exigua, one of the most devastating pests worldwide. However, due to limited information, the molecular mechanisms underlying interactions between S. exigua and SeNPV remain to be elucidated. In this study, RNA-Seq and differentially expressed gene (DEG) analysis of the S. exigua larva midgut were performed to explore molecular responses to SeNPV infection. A total of 1785 DEGs, including 935 upregulated and 850 downregulated genes, were identified in the midgut of SeNPV-infected S. exigua larvae. Ultrastructural observations showed that after SeNPV infection, the peritrophic matrix (PM) became a loose and highly porous surface with many clear ruptures; these changes were most likely associated with upregulation of chitin deacetylases. In addition, 124 putative innate immunity-related DEGs were identified and divided into several groups, including pattern recognition proteins, signaling pathways, signal modulation, antimicrobial peptides and detoxification. Interestingly, upregulation of some pattern recognition proteins, induction of the JAK/STAT signaling pathway and promotion of REPAT synthesis might be the main innate immunity responses occurring in the S. exigua larva midgut after SeNPV infection. According to quantitative real-time PCR, the expression profiles of 19 random DEGs were consistent with those obtained by RNA-Seq. These findings provide important basic information for understanding the molecular mechanisms of SeNPV invasion and the anti-SeNPV responses of the S. exigua midgut, promoting the utility of SeNPV as a bioinsecticide for the effective control of S. exigua and related pests.

The unique antimicrobial peptide repertoire of stick insects

The comparative analysis of innate immunity across different insect taxa has revealed unanticipated evolutionary plasticity, providing intriguing examples of immunity-related effector gene expansion and loss. Phasmatodea, the stick and leaf insects, is an order of hemimetabolous insects that can provide insight into ancestral innate immunity genes lost by later insect clades. We injected the stick insect Peruphasma schultei with a mixture of microbial elicitors to activate a strong immune response, followed by RNA-Seq analysis to screen for induced immunity-related effector genes. This revealed a highly diverse spectrum of antimicrobial peptides (AMPs) belonging to the attacin, coleoptericin, defensin, thaumatin, and tachystatin families. In addition, we identified a large group of short, cysteine-rich putative AMPs, some of which were strongly elicited. The immunity-related effector gene repertoire also included c-type and i-type lysozymes and several pattern-recognition proteins, such as proteins that recognize Gram-negative bacteria and peptidoglycans. Finally, we identified 45 hemolymph lipopolysaccharide-binding protein sequences, an unusually large number for insects. Taken together, our results indicate that at least some phasmids synthesize a broad spectrum of diverse AMPs that deserve further in-depth analysis.

Insect Defense Proteins and Peptides

The composition of insect hemolymph can change depending on many factors, e.g. access to nutrients, stress conditions, and current needs of the insect. In this chapter, insect immune-related polypeptides, which can be permanently or occasionally present in the hemolymph, are described. Their division into peptides or low-molecular weight proteins is not always determined by the length or secondary structure of a given molecule but also depends on the mode of action in insect immunity and, therefore, it is rather arbitrary. Antimicrobial peptides (AMPs) with their role in immunity, modes of action, and classification are presented in the chapter, followed by a short description of some examples: cecropins, moricins, defensins, proline- and glycine-rich peptides. Further, we will describe selected immune-related proteins that may participate in immune recognition, may possess direct antimicrobial properties, or can be involved in the modulation of insect immunity by both abiotic and biotic factors. We briefly cover Fibrinogen-Related Proteins (FREPs), Down Syndrome Cell Adhesion Molecules (Dscam), Hemolin, Lipophorins, Lysozyme, Insect Metalloproteinase Inhibitor (IMPI), and Heat Shock Proteins. The reader will obtain a partial picture presenting molecules participating in one of the most efficient immune strategies found in the animal world, which allow insects to inhabit all ecological land niches in the world.

Mini Review on Antimicrobial Peptides, Sources, Mechanism and Recent Applications

Antimicrobial peptides in recent years have gained increased interest among scientists, health professionals and the pharmaceutical companies owing to their therapeutic potential. These are low molecular weight proteins with broad range antimicrobial and immuno modulatory activities against infectious bacteria (Gram positive and Gram negative), viruses and fungi. Inability of micro-organisms to develop resistance against most of the antimicrobial peptide has made them as an efficient product which can greatly impact the new era of antimicrobials. In addition to this these peptides also demonstrates increased efficacy, high specificity, decreased drug interaction, low toxicity, biological diversity and direct attacking properties. Pharmaceutical industries are therefore conducting appropriate clinical trials to develop these peptides as potential therapeutic drugs. More than 60 peptide drugs have already reached the market and several hundreds of novel therapeutic peptides are in preclinical and clinical development. Rational designing can be used further to modify the chemical and physical properties of existing peptides. This mini review will discuss the sources, mechanism and recent therapeutic applications of antimicrobial peptides in treatment of infectious diseases.

Antibacterial immune functions of subadults and adults in a semelparous spider

Although capacity to mount an efficient immune response plays a critical role in individuals’ survival, its dynamics across ontogenetic stages is still largely unexplored. Life stage-dependent variation in the encountered diversity and prevalence of parasites were proposed to contribute to stage-dependent changes in immunity, but differences in life history objectives between developmental stages may also lead to stage-specific changes in efficiency of given immune mechanisms. The reason for this is that juveniles and subadults are unable to reproduce, therefore they invest resources mainly into survival, while adults have to partition their resources between survival and reproduction. The general trade-off between somatic maintenance and reproductive effort is expected to impair immune function. Especially so in semelparous organisms that only reproduce once throughout their lifetime, hence they do not face the trade-off between current and future reproduction. We hypothesised that in a semelparous species individuals would be characterised by decreased investment into somatic maintenance after maturation, in order to maximise their reproductive output. Accordingly, we predicted that (1) elements of somatic maintenance, such as immunity, should be relatively weaker in adults in comparison to subadults, and (2) increased reproductive investment in adults should be associated with lower immune efficiency. We quantified two markers of immunity in subadult and adult individuals of the semelparous wolf spider Pardosa agrestis (Westring, 1861), namely bacterial growth inhibition power and bacterial cell wall lytic activity. We found that subadults showed significantly higher cell wall lytic activity than adults, but the two life stages did not differ in their capacity to inhibit bacterial growth. Also, we found weaker immune measures in mated females compared to virgins. Furthermore, in mated females bacterial growth inhibition power was negatively associated with fecundity.

RNA-seq profiles of putative genes involved in specific immune priming in Bombyx mori haemocytes

BACKGROUND

The immune system of many invertebrates, including insects, has been shown to comprise memory, or specific immune priming. However, knowledge of the molecular mechanisms especially the candidate immune-related genes mediated the specificity of the immune priming are still very scarce and fragmentary. We therefore used two closely related Gram-negative pathogenic bacteria (Photorhabdus luminescens TT01 and P. luminescens H06) as the priming agents and employed Illumina/Solexa platform to investigate the transcriptional changes of the haemocytes of Bombyx mori larvae after priming.

RESULTS

In total, 23.0 Gbp of sequence data and 153,331,564 reads were generated, representing 10,496 genes. Approximately 89% of the genes or sequenced reads could be aligned to the silkworm reference genome. The differentially expressed genes (DEGs) of PBS-vs-TT01 (up-regulated expression of TT01 relative to PBS), PBS-vs-H06 (up-regulated expression of H06 relative to PBS) and TT01-vs-H06 (up-regulated expression of H06 relative to TT01) were 707, 159 and 461 respectively. In addition, expression patterns of 25 selected DEGs derived from quantitative real-time polymerase chain reaction (qRT-PCR) were consistent with their transcript abundance changes obtained by transcriptomic analyses. The DEGs are mainly related to pattern recognition receptors (PRRs), antimicrobial peptides (AMPs), signaling molecular, effector molecules, phagosome and spliceosome, indicating that they have participated in the regulation of the specific immune priming in the B. mori larvae.

CONCLUSIONS

The transcriptome profiling data sets from this study will provide valuable resources to better understand the molecular and biological mechanisms regulating the specificity of invertebrates' immune priming. All these will shed light on controlling insect pests or preventing epidemic of infectious diseases in economic invertebrates.

Pectin Digestion in Herbivorous Beetles: Impact of Pseudoenzymes Exceeds That of Their Active Counterparts

Many protein families harbor pseudoenzymes that have lost the catalytic function of their enzymatically active counterparts. Assigning alternative function and importance to these proteins is challenging. Because the evolution toward pseudoenzymes is driven by gene duplication, they often accumulate in multigene families. Plant cell wall-degrading enzymes (PCWDEs) are prominent examples of expanded gene families. The pectolytic glycoside hydrolase family 28 (GH28) allows herbivorous insects to break down the PCW polysaccharide pectin. GH28 in the Phytophaga clade of beetles contains many active enzymes but also many inactive counterparts. Using functional characterization, gene silencing, global transcriptome analyses, and recordings of life history traits, we found that not only catalytically active but also inactive GH28 proteins are part of the same pectin-digesting pathway. The robustness and plasticity of this pathway and thus its importance for the beetle is supported by extremely high steady-state expression levels and counter-regulatory mechanisms. Unexpectedly, the impact of pseudoenzymes on the pectin-digesting pathway in Phytophaga beetles exceeds even the influence of their active counterparts, such as a lowered efficiency of food-to-energy conversion and a prolongation of the developmental period.

Vitellogenin expression in the oilseed rape pest Meligethes aeneus

When investigating insecticide resistance of pest insects, for example, the pollen beetle Meligethes aeneus, it is relevant to differentiate toxicological and molecular genetic data between male and female specimens. A molecular sex determination method would allow resistance testing to be run without prior sorting of the samples. A one-step quantitative RT-PCR method for quantification of the yolk protein vitellogenin expression in the pollen beetle was established. The expression level of vitellogenin relative to tubulin was determined. Pollen beetles were tested at different time points during their development to determine if vitellogenin is a reliable molecular marker for detection of sexually mature females. The differentiation between females and males by relative expression of vitellogenin to tubulin is conditional regarding the life cycle. Sexually mature females and males could easily be distinguished, whereas immature specimens could not be seperated. Vitellogenin expression is a successful marker for identification of sexually mature pollen beetles. Females from the spring populations showed vitellogenin expression when the temperature was above 10.2°C. Further, detailed observations of vitellogenin throughout the spring indicated a strong relationship between daily temperatures and vitellogenin expression, which is an indicator of oviposition ability.

Beetles as Model Organisms in Physiological, Biomedical and Environmental Studies - A Review

Model organisms are often used in biological, medical and environmental research. Among insects, Drosophila melanogaster, Galleria mellonella, Apis mellifera, Bombyx mori, Periplaneta americana, and Locusta migratoria are often used. However, new model organisms still appear. In recent years, an increasing number of insect species has been suggested as model organisms in life sciences research due to their worldwide distribution and environmental significance, the possibility of extrapolating research studies to vertebrates and the relatively low cost of rearing. Beetles are the largest insect order, with their representative - Tribolium castaneum - being the first species with a completely sequenced genome, and seem to be emerging as new potential candidates for model organisms in various studies. Apart from T. castaneum, additional species representing various Coleoptera families, such as Nicrophorus vespilloides, Leptinotarsa decemlineata, Coccinella septempunctata, Poecilus cupreus, Tenebrio molitor and many others, have been used. They are increasingly often included in two major research aspects: biomedical and environmental studies. Biomedical studies focus mainly on unraveling mechanisms of basic life processes, such as feeding, neurotransmission or activity of the immune system, as well as on elucidating the mechanism of different diseases (neurodegenerative, cardiovascular, metabolic, or immunological) using beetles as models. Furthermore, pharmacological bioassays for testing novel biologically active substances in beetles have also been developed. It should be emphasized that beetles are a source of compounds with potential antimicrobial and anticancer activity. Environmental-based studies focus mainly on the development and testing of new potential pesticides of both chemical and natural origin. Additionally, beetles are used as food or for their valuable supplements. Different beetle families are also used as bioindicators. Another important research area using beetles as models is behavioral ecology studies, for instance, parental care. In this paper, we review the current knowledge regarding beetles as model organisms and their practical application in various fields of life science.

The transcriptome analysis of Protaetia brevitarsis Lewis larvae

Larvae of the pest Protaetia brevitarsis are used to treat infections in traditional Chinese medicine. However, genomic information about this non-model species is currently lacking. To better understand the fundamental biology of this non-model species, its transcriptome was obtained using next generation sequencing and then analyzed. A total of 7.62 Gb of clean reads were obtained, which were assembled into 169,087 transcripts corresponding to 142,000 annotated unigenes. These unigenes were functionally classified according to Gene Ontology (GO), euKaryotic Ortholog Groups of proteins (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations. A total of 41,921 unigenes were assigned to 56 GO terms, 21,454 unigenes were divided among 26 KOG categories, and 16,368 unigenes were assigned to 32 KEGG pathways. In addition, 19,144 simple sequence repeats (SSRs) were identified. Furthermore, several kinds of natural antimicrobial peptides and proteins, 4 histones with potential antimicrobial activity, and 41 potential antimicrobial peptide sequences were identified. These data are the first reported whole transcriptome sequence of P. brevitarsis larvae, which represents a valuable genomic resource for studying this species, thus promoting the utilization of its medical potential.

Functional identification of lncRNAs in sweet cherry (Prunus avium) pollen tubes via transcriptome analysis using single-molecule long-read sequencing

Sweet cherry (Prunus avium) is a popular fruit with high nutritional value and excellent flavor. Although pollen plays an important role in the double fertilization and subsequent fruit production of this species, little is known about its pollen tube transcriptome. In this study, we identified 16,409 transcripts using single-molecule sequencing. After filtering 292 transposable elements, we conducted further analyses including mRNA classification, gene function prediction, alternative splicing (AS) analysis, and long noncoding RNA (lncRNA) identification to gain insight into the pollen transcriptome. The filtered transcripts could be matched with 3,438 coding region sequences from the sweet cherry genome. GO and KEGG analyses revealed complex biological processes during pollen tube elongation. A total of 2043 AS events were predicted, 7 of which were identified in different organs, such as the leaf, pistil and pollen tube. Using BLASTnt and the Coding-Potential Assessment Tool (CPAT), we distinguished a total of 284 lncRNAs, among which 154 qualified as natural antisense transcripts (NATs). As the NATs could be the reverse complements of coding mRNA sequences, they might bind to coding sequences. Antisense transfection assays showed that the NATs could regulate the expression levels of their complementary sequences and even affect the growth conditions of pollen tubes. In summary, this research characterizes the transcripts of P. avium pollen and lays the foundation for elucidating the physiological and biochemical mechanisms underlying sexual reproduction in the male gametes of this species.

Novel Factors of Viral Origin Inhibit TOR Pathway Gene Expression

Polydnaviruses (PDVs) are obligate symbionts of endoparasitoid wasps, which exclusively attack the larval stages of their lepidopteran hosts. The Polydnavirus is injected by the parasitoid female during oviposition to selectively infect host tissues by the expression of viral genes without undergoing replication. Toxoneuron nigriceps bracovirus (TnBV) is associated with Toxoneuron nigriceps (Hymenoptera: Braconidae) wasp, an endoparasitoid of the tobacco budworm larval stages, Heliothis virescens (Lepidoptera: Noctuidae). Previous studies showed that TnBV is responsible for alterations in host physiology. The arrest of ecdysteroidogenesis is the main alteration which occurs in last (fifth) instar larvae and, as a consequence, prevents pupation. TnBV induces the functional inactivation of H. virescens prothoracic glands (PGs), resulting in decreased protein synthesis and phosphorylation. Previous work showed the involvement of the PI3K/Akt/TOR pathway in H. virescens PG ecdysteroidogenesis. Here, we demonstrate that this cellular signaling is one of the targets of TnBV infection. Western blot analysis and enzyme immunoassay (EIA) showed that parasitism inhibits ecdysteroidogenesis and the phosphorylation of the two targets of TOR (4E-BP and S6K), despite the stimulation of PTTH contained in the brain extract. Using a transcriptomic approach, we identified viral genes selectively expressed in last instar H. virescens PGs, 48 h after parasitization, and evaluated expression levels of PI3K/Akt/TOR pathway genes in these tissues. The relative expression of selected genes belonging to the TOR pathway (tor, 4e-bp, and s6k) in PGs of parasitized larvae was further confirmed by qRT-PCR. The down-regulation of these genes in PGs of parasitized larvae supports the hypothesis of TnBV involvement in blocking ecdysteroidogenesis, through alterations of the PI3K/Akt/TOR pathway at the transcriptional level.

UDP-glycosyltransferase family in Haemonchus contortus: Phylogenetic analysis, constitutive expression, sex-differences and resistance-related differences

UDP-glycosyltransferases (UGT), catalysing conjugation of UDP-activated sugar donors to small lipophilic chemicals, are widespread in living organisms from bacteria to fungi, plant, or animals. The progress of genome sequencing has enabled an assessment of the UGT multigene family in Haemonchus contortus (family Trichostrongylidae, Nematoda), a hematophagous gastrointestinal parasite of small ruminants. Here we report 32 putative UGT genes divided into 15 UGT families. Phylogenetic analysis in comparison with UGTs from Caenorhabditis elegans, a free-living model nematode, revealed several single member homologues, a lack of the dramatic gene expansion seen in C. elegans, but also several families (UGT365, UGT366, UGT368) expanded in H. contortus only. The assessment of constitutive UGT mRNA expression in H. contortus adults identified significant differences between females and males. In addition, we compared the expression of selected UGTs in the drug-sensitive ISE strain to two benzimidazole-resistant strains, IRE and WR, with different genetic backgrounds. Constitutive expression of UGT368B2 was significantly higher in both resistant strains than in the sensitive strain. As resistant strains were able to deactivate benzimidazole anthelmintics via glycosylation more effectively then the sensitive strain, UGT368B2 enhanced constitutive expression might contribute to drug resistance in H. contortus.

Comparative analysis of the Monochamus alternatus immune system

The pine sawyer beetle, Monochamus alternatus, is regarded as a notorious forest pest in Asia, vectoring an invasive pathogenic nematode, Bursaphelenchus xylophilus, which is known to cause pine wilt disease. However, little sequence information is available for this vector beetle. This hampered the research on its immune system. Based on the transcriptome of M. alternatus, we have identified and characterized 194 immunity-related genes in M. alternatus, and compared them with homologues molecules from other species known to exhibit immune responses against invading microbes. The lower number of putative immunity-related genes in M. alternatus were attributed to fewer C-type lectin, serine protease (SP) and anti-microbial peptide (AMP) genes. Phylogenetic analysis revealed that M. alternatus had a unique recognition gene, galectin3, orthologues of which were not identified in Tribolium castaneum, Drosophila melanogastor, Anopheles gambiae and Apis mellifera. This suggested a lineage-specific gene evolution for coleopteran insects. Our study provides the comprehensive sequence resources of the immunity-related genes of M. alternatus, presenting valuable information for better understanding of the molecular mechanism of innate immunity processes in M. alternatus against B. xylophilus.

Molecular identification and characterization of rhodaneses from the insect herbivore Pieris rapae

The association of cabbage white butterflies (Pieris spec., Lepidoptera: Pieridae) with their glucosinolate-containing host plants represents a well-investigated example of the sequential evolution of plant defenses and insect herbivore counteradaptations. The defensive potential of glucosinolates, a group of amino acid-derived thioglucosides present in plants of the Brassicales order, arises mainly from their rapid breakdown upon tissue disruption resulting in formation of toxic isothiocyanates. Larvae of P. rapae are able to feed exclusively on glucosinolate-containing plants due to expression of a nitrile-specifier protein in their gut which redirects glucosinolate breakdown to the formation of nitriles. The release of equimolar amounts of cyanide upon further metabolism of the benzylglucosinolate-derived nitrile suggests that the larvae are also equipped with efficient means of cyanide detoxification such as β-cyanoalanine synthases or rhodaneses. While insect β-cyanoalanine synthases have recently been identified at the molecular level, no sequence information was available of characterized insect rhodaneses. Here, we identify and characterize two single-domain rhodaneses from P. rapae, PrTST1 and PrTST2. The enzymes differ in their kinetic properties, predicted subcellular localization and expression in P. rapae indicating different physiological roles. Phylogenetic analysis together with putative lepidopteran rhodanese sequences indicates an expansion of the rhodanese family in Pieridae.

Metabolomic and transcriptomic data on major metabolic/biosynthetic pathways in workers and soldiers of the termite Prorhinotermes simplex (Isoptera: Rhinotermitidae) and chemical synthesis of intermediates of defensive (E)-nitropentadec-1-ene biosynthesis

Production of nitro compounds has only seldom been recorded in arthropods. The aliphatic nitroalkene (E)-nitropentadec-1-ene (NPD), identified in soldiers of the termite genus Prorhinotermes, was the first case documented in insects in early seventies. Yet, the biosynthetic origin of NPD has long remained unknown. We previously proposed that NPD arises through the condensation of amino acids glycine and/or l-serine with tetradecanoic acid along a biosynthetic pathway analogous to the formation of sphingolipids. Here, we provide a metabolomics and transcriptomic data of the Prorhinotermes simplex termite workers and soldiers. Data are related to NPD biosynthesis in P. simplex soldiers. Original metabolomics data were deposited in MetaboLights metabolomics database and are become publicly available after publishing the original article. Additionally, chemical synthesis of biosynthetic intermediates of NPD in nonlabeled and stable labeled forms are reported. Data extend our poor knowledge of arthropod metabolome and transcriptome and would be useful for comparative study in termites or other arthropods. The data were used for de-replication of NPD biosynthesis and published separately (Jirošová et al., 2017) [1].

Symbiotic polydnavirus and venom reveal parasitoid to its hyperparasitoids

Symbiotic relationships may provide organisms with key innovations that aid in the establishment of new niches. For example, during oviposition, some species of parasitoid wasps, whose larvae develop inside the bodies of other insects, inject polydnaviruses into their hosts. These symbiotic viruses disrupt host immune responses, allowing the parasitoid's progeny to survive. Here we show that symbiotic polydnaviruses also have a downside to the parasitoid's progeny by initiating a multitrophic chain of interactions that reveals the parasitoid larvae to their enemies. These enemies are hyperparasitoids that use the parasitoid progeny as host for their own offspring. We found that the virus and venom injected by the parasitoid during oviposition, but not the parasitoid progeny itself, affected hyperparasitoid attraction toward plant volatiles induced by feeding of parasitized caterpillars. We identified activity of virus-related genes in the caterpillar salivary gland. Moreover, the virus affected the activity of elicitors of salivary origin that induce plant responses to caterpillar feeding. The changes in caterpillar saliva were critical in inducing plant volatiles that are used by hyperparasitoids to locate parasitized caterpillars. Our results show that symbiotic organisms may be key drivers of multitrophic ecological interactions. We anticipate that this phenomenon is widespread in nature, because of the abundance of symbiotic microorganisms across trophic levels in ecological communities. Their role should be more prominently integrated in community ecology to understand organization of natural and managed ecosystems, as well as adaptations of individual organisms that are part of these communities.

Ecdysteroidogenesis and development in Heliothis virescens (Lepidoptera: Noctuidae): Focus on PTTH-stimulated pathways

Post-embryonic development and molting in insects are regulated by endocrine changes, including prothoracicotropic hormone (PTTH)-stimulated ecdysone secretion by the prothoracic glands (PGs). In Lepidoptera, two pathways are potentially involved in PTTH-stimulated ecdysteroidogenesis, mitogen-activated protein kinase (MAPK) and phosphoinositide 3-kinase/protein kinase B/target of rapamycin (PI3K/Akt/TOR). We investigated the potential roles of both these pathways in Heliothis virescens ecdysteroidogenesis. We identified putative proteins belonging to MAPK and PI3K/Akt/TOR signaling cascades, using transcriptomic analyses of PGs from last (fifth) instar larvae. Using western blots, we measured the phosphorylation of 4E-BP and S6K proteins, the main targets of TOR, following the in vitro exposure of PGs to brain extract containing PTTH (hereafter referred to as PTTH) and/or the inhibitors of MAPK (U0126), PI3K (LY294002) or TOR (rapamycin). Next, we measured ecdysone production, under the same experimental conditions, by enzyme immunoassay (EIA). We found that in Heliothis virescens last instar larvae, both pathways modulated PTTH-stimulated ecdysteroidogenesis. Finally, we analyzed the post-embryonic development of third and fourth instar larvae fed on diet supplemented with rapamycin, in order to better understand the role of the TOR pathway in larval growth. When rapamycin was added to the diet of larvae, the onset of molting was delayed, the growth rate was reduced and abnormally small larvae/pupae with high mortality rates resulted. In larvae fed on diet supplemented with rapamycin, the growth of PGs was suppressed, and ecdysone production and secretion were inhibited. Overall, the in vivo and in vitro results demonstrated that, similarly to Bombyx mori, MAPK and PI3K/Akt/TOR pathways are involved in PTTH signaling-stimulated ecdysteroidogenesis, and indicated the important role of TOR protein in H. virescens systemic growth.

Gene silencing in Tribolium castaneum as a tool for the targeted identification of candidate RNAi targets in crop pests

RNAi shows potential as an agricultural technology for insect control, yet, a relatively low number of robust lethal RNAi targets have been demonstrated to control insects of agricultural interest. In the current study, a selection of lethal RNAi target genes from the iBeetle (Tribolium castaneum) screen were used to demonstrate efficacy of orthologous targets in the economically important coleopteran pests Diabrotica virgifera virgifera and Meligethes aeneus. Transcript orthologs of 50 selected genes were analyzed in D. v. virgifera diet-based RNAi bioassays; 21 of these RNAi targets showed mortality and 36 showed growth inhibition. Low dose injection- and diet-based dsRNA assays in T. castaneum and D. v. virgifera, respectively, enabled the identification of the four highly potent RNAi target genes: Rop, dre4, ncm, and RpII140. Maize was genetically engineered to express dsRNA directed against these prioritized candidate target genes. T₀ plants expressing Rop, dre4, or RpII140 RNA hairpins showed protection from D. v. virgifera larval feeding damage. dsRNA targeting Rop, dre4, ncm, and RpII140 in M. aeneus also caused high levels of mortality both by injection and feeding. In summary, high throughput systems for model organisms can be successfully used to identify potent RNA targets for difficult-to-work with agricultural insect pests.

RNA-sequencing analysis of fungi-induced transcripts from the bamboo wireworm Melanotus cribricollis (Coleoptera: Elateridae) larvae

Larvae of Melanotus cribricollis, feed on bamboo shoots and roots, causing serious damage to bamboo in Southern China. However, there is currently no effective control measure to limit the population of this underground pest. Previously, a new entomopathogenic fungal strain isolated from M. cribricollis larvae cadavers named Metarhizium pingshaense WP08 showed high pathogenic efficacy indoors, indicated that the fungus could be used as a bio-control measure. So far, the genetic backgrounds of both M. cribricollis and M. pingshaense WP08 were blank. Here, we analyzed the whole transcriptome of M. cribricollis larvae, infected with M. pingshaense WP08 or not, using high-throughput next generation sequencing technology. In addition, the transcriptome sequencing of M. pingshaense WP08 was also performed for data separation of those two non-model species. The reliability of the RNA-Seq data was also validated through qRT-PCR experiment. The de novo assembly, functional annotation, sequence comparison of four insect species, and analysis of DEGs, enriched pathways, GO terms and immune related candidate genes were operated. The results indicated that, multiple defense mechanisms of M. cribricollis larvae are initiated to protect against the more serious negative effects caused by fungal infection. To our knowledge, this was the first report of transcriptome analysis of Melanotus spp. infected with a fungus, and it could provide insights to further explore insect-fungi interaction mechanisms.

Rust Infection of Black Poplar Trees Reduces Photosynthesis but Does Not Affect Isoprene Biosynthesis or Emission

Poplar (Populus spp.) trees are widely distributed and play an important role in ecological communities and in forestry. Moreover, by releasing high amounts of isoprene, these trees impact global atmospheric chemistry. One of the most devastating diseases for poplar is leaf rust, caused by fungi of the genus Melampsora. Despite the wide distribution of these biotrophic pathogens, very little is known about their effects on isoprene biosynthesis and emission. We therefore infected black poplar (P. nigra) trees with the rust fungus M. larici-populina and monitored isoprene emission and other physiological parameters over the course of infection to determine the underlying mechanisms. We found an immediate and persistent decrease in photosynthesis during infection, presumably caused by decreased stomatal conductance mediated by increased ABA levels. At the same time, isoprene emission remained stable during the time course of infection, consistent with the stability of its biosynthesis. There was no detectable change in the levels of intermediates or gene transcripts of the methylerythritol 4-phosphate (MEP) pathway in infected compared to control leaves. Rust infection thus does not affect isoprene emission, but may still influence the atmosphere via decreased fixation of CO₂.

Nutritional immunology: Diversification and diet-dependent expression of antimicrobial peptides in the black soldier fly Hermetia illucens

The black soldier fly Hermetia illucens is used for the bioconversion of organic waste into feed for livestock and aquaculture, and is economically among the most important farmed insects in the world. The larvae can be fed on agricultural waste and even liquid manure, resulting in highly unpredictable pathogen levels and dietary conditions. Here we show that H. illucens larvae express a remarkably expanded spectrum of antimicrobial peptides (AMPs), many of which are induced by feeding on a diet containing high bacterial loads. The addition of sulfonated lignin, cellulose, chitin, brewer's grains or sunflower oil revealed the diet-dependent expression profiles of AMPs in the larvae. The highest number of AMPs and the highest levels of AMP expression were induced by feeding larvae on diets supplemented with protein or sunflower oil. Strikingly, the diet-dependent expression of AMPs translated into diet-dependent profiles of inhibitory activities against a spectrum of bacteria, providing an intriguing example for the emerging field of nutritional immunology. We postulate that the fine-tuned expression of the expanded AMP repertoire mediates the adaptation of the gut microbiota to the digestion of unusual diets, and this feature could facilitate the use of H. illucens for the bioconversion of organic waste.

Drastic Genome Reduction in an Herbivore's Pectinolytic Symbiont

Pectin, an integral component of the plant cell wall, is a recalcitrant substrate against enzymatic challenges by most animals. In characterizing the source of a leaf beetle's (Cassida rubiginosa) pectin-degrading phenotype, we demonstrate its dependency on an extracellular bacterium housed in specialized organs connected to the foregut. Despite possessing the smallest genome (0.27 Mb) of any organism not subsisting within a host cell, the symbiont nonetheless retained a functional pectinolytic metabolism targeting the polysaccharide's two most abundant classes: homogalacturonan and rhamnogalacturonan I. Comparative transcriptomics revealed pectinase expression to be enriched in the symbiotic organs, consistent with enzymatic buildup in these structures following immunostaining with pectinase-targeting antibodies. Symbiont elimination results in a drastically reduced host survivorship and a diminished capacity to degrade pectin. Collectively, our findings highlight symbiosis as a strategy for an herbivore to metabolize one of nature's most complex polysaccharides and a universal component of plant tissues.

RNA-seq of Rice Yellow Stem Borer Scirpophaga incertulas Reveals Molecular Insights During Four Larval Developmental Stages

The yellow stem borer (YSB), Scirpophaga incertulas, is a prominent pest in rice cultivation causing serious yield losses. The larval stage is an important stage in YSB, responsible for maximum infestation. However, limited knowledge exists on the biology and mechanisms underlying the growth and differentiation of YSB. To understand and identify the genes involved in YSB development and infestation, so as to design pest control strategies, we performed de novo transcriptome analysis at the first, third, fifth, and seventh larval developmental stages employing Illumina Hi-seq. High-quality reads (HQR) of ∼229 Mb were assembled into 24,775 transcripts with an average size of 1485 bp. Genes associated with various metabolic processes, i.e., detoxification mechanism [CYP450, GSTs, and carboxylesterases (CarEs)], RNA interference (RNAi) machinery (Dcr-1, Dcr-2, Ago-1, Ago-2, Sid-1, Sid-2, Sid-3, and Sid-1-related gene), chemoreception (CSPs, GRs, OBPs, and ORs), and regulators [transcription factors (TFs) and hormones] were differentially regulated during the developmental stages. Identification of stage-specific transcripts made it possible to determine the essential processes of larval development. Comparative transcriptome analysis revealed that YSB has not evolved much with respect to the detoxification mechanism, but showed the presence of distinct RNAi machinery. The presence of strong specific visual recognition coupled with chemosensory mechanisms supports the monophagous nature of YSB. Designed expressed sequenced tags-simple-sequence repeats (EST-SSRs) will facilitate accurate estimation of the genetic diversity of YSB. This is the first report on characterization of the YSB transcriptome and the identification of genes involved in key processes, which will help researchers and industry to devise novel pest control strategies. This study also opens up a new avenue to develop next-generation resistant rice using RNAi or genome editing approaches.

Arabidopsis glucosinolates trigger a contrasting transcriptomic response in a generalist and a specialist herbivore

Phytophagous insects have to deal with toxic defense compounds from their host plants. Although it is known that insects have evolved genes and mechanisms to detoxify plant allochemicals, how specialist and generalist precisely respond to specific secondary metabolites at the molecular level is less understood. Here we studied the larval performance and transcriptome of the generalist moth Heliothis virescens and the specialist butterfly Pieris brassicae feeding on Arabidopsis thaliana genotypes with different glucosinolate (GS) levels. H. virescens larvae gained significantly more weight on the GS-deficient mutant quadGS compared to wild-type (Col-0) plants. On the contrary, P. brassicae was unaffected by the presence of GS and performed equally well on both genotypes. Strikingly, there was a considerable differential gene expression in H. virescens larvae feeding on Col-0 compared to quadGS. In contrast, compared to H. virescens, P. brassicae displayed a much-reduced transcriptional activation when fed on both plant genotypes. Transcripts coding for putative detoxification enzymes were significantly upregulated in H. virescens, along with digestive enzymes and transposable elements. These data provide an unprecedented view on transcriptional changes that are specifically activated by GS and illustrate differential molecular responses that are linked to adaptation to diet in lepidopteran herbivores.

The digestive and defensive basis of carcass utilization by the burying beetle and its microbiota

Insects that use ephemeral resources must rapidly digest nutrients and simultaneously protect them from competitors. Here we use burying beetles (Nicrophorus vespilloides), which feed their offspring on vertebrate carrion, to investigate the digestive and defensive basis of carrion utilization. We characterize gene expression and microbiota composition in the gut, anal secretions, and on carcasses used by the beetles. We find a strict functional compartmentalization of the gut involving differential expression of immune effectors (antimicrobial peptides and lysozymes), as well as digestive and detoxifying enzymes. A distinct microbial community composed of Firmicutes, Proteobacteria and a clade of ascomycetous yeasts (genus Yarrowia) is present in larval and adult guts, and is transmitted to the carcass via anal secretions, where the yeasts express extracellular digestive enzymes and produce antimicrobial compounds. Our results provide evidence of potential metabolic cooperation between the host and its microbiota for digestion, detoxification and defence that extends from the beetle's gut to its nutritional resource.