Insight into the Salivary Gland Transcriptome of Lygus lineolaris (Palisot de Beauvois)

Exploring the contribution of the salivary gland and midgut to digestion in the swede midge (Contarinia nasturtii) through a genomics-guided approach

The larvae of Contarinia nasturtii (Kieffer) (Diptera: Cecidomyiidae), the swede midge, targets the meristem of brassica crops where they induce the formation of galls and disrupt seed and vegetable production. Previously, we examined the salivary gland transcriptome of newly-hatched first instar larvae as they penetrated the host and initiated gall formation. Here we examine the salivary gland and midgut transcriptome of third instar larvae and provide evidence for cooperative nutrient acquisition beginning with secretion of enzymes and feeding facilitators followed by gastrointestinal digestion. Sucrose, presumably obtained from the phloem, appeared to be a major nutrient source as several α-glucosidases (sucrases, maltases) and β-fructofuranosidases (invertases) were identified. Genes encoding β-fructofuranosidases/invertases were among the most highly expressed in both tissues and represented two distinct gene families that may have originated via horizontal gene transfer from bacteria. The importance of the phloem as a nutrient source is underscored by the expression of genes encoding regucalcin and ARMET (arginine-rich mutated in early stages of tumor) which interfere with calcium signalling and prevent sieve tube occlusion. Lipids, proteins, and starch appear to serve as a secondary nutrient sources. Genes encoding enzymes involved in the detoxification of glucosinolates (myrosinases, arylsulfatases, and glutathione-S-transferases) were expressed indicative of Brassicaceae host specialization. The midgut expressed simple peritrophins and mucins typical of those found in Type II peritrophic matrices, the first such description for a gall midge.

Comparative transcriptomic analysis of salivary glands between the zoophytophagous Cyrtorhinus lividipennis and the phytozoophagous Apolygus lucorum

BACKGROUND

Saliva plays a crucial role in shaping the feeding behavior of insects, involving processes such as food digestion and the regulation of interactions between insects and their hosts. Cyrtorhinus lividipennis serves as a predominant natural enemy of rice pests, while Apolygus lucorum, exhibiting phytozoophagous feeding behavior, is a destructive agricultural pest. In this study, a comparative transcriptome analysis, incorporating the published genomes of C.lividipennis and A.lucorum, was conducted to reveal the role of salivary secretion in host adaptation.

RESULTS

In contrast to A.lucorum, C.lividipennis is a zoophytophagous insect. A de novo genome analysis of C.lividipennis yielded 19,706 unigenes, including 16,217 annotated ones. On the other hand, A.lucorum had altogether 20,111 annotated genes, as obtained from the published official gene set (20,353 unigenes). Functional analysis of the top 1,000 salivary gland (SG)-abundant genes in both insects revealed that the SG was a dynamically active tissue engaged in protein synthesis and secretion. Predictions of other tissues and signal peptides were compared. As a result, 94 and 157 salivary proteins were identified in C.lividipennis and A.lucorum, respectively, and were categorized into 68 and 81 orthogroups. Among them, 26 orthogroups were shared, potentially playing common roles in digestion and detoxification, including several venom serine proteases. Furthermore, 42 and 55 orthogroups were exclusive in C.lividipennis and A.lucorum, respectively, which were exemplified by a hyaluronidase in C.lividipennis that was associated with predation, while polygalacturonases in A.lucorum were involved in mesophyll-feeding patterns.

CONCLUSIONS

Findings in this study provide a comprehensive insight into saliva secretions in C.lividipennis and A.lucorum via a transcriptome approach, reflecting the intricate connections between saliva secretions and feeding behaviors. It is found that conserved salivary secretions are involved in shaping the overlapping feeding patterns, while a plethora of unique salivary secretions may drive the evolution of specific feeding behaviors crucial for their survival. These results enhance our understanding of the feeding mechanisms in different insects from the perspective of saliva and contribute to future environmentally friendly pest control by utilizing predatory insects.

The resistance of the jujube (Ziziphus jujuba) to the devastating insect pest Apolygus lucorum (Hemiptera, Insecta) involves the jasmonic acid signaling pathway

