Gene expression profiling in the thiamethoxam resistant and susceptible B-biotype sweetpotato whitefly, Bemisia tabaci

Transcriptional analysis of Bemisia tabaci MEAM1 cryptic species under the selection pressure of neonicotinoids imidacloprid, acetamiprid and thiamethoxam

BACKGROUND

Neonicotinoids are widely applied in the control of the destructive agricultural pest Bemisia tabaci, and resistance against these chemicals has become a common, severe problem in the control of whiteflies. To investigate the molecular mechanism underlying resistance against nenonicotinoids in whiteflies, RNA-seq technology was applied, and the variation in the transcriptomic profiles of susceptible whiteflies and whiteflies selected by imidacloprid, acetamiprid and thiamethoxam treatment was characterized.

RESULTS

A total of 90.86 GB of clean sequence data were obtained from the 4 transcriptomes. Among the 16,069 assembled genes, 584, 110 and 147 genes were upregulated in the imidacloprid-selected strain (IMI), acetamiprid-selected strain (ACE), and thiamethoxam (THI)-selected strain, respectively, relative to the susceptible strain. Detoxification-related genes including P450s, cuticle protein genes, GSTs, UGTs and molecular chaperone HSP70s were overexpressed in the selected resistant strains, especially in the IMI strain. Five genes were downregulated in all three selected resistant strains, including 2 UDP-glucuronosyltransferase 2B18-like genes (LOC 109030370 and LOC 109032577).

CONCLUSIONS

Ten generations of selection with the three neonicotinoids induced different resistance levels and gene expression profiles, mainly involving cuticle protein and P450 genes, in the three selected resistant whitefly strains. The results provide a reference for research on resistance and cross-resistance against neonicotinoids in B. tabaci.

The neonicotinoid thiamethoxam impairs male fertility in solitary bees, Osmia cornuta

The ongoing loss of global biodiversity is endangering ecosystem functioning and human food security. While environmental pollutants are well known to reduce fertility, the potential effects of common neonicotinoid insecticides on insect fertility remain poorly understood. Here, we show that field-realistic neonicotinoid exposure can drastically impact male insect fertility. In the laboratory, male and female solitary bees Osmia cornuta were exposed to four concentrations of the neonicotinoid thiamethoxam to measure survival, food consumption, and sperm traits. Despite males being exposed to higher dosages of thiamethoxam, females revealed an overall increased hazard rate for survival; suggesting sex-specific differences in toxicological sensitivity. All tested sublethal concentrations (i.e., 1.5, 4.5 and 10 ng g^-1) reduced sperm quantity by 57% and viability by 42% on average, with the lowest tested concentration leading to a reduction in total living sperm by 90%. As the tested sublethal concentrations match estimates of global neonicotinoid pollution, this reveals a plausible mechanism for population declines, thereby reflecting a realistic concern. An immediate reduction in environmental pollutants is required to decelerate the ongoing loss of biodiversity.

Role of several cytochrome P450s in the resistance and cross-resistance against imidacloprid and acetamiprid of Bemisia tabaci (Hemiptera: Aleyrodidae) MEAM1 cryptic species in Xinjiang, China

Neonicotinoids are commonly used for the control of the whitefly Bemisia tabaci in cotton field. Laboratory test and field experiments have found that whitefly has a high risk of developing resistance and cross-resistance to the pesticide. Over expression of cytochrome P450 is one of the main mechanism that controls pesticide resistance in many insects. In this study we use MEAM1 whitefly, the dominant cryptic species of B. tabaci in Xinjiang cotton field, to investigate the possible resistance and cross-resistance mechanism controlled by cytochrome P450 enzymes. The P450 enzyme activity was higher in both selected strains of imidacloprid and acetamipird than that of susceptible strain. Synergism test showed that piperonyl butoxide (PBO) distinctly increased the control efficiency of imidacloprid and acetamiprid to the two resistance selected strains. Four out of 13 cytochrome genes, CYP4CS3, CYP6CX5, CYP6DW2 and CYP6CM1 were significantly up-regulated in the two selected strains based on real-time fluorescence quantitative PCR results. Other 3 genes, CYP6CX2, CYP6CX4 and CYP6DW3 were only highly expressed in the acetamiprid selected strain instead of the susceptible strain and imidacloprid selected strain. CYP6CM1 showed the highest expression level among all the 13 tested genes. No functional mutation of CYP6CM1 was found by sequence analysis. The possible role of these genes involving the resistance and cross-resistance of the whitefly MEAM1 cryptic species against neonicotinoids was discussed.

