Generation of a Transcriptome in a Model Lepidopteran Pest, Heliothis virescens, Using Multiple Sequencing Strategies for Profiling Midgut Gene Expression

Hemocyte response to treatment of susceptible and resistant Asian corn borer (Ostrinia furnacalis) larvae with Cry1F toxin from Bacillus thuringiensis

Midgut receptors have been recognized as the major mechanism of resistance to Cry proteins in lepidopteran larvae, while there is a dearth of data on the role of hemocyte's response to Cry intoxication and resistance development. We aimed at investigating the role of circulating hemocytes in the intoxication of Cry1F toxin in larvae from susceptible (ACB-BtS) and resistant (ACB-FR) strains of the Asian corn borer (ACB), Ostrinia furnacalis. Transcriptome and proteome profiling identified genes and proteins involved in immune-related (tetraspanin and C-type lectins) and detoxification pathways as significantly up-regulated in the hemocytes of Cry1F treated ACB-FR. High-throughput in vitro assays revealed the binding affinity of Cry1F with the tetraspanin and C-type lectin family proteins. We found significant activation of MAPKinase (ERK 1/2, p38α, and JNK 1/2) in the hemocytes of Cry1F treated ACB-FR. In testing plausible crosstalk between a tetraspanin (CD63) and downstream MAPK signaling, we knocked down CD63 expression by RNAi and detected an alteration in JNK 1/2 level but a significant increase in susceptibility of ACB-FR larvae to Cry1F toxin. Information from this study advances a change in knowledge on the cellular immune response to Cry intoxication and its potential role in resistance in a lepidopteran pest.

Comparative gut transcriptome analysis of Diatraea saccharalis in response to the dietary source

The sugarcane borer (Diatraea saccharalis, Fabricius, 1794) is a devastating pest that causes millions of dollars of losses each year to sugarcane producers by reducing sugar and ethanol yields. The control of this pest is difficult due to its endophytic behavior and rapid development. Pest management through biotechnological approaches has emerged in recent years as an alternative to currently applied methods. Genetic information about the target pests is often required to perform biotechnology-based management. The genomic and transcriptomic data for D. saccharalis are very limited. Herein, we report a tissue-specific transcriptome of D. saccharalis larvae and a differential expression analysis highlighting the physiological characteristics of this pest in response to two different diets: sugarcane and an artificial diet. Sequencing was performed on the Illumina HiSeq 2000 platform, and a de novo assembly was generated. A total of 27,626 protein-coding unigenes were identified, among which 1,934 sequences were differentially expressed between treatments. Processes such as defence, digestion, detoxification, signaling, and transport were highly represented among the differentially expressed genes (DEGs). Furthermore, seven aminopeptidase genes were identified as candidates to encode receptors of Cry proteins, which are toxins of Bacillus thuringiensis used to control lepidopteran pests. Since plant-insect interactions have produced a considerable number of adaptive responses in hosts and herbivorous insects, the success of phytophagous insects relies on their ability to overcome challenges such as the response to plant defences and the intake of nutrients. In this study, we identified metabolic pathways and specific genes involved in these processes. Thus, our data strongly contribute to the knowledge advancement of insect transcripts, which can be a source of target genes for pest management.

Transcriptome analysis of the almond moth, Cadra cautella, female abdominal tissues and identification of reproduction control genes

BACKGROUND

The almond moth, Cadra cautella is a destructive pest of stored food commodities including dates that causes severe economic losses for the farming community worldwide. To date, no genetic information related to the molecular mechanism/strategies of its reproduction is available. Thus, transcriptome analysis of C. cautella female abdominal tissues was performed via next-generation sequencing (NGS) to recognize the genes responsible for reproduction.

RESULTS

