Transcriptome analysis of the salivary glands of potato leafhopper, Empoasca fabae

Transcriptomic response of citrus psyllid salivary glands to the infection of citrus Huanglongbing pathogen

The Asian citrus psyllid, Diaphorina citri Kuwayama (Hemiptera: Psyllidae), is the key vector insect transmitting the Candidatus Liberibacter asiaticus (CLas) bacterium that causes the devastating citrus greening disease (Huanglongbing, HLB) worldwide. The D. citri salivary glands (SG) exhibit an important barrier against the transmission of HLB pathogen. However, knowledge on the molecular mechanism of SG defence against CLas infection is still limited. In the present study, we compared the SG transcriptomic response of CLas-free and CLas-infected D. citri using an illumine paired-end RNA sequencing. In total of 861 differentially expressed genes (DEGs) in the SG upon CLas infection, including 202 upregulated DEGs and 659 downregulated DEGs were identified. Functional annotation analysis showed that most of the DEGs were associated with cellular processes, metabolic processes, and the immune response. Gene ontology and Kyoto Encyclopaedia of Genes and Genomes enrichment analyses revealed that these DEGs were enriched in pathways involving carbohydrate metabolism, amino acid metabolism, the immune system, the digestive system, the lysosome, and endocytosis. A total of 16 DEGs were randomly selected to further validate the accuracy of RNA-Seq dataset by reverse-transcription quantitative polymerase chain reaction. This study provides substantial transcriptomic information regarding the SG of D. citri in response to CLas infection, which may shed light on the molecular interaction between D. citri and CLas, and provides new ideas for the prevention and control of citrus psyllid.

Identification of salivary proteins of the cowpea aphid Aphis craccivora by transcriptome and LC-MS/MS analyses

Aphid salivary proteins mediate the interaction between aphids and their host plants. Moreover, these proteins facilitate digestion, detoxification of secondary metabolites, as well as activation and suppression of plant defenses. The cowpea aphid, Aphis craccivora, is an important sucking pest of leguminous crops worldwide. Although aphid saliva plays an important role in aphid plant interactions, knowledge of the cowpea aphid salivary proteins is limited. In this study, we performed transcriptomic and LC-MS/MS analyses to identify the proteins present in the salivary glands and saliva of A. craccivora. A total of 1,08,275 assembled transcripts were identified in the salivary glands of aphids. Of all these assembled transcripts, 53,714 (49.11%) and 53,577 (49.48%) transcripts showed high similarity to known proteins in the Nr and UniProt databases, respectively. A total of 2159 proteins were predicted as secretory proteins from the salivary gland transcriptome dataset, which contain digestive enzymes, detoxification enzymes, previously known effectors and elicitors, and potential proteins whose functions have yet to be determined. The proteomic analysis of aphid saliva resulted in the identification of 171 proteins. Tissue-specific expression of selected genes using RT-PCR showed that three genes were expressed only in the salivary glands. Overall, our results provide a comprehensive repertoire of cowpea aphid salivary proteins from the salivary gland and saliva, which will be a good resource for future effector functional studies and might also be useful for sustainable aphid management.

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.

Screening Candidate Effectors of the Bean Bug Riptortus pedestris by Proteomic and Transcriptomic Analyses

The damage of Riptortus pedestris is exceptional by leading soybean plants to keep green in late autumn. Identification of the salivary proteins is essential to understand how the pest-plant interaction occurs. Here, we have tried to identify them by a combination of proteomic and transcriptomic analyses. The transcriptomes of salivary glands from R. pedestris males, females and nymphs showed about 28,000 unigenes, in which about 40% had open reading frames (ORFs). Therefore, the predicted proteins in the transcriptomes with secretion signals were obtained. Many of the top 1,000 expressed transcripts were involved in protein biosynthesis and transport, suggesting that the salivary glands produce a rich repertoire of proteins. In addition, saliva of R. pedestris males, females and nymphs was collected and proteins inside were identified. In total, 155, 20, and 11 proteins were, respectively, found in their saliva. We have tested the tissue-specific expression of 68 genes that are likely to be effectors, either because they are homologs of reported effectors of other sap-feeding arthropods, or because they are within the top 1,000 expressed genes or found in the salivary proteomes. Their potential functions in regulating plant defenses were discussed. The datasets reported here represent the first step in identifying effectors of R. pedestris.

Detecting Ingested Host Plant DNA in Potato Leafhopper (Hemiptera: Cicadellidae): Potential Use of Molecular Markers for Gut Content Analysis

Detection of host plant DNA from sap-feeding insects can be challenging due to potential low concentration of ingested plant DNA. Although a few previous studies have demonstrated the possibility of detecting various fragments of plant DNA from some sap-feeders, there are no protocols available for potato leafhopper, Empoasca fabae (Harris) (Hemiptera: Cicadellidae), a significant agricultural pest. In this study we focused on optimizing a DNA-based method for host plant identification of E. fabae and investigating the longevity of the ingested plant DNA as one of the potential applications of the protocol. We largely utilized and modified our previously developed PCR-based method for detecting host plant DNA from grasshopper and the spotted lanternfly gut contents. We have demonstrated that the trnL (UAA) gene can be successfully utilized for detecting ingested host plant DNA from E. fabae and determining plant DNA longevity. The developed protocol is a relatively quick and low-cost method for detecting plant DNA from E. fabae. It has a number of important applications-from determining host plants and dispersal of E. fabae to developing effective pest management strategies.

Food Resource Sharing of Alder Leaf Beetle Specialists (Coleoptera: Chrysomelidae) as Potential Insect-Plant Interface for Horizontal Transmission of Endosymbionts

Recent studies suggest that endosymbionts of herbivore insects can be horizontally transferred to other herbivores feeding on the same host plants, whereby the plant acts as an intermediate stage in the chain of transmission. If this mechanism operates, it is also expected that insect communities sharing the same host plant will have higher chances to share their endosymbionts. In this study, we use a high-throughput 16S rRNA metabarcoding approach to investigate the presence, diversity, and potential sharing of endosymbionts in several species of leaf beetles (Coleoptera: Chrysomelidae) of a local community specialized on an alder diet in North America. Rickettsia and Wolbachia were predominant in the sample, with strong evidence for each species having their own dominant infection, of either or both types of bacteria. However, all species shared a much lower proportion of a particular Wolbachia type, compatible with the same strain dominant in one of the species of leaf beetles. Crucially, the same 16S rRNA haplotype of Wolbachia was found on alder leaf extracts. The combined evidence and the absence of this strain in a syntopic species of leaf beetle feeding on a different host plant support the hypothesis that at least the initial stages of the mechanism that would allow horizontal transmission of endosymbionts across species feeding on the same plant is possible. The accessibility and characteristics of endosymbiont associations of this system make it suitable for deeper analyses of their diversity and transmission in natural conditions.

