Transcriptome Comparison Analysis of Ostrinia furnacalis in Four Developmental Stages

Analysis of Genes Associated with Feeding Preference and Detoxification in Various Developmental Stages of Aglais urticae

Herbivorous insects and host plants have developed a close and complex relationship over a long period of co-evolution. Some plants provide nutrients for insects, but plants' secondary metabolites also influence their growth and development. Urtica cannabina roots and leaves are poisonous, yet Aglais urticae larvae feed on them, so we aimed to clarify the mechanism enabling this interaction. At present, studies on the detoxification mechanism of the A. urticae are rare. In our study, first, we used the A. urticae larval odor selection behavior bioassay and choice feeding preference assay to analyze the feeding preferences of A. urticae on its host plant, U. cannabina. Next, we used transcriptome sequencing to obtain the unigenes annotated and classified by various databases, such as KEGG and GO. In this study, we found that U. cannabina could attract A. urticae larvae to feed via scent, and the feeding preference assay confirmed that larvae preferred U. cannabina leaves over three other plants: Cirsium japonicum, Cannabis sativa, and Arctium lappa. The activity of detoxifying enzymes GST and CarE changed in larvae that had consumed U. cannabina. Furthermore, through transcriptomic sequencing analysis, 77,624 unigenes were assembled from raw reads. The numbers of differentially expressed genes were calculated using pairwise comparisons of all life stages; the expression of detoxification enzyme genes was substantially higher in larvae than in the pupal and adult stages. Finally, we identified and summarized 34 genes associated with detoxification enzymes, such as UDP-glucose 4-epimerase gene, 5 Glutathione S-transferase genes, 4 Carboxylesterase genes, 4 Cytochrome P450 genes, 10 ATP-binding cassette genes, 4 Superoxide dismutase, and Peroxidase. Moreover, we identified 28 genes associated with the development of A. urticae. The qRT-PCR results were nearly consistent with the transcriptomic data, showing an increased expression level of four genes in larvae. Taken together, this study examines the correlation between A. urticae and host plants U. cannabina, uncovering a pronounced preference for A. urticae larvae toward host plants. Consistent with RNA-seq, we investigated the mechanism of A. urticae's interaction with host plants and identified detoxification-related genes. The present study provides theoretical support for studying insect adaptation mechanisms and biological control.

Transcriptome profiling for developmental stages Protaetia brevitarsis seulensis with focus on wing development and metamorphosis

A white-spotted flower chafer Protaetia brevitarsis seulensis widely distributed in Asian countries is traditionally used in oriental medicine. This study explored gene expression abundance with respect to wing development and metamorphosis in P. b. seulensis based on the large-scale RNA-seq data. The transcriptome assembly consists of 23,551 high-quality transcripts which are approximately 96.7% covered. We found 265 wing development genes, 19 metamorphosis genes, and 1,314 candidates. Of the 1,598 genes, 1,594 are included exclusively in cluster 4 with similar gene co-expression patterns. The network centrality analyses showed that wing development- and metamorphosis-related genes have a high degree of betweenness centrality and are expressed most highly in eggs, moderately in pupa and adults, and lowest in larva. This study provides some meaningful clues for elucidating the genetic modulation mechanism of wing development and metamorphosis in P. b. seulensis.

Recapitulation of the embryonic transcriptional program in holometabolous insect pupae

Holometabolous insects are predominantly motionless during metamorphosis, when no active feeding is observed and the body is enclosed in a hardened cuticle. These physiological properties as well as undergoing processes resemble embryogenesis, since at the pupal stage organs and systems of the imago are formed. Therefore, recapitulation of the embryonic expression program during metamorphosis could be hypothesized. To assess this hypothesis at the transcriptome level, we have performed a comprehensive analysis of the developmental datasets available in the public domain. Indeed, for most datasets, the pupal gene expression resembles the embryonic rather than the larval pattern, interrupting gradual changes in the transcriptome. Moreover, changes in the transcriptome profile during the pupa-to-imago transition are positively correlated with those at the embryo-to-larvae transition, suggesting that similar expression programs are activated. Gene sets that change their expression level during the larval stage and revert it to the embryonic-like state during the metamorphosis are enriched with genes associated with metabolism and development.

Ostrinia furnacalis PBP2 solution NMR structure: Insight into ligand binding and release mechanisms

