Insights into the transcriptomics of polyphagy: Bemisia tabaci adaptability to phenylpropanoids involves coordinated expression of defense and metabolic genes

Integrative transcriptomics reveals association of abscisic acid and lignin pathways with cassava whitefly resistance

BACKGROUND

Whiteflies are a global threat to crop yields, including the African subsistence crop cassava (Manihot esculenta). Outbreaks of superabundant whitefly populations throughout Eastern and Central Africa in recent years have dramatically increased the pressures of whitefly feeding and virus transmission on cassava. Whitefly-transmitted viral diseases threaten the food security of hundreds of millions of African farmers, highlighting the need for developing and deploying whitefly-resistant cassava. However, plant resistance to whiteflies remains largely poorly characterized at the genetic and molecular levels. Knowledge of cassava-defense programs also remains incomplete, limiting characterization of whitefly-resistance mechanisms. To better understand the genetic basis of whitefly resistance in cassava, we define the defense hormone- and Aleurotrachelus socialis (whitefly)-responsive transcriptome of whitefly-susceptible (COL2246) and whitefly-resistant (ECU72) cassava using RNA-seq. For broader comparison, hormone-responsive transcriptomes of Arabidopsis thaliana were also generated.

RESULTS

Whitefly infestation, salicylic acid (SA), jasmonic acid (JA), ethylene (ET), and abscisic acid (ABA) transcriptome responses of ECU72 and COL2246 were defined and analyzed. Strikingly, SA responses were largely reciprocal between the two cassava genotypes and we suggest candidate regulators. While susceptibility was associated with SA in COL2246, resistance to whitefly in ECU72 was associated with ABA, with SA-ABA antagonism observed. This was evidenced by expression of genes within the SA and ABA pathways and hormone levels during A. socialis infestation. Gene-enrichment analyses of whitefly- and hormone-responsive genes suggest the importance of fast-acting cell wall defenses (e.g., elicitor recognition, lignin biosynthesis) during early infestation stages in whitefly-resistant ECU72. A surge of ineffective immune and SA responses characterized the whitefly-susceptible COL2246's response to late-stage nymphs. Lastly, in comparison with the model plant Arabidopsis, cassava's hormone-responsive genes showed striking divergence in expression.

CONCLUSIONS

This study provides the first characterization of cassava's global transcriptome responses to whitefly infestation and defense hormone treatment. Our analyses of ECU72 and COL2246 uncovered possible whitefly resistance/susceptibility mechanisms in cassava. Comparative analysis of cassava and Arabidopsis demonstrated that defense programs in Arabidopsis may not always mirror those in crop species. More broadly, our hormone-responsive transcriptomes will also provide a baseline for the cassava community to better understand global responses to other yield-limiting pests/pathogens.

Quantity versus quality: Effects of diet protein-carbohydrate ratios and amounts on insect herbivore gene expression

Dietary protein and digestible carbohydrates are two key macronutrients for insect herbivores, but the amounts and ratios of these two macronutrients in plant vegetative tissues can be highly variable. Typically, insect herbivores regulate their protein-carbohydrate intake by feeding selectively on nutritionally complementary plant tissues, but this may not always be possible. Interestingly, lab experiments consistently demonstrate that performance - especially growth and survival - does not vary greatly when caterpillars and nymphal grasshoppers are reared on diets that differ in their protein-carbohydrate content. This suggests insect herbivores employ post-ingestive physiological mechanisms to compensate for variation in diet protein-carbohydrate profile. However, the molecular mechanisms that underlie this compensation are not well understood. Here we explore, for the first time in an insect herbivore, the transcriptional effects of two dietary factors: protein-to-carbohydrate ratio (p:c) and total macronutrient (p + c) content. Specifically, we reared Helicoverpa zea caterpillars on three diets that varied in diet p:c ratio and one diet that varied in total p + c concentration, all within an ecologically-relevant range. We observed two key findings. Caterpillars reared on diets with elevated total p + c content showed large differences in gene expression. In contrast, only small differences in gene expression were observed when caterpillars were reared on diets with different p:c ratios (spanning from protein-biased to carbohydrate-biased). The invariable expression of many metabolic genes across these variable diets suggests that H. zea caterpillars employ a strategy of constitutive expression to deal with protein-carbohydrate imbalances rather than diet-specific changes. This is further supported by two findings. First, few genes were uniquely associated with feeding on a protein- and carbohydrate-biased diet. Second, many differentially-expressed genes were shared across protein-biased, carbohydrate-biased, and concentrated diet treatments. Our study provides insights into the post-ingestive physiological mechanisms insect herbivores employ to regulate protein-carbohydrate intake. Most notably, it suggests that H. zea, and perhaps other generalist species, use similar post-ingestive mechanisms to deal with protein-carbohydrate imbalances - regardless of the direction of the imbalance.