In silico identification of effector proteins from generalist herbivore Spodoptera litura

BACKGROUND

The common cutworm, Spodoptera litura Fabricius is a leaf and fruit feeding generalist insect of the order Lepidoptera and a destructive agriculture pest. The broad host range of the herbivore is due to its ability to downregulate plant defense across different plants. The identity of Spodoptera litura released effectors that downregulate plant defense are largely unknown. The current study aims to identify genes encoding effector proteins from salivary glands of S. litura (Fab.).

RESULTS

Head and salivary glands of Spodoptera litura were used for de-novo transcriptome analysis and effector prediction. Eight hundred ninety-nine proteins from the head and 330 from salivary gland were identified as secretory proteins. Eight hundred eight proteins from the head and 267 from salivary gland proteins were predicted to be potential effector proteins.

CONCLUSIONS

This study is the first report on identification of potential effectors from Spodoptera litura salivary glands.

A Transcriptomic and Proteomic Analysis of the Diaphorina citri Salivary Glands Reveals Genes Responding to Candidatus Liberibacter asiaticus

The Asian citrus psyllid (ACP), Diaphorina citri Kuwayama, is the principal vector of the Candidatus Liberibacter asiaticus (CLas) bacterium that causes Huanglongbing (HLB) disease. The D. citri salivary glands (SG) is an important barrier to the transmission of CLas. Despite its importance, the transcriptome and proteome of SG defense against CLas are unstudied in D. citri. In the present study, we generated a comparative transcriptome dataset of the SG in infected and uninfected D. citri using an Illumina RNA-Seq technology. We obtained 407 differentially expressed genes (DEGs), including 159 upregulated DEGs and 248 downregulated DEGs. Functional categories showed that many DEGs were associated with the ribosome, the insecticide resistance, the immune response and the digestion in comparison with CLas-infected SG and CLas-free SG. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases confirmed that metabolism and immunity were important functions in the SG. Among the DEGs, 68 genes (35 upregulated and 33 downregulated) encoding putative-secreted proteins were obtained with a signal peptide, suggesting that these genes may play important roles in CLas infection. A total of 673 SG proteins were identified in uninfected D. citri by liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) analysis, and 30 DEGs (15 upregulated and 15 downregulated) were found using the local tBLASTP programs. Among the 30 DEGs, many DEGs mainly involved in the metabolism and cellular processes pathways. This study provides basic transcriptome and proteome information for the SG in D. citri, and helps illuminate the molecular interactions between CLas and D. citri.

Brown marmorated stink bug, Halyomorpha halys (Stål), genome: putative underpinnings of polyphagy, insecticide resistance potential and biology of a top worldwide pest

BACKGROUND

Halyomorpha halys (Stål), the brown marmorated stink bug, is a highly invasive insect species due in part to its exceptionally high levels of polyphagy. This species is also a nuisance due to overwintering in human-made structures. It has caused significant agricultural losses in recent years along the Atlantic seaboard of North America and in continental Europe. Genomic resources will assist with determining the molecular basis for this species' feeding and habitat traits, defining potential targets for pest management strategies.

RESULTS

Analysis of the 1.15-Gb draft genome assembly has identified a wide variety of genetic elements underpinning the biological characteristics of this formidable pest species, encompassing the roles of sensory functions, digestion, immunity, detoxification and development, all of which likely support H. halys' capacity for invasiveness. Many of the genes identified herein have potential for biomolecular pesticide applications.

CONCLUSIONS

Availability of the H. halys genome sequence will be useful for the development of environmentally friendly biomolecular pesticides to be applied in concert with more traditional, synthetic chemical-based controls.

Phenotypic and transcriptomic responses of two Nilaparvata lugens populations to the Mudgo rice containing Bph1

The Bph1 gene was the first reported brown planthopper (BPH, Nilaparvata lugens) resistance gene in Mudgo rice and was widely used as a commercial cultivar for controlling BPH infestations. However, rapid adaptations of BPH on the Mudgo rice resulted in its resistance breakdown and the emergence of virulent BPH populations. Thus, specific BPH populations and rice varieties can serve as good model systems for studying the roles of different bio-compounds and proteins in the insect-plant interactions. Although our understandings have been improved on the complexity of BPH and rice interactions, the underlying molecular mechanisms remain largely unknown. Here we analyzed the feeding performances and the transcriptomic responses of two BPH populations (Mugdo-BPH and TN1-BPH) during compatible (Mudog-BPH feeding on Mudgo rice) and incompatible (TN1-BPH feeding on Mudgo rice) interactions. The electrical penetration graph (EPG) results indicated that the BPH feeding and performances during the incompatible interaction are significantly affected in terms of decreased honeydew, loss of weight, decreased phloem sap ingestion (N4 waveform), but increased non-penetration (NP waveform) phase. Abundance of glucose and trehalose was reduced in BPH during the incompatible interaction. Transcriptomic surveys of insects in both interactions revealed that genes involved in cuticle formation, detoxification, metabolite transport, digestion, RNA processing, lipid or fatty acid metabolism, and proteolysis were significantly down-regulated during the incompatible interaction, whereas genes involved in insulin signaling were significantly upregulated. Knockdown of four genes, including the sugar transporter NlST45, the serine and arginine-rich protein NlSRp54, the cytochrome P450 gene NlCYP6AY1, and the cuticle protein NlCPR70 through RNA-interference revealed thess genes are important for BPH survival. Overall, the results of this study will be helpful for the future researches on BPH virulence shifts.

Exposure to spinosad induces histopathological and cytotoxic effects on the salivary complex of the non-target predator Podisus nigrispinus

In integrated pest management systems, biological and chemical controls must be compatible. The insecticide spinosad affects some non-target insects and might compromise their fitness. The objective of this study was to evaluate the histopathological and cytotoxic effects of spinosad on the salivary complex of the predatory bug Podisus nigrispinus (Heteroptera: Pentatomidae). Spinosad toxicity and insect survival were determined using six concentrations of insecticide. Ultrastructural changes and cell death of salivary glands were analyzed after P. nigrispinus exposure to spinosad LC₅₀ (3.15 μg L^-1). The insecticide caused toxicity to P. nigrispinus; survival was 32% after 48 h of exposure to LC₅₀. The main histological changes in the salivary complex were disorganization of the epithelium, cytoplasmic vacuolization, and apocrine secretion into the gland lumen. Cytotoxic effects, such as release of granules and vacuoles into the lumen, presence of autophagosomes, enlargement of basal plasma membrane infoldings, and apoptosis, were observed. Spinosad causes toxicity, decreases survival, and changes the histology and cytology of the P. nigrispinus salivary complex. The results suggest that the cellular stress induced by the insecticide affects extra-oral digestion, compromising the potential of P. nigrispinus as a biological pest control agent.