The NGS was performed with an Illumina Hiseq 2000 sequencer (Beijing Genomics Institute: BGI). From the transcriptome data, 9,804,804,120 nucleotides were generated and their assemblage resulted in 62,687 unigenes. The functional annotation analyses done by different databases, annotated, 27,836 unigenes in total. The transcriptome data of C. cautella female abdominal tissue was submitted to the National Center for Biotechnology Information (accession no: PRJNA484692). The transcriptome analysis yielded several genes responsible for C. cautella reproduction including six Vg gene transcripts. Among the six Vg gene transcripts, only one was highly expressed with 3234.95 FPKM value (fragments per kilobase per million mapped reads) that was much higher than that of the other five transcripts. Higher differences in the expression level of the six Vg transcripts were confirmed by running the RT-PCR using gene specific primers, where the expression was observed only in one transcript it was named as the CcVg.

CONCLUSIONS

This is the first study to explore C. cautella reproduction control genes and it might be supportive to explore the reproduction mechanism in this pest at the molecular level. The NGS based transcriptome pool is valuable to study the functional genomics and will support to design biotech-based management strategies for C. cautella.

Transcriptome and differentially expressed genes of Busseola fusca (Lepidoptera: Noctuidae) larvae challenged with Cry1Ab toxin

Busseola fusca (Fuller) (Lepidoptera: Noctuidae), a major insect pest of maize in sub-Saharan Africa, has developed high levels of non-recessive resistance to Cry1Ab toxin expressed in genetically modified Bt maize. Multiple resistance mechanisms to various Cry toxins have been identified in Lepidoptera, but no study has yet been done to determine the mechanism of Cry1Ab resistance in B. fusca. Therefore, the larval transcriptome of B. fusca was sequenced, de novo assembled and characterized. Differential expression analysis was performed to compare gene expression profiles of Cry toxin challenged and unchallenged neonate larvae to assess the molecular basis of the defence mechanism employed by this insect. Several genes associated with Cry toxin resistance in other lepidopteran pests were detected in B. fusca. Results suggest that differential expression of metabolic and immune-related genes might explain Cry1Ab toxin defence in this pest (supplemental file). Transcript expression profiles of neonates demonstrated that 33.59% and 60.31% of the 131 differentially expressed genes were upregulated and downregulated in the toxin-challenged neonate larvae, respectively. Transcripts were grouped into two subclusters according to the similarity of their expression patterns. Transcripts in subcluster 1 were moderately upregulated in the toxin-challenged neonate larvae, and, conversely, downregulated in the unchallenged neonate larvae. The solute carrier organic anion transporter, which is involved in insecticide detoxification, was upregulated in the toxin-challenged neonate larvae. Conversely, most of the transcripts in subcluster 2 were moderately downregulated in the toxin-challenged neonate larvae, and upregulated for neonates feeding on non-challenged maize. Four unidentified transcripts were extremely down-regulated in the toxin-challenged neonate larvae, and upregulated in the unchallenged neonate larvae. Further studies are recommended to establish if there is a direct correlation between these differentially expressed genes and the observed resistance. Elucidation of such defence mechanisms is crucial for developing insect resistance management strategies to ensure sustainable use of genetically modified maize in Africa. Nevertheless, this is the first study on gene expression profiles of B. fusca strains challenged with Cry toxin. The transcriptome characterized in this study provides a significant resource base for future studies on B. fusca and contributes to understanding some of the gene regulation and signalling networks involved in the defence of B. fusca against Cry1Ab toxin.

HT-SuperSAGE of the gut tissue of a Vip3Aa-resistant Heliothis virescens (Lepidoptera: Noctuidae) strain provides insights into the basis of resistance