Apolygus lucorum (Hemiptera, Insecta), cosmopolitan true bug, is a major pest of the Chinese jujube (Ziziphus jujuba). To propose control measures of A. lucorum, we investigated the molecular mechanisms of resistance in two varieties of jujube (wild jujube and winter jujube) with different sensitivities to this pest. We monitored changes of two species of jujube in the transcriptome, jasmonic acid (JA) and salicylic acid (SA) content, and the expression of genes involved in signaling pathways. The preference of A. lucorum for jujube with exogenous SA and methyl jasmonate (MeJA) were also examined. The results showed that wild jujube leaves infested by A. lucorum showed stronger resistance and non-selectivity to A. lucorum than winter jujube. By comparing data from the A. lucorum infested plants with the control, A total of 438 and 796 differentially expressed genes (DEGs) were found in winter and wild jujube leaves, respectively. GO analysis revealed that biological process termed "plant-pathogen interactions", "plant hormone transduction" and "phenylpropanoid biosynthesis". Most of DEGs enriched in JA pathways were upregulated, while most DEGs of SA pathways were downregulated. A. lucorum increased the JA content but decreased the SA content in jujube. Consistently, the JA and SA contents in winter jujube were lower than those in wild jujube leaves. The key genes ZjFAD3, ZjLOX, ZjAOS, ZjAOC3 and ZjAOC4 involved in JA synthesis of jujube leaves were significantly up-regulated after A. lucorum infestation, especially the expression and up-regulation ratio of ZjFAD3, ZjLOX and ZjAOS in wild jujube were significantly higher than those in winter jujube. MeJA-treated jujube showed an obvious repellent effect on A. lucorum. Based on these findings, we conclude that A. lucorum infestation of jujube induced the JA pathway and suppressed the SA pathway. In jujube leaves the ZjFAD3, ZjLOX and ZjAOS played important roles in increasing of JA content in jujube leaves. Thus, JA played an important role in repelling and resisting against A. lucorum in jujube.

A chromosome scale assembly of the tarnished plant bug, Lygus lineolaris (Palisot de Beauvois), genome

OBJECTIVE

The tarnished plant bug (TPB), Lygus lineolaris (Palisot de Beauvois) (Hemiptera: Miridae), is a pest damaging many cultivated crops in North America. Although partial transcriptome data are available for this pest, a genome assembly was not available for this species. This assembly of a high-quality chromosome-length genome of TPB is aimed to develop the genetic resources that can provide the foundation required for advancing research on this species.

RESULTS

The initial genome of TPB assembled with paired-end nucleotide sequences generated with Illumina technology was scaffolded with Illumina HiseqX reads generated from a proximity ligated (HiC) library to obtain a high-quality genome assembly. The final assembly contained 3963 scaffolds longer than 1 kbp to yield a genome of 599.96 Mbp. The N50 of the TPB genome assembly was 35.64 Mbp and 98.68% of the genome was assembled into 17 scaffolds larger than 1 Mbp. This megabase scaffold number is the same as the number of chromosomes observed in karyotyping of this insect. The TPB genome is known to have high repetitive DNA content, and the reduced assembled genome size compared to flowcytometric estimates of approximately 860 Mbp may be due to the collapsed assembly of highly similar regions.

High-quality chromosome-level scaffolds of the plant bug Pachypeltis micranthus provide insights into the availability of Mikania micrantha control

BACKGROUND

The plant bug, Pachypeltis micranthus Mu et Liu (Hemiptera: Miridae), is an effective potential biological control agent for Mikania micrantha H.B.K. (Asteraceae; one of the most notorious invasive weeds worldwide). However, limited knowledge about this species hindered its practical application and research. Accordingly, sequencing the genome of this mirid bug holds great significance in controlling M. micrantha.

RESULTS

Here, 712.72 Mb high-quality chromosome-level scaffolds of P. micranthus were generated, of which 707.51 Mb (99.27%) of assembled sequences were anchored onto 15 chromosome-level scaffolds with contig N50 of 16.84 Mb. The P. micranthus genome had the highest GC content (42.43%) and the second highest proportion of repetitive sequences (375.82 Mb, 52.73%) than the three other mirid bugs (i.e., Apolygus lucorum, Cyrtorhinus lividipennis, and Nesidiocoris tenuis). Phylogenetic analysis showed that P. micranthus clustered with other mirid bugs and diverged from the common ancestor approximately 200 million years ago. Gene family expansion and/or contraction were analyzed, and significantly expanded gene families associated with P. micranthus feeding and adaptation to M. micrantha were manually identified. Compared with the whole body, transcriptome analysis of the salivary gland revealed that most of the upregulated genes were significantly associated with metabolism pathways and peptidase activity, particularly among cysteine peptidase, serine peptidase, and polygalacturonase; this could be one of the reasons for precisely and highly efficient feeding by the oligophagous bug P. micranthus on M. micrantha.