Short-Term Physiological Response of a Native Hawaiian Plant, Hibiscus arnottianus, to Injury by the Exotic Leafhopper, Sophonia orientalis (Hemiptera: Cicadellidae)

Sophonia orientalis (Matsumura), also known as the two-spotted leafhopper, is a widespread exotic pest of many economically important crop plants and ornamental plants in Hawaii. Sophonia orientalis is highly polyphagous and is a major threat to some of the native endemic plants. Despite the successful establishment in Hawaii, interactions of S. orientalis with its host plants remain poorly understood. Previous studies primarily focused on distribution, parasitism, and oviposition of S. orientalis in Hawaii, whereas plant physiological responses to the leafhopper's injury, and, specifically, gas exchange rates in plants, have not yet been described. In this study, we examined a short-term physiological response of a native Hawaiian plant, Hibiscus arnottianus (A. Gray), to injury by S. orientalis. We also explored whether Camellia sinensis (L.) Kuntze, a native host plant of S. orientalis in Asia, exhibits a similar response. We found that H. arnottianus plants demonstrated a rapid (2-d) physiological response to injury accompanied by 40% reduction in rate of photosynthesis and 42% reduction in rate of transpiration, whereas C. sinensis did not exhibit any reduction in gas exchange rates. We did not record any changes in plant chlorophyll levels after plant injury in either species. Our results suggest that H. arnottianus is responding to the leafhopper feeding with a generalized wound response predicted for novel plant-insect herbivore associations. We discuss potential future directions for studies which might focus on host plant responses to S. orientalis in its native versus introduced range.

The Dual Prey-Inactivation Strategy of Spiders-In-Depth Venomic Analysis of Cupiennius salei

Most knowledge of spider venom concerns neurotoxins acting on ion channels, whereas proteins and their significance for the envenomation process are neglected. The here presented comprehensive analysis of the venom gland transcriptome and proteome of Cupiennius salei focusses on proteins and cysteine-containing peptides and offers new insight into the structure and function of spider venom, here described as the dual prey-inactivation strategy. After venom injection, many enzymes and proteins, dominated by α-amylase, angiotensin-converting enzyme, and cysteine-rich secretory proteins, interact with main metabolic pathways, leading to a major disturbance of the cellular homeostasis. Hyaluronidase and cytolytic peptides destroy tissue and membranes, thus supporting the spread of other venom compounds. We detected 81 transcripts of neurotoxins from 13 peptide families, whereof two families comprise 93.7% of all cysteine-containing peptides. This raises the question of the importance of the other low-expressed peptide families. The identification of a venom gland-specific defensin-like peptide and an aga-toxin-like peptide in the hemocytes offers an important clue on the recruitment and neofunctionalization of body proteins and peptides as the origin of toxins.

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.

De novo transcriptome assembly of the Chinese pearl barley, adlay, by full-length isoform and short-read RNA sequencing

Adlay (Coix lacryma-jobi) is a tropical grass that has long been used in traditional Chinese medicine and is known for its nutritional benefits. Recent studies have shown that vitamin E compounds in adlay protect against chronic diseases such as cancer and heart disease. However, the molecular basis of adlay's health benefits remains unknown. Here, we generated adlay gene sets by de novo transcriptome assembly using long-read isoform sequencing (Iso-Seq) and short-read RNA-Sequencing (RNA-Seq). The gene sets obtained from Iso-seq and RNA-seq contained 31,177 genes and 57,901 genes, respectively. We confirmed the validity of the assembled gene sets by experimentally analyzing the levels of prolamin and vitamin E biosynthesis-associated proteins in adlay plant tissues and seeds. We compared the screened adlay genes with known gene families from closely related plant species, such as rice, sorghum and maize. We also identified tissue-specific genes from the adlay leaf, root, and young and mature seed, and experimentally validated the differential expression of 12 randomly-selected genes. Our study of the adlay transcriptome will provide a valuable resource for genetic studies that can enhance adlay breeding programs in the future.

Transcriptome-wide responses of adult melon thrips (Thrips palmi) associated with capsicum chlorosis virus infection

Thrips palmi is a widely distributed major agricultural pest in the tropics and subtropics, causing significant losses in cucurbit and solanaceous crops through feeding damage and transmission of tospoviruses. Thrips palmi is a vector of capsicum chlorosis virus (CaCV) in Australia. The present understanding of transmission biology and potential effects of CaCV on T. palmi is limited. To gain insights into molecular responses to CaCV infection, we performed RNA-Seq to identify thrips transcripts that are differentially-abundant during virus infection of adults. De-novo assembly of the transcriptome generated from whole bodies of T. palmi adults generated 166,445 contigs, of which ~24% contained a predicted open reading frame. We identified 1,389 differentially-expressed (DE) transcripts, with comparable numbers up- (708) and down-regulated (681) in virus-exposed thrips compared to non-exposed thrips. Approximately 59% of these DE transcripts had significant matches to NCBI non-redundant proteins (Blastx) and Blast2GO identified provisional functional categories among the up-regulated transcripts in virus-exposed thrips including innate immune response-related genes, salivary gland and/or gut-associated genes and vitellogenin genes. The majority of the immune-related proteins are known to serve functions in lysosome activity and melanisation in insects. Most of the up-regulated oral and extra-oral digestion-associated genes appear to be involved in digestion of proteins, lipids and plant cell wall components which may indirectly enhance the likelihood or frequency of virus transmission or may be involved in the regulation of host defence responses. Most of the down-regulated transcripts fell into the gene ontology functional category of 'structural constituent of cuticle'. Comparison to DE genes responsive to tomato spotted wilt virus in Frankliniella occidentalis indicates conservation of some thrips molecular responses to infection by different tospoviruses. This study assembled the first transcriptome in the genus Thrips and provides important data to broaden our understanding of networks of molecular interactions between thrips and tospoviruses.

The green rice leafhopper, Nephotettix cincticeps (Hemiptera: Cicadellidae), salivary protein NcSP75 is a key effector for successful phloem ingestion