Transcriptomic analysis of interactions between Lymantria dispar larvae and carvacrol

Due to its biological activity, carvacrol (CAR) is widely used in medicine, agriculture, and forestry. Our previous studies showed that in Lymantria dispar larvae, CAR treatment can induce the production of antifeedants and lead to growth inhibition and death of larvae. However, the effect CAR exerts on RNA levels in L. dispar larvae remains unclear. In this study, the Illumina HiSeq4000 sequencing platform was used to sequence the total RNA of L. dispar larvae. A total of six cDNA libraries (three treatments and three controls) were established and 39,807 genes were generated. Compared with the control group, 296 differentially expressed genes (DEGs) (142 up-regulated and 154 down-regulated) were identified after CAR treatment. GO and KEGG enrichment analyses showed that these DEGs mainly clustered in the metabolism of xenobiotics, carbohydrates, and lipids. Furthermore, 12 DEGs were found to be involved in detoxification, including six cytochrome P450s, two esterases, one glutathione peroxidase, one UDP-glycosyltransferase gene, and two genes encoding heat shock proteins. The expression levels of detoxification genes changed under CAR treatment (especially P450s), which further yielded candidate genes for explorations of the insecticidal mechanism of CAR. The reliability of transcriptome data was verified by qRT-PCR. The enzyme activities of CYP450 and acid phosphatase significantly increased (by 38.52 U/mg·prot and 0.12 μmol/min·mg, respectively) 72 h after CAR treatment. However, the activity of alkaline phosphatase did not change significantly. These changes in enzyme activity corroborated the reliability of the transcriptome data at the protein level. The results of GO enrichment analysis of DEGs indicated that CAR influenced the oxidation-reduction process in L. dispar larvae. Furthermore, CAR can cause oxidative stress in L. dispar larvae, identified through the determination of peroxidase and polyphenol oxidase activities, total antioxidant capacity, and hydrogen peroxide content. This study provides useful insight into the insecticidal mechanism of CAR.

The Chorion Proteome of Diaphorina citri, the Vector of Huanglongbing Disease in Citrus

Nowadays, the Asian citrus psyllid, Diaphorina citri (Kuwayama) (Hemiptera: Liviidae) is considered the most devastating pest of citrus because it transmits "Candidatus Liberibacter asiaticus", the putative causal agent of huanglongbing (HLB) or citrus greening. Controlling the vector is the main strategy used to mitigate HLB. Targeting D. citri at the very early stages of its development may offer an effective control strategy. Identifying chorion proteins will contribute to a better understanding of embryo development and egg hatching and thus could lead to valuable targets to better control psyllid populations. Herein, we analyze the chorion proteins of D. citri. Mass spectrometry-based bottom-up/shotgun proteomics and databases were queried to achieve protein identification. Fifty-one proteins were identified in D. citri chorion. The D. citri chorion proteins were divided into eight categories according to their biological or molecular function: i-enzymes (25%); ii-binding proteins (10%); iii-structural proteins (8%); iv-homeostasis-related proteins, mostly vitellogenins (8%); v-proteins related to gene expression (6%); vi-immune system proteins (6%); vii-other proteins (16%); and viii-uncharacterized proteins (21%). The composition of the chorion proteome suggested that the hatching rate could be reduced by silencing chorion-related genes. The proteomic analysis of D. citri chorion tissue allowed us to identify its proteins, providing promising new targets for D. citri control through RNA interference technology.

Host plant diet affects growth and induces altered gene expression and microbiome composition in the wood white (Leptidea sinapis) butterfly

