The survival of a species depends on its capacity to adjust to changing environmental conditions, and new stressors. Such new, anthropogenic stressors include the neonicotinoid class of crop-protecting agents, which have been implicated in the population declines of pollinating insects, including honeybees (Apis mellifera). The low-dose effects of these compounds on larval development and physiological responses have remained largely unknown. Over a period of 15 days, we provided syrup tainted with low levels (2 µg/L −1 ) of the neonicotinoid insecticide imidacloprid to beehives located in the field. We measured transcript levels by RNA sequencing and established lipid profiles using liquid chromatography coupled with mass spectrometry from worker-bee larvae of imidacloprid-exposed (IE) and unexposed, control (C) hives. Within a catalogue of 300 differentially expressed transcripts in larvae from IE hives, we detect significant enrichment of genes functioning in lipid-carbohydrate-mitochondrial metabolic networks. Myc-involved transcriptional response to exposure of this neonicotinoid is indicated by overrepresentation of E-box elements in the promoter regions of genes with altered expression. RNA levels for a cluster of genes encoding detoxifying P450 enzymes are elevated, with coordinated downregulation of genes in glycolytic and sugar-metabolising pathways. Expression of the environmentally responsive Hsp90 gene is also reduced, suggesting diminished buffering and stability of the developmental program. The multifaceted, physiological response described here may be of importance to our general understanding of pollinator health. Muscles, for instance, work at high glycolytic rates and flight performance could be impacted should low levels of this evolutionarily novel stressor likewise induce downregulation of energy metabolising genes in adult pollinators.

Competing interests: This project was supported by The Co-operative Group, UK. Alessandro Guffanti, Paolo Pavan and Anna Moles are employed by Genomnia srl, (Milan/Italy). Thomas Ryder is owner of Parks Apiaries (Nottingham/UK). There are no patents, products in development or marketed products to declare. This does not alter the authors' adherence to all the PLOS ONE policies on sharing data and materials.

Copyright: © 2013 Derecka et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Here we present and analyse molecular profiles obtained from worker-honeybee larvae after hives in the field were given access to syrup tainted with a low dose of the widely used neonicotinoid insecticide, imidacloprid. Genome-wide RNA transcriptional responses and lipid profiles provide insight into the physiological responses of developing organisms confronted with a new, real-world environmental threat.

Recent field studies showed that even low doses of neonicotinoid-contaminated food disturb the population dynamics of bumblebee (Bombus terrestris) colonies [17] , [18] , and weaken flight and homing performance of honeybees [19] , [20] . The brood of honeybees fully depends on food collected by foraging workers. It is conceivable, then, that neonicotinoids could be imported into the hive via contaminated nectar and pollen from surrounding insecticide-treated crops. The consequences of larval exposure to low levels of this novel stressor have not been studied under field-realistic conditions.

Because neonicotinoid compounds can be present in all parts of a plant, it is plausible that nectar-collecting and pollen-harvesting insects may ingest low doses of the insecticide, potentially affecting their health and behaviour [15] , [16] . Indeed, residual levels of neonicotinoids, as well as other synthetic compounds, have been detected in the nectar and pollen of pesticide-treated, flowering crops (reviewed in [15] ).

Entirely new environmental hazards, and those not encountered previously by a species present a conundrum. The organism must rapidly deal with the novel stressor in order to survive, but the genome-encoded options for responding are often limited. The synthetic class of neonicotinoid insecticides represents such a new type of environmental stressor. Registered for use in more than 120 countries and accounting for around 24% of the global market, neonicotinoids have been the most important class of agricultural pesticides in recent years [13] . The success of these synthetic insecticides in agricultural applications largely derives from two factors. First, neonicotinoids selectively target nicotinic acetylcholine receptors (nAChRs) of insects and as a result have replaced older, more hazardous organophosphate and organochlorine pesticides, which have a much higher toxicity for vertebrates [14] . Second, neonicotinoids act systemically; they are absorbed, distributed and maintained throughout the plant upon seed or soil treatment. The necessity of repeated crop spraying is reduced. Although development and application of neonicotinoids can be viewed favourably insofar as they contribute to an overall improvement to modern commercial agricultural practice and global food security, the systemic nature of their action has raised concern, as this very property may inadvertently harm insects considered beneficial to humans.