Nephotettix cincticeps, a prevalent rice pest, injects gelling and watery saliva into plant tissues during the sucking process. Certain components within the saliva are believed to interact with plant cellular constituents and play important roles in overcoming host plant defense responses. Based on our previous analysis of the salivary gland transcriptome and secreted saliva proteome of N. cincticeps, in this study, we analyzed the biological functions of salivary protein, NcSP75 (N. cincticepssalivary protein 75 kD). NcSP75, a salivary glands-specific gene, showed low similarities to any previously reported sequences. Knockdown of NcSP75 by RNA interference (RNAi) reduced the longevity of treated nymphs to approximately half of the longevity of controls and caused severe developmental retardation. Furthermore, the knockdown of NcSP75 decreased the survival rate of adults, and reduced the number of deposited eggs and hatched nymphs. Thus, the adverse effects caused by the knockdown of NcSP75 were observed throughout the lifetime of N. cincticeps, when feeding on rice plants. In contrast, no reduction was observed in the survival rate of the knockdown of NcSP75 adults fed on an artificial diet. Electrical penetration graph measurements taken from adult females feeding on rice plants showed a significantly shorter duration of phloem ingestion associated with the knockdown of NcSP75 than the knockdown of the enhanced green fluorescent protein (EGFP). Furthermore, the total sugar content of the honeydew was lower when NcSP75 was knocked down. These results suggest that the NcSP75 protein contribute to successful and sustainable ingestion from the sieve elements of rice plants. The NcSP75 protein of N. cincticeps can, accordingly, be considered as a key effector for establishing compatible interaction with rice plants and could be a potential target for controlling this species.

Two Phytoplasmas Elicit Different Responses in the Insect Vector Euscelidius variegatus Kirschbaum

Phytoplasmas are plant-pathogenic bacteria transmitted by hemipteran insects. The leafhopper Euscelidius variegatus is a natural vector of chrysanthemum yellows phytoplasma (CYp) and a laboratory vector of flavescence dorée phytoplasma (FDp). The two phytoplasmas induce different effects on this species: CYp slightly improves whereas FDp negatively affects insect fitness. To investigate the molecular bases of these different responses, transcriptome sequencing (RNA-seq) analysis of E. variegatus infected with either CYp or FDp was performed. The sequencing provided the first de novo transcriptome assembly for a phytoplasma vector and a starting point for further analyses on differentially regulated genes, mainly related to immune system and energy metabolism. Insect phenoloxidase activity, immunocompetence, and body pigmentation were measured to investigate the immune response, while respiration and movement rates were quantified to confirm the effects on energy metabolism. The activation of the insect immune response upon infection with FDp, which is not naturally transmitted by E. variegatus, confirmed that this bacterium is mostly perceived as a potential pathogen. Conversely, the acquisition of CYp, which is naturally transmitted by E. variegatus, seems to increase the insect fitness by inducing a prompt response to stress. This long-term relationship is likely to improve survival and dispersal of the infected insect, thus enhancing the opportunity of phytoplasma transmission.

Current state of knowledge on Wolbachia infection among Coleoptera: a systematic review

Background

Despite great progress in studies on Wolbachia infection in insects, the knowledge about its relations with beetle species, populations and individuals, and the effects of bacteria on these hosts, is still unsatisfactory. In this review we summarize the current state of knowledge about Wolbachia occurrence and interactions with Coleopteran hosts.

Methods

An intensive search of the available literature resulted in the selection of 86 publications that describe the relevant details about Wolbachia presence among beetles. These publications were then examined with respect to the distribution and taxonomy of infected hosts and diversity of Wolbachia found in beetles. Sequences of Wolbachia genes (16S rDNA, ftsZ) were used for the phylogenetic analyses.

Results

The collected publications revealed that Wolbachia has been confirmed in 204 beetle species and that the estimated average prevalence of this bacteria across beetle species is 38.3% and varies greatly across families and genera (0–88% infected members) and is much lower (c. 13%) in geographic studies. The majority of the examined and infected beetles were from Europe and East Asia. The most intensively studied have been two groups of herbivorous beetles: Curculionidae and Chrysomelidae. Coleoptera harbor Wolbachia belonging to three supergroups: F found in only three species, and A and B found in similar numbers of beetles (including some doubly infected); however the latter two were most prevalent in different families. A total of 59% of species with precise data were found to be totally infected. Single infections were found in 69% of species and others were doubly- or multiply-infected. Wolbachia caused numerous effects on its beetle hosts, including selective sweep with host mtDNA (found in 3% of species), cytoplasmic incompatibility (detected in c. 6% of beetles) and other effects related to reproduction or development (like male-killing, possible parthenogenesis or haplodiploidy induction, and egg development). Phylogenetic reconstructions for Wolbachia genes rejected cospeciation between these bacteria and Coleoptera, with minor exceptions found in some Hydraenidae, Curculionidae and Chrysomelidae. In contrast, horizontal transmission of bacteria has been suspected or proven in numerous cases (e.g., among beetles sharing habitats and/or host plants).

Discussion

The present knowledge about Wolbachia infection across beetle species and populations is very uneven. Even the basic data about infection status in species and frequency of infected species across genera and families is very superficial, as only c. 0.15% of all beetle species have been tested so far. Future studies on Wolbachia diversity in Coleoptera should still be based on the Multi-locus Sequence Typing system, and next-generation sequencing technologies will be important for uncovering Wolbachia relations with host evolution and ecology, as well as with other, co-occurring endosymbiotic bacteria.

Study of the Metatranscriptome of Eight Social and Solitary Wild Bee Species Reveals Novel Viruses and Bee Parasites

Bees are associated with a remarkable diversity of microorganisms, including unicellular parasites, bacteria, fungi, and viruses. The application of next-generation sequencing approaches enables the identification of this rich species composition as well as the discovery of previously unknown associations. Using high-throughput polyadenylated ribonucleic acid (RNA) sequencing, we investigated the metatranscriptome of eight wild bee species (Andrena cineraria, Andrena fulva, Andrena haemorrhoa, Bombus terrestris, Bombus cryptarum, Bombus pascuorum, Osmia bicornis, and Osmia cornuta) sampled from four different localities in Belgium. Across the RNA sequencing libraries, 88–99% of the taxonomically informative reads were of the host transcriptome. Four viruses with homology to insect pathogens were found including two RNA viruses (belonging to the families Iflaviridae and Tymoviridae that harbor already viruses of honey bees), a double stranded DNA virus (family Nudiviridae) and a single stranded DNA virus (family Parvoviridae). In addition, we found genomic sequences of 11 unclassified arthropod viruses (related to negeviruses, sobemoviruses, totiviruses, rhabdoviruses, and mononegaviruses), seven plant pathogenic viruses, and one fungal virus. Interestingly, nege-like viruses appear to be widespread, host-specific, and capable of attaining high copy numbers inside bees. Next to viruses, three novel parasite associations were discovered in wild bees, including Crithidia pragensis and a tubulinosematid and a neogregarine parasite. Yeasts of the genus Metschnikowia were identified in solitary bees. This study gives a glimpse of the microorganisms and viruses associated with social and solitary wild bees and demonstrates that their diversity exceeds by far the subset of species first discovered in honey bees.

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.

