Transcriptomic immune response of Tenebrio molitor pupae to parasitization by Scleroderma guani

Immune system modulation & virus transmission during parasitism identified by multi-species transcriptomics of a declining insect biocontrol system

BACKGROUND

The Argentine stem weevil (ASW, Listronotus bonariensis) is a significant pasture pest in Aotearoa New Zealand, primarily controlled by the parasitoid biocontrol agent Microctonus hyperodae. Despite providing effective control of ASW soon after release, M. hyperodae parasitism rates have since declined significantly, with ASW hypothesised to have evolved resistance to its biocontrol agent. While the parasitism arsenal of M. hyperodae has previously been investigated, revealing many venom components and an exogenous novel DNA virus Microctonus hyperodae filamentous virus (MhFV), the effects of said arsenal on gene expression in ASW during parasitism have not been examined. In this study, we performed a multi-species transcriptomic analysis to investigate the biology of ASW parasitism by M. hyperodae, as well as the decline in efficacy of this biocontrol system.

RESULTS

The transcriptomic response of ASW to parasitism by M. hyperodae involves modulation of the weevil's innate immune system, flight muscle components, and lipid and glucose metabolism. The multispecies approach also revealed continued expression of venom components in parasitised ASW, as well as the transmission of MhFV to weevils during parasitism and some interrupted parasitism attempts. Transcriptomics did not detect a clear indication of parasitoid avoidance or other mechanisms to explain biocontrol decline.

CONCLUSIONS

This study has expanded our understanding of interactions between M. hyperodae and ASW in a biocontrol system of critical importance to Aotearoa-New Zealand's agricultural economy. Transmission of MhFV to ASW during successful and interrupted parasitism attempts may link to a premature mortality phenomenon in ASW, hypothesised to be a result of a toxin-antitoxin system. Further research into MhFV and its potential role in ASW premature mortality is required to explore whether manipulation of this viral infection has the potential to increase biocontrol efficacy in future.

Molecular Characterization and Functional Analysis of the Dipeptidyl Peptidase IV from Venom of the Ectoparasitoid Scleroderma guani

Dipeptidyl peptidase IV (DPPIV) is a proline-specific serine peptidase that remains poorly investigated in terms of venom composition. Here, we describe the molecular characteristics and possible functions of DPPIV as a major venom component of the ant-like bethylid ectoparasitoid, Scleroderma guani, named SgVnDPPIV. The SgVnDPPIV gene was cloned, which encodes a protein with the conserved catalytic triads and substrate binding sites of mammalian DPPIV. This venom gene is highly expressed in the venom apparatus. Recombinant SgVnDPPIV, produced in Sf9 cells using the baculovirus expression system, has high enzymatic activity, which can be efficiently inhibited by vildagliptin and sitagliptin. Functional analysis revealed that SgVnDPPIV affects genes related to detoxification, lipid synthesis and metabolism, response to stimuli, and ion exchange in pupae of Tenebrio molitor, an envenomated host of S. guani. The present work contributes towards understanding the role of venom DPPIV involved in the interaction between parasitoid wasp and its host.

The Wnt gene family in Tenebrio molitor and other coleopterans

The Wnt gene family is involved in a wide range of developmental processes. Despite its significance, the evolution and function of Wnt genes remain largely unclear. Here, an exhaustive survey of Wnt genes was conducted in Tenebrio molitor and 17 other beetle genomes. A total of 146 Wnt genes were identified, creating a comprehensive coleopteran Wnt gene catalog. Comparative genomics indicates that dynamic evolutionary patterns of Wnt gene loss and duplication occurred in Coleoptera, leading to the diverse Wnt gene repertoire in various beetles. A striking loss of particular Wnt gene subfamilies occurs in Coleoptera. Remarkably, Wnt gene duplication was discovered for the first time in insects. Further analysis of Wnt gene expression in T. molitor indicates that each Wnt gene, including the duplicated ones, has a unique spatial or temporal expression pattern. The current study provides valuable insight into the evolution and functional validation of Wnt genes in Coleoptera.

Fossoriality in desert-adapted tenebrionid (Coleoptera) larvae

In many extreme arid ecosystems, insects constitute major faunal components and are key contributors in nutrient cycling. Previous research on xerophily in insects has focused on adult forms. This study investigates skeletomuscular and behavioural adaptations of the Kalahari sandworm beetle larvae (Gonopus tibialis Fabricius) for dwelling in the sand. Microcomputed tomography enabled cuticle thickness distribution analysis, revealing structural reinforcements of the mandibular edge, the middle part of the head, and the ventral side of the front legs. Laboratory observations and the analysis of muscular system allowed for the definition and functional description of the elements of the digging apparatus of the sandworm larvae. Obtained results point to the crucial role of the head and mandibles in the digging process. These observations are important for understanding desert ecology and pose a challenge to develop newer excavation techniques.

Analysis of Genetic Diversity of Banana Weevils (Cosmopolites sordidus) (Coleoptera: Curculionidae) Using Transcriptome-Derived Simple Sequence Repeat Markers

The banana weevil, Cosmopolites sordidus (Germar) (Coleoptera: Curculionidae) is an economically important insect pest of bananas. It causes up to 100% yield losses and substantial lifespan reduction in bananas. Advances in genomics, proteomics, and sequencing technologies have provided powerful pathways to genotyping disastrous pests such as C. sordidus. However, such technologies are often not available to the majority of rural subtropical African banana growers and pest control managers. This study was therefore motivated by the need to create cheap and easily accessible C. sordidus genotyping methods that could be deployed by banana pest control managers to the benefit of C. sordidus control programs in the tropics where such advanced technologies are not readily accessible. We used an in-house C. sordidus transcriptome from the an-ongoing study from which we mined an array of simple sequence repeat (SSR) markers. Of these, six highly polymorphic transcriptome-derived SSR markers were used to successfully genotype within and among banana weevil population genetic diversity of 12 C. sordidus populations collected from four banana-growing agro-ecological zones (AEZs) in Uganda. The developed transcriptome-derived SSR markers can be used by researchers in population genetics for characterization of the C. sordidus and identification of new genes that are linked to traits of particular interest. The significant genetic diversity revealed in C. sordidus provides pertinent information for integrated pest management strategies.

Morphological Characterisation of Haemocytes in the Mealworm Beetle Tenebrio molitor (Coleoptera, Tenebrionidae)

The immunocompetence of the mealworm beetle Tenebrio molitor has been well investigated at molecular and physiological levels, but information on morphological and functional characteristics of its immune cells (haemocytes) is still scarce and fragmentary. This study provides an updated overview of the morphology of circulating immune cells from mealworm beetle adults, using light and transmission electron microscopy. Based on their affinities for May-Grünwald Giemsa stain, haemocytes were defined as either eosinophilic, basophilic or neutral. Ultrastructural descriptions allowed to detect four main cell types in the haemolymph: prohaemocytes, plasmatocytes, granular cells and oenocytoids. The morphological plasticity of haemocytes and the evidence of mitotic circulating cells, intermediate cell stages, as well as autophagic activities suggest haemocyte proliferation, turnover and transdifferentiation as constantly active processes in the haemolymph. Cytochemical tests revealed differences in the distribution of carbohydrates among cell types underling the great plasticity of the immune response and the direct involvement of circulating immune cells in the resource allocation. In addition, our results provide a detailed morphological description of vesicle trafficking, macro- and microautophagy, apoptotic and necrotic processes, confirming the suitability of T. molitor haemocytes as a model for studying evolutionarily conserved cellular mechanisms.