Taxing conditions often require energy-consuming physiological adjustments: altered metabolic demands must be coordinated with organismal growth to ensure completion of the developmental program under adverse conditions. Hormones participate in regulating metabolic aspects of the stress response [7] . In the fruit fly Drosophila, hormonal and nutritional signals are integrated by the transcription factor Myc (dMyc), which functions as an intermediate regulator of growth during development [8] . Insulin indirectly regulates dMyc in the larval fat body and muscle through FOXO and TOR complex 1 (TORC1) signaling pathways, causing changes in expression of a large metabolic gene network [9] , [10] . Both TORC1 and FOXO, which act upstream of dMyc, also respond to additional environmental cues, including immune challenges and changes in oxygen levels [11] , [12] . Thus, there is evidence for the evolution of developmental buffering systems to meet the challenges of environmental stressors recurrently experienced during life histories.

Species populating today’s Earth all derive from ancestors that prevailed in changeable environments. Waddington described the tendency of an organism to produce definitive phenotypes under variable conditions as ‘canalisation’, which implies some robustness of the genetically encoded developmental program [1] . Such canalisation is not absolute. With too great of an environmental stimulus, developmental buffering systems can fail, causing abnormities or death during development. A number of molecular mechanisms have been proposed and identified as contributing to canalisation. These mechanisms include microRNAs [2] , [3] , [4] , [5] , small non-coding RNAs that modulate networks of protein-coding genes and the heat-shock protein 90 (Hsp90), a molecular chaperone that stabilises developmental programs [6] . A characteristic property for both Hsp90 and microRNA buffering systems is their ability to facilitate rapid adjustments.

Results and Discussion

Field-realistic Feeding Trials of Control and Imidacloprid-exposed Hives To mimic honey flow of a neonicotinoid-treated nectar source in bloom, we provided colonies of free-foraging honeybees in the field an additional, imidacloprid-tainted source of food. Over a 15-day period, three experimental, imidacloprid-exposed (IE) hives received a daily ration of 100 mls syrup containing imidacloprid, while three control (C) hives received 100 mls of untainted syrup. We provided a concentration of imidacloprid (2 µg imidacloprid/L−1 syrup/≈2 parts per billion) that lies within the range (0.5 ppb–10 pbb) detected in contaminated pollen and nectar of a variety of crops (reviewed in [15]). During these 15-day feeding trials, all individuals within IE colonies, including egg-laying queens and the developing brood, will have come into contact with this insecticide through communal food exchange, a behavioural trait of social insects termed trophollaxis [21]. Following the 15-day feeding trials, we collected worker bee larvae for RNA transcriptome analysis and lipid profiling (Figure S1).

Altered Expression of 15 microRNAs MicroRNAs have emerged as modulators of cell physiological processes. We asked whether altered microRNA levels are detectable in imidacloprid-exposed larvae. From each of the three C and the three IE hives, respectively. RNA samples were isolated and pooled from whole larvae to establish six independent libraries for deep sequencing. Using statistical methodology implemented by edgeR and DESeq analysis [22] on the tag numbers of identified microRNAs, we established a list of 15 candidate microRNAs for which expression levels differed significantly (≥0.5 log 2 -fold change/p≤0.05; [23]) between the IE groups and C groups (Table 1, Data S1). Six of these microRNAs have been detected and annotated only in honeybees; their gene targets are presently unknown. The other nine microRNAs on the list are conserved in the evolutionarily distant fruit fly Drosophila. One of these genes is mir-14, which has diminished levels in the IE data set (−0.52 log 2 -fold change) (Table 1). In Drosophila, mir-14 deficiency is associated with a variety of phenotypes, including a lower probability of larvae to survive development to adulthood, reduced adult lifespan, stress-sensitivity and abnormal energy metabolism [24], [25]. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Table 1. DESeq-based expression levels identified, 15 annotated, mature candidate microRNAs (miRNAs), which show subtle differences in abundance between larval samples collected from imidacloprid-exposed (IE) hives and unexposed, control (C) hives. https://doi.org/10.1371/journal.pone.0068191.t001