Two Fungicides Alter Reproduction of the Small Brown Planthopper Laodelphax striatellus by Influencing Gene and Protein Expression

Aside from their intended actions, fungicides can drive pest insect outbreaks due to virtually continuous use and pest evolution. Small brown planthopper (SBPH), Laodelphax striatellus, outbreaks occurred recently in many provinces in China, with devastating rice losses. Because exposure to the fungicide jinggangmycin (JGM) increased reproduction of the brown plant hopper, Nilaparvata lugens, via its influence on fatty acid synthase, we posed the hypothesis that JGM and carbendazim (CBM) influence SBPH reproduction via their influence on enzymes involved in other aspects of lipid metabolism. Exposure to the fungicide CBM stimulated SBPH reproduction (egg-laying up by 78%) and to another fungicide, JGM, led to decreased egg-laying (down by 47.3%). These inverse effects are mediated by down-regulated expression of l-3-hydroxyacyl-coenzyme A dehydrogenase (LCHAD) in JGM-treated females and up-regulated expression of hydroxysteroid dehydrogenase-like protein 2-like (HSD) in CBM-treated females. RNAi knockdown of, separately, LCHAD and HSD led to reduced egg-laying (down by 52% for dsLCHAD and by 73% for dsHSD). dsLCHAD, dsHSD, and JGM treatments also led to severely reduced ovarian development in experimental SBPH, with shorted and thinned valvula and lack of egg cells in ovaries. Valvula of CBM-treated females enlarged, with banana-shaped eggs in ovaries. These data strongly support our hypothesis.

The invasive MED/Q Bemisia tabaci genome: a tale of gene loss and gene gain

BACKGROUND

Sweetpotato whitefly, Bemisia tabaci MED/Q and MEAM1/B, are two economically important invasive species that cause considerable damages to agriculture crops through direct feeding and indirect vectoring of plant pathogens. Recently, a draft genome of B. tabaci MED/Q has been assembled. In this study, we focus on the genomic comparison between MED/Q and MEAM1/B, with a special interest in MED/Q's genomic signatures that may contribute to the highly invasive nature of this emerging insect pest.

RESULTS

The genomes of both species share similarity in syntenic blocks, but have significant divergence in the gene coding sequence. Expansion of cytochrome P450 monooxygenases and UDP glycosyltransferases in MED/Q and MEAM1/B genome is functionally validated for mediating insecticide resistance in MED/Q using in vivo RNAi. The amino acid biosynthesis pathways in MED/Q genome are partitioned among the host and endosymbiont genomes in a manner distinct from other hemipterans. Evidence of horizontal gene transfer to the host genome may explain their obligate relationship. Putative loss-of-function in the immune deficiency-signaling pathway due to the gene loss is a shared ancestral trait among hemipteran insects.

CONCLUSIONS

The expansion of detoxification genes families, such as P450s, may contribute to the development of insecticide resistance traits and a broad host range in MED/Q and MEAM1/B, and facilitate species' invasions into intensively managed cropping systems. Numerical and compositional changes in multiple gene families (gene loss and gene gain) in the MED/Q genome sets a foundation for future hypothesis testing that will advance our understanding of adaptation, viral transmission, symbiosis, and plant-insect-pathogen tritrophic interactions.

RNA interference-mediated knockdown of the hydroxyacid-oxoacid transhydrogenase gene decreases thiamethoxam resistance in adults of the whitefly Bemisia tabaci

Bemisia tabaci has developed a high level of resistance to thiamethoxam, a second generation neonicotinoid insecticide that has been widely used to control this pest. In this study, we investigated whether hydroxyacid-oxoacid transhydrogenase (HOT) is involved in resistance to the neonicotinoid insecticide thiamethoxam in the whitefly. We cloned the full-length gene that encodes HOT in B. tabaci. Its cDNA contains a 1428-bp open reading frame encoding 475 amino acid residues. Then we evaluated the mRNA expression level of HOT in different developmental stages, and found HOT expression was significantly greater in thiamethoxam resistance adults than in thiamethoxam susceptible adults. Subsequently, seven field populations of B. tabaci adults were sampled, the expression of mRNA level of HOT significant positive correlated with thiamethoxam resistance level. At last, we used a modified gene silencing system to knock-down HOT expression in B. tabaci adults. The results showed that the HOT mRNA levels decreased by 57% and thiamethoxam resistance decreased significantly after 2 days of feeding on a diet containing HOT dsRNA. The results indicated that down-regulation of HOT expression decreases thiamethoxam resistance in B. tabaci adults.