In a time with decreasing biodiversity, especially among insects, a detailed understanding about specific resource utilization strategies is crucial. The physiological and behavioural responses to host switches in phytophagous insects are poorly understood. Earlier studies indicate that a host plant switch might be associated with distinctive molecular and physiological responses in different lineages. Expanding the assessment of such associations across Lepidoptera will reveal if there are general patterns in adaptive responses, or if each switch event is more of a unique character. We investigated host plant preference, fitness consequences, effects on expression profiles and gut microbiome composition in two common wood white (Leptidea sinapis) populations with different host plant preferences from the extremes of the species distribution area (Sweden and Catalonia). Our results show that female Catalonian wood whites lack preference for either host plant (Lotus corniculatus or L. dorycnium), while Swedish females laid significantly more eggs on L. corniculatus. Individuals from both populations reared on L. dorycnium had longer developmental times and smaller body size as adults. This indicates that both environmental and genetic factors determine the choice to use a specific host plant. Gene expression analysis revealed a more pronounced response to host plant in the Catalonian compared to the Swedish population. In addition, host plant treatment resulted in a significant shift in microbiome community structure in the Catalonian population. Together, this suggests that population specific plasticity associated with local conditions underlies host plant utilisation in wood whites.

Comparative analysis of diet-associated responses in two rice planthopper species

Background

Host adaptation is the primary determinant of insect diversification. However, knowledge of different host ranges in closely related species remains scarce. The brown planthopper (Nilaparvata lugens, BPH) and the small brown planthopper (Laodelphax striatellus, SBPH) are the most destructive insect pests within the family Delphacidae. These two species differ in their host range (SBPH can well colonize rice and wheat plants, whereas BPH survives on only rice plants), but the underlying mechanism of this difference remains unknown. High-throughput sequencing provides a powerful approach for analyzing the association between changes in gene expression and the physiological responses of insects. Therefore, gut transcriptomes were performed to elucidate the genes associated with host adaptation in planthoppers. The comparative analysis of planthopper responses to different diets will improve our knowledge of host adaptation regarding herbivorous insects.

Results

In the present study, we analyzed the change in gene expression of SBPHs that were transferred from rice plants to wheat plants over the short term (rSBPH vs tSBPH) or were colonized on wheat plants over the long term (rSBPH vs wSBPH). The results showed that the majority of differentially expressed genes in SBPH showed similar changes in expression for short-term transfer and long-term colonization. Based on a comparative analysis of BPH and SBPH after transfer, the genes associated with sugar transporters and heat-shock proteins showed similar variation. However, most of the genes were differentially regulated between the two species. The detoxification-related genes were upregulated in SBPH after transfer from the rice plants to the wheat plants, but these genes were downregulated in BPH under the same conditions. In contrast, ribosomal-related genes were downregulated in SBPH after transfer, but these genes were upregulated in BPH under the same conditions.

Conclusions

The results of this study provide evidence that host plants played a dominant role in shaping gene expression and that the low fitness of BPH on wheat plants might be determined within 24 h after transfer. This study deepens our understanding of different host ranges for the two planthopper species, which may provide a potential strategy for pest management.

Molecular Evolution of the Glutathione S-Transferase Family in the Bemisia tabaci Species Complex

The glutathione S-transferase (GST) family plays an important role in the adaptation of herbivorous insects to new host plants and other environmental constrains. The family codes for enzymes that neutralize reactive oxygen species and phytotoxins through the conjugation of reduced glutathione. Here, we studied the molecular evolution of the GST family in Bemisia tabaci, a complex of >35 sibling species, differing in their geographic and host ranges. We tested if some enzymes evolved different functionality, by comparing their sequences in six species, representing five of the six major genetic clades in the complex. Comparisons of the nonsynonymous to synonymous substitution ratios detected positive selection events in 11 codons of 5 cytosolic GSTs. Ten of them are located in the periphery of the GST dimer, suggesting a putative involvement in interactions with other proteins. Modeling the tertiary structure of orthologous enzymes, identified additional 19 mutations in 9 GSTs, likely affecting the enzymes' functionality. Most of the mutation events were found in the environmentally responsive classes Delta and Sigma, indicating a slightly different delta/sigma tool box in each species. At a broader genomic perspective, our analyses indicated a significant expansion of the Delta GST class in B. tabaci and a general association between the diet breadth of hemipteran species and their total number of GST genes. We raise the possibility that at least some of the identified changes improve the fitness of the B. tabaci species carrying them, leading to their better adaptation to specific environments.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

Gossypol-induced fitness gain and increased resistance to deltamethrin in beet armyworm, Spodoptera exigua (Hübner)

BACKGROUND