Identification of Differentially Expressed Genes We next performed whole transcriptome sequencing (RNA-Seq) to explore whether steady-state levels of certain protein-encoding RNAs differ between worker larvae from C and IE hives. Following differential expression analysis using DEGseq (Table S1), we collated a list of 300 genes according to the DEGseq statistical tests FET, LRT, and MARS (p-value ≤0.001) and FC ≥0.5 (log 2 normalised fold change) representing differences in RNA expression between samples C and IE (Table S2). Of this list, 65% (195/300) of IE genes have reduced RNA levels, while 35% of the IE genes (105/300) have elevated RNA levels relative to the same genes in the C group (Figure 1A, Table S2). PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Figure 1. Functional annotation of differentially expressed genes. (A) Based on RNA-Seq data, expression levels for 300 genes were found to be significantly changed in imidacloprid-exposed larvae (compared with data obtained from non-exposed larvae). Expression is reduced for 195 genes (blue); expression is increased for 105 genes (red). The group of over-expressed genes includes nine cytochrome P450s; their gene IDs and normalised fold-change values are shown in the table. (B) Selection of non-redundant Gene Ontology (GO) terms for the three ontologies: Biological Process (BP), Molecular Function (MF) and Cellular Component (CC), and Interpro (IP) protein domains that are overrepresented in the 105 up (red) and 195 down (blue) regulated genes. Significance of enrichment is based on Fisher’s Exact Test. GO terms for the most significant enrichment groups are indicated in bold letters. https://doi.org/10.1371/journal.pone.0068191.g001

Increased Expression of Lipid and Xenobiotic Metabolising Cytochrome P450 Genes Genes are conveniently categorised by their known functions and predicted biological roles. Such gene ontology (GO) analyses revealed that our list of differentially expressed RNAs is significantly enriched for genes operating in a lipid-carbohydrate-mitochondrial metabolic network (p<10−5/Figure 1B). Within the overexpressed group of 105 transcripts, we find enrichment for nine genes of the cytochrome P450 monooxygenase superfamily (Figure 1A; Table S2). P450 (CYP) enzymes are oxidation catalysts of many cellular compounds, including lipids, steroid hormones and arachidonic acid metabolites. These enzymes also metabolise xenobiotic compounds and catalyse the breakdown of a wide range of structurally different toxins and synthetic insecticides. Overexpression of genes coding for the P450 clades (CYP4, CYP6 and CYP9) contribute considerably to insecticide-resistance [26]. In particular, filed and laboratory studies have shown a causal link between Cyp6g1 overexpression and resistance to DTT and neonicotinoids [27]. Consistent with those results, RNAi mediated knockdown of Cyp6g1 renders adult Drosophila more susceptible to imidacloprid [28]. Many insect CYP/P450 genes appear to be under direct or indirect control of the nuclear receptor dHR96 in response to structurally diverse xenobiotics [28], [29], [30]. Knockdown of dHR96 increases tolerance of adult Drosophila to imidacloprid exposure [28]. As dHR96 influences both gene activation and repression, imidacloprid tolerance could be mediated, in part, by upregulation of certain CYP/P450 genes in these knockdown flies [28]. The nine upregulated honeybee CYP/P450 genes in IE larval samples belong to the CYP4, CYP6 and CYP9 clades. It remains to be shown if altered transcription of these CYP/P450 genes is a specific detoxification-response and whether the encoded enzymes are capable of metabolising imidacloprid. Further experimental studies will be required to establish whether honeybees have a functional orthologue of the Drosophila Cyp6g1, the CYP/P450 enzyme capable of protecting insects against the neonicotinoid insecticide. Compared with Drosophila, the honeybee genome carries s only around half the number of CYP/P450 genes [31]. The shortfall of honeybee CYP/P450 genes has been suggested to decrease the efficiency of detoxification processes, explaining the vulnerability of this pollinating insect to pesticides [32].