Ostrinia furnacalis is an invasive lepidopteran agricultural pest that relies on olfaction for mating and reproduction. Male moths have an extremely sensitive olfactory system that can detect the sex pheromones emitted by females over a great distance. Pheromone-binding proteins present in the male moth antenna play a key role in the pheromone uptake, transport, and release at the dendritic membrane of the olfactory neuron. Here, we report the first high-resolution NMR structure of a pheromone-binding protein from an Ostrinia species at pH 6.5. The core of the Ostrinia furnacalis PBP2 (OfurPBP2) consists of six helices, α1a (2-14), α1b (16-22), α2 (27-37), α3 (46-60), α4 (70-80), α5 (84-100), and α6 (107-124) surrounding a large hydrophobic pocket. The structure is stabilized by three disulfide bridges, 19-54, 50-108, and 97-117. In contrast to the unstructured C-terminus of other lepidopteran PBPs, the C-terminus of OfurPBP2 folds into an α-helix (α7) at pH 6.5. The protein has nanomolar affinity towards both pheromone isomers. Molecular docking of both pheromones, E-12 and Z-12-tetradecenyl acetate, to OfurPBP2 revealed that the residues Met5, Lys6, Met8, Thr9, Phe12, Phe36, Trp37, Phe76, Ser115, Phe118, Lys119, Ile122, His123, and Ala128 interact with both isomers, while Thr9 formed a hydrogen bond with the acetate head group. NMR structure and thermal unfolding studies with CD suggest that ligand release at pH 4.5 is likely due to the partial unfolding of the protein.

Antennal Transcriptome of the Fruit-Sucking Moth Eudocima materna: Identification of Olfactory Genes and Preliminary Evidence for RNA-Editing Events in Odorant Receptors

Unappealing shriveled fruits are a characteristic of one of the most elusive fruit pests. The perpetrator, Eudocima materna, attacks the fruit at a fully formed stage and, therefore, the antennal transcriptome for this insect was deduced to identify the molecular elicitors involved in the attraction to its host plants. A total of 260 olfactory genes, including 16 odorant-binding proteins (OBPs), four pheromone-binding proteins (PBPs), 40 antennal-binding proteins (ABPs), 178 odorant receptors (ORs), 17 chemosensory proteins (CSPs) and five sensory neuron membrane proteins (SNMPs) were identified. Phylogenetic analysis shows the divergence of E. materna proteins from closely related lepidopterans and provides insights on genes that have exclusively evolved in this insect. STRING network analysis revealed interactions of olfactory proteins among themselves and the proteins of other groups. Interestingly, online tools predicted RNA-editing events in the odorant receptor sequences, suggesting the possibility of multiple protein forms. Transcripts matching transposable element sequences were also detected in the dataset. Thus, the work reported here provides a valuable resource to design molecular methods for pest control.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

Background

Methods

Results

Conclusion

Transcriptome and Proteome Alternation With Resistance to Bacillus thuringiensis Cry1Ah Toxin in Ostrinia furnacalis

Background: Asian corn borer (ACB), Ostrinia furnacalis can develop resistance to transgenic Bacillus thuringiensis (Bt) maize expressing Cry1Ah-toxin. However, the mechanisms that regulate the resistance of ACB to Cry1Ah-toxin are unknown.

Objective: In order to understand the molecular basis of the Cry1Ah-toxin resistance in ACB, “omics” analyses were performed to examine the difference between Cry1Ah-resistant (ACB-AhR) and susceptible (ACB-BtS) strains of ACB at both transcriptional and translational levels.

Results: A total of 7,007 differentially expressed genes (DEGs) and 182 differentially expressed proteins (DEPs) were identified between ACB-AhR and ACB-BtS and 90 genes had simultaneous transcription and translation profiles. Down-regulated genes associated with Cry1Ah resistance included aminopeptidase N, ABCC3, DIMBOA-induced cytochrome P450, alkaline phosphatase, glutathione S-transferase, cadherin-like protein, and V-ATPase. Whereas, anti-stress genes, such as heat shock protein 70 and carboxylesterase were up-regulated in ACB-AhR, displaying that a higher proportion of genes/proteins related to resistance was down-regulated compared to up-regulated. The Kyoto encyclopedia of genes and genomes (KEGG) analysis mapped 578 and 29 DEGs and DEPs, to 27 and 10 pathways, respectively (P < 0.05). Furthermore, real-time quantitative (qRT-PCR) results based on relative expression levels of randomly selected genes confirmed the “omics” response.