CONCLUSION

Collectively, this work provides a crucial chromosome-level scaffolds resource to study the evolutionary adaptation between mirid bug and their host. It is also helpful in searching for novel environment-friendly biological strategies to control M. micrantha.

Comparative analyses of the venom components in the salivary gland transcriptomes and saliva proteomes of some heteropteran insects

Salivary gland-specific transcriptomes of nine heteropteran insects with distinct feeding strategies (predaceous, hematophagous, and phytophagous) were analyzed and annotated to compare and identify the venom components as well as their expression profiles. The transcriptional abundance of venom genes was verified via quantitative real-time PCR. Hierarchical clustering of 30 representative differentially expressed venom genes from the nine heteropteran species revealed unique groups of salivary gland-specific genes depending on their feeding strategy. The commonly transcribed genes included a paralytic neurotoxin (arginine kinase), digestive enzymes (cathepsin and serine protease), an anti-inflammatory protein (cystatin), hexamerin, and an odorant binding protein. Both predaceous and hematophagous (bed bug) heteropteran species showed relatively higher transcription levels of genes encoding proteins involved in proteolysis and cytolysis, whereas phytophagous heteropterans exhibited little or no expression of these genes, but had a high expression of vitellogenin, a multifunctional allergen. Saliva proteomes from four representative species were also analyzed. All venom proteins identified via saliva proteome analysis were annotated using salivary gland transcriptome data. The proteomic expression profiles of venom proteins were in good agreement with the salivary gland-specific transcriptomic profiles. Our results indicate that profiling of the salivary gland transcriptome provides important information on the composition and evolutionary features of venoms depending on their feeding strategy.

Stomata-mediated interactions between plants, herbivores, and the environment

Stomata play a central role in plant responses to abiotic and biotic stresses. Existing knowledge regarding the roles of stomata in plant stress is centered on abiotic stresses and plant-pathogen interactions, but how stomata influence plant-herbivore interactions remains largely unclear. Here, we summarize the functions of stomata in plant-insect interactions and highlight recent discoveries of how herbivores manipulate plant stomata. Because stomata are linked to interrelated physiological processes in plants, herbivory-induced changes in stomatal dynamics might have cellular, organismic, and/or even community-level impacts. We summarize our current understanding of how stomata mediate plant responses to herbivory and environmental stimuli, propose how herbivores may influence these responses, and identify key knowledge gaps in plant-herbivore interactions.

Biology, Ecology, and Pest Management of the Tarnished Plant Bug, Lygus lineolaris (Palisot de Beauvois) in Southern Row Crops

Simple Summary

The tarnished plant bug, Lygus lineolaris, is a polyphagous, sap-feeder that causes significant economic damage in several field crops, especially cotton (Gossypium hirsutum L.) in the mid-southern United States. In 2020, it was reported that 4.8 million acres of cotton were infested by Lygus spp. in the United States. A broad host range, polyphagous feeding behavior and high mobility of this pest along with resistance development to conventional pesticides helped them establish as a significant pest of concern for cotton growers in the mid-south. Since the publication of a review by Layton (2000) on damage caused by Lygus lineolaris, many new research studies have been published on the Lygus biology, ecology, and integrated pest management strategies. A comprehensive review paper that summarizes these latest research developments and Lygus management strategies will be useful for researchers and cotton growers. In this review, we report and discuss the latest developments in Lygus research and the new control strategies that have been developed in the last two decades.

Abstract

The tarnished plant bug, Lygus lineolaris (Palisot de Beauvois), (Hemiptera: Miridae) is considered the most damaging pest of cotton (Gossypium hirsutum L.) in the mid-southern United States, although it is established throughout the United States, southern Canada, and northern Mexico. The introduction of transgenic crops for the control of moths in the Heliothine complex and eradication of the boll weevil, Anthonomus grandis, from much of the United States led to greatly reduced pesticide use in cotton fields, which allowed L. lineolaris to emerge as a new primary pest of cotton in the mid-southern United States. Since the publication of a review by Layton (2000) on damage caused by Lygus lineolaris, many new studies have been published on the changes in host range, population dynamics, sampling methods and thresholds, cultural practices, sex pheromones and attractant blends, novel pesticides and insecticide resistance mechanisms, olfactory and feeding behaviors, introduction of biological control agents, host-plant resistance mechanisms, and new molecular and genetic tools for integrated pest management of Lygus species in cotton and other important crops. Here, we review and discuss the latest developments in L. lineolaris research in the last two decades.