In plant-insect interactions, phytotoxins such as gossypol, exert a defensive role on behalf of the plant by interfering with the essential metabolic, biochemical and physiological pathways of herbivorous insects. The beet armyworm, Spodoptera exigua (Hübner), is a key pest for many important crops including a range of vegetables, ornamentals, and cotton. In this study, we investigated how resistance to deltamethrin relates to enzyme activity in gossypol-pretreated larvae of S. exigua.

RESULTS

Following selection with deltamethrin insecticides on gossypol-pretreated larvae for 10 generations, the Gos-SEL population developed a 113.29-fold resistance. Under the same conditions, the Delta-SEL selected population showed a 69.76-fold increase in resistance along with corresponding levels of xenobiotic defense enzyme activity. Similarly, the fecundity of the Delta-SEL population together with male and female longevity were found to be significantly lower when compared with the Gos-SEL population and the laboratory susceptible-strain group (SS-Strain). In addition, the activities of cytochrome P450s in S. exigua were significantly enhanced when the insects were fed on a deltamethrin and gossypol-pretreated diet compared with being fed on deltamethrin alone.

CONCLUSION

The reproductive capacity of S. exigua is significantly reduced in Delta-SEL and Gos-SEL populations compared with the control group (SS-Strain). Elevation of the major detoxification enzyme cytochrome P450 monooxygenase and esterase might have an important role in inducing tolerance to deltamethrin in gossypol-fed S. exigua populations. This study enhances our understanding of detoxification enzyme pathways for S. exigua gene expression and their role in responses to insecticides and plant secondary metabolites. © 2018 Society of Chemical Industry.

Enhanced effects of dietary tannic acid with chlorantraniliprole on life table parameters and nutritional physiology of Spodoptera exigua (Hübner)

The beet armyworm, Spodoptera exigua (Hübner) (Lepidoptera: Noctuidae) is a highly polyphagous pest which causes considerable economic losses to cotton and many vegetable crops. Tannins are among the most important secondary metabolites in cotton plants. We show that tannic acid enhances the toxic effect of chlorantraniliprole on S. exigua when presented in combination. Bioassays using third-instar S. exigua larvae on an artificial diet showed that consumption of tannic acid with chlorantraniliprole at the concentration of (2 mg/g and LC₅₀ 0.018 mg/L) had higher toxicity when compared to either chlorantraniliprole or tannic acid alone (LC₅₀ 0.027 mg/L). The diet containing tannic acid with chlorantraniliprole significantly prolonged larval and pupal developmental time and extended mean generation time and total pre-oviposition period compared to either chemical alone. Moreover, fecundity, survival rate, reproductive value, intrinsic rate of increase, finite rate of increase and net reproduction rate declined significantly when exposed to the combined treatment. No difference was observed between tannic acid and the control. Meanwhile, tannic acid with chlorantraniliprole had markedly antifeedant effects; causing significant decline in the relative growth rate (RGR), the relative consumption rate (RCR), the efficiency of conversion of ingested food (ECI), the efficiency of conversion of digested food and an increase in the approximate digestibility (AD) compared to either chemical alone. Tannic acid with chlorantraniliprole also decreased the insect's carbohydrate, lipid and protein contents significantly. The results showed that the interaction between tannic acid and chlorantraniliprole on the growth inhibition of larvae was additive and tannic acid increased the toxicity of chlorantraniliprole to insects. The results of this study provide information useful in integrated pest management programs for S. exigua and show that tannic acid combined with chlorantraniliprole may be a route to reducing the use of synthetic pesticides.

Comparative transcriptome analysis reveals host-associated differentiation in Chilo suppressalis (Lepidoptera: Crambidae)

The striped stem borer, Chilo suppressalis Walker (Lepidoptera: Crambidae), is one of the most serious rice pests. Besides attacking rice, it also feeds on an economically important vegetable crop, water-oat Zizania latifolia. The species feeding on water-oat has higher growth and survival rate than those on rice, suggesting their success in adaptation to the new host plant. However, little is known about the molecular mechanisms of host plant adaptation. Here we investigated the midgut transcriptome responses of C. suppressalis larvae reared on rice and water-oat. A total of 1,633 differentially expressed genes were identified, with a greater number up-regulated on the more delicious new host. The up-regulation of most digestive and detoxification-related genes may be the result of adaptation to the changes in nutritional requirements and toxic chemicals during host shift. In contrast, down-regulation of ribosomal genes may be related to their better development performance when feeding on the new host. In conclusion, our results suggest that transcriptional regulation of genes related to digestion, detoxification and ribosome may play an important role in adaptation of C. suppressalis to a new host plant.