Recycling of spent mushroom substrate: Utilization as feed material for the larvae of the yellow mealworm Tenebrio molitor (Coleoptera: Tenebrionidae)

Spent mushroom substrate is made from the waste remaining after the harvest of mushrooms. Here, we evaluated the potential of five spent edible fungi (Auricularia cornea, Lentinus edodes, Pleurotus eryngii, P. citrinopileatus and P. ostreatus) substrates as feed sources for Tenebrio molitor larvae. Young larvae did not survive on any substrate except the spent L. edodes substrate (36.7%). The survival rates in young larvae were similar among the different diets in which wheat bran or rice bran was replaced with 0, 20, 30, 40, 50, or 60% spent L. edodes substrate. The weights of the surviving larvae were decreased only when 70% of wheat bran and > 40% of rice bran was replaced with spent L. edodes substrate. In addition, the middle-aged larvae fed wheat bran only were significantly larger than those fed diets with 30~60% spent L. edodes substrate in dry feed, but the larvae of all treatments failed to pupate. Whereas the green feed was added in dry feed, there were no significant differences in pupal weight, pupation rate, pupal duration, adult emergence, or deformed adults among the three treatments in middle-aged larvae that were fed on diets containing 0, 30, or 40% spent L. edodes substrate. Collectively, these results suggest that spent L. edodes substrate has considerable potential to be used as a partial replacement (< 40%) of conventional feed for T. molitor, and spent mushroom substrate waste may be recycled as feed material for resource insects.

Exposure of Agaricus bisporus to Trichoderma aggressivum f. europaeum leads to growth inhibition and induction of an oxidative stress response

Green mould disease of mushroom, Agaricus bisporus,is caused by Trichodermaspecies and can result in substantial crop losses.Label free proteomic analysis of changes in the abundance of A. bisporusproteins following exposure to T. aggressivumsupernatantin vitroindicated increased abundance of proteins associated with an oxidative stress response (zinc ion binding (+6.6 fold); peroxidase activity (5.3-fold); carboxylic ester hydrolase (+2.4 fold); dipeptidase (+3.2 fold); [2Fe-2S] cluster assembly (+3.3 fold)). Proteins that decreased in relative abundance were associated with growth: structural constituent of ribosome, translation (-12 fold), deadenylation-dependent decapping of nuclear-transcribed mRNA (-3.4 fold), and small GTPase mediated signal transduction (-2.6 fold). In vivoanalysis revealed that 10^-4 T. aggressivuminoculum decreased the mushroom yield by 29% to 56% and 10^-3 T. aggressivuminoculum decreased the mushroom yield by 68% to 100%. Proteins that increased in abundance in A. bisporusin vivofollowing exposure to T. aggressivumindicated an oxidative stress response and included proteins with pyruvate kinase activity (+2.6 fold) and hydrolase activity (+2.1 fold)). The results indicate that exposure of A. bisporusmycelium to T. aggressivum in vitroand in vivoresulted in an oxidative stress response and reduction in growth.

Molecular Evidence that Lysiphlebia japonica Regulates the Development and Physiological Metabolism of Aphis gossypii

Lysiphlebia japonica Ashmead (Hymenoptera, Braconidae) is an endophagous parasitoid and Aphis gossypii Glover (Hemiptera, Aphididae) is a major pest in cotton. The relationship between insect host-parasitoids and their hosts involves complex physiological, biochemical and genetic interactions. This study examines changes in the development and physiological metabolism of A. gossypii regulated by L. japonica. Our results demonstrated that both the body length and width increased compared to non-parasitized aphids. We detected significantly increases in the developmental period as well as severe reproductive castration following parasitization by L. japonica. We then used proteomics to characterize these biological changes, and when combined with transcriptomes, this analysis demonstrated that the differential expression of mRNA (up or downregulation) captured a maximum of 48.7% of the variations of protein expression. We assigned these proteins to functional categories that included immunity, energy metabolism and transport, lipid metabolism, and reproduction. We then verified the contents of glycogen and 6-phosphate glucose, which demonstrated that these important energy sources were significantly altered following parasitization. These results uncover the effects on A. gossypii following parasitization by L. japonica, additional insight into the mechanisms behind insect-insect parasitism, and a better understanding of host-parasite interactions.

Venom serine proteinase homolog of the ectoparasitoid Scleroderma guani impairs host phenoloxidase cascade

The ant-like bethylid ectoparasitoid Scleroderma guani (Hymenoptera: Bethylidae) envenomates host to suppress immune response. Yet, the roles of its venom in inhibiting melanization of the host hemolymph have not been fully characterized. Here, we demonstrated that S. guani envenomation induced strong inhibition of melanization of the hemolymph from Tenebrio molitor (Coleoptera: Tenebrionidae), permitting the successful development of parasitoid offspring. To reveal venom component associated with such function, a serine proteinase homolog (SguaSPH) rich in the venom of S. guani was characterized. It was found that one of the catalytic triad residues for serine proteinase is absent in the amino acid sequence of SguaSPH. This venom component was abundantly expressed in venom apparatus and adult stages. By enzymatic assays, SguaSPH displayed low trypsin and no chymotrypsin activity, and was able to inhibit phenoloxidase activity in the hemolymph of Ostrinia furnacalis (Lepidoptera: Crambidae). The findings suggest that SguaSPH is essential for interfering with hemolymph melanization of S. guani envenomated host via phenoloxidase cascade disruption.

Panaxadiol saponins treatment caused the subtle variations in the global transcriptional state of Asiatic corn borer, Ostrinia furnacalis

Background

Methods

Results

Conclusion

Effects of targeting eye color in Tenebrio molitor through RNA interference of tryptophan 2,3-dioxygenase (vermilion): Implications for insect farming

The gene vermilion encodes tryptophan 2,3-dioxygenase, part of the ommochrome pathway, and is responsible for the dark pigmented eyes in some insects, including beetles. Using RNA interference, we targeted the vermilion gene ortholog in embryos and pupae of the yellow mealworm, Tenebrio molitor, resulting in larvae and adults, respectively, that lacked eye pigment. RNA-Seq was used to analyze the impact of vermilion-specific RNA interference on gene expression. There was a 425-fold reduction in vermilion gene expression (p = 0.0003), as well as significant (p < 0.05) differential expression of 109 other putative genes, most of which were downregulated. Enrichment analysis of Gene Ontology terms found in the differentially expressed data set included genes known to be involved in the ommochrome pathway. However, enrichment analysis also revealed the influence of vermilion expression on genes involved in protein translocation to the endoplasmic reticulum, signal transduction, G-protein-coupled receptor signaling, cell-cycle arrest, mannose biosynthesis, and vitamin transport. These data demonstrate that knockdown of vermilion in T. molitor results in complete loss of eye color (white-eyed phenotype) and identify other interrelated genes in the vermilion metabolic pathway. Therefore, a dominant marker system based on eye color can be developed for the genetic manipulation of T. molitor to increase the value of mealworms as an alternative food source by decreasing negative traits, such as disease susceptibility, and increasing desired traits, such as protein content and vitamin production.