Multitoxin Bt-crops expressing insecticidal toxins with different modes of action, for example, Cry and Vip, are expected to improve resistance management in target pests. While Cry1A resistance has been relatively well characterized in some insect species, this is not the case for Vip3A, for which no mechanism of resistance has yet been identified. Here we applied HT-SuperSAGE to analyze the transcriptome of the gut tissue of tobacco budworm Heliothis virescens (F.) laboratory-selected for Vip3Aa resistance. From a total of 1 324 252 sequence reads, 5 895 126-bp tags were obtained representing 17 751 nonsingleton unique transcripts (UniTags) from genetically similar Vip3Aa-resistant (Vip-Sel) and susceptible control (Vip-Unsel) strains. Differential expression was significant (≥2.5 fold or ≤0.4; P < 0.05) for 1989 sequences (11.2% of total UniTags), where 420 represented overexpressed (OE) and 1569 underexpressed (UE) genes in Vip-Sel. BLASTN searches mapped 419 UniTags to H. virescens sequence contigs, of which, 416 (106 OE and 310 UE) were unambiguously annotated to proteins in NCBI nonredundant protein databases. Gene Ontology distributed 345 of annotated UniTags in 14 functional categories with metabolism (including serine-type hydrolases) and translation/ribosome biogenesis being the most prevalent. A UniTag homologous to a particular member of the REsponse to PAThogen (REPAT) family was found among most overexpressed, while UniTags related to the putative Vip3Aa-binding ribosomal protein S2 (RpS2) were underexpressed. qRT-PCR of a subset of UniTags validated the HT-SuperSAGE data. This study is the first providing lepidopteran gut transcriptome associated with Vip3Aa resistance and a foundation for future attempts to elucidate the resistance mechanism.

Transcriptomic response of female adult moths to host and non-host plants in two closely related species

BACKGROUND

Divergent selection has been shown to promote speciation in many taxa and especially in phytophagous insects. In the Ostrinia species complex, the European corn borer (ECB) and adzuki bean borer (ABB) are two sibling species specialized to different host plants. The first is a well-known maize pest, whereas the second is a polyphagous species associated with various dicotyledons. Their specialization to host plants is driven by morphological, behavioral and physiological adaptations. In particular, previous studies have shown that ECB and ABB display marked behavior with regard to plant choice during oviposition, involving specific preference and avoidance mechanisms. In this study, our goal was to identify the mechanisms underlying this host-plant specialization in adult females through an analysis of their gene expression. We assembled and annotated a de novo reference transcriptome and measured differences in gene expression between ECB and ABB females, and between environments. We related differentially expressed genes to host preference behavior, and highlighted the functional categories involved. We also conducted a specific analysis of chemosensory genes, which are considered to be good candidates for host recognition before oviposition.

RESULTS

We recorded more differentially expressed genes in ECB than in ABB samples, and noticed that the majority of genes potentially involved in the host preference were different between the two species. At the functional level, the response to plant environment in adult females involved many processes, including the metabolism of carbohydrates, lipids, proteins, and amino acids; detoxification mechanisms and immunity; and the chemosensory repertoire (as expected). Until now, most of the olfactory receptors described in Ostrinia spp. had been tested for their putative role in pheromone recognition by males. Here we observed that one specific olfactory receptor was clearly associated with ECB's discrimination between maize and mugwort conditions, highlighting a potential new candidate involved in plant odor discrimination in adult females.

CONCLUSIONS

Our results are a first step toward the identification of candidate genes and functions involved in chemosensory processes, carbohydrate metabolism, and virus and retrovirus dynamics. These candidates provide new avenues for research into understanding the role of divergent selection between different environments in species diversification.

Nucleases as a barrier to gene silencing in the cotton boll weevil, Anthonomus grandis

RNA interference (RNAi) approaches have been applied as a biotechnological tool for controlling plant insect pests via selective gene down regulation. However, the inefficiency of RNAi mechanism in insects is associated with several barriers, including dsRNA delivery and uptake by the cell, dsRNA interaction with the cellular membrane receptor and dsRNA exposure to insect gut nucleases during feeding. The cotton boll weevil (Anthonomus grandis) is a coleopteran in which RNAi-mediated gene silencing does not function efficiently through dsRNA feeding, and the factors involved in the mechanism remain unknown. Herein, we identified three nucleases in the cotton boll weevil transcriptome denoted AgraNuc1, AgraNuc2, and AgraNuc3, and the influences of these nucleases on the gene silencing of A. grandis chitin synthase II (AgraChSII) were evaluated through oral dsRNA feeding trials. A phylogenetic analysis showed that all three nucleases share high similarity with the DNA/RNA non-specific endonuclease family of other insects. These nucleases were found to be mainly expressed in the posterior midgut region of the insect. Two days after nuclease RNAi-mediated gene silencing, dsRNA degradation by the gut juice was substantially reduced. Notably, after nucleases gene silencing, the orally delivered dsRNA against the AgraChSII gene resulted in improved gene silencing efficiency when compared to the control (non-silenced nucleases). The data presented here demonstrates that A. grandis midgut nucleases are effectively one of the main barriers to dsRNA delivery and emphasize the need to develop novel RNAi delivery strategies focusing on protecting the dsRNA from gut nucleases and enhancing its oral delivery and uptake to crop insect pests.