Deep sequencing and transcriptome analyses to identify genes involved in secoiridoid biosynthesis in the Tibetan medicinal plant Swertia mussotii

Swertia mussotii Franch. is an important traditional Tibetan medicinal plant with pharmacological properties effective in the treatment of various ailments including hepatitis. Secoiridoids are the major bioactive compounds in S. mussotii. To better understand the secoiridoid biosynthesis pathway, we generated transcriptome sequences from the root, leaf, stem, and flower tissues, and performed de novo sequence assembly, yielding 98,613 unique transcripts with an N50 of 1,085 bp. Putative functions could be assigned to 35,029 transcripts (35.52%) based on BLAST searches against annotation databases including GO and KEGG. The expression profiles of 39 candidate transcripts encoding the key enzymes for secoiridoid biosynthesis were examined in different S. mussotii tissues, validated by qRT-PCR, and compared with the homologous genes from S. japonica, a species in the same family, unveiling the gene expression, regulation, and conservation of the pathway. The examination of the accumulated levels of three bioactive compounds, sweroside, swertiamarin, and gentiopicroside, revealed their considerable variations in different tissues, with no significant correlation with the expression profiles of key genes in the pathway, suggesting complex biological behaviours in the coordination of metabolite biosynthesis and accumulation. The genomic dataset and analyses presented here lay the foundation for further research on this important medicinal plant.

Gene Duplication and Gene Expression Changes Play a Role in the Evolution of Candidate Pollen Feeding Genes in Heliconius Butterflies

Heliconius possess a unique ability among butterflies to feed on pollen. Pollen feeding significantly extends their lifespan, and is thought to have been important to the diversification of the genus. We used RNA sequencing to examine feeding-related gene expression in the mouthparts of four species of Heliconius and one nonpollen feeding species, Eueides isabella. We hypothesized that genes involved in morphology and protein metabolism might be upregulated in Heliconius because they have longer proboscides than Eueides, and because pollen contains more protein than nectar. Using de novo transcriptome assemblies, we tested these hypotheses by comparing gene expression in mouthparts against antennae and legs. We first looked for genes upregulated in mouthparts across all five species and discovered several hundred genes, many of which had functional annotations involving metabolism of proteins (cocoonase), lipids, and carbohydrates. We then looked specifically within Heliconius where we found eleven common upregulated genes with roles in morphology (CPR cuticle proteins), behavior (takeout-like), and metabolism (luciferase-like). Closer examination of these candidates revealed that cocoonase underwent several duplications along the lineage leading to heliconiine butterflies, including two Heliconius-specific duplications. Luciferase-like genes also underwent duplication within lepidopterans, and upregulation in Heliconius mouthparts. Reverse-transcription PCR confirmed that three cocoonases, a peptidase, and one luciferase-like gene are expressed in the proboscis with little to no expression in labial palps and salivary glands. Our results suggest pollen feeding, like other dietary specializations, was likely facilitated by adaptive expansions of preexisting genes—and that the butterfly proboscis is involved in digestive enzyme production.

Transcriptomic and Expression Analysis of the Salivary Glands in White-Backed Planthoppers, Sogatella furcifera

The white-backed planthopper (WBPH), Sogatella furcifera (Horváth), is one of the serious rice pests because of its destructive feeding. The salivary glands of the WBPH play an important role in the feeding behaviour. Currently, however, very little is known about the salivary glands at the molecular level. We sequenced the salivary gland transcriptome (sialotranscripome) of adult WBPHs using the Illumina sequencing. A total of 65,595 transcripts and 51,842 unigenes were obtained from salivary glands. According to annotations against the Nr database, many of the unigenes identified were associated with the most studied enzymes in hemipteran saliva. In the present study, we identified 32 salivary protein genes from the WBPH sialotranscripome, which were categorized as those involved in sugar metabolism, detoxification, suppression of plant defense responses, immunity-related responses, general digestion, and other phytophagy processes. Tissue expression profiles analysis revealed that four of 32 salivary protein genes (multicopper oxidase 4, multicopper oxidase 6, carboxylesterase and uridine phosphorylase 1 isform X2) were primarily expressed in the salivary gland, suggesting that they played putative role in insect-rice interactions. 13 of 32 salivary protein genes were primarily expressed in gut, which might play putative role in digestive and detoxify mechanism. Development expression profiles analysis revealed that the expression level of 26 of 32 salivary protein genes had no significant difference, suggesting that they may play roles in every developmental stages of salivary gland of WBPH. The other six genes have a high expression level in the salivary gland of adult. 31 of 32 genes (except putative acetylcholinesterase 1) have no significant difference in male and female adult, suggesting that their expression level have no difference between sexes. This report analysis of the sialotranscripome for the WBPH, and the transcriptome provides a foundational list of the genes involved in feeding. Our data will be useful to investigate the mechanisms of interaction between the WBPH and the host plant.

Comparative transcriptomics and proteomics of three different aphid species identifies core and diverse effector sets

BACKGROUND

Aphids are phloem-feeding insects that cause significant economic losses to agriculture worldwide. While feeding and probing these insects deliver molecules, called effectors, inside their host to enable infestation. The identification and characterization of these effectors from different species that vary in their host range is an important step in understanding the infestation success of aphids and aphid host range variation. This study employs a multi-disciplinary approach based on transcriptome sequencing and proteomics to identify and compare effector candidates from the broad host range aphid Myzus persicae (green peach aphid) (genotypes O, J and F), and narrow host range aphids Myzus cerasi (black cherry aphid) and Rhopalosiphum padi (bird-cherry oat aphid).

RESULTS

