A De Novo Transcriptomics Approach Reveals Genes Involved in Thrips Tabaci Resistance to Spinosad

Role of CYP9E2 and a long non-coding RNA gene in resistance to a spinosad insecticide in the Colorado potato beetle, Leptinotarsa decemlineata

Spinosads are insecticides used to control insect pests, especially in organic farming where limited tools for pest management exist. However, resistance has developed to spinosads in economically important pests, including Colorado potato beetle (CPB), Leptinotarsa decemlineata. In this study, we used bioassays to determine spinosad sensitivity of two field populations of CPB, one from an organic farm exposed exclusively to spinosad and one from a conventional farm exposed to a variety of insecticides, and a reference insecticide naïve population. We found the field populations exhibited significant levels of resistance compared with the sensitive population. Then, we compared transcriptome profiles between the two field populations to identify genes associated primarily with spinosad resistance and found a cytochrome P450, CYP9E2, and a long non-coding RNA gene, lncRNA-2, were upregulated in the exclusively spinosad-exposed population. Knock-down of these two genes simultaneously in beetles of the spinosad-exposed population using RNA interference (RNAi) resulted in a significant increase in mortality when gene knock-down was followed by spinosad exposure, whereas single knock-downs of each gene produced smaller effects. In addition, knock-down of the lncRNA-2 gene individually resulted in significant reduction in CYP9E2 transcripts. Finally, in silico analysis using an RNA-RNA interaction tool revealed that CYP9E2 mRNA contains multiple binding sites for the lncRNA-2 transcript. Our results imply that CYP9E2 and lncRNA-2 jointly contribute to spinosad resistance in CPB, and lncRNA-2 is involved in regulation of CYP9E2 expression. These results provide evidence that metabolic resistance, driven by overexpression of CYP and lncRNA genes, contributes to spinosad resistance in CPB.

Pest status, molecular evolution, and epigenetic factors derived from the genome assembly of Frankliniella fusca, a thysanopteran phytovirus vector

BACKGROUND

The tobacco thrips (Frankliniella fusca Hinds; family Thripidae; order Thysanoptera) is an important pest that can transmit viruses such as the tomato spotted wilt orthotospovirus to numerous economically important agricultural row crops and vegetables. The structural and functional genomics within the order Thysanoptera has only begun to be explored. Within the > 7000 known thysanopteran species, the melon thrips (Thrips palmi Karny) and the western flower thrips (Frankliniella occidentalis Pergrande) are the only two thysanopteran species with assembled genomes.

RESULTS

A genome of F. fusca was assembled by long-read sequencing of DNA from an inbred line. The final assembly size was 370 Mb with a single copy ortholog completeness of ~ 99% with respect to Insecta. The annotated genome of F. fusca was compared with the genome of its congener, F. occidentalis. Results revealed many instances of lineage-specific differences in gene content. Analyses of sequence divergence between the two Frankliniella species' genomes revealed substitution patterns consistent with positive selection in ~ 5% of the protein-coding genes with 1:1 orthologs. Further, gene content related to its pest status, such as xenobiotic detoxification and response to an ambisense-tripartite RNA virus (orthotospovirus) infection was compared with F. occidentalis. Several F. fusca genes related to virus infection possessed signatures of positive selection. Estimation of CpG depletion, a mutational consequence of DNA methylation, revealed that F. fusca genes that were downregulated and alternatively spliced in response to virus infection were preferentially targeted by DNA methylation. As in many other insects, DNA methylation was enriched in exons in Frankliniella, but gene copies with homology to DNA methyltransferase 3 were numerous and fragmented. This phenomenon seems to be relatively unique to thrips among other insect groups.

CONCLUSIONS

The F. fusca genome assembly provides an important resource for comparative genomic analyses of thysanopterans. This genomic foundation allows for insights into molecular evolution, gene regulation, and loci important to agricultural pest status.

Antennal Transcriptome Analysis of Olfactory Genes and Characterization of Odorant Binding Proteins in Odontothrips loti (Thysanoptera: Thripidae)

To identify odors in complex environments accurately, insects have evolved multiple olfactory proteins. In our study, various olfactory proteins of Odontothrips loti Haliday, an oligophagous pest that primarily affects Medicago sativa (alfalfa), were explored. Specifically, 47 putative olfactory candidate genes were identified in the antennae transcriptome of O. loti, including seven odorant-binding proteins (OBPs), nine chemosensory proteins (CSPs), seven sensory neuron membrane proteins (SNMPs), eight odorant receptors (ORs), and sixteen ionotropic receptors (IRs). PCR analysis further confirmed that 43 out of 47 genes existed in O. loti adults, and O.lotOBP1, O.lotOBP4, and O.lotOBP6 were specifically expressed in the antennae with a male-biased expression pattern. In addition, both the fluorescence competitive binding assay and molecular docking showed that p-Menth-8-en-2-one, a component of the volatiles of the host, had strong binding ability to the O.lotOBP6 protein. Behavioral experiments showed that this component has a significant attraction to both female and male adults, indicating that O.lotOBP6 plays a role in host location. Furthermore, molecular docking reveals potential active sites in O.lotOBP6 that interact with most of the tested volatiles. Our results provide insights into the mechanism of O. loti odor-evoked behavior and the development of a highly specific and sustainable approach for thrip management.