Immune Defenses of a Beneficial Pest: The Mealworm Beetle, Tenebrio molitor

The mealworm beetle, Tenebrio molitor, is currently considered as a pest when infesting stored grains or grain products. However, mealworms are now being promoted as a beneficial insect because their high nutrient content makes them a viable food source and because they are capable of degrading polystyrene and plastic waste. These attributes make T. molitor attractive for mass rearing, which may promote disease transmission within the insect colonies. Disease resistance is of paramount importance for both the control and the culture of mealworms, and several biotic and abiotic environmental factors affect the success of their anti-parasitic defenses, both positively and negatively. After providing a detailed description of T. molitor's anti-parasitic defenses, we review the main biotic and abiotic environmental factors that alter their presentation, and we discuss their implications for the purpose of controlling the development and health of this insect.

Analysis of the Gynaephora qinghaiensis pupae immune transcriptome in response to parasitization by Thektogaster sp

As a pest on the Qinghai-Tibet Plateau, Gynaephora qinghaiensis causes severe damage to grassland vegetation and its pupae are also natural hosts of Thektogaster sp. To successfully parasitize, endoparasitoids generally introduce or secrete multiple parasitic factors into the host body during the spawning stage to suppress the host immune response. To study the parasitic effects of Thektogaster sp. on G. qinghaiensis, a transcriptome analysis of immune-related genes in parasitized and nonparasitized G. qinghaiensis pupae was performed. A total of 371,260,704 clean reads were assembled into 118,144 unigenes with an average length of 884.33 base pairs. Of these, 23,660 unigenes were annotated in at least one database and 94,484 unigenes were not annotated in any databases. These findings indicated that the majority of the genetic resources (79.97% of all unigenes) in Gynaephora should be further explored. Parasitization significantly affected the transcriptional profile of G. qinghaiensis pupae. The present study identified 12,322 differentially expressed genes and 57 immune-related genes were identified in parasitized G. qinghaiensis pupae. Most immune-related genes were downregulated, potentially resulting from the inhibitory effect of Thektogaster sp. on G. qinghaiensis pupae after parasitization. Overall, the transcriptome analysis sheds valuable light on the molecular mechanisms of G. qinghaiensis parasitization by Thektogaster sp. and promotes the development of novel biocontrol strategies for Gynaephora based on immune defense.

Optimized pupal age of Tenebrio molitor L. (Coleoptera: Tenebrionidae) enhanced mass rearing efficiency of Chouioia cunea Yang (Hymenoptera: Eulophidae)

Chouioia cunea Yang (Hymenoptera: Eulophidae) has been widely used for biological control of the fall webworm, Hyphantria cunea (Drury) (Lepidoptera: Arctiidae), in China. The yellow mealworm, Tenebrio molitor L. (Coleoptera: Tenebrionidae), an important resource insect species distributed worldwide, is considered to be a potential alternative host for mass rearing of C. cunea to the Chinese oak silkworm, Antheraea pernyi (Guerin-Meneville) (Lepidoptera: Saturniidae), which is currently used. In this study, we investigated the effects of host age on C. cunea mass rearing by measuring parasitism, development and adult fertility of C. cunea on T. molitor pupae of different ages. The results showed no significant differences in the percentage of parasitized hosts and developmental time of C. cunea in pupae of different ages. However, the number of C. cunea adults (137.2–154.7 adults per host) that emerged from 0, 1, and 2-day-old pupae was significantly higher than that from 4-day-old pupae. The lowest percentages of unemerged adults were found in 2-day-old (1.2%) and 3-day-old (1.4%) pupae, which were significantly lower than that of 4-day-old pupae (10.3%). The emergence of adult females from 0 to 2-day-old pupae (120.2–142.3 per pupa) was significantly higher than that from 4-day-old hosts (64.6). Adult females emerging from 2-day-old pupae carried significantly more eggs (258.2 eggs/female) than those from 0 and 1-day-old pupae (178.4–178.9 eggs/female). Our findings indicated that 2-day-old pupae of T. molitor were most suitable to rear C. cunea. Overall, this research provided valuable information to optimize pupae for the mass rearing of C. cunea on host T. molitor.

Superoxide dismutase from venom of the ectoparasitoid Scleroderma guani inhibits melanization of hemolymph

Superoxide dismutase (SOD) known as an important antioxidative stress protein has been recently found in venoms of several parasitoid wasps. However, its functions and characteristics as a virulent factor remain scarcely described. Here, we report the characterization of two venomous SOD genes (SguaSOD1 and SguaSOD3) from the ectoparasitoid, Scleroderma guani. The metal binding sites, cysteine amino acid positions and signature sequences of the SOD family were conserved within SguaSOD1 and SguaSOD3. Relatively high levels of their transcripts were observed in pupae followed a decrease in early adults, after which they had the highest transcriptions, indicating that their productions would be regulated in venom apparatus. Although the two genes showed lower expression in venom apparatus compared to head and thorax, the enzymatic assay revealed that SOD indeed had activity in venom. Further, we showed that recombinant SguaSOD3 suppressed melanization of host hemolymph, implying that this protein used as a virulent factor uniquely impacts the prophenoloxidase cascade.

Parasitism and venom of ectoparasitoid Scleroderma guani impairs host cellular immunity

Venom is a prominently maternal virulent factor utilized by parasitoids to overcome hosts immune defense. With respect to roles of this toxic mixture involved in manipulating hosts immunity, great interest has been mostly restricted to Ichneumonoidea parasitoids associated with polydnavirus (PDV), of which venom is usually considered as a helper component to enhance the role of PDV, and limited Chalcidoidea species. In contrast, little information is available in other parasitoids, especially ectoparasitic species not carrying PDV. The ectoparasitoid Scleroderma guani injects venom into its host, Tenebrio molitor, implying its venom was involved in suppression of hosts immune response for successful parasitism. Thus, we investigated the effects of parasitism and venom of this parasitoid on counteracting the cellular immunity of its host by examining changes of hemocyte counts, and hemocyte spreading and encapsulation ability. Total hemocyte counts were elevated in parasitized and venom-injected pupae. The spreading behavior of both granulocytes and plasmatocytes was impaired by parasitization and venom. High concentration of venom led to more severely increased hemocyte counts and suppression of hemocyte spreading. The ability of hemocyte encapsulation was inhibited by venom in vitro. In addition to immediate effects observed, venom showed persistent interference in hosts cellular immunity. These results indicate that venom alone from S. guani plays a pivotal role in blocking hosts cellular immune response, serving as a regulator that guarantees the successful development of its progenies. The findings provide a foundation for further investigation of the underlying mechanisms in immune inhibitory action of S. guani venom.

De novo Assembly of the Burying Beetle Nicrophorus orbicollis (Coleoptera: Silphidae) Transcriptome Across Developmental Stages with Identification of Key Immune Transcripts