Transcript expression plasticity as a response to alternative larval host plants in the speciation process of corn and rice strains of Spodoptera frugiperda

BACKGROUND

Our main purpose was to evaluate the expression of plastic and evolved genes involved in ecological speciation in the noctuid moth Spodoptera frugiperda, the fall armyworm (FAW); and to demonstrate how host plants might influence lineage differentiation in this polyphagous insect. FAW is an important pest of several crops worldwide, and it is differentiated into host plant-related strains, corn (CS) and rice strains (RS). RNA-Seq and transcriptome characterization were applied to evaluate unbiased genetic expression differences in larvae from the two strains, fed on primary (corn) and alternative (rice) host plants. We consider that genes that are differently regulated by the same FAW strain, as a response to different hosts, are "plastic". Otherwise, differences in gene expression between the two strains fed on the same host are considered constitutive differences.

RESULTS

Individual performance parameters (larval and pupal weight) varied among conditions (strains vs. hosts). A total of 3657 contigs was related to plastic response, and 2395 contigs were differentially regulated in the two strains feeding on preferential and alternative hosts (constitutive contigs). Three molecular functions were present in all comparisons, both down- and up-regulated: oxidoreductase activity, metal-ion binding, and hydrolase activity.

CONCLUSIONS

Metabolization of foreign chemicals is among the key functions involved in the phenotypic variation of FAW strains. From an agricultural perspective, high plasticity in families of detoxifying genes indicates the capacity for a rapid response to control compounds such as insecticides.

Alpha-arylphorin is a mitogen in the Heliothis virescens midgut cell secretome upon Cry1Ac intoxication

Insecticidal crystal (Cry) proteins produced by the bacterium Bacillus thuringiensis (Bt) target cells in the midgut epithelium of susceptible larvae. While the mode of action of Cry toxins has been extensively investigated, the midgut response to Cry intoxication and its regulation are not well characterized. In this work, we describe the secreted proteome (secretome) of primary mature midgut cell cultures from Heliothis virescens larvae after exposure to Cry1Ac toxin compared to control buffer treatment. The Cry1Ac-induced secretome caused higher proliferation and differentiation and an overall reduction in total cell mortality over time in primary H. virescens midgut stem cell cultures when compared to treatment with control buffer secretome. Differential proteomics identified four proteins with significant differences in abundance comparing Cry1Ac-treated and control secretomes. The most significant difference detected in the Cry1Ac secretome was an arylphorin subunit alpha protein not detected in the control secretome. Feeding of purified alpha-arylphorin to H. virescens larvae resulted in midgut hyperplasia and significantly reduced susceptibility to Cry1Ac toxin compared to controls. These data identify alpha-arylphorin as a protein with a new putative role in the midgut regeneration process in response to Cry1Ac intoxication and possibly pathogen/abiotic stress, identifying alpha-arylphorin as a potential gene to target with insecticidal gene silencing for pest control.

A transcriptomic survey of Migdolus fryanus (sugarcane rhizome borer) larvae