Transcription analysis of neonicotinoid resistance in Mediterranean (MED) populations of B. tabaci reveal novel cytochrome P450s, but no nAChR mutations associated with the phenotype

Background

Bemisia tabaci is one of the most damaging agricultural pests world-wide. Although its control is based on insecticides, B. tabaci has developed resistance against almost all classes of insecticides, including neonicotinoids.

Results

We employed an RNA-seq approach to generate genome wide expression data and identify genes associated with neonicotinoid resistance in Mediterranean (MED) B. tabaci (Q1 biotype). Twelve libraries from insecticide resistant and susceptible whitefly populations were sequenced on an Illumina Next-generation sequencing platform, and genomic sequence information of approximately 73 Gbp was generated.

A reference transcriptome was built by de novo assembly and functionally annotated. A total of 146 P450s, 18 GSTs and 23 CCEs enzymes (unigenes) potentially involved in the detoxification of xenobiotics were identified, along with 78 contigs encoding putative target proteins of six different insecticide classes. Ten unigenes encoding nicotinic Acetylcholine Receptors (nAChR), the target of neoinicotinoids, were identified and phylogenetically classified. No nAChR polymorphism potentially related with the resistant phenotypes, was observed among the studied strains.

DE analysis revealed that among the 550 differentially (logFC > 1) over-transcribed unigenes, 52 detoxification enzymes were over expressed including unigenes with orthologues in P450s, GSTs, CCE and UDP-glucuronosyltransferases.

Eight P450 unigenes belonging to clades CYP2, CYP3 and CYP4 were highly up-regulated (logFC > 2) including CYP6CM1, a gene already known to confer imidacloprid resistance in B. tabaci. Using quantitative qPCRs, a larger screening of field MED B. tabaci from Crete with known neonicotinoid phenotype was performed to associate expression levels of P450s with resistance levels. Expression levels of five P450s, including CYP6CM1, were found associated with neonicotinoid resistance. However, a significant correlation was found only in CYP303 and CYP6CX3, with imidacloprid and acetamiprid respectively.

Conclusion

Our work has generated new toxicological data and genomic resources which will significantly enrich the available dataset and substantially facilitate the molecular studies in MED B. tabaci. No evidence of target site neonicotinoid resistance has been found. Eight P450 unigenes, including CYP6CM1, were found significantly over-expressed in resistant B. tabaci. This study suggests at least two novel P450s (CYP303 and CYP6CX3) as candidates for their functional characterization as detoxification mechanisms of neonicotinoid resistance in B. tabaci.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2161-5) contains supplementary material, which is available to authorized users.

The global status of insect resistance to neonicotinoid insecticides

The first neonicotinoid insecticide, imidacloprid, was launched in 1991. Today this class of insecticides comprises at least seven major compounds with a market share of more than 25% of total global insecticide sales. Neonicotinoid insecticides are highly selective agonists of insect nicotinic acetylcholine receptors and provide farmers with invaluable, highly effective tools against some of the world's most destructive crop pests. These include sucking pests such as aphids, whiteflies, and planthoppers, and also some coleopteran, dipteran and lepidopteran species. Although many insect species are still successfully controlled by neonicotinoids, their popularity has imposed a mounting selection pressure for resistance, and in several species resistance has now reached levels that compromise the efficacy of these insecticides. Research to understand the molecular basis of neonicotinoid resistance has revealed both target-site and metabolic mechanisms conferring resistance. For target-site resistance, field-evolved mutations have only been characterized in two aphid species. Metabolic resistance appears much more common, with the enhanced expression of one or more cytochrome P450s frequently reported in resistant strains. Despite the current scale of resistance, neonicotinoids remain a major component of many pest control programmes, and resistance management strategies, based on mode of action rotation, are of crucial importance in preventing resistance becoming more widespread. In this review we summarize the current status of neonicotinoid resistance, the biochemical and molecular mechanisms involved, and the implications for resistance management.