Burying beetles (Nicrophorus spp.) are among the relatively few insects that provide parental care while not belonging to the eusocial insects such as ants or bees. This behavior incurs energy costs as evidenced by immune deficits and shorter life-spans in reproducing beetles. In the absence of an assembled transcriptome, relatively little is known concerning the molecular biology of these beetles. This work details the assembly and analysis of the Nicrophorus orbicollis transcriptome at multiple developmental stages. RNA-Seq reads were obtained by next-generation sequencing and the transcriptome was assembled using the Trinity assembler. Validation of the assembly was performed by functional characterization using Gene Ontology (GO), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Differential expression analysis highlights developmental stage-specific expression patterns, and immunity-related transcripts are discussed. The data presented provides a valuable molecular resource to aid further investigation into immunocompetence throughout this organism's sexual development.

Identification and Characterization of C-type Lectins in Ostrinia furnacalis (Lepidoptera: Pyralidae)

C-type lectins (CTLs) are a large family of calcium-dependent carbohydrate-binding proteins. They function primarily in cell adhesion and immunity by recognizing various glycoconjugates. We identified 14 transcripts encoding proteins with one or two CTL domains from the transcriptome from Asian corn borer, Ostrinia furnacalis (Guenée; Lepidoptera: Pyralidae). Among them, five (OfCTL-S1 through S5) only contain one CTL domain, the remaining nine (OfIML-1 through 9) have two tandem CTL domains. Five CTL-Ss and six OfIMLs have a signal peptide are likely extracellular while another two OfIMLs might be cytoplasmic. Phylogenetic analysis indicated that OfCTL-Ss had 1:1 orthologs in Lepidoptera, Diptera, Coleoptera and Hymenoptera species, but OfIMLs only clustered with immulectins (IMLs) from Lepidopteran. Structural modeling revealed that the 22 CTL domains adopt a similar double-loop fold consisting of β-sheets and α-helices. The key residues for calcium-dependent or independent binding of specific carbohydrates by CTL domains were predicted with homology modeling. Expression profiles assay showed distinct expression pattern of 14 CTLs: the expression and induction were related to the developmental stages and infected microorganisms. Overall, our work including the gene identification, sequence alignment, phylogenetic analysis, structural modeling, and expression profile assay would provide a valuable basis for the further functional studies of O. furnacalis CTLs.

De Novo Assembly and Characterization of the Xenocatantops brachycerus Transcriptome

Grasshoppers are common pests but also have high nutritional and commercial potential. Xenocatantops brachycerus Willemse (Orthoptera: Acrididae) is an economically important grasshopper species that is reared in China. Using the IlluminaHiSeqTM 4000 platform, three transcriptomes of the adult male, adult female, and nymph of X. brachycerus were sequenced. A total of 133,194,848 clean reads were obtained and de novo assembled into 43,187 unigenes with an average length of 964 bp (N50 of 1799 bp); of these, 24,717 (57.23%) unigenes matched known proteins. Based on these annotations, many putative transcripts related to X. brachycerus growth, development, environmental adaptability, and metabolism of nutritional components and bioactive components were identified. In addition, the expression profiles of all three transcriptome datasets were analyzed, and many differentially expressed genes were detected using RSEM and PossionDis. Unigenes. Unigenes with functions associated with growth and development exhibited higher transcript levels at the nymph stage, and unigenes associated with environmental adaptability showed increased transcription in the adults. These comprehensive X. brachycerus transcriptomic data will provide a useful molecular resource for gene prediction, molecular marker development, and studies on signaling pathways in this species and will serve as a reference for the efficient use of other grasshoppers.

Induction of stress- and immune-associated genes in the Indian meal moth Plodia interpunctella against envenomation by the ectoparasitoid Bracon hebetor

Envenomation is an important process in parasitism by parasitic wasps; it suppresses the immune and development of host insects. However, the molecular mechanisms of host responses to envenomation are not yet clear. This study aimed to determine the transcription-level responses of the Indian meal moth Plodia interpunctella against envenomation of the ectoparasitoid Bracon hebetor. Quantitative real-time reverse-transcription PCR was used to determine the transcriptional changes of 13 selected genes, which are associated with development, metabolism, stress, or immunity, in the feeding and wandering fifth instar larvae over a 4-day period after envenomation. The effects of envenomation on the feeding-stage larvae were compared with those of starvation in the transcriptional levels of the 13 genes. Most selected genes were altered in their expression by either envenomation or starvation. In particular, a heat shock protein, hsp70, was highly upregulated in envenomated larvae in both the feeding and wandering stages as well as in starved larvae. Further, some genes were upregulated by envenomation in a stage-specific manner. For example, hsp25 was upregulated after envenomation in the feeding larvae, but hsp90 and an immune-associated gene, hemolin, were upregulated in the wandering larvae. However, both envenomation and starvation resulted in the downregulation of genes associated with development and metabolism. Taken together, P. interpunctella upregulated stress- and immune-responsive genes, but downregulated genes associated with development and metabolism after envenomation. This study provides important information for understanding the molecular mechanisms of host responses to parasitism.

Re-evaluation of insect melanogenesis research: Views from the dark side

Melanins (eumelanin and pheomelanin) are synthesized in insects for several purposes including cuticle sclerotization and color patterning, clot formation, organogenesis, and innate immunity. Traditional views of insect immunity detail the storage of pro-phenoloxidases inside specialized blood cells (hemocytes) and their release upon recognition of foreign bodies. Activated phenoloxidases convert monophenols into reactive quinones in a two-step enzymatic reaction, and until recently, the mechanism of tyrosine hydroxylation remained a mystery. Herein, we present our interpretations of these enzyme-substrate complexes. The resultant melanins are deposited onto the surface of microbes to immobilize, agglutinate, and suffocate them. Phenoloxidase activity and melanin production are not limited to the blood (hemolymph) or cuticle, as recent evidence points to more diverse, sophisticated interactions in the gut and with the resident symbionts. This review offers insight into the somewhat neglected areas of insect melanogenesis research, particularly in innate immunity, its role in beneficial insects such as pollinators, the functional versatility of phenoloxidases, and the limitations of common experimental approaches that may impede progress inadvertently.

Endogenous egg immune defenses in the yellow mealworm beetle (Tenebrio molitor)

In order to survive microbe encounters, insects rely on both physical barriers as well as local and systemic immune responses. Most research focusses on adult or larval defenses however, whereas insect eggs are also in need of protection. Lately, the defense of eggs against microbes has received an increasing amount of attention, be it through endogenous egg defenses, trans-generational immune priming (TGIP) or parental investment. Here we studied the endogenous immune response in eggs and adults of Tenebrio molitor. We show that many immune genes are induced in both adults and eggs. Furthermore, we show that eggs reach comparable levels of immune gene expression as adults. These findings show that the eggs of Tenebrio are capable of an impressive endogenous immune response, and indicate that such inducible egg defenses are likely common in insects.

Oral immune priming with Bacillus thuringiensis induces a shift in the gene expression of Tribolium castaneum larvae

Background

The phenomenon of immune priming, i.e. enhanced protection following a secondary exposure to a pathogen, has now been demonstrated in a wide range of invertebrate species. Despite accumulating phenotypic evidence, knowledge of its mechanistic underpinnings is currently very limited. Here we used the system of the red flour beetle, Tribolium castaneum and the insect pathogen Bacillus thuringiensis (Bt) to further our molecular understanding of the oral immune priming phenomenon. We addressed how ingestion of bacterial cues (derived from spore supernatants) of an orally pathogenic and non-pathogenic Bt strain affects gene expression upon later challenge exposure, using a whole-transcriptome sequencing approach.