Multiple-P450 Gene Co-Up-Regulation in the Development of Permethrin Resistance in the House Fly, Musca domestica

This paper reports a study conducted at the whole transcriptome level to characterize the P450 genes involved in the development of pyrethroid resistance, utilizing expression profile analyses of 86 cytochrome P450 genes in house fly strains with different levels of resistance to pyrethroids/permethrin. Interactions among the up-regulated P450 genes and possible regulatory factors in different autosomes were examined in house fly lines with different combinations of autosomes from a resistant house fly strain, ALHF. Eleven P450 genes that were significantly up-regulated, with levels > 2-fold those in the resistant ALHF house flies, were in CYP families 4 and 6 and located on autosomes 1, 3 and 5. The expression of these P450 genes was regulated by trans- and/or cis-acting factors, especially on autosomes 1 and 2. An in vivo functional study indicated that the up-regulated P450 genes also conferred permethrin resistance in Drosophila melanogaster transgenic lines. An in vitro functional study confirmed that the up-regulated P450 genes are able to metabolize not only cis- and trans-permethrin, but also two metabolites of permethrin, PBalc and PBald. In silico homology modeling and the molecular docking methodology further support the metabolic capacity of these P450s for permethrin and substrates. Taken together, the findings of this study highlight the important function of multi-up-regulated P450 genes in the development of insecticide resistance in house flies.

Groundnut Bud Necrosis Virus Modulates the Expression of Innate Immune, Endocytosis, and Cuticle Development-Associated Genes to Circulate and Propagate in Its Vector, Thrips palmi

Thrips palmi (Thysanoptera: Thripidae) is the predominant tospovirus vector in Asia-Pacific region. It transmits economically damaging groundnut bud necrosis virus (GBNV, family Tospoviridae) in a persistent propagative manner. Thrips serve as the alternate host, and virus reservoirs making tospovirus management very challenging. Insecticides and host plant resistance remain ineffective in managing thrips–tospoviruses. Recent genomic approaches have led to understanding the molecular interactions of thrips–tospoviruses and identifying novel genetic targets. However, most of the studies are limited to Frankliniella species and tomato spotted wilt virus (TSWV). Amidst the limited information available on T. palmi–tospovirus relationships, the present study is the first report of the transcriptome-wide response of T. palmi associated with GBNV infection. The differential expression analyses of the triplicate transcriptome of viruliferous vs. nonviruliferous adult T. palmi identified a total of 2,363 (1,383 upregulated and 980 downregulated) significant transcripts. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses showed the abundance of differentially expressed genes (DEGs) involved in innate immune response, endocytosis, cuticle development, and receptor binding and signaling that mediate the virus invasion and multiplication in the vector system. Also, the gene regulatory network (GRN) of most significant DEGs showed the genes like ABC transporter, cytochrome P450, endocuticle structural glycoprotein, gamma-aminobutyric acid (GABA) receptor, heat shock protein 70, larval and pupal cuticle proteins, nephrin, proline-rich protein, sperm-associated antigen, UHRF1-binding protein, serpin, tyrosine–protein kinase receptor, etc., were enriched with higher degrees of interactions. Further, the expression of the candidate genes in response to GBNV infection was validated in reverse transcriptase-quantitative real-time PCR (RT-qPCR). This study leads to an understanding of molecular interactions between T. palmi and GBNV and suggests potential genetic targets for generic pest control.

A simple guide to de novo transcriptome assembly and annotation

A transcriptome constructed from short-read RNA sequencing (RNA-seq) is an easily attainable proxy catalog of protein-coding genes when genome assembly is unnecessary, expensive or difficult. In the absence of a sequenced genome to guide the reconstruction process, the transcriptome must be assembled de novo using only the information available in the RNA-seq reads. Subsequently, the sequences must be annotated in order to identify sequence-intrinsic and evolutionary features in them (for example, protein-coding regions). Although straightforward at first glance, de novo transcriptome assembly and annotation can quickly prove to be challenging undertakings. In addition to familiarizing themselves with the conceptual and technical intricacies of the tasks at hand and the numerous pre- and post-processing steps involved, those interested must also grapple with an overwhelmingly large choice of tools. The lack of standardized workflows, fast pace of development of new tools and techniques and paucity of authoritative literature have served to exacerbate the difficulty of the task even further. Here, we present a comprehensive overview of de novo transcriptome assembly and annotation. We discuss the procedures involved, including pre- and post-processing steps, and present a compendium of corresponding tools.