A stationary tweezer platform for high throughput dissections of minute arthropods and extirpation of their minute organs

A homemade platform satisfied the need for fast, efficient dissection of minute arthropods and extirpation of their key organs, such as salivary glands and midguts, involved in agricultural disease transmission pathways. With its implementation, ~200 organs could be extirpated per 8 h workday while the subjects are submerged in protein or transcript protectant. A vacuum wand is used to capture insects and position them in the field of view. Two stationary tweezers are positioned on an adjustable scaffold that spans the microscope stage transversely such that their tips, and the insects they immobilize, can be submerged in select dissection media. High tensile strength fishlines are attached to the stationary tweezers for opening and closing with the 5th fingers while hand-held dissection tweezers load insects from the wand to their tines, then extirpate the target organs. Organs are lifted out with glass splints or plastic toothpicks into a final tube of select preservation media for freezing at session end.•Constructed from common retail materials•Adjustable design fits many microscopes•Can also be used in a wide variety of applications, including materials science.

Expression Profiles of Digestive Genes in the Gut and Salivary Glands of Tarnished Plant Bug (Hemiptera: Miridae)

Host plant preference of agricultural pests may shift throughout the growing season, allowing the pests to persist on wild hosts when crops are not available. Lygus Hahn (Hemiptera: Miridae) bugs are severe pests of cotton during flowering and fruiting stages, but can persist on alternative crops, or on weed species. Diversity of digestive enzymes produced by salivary glands and gut tissues play a pivotal role in an organism's ability to utilize various food sources. Polyphagous insects produce an array of enzymes that can process carbohydrates, lipids, and proteins. In this study, the digestive enzyme repertoire of the tarnished plant bug, Lygus lineolaris (Palisot de Beauvois), was identified by high-throughput sequencing followed by cDNA cloning and sequencing. This study identified 87 digestive genes, including 30 polygalacturonases (PG), one β-galactosidase, three α-glucosidases, six β-glucosidases, 28 trypsin-like proteases, three serine proteases, one apyrase-like protease, one cysteine protease, 12 lipases, and two transcripts with low similarity to a xylanase A-like genes. RNA-Seq expression profiles of these digestive genes in adult tarnished plant bugs revealed that 57 and 12 genes were differentially expressed in the salivary gland and gut (≥5-fold, P ≤ 0.01), respectively. All polygalacturonase genes, most proteases, and two xylanase-like genes were differentially expressed in salivary glands, while most of the carbohydrate and lipid processing enzymes were differentially expressed in the gut. Seven of the proteases (KF208689, KF208697, KF208698, KF208699, KF208700, KF208701, and KF208702) were not detected in either the gut or salivary glands.

Extraoral digestion: outsourcing the role of the hemipteran midgut

Extraoral digestion allows for breakdown of dietary components before they reach the midgut for final enzymatic degradation and absorption. In the Hemiptera, this is achieved by the secretion of enzyme-rich fluids from the salivary gland, with the combination of protein and mRNA from these tissues termed the sialome. Separate channels within the hemipteran stylets allow for secretion of saliva and ingestion of predigested material in a non-reflux mechanism. Both feeding mode and diet type influence the composition of the hemipteran sialome, as illustrated by 1) differences in protease abundance between hematophagous and predatory heteropteran sialomes, 2) diet specific aminopeptidase-N genes among aphid biotypes, and 3) adaptation-induced sialome variation in related cicada populations. Despite challenges associated with incomplete genome annotation, -omics analysis of the sialomes of diverse hemipteran species will enhance understanding of both sialome function and the evolution of extraoral digestion within the order.

Comparative sialotranscriptome analysis of the rare Chinese cicada Subpsaltria yangi, with identification of candidate genes related to host-plant adaptation