Elevated carboxylesterase activity contributes to the lambda-cyhalothrin insensitivity in quercetin fed Helicoverpa armigera (Hübner)

Quercetin as one of the key plant secondary metabolite flavonol is ubiquitous in terrestrial plants. In this study, the decrease in sensitivity to lambda-cyhalothrin was observed in quercetin-fed Helicoverpa armigera larvae. In order to figure out the mechanisms underlying the decreased sensitivity of H. armigera larvae to lambda-cyhalothrin by quercetin induction, the changes in carboxylesterase activity and in-vitro hydrolytic metabolic capacity to lambda-cyhalothrin were examined. The LC50 value of quercetin-fed H. armigera larvae to lambda-cyhalothrin showed 2.41-fold higher than that of the control. S, S, S-Tributyl phosphorotrithioate (DEF) treatment showed a synergism effect on lambda-cyhalothrin toxicity to quercetin-fed H. armigera. Moreover, the activity of carboxylesterase was significantly higher in quercetin-fed H. armigera larvae after fed on quercetin for 48 h. The in-vitro hydrolytic metabolic capacity to lambda-cyhalothrin in quercetin-fed H. armigera larvae midgut was 289.82 nmol 3-PBA/mg protein/min, which is significant higher than that in the control group (149.60 nmol 3-PBA/mg protein/min). The elevated CarE enzyme activity and corresponding increased hydrolytic metabolic capacity to lambda-cyhalothrin in quercetin-fed H. armigera contributed to the enhanced tolerance to lambda-cyhalothrin.

Gene expression plasticity across hosts of an invasive scale insect species

For plant-eating insects, we still have only a nascent understanding of the genetic basis of host-use promiscuity. Here, to improve that situation, we investigated host-induced gene expression plasticity in the invasive lobate lac scale insect, Paratachardina pseudolobata (Hemiptera: Keriidae). We were particularly interested in the differential expression of detoxification and effector genes, which are thought to be critical for overcoming a plant's chemical defenses. We collected RNA samples from P. pseudolobata on three different host plant species, assembled transcriptomes de novo, and identified transcripts with significant host-induced gene expression changes. Gene expression plasticity was pervasive, but the expression of most detoxification and effector genes was insensitive to the host environment. Nevertheless, some types of detoxification genes were more differentially expressed than expected by chance. Moreover, we found evidence of a trade-off between expression of genes involved in primary and secondary metabolism; hosts that induced lower expression of genes for detoxification induced higher expression of genes for growth. Our findings are largely consonant with those of several recently published studies of other plant-eating insect species. Thus, across plant-eating insect species, there may be a common set of gene expression changes that enable host-use promiscuity.

Expression and RNA Interference of Ribosomal Protein L5 Gene in Nilaparvata lugens (Hemiptera: Delphacidae)

The ribosomal proteins play important roles in the growth and development of organisms. This study aimed to explore the function of NlRPL5 (GenBank KX379234), a ribosomal protein L5 gene, in the brown planthopper Nilaparvata lugens. The open reading frame of NlRPL5 was cloned from N. lugens based on a previous transcriptome analysis. The results revealed that the open reading frame of NlRPL5 is of 900 bp, encoding 299 amino acid residues. The reverse transcription quantitative PCR results suggested that the expression of NlRPL5 gene was stronger in gravid females, but was relatively low in nymphs, males, and newly emerged females. The expression level of NlRPL5 in the ovary was about twofolds of that in the head, thorax, or fat body. RNAi of dsNlRPL5 resulted in a significant reduction of mRNA levels, ∼50% decrease in comparison with the dsGFP control at day 6. Treatment of dsNlRPL5 significantly restricted the ovarian development, and decreased the number of eggs laid on the rice (Oryza sativa) plants. This study provided a new clue for further study on the function and regulation mechanism of NlRPL5 in N. lugens.

The Involvement of Heat Shock Proteins in the Establishment of Tomato Yellow Leaf Curl Virus Infection