Results

Whereas gene expression of individuals primed with the orally non-pathogenic strain showed minor changes to controls, we found that priming with the pathogenic strain induced regulation of a large set of distinct genes, many of which are known immune candidates. Intriguingly, the immune repertoire activated upon priming and subsequent challenge qualitatively differed from the one mounted upon infection with Bt without previous priming. Moreover, a large subset of priming-specific genes showed an inverse regulation compared to their regulation upon challenge only.

Conclusions

Our data demonstrate that gene expression upon infection is strongly affected by previous immune priming. We hypothesise that this shift in gene expression indicates activation of a more targeted and efficient response towards a previously encountered pathogen, in anticipation of potential secondary encounter.

Electronic supplementary material

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

Identification and expression of cuticular protein genes based on Locusta migratoria transcriptome

Many types of cuticular proteins are found in a single insect species, and their number and features are very diversified among insects. The cuticle matrix consists of many different proteins that confer the physical properties of the exoskeleton. However, the number and properties of cuticle proteins in Locusta migratoria remain unclear. In the present study, Illumina sequencing and de novo assembly were combined to characterize the transcriptome of L. migratoria. Eighty-one cuticular protein genes were identified and divided into five groups: the CPR family (51), Tweedle (2), CPF/CPFLs (9), CPAP family (9), and other genes (10). Based on the expression patterns in different tissues and stages, most of the genes as a test were distributed in the integument, pronotum and wings, and expressed in selected stages with different patterns. The results showed no obvious correlation between the expression patterns and the conservative motifs. Additionally, each cluster displayed a different expression pattern that may possess a different function in the cuticle. Furthermore, the complexity of the large variety of genes displayed differential expression during the molting cycle may be associated with cuticle formation and may provide insights into the gene networks related to cuticle formation.

De novo transcriptome analysis and microsatellite marker development for population genetic study of a serious insect pest, Rhopalosiphum padi (L.) (Hemiptera: Aphididae)

The bird cherry-oat aphid, Rhopalosiphum padi (L.), is one of the most abundant aphid pests of cereals and has a global distribution. Next-generation sequencing (NGS) is a rapid and efficient method for developing molecular markers. However, transcriptomic and genomic resources of R. padi have not been investigated. In this study, we used transcriptome information obtained by RNA-Seq to develop polymorphic microsatellites for investigating population genetics in this species. The transcriptome of R. padi was sequenced on an Illumina HiSeq 2000 platform. A total of 114.4 million raw reads with a GC content of 40.03% was generated. The raw reads were cleaned and assembled into 29,467 unigenes with an N50 length of 1,580 bp. Using several public databases, 82.47% of these unigenes were annotated. Of the annotated unigenes, 8,022 were assigned to COG pathways, 9,895 were assigned to GO pathways, and 14,586 were mapped to 257 KEGG pathways. A total of 7,936 potential microsatellites were identified in 5,564 unigenes, 60 of which were selected randomly and amplified using specific primer pairs. Fourteen loci were found to be polymorphic in the four R. padi populations. The transcriptomic data presented herein will facilitate gene discovery, gene analyses, and development of molecular markers for future studies of R. padi and other closely related aphid species.

Differential expression of genes in the alate and apterous morphs of the brown citrus aphid, Toxoptera citricida

Winged and wingless morphs in insects represent a trade-off between dispersal ability and reproduction. We studied key genes associated with apterous and alate morphs in Toxoptera citricida (Kirkaldy) using RNAseq, digital gene expression (DGE) profiling, and RNA interference. The de novo assembly of the transcriptome was obtained through Illumina short-read sequencing technology. A total of 44,199 unigenes were generated and 27,640 were annotated. The transcriptomic differences between alate and apterous adults indicated that 279 unigenes were highly expressed in alate adults, whereas 5,470 were expressed at low levels. Expression patterns of the top 10 highly expressed genes in alate adults agreed with wing bud development trends. Silencing of the lipid synthesis and degradation gene (3-ketoacyl-CoA thiolase, mitochondrial-like) and glycogen genes (Phosphoenolpyruvate carboxykinase [GTP]-like and Glycogen phosphorylase-like isoform 2) resulted in underdeveloped wings. This suggests that both lipid and glycogen metabolism provide energy for aphid wing development. The large number of sequences and expression data produced from the transcriptome and DGE sequencing, respectively, increases our understanding of wing development mechanisms.

Deep Sequencing-Based Transcriptome Analysis Reveals the Regulatory Mechanism of Bemisia tabaci (Hemiptera: Aleyrodidae) Nymph Parasitized by Encarsia sophia (Hymenoptera: Aphelinidae)

The whitefly Bemisia tabaci is a genetically diverse complex with multiple cryptic species, and some are the most destructive invasive pests of many ornamentals and crops worldwide. Encarsia sophia is an autoparasitoid wasp that demonstrated high efficiency as bio-control agent of whiteflies. However, the immune mechanism of B. tabaci parasitization by E. sophia is unknown. In order to investigate immune response of B. tabaci to E. Sophia parasitization, the transcriptome of E. sophia parasitized B. tabaci nymph was sequenced by Illumina sequencing. De novo assembly generated 393,063 unigenes with average length of 616 bp, in which 46,406 unigenes (15.8% of all unigenes) were successfully mapped. Parasitization by E. sophia had significant effects on the transcriptome profile of B. tabaci nymph. A total of 1482 genes were significantly differentially expressed, of which 852 genes were up-regulated and 630 genes were down-regulated. These genes were mainly involved in immune response, development, metabolism and host signaling pathways. At least 52 genes were found to be involved in the host immune response, 33 genes were involved in the development process, and 29 genes were involved in host metabolism. Taken together, the assembled and annotated transcriptome sequences provided a valuable genomic resource for further understanding the molecular mechanism of immune response of B. tabaci parasitization by E. sophia.

Deciphering the main venom components of the ectoparasitic ant-like bethylid wasp, Scleroderma guani

Similar to venom found in most venomous animals, parasitoid venoms contain a complex cocktail of proteins with potential agrichemical and pharmaceutical use. Even though parasitoids are one of the largest group of venomous animals, little is known about their venom composition. Recent few studies revealed high variated venom composition existing not only in different species but also between closely related strains, impling that increasing information on the venom proteins from more greater diversity of species of different taxa is key to comprehensively uncover the complete picture of parasitoid venom. Here, we explored the major protein components of the venom of ectoparasitic ant-like bethylid wasp, Scleroderma guani by an integrative transcriptomic-proteomic approach. Illumina deep sequencing of venom apparatus cDNA produced 49,873 transcripts. By mapping the peptide spectral data derived from venom reservoir against these transcripts, mass spectrometry analysis revealed ten main venom proteins, including serine proteinase, metalloprotease, dipeptidyl peptidase IV, esterase, antithrombin-III, acid phosphatase, neural/ectodermal development factor IMP-L2 like protein, venom allergen 3, and unknown protein. Interestingly, one serine proteinase was firstly identified with rarely high molecular weight about 200 kDa in parasitoid venom. The occurrence of abundant acid phosphatase, antithrombin-III and venom allergen 3 demonstrated that S. guani venom composition is similar to that of social wasp venoms. All identified venom genes showed abundantly biased expression in venom apparatus, indicating their virulent functions involved in parasitization. This study shed light on the more better understanding of parasitoid venom evolution across species and will facilitate the further elucidation of function and toxicity of these venom proteins.