Genes Involved in Lipid Metabolism Our finding of changes in mRNA levels for fatty acid metabolising enzymes suggests that lipid homeostasis is altered in imidacloprid-exposed larvae (Figure 1B). For example, reduced transcriptional activities are seen for fatty acid synthase (FAS), a glycerol-3-phosphate acyltransferase (GPAT) and the ATP citrate lyase, ATPCL (Table S2; Data S3). These three genes are downregulated in adult Drosophila as part of an immune response to infection; their lowered expression levels reflects energy reallocation, a process proposed to be controlled by the FOXO signalling pathway [33]. FOXO might be involved in regulating the aforementioned nuclear receptor dHR96 in Drosophila [34]. Besides responding to xenobiotic stresses, dHR96 also plays a key role in lipid and cholesterol metabolism, in adult Drosophila [30], [35], [36]. We examined if some of the 300 differentially expressed genes in IE honeybee larvae could be orthologues of lipid-metabolising genes found to be dHR96-regulated in Drosophila. Potential dHR96-regulated genes are GB12567, GB17220 and GB10584, encoding a predicted long-chain fatty acid-CoA ligase, a Lip3-like lipase and an enzyme involved in sphingolipid metabolic processes, respectively. That is, we only observe a limited overlap between our honeybee data set of genes taking part in lipid-metabolism and that reported for dHR96-regulated lipid genes in Drosophila. However, there was an overrepresentation of Drosophila orthologues that are upregulated (P<.004) and downregulated (P<10–9) by CncC, the cap ‘n’ collar isoform-C regulator of xenobiotic detoxification responses in Drosophila [37]. Another transcription factor that may be involved in the response to imidacloprid-exposure is the Hepatic Nuclear Factor 4 (HNF4) as we identified 15 orthologues of transcripts that show altered expression in Drosophila dHNF4 mutants [38]. Amongst them are the downregulated AcetylCoenzyme A synthase (AcCoAS) and an acetyl-CoA carboxylase (ACC) and two upregulated Lip3-like lipases (GB17220, GB17745) (Table S2). Elevated levels of Lip3, encoding a predicted cholesterol ester hydrolase, have been observed in starved Drosophila larva [39]. The dHNF4 transcription factor is activated in response to starvation or exogenous long-chain fatty acids and it has been suggested that dHNF4 senses free fatty acid levels in the fly larva and drives lipid mobilization and catabolism through lipolysis and β-oxidation in the mitochondria [38]. Together, these findings indicate that imidacloprid-exposure triggers changes in gene expression affecting various facets of larval lipid metabolism.

Concerted Effect on Glycolytic and Carbohydrate-metabolising genes Balancing energy requirements under varying conditions necessitates coordination of lipid and glucose-metabolising pathways. We found that genes of the gluconeogenesis and glycolysis pathways have altered expression upon imidacloprid exposure (Figure 2). Increased expression is detected for the stress-responsive phosphoenolpyruvate carboxykinase (PEPCK), a gene that regulates glucose levels by catalysing the first step in gluconeogenesis. In Drosophila, upregulation of PEPCK is induced by xenobiotic insults [29] and under fasting conditions [39], the latter involving the FOXO signaling pathway [40]. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Figure 2. Expression of genes encoding carbohydrate-metabolising enzymes is affected in imidacloprid-exposed worker bee larvae. Genes/enzymes, including paralogues, and their positions (coloured/grey boxes) in the glycolytic and related carbohydrate pathways are placed with reference to honeybee-specific pathway variations [72], [73]. Based on DEGseq analysis, ten genes are downregulated (blue), including PGI, PGK, PGLYM, ENO and PYK, of the glycolytic pathway; PEPCK is upregulated (red). A key for gene names is provided; fold-changes in expression as determined by RNA-Seq (IE relative to C data set). https://doi.org/10.1371/journal.pone.0068191.g002 Within the group of 195 downregulated mRNAs in the IE data set, we identified a cluster of sugar metabolising genes; five of them encode core enzymes of the glycolytic pathway (Figure 2). To better understand the potential interactions between the up- and downregulated genes, we used the gene lists as input in the online database resource Search Tool for the Retrieval of Interacting Genes (STRING) program. Analysis of the downregulated genes revealed a highly connected network with more than twice the number of putative interactions expected by chance (377 observed, 153 expected, P = 0). Our analysis showed strong links between the down regulation of genes involved in sugar metabolism and RNA translation in both the cell and the mitochondria (Figure S2). This STRING analysis further supports our discovery of a concerted, non-random down regulation of metabolic processes in IE larvae. That is, the reported transcriptional changes for the genes in our list cannot be explained by normal variation in expression found among bee colonies. While transcriptional changes of individual glycolytic genes have little effect on glucose flux [41], several reports indicate that coordinated changes in their expression levels occur as part of the physiological response to starvation. Diminished expression of genes involved in glycolysis has been observed in nutrient deprived Drosophila larvae [38] and can be modulated by at least two different transcription factors: the ecdysone receptor (EcR) and the estrogen-related receptor (dERR), respectively [42], [43]. The Drosophila ERR is an essential regulator of carbohydrate metabolism and a deficiency in mutant Drosophila larvae is associated with diminished ATP and triacylglyceride (TAG) levels [43]. The overall altered transcription profile of the IE data set (Figure 1), echoes the observations reported for xenobiotic stress responses in Drosophila [29]. Treatment of fruit flies with the drug phenobarbital is also associated with deregulation of genes involved in energy and sugar metabolism [29].