Tomato yellow leaf curl virus (TYLCV), a begomovirus, induces protein aggregation in infected tomatoes and in its whitefly vector Bemisia tabaci. The interactions between TYLCV and HSP70 and HSP90 in plants and vectors are necessity for virus infection to proceed. In infected host cells, HSP70 and HSP90 are redistributed from a soluble to an aggregated state. These aggregates contain, together with viral DNA/proteins and virions, HSPs and components of the protein quality control system such as ubiquitin, 26S proteasome subunits, and the autophagy protein ATG8. TYLCV CP can form complexes with HSPs in tomato and whitefly. Nonetheless, HSP70 and HSP90 play different roles in the viral cell cycle in the plant host. In the infected host cell, HSP70, but not HSP90, participates in the translocation of CP from the cytoplasm into the nucleus. Viral amounts decrease when HSP70 is inhibited, but increase when HSP90 is downregulated. In the whitefly vector, HSP70 impairs the circulative transmission of TYLCV; its inhibition increases transmission. Hence, the efficiency of virus acquisition by whiteflies depends on the functionality of both plant chaperones and their cross-talk with other protein mechanisms controlling virus-induced aggregation.

Whitefly interactions with plants

Whiteflies are important pests of many crops worldwide. They are polyphagous and effectively feed on phloem sap using mouthparts modified into long, flexible stylets. Plants respond to whitefly attack by activating defense genes leading to production of toxic compounds. To reach plant phloem and survive on host plants, whiteflies secret effectors in the saliva to regulate plant responses and activate detoxification system to cope with plant defenses. Additionally, whitefly-transmitted viruses may exert substantial effects on host plants and in turn the performance of whiteflies. Understanding the interactions between whiteflies and host plants will promote the development of novel strategies for controlling whiteflies. Here, we summarize the genetics, molecular genetics and genomics of the whitefly's interactions with plants.

The draft genome of whitefly Bemisia tabaci MEAM1, a global crop pest, provides novel insights into virus transmission, host adaptation, and insecticide resistance

Background

The whitefly Bemisia tabaci (Hemiptera: Aleyrodidae) is among the 100 worst invasive species in the world. As one of the most important crop pests and virus vectors, B. tabaci causes substantial crop losses and poses a serious threat to global food security.

Results

We report the 615-Mb high-quality genome sequence of B. tabaci Middle East-Asia Minor 1 (MEAM1), the first genome sequence in the Aleyrodidae family, which contains 15,664 protein-coding genes. The B. tabaci genome is highly divergent from other sequenced hemipteran genomes, sharing no detectable synteny. A number of known detoxification gene families, including cytochrome P450s and UDP-glucuronosyltransferases, are significantly expanded in B. tabaci. Other expanded gene families, including cathepsins, large clusters of tandemly duplicated B. tabaci-specific genes, and phosphatidylethanolamine-binding proteins (PEBPs), were found to be associated with virus acquisition and transmission and/or insecticide resistance, likely contributing to the global invasiveness and efficient virus transmission capacity of B. tabaci. The presence of 142 horizontally transferred genes from bacteria or fungi in the B. tabaci genome, including genes encoding hopanoid/sterol synthesis and xenobiotic detoxification enzymes that are not present in other insects, offers novel insights into the unique biological adaptations of this insect such as polyphagy and insecticide resistance. Interestingly, two adjacent bacterial pantothenate biosynthesis genes, panB and panC, have been co-transferred into B. tabaci and fused into a single gene that has acquired introns during its evolution.

Conclusions

The B. tabaci genome contains numerous genetic novelties, including expansions in gene families associated with insecticide resistance, detoxification and virus transmission, as well as numerous horizontally transferred genes from bacteria and fungi. We believe these novelties likely have shaped B. tabaci as a highly invasive polyphagous crop pest and efficient vector of plant viruses. The genome serves as a reference for resolving the B. tabaci cryptic species complex, understanding fundamental biological novelties, and providing valuable genetic information to assist the development of novel strategies for controlling whiteflies and the viruses they transmit.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0321-y) contains supplementary material, which is available to authorized users.

Variation in life history traits and transcriptome associated with adaptation to diet shifts in the ladybird Cryptolaemus montrouzieri

BACKGROUND

Despite the broad diet range of many predatory ladybirds, the mechanisms involved in their adaptation to diet shifts are not completely understood. Here, we explored how a primarily coccidophagous ladybird Cryptolaemus montrouzieri adapts to feeding on aphids.

RESULTS

Based on the lower survival rate, longer developmental time, and lower adult body weight and reproduction rate of the predator, the aphid Megoura japonica proved being less suitable to support C. montrouzieri as compared with the citrus mealybug Planococcus citri. The results indicated up-regulation of genes related to ribosome and translation in fourth instars, which may be related to their suboptimal development. Also, several genes related to biochemical transport and metabolism, and detoxification were up-regulated as a result of adaptation to the changes in nutritional and non-nutritional (toxic) components of the prey.