A de novo transcriptomic analysis to reveal functional genes in Apolygus lucorum

The widespread planting of genetically engineered cotton producing the Cry1Ac toxin has led to significantly reduced pesticide applications since 1997. However, consequently, the number of green mirid bugs (GMB), Apolygus lucorum (Meyer-Dür) has increased. So far the GMB, instead of the cotton bollworm Helicoverpa armigera (Hübner), has become the major pest in the transgenic Bt cotton field and has influenced cotton yield. Disproportionately, only a few studies on GMB at molecular level have been reported. Libraries from both third instar nymphs and adults were sequenced using Illumina technology, producing more than 106 million short reads and assembled into 63 029 unigenes of mean length 597 nt and N50 813 nt, ranging from 300 nt to 9771 nt. BLASTx analysis against Nr, Swissprot, GO and COG was performed to annotate these unigenes. As a result, 26 478 unigenes (42.01%) matched to known proteins and 107 immune-related, 320 digestive-related and 53 metamorphosis-related genes were detected in these annotated unigenes. Additionally, we profiled gene expression using mapping based differentially expressed genes (DEGs) strategy between the two developmental stages: nymph and adult. The results demonstrated that thousands of genes were significantly differentially expressed at different developmental stages. The transcriptome and gene expression data provided comprehensive and global gene resources of GMB. This transcriptome would improve our understanding of the molecular mechanisms of various underlying biological characteristics, including development, digestion and immunity in GMB. Therefore, these findings could help elucidate the intrinsic factors of the GMB resurgence, offering novel pest management targets for future transgenic cotton breeding.

Ixodes scapularis and Ixodes ricinus tick cell lines respond to infection with tick-borne encephalitis virus: transcriptomic and proteomic analysis

BACKGROUND

Ixodid ticks are important vectors of a wide variety of viral, bacterial and protozoan pathogens of medical and veterinary importance. Although several studies have elucidated tick responses to bacteria, little is known about the tick response to viruses. To gain insight into the response of tick cells to flavivirus infection, the transcriptomes and proteomes of two Ixodes spp cell lines infected with the flavivirus tick-borne encephalitis virus (TBEV) were analysed.

METHODS

RNA and proteins were isolated from the Ixodes scapularis-derived cell line IDE8 and the Ixodes ricinus-derived cell line IRE/CTVM19, mock-infected or infected with TBEV, on day 2 post-infection (p.i.) when virus production was increasing, and on day 6 p.i. when virus production was decreasing. RNA-Seq and mass spectrometric technologies were used to identify changes in abundance of, respectively, transcripts and proteins. Functional analyses were conducted on selected transcripts using RNA interference (RNAi) for gene knockdown in tick cells infected with the closely-related but less pathogenic flavivirus Langat virus (LGTV).

RESULTS

Differential expression analysis using DESeq resulted in totals of 43 and 83 statistically significantly differentially-expressed transcripts in IDE8 and IRE/CTVM19 cells, respectively. Mass spectrometry detected 76 and 129 statistically significantly differentially-represented proteins in IDE8 and IRE/CTVM19 cells, respectively. Differentially-expressed transcripts and differentially-represented proteins included some that may be involved in innate immune and cell stress responses. Knockdown of the heat-shock proteins HSP90, HSP70 and gp96, the complement-associated protein Factor H and the protease trypsin resulted in increased LGTV replication and production in at least one tick cell line, indicating a possible antiviral role for these proteins. Knockdown of RNAi-associated proteins Argonaute and Dicer, which were included as positive controls, also resulted in increased LGTV replication and production in both cell lines, confirming their role in the antiviral RNAi pathway.

CONCLUSIONS

This systems biology approach identified several molecules that may be involved in the tick cell innate immune response against flaviviruses and highlighted that ticks, in common with other invertebrate species, have other antiviral responses in addition to RNAi.

THE ENDOPARASITOID Pteromalus puparum INFLUENCES HOST GENE EXPRESSION WITHIN FIRST HOUR OF PARASITIZATION

The small cabbage butterfly, Pieris rapae, is an important pest of cruciferous corps, and Pteromalus puparum is a predominant pupal endoparasitoid wasp of this butterfly. For successful development of parasitoid offspring, female parasitoids usually introduce one or several kinds of maternal factors into the hemocoels during oviposition to suppress host immunity. To investigate the early changes in host immune-related genes following parasitization, we analyzed transcriptomes of parasitized and unparasitized, control, host pupae. Approximately 17.7 and 19.3 million paired-end reads were generated from nonparasitized and parasitized host pupae, and assembled de novo into 45,639 transcripts and 27,659 nonredundant unigenes. The average unigene length was 790 bp. A total 18,377 of 27,659 unigenes were annotated and we identified 557 differentially expressed unigenes in host pupae at 1 h after parasitization, of which 21 were immune-related. Parasitization led to downregulation of most pattern recognition receptors and upregulation of all serine protease inhibitors. The transcirptomic profile of P. rapae is considerably affected by parasitization. This study provides valuable sources for future investigations of the molecular interaction between P. puparum and its host P. rapae.

The Last Breath: A μCT-based method for investigating the tracheal system in Hexapoda

In recent years, μCT-based studies of the insect tracheal system have become an increasingly important area of research. Nevertheless, the methods proposed in previous research for investigating the respiratory system in the three-dimensional space were described and tested based on a relatively small number of specimens. Additionally, the individuals studied in all these cases represented only a single post-embryonic stadium - pupa or imago - of a particular insect species. Therefore, in the current situation it is difficult to predict the reliability and possible limitations of these methods. To address this problem we conducted a methodological study, during which we used 65 individuals representing larvae, pupae and imagines of the mealworm beetle (Tenebrio molitor). In addition to the protocol previously described, which implicated freezing as a killing technique, we also tested a novel one, which was based on ethyl acetate fumigation of the specimens studied. We included step-by-step guides for the manual and semiautomatic approaches in order to facilitate the digital visualization of the tracheal system. Our investigations enabled us to generate multiple models of the tracheal system of all post-embryonic stages of the mealworm beetle. The methods used proved to be minimally invasive, thus allowing for the application of post-scanning manipulations, such as drying with critical point dryer (CPD). This approach enabled us to merge different three-dimensional models into a single picture and analyse the relationship of the tracheal system with other tissues (e.g., muscles, nervous system). We comprehensively discuss the advantages and possible limitations of the tested methods and provide practical suggestions for conducting the analyses on a wider scale. The visualizations presented in this publication are the first three-dimensional models of the respiratory system using a representative of the extremely diverse order Coleoptera.

The extraembryonic serosa is a frontier epithelium providing the insect egg with a full-range innate immune response

Drosophila larvae and adults possess a potent innate immune response, but the response of Drosophila eggs is poor. In contrast to Drosophila, eggs of the beetle Tribolium are protected by a serosa, an extraembryonic epithelium that is present in all insects except higher flies. In this study, we test a possible immune function of this frontier epithelium using Tc-zen1 RNAi-mediated deletion. First, we show that bacteria propagate twice as fast in serosa-less eggs. Then, we compare the complete transcriptomes of wild-type, control RNAi, and Tc-zen1 RNAi eggs before and after sterile or septic injury. Infection induces genes involved in Toll and IMD-signaling, melanisation, production of reactive oxygen species and antimicrobial peptides in wild-type eggs but not in serosa-less eggs. Finally, we demonstrate constitutive and induced immune gene expression in the serosal epithelium using in situ hybridization. We conclude that the serosa provides insect eggs with a full-range innate immune response.