Using a combination of aphid transcriptome sequencing on libraries derived from head versus body tissues as well as saliva proteomics we were able to predict candidate effectors repertoires from the different aphid species and genotypes. Among the identified conserved or core effector sets, we identified a significant number of previously identified aphid candidate effectors indicating these proteins may be involved in general infestation strategies. Moreover, we identified aphid candidate effector sequences that were specific to one species, which are interesting candidates for further validation and characterization with regards to species-specific functions during infestation. We assessed our candidate effector repertoires for evidence of positive selection, and identified 49 candidates with DN/DS ratios >1. We noted higher rates of DN/DS ratios in predicted aphid effectors than non-effectors. Whether this reflects positive selection due to co-evolution with host plants, or increased neofunctionalization upon gene duplication remains to be investigated.

CONCLUSION

Our work provides a comprehensive overview of the candidate effector repertoires from three different aphid species with varying host ranges. Comparative analyses revealed candidate effectors that are most likely are involved in general aspects of infestation, whereas others, that are highly divergent, may be involved in specific processes important for certain aphid species. Insights into the overlap and differences in aphid effector repertoires are important in understanding how different species successfully infest different ranges of plant species.

Organ-specific transcriptome response of the small brown planthopper toward rice stripe virus

Rice stripe virus (RSV) causes rice stripe disease and is transmitted by the small brown planthopper (Laodelphax striatellus, SBPH) in a persistent, circulative, and propagative manner. The alimentary canal and salivary gland of SBPH play important roles in viral replication and transmission. However, little is known about the underlying molecular functions of these two organs in the interaction between RSV and SBPH. In this study, organ-specific transcriptomes of the alimentary canal and salivary gland were analyzed in viruliferous and naïve SBPH. The number of differentially expressed unigenes in the alimentary canal was considerably greater than that in the salivary gland after RSV infection, and only 23 unigenes were co-regulated in the two organs. In the alimentary canal, genes involved in lysosome, digestion and detoxification were activated upon RSV infection, whereas the genes related to DNA replication and repair were suppressed. RSV activated RNA transport and repressed the MAPK, mTOR, Wnt, and TGF-beta signaling pathways in the salivary gland. The overall immune reaction toward RSV was much stronger in the salivary gland than in the alimentary canal. RSV activated the pattern recognition molecules and Toll pathway in the salivary gland but inhibited these two reactions in the alimentary canal. The responses from reactive oxygen and the immune-responsive effectors were stronger in the salivary gland than in the alimentary canal after RSV infection. These findings provide clues on the roles of the two organs in confronting RSV infection and aid in the understanding of the interaction between RSV and SBPHs.

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

The tarnished plant bug (TPB), Lygus lineolaris (Palisot de Beauvois) is a polyphagous, phytophagous insect that has emerged as a major pest of cotton, alfalfa, fruits, and vegetable crops in the eastern United States and Canada. Using its piercing-sucking mouthparts, TPB employs a “lacerate and flush” feeding strategy in which saliva injected into plant tissue degrades cell wall components and lyses cells whose contents are subsequently imbibed by the TPB. It is known that a major component of TPB saliva is the polygalacturonase enzymes that degrade the pectin in the cell walls. However, not much is known about the other components of the saliva of this important pest. In this study, we explored the salivary gland transcriptome of TPB using Illumina sequencing. After in silico conversion of RNA sequences into corresponding polypeptides, 25,767 putative proteins were discovered. Of these, 19,540 (78.83%) showed significant similarity to known proteins in the either the NCBI nr or Uniprot databases. Gene ontology (GO) terms were assigned to 7,512 proteins, and 791 proteins in the sialotranscriptome of TPB were found to collectively map to 107 Kyoto Encyclopedia of Genes and Genomes (KEGG) database pathways. A total of 3,653 Pfam domains were identified in 10,421 sialotranscriptome predicted proteins resulting in 12,814 Pfam annotations; some proteins had more than one Pfam domain. Functional annotation revealed a number of salivary gland proteins that potentially facilitate degradation of host plant tissues and mitigation of the host plant defense response. These transcripts/proteins and their potential roles in TPB establishment are described.

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.

Morphology and Proteome Characterization of the Salivary Glands of the Western Chinch Bug (Hemiptera: Blissidae)

The western chinch bug, Blissus occiduus Barber, is a serious pest of buffalograss, Buchloe dactyloides (Nuttall) due to physical and chemical damage caused during the feeding process. Although previous work has investigated the feeding behaviors of chinch bugs in the Blissus complex, no study to date has explored salivary gland morphology and the associated salivary complex of this insect. Whole and sectioned B. occiduus salivary glands were visualized using light and scanning electron microscopy to determine overall structure and cell types of the salivary glands and their individual lobes. Microscopy revealed a pair of trilobed principal glands and a pair of tubular accessory glands of differing cellular types. To link structure with function, the salivary gland proteome was characterized using liquid chromatography tandem mass spectrometry. The salivary proteome analysis resulted in B. occiduus sequences matching 228 nonhomologous protein sequences of the pea aphid, Acyrthosiphon pisum (Harris), with many specific to the proteins present in the salivary proteome of A. pisum. A number of sequences were assigned the molecular function of hydrolase and oxido-reductase activity, with one specific protein sequence revealing a peroxidase-like function. This is the first study to characterize the salivary proteome of B. occiduus and the first of any species in the family Blissidae.

Mechanisms and ecological consequences of plant defence induction and suppression in herbivore communities

BACKGROUND

Plants are hotbeds for parasites such as arthropod herbivores, which acquire nutrients and energy from their hosts in order to grow and reproduce. Hence plants are selected to evolve resistance, which in turn selects for herbivores that can cope with this resistance. To preserve their fitness when attacked by herbivores, plants can employ complex strategies that include reallocation of resources and the production of defensive metabolites and structures. Plant defences can be either prefabricated or be produced only upon attack. Those that are ready-made are referred to as constitutive defences. Some constitutive defences are operational at any time while others require activation. Defences produced only when herbivores are present are referred to as induced defences. These can be established via de novo biosynthesis of defensive substances or via modifications of prefabricated substances and consequently these are active only when needed. Inducibility of defence may serve to save energy and to prevent self-intoxication but also implies that there is a delay in these defences becoming operational. Induced defences can be characterized by alterations in plant morphology and molecular chemistry and are associated with a decrease in herbivore performance. These alterations are set in motion by signals generated by herbivores. Finally, a subset of induced metabolites are released into the air as volatiles and function as a beacon for foraging natural enemies searching for prey, and this is referred to as induced indirect defence.

SCOPE