CONCLUSION

Our results indicated that C. montrouzieri succeeded in feeding on aphids by regulation of genes related to development, digestion and detoxification. Thus, we argue that these candidate genes are valuable for further studies of the functional evolution of ladybirds led by diet shifts.

Only a minority of broad-range detoxification genes respond to a variety of phytotoxins in generalist Bemisia tabaci species

Generalist insect can utilize two different modes for regulating their detoxification genes, the constitutive mode and the induced mode. Here, we used the Bemisia tabaci sibling species MEAM1 and MED, as a model system for studying constitutive and induced detoxification resistance and their associated tradeoffs. B. tabaci adults were allowed to feed through membranes for 24 h on diet containing only sucrose or sucrose with various phytotoxins. Quantitative real-time PCR analyses of 18 detoxification genes, indicated that relatively few transcripts were changed in both the MEAM1 and MED species, in response to the addition of phytotoxins to the diet. Induced transcription of detoxification genes only in the MED species, in response to the presence of indole-3-carbinol in the insect’s diet, was correlated with maintenance of reproductive performance in comparison to significant reduction in performance of the MEAM1 species. Three genes, COE2, CYP6-like 5 and BtGST2, responded to more than one compound and were highly transcribed in the insect gut. Furthermore, functional assays showed that the BtGST2 gene encodes a protein capable of interacting with both flavonoids and glucosinolates. In conclusion, several detoxification genes were identified that could potentially be involved in the adaptation of B. tabaci to its host plants.

Differential expression of the chemosensory transcriptome in two populations of the stemborer Sesamia nonagrioides

Among the proposed mechanisms of local adaptation to different ecological environments, transcriptional changes may play an important role. In this study, we investigated whether such variability occurred within the chemosensory organs of a herbivorous insect, for which chemosensation guides most of its host preferences. A European and an African population of the noctuid Sesamia nonagrioides that display significant differences in their ecological preferences were collected on Zea mays and Typha domingensis, respectively. RNAseq were used between the two populations for digital expression profiling of chemosensory organs from larval antennae and palps. Preliminary data on adult female antennae and ovipositors were also collected. We found 6,550 differentially expressed transcripts in larval antennae and palps. Gene ontology enrichment analyses suggested that transcriptional activity was overrepresented in the French population and that virus and defense activities were overrepresented in the Kenyan population. In addition, we found differential expression of a variety of cytochrome P450s, which may be linked to the different host-plant diets. Looking at olfactory genes, we observed differential expression of numerous candidate odorant-binding proteins, chemosensory proteins, and one olfactory receptor, suggesting that differences in olfactory sensitivity participate in insect adaptation.

We can't all be supermodels: the value of comparative transcriptomics to the study of non-model insects

Insects are the most diverse group of organisms on the planet. Variation in gene expression lies at the heart of this biodiversity and recent advances in sequencing technology have spawned a revolution in researchers' ability to survey tissue-specific transcriptional complexity across a wide range of insect taxa. Increasingly, studies are using a comparative approach (across species, sexes and life stages) that examines the transcriptional basis of phenotypic diversity within an evolutionary context. In the present review, we summarize much of this research, focusing in particular on three critical aspects of insect biology: morphological development and plasticity; physiological response to the environment; and sexual dimorphism. A common feature that is emerging from these investigations concerns the dynamic nature of transcriptome evolution as indicated by rapid changes in the overall pattern of gene expression, the differential expression of numerous genes with unknown function, and the incorporation of novel, lineage-specific genes into the transcriptional profile.

Host plant-specific remodeling of midgut physiology in the generalist insect herbivore Trichoplusia ni

