Analysis of expressed sequence tags for Frankliniella occidentalis, the western flower thrips

Genome-enabled insights into the biology of thrips as crop pests

BACKGROUND

The western flower thrips, Frankliniella occidentalis (Pergande), is a globally invasive pest and plant virus vector on a wide array of food, fiber, and ornamental crops. The underlying genetic mechanisms of the processes governing thrips pest and vector biology, feeding behaviors, ecology, and insecticide resistance are largely unknown. To address this gap, we present the F. occidentalis draft genome assembly and official gene set.

RESULTS

We report on the first genome sequence for any member of the insect order Thysanoptera. Benchmarking Universal Single-Copy Ortholog (BUSCO) assessments of the genome assembly (size = 415.8 Mb, scaffold N50 = 948.9 kb) revealed a relatively complete and well-annotated assembly in comparison to other insect genomes. The genome is unusually GC-rich (50%) compared to other insect genomes to date. The official gene set (OGS v1.0) contains 16,859 genes, of which ~ 10% were manually verified and corrected by our consortium. We focused on manual annotation, phylogenetic, and expression evidence analyses for gene sets centered on primary themes in the life histories and activities of plant-colonizing insects. Highlights include the following: (1) divergent clades and large expansions in genes associated with environmental sensing (chemosensory receptors) and detoxification (CYP4, CYP6, and CCE enzymes) of substances encountered in agricultural environments; (2) a comprehensive set of salivary gland genes supported by enriched expression; (3) apparent absence of members of the IMD innate immune defense pathway; and (4) developmental- and sex-specific expression analyses of genes associated with progression from larvae to adulthood through neometaboly, a distinct form of maturation differing from either incomplete or complete metamorphosis in the Insecta.

CONCLUSIONS

Analysis of the F. occidentalis genome offers insights into the polyphagous behavior of this insect pest that finds, colonizes, and survives on a widely diverse array of plants. The genomic resources presented here enable a more complete analysis of insect evolution and biology, providing a missing taxon for contemporary insect genomics-based analyses. Our study also offers a genomic benchmark for molecular and evolutionary investigations of other Thysanoptera species.

Discovery of Novel Thrips Vector Proteins That Bind to the Viral Attachment Protein of the Plant Bunyavirus Tomato Spotted Wilt Virus

The plant-pathogenic virus tomato spotted wilt virus (TSWV) encodes a structural glycoprotein (GN) that, like with other bunyavirus/vector interactions, serves a role in viral attachment and possibly in entry into arthropod vector host cells. It is well documented that Frankliniella occidentalis is one of nine competent thrips vectors of TSWV transmission to plant hosts. However, the insect molecules that interact with viral proteins, such as GN, during infection and dissemination in thrips vector tissues are unknown. The goals of this project were to identify TSWV-interacting proteins (TIPs) that interact directly with TSWV GN and to localize the expression of these proteins in relation to virus in thrips tissues of principal importance along the route of dissemination. We report here the identification of six TIPs from first-instar larvae (L1), the most acquisition-efficient developmental stage of the thrips vector. Sequence analyses of these TIPs revealed homology to proteins associated with the infection cycle of other vector-borne viruses. Immunolocalization of the TIPs in L1 revealed robust expression in the midgut and salivary glands of F. occidentalis, the tissues most important during virus infection, replication, and plant inoculation. The TIPs and GN interactions were validated using protein-protein interaction assays. Two of the thrips proteins, endocuticle structural glycoprotein and cyclophilin, were found to be consistent interactors with GN These newly discovered thrips protein-GN interactions are important for a better understanding of the transmission mechanism of persistent propagative plant viruses by their vectors, as well as for developing new strategies of insect pest management and virus resistance in plants.IMPORTANCE Thrips-transmitted viruses cause devastating losses to numerous food crops worldwide. For negative-sense RNA viruses that infect plants, the arthropod serves as a host as well by supporting virus replication in specific tissues and organs of the vector. The goal of this work was to identify thrips proteins that bind directly to the viral attachment protein and thus may play a role in the infection cycle in the insect. Using the model plant bunyavirus tomato spotted wilt virus (TSWV), and the most efficient thrips vector, we identified and validated six TSWV-interacting proteins from Frankliniella occidentalis first-instar larvae. Two proteins, an endocuticle structural glycoprotein and cyclophilin, were able to interact directly with the TSWV attachment protein, GN, in insect cells. The TSWV GN-interacting proteins provide new targets for disrupting the viral disease cycle in the arthropod vector and could be putative determinants of vector competence.