The objective of this review is to evaluate (1) which strategies plants have evolved to cope with herbivores and (2) which traits herbivores have evolved that enable them to counter these defences. The primary focus is on the induction and suppression of plant defences and the review outlines how the palette of traits that determine induction/suppression of, and resistance/susceptibility of herbivores to, plant defences can give rise to exploitative competition and facilitation within ecological communities "inhabiting" a plant.

CONCLUSIONS

Herbivores have evolved diverse strategies, which are not mutually exclusive, to decrease the negative effects of plant defences in order to maximize the conversion of plant material into offspring. Numerous adaptations have been found in herbivores, enabling them to dismantle or bypass defensive barriers, to avoid tissues with relatively high levels of defensive chemicals or to metabolize these chemicals once ingested. In addition, some herbivores interfere with the onset or completion of induced plant defences, resulting in the plant's resistance being partly or fully suppressed. The ability to suppress induced plant defences appears to occur across plant parasites from different kingdoms, including herbivorous arthropods, and there is remarkable diversity in suppression mechanisms. Suppression may strongly affect the structure of the food web, because the ability to suppress the activation of defences of a communal host may facilitate competitors, whereas the ability of a herbivore to cope with activated plant defences will not. Further characterization of the mechanisms and traits that give rise to suppression of plant defences will enable us to determine their role in shaping direct and indirect interactions in food webs and the extent to which these determine the coexistence and persistence of species.

Differentially regulated genes in the salivary glands of brown planthopper after feeding in resistant versus susceptible rice varieties

Brown planthopper (BPH) is a damaging insect pest of rice. We used suppression subtractive hybridization (SSH) and mirror orientation selection to identify differentially regulated genes in salivary glands of BPH after feeding on resistant and susceptible varieties. The forward SSH library included 768 clones with insertions ranging from 250 to 1000 bp. After differential screening, a total of 112 transcripts were identified, which included 27 upregulated genes and seven downregulated genes. Several of these transcripts showed sequence homology to known proteins such as trehalase, mucin-like protein, vitellogenin, calcium ion binding protein, and eukaryotic initiation factor-like protein. About half of the transcripts, however, did not match to any sequences in the protein databases currently available. Functional annotation of the transcripts showed gene ontology association with metabolism, signal transduction, and regulatory responses. Notably, many known functional genes were predicted to be secreted proteins. Also, gene expression profiles of the salivary glands of BPH feeding on resistant rice (B5) and susceptible rice (TN1) varieties were compared. Our data provide a molecular resource for future functional studies on salivary glands and will be useful for elucidating the molecular mechanisms between BPH feeding and rice varieties with BPH resistance differences.

Proteome Analysis of Watery Saliva Secreted by Green Rice Leafhopper, Nephotettix cincticeps

The green rice leafhopper, Nephotettix cincticeps, is a vascular bundle feeder that discharges watery and gelling saliva during the feeding process. To understand the potential functions of saliva for successful and safe feeding on host plants, we analyzed the complexity of proteinaceous components in the watery saliva of N. cincticeps. Salivary proteins were collected from a sucrose diet that adult leafhoppers had fed on through a membrane of stretched parafilm. Protein concentrates were separated using SDS-PAGE under reducing and non-reducing conditions. Six proteins were identified by a gas-phase protein sequencer and two proteins were identified using LC-MS/MS analysis with reference to expressed sequence tag (EST) databases of this species. Full -length cDNAs encoding these major proteins were obtained by rapid amplification of cDNA ends-PCR (RACE-PCR) and degenerate PCR. Furthermore, gel-free proteome analysis that was performed to cover the broad range of salivary proteins with reference to the latest RNA-sequencing data from the salivary gland of N. cincticeps, yielded 63 additional protein species. Out of 71 novel proteins identified from the watery saliva, about 60 % of those were enzymes or other functional proteins, including GH5 cellulase, transferrin, carbonic anhydrases, aminopeptidase, regucalcin, and apolipoprotein. The remaining proteins appeared to be unique and species- specific. This is the first study to identify and characterize the proteins in watery saliva of Auchenorrhyncha species, especially sheath-producing, vascular bundle-feeders.

Differential expression of candidate salivary effector proteins in field collections of Hessian fly, Mayetiola destructor

Evidence is emerging that some proteins secreted by gall-forming parasites of plants act as effectors responsible for systemic changes in the host plant, such as galling and nutrient tissue formation. A large number of secreted salivary gland proteins (SSGPs) that are the putative effectors responsible for the physiological changes elicited in susceptible seedling wheat by Hessian fly, Mayetiola destructor (Say), larvae have been documented. However, how the genes encoding these candidate effectors might respond under field conditions is unknown. The goal of this study was to use microarray analysis to investigate variation in SSGP transcript abundance amongst field collections from different geographical regions (southeastern USA, central USA, and the Middle East). Results revealed significant variation in SSGP transcript abundance amongst the field collections studied. The field collections separated into three distinct groups that corresponded to the wheat classes grown in the different geographical regions as well as to recently described Hessian fly populations. These data support previous reports correlating Hessian fly population structure with micropopulation differences owing to agro-ecosystem parameters such as cultivation of regionally adapted wheat varieties, deployment of resistance genes and variation in climatic conditions.

Transcriptome profiling of the testis reveals genes involved in spermatogenesis and marker discovery in the oriental fruit fly, Bactrocera dorsalis

The testis is a highly specialized tissue that plays a vital role in ensuring fertility by producing spermatozoa, which are transferred to the female during mating. Spermatogenesis is a complex process, resulting in the production of mature sperm, and involves significant structural and biochemical changes in the seminiferous epithelium of the adult testis. The identification of genes involved in spermatogenesis of Bactrocera dorsalis (Hendel) is critical for a better understanding of its reproductive development. In this study, we constructed a cDNA library of testes from male B. dorsalis adults at different ages, and performed de novo transcriptome sequencing to produce a comprehensive transcript data set, using Illumina sequencing technology. The analysis yielded 52 016 732 clean reads, including a total of 4.65 Gb of nucleotides. These reads were assembled into 47 677 contigs (average 443 bp) and then clustered into 30 516 unigenes (average 756 bp). Based on BLAST hits with known proteins in different databases, 20 921 unigenes were annotated with a cut-off E-value of 10(-5). The transcriptome sequences were further annotated using the Clusters of Orthologous Groups, Gene Orthology and the Kyoto Encyclopedia of Genes and Genomes databases. Functional genes involved in spermatogenesis were analysed, including cell cycle proteins, metalloproteins, actin, and ubiquitin and antihyperthermia proteins. Several testis-specific genes were also identified. The transcripts database will help us to understand the molecular mechanisms underlying spermatogenesis in B. dorsalis. Furthermore, 2913 simple sequence repeats and 151 431 single nucleotide polymorphisms were identified, which will be useful for investigating the genetic diversity of B. dorsalis in the future.