Species diversity in terrestrial ecosystems is influenced by plant defense compounds that alter the behavior, physiology, and host preference of insect herbivores. Although it is established that insects evolved the ability to detoxify specific allelochemicals, the mechanisms by which polyphagous insects cope with toxic compounds in diverse host plants are not well understood. Here, we used defended and non-defended plant genotypes to study how variation in chemical defense affects midgut responses of the lepidopteran herbivore Trichoplusia ni, which is a pest of a wide variety of native and cultivated plants. The genome-wide midgut transcriptional response of T. ni larvae to glucosinolate-based defenses in the crucifer Arabidopsis thaliana was characterized by strong induction of genes encoding Phase I and II detoxification enzymes. In contrast, the response of T. ni to proteinase inhibitors and other jasmonate-regulated defenses in tomato (Solanum lycopersicum) was dominated by changes in the expression of digestive enzymes and, strikingly, concomitant repression of transcripts encoding detoxification enzymes. Unbiased proteomic analyses of T. ni feces demonstrated that tomato defenses remodel the complement of T.ni digestive enzymes, which was associated with increased amounts of serine proteases and decreased lipase protein abundance upon encountering tomato defense chemistry. These collective results indicate that T. ni adjusts its gut physiology to the presence of host plant-specific chemical defenses, and further suggest that plants may exploit this digestive flexibility as a defensive strategy to suppress the production of enzymes that detoxify allelochemicals.

Activation of the Phenylpropanoid pathway in Nicotiana tabacum improves the performance of the whitefly Bemisia tabaci via reduced jasmonate signaling

BACKGROUND

Phloem-feeding insects can manipulate plant-induced resistance and are able to suppress effective jasmonic acid/ethylene (JA/ET) defenses by the induction of inefficient salicylic acid (SA) based responses. As a result, activation of the phenylpropanoid biosynthesis pathway in transgenic plants is anticipated to cause complex interactions between phloem-feeding insects and their host plants due to predicted contradiction between two defense forces: the toxicity of various phenylpropanoids and the accumulation of SA via a branch of the activated pathway.

METHODOLOGY/PRINCIPAL FINDINGS

Here, we investigated the effect of activating the phenylpropanoids pathway in Nicotiana tabacum, by over-expression of the PAP1 transcription factor, on the whitefly Bemisia tabaci, a phloem-feeding insect model. Our performance assays indicated that the over-expression made the transgenic plants a more suitable host for B. tabaci than wild-type (WT) plants, although these plants accumulated significantly higher levels of flavonoids. Transcription analyses of indicator genes in the SA (PR1a) and JA/ET (ERF1, COI1 and AOC) pathways followed by quantification of the SA and JA hormone levels, indicated that B. tabaci infestation periods longer than 8 hours, caused higher levels of activity of SA signaling in transgenic plants and higher levels of JA/ET signaling in WT plants.

CONCLUSIONS/SIGNIFICANCE

Taken together, these results emphasize the important role JA/ET-induced defenses play in protecting plants from successful infestation by B. tabaci and likely other phloem-feeding insects. It also indicates the necessity of phloem feeders to suppress these defenses for efficient utilization of plant hosts. Our data also indicate that the defensive chemistry produced by the phenylpropanoids pathway has only a minor effect on the insect fitness.

Genes involved in the evolution of herbivory by a leaf-mining, Drosophilid fly

Herbivorous insects are among the most successful radiations of life. However, we know little about the processes underpinning the evolution of herbivory. We examined the evolution of herbivory in the fly, Scaptomyza flava, whose larvae are leaf miners on species of Brassicaceae, including the widely studied reference plant, Arabidopsis thaliana (Arabidopsis). Scaptomyza flava is phylogenetically nested within the paraphyletic genus Drosophila, and the whole genome sequences available for 12 species of Drosophila facilitated phylogenetic analysis and assembly of a transcriptome for S. flava. A time-calibrated phylogeny indicated that leaf mining in Scaptomyza evolved between 6 and 16 million years ago. Feeding assays showed that biosynthesis of glucosinolates, the major class of antiherbivore chemical defense compounds in mustard leaves, was upregulated by S. flava larval feeding. The presence of glucosinolates in wild-type (WT) Arabidopsis plants reduced S. flava larval weight gain and increased egg-adult development time relative to flies reared in glucosinolate knockout (GKO) plants. An analysis of gene expression differences in 5-day-old larvae reared on WT versus GKO plants showed a total of 341 transcripts that were differentially regulated by glucosinolate uptake in larval S. flava. Of these, approximately a third corresponded to homologs of Drosophila melanogaster genes associated with starvation, dietary toxin-, heat-, oxidation-, and aging-related stress. The upregulated transcripts exhibited elevated rates of protein evolution compared with unregulated transcripts. The remaining differentially regulated transcripts also contained a higher proportion of novel genes than the unregulated transcripts. Thus, the transition to herbivory in Scaptomyza appears to be coupled with the evolution of novel genes and the co-option of conserved stress-related genes.