Protein profiles of Chinese white wax scale, Ericerus pela, at the male pupal stage by high-throughput proteomics

The Chinese white wax scale insect (Ericerus pela) is sexually dimorphic with holometabolous males and hemimetabolous females. Holometabolous insects were assumed to originate from hemimetabolous ancestors. Therefore, the male pupal stage is a major innovation compared with hemimetabolous female insects. Here, the protein profiles of the male pupae were obtained by high-throughput proteomics and analyzed using bioinformatics methods. A total of 1,437 peptides were identified and assigned to 677 protein groups. Most of the proteins had molecular weights below 40 kDa and isoelectric points from 4 to 7. Gene Ontology terms were assigned to 331 proteins, including metabolic process, developmental process, and cellular process. Kyoto Encyclopedia of Genes and Genomes annotations identified 142 pathways and most proteins were assigned to metabolism events. Pathways involved in cell growth and death, signal transduction, folding, and sorting and degradation were also identified. Six proteins that had undergone positive selection were classified into four groups, protein biosynthesis, protein degeneration, signal transduction, and detoxification. Many of the high-abundance proteins were enzymes involved in carbohydrate, lipid, and amino acid metabolism; signal transduction; degradation; and immunization, which indicated that metabolism, disruption, and development occurred intensely at the pupal stage. These processes are closely related to the physiological status of pupae. The results also suggested that these related proteins may be fundamental factors in the formation of pupae. This study describes pupal characterization at the molecular level and provides a basis for further physiological studies.

Complete mitochondrial genome of yellow meal worm (Tenebrio molitor)

The yellow meal worm (Tenebrio molitor L.) is an important resource insect typically used as animal feed additive. It is also widely used for biological research. The first complete mitochondrial genome of T. molitor was determined for the first time by long PCR and conserved primer walking approaches. The results showed that the entire mitogenome of T. molitor was 15 785 bp long, with 72.35% A+T content [deposited in GenBank with accession number KF418153]. The gene order and orientation were the same as the most common type suggested as ancestral for insects. Two protein-coding genes used atypical start codons (CTA in ND2 and AAT in COX1), and the remaining 11 protein-coding genes started with a typical insect initiation codon ATN. All tRNAs showed standard clover-leaf structure, except for tRNA(Ser) (AGN), which lacked a dihydrouridine (DHU) arm. The newly added T. molitor mitogenome could provide information for future studies on yellow meal worm.

The parasitoid, Cotesia flavipes (Cameron) (Hymenoptera: Braconidae), influences food consumption and utilization by larval Diatraea saccharalis (F.) (Lepidoptera: Crambidae)

Parasitoids exploit host insects for food and other resources; they alter host development and physiology to optimize conditions to favor parasitoid development. Parasitoids influence their hosts by injecting eggs, along with a variety of substances, including venoms, polydnaviruses, ovarian fluids, and other maternal factors, into hosts. These factors induce profound changes in hosts, such as behavior, metabolism, endocrine events, and immune defense. Because endoparasitoids develop and consume tissues from within their hosts, it is reasonable to suggest that internal parasitization would also influence host food consumption and metabolism. We report on the effects of parasitism by Cotesia flavipes on the food consumption and utilization of its host, Diatraea saccharalis. Cotesia flavipes reduces the host food consumption, but parasitized larvae considered a unit with their parasitoid's attained the same final weight as the nonparasitized larvae. Nutritional indices, midgut activities of carbohydrases, and trypsin of parasitized and nonparasitized D. saccharalis were assessed. Parasitized larvae had reduced relative food consumption, metabolic and growth rates, coupled with higher efficiency for conversion of the digested, but not ingested, food into body mass. Parasitism also affected food flux through the gut and protein contents in the midgut of parasitized larvae. The activity of α-amylase and trehalase in parasitized host was enhanced in the first day after parasitism relative to control larvae. Saccharase activity remained unchanged during larval development. Trypsin activity was reduced from the fifth to ninth day after parasitism. We argue on the mechanisms involved in host food processing after parasitism.

De novo transcriptome analysis of wing development-related signaling pathways in Locusta migratoria manilensis and Ostrinia furnacalis (Guenée)

Background

Orthopteran migratory locust, Locusta migratoria, and lepidopteran Asian corn borer, Ostrinia furnacalis, are two types of insects undergoing incomplete and complete metamorphosis, respectively. Identification of candidate genes regulating wing development in these two insects would provide insights into the further study about the molecular mechanisms controlling metamorphosis development. We have sequenced the transcriptome of O. furnacalis larvae previously. Here we sequenced and characterized the transcriptome of L. migratoria wing discs with special emphasis on wing development-related signaling pathways.

Methodology/Principal Findings

Illumina Hiseq2000 was used to sequence 8.38 Gb of the transcriptome from dissected nymphal wing discs. De novo assembly generated 91,907 unigenes with mean length of 610 nt. All unigenes were searched against five databases including Nt, Nr, Swiss-Prot, COG, and KEGG for annotations using blastn or blastx algorithm with an cut-off E-value of 10⁻⁵. A total of 23,359 (25.4%) unigenes have homologs within at least one database. Based on sequence similarity to homologs known to regulate Drosophila melanogaster wing development, we identified 50 and 46 potential wing development-related unigenes from L. migratoria and O. furnacalis transcriptome, respectively. The identified unigenes encode putative orthologs for nearly all components of the Hedgehog (Hh), Decapentaplegic (Dpp), Notch (N), and Wingless (Wg) signaling pathways, which are essential for growth and pattern formation during wing development. We investigated the expression profiles of the component genes involved in these signaling pathways in forewings and hind wings of L. migratoria and O. furnacalis. The results revealed the tested genes had different expression patterns in two insects.

Conclusions/Significance

This study provides the comprehensive sequence resource of the wing development-related signaling pathways of L. migratoria. The obtained data gives an insight into better understanding the molecular mechanisms involved in the wing development in L. migratoria and O. furnacalis, two insect species with different metamorphosis types.

Parasitization by Scleroderma guani influences expression of superoxide dismutase genes in Tenebrio molitor

Superoxide dismutase (SOD) is an antioxidant enzyme involved in detoxifying reactive oxygen species. In this study, we identified genes encoding the extracellular and intracellular copper-zinc SODs (ecCuZnSOD and icCuZnSOD) and a manganese SOD (MnSOD) in the yellow mealworm beetle, Tenebrio molitor. The cDNAs for ecCuZnSOD, icCuZnSOD, and MnSOD, respectively, encode 24.55, 15.81, and 23.14 kDa polypeptides, which possess structural features typical of other insect SODs. They showed 20-94% identity to other known SOD sequences from Bombyx mori, Musca domestica, Nasonia vitripennis, Pediculus humanus corporis, and Tribolium castaneum. Expression of these genes was analyzed in selected tissues and developmental stages, and following exposure to Escherichia coli and parasitization by Scleroderma guani. We recorded expression of all three SODs in cuticle, fat body, and hemocytes and in the major developmental stages. Relatively higher expressions were detected in late-instar larvae and pupae, compared to other developmental stages. Transcriptional levels were upregulated following bacterial infection. Analysis of pupae parasitized by S. guani revealed that expression of T. molitor SOD genes was significantly induced following parasitization. We infer that these genes act in immune response and in host-parasitoid interactions.