Possible Participation of Myc-mediated Gene Expression Coordinated responses of metabolic gene networks to a variety of different environmental stimuli suggest genome-wide, synchronised regulation by one or more transcription factors. We applied ‘Clover’, a program developed to identify statistical overrepresentation of functional DNA sequence motives [44] to the promoters of our 300 differentially expressed protein-coding genes (Data S2). A significant proportion of genes (19% [56/300]) contain the canonical enhancer box (E-box) sequence CACGTG, a DNA element recognised by the Myc family of transcription factors [45]. For the majority of genes containing this sequence motive (73% [41/56]), expression levels are diminished in the IE group (Data S2). We note that six of the ‘E-box genes’ have Drosophila orthologues that carry the same promoter DNA element; their expression has been shown to be TORC1-regulated and dMyc-dependent (Data S2) [9], [46]). One of the factors influencing the activity of dMyc in the fat body is ecdysone [47]. Elevated levels of this steroid hormone and its receptor repress dMyc; this leads to altered expression of Myc-target genes and a slowed larval growth rate [47]. Among the many potential genomic targets, Myc has been implicated in the control of nuclear-encoded genes affecting mitochondrial function and activity [48]. We find six of these genes to be downregulated in the IE-Myc data set; they encode mitochondrial membrane proteins Tim8, Tim9a, Tim13 and bor/dATAD3A, as well as the mitochondrial ribosomal proteins mRpL9 and mRpL12 (Data S2). Cells of a fruit fly lacking the protein mRpL12 have an impaired growth rate and reduced mitochondrial activity [49]. The cluster of downregulated mitochondrial genes in our ‘E-box list’ suggests that the function of this energy-producing organelle could be affected in imidacloprid-exposed larvae.