Identification of salt gland-associated genes and characterization of a dehydrin from the salt secretor mangrove Avicennia officinalis

BACKGROUND

Salt stress is a major challenge for growth and development of plants. The mangrove tree Avicennia officinalis has evolved salt tolerance mechanisms such as salt secretion through specialized glands on its leaves. Although a number of structural studies on salt glands have been done, the molecular mechanism of salt secretion is not clearly understood. Also, studies to identify salt gland-specific genes in mangroves have been scarce.

RESULTS

By subtractive hybridization (SH) of cDNA from salt gland-rich cell layers (tester) with mesophyll tissues as the driver, several Expressed Sequence Tags (ESTs) were identified. The major classes of ESTs identified include those known to be involved in regulating metabolic processes (37%), stress response (17%), transcription (17%), signal transduction (17%) and transport functions (12%). A visual interactive map generated based on predicted functional gene interactions of the identified ESTs suggested altered activities of hydrolase, transmembrane transport and kinases. Quantitative Real-Time PCR (qRT-PCR) was carried out to validate the expression specificity of the ESTs identified by SH. A Dehydrin gene was chosen for further experimental analysis, because it is significantly highly expressed in salt gland cells, and dehydrins are known to be involved in stress remediation in other plants. Full-length Avicennia officinalis Dehydrin1 (AoDHN1) cDNA was obtained by Rapid Amplification of cDNA Ends. Phylogenetic analysis and further characterization of this gene suggested that AoDHN1 belongs to group II Late Embryogenesis Abundant proteins. qRT-PCR analysis of Avicennia showed up-regulation of AoDHN1 in response to salt and drought treatments. Furthermore, some functional insights were obtained by growing E. coli cells expressing AoDHN1. Growth of E. coli cells expressing AoDHN1 was significantly higher than that of the control cells without AoDHN1 under salinity and drought stresses, suggesting that the mangrove dehydrin protein helps to mitigate the abiotic stresses.

CONCLUSIONS

Thirty-four ESTs were identified to be enriched in salt gland-rich tissues of A. officinalis leaves. qRT-PCR analysis showed that 10 of these were specifically enriched in the salt gland-rich tissues. Our data suggest that one of the selected genes, namely, AoDHN1 plays an important role to mitigate salt and drought stress responses.

Identification and validation of reference genes for quantitative real-time PCR in Drosophila suzukii (Diptera: Drosophilidae)

To accurately evaluate gene expression levels and obtain more accurate quantitative real-time RT-PCR (qRT-PCR) data, normalization relative to reliable reference gene(s) is required. Drosophila suzukii, is an invasive fruit pest native to East Asia, and recently invaded Europe and North America, the stability of its reference genes have not been previously investigated. In this study, ten candidate reference genes (RPL18, RPS3, AK, EF-1β, TBP, NADH, HSP22, GAPDH, Actin, α-Tubulin), were evaluated for their suitability as normalization genes under different biotic (developmental stage, tissue and population), and abiotic (photoperiod, temperature) conditions. The three statistical approaches (geNorm, NormFinder and BestKeeper) and one web-based comprehensive tool (RefFinder) were used to normalize analysis of the ten candidate reference genes identified α-Tubulin, TBP and AK as the most stable candidates, while HSP22 and Actin showed the lowest expression stability. We used three most stable genes (α-Tubulin, TBP and AK) and one unstably expressed gene to analyze the expression of P-glycoprotein in abamectin-resistant and sensitive strains, and the results were similar to reference genes α-Tubulin, TBP and AK, which show good stability, while the result of HSP22 has a certain bias. The three validated reference genes can be widely used for quantification of target gene expression with qRT-PCR technology in D.suzukii.

Di-(2-ethylhexyl) phthalate and bisphenol A exposure impairs mouse primordial follicle assembly in vitro