Parasitization by Scleroderma guani influences protein expression in Tenebrio molitor pupae

Ectoparasitoid wasps deposit their eggs onto the surface and inject venom into their hosts. Venoms are chemically complex and they exert substantial impact on hosts, including permanent or temporary paralysis and developmental arrest. These visible venom effects are due to changes in expression of genes encoding physiologically relevant proteins. While the influence of parasitization on gene expression in several lepidopterans has been reported, the molecular details of parasitoid/beetle relationships remain mostly unknown. This shortcoming led us to pose the hypothesis that envenomation by the ectoparasitic ant-like bethylid wasp Scleroderma guani leads to changes in protein expression in the yellow mealworm beetle Tenebrio molitor. We tested our hypothesis by comparing the proteomes of non-parasitized and parasitized host pupae using iTRAQ-based proteomics. We identified 41 proteins that were differentially expressed (32↑- and 9↓-regulated) in parasitized pupae. We assigned these proteins to functional categories, including immunity, stress and detoxification, energy metabolism, development, cytoskeleton, signaling and others. We recorded parallel changes in mRNA levels and protein abundance in 14 selected proteins following parasitization. Our findings support our hypothesis by documenting changes in protein expression in parasitized hosts.

Genomics: moving behavioural ecology beyond the phenotypic gambit

Researchers studying the adaptive significance of behaviour typically assume that genetic mechanisms will not inhibit evolutionary trajectories, an assumption commonly known as the 'phenotypic gambit'. Although the phenotypic gambit continues to be a useful heuristic for behavioural ecology, here we discuss how genomic methods provide new tools and conceptual approaches that are relevant to behavioural ecology. We first describe how the concept of a genetic toolkit for behaviour can allow behavioural ecologists to synthesize both genomic and ecological information when assessing behavioural adaptation. Then we show how gene expression profiles can be viewed as complex phenotypic measurements, used to (1) predict behaviour, (2) evaluate phenotypic plasticity and (3) devise methods to manipulate behaviour in order to test adaptive hypotheses. We propose that advances in genomics and bioinformatics may allow researchers to overcome some of the logistical obstacles that motivated the inception of the phenotypic gambit. Behavioural ecology and genomics are mutually informative, providing potential synergy that could lead to powerful advances in the field of animal behaviour.

Transcriptome immune analysis of the invasive beetle Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae) parasitized by Tetrastichus brontispae Ferrière (Hymenoptera: Eulophidae)

The beetle Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae) is a serious invasive insect pest of palm plants in southern China, and the endoparasitoid Tetrastichus brontispae Ferrière (Hymenoptera: Eulophidae) is a natural enemy of this pest that exhibits great ability in the biocontrol of O. nipae. For successful parasitism, endoparasitoids often introduce or secrete various virulence factors to suppress host immunity. To investigate the effects of parasitization by T. brontispae on the O. nipae immune system, the transcriptome of O. nipae pupae was analyzed with a focus on immune-related genes through Illumina sequencing. De novo assembly generated 49,919 unigenes with a mean length of 598 bp. Of these genes, 27,490 unigenes (55.1% of all unigenes) exhibited clear homology to known genes in the NCBI nr database. Parasitization had significant effects on the transcriptome profile of O. nipae pupae, and most of these differentially expressed genes were down-regulated. Importantly, the expression profiles of immune-related genes were significantly regulated after parasitization. Taken together, these transcriptome sequencing efforts shed valuable light on the host (O. nipae) manipulation mechanisms induced by T. brontispae, which will pave the way for the development of novel immune defense-based management strategies of O. nipae, and provide a springboard for further molecular analyses, particularly of O. nipae invasion.

De novo assembly and characterization of the global transcriptome for Rhyacionia leptotubula using Illumina paired-end sequencing

Background

The pine tip moth, Rhyacionia leptotubula (Lepidoptera: Tortricidae) is one of the most destructive forestry pests in Yunnan Province, China. Despite its importance, less is known regarding all aspects of this pest. Understanding the genetic information of it is essential for exploring the specific traits at the molecular level. Thus, we here sequenced the transcriptome of R. leptotubula with high-throughput Illumina sequencing.

Methodology/Principal Findings

In a single run, more than 60 million sequencing reads were generated. De novo assembling was performed to generate a collection of 46,910 unigenes with mean length of 642 bp. Based on Blastx search with an E-value cut-off of 10⁻⁵, 22,581 unigenes showed significant similarities to known proteins from National Center for Biotechnology Information (NCBI) non-redundant (Nr) protein database. Of these annotated unigenes, 10,360, 6,937 and 13,894 were assigned to Gene Ontology (GO), Clusters of Orthologous Group (COG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. A total of 5,926 unigenes were annotated with domain similarity derived functional information, of which 55 and 39 unigenes respectively encoding the insecticide resistance related enzymes, cytochrome P450 and carboxylesterase. Using the transcriptome data, 47 unigenes belonging to the typical “stress” genes of heat shock protein (Hsp) family were retrieved. Furthermore, 1,450 simple sequence repeats (SSRs) were detected; 3.09% of the unigenes contained SSRs. Large numbers of SSR primer pairs were designed and out of randomly verified primer pairs 80% were successfully yielded amplicons.

Conclusions/Significance

A large of putative R. leptotubula transcript sequences has been obtained from the deep sequencing, which extensively increases the comprehensive and integrated genomic resources of this pest. This large-scale transcriptome dataset will be an important information platform for promoting our investigation of the molecular mechanisms from various aspects in this species.

Molecular cloning, sequence characterization and expression analysis of a CD63 homologue from the coleopteran beetle, Tenebrio molitor

CD63, a member of the tetraspanin membrane protein family, plays a pivotal role in cell growth, motility, signal transduction, host-pathogen interactions and cancer. In this work, the cDNA encoding CD63 homologue (TmCD63) was cloned from larvae of a coleopteran beetle, Tenebrio molitor. The cDNA is comprised of an open reading frame of 705 bp, encoding putative protein of 235 amino acid residues. In silico analysis shows that the protein has four putative transmembrane domains and one large extracellular loop. The characteristic “Cys-Cys-Gly” motif and “Cys188” residues are highly conserved in the large extracellular loop. Phylogenetic analysis of TmCD63 revealed that they belong to the insect cluster with 50%–56% identity. Analysis of spatial expression patterns demonstrated that TmCD63 mRNA is mainly expressed in gut and Malphigian tubules of larvae and the testis of the adult. Developmental expression patterns of CD63 mRNA showed that TmCD63 transcripts are detected in late larval, pupal and adult stages. Interestingly, TmCD63 transcripts are upregulated to the maximum level of 4.5 fold, in response to DAP-type peptidoglycan during the first 6 h, although other immune elicitors also caused significant increase to the transcript level at later time-points. These results suggest that CD63 might contribute to T. molitor immune response against various microbial pathogens.