Plant Virus-Insect Vector Interactions: Current and Potential Future Research Directions

Acquisition and transmission by an insect vector is central to the infection cycle of the majority of plant pathogenic viruses. Plant viruses can interact with their insect host in a variety of ways including both non-persistent and circulative transmission; in some cases, the latter involves virus replication in cells of the insect host. Replicating viruses can also elicit both innate and specific defense responses in the insect host. A consistent feature is that the interaction of the virus with its insect host/vector requires specific molecular interactions between virus and host, commonly via proteins. Understanding the interactions between plant viruses and their insect host can underpin approaches to protect plants from infection by interfering with virus uptake and transmission. Here, we provide a perspective focused on identifying novel approaches and research directions to facilitate control of plant viruses by better understanding and targeting virus-insect molecular interactions. We also draw parallels with molecular interactions in insect vectors of animal viruses, and consider technical advances for their control that may be more broadly applicable to plant virus vectors.

Stable Reference Gene Selection for RT-qPCR Analysis in Nonviruliferous and Viruliferous Frankliniella occidentalis

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a reliable technique for measuring and evaluating gene expression during variable biological processes. To facilitate gene expression studies, normalization of genes of interest relative to stable reference genes is crucial. The western flower thrips Frankliniella occidentalis (Pergande) (Thysanoptera: Thripidae), the main vector of tomato spotted wilt virus (TSWV), is a destructive invasive species. In this study, the expression profiles of 11 candidate reference genes from nonviruliferous and viruliferous F. occidentalis were investigated. Five distinct algorithms, geNorm, NormFinder, BestKeeper, the ΔC_t method, and RefFinder, were used to determine the performance of these genes. geNorm, NormFinder, BestKeeper, and RefFinder identified heat shock protein 70 (HSP70), heat shock protein 60 (HSP60), elongation factor 1 α, and ribosomal protein l32 (RPL32) as the most stable reference genes, and the ΔC_t method identified HSP60, HSP70, RPL32, and heat shock protein 90 as the most stable reference genes. Additionally, two reference genes were sufficient for reliable normalization in nonviruliferous and viruliferous F. occidentalis. This work provides a foundation for investigating the molecular mechanisms of TSWV and F. occidentalis interactions.

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

Heliothine pests such as the tobacco budworm, Heliothis virescens (F.), pose a significant threat to production of a variety of crops and ornamental plants and are models for developmental and physiological studies. The efforts to develop new control measures for H. virescens, as well as its use as a relevant biological model, are hampered by a lack of molecular resources. The present work demonstrates the utility of next-generation sequencing technologies for rapid molecular resource generation from this species for which lacks a sequenced genome. In order to amass a de novo transcriptome for this moth, transcript sequences generated from Illumina, Roche 454, and Sanger sequencing platforms were merged into a single de novo transcriptome assembly. This pooling strategy allowed a thorough sampling of transcripts produced under diverse environmental conditions, developmental stages, tissues, and infections with entomopathogens used for biological control, to provide the most complete transcriptome to date for this species. Over 138 million reads from the three platforms were assembled into the final set of 63,648 contigs. Of these, 29,978 had significant BLAST scores indicating orthologous relationships to transcripts of other insect species, with the top-hit species being the monarch butterfly (Danaus plexippus) and silkworm (Bombyx mori). Among identified H. virescens orthologs were immune effectors, signal transduction pathways, olfactory receptors, hormone biosynthetic pathways, peptide hormones and their receptors, digestive enzymes, and insecticide resistance enzymes. As an example, we demonstrate the utility of this transcriptomic resource to study gene expression profiling of larval midguts and detect transcripts of putative Bacillus thuringiensis (Bt) Cry toxin receptors. The substantial molecular resources described in this study will facilitate development of H. virescens as a relevant biological model for functional genomics and for new biological experimentation needed to develop efficient control efforts for this and related Noctuid pest moths.

Evaluation and validation of reference genes for qRT-PCR normalization in Frankliniella occidentalis (Thysanoptera: Thripidae)