Sugarcane, a major crop grown in the tropical and subtropical areas of the world, is produced mainly for sucrose, which is used as a sweetener or for the production of bioethanol. Among the numerous pests that significantly affect the yield of sugarcane, the sugarcane rhizome borer (Migdolus fryanus, a cerambycidae beetle) is known to cause severe damage to the crops in Brazil. The absence of molecular information about this insect reinforces the need for studies and an effective method to control this pest. In this study, RNA-Seq technology was employed to study different parts of M. fryanus larvae. The generated data will help in further investigations about the taxonomy, development, and adaptation of this insect. RNA was extracted from six different parts (head, fat body, integument, hindgut, midgut, and foregut) using Trizol methodology. Using Illumina paired-end sequencing technology and the Trinity platform, trimming and de novo assembly was performed, resulting in 44,567 contigs longer than 200 nt for a reunion of data from all transcriptomes, with a mean length of 1,095.27 nt. Transcripts were annotated using BLAST against different protein databanks (Uniprot/Swissprot, PFAM, KEEG, SignalP 4.1, Gene Ontology, and CAZY) and were compared for similarity using a Venn diagram. Differential expression patterns were studied for select genes through qPCR and FPKM comprising important protein families (digestive peptidases, glucosyl hydrolases, serine protease inhibitors and otopetrin), which allowed a better understanding of the insect’s digestion, immunity and gravity sensorial mechanisms.

Exploring developmental gene toolkit and associated pathways in a potential new model crustacean using transcriptomic analysis

The crustaceans are one of the largest, most diverse, and most successful groups of invertebrates. The diversity among the crustaceans is also reflected in embryonic development models. However, the molecular genetics that regulates embryonic development is not known in those crustaceans that have a short germ-band development with superficial cleavage, such as Macrobrachium olfersi. This species is a freshwater decapod and has great potential to become a model for developmental biology, as well as for evolutionary and environmental studies. To obtain sequence data of M. olfersi from an embryonic developmental perspective, we performed de novo assembly and annotation of the embryonic transcriptome. Using a pooling strategy of total RNA, paired-end Illumina sequencing, and assembly with multiple k-mers, a total of 25,636,097 pair reads were generated. In total, 99,751 unigenes were identified, and 20,893 of these returned a Blastx hit. KEGG pathway analysis mapped a total of 6866 unigenes related to 129 metabolic pathways. In general, 21,845 unigenes were assigned to gene ontology (GO) categories: molecular function (19,604), cellular components (10,254), and biological processes (13,841). Of these, 2142 unigenes were assigned to the developmental process category. More specifically, a total of 35 homologs of embryonic development toolkit genes were identified, which included maternal effect (one gene), gap (six), pair-rule (six), segment polarity (seven), Hox (four), Wnt (eight), and dorsoventral patterning genes (three). In addition, genes of developmental pathways were found, including TGF-β, Wnt, Notch, MAPK, Hedgehog, Jak-STAT, VEGF, and ecdysteroid-inducible nuclear receptors. RT-PCR analysis of eight genes related to embryonic development from gastrulation to late morphogenesis/organogenesis confirmed the applicability of the transcriptome analysis.

Sequencing, de novo assembly and annotation of a pink bollworm larval midgut transcriptome

BACKGROUND

The pink bollworm Pectinophora gossypiella (Saunders) (Lepidoptera: Gelechiidae) is one of the world's most important pests of cotton. Insecticide sprays and transgenic cotton producing toxins of the bacterium Bacillus thuringiensis (Bt) are currently used to manage this pest. Bt toxins kill susceptible insects by specifically binding to and destroying midgut cells, but they are not toxic to most other organisms. Pink bollworm is useful as a model for understanding insect responses to Bt toxins, yet advances in understanding at the molecular level have been limited because basic genomic information is lacking for this cosmopolitan pest. Here, we have sequenced, de novo assembled and annotated a comprehensive larval midgut transcriptome from a susceptible strain of pink bollworm.

FINDINGS

A de novo transcriptome assembly for the midgut of P. gossypiella was generated containing 46,458 transcripts (average length of 770 bp) derived from 39,874 unigenes. The size of the transcriptome is similar to published midgut transcriptomes of other Lepidoptera and includes up to 91 % annotated contigs. The dataset is publicly available in NCBI and GigaDB as a resource for researchers.

CONCLUSIONS

Foundational knowledge of protein-coding genes from the pink bollworm midgut is critical for understanding how this important insect pest functions. The transcriptome data presented here represent the first large-scale molecular resource for this species, and may be used for deciphering relevant midgut proteins critical for xenobiotic detoxification, nutrient digestion and allocation, as well as for the discovery of protein receptors important for Bt intoxication.