Little is known about the mechanisms underlying the relationship between changes in salivary proteins and the adaptation of insects to different host-plants. To address this knowledge gap, the transcriptional profiles of salivary glands were compared among three populations of the rare cicada Subpsaltria yangi, in which two populations specialize on Zizyphus jujuba var. spinosa, but the population occurring in the Helan (HL) Mountains is locally specialized on the endemic plant Ephedra lepidosperma. The comparisons indicate that genes related to digestion and detoxification are differentially regulated in populations feeding on different plants, possibly reflecting adaptative changes in salivary proteins of S. yangi in response to different host chemistries. In detail, 38 differentially expressed genes and 21 up-regulated genes related to digestion and detoxification were identified respectively in two pairwise comparisons among the populations using different hosts, with some genes exclusively expressed in the HL population. Our results are consistent with the hypothesis that the host plant shift in the HL population was facilitated by differential regulation of genes related to digestion and detoxification. This study provides new information for elucidating the molecular mechanisms underlying the relationship between changed salivary proteins and the adaptability of plant-feeding insects to novel host plants.

Transcriptomic and Expression Analysis of the Salivary Glands in Brown Planthoppers, Nilaparvata lugens (Hemiptera: Delphacidae)

The brown planthopper (BPH), Nilaparvata lugens Stål (Hemiptera: Delphacidae), is a serious rice pest because of its destructive feeding. The salivary glands of the BPH play a key role in successful feeding. In this study, we explored the salivary gland transcriptome (sialotranscriptome) of adult BPHs using Illumina sequencing and a total of 55,913 transcripts and 45,421 unigenes were obtained. We identified one reference gene RPL9 (Ribosomal protein L9) and 19 salivary protein genes from the BPH sialotranscripome, which were categorized as those involved in sugar metabolism, extra-oral digestion of cell wall components, detoxification, and suppression of plant defenses. Tissue expression profiles of 19 salivary protein genes analysis revealed that the expression level of alpha-glucosidase family 31 had no difference in five tissues, suggesting that it may have functions in the whole-body parts. Glucose dehydrogenase (flavine adenine dinucleotide, quinone)-like was expressed highly in the salivary gland, which might play putative role in insect feeding. Glucose dehydrogenase (acceptor) was expressed the highest level in head without salivary gland. Other salivary protein genes were highly expressed in gut or malpighian tubule, suggesting that they may play roles in digestive and detoxification mechanism. Moreover, we detected RPL9 is one of the genes that is most consistently present for acquisition of gene expression in different tissues. Thus, RPL9 can be a new reference for expression studies of BPH. The obtained BPH sialotranscripome provides a list of genes that have potential roles in feeding and interaction between BPHs and rice plants.

By land, air, and sea: hemipteran diversity through the genomic lens

Thanks to a recent spate of sequencing projects, the Hemiptera are the first hemimetabolous insect order to achieve a critical mass of species with sequenced genomes, establishing the basis for comparative genomics of the bugs. However, as the most speciose hemimetabolous order, there is still a vast swathe of the hemipteran phylogeny that awaits genomic representation across subterranean, terrestrial, and aquatic habitats, and with lineage-specific and developmentally plastic cases of both wing polyphenisms and flightlessness. In this review, we highlight opportunities for taxonomic sampling beyond obvious pest species candidates, motivated by intriguing biological features of certain groups as well as the rich research tradition of ecological, physiological, developmental, and particularly cytogenetic investigation that spans the diversity of the Hemiptera.

Transcriptome analysis of the salivary glands of the grain aphid, Sitobion avenae

Aphid saliva plays important roles in aphid-host interactions, such as assisting aphid digestion, detoxification, activating or suppressing plant defenses. The grain aphid, Sitobion avenae, is one of the most devastating pests of cereals worldwide. In this study, we performed the transcriptome analysis of salivary glands of S. avenae. A total of 33,079 assembled unigenes were identified in the salivary glands of aphids. Of the all obtained unigenes, 15,833(47.86%) and 10,829(32.73%) unigenes showed high similarity to known proteins in Nr and Swiss-Prot databases respectively. 526 unigenes were predicted to encode secretory proteins, including some digestive and detoxifying enzymes and potential effectors. The RT-PCR and RT-qPCR results showed that all of the 15 most highly expressed putative secretory proteins specifically expressed in salivary glands. Interestingly, 11 of the 15 most highly expressed putative secretory proteins were still not matched to function-known proteins. We also detected the expression of 9 interested putative secretory proteins in aphid different tissues, including some digestive and detoxifying enzymes, effectors and Ca2+ binding proteins. The results showed that only glutathione-S-transferase 1 was specifically expressed in salivary glands. These findings provide a further insight into the identification of potential effectors involving in aphid-cereals interactions.