Quantitative real time PCR (qRT-PCR) has emerged as a reliable and reproducible technique for studying gene expression analysis. For accurate results, the normalization of data with reference genes is particularly essential. Once the transcriptome sequencing of Frankliniella occidentalis was completed, numerous unigenes were identified and annotated. Unfortunately, there are no studies on the stability of reference genes used in F. occidentalis. In this work, seven candidate reference genes, including actin, 18S rRNA, H3, tubulin, GAPDH, EF-1 and RPL32, were evaluated for their suitability as normalization genes under different experimental conditions using the statistical software programs BestKeeper, geNorm, Normfinder and the comparative ΔCt method. Because the rankings of the reference genes provided by each of the four programs were different, we chose a user-friendly web-based comprehensive tool RefFinder to get the final ranking. The result demonstrated that EF-1 and RPL32 displayed the most stable expression in different developmental stages; RPL32 and GAPDH showed the most stable expression at high temperatures, while 18S and EF-1 exhibited the most stable expression at low temperatures. In this study, we validated the suitable reference genes in F. occidentalis for gene expression profiling under different experimental conditions. The choice of internal standard is very important in the normalization of the target gene expression levels, thus validating and selecting the best genes will help improve the quality of gene expression data of F. occidentalis. What is more, these validated reference genes could serve as the basis for the selection of candidate reference genes in other insects.

Developing conversed microsatellite markers and their implications in evolutionary analysis of the Bemisia tabaci complex

The study of population genetics among the Bemisia tabaci complex is limited due to the lack of conserved molecular markers. In this study, 358, 433 and 322 new polynucleotide microsatellites are separately identified from the transcriptome sequences of three cryptic species of the B. tabaci complex. The cross species transferability of 57 microsatellites was then experimentally validated. The results indicate that these markers are conserved and have high inter-taxon transferability. Thirteen markers were employed to assess the genetic relationships among six cryptic species of the B. tabaci complex. To our surprise, the inferred phylogeny was consistent with that of mitochondrial COI sequences, indicating that microsatellites have the potential to distinguish species of the B. tabaci complex. Our results demonstrate that development of microsatellites from transcriptome data is a fast and cost-effective approach. These markers can be used to analyze the population genetics and evolutionary patterns of the B. tabaci complex.

The NSs protein of tomato spotted wilt virus is required for persistent infection and transmission by Frankliniella occidentalis

Tomato spotted wilt virus (TSWV) is the type member of tospoviruses (genus Tospovirus), plant-infecting viruses that cause severe damage to ornamental and vegetable crops. Tospoviruses are transmitted by thrips in the circulative propagative mode. We generated a collection of NSs-defective TSWV isolates and showed that TSWV coding for truncated NSs protein could not be transmitted by Frankliniella occidentalis. Quantitative reverse transcription (RT)-PCR and immunostaining of individual insects detected the mutant virus in second-instar larvae and adult insects, demonstrating that insects could acquire and accumulate the NSs-defective virus. Nevertheless, adults carried a significantly lower viral load, resulting in the absence of transmission. Genome sequencing and analyses of reassortant isolates showed genetic evidence of the association between the loss of competence in transmission and the mutation in the NSs coding sequence. Our findings offer new insight into the TSWV-thrips interaction and Tospovirus pathogenesis and highlight, for the first time in the Bunyaviridae family, a major role for the S segment, and specifically for the NSs protein, in virulence and efficient infection in insect vector individuals.

IMPORTANCE

Our work is the first to show a role for the NSs protein in virus accumulation in the insect vector in the Bunyaviridae family: demonstration was obtained for the system TSWV-F. occidentalis, arguably one of the most damaging combination for vegetable crops. Genetic evidence of the involvement of the NSs protein in vector transmission was provided with multiple approaches.

Analysis of the salivary gland transcriptome of Frankliniella occidentalis