Decreased Expression Levels of Hsp90 The developmental Hsp90 buffering system is versatile. It stabilises misfolded and metabstable protein complexes [50], suppresses the activity of transposable elements [51] and is involved in establishing altered chromatin states [52]. Levels of the ubiquitously expressed Hsp90 gene are normally high. This ensures immediate buffering capacity if demand exceeds the required basal level [6], [53]. In Drosophila, Hsp90 is one of the TORC1/dMyc-regulated genes that contain a canonical E-box element in the promoter region [9]. The honeybee Hsp90 promoter carries the variant E-Box sequence CACATG, which may also be responsive to Myc. Our RNA-Seq expression data indicate reduced levels of Hsp90 transcripts in IE larvae (Table S2). Overall diminished Hsp90 RNA levels in IE larvae was also detected by multiplexed RT-PCR using additional larval samples (Figure 3). Reduction of Hsp90 levels would be a strategically sound response to a new environmental stressor. It allows the release from – and expression of – pre-existing cryptic genetic and epigenetic variants. Some of these variants may improve the chances of normal larval development, while others will increase imidacloprid sensitivity, potentially leading to larval death. If Hsp90 downregulation upon imidacloprid exposure is a common response, the genetic background may prove to be an important factor in determining the fate of a bee colony. Inbred populations of fruit flies, for example, with slightly decreased Hsp90 expression levels are less fecund and shorter-lived than populations with normal Hsp90 levels [54]. Together, our finding of reduced Hsp90 expression implies that imidacloprid-exposure affects the developmental buffering system. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Figure 3. Multiplex RT-PCR validation of HSP90 and a P450 gene (CYP9-clade). Variance in measured RNA expression levels was assessed between treatment groups (left), among hives (middle) and between the two measurements taken for individual larvae (right/larval samples 1–34 came from control hives C1–C3; Larval samples 35–74 came from imidacloprid-exposed hives IE1–IE3. As two measurements were made for each larva, the height of each bar in the larva plot gives the range of values recorded; the mean value for each larva is indicated in the middle of this range. For those larvae that yielded only a single successful measurement for a given gene, only one point is shown. Box and whisker plots show 25th percentile (bottom of box), 50th percentile (middle of box), 75th percentile (top of box), median (line in middle of box), and full range of data (low bar to high bar). Relative expression levels are compared to a standardisation factor. The scale of the y-axis differs between plots for the two genes and indicates fold-difference for the expression level relative to the standardisation factor. https://doi.org/10.1371/journal.pone.0068191.g003

Measuring RNA Levels in Samples of Individual Larvae We sought to validate results that had been obtained by RNA-Seq. From the list of 300 differentially expressed genes, a set of 17 genes was selected for analysis, including Hsp90 and members of the cytochrome P450 clades CYP6 and CYP9 (Figure 3; Data S3). We employed a multiplexed, quantitative gene expression analysis system (GenomeLab GeXP/Beckman Coulter [55]. This approach allowed us to simultaneously measure relative mRNA levels of the 17 genes in a given RNA sample and perform ANOVA on expression data. RNA-expression levels were standardised by the geometric mean of three genes (Actin, Ubiquitin and GB19767), whose expression levels were relatively consistent. Altogether we analysed 74 RNA samples, each isolated from an individual larva. Of these, 34 larvae were collected from the three C hives, and 40 from the IE hives. We asked about the significance of variation in RNA expression levels found to exist between repeated measurements for the same larva, among hives, and between IE and C groups. Three genes, Hsp90, FAS and G6PD, respectively, were found, following the conservative Bonferroni correction [56] to exhibit marginal significance for variation among C and IE larvae. Thus, differences between measures of individual genes from C and IE groups are subtle, but overall in good agreement with the results obtained by RNA-Seq (Figure 3; Data S3). Given that even genetically identical organisms can have pronounced fluctuations in the expression of individual genes [57], we were not surprised to detect varying expression levels among individual larvae from heterogeneous, wild-type bee populations. The genetic background – reflected by hive-impact – contributes to variations in expression levels of some genes (i.e. Cyp6/GB19113 in Data S3).

Changes in the Composition of Lipid Compounds Analysis of our RNA-Seq data point towards altered lipid metabolism in imidacloprid-exposed larvae (Figure 1). To find possible unbiased differences in lipid composition, 20 individual lipid profiles were established from extracts of ten C and IE larvae, respectively. For this, we used liquid chromatography coupled with mass spectrometry (LC-MS) analysis. We detected a total of 1638 lipid metabolites with strong MS signals (threshold >1000 cps) that were present both in C and IE samples. Of these, 247 metabolites (15%) showed significant differences in ratios between C and IE samples (P<0.001, Student’s T) (Figure 4; Table S3). While identification of lipid species on the basis of accurate mass of the ion is not definitive, we classified around 27% of the compounds (68/247) in this list. Species from many lipid classes showed altered levels, including diacylglycerols (DAG), triacylglycerols (TAG), ceramides (CE), phosphoethanolamines (PE), phosphocholines (PC), phosphoserines (PS), phosphoinositol (PI) and free fatty acids (FFA) (Table S3). Thus, altered abundance is not restricted to lipid species directly involved with energy metabolism. Constituents of cell membranes, such as PC and PE lipid species, also differ in relative quantities between IE and C larvae. Small changes in lipid composition of cell membranes could alter the structure and function of membrane-embedded proteins [58]. Functional nAChRs – the molecular target site of neonicotinoids – reside in the plasma membrane as ligand-gated ion channels, composed of five protein-subunits [59]. In vitro studies demonstrate that drug action at nAChRs is influenced by the composition of membrane lipids [58]. Altered abundance of certain membrane lipid species may therefore be a physiological adjustment that counters nAChR-mediated effects of imidacloprid. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Figure 4. Exemplary data of high resolution LC-MS analysis of larval lipid extracts. Shown are graphic representations of summed spectra (positive electrospray mode) from a larval sample of an imidacloprid-exposed hive (A) and a larval sample of an unexposed, control hive (B). Y-axes (ion intensity) are normalised to the most intense ion species; x-axes indicate the specific mass-to-charge ratio (m/z). Smaller insets in A and B are zoomed-in portions of the graph axis scale (m/z 840–855) and provide an example of a typical change at m/z 850.587, which is elevated in imidacloprid-exposed samples. While some differences between the spectra can be readily observed, data processing enabled more detailed analysis. https://doi.org/10.1371/journal.pone.0068191.g004