The genome of the cotton bacterial blight pathogen Xanthomonas citri pv. malvacearum strain MSCT1

Xanthomonas citri pv. malvacearum is a major pathogen of cotton, Gossypium hirsutum L.. In this study we report the complete genome of the X. citri pv. malvacearum strain MSCT1 assembled from long read DNA sequencing technology. The MSCT1 genome is the first X. citri pv. malvacearum genome with complete coding regions for X. citri pv. malvacearum transcriptional activator-like effectors. In addition functional and structural annotations are presented in this study that will provide a foundation for future pathogenesis studies with MSCT1.

Transcriptome analysis of three cotton pests reveals features of gene expressions in the mesophyll feeder Apolygus lucorum

The green mirid bug Apolygus lucorum is an agricultural pest that is known to cause damage to more than 150 plant species. Here, we report the transcriptomes of A. lucorum at three different developmental stages (the second and fifth instar nymphs and adults). A total of 98,236 unigenes with an average length of 1,335 nt was obtained, of which 50,640 were annotated, including those encoding digestive enzymes and cytochrome P450s. Comparisons with cotton bollworm and cotton aphid transcriptomes revealed distinct features of A. lucorum as a mesophyll feeder. The gene expression dynamics varied during development from young nymphs to adults. The high-quality transcriptome data and the gene expression dynamics reported here provide valuable data for a more comprehensive understanding of the physiology and development of mirid bugs, and for mining targets for their control.

Functional analysis of two polygalacturonase genes in Apolygus lucorum associated with eliciting plant injury using RNA interference

Salivary enzymes of many piercing-sucking insects lead to host plant injury. The salivary enzymes, polygalacturonase (PGs), act in insect feeding. PG family genes have been cloned from the mirid bug Apolygus lucorum, a pest of cotton and other host crops in China. We investigated the function of two PG genes that are highly expressed in A. lucorum nymphs (PG3-4) and adults (PG3-5), using siRNA injection-based RNA interference (RNAi). Accumulation of mRNA encoding both genes and their cognate proteins was significantly reduced (>60%) in experimental compared control green fluorescent protein (GFP) siRNA-treated mirids at 48 h post injection. Injury levels of cotton buds were also significantly reduced after injecting saliva isolated from PG3-4 and PG3-5 siRNA-treated A. lucorum. These results demonstrate that these two PG act in A. lucorum elicitation of plant injury.

Identification of Genes Potentially Responsible for extra-Oral Digestion and Overcoming Plant Defense from Salivary Glands of the Tarnished Plant Bug (Hemiptera: Miridae) Using cDNA Sequencing

Saliva is known to play a crucial role in tarnished plant bug (TPB, Lygus lineolaris [Palisot de Beauvois]) feeding. By facilitating the piercing, the enzyme-rich saliva may be used for extra-oral digestion and for overcoming plant defense before the plant fluids are ingested by TPBs. To identify salivary gland genes, mRNA was extracted from salivary glands and cDNA library clones were sequenced. A de novo-assembling of 7,000 Sanger sequences revealed 666 high-quality unique cDNAs with an average size of 624 bp, in which the identities of 347 cDNAs were determined using Blast2GO. Kyoto Encyclopedia of Genes and Genomes analysis indicated that these genes participate in eighteen metabolic pathways. Identifications of large number of enzyme genes in TPB salivary glands evidenced functions for extra-oral digestion and feeding damage mechanism, including 45 polygalacturonase, two α- amylase, one glucosidase, one glycan enzyme, one aminopeptidase, four lipase, and many serine protease cDNAs. The presence of multiple transcripts, multigene members, and high abundance of cell wall degradation enzymes (polygalacturonases) indicated that the enzyme-rich saliva may cause damage to plants by breaking down plant cell walls to make nutrients available for feeding. We also identified genes potentially involved in insect adaptation and detoxifying xenobiotics that may allow insects to overcome plant defense responses, including four glutathione S-transferases, three esterases, one cytochrome P450, and several serine proteases. The gene profiles of TPB salivary glands revealed in this study provides a foundation for further understanding and potential development of novel enzymatic inhibitors, or other RNAi approaches that may interrupt or minimize TPB feeding damage.