Bisphenol-A (BPA) and diethylhexyl phthalate (DEHP) are estrogenic compounds widely used in commercial plastic products. Previous studies have shown that exposure to such compounds have adverse effects on various aspects of mammalian reproduction including folliculogenesis. The objective of this study was to examine the effects of BPA and DEHP exposure on primordial follicle formation. We found that germ cell nest breakdown and primordial follicle assembly were significantly reduced when newborn mouse ovaries were exposed to 10 or 100 μM BPA and DEHP in vitro. Moreover, BPA and DEHP exposure increased the number of TUNEL positive oocytes and the mRNA level of the pro-apoptotic gene Bax in oocytes. These effects were associated with decreased expression of oocyte specific genes such as LIM homeobox 8 (Lhx8), factor in the germline alpha (Figla), spermatogenesis and oogenesis helix-loop-helix (Sohlh2), and newborn ovary homeobox (Nobox). Interestingly, BPA and DEHP exposure also prevented DNA demethylation of CpG sites of the Lhx8 gene in oocytes, a process normally associated with folliculogenesis. Finally, folliculogenesis was severely impaired in BPA and DEHP exposed ovaries after transplantation into the kidney capsules of immunodeficient mice. In conclusion, BPA and DEHP exposures impair mouse primordial follicle assembly in vitro.

Transcriptome profiling of the whitefly Bemisia tabaci reveals stage-specific gene expression signatures for thiamethoxam resistance

Bemisia tabaci has developed high levels of resistance to many insecticides including the neonicotinoids and there is strong evidence that for some compounds resistance is stage-specific. To investigate the molecular basis of B. tabaci resistance to the neonicotinoid thiamethoxam we used a custom whitefly microarray to compare gene expression in the egg, nymph and adult stages of a thiamethoxam-resistant strain (TH-R) with a susceptible strain (TH-S). Gene ontology and bioinformatic analyses revealed that in all life stages many of the differentially expressed transcripts encoded enzymes involved in metabolic processes and/or metabolism of xenobiotics. Several of these are candidate resistance genes and include the cytochrome P450 CYP6CM1, which has been shown to confer resistance to several neonicotinoids previously, a P450 belonging to the Cytochrome P450s 4 family and a glutathione S-transferase (GST) belonging to the sigma class. Finally several ATP-binding cassette transporters of the ABCG subfamily were highly over-expressed in the adult stage of the TH-R strain and may play a role in resistance by active efflux. Here, we evaluated both common and stage-specific gene expression signatures and identified several candidate resistance genes that may underlie B. tabaci resistance to thiamethoxam.

Transcriptomic and proteomic responses of sweetpotato whitefly, Bemisia tabaci, to thiamethoxam

BACKGROUND

The sweetpotato whitefly, Bemisia tabaci (Hemiptera: Aleyrodidae), is one of the most widely distributed agricultural pests. Although it has developed resistance to many registered insecticides including the neonicotinoid insecticide thiamethoxam, the mechanisms that regulate the resistance are poorly understood. To understand the molecular basis of thiamethoxam resistance, "omics" analyses were carried out to examine differences between resistant and susceptible B. tabaci at both transcriptional and translational levels.

RESULTS

A total of 1,338 mRNAs and 52 proteins were differentially expressed between resistant and susceptible B. tabaci. Among them, 11 transcripts had concurrent transcription and translation profiles. KEGG analysis mapped 318 and 35 differentially expressed genes and proteins, respectively, to 160 and 59 pathways (p<0.05). Thiamethoxam treatment activated metabolic pathways (e.g., drug metabolism), in which 118 transcripts were putatively linked to insecticide resistance, including up-regulated glutathione-S-transferase, UDP glucuronosyltransferase, glucosyl/glucuronosyl transferase, and cytochrome P450. Gene Ontology analysis placed these genes and proteins into protein complex, metabolic process, cellular process, signaling, and response to stimulus categories. Quantitative real-time PCR analysis validated "omics" response, and suggested a highly overexpressed P450, CYP6CX1, as a candidate molecular basis for the mechanistic study of thiamethoxam resistance in whiteflies. Finally, enzymatic activity assays showed elevated detoxification activities in the resistant B. tabaci.

CONCLUSIONS

This study demonstrates the applicability of high-throughput omics tools for identifying molecular candidates related to thiamethoxam resistance in an agricultural important insect pest. In addition, transcriptomic and proteomic analyses provide a solid foundation for future functional investigations into the complex molecular mechanisms governing the neonicotinoid resistance in whiteflies.