Conclusion: Despite the previous studies, this is the first combination of a study using RNA-Seq and iTRAQ approaches on Cry1Ah-toxin binding, which led to the identification of longer length of unigenes in ACB. The DEGs and DEPs results are valuable for further clarifying Cry1Ah-mediated resistance.

The Asian corn borer Ostrinia furnacalis feeding increases the direct and indirect defence of mid-whorl stage commercial maize in the field

The Asian corn borer (Ostrinia furnacalis Guenée) is a destructive pest of maize (Zea mays L.). Despite large-scale commercial maize production, little is known about the defensive responses of field-grown commercial maize to O. furnacalis herbivory, and how these responses result in direct and indirect defence against this pest. To elucidate the maize transcriptome response to O. furnacalis feeding, leaves of maize hybrid Jingke968 were infested with O. furnacalis for 0, 2, 4, 12 and 24 h. Ostrinia furnacalis feeding elicited stronger and more rapid changes in the defence-related gene expression (i.e. after 2 h), and more differentially expressed genes (DEGs) were up-regulated than down-regulated at all times post-induction (i.e. 2, 4, 12 and 24 h) in the O. furnacalis pre-infested maize plants. KEGG pathway analysis indicated that the DEGs in the O. furnacalis pre-infested maize are involved in benzoxazinoids, phytohormones, volatiles, and other metabolic pathways related to maize resistance to herbivores. In addition, the maize leaves previously infested by O. furnacalis for 24 h showed an obvious inhibition of the subsequent O. furnacalis performance, and maize volatiles induced by O. furnacalis feeding for 24 and 48 h attracted the parasitic wasp, Macrocentrus cingulum Brischke. The increased direct and indirect defences induced by O. furnacalis feeding were correlated with O. furnacalis-induced phytohormones, benzoxazinoids, and volatiles. Together, our findings provide new insights into how commercial maize orchestrates its transcriptome and metabolome to directly and indirectly defend against O. furnacalis at the mid-whorl stage in the field.

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.

De Novo Assembly and Analysis of the White-Backed Planthopper (Sogatella furcifera) Transcriptome

The white-backed planthopper Sogatella furcifera (Horváth) has become an important pest on rice in China and Southeast Asian countries. White-backed planthopper wing bud length is in relation to adult wing length, but little is known about the development and differentiation of wing buds at the molecular level. Using Illumina HiSeq high-throughput sequencing technology, we sequenced four cDNA libraries, two biological replicates of long-winged female fifth-instar nymphs (LW), and two of short-winged nymphs (SW). In total, 62,154 unigenes with an average length of 984 bp and N50 length of 1,878 bp were obtained by de novo transcriptome assembly. A total of 18,416 open reading frames (ORFs) were predicted based on the unigenes. Ninety-three percentage of these ORFs could be annotated by searching for homology in six protein databases. A total of 184 differentially expressed genes (DEGs) with 129 upregulated and 55 downregulated were found in SW compared to LW. Gene Ontology and euKaryotic Orthologous Group classification provided comprehensive information about the function of each gene. Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed five enriched pathways including three metabolic pathways. In addition, we found that some DEGs were relevant to muscle movement and cuticle and likely involved in development and differentiation of wing buds. This study provided transcriptome resource of female fifth-instar nymphs of white-backed planthopper including long-winged and short-winged nymphs, and different molecular features between them lay the foundation for adult wing morph prediction, promoting further studies on planthopper population management.

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.

Transcriptomic analysis of gene expression in mice treated with troxerutin

Troxerutin, a semi-synthetic derivative of the natural bioflavanoid rutin, has been reported to possess many beneficial effects in human bodies, such as vasoprotection, immune support, anti-inflammation and anti-aging. However, the effects of troxerutin on genome-wide transcription in blood cells are still unknown. In order to find out effects of troxerutin on gene transcription, a high-throughput RNA sequencing was employed to analysis differential gene expression in blood cells consisting of leucocytes, erythrocytes and platelets isolated from the mice received subcutaneous injection of troxerutin. Transcriptome analysis demonstrated that the expression of only fifteen genes was significantly changed by the treatment with troxerutin, among which 5 genes were up-regulated and 10 genes were down-regulated. Bioinformatic analysis of the fifteen differentially expressed genes was made by utilizing the Gene Ontology (GO), and the differential expression induced by troxerutin was further evaluated by real-time quantitative PCR (Q-PCR).

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

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