Transcriptome analysis of the salivary glands of Nephotettix cincticeps (Uhler)

The green rice leafhopper (GRH), Nephotettix cincticeps, is one of the most important pests of rice in temperate Asian countries. GRH, a vascular feeder, secretes watery and gelling saliva in the process of feeding on phloem and xylem sap. It is known that GRH saliva contains several bioactive proteins, including enzymes such as laccase and beta-glucosidase. In this study, we performed transcriptome analysis of salivary glands of GRH using Illumina paired-end sequencing. Of 51,788 assembled contigs, 16,017 (30.9%) showed significant similarity to known proteins in the NCBI nr database, while 34,978 (67.5%) could not be annotated by similarity search, Pfam, or gene ontology (GO). Contigs (905) with predicted signal peptides and no putative transmembrane domains are suggested to represent secreted protein coding genes. Among the 76 most highly expressed putative secretory protein contigs, 68 transcripts were found to be salivary gland-specific or at least -dominant, but not expressed in stomach or Malpighian tubules. However, 45 of the 68 transcripts were unknown proteins. These findings suggest that most of the GRH transcripts encoding secreted proteins expressed in salivary glands are species and/or tissue specific. Our results provide a fundamental list of genes involved in GRH-Poaceae host plant interactions including successful feeding and plant pathogen transmission.

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.

Salivary proteins of plant-feeding hemipteroids - implication in phytophagy

Many hemipteroids are major pests and vectors of microbial pathogens, infecting crops. Saliva of the hemipteroids is critical in enabling them to be voracious feeders on plants, including the economically important ones. A plethora of hemipteroid salivary enzymes is known to inflict stress in plants, either by degrading the plant tissue or by affecting their normal metabolism. Hemipteroids utilize one of the following three strategies of feeding behaviour: salivary sheath feeding, osmotic-pump feeding and cell-rupture feeding. The last strategy also includes several different tactics such as lacerate-and-flush, lacerate-and-sip and macerate-and-flush. Understanding hemipteroid feeding mechanisms is critical, since feeding behaviour directs salivary composition. Saliva of the Heteroptera that are specialized as fruit and seed feeders, includes cell-degrading enzymes, auchenorrhynchan salivary composition also predominantly consists of cell-degrading enzymes such as amylase and protease, whereas that of the Sternorhyncha includes a variety of allelochemical-detoxifying enzymes. Little is known about the salivary composition of the Thysanoptera. Cell-degrading proteins such as amylase, pectinase, cellulase and pectinesterase enable stylet entry into the plant tissue. In contrast, enzymes such as glutathione peroxidase, laccase and trehalase detoxify plant chemicals, enabling the circumvention of plant-defence mechanisms. Salivary enzymes such as M1-zinc metalloprotease and CLIP-domain serine protease as in Acyrthosiphon pisum (Aphididae), and non-enzymatic proteins such as apolipophorin, ficolin-3-like protein and 'lava-lamp' protein as in Diuraphis noxia (Aphididae) have the capacity to alter host-plant-defence mechanisms. A majority of the hemipteroids feed on phloem, hence Ca++-binding proteins such as C002 protein, calreticulin-like isoform 1 and calmodulin (critical for preventing sieve-plate occlusion) are increasingly being recognized in hemipteroid-plant interactions. Determination of a staggering variety of proteins shows the complexity of hemipteroid saliva: effector proteins localized in hemipteran saliva suggest a similarity to the physiology of pathogen-plant interactions.

Comparative transcriptome analysis of salivary glands of two populations of rice brown planthopper, Nilaparvata lugens, that differ in virulence

Background

The brown planthopper (BPH), Nilaparvata lugens (Stål), a destructive rice pest in Asia, can quickly overcome rice resistance by evolving new virulent populations. Herbivore saliva plays an important role in plant–herbivore interactions, including in plant defense and herbivore virulence. However, thus far little is known about BPH saliva at the molecular level, especially its role in virulence and BPH–rice interaction.

Methodology/Principal Findings

Using cDNA amplification in combination with Illumina short-read sequencing technology, we sequenced the salivary-gland transcriptomes of two BPH populations with different virulence; the populations were derived from rice variety TN1 (TN1 population) and Mudgo (M population). In total, 37,666 and 38,451 unigenes were generated from the salivary glands of these populations, respectively. When combined, a total of 43,312 unigenes were obtained, about 18 times more than the number of expressed sequence tags previously identified from these glands. Gene ontology annotations and KEGG orthology classifications indicated that genes related to metabolism, binding and transport were significantly active in the salivary glands. A total of 352 genes were predicted to encode secretory proteins, and some might play important roles in BPH feeding and BPH–rice interactions. Comparative analysis of the transcriptomes of the two populations revealed that the genes related to ‘metabolism,’ ‘digestion and absorption,’ and ‘salivary secretion’ might be associated with virulence. Moreover, 67 genes encoding putative secreted proteins were differentially expressed between the two populations, suggesting these genes may contribute to the change in virulence.

Conclusions/Significance

This study was the first to compare the salivary-gland transcriptomes of two BPH populations having different virulence traits and to find genes that may be related to this difference. Our data provide a rich molecular resource for future functional studies on salivary glands and will be useful for elucidating the molecular mechanisms underlying BPH feeding and virulence differences.

How phloem-feeding insects face the challenge of phloem-located defenses

Due to the high content of nutrient, sieve tubes are a primary target for pests, e.g., most phytophagous hemipteran. To protect the integrity of the sieve tubes as well as their content, plants possess diverse chemical and physical defense mechanisms. The latter mechanisms are important because they can potentially interfere with the food source accession of phloem-feeding insects. Physical defense mechanisms are based on callose as well as on proteins and often plug the sieve tube. Insects that feed from sieve tubes are potentially able to overwhelm these defense mechanisms using their saliva. Gel saliva forms a sheath in the apoplast around the stylet and is suggested to seal the stylet penetration site in the cell plasma membrane. In addition, watery saliva is secreted into penetrated cells including sieve elements; the presence of specific enzymes/effectors in this saliva is thought to interfere with plant defense responses. Here we detail several aspects of plant defense and discuss the interaction of plants and phloem-feeding insects. Recent agro-biotechnological phloem-located aphid control strategies are presented.