Saliva is known to play a crucial role in insect feeding behavior and virus transmission. Currently, little is known about the salivary glands and saliva of thrips, despite the fact that Frankliniella occidentalis (Pergande) (the western flower thrips) is a serious pest due to its destructive feeding, wide host range, and transmission of tospoviruses. As a first step towards characterizing thrips salivary gland functions, we sequenced the transcriptome of the primary salivary glands of F. occidentalis using short read sequencing (Illumina) technology. A de novo-assembled transcriptome revealed 31,392 high quality contigs with an average size of 605 bp. A total of 12,166 contigs had significant BLASTx or tBLASTx hits (E≤1.0E-6) to known proteins, whereas a high percentage (61.24%) of contigs had no apparent protein or nucleotide hits. Comparison of the F. occidentalis salivary gland transcriptome (sialotranscriptome) against a published F. occidentalis full body transcriptome assembled from Roche-454 reads revealed several contigs with putative annotations associated with salivary gland functions. KEGG pathway analysis of the sialotranscriptome revealed that the majority (18 out of the top 20 predicted KEGG pathways) of the salivary gland contig sequences match proteins involved in metabolism. We identified several genes likely to be involved in detoxification and inhibition of plant defense responses including aldehyde dehydrogenase, metalloprotease, glucose oxidase, glucose dehydrogenase, and regucalcin. We also identified several genes that may play a role in the extra-oral digestion of plant structural tissues including β-glucosidase and pectin lyase; and the extra-oral digestion of sugars, including α-amylase, maltase, sucrase, and α-glucosidase. This is the first analysis of a sialotranscriptome for any Thysanopteran species and it provides a foundational tool to further our understanding of how thrips interact with their plant hosts and the viruses they transmit.

Proteomic analysis of Frankliniella occidentalis and differentially expressed proteins in response to tomato spotted wilt virus infection

Tomato spotted wilt virus (TSWV) is transmitted by Frankliniella occidentalis in a persistent propagative manner. Despite the extensive replication of TSWV in midgut and salivary glands, there is little to no pathogenic effect on F. occidentalis. We hypothesize that the first-instar larva (L1) of F. occidentalis mounts a response to TSWV that protects it from pathogenic effects caused by virus infection and replication in various insect tissues. A partial thrips transcriptome was generated using 454-Titanium sequencing of cDNA generated from F. occidentalis exposed to TSWV. Using these sequences, the L1 thrips proteome that resolved on a two-dimensional gel was characterized. Forty-seven percent of the resolved protein spots were identified using the thrips transcriptome. Real-time quantitative reverse transcriptase PCR (RT-PCR) analysis of virus titer in L1 thrips revealed a significant increase in the normalized abundance of TSWV nucleocapsid RNA from 2 to 21 h after a 3-h acquisition access period on virus-infected plant tissue, indicative of infection and accumulation of virus. We compared the proteomes of infected and noninfected L1s to identify proteins that display differential abundances in response to virus. Using four biological replicates, 26 spots containing 37 proteins were significantly altered in response to TSWV. Gene ontology assignments for 32 of these proteins revealed biological roles associated with the infection cycle of other plant- and animal-infecting viruses and antiviral defense responses. Our findings support the hypothesis that L1 thrips display a complex reaction to TSWV infection and provide new insights toward unraveling the molecular basis of this interaction.

Development and characterization of 18 novel EST-SSRs from the western flower Thrips, Frankliniella occidentalis (Pergande)

The western flower thrips, Frankliniella occidentalis (Pergande), is an invasive species and the most economically important pest within the insect order Thysanoptera. For a better understanding of the genetic makeup and migration patterns of F. occidentalis throughout the world, we characterized 18 novel polymorphic EST-derived microsatellites. The mutational mechanism of these EST-SSRs was also investigated to facilitate the selection of appropriate combinations of markers for population genetic studies. Genetic diversity of these novel markers was assessed in 96 individuals from three populations in China (Harbin, Dali, and Guiyang). The results showed that all these 18 loci were highly polymorphic; the number of alleles ranged from 2 to 15, with an average of 5.50 alleles per locus. The observed (H_O) and expected (H_E) heterozygosities ranged from 0.072 to 0.707 and 0.089 to 0.851, respectively. Furthermore, only two locus/population combinations (WFT144 in Dali and WFT50 in Guiyang) significantly deviated from Hardy–Weinberg equilibrium (HWE). Pairwise F_ST analysis showed a low but significant differentiation (0.026 < F_ST < 0.032) among all three pairwise population comparisons. Sequence analysis of alleles per locus revealed a complex mutational pattern of these EST-SSRs. Thus, these EST-SSRs are useful markers but greater attention should be paid to the mutational characteristics of these microsatellites when they are used in population genetic studies.