BACKGROUND
because organophosphate  pesticides are widely used for the control of agricultural pests and of arthropod disease vectors commercially, residentially, and institutionally, large numbers of people throughout the world are routinely exposed to ops, many at toxic levels.
it has been estimated that there may be as many as three million poisonings and  <dig>  deaths per year from op exposure worldwide
<cit> .
all ops exert their primary effects as acetylcholinesterase  inhibitors.
both short term and long-term adverse effects of ache inhibitor exposure have been described, ranging from acutely life-threatening conditions to subtle, long term behavioral deficits.
however, the consequences of low level exposure and the nature and mechanism of persistent effects are poorly understood and, in many cases, controversial .
without intervention, high-dose, acute exposures to ops and other ache inhibitors result in death from respiratory failure; less severe exposures may cause excessive salivation, lacrimation, and urination, as well as diaphoresis, gastrointestinal motility, and emesis followed by paralysis
<cit> .
people and animals that recover from acute exposures to many of these compounds may present with a delayed syndrome; organophosphate induced delayed polyneuropathy , characterized by the appearance of numbness, weakness, and parathesia in the limbs, and degeneration of peripheral nerve and central nervous system myelin sheaths 7– <dig> days after exposure
<cit> .
a so-called intermediate syndrome, which presents 24– <dig> h after exposure, has been identified and is characterized by weakness of the neck, proximal limb, and respiratory musculature.
this syndrome is believed to result from acetylcholine receptor desensitization
<cit> .
the occurrence of persistent neurological and neuropsychiatric effects after low level exposure and of developmental neuro-behavioral effects has also been described
<cit> .
these effects remain controversial, however,
<cit>  because of the difficulty in obtaining persuasive evidence of low-level and transient exposures to these compounds
<cit> .
in an effort to resolve some of the issues surrounding acute, transient, and low-level exposures to ops, we have undertaken studies tracking both development of and recovery from op intoxication at the global gene expression level using the genomic model organism caenorhabditis elegans and affymetrix whole genome c. elegans genechip microarrays.
c. elegans shows substantial similarity to mammals in the relevant biochemistry and genomics.
the acute toxicity of the commercial ops, including dichlorvos, results from inhibition of acetylcholinesterase in vertebrates
<cit> , and this mechanism appears to be the same in worms as well
<cit> .
in contrast to mammals, however, the nematode possesses four ache genes rather than a single alternatively spliced gene
<cit> .
mutations in at least  <dig> c. elegans genes have been shown to confer resistance to the op-like carbamate pesticide aldicarb, and all these resistance genes are known—or plausibly believed—to be involved in acetylcholine metabolism, secretion, or recycling
<cit> .
the c. elegans genome also contains a close homolog of the vertebrate secondary op target, neuropathy target esterase
<cit> .
impairment of nte function is thought to underlie opidp
<cit> .
in previous work, we found that c. elegans is relatively resistant to op lethality
<cit>  but incurred developmental abnormalities after  <dig> h in high concentrations of the ache inhibitors aldicarb , fenamiphos, and dichlorvos .
the abnormalities included malformed cuticles, severe anatomical disorganization, and protruding vulvae, suggesting that the ops also target molecules other than ache.
a large number of non-ache and non-nte targets have been proposed in humans and other vertebrates, including nicotinic, muscarinic, and cannabinoid receptors, kinases, and carboxylesterases in addition to ache and nte
<cit> .
the potent neurotoxicity of the ops, however, has made it difficult to study “off-target” effects.
we chose to perform the present studies with dichlorvos in order to maintain comparability with previous work and also because dichlorvos is known to reversibly inactivate ache
<cit> .
we reasoned that it might be possible to wash dichlorvos out of the worm culture so that the worms’ recovery from acute intoxication could be followed.
while it is unlikely that c. elegans experiences a syndrome like opidp because of the brevity of its life-span and the shortness of its neurons
<cit> , opidp has been described in humans following dichlorvos exposure
<cit> .
this raises the possibility that the consequences of nte inhibition by dichlorvos might nevertheless be observable in c. elegans at the gene expression level.
as noted above, we have observed developmental abnormalities in worms exposed to dichlorvos suggesting that there could be undiscovered targets for ops in c. elegans.
because c. elegans is relatively resistant to op lethality, yet shows substantial similarity to mammals in the principal biochemical systems targeted by ops, we conjectured that it might be possible to highlight non-cholinergic effects from op exposures in the worm that are difficult to discern using animal models less resistant to ops.
we either exposed synchronized l3/l <dig> stage c. elegans cultures to dichlorvos or washed the dichlorvos out of the culture to allow the worms to recover .
we harvested worms at several time points and performed an analysis of global gene expression.
in these studies, we have elucidated perturbations in key pathways related to energy metabolism, innate immunity, and cellular damage resulting from op exposure.
c. elegans
following dichlorvos exposure
worms were exposed to various concentrations of dichlorvos either continuously for  <dig> h or for  <dig> h followed by removal of the dichlorvos and recovery in fresh medium for  <dig> or  <dig> h. worms were examined microscopically, and the amount of movement compared with controls was recorded.
RESULTS
organophosphate  pesticides are commonly used and may account for as many as three million poisonings and  <dig>  deaths per year worldwide
<cit> .
the principal mechanism of acute toxicity, inhibition of acetylcholinesterase , has been well studied, yet little is known about the cause of “off-target” effects, including persistent and delayed neurological effects.
to investigate off-target effects, we chose the nematode, c. elegans, because of its resistance to lethality from ache inhibition.
range-finding experiments were performed to determine suitable exposure concentrations, and worms were exposed to two concentrations dichlorvos.
in a subset of the cultures, the toxicant was removed to examine recovery specific mechanisms.
gene expression changes were monitored throughout the course of the experiment using whole genome microarrays.
a striking consequence of the exposure was delayed development , which complicated analysis because of the interplay between developmentally-regulated events and toxicant-induced expression changes.
in general, the worms showed time-dependent responses to the toxicant, which progressively intensified.
the major effects of the exposure included metabolic disruption, and innate immunity-like and cellular repair responses.
effects of dichlorvos on worm motility
we performed range-finding experiments by observing the effects of varying concentrations of dichlorvos on motility since impairment of worm movement is a readily observable consequence of op exposure in worms
<cit> .
for the definitive experiments, we chose concentrations of dichlorvos that produced observable effects on motility yet allowed the worms to recover following washout.
commencing at the early l <dig> stage , worms were exposed for  <dig> h to a series of 5-fold dilutions of dichlorvos beginning with  <dig> μm, a concentration with effects on gene expression at  <dig> h of exposure, as we have previously described
<cit> .
at the end of the exposure, the worms cultured with  <dig>  and  <dig>  μm dichlorvos were immobile.
immobility was associated with the worms assuming a static form with multiple sharp bends , consistent with hypercontracted paralysis resulting from excess acetylcholine at the neuromuscular junction
<cit>  and resembling the sharp bends observed by mahoney et al.
in worms exposed to aldicarb, a carbamate inhibitor of acetyl cholinesterase
<cit> .
although the rate at which worms writhed when exposed to  <dig>  μm dichlorvos was visibly slowed at  <dig> h post-exposure, few, if any, worms became sharply bent or completely immobile even after  <dig> h of exposure.
there was no obvious difference between control cultures and cultures exposed to  <dig>  μm even after  <dig> h of exposure.
these observations are consistent with those of cole and coworkers
<cit> , who developed a software package to quantify worm mobility which demonstrated that a  <dig> h exposure to  <dig>  μm dichlorvos reduced worm motion by 50%.
to verify the reversibility of the effects of dichlorvos exposure, time-course experiments were performed in which early l <dig>  worms were exposed to varying concentrations of dichlorvos for  <dig> h, and then the culture medium was changed to remove the dichlorvos.
by visual inspection, worms exposed to  <dig>  μm or  <dig>  μm dichlorvos showed an incomplete recovery over  <dig> h. in contrast, worms exposed to  <dig>  μm dichlorvos appeared to be fully recovered by  <dig> h. by  <dig> h, we could not readily distinguish the behavior of the worms in any of the exposures from control.
there was no discernible effect of  <dig>  μm dichlorvos on the worms under any condition.
although the worms resumed normal behavior based on visual observations, ache activity may still have been inhibited.
in studies of the recovery of adult c. elegans from a  <dig> h exposure to the ache inhibitory pesticides carbofuran and fenamiphos, the recovery of ache activity was incomplete even when behavioral evidence indicated complete recovery
<cit> .
gene expression analysis of dichlorvos exposure
for the definitive study, worms were exposed to two concentrations of dichlorvos;  <dig>  μm—the loel  determined in the motility assays—and  <dig> μm—a concentration previously shown to elicit changes in gene expression
<cit> .
the exposures began at the l3/l <dig> stage cuticular molt  and continued for  <dig> h, which permitted a full  <dig> h recovery after an  <dig> h exposure but did not extend into the egg-laying period .
three exposure regimes were utilized .
in one, there were no changes in the medium prior to harvesting the samples.
in the other two, the medium was washed out and replaced with fresh medium after  <dig> or  <dig> h to remove the dichlorvos, allowing the worms to recover; control flasks were also subjected to the washout procedure .
samples were harvested at the zero time point,  <dig> h, and then every  <dig> h until the conclusion of the experiment at  <dig> h. the entire exposure protocol was completed four times on separate days to be used as biological replicates for the microarray analysis.
levels of expression were determined for approximately  <dig>  genes represented on the whole genome microarray for each of the  <dig> samples; two samples failed to provide data that passed quality control standards.
a total of  <dig>  probe sets were detected above background levels in all the replicates of at least one condition.
we considered only these probe sets in subsequent analyses.
in the initial evaluation of the data, it was clear that the  <dig> h washout regime and the low dose exposures produced few detectable effects; therefore, we focused our investigation on the other conditions.
dichlorvos-dependent delayed development
previous work
<cit>  and preliminary data exploration clearly indicated that exposure to dichlorvos causes a delay of development in the l <dig> larvae as judged by gene expression.
a principal components analysis  of the samples from the continuous exposure protocol  shows that the largest effect in this experiment is associated with the time of harvest, presumably resulting from changes in gene expression throughout normal development.
the spread of the clusters of worms harvested at different times in the high concentration exposure regime is compressed along the first principal component  toward the early harvest times for the unexposed samples, illustrating the dichlorvos-dependent developmental delay that complicates the analysis of these data.
the predicted age of worms cultured under the continuous or the  <dig> h washout conditions described in the text was calculated as described in methods; the age is calculated from the initiation of the exposure.
the number of probe sets used for predicting each age is shown, and predicted ages are reported ± sem.
statistically significant differences from control are indicated by asterisks.
the data presented in this table are also plotted in figure
<dig>
in an attempt to better understand the effects of exposure on development and to predict the amount of developmental delay, we developed a linear regression-based method for computing a predicted age  for comparison purposes based on the expression level of developmentally-regulated genes and the time of harvest .
for the pa, the beginning of the exposure is set as the zero time.
as expected, the unexposed samples in both the continuous regime and washout regime have pas that are similar to their harvest times .
the washout procedure itself somewhat delayed development, but the effect is only significant in the  <dig> h-harvest condition.
however, the exposed samples show significant differences from the harvest time for all times 8– <dig> h post-exposure.
in the  <dig> h- and  <dig> h-harvest samples, there is approximately a 5– <dig> h delay in development in the exposed worms.
surprisingly, the worms in the continuous exposure and  <dig> h washout protocols showed similar developmental delays , suggesting that some key developmental process in early l <dig> larvae are susceptible to interference by dichlorvos or that the cause of the delay persists at least for a time after the dichlorvos is removed.
while the profiles of the majority of the developmentally regulated genes are consistent with a well-regulated but delayed developmental process, we identified  <dig> developmentally-regulated genes  that do not appear to be affected by the developmental delay.
these genes were identified by filtering the list of  <dig>  developmentally regulated genes  using the following inclusion criteria applied to data from the continuous exposure protocol: 1) a ≤  <dig> -fold difference in expression level between exposed and control conditions at each time point; 2) a ≥  <dig> -fold difference between at least one pair of consecutive harvest times in the 8– <dig> h time range.
the only significantly enriched biological classification identified for this set was the presence of two genes that encode proteins involved in amino acid transport and metabolism .
based on the limited annotation available for these genes, we cannot provide a clear biological interpretation of the role this small set of genes plays in development.
it is also possible that the genes are regulated by both dichlorvos exposure and development in such a way that some interaction between the two makes it appear that the expression profile of these genes is not delayed.
dichlorvos exposure significantly delayed worm development—by more than  <dig> h at the  <dig> h-harvest point as judged by the pa value.
overall, the coordinated progression of development does not appear to be disrupted, but we are unable to rule out desynchronization of a small number of developmental processes.
biological responses to dichlorvos exposure
in order to understand the mechanism of dichlorvos toxicity, we attempted to identify gene expression changes that were associated with dichlorvos exposure but were not the consequence of developmental delay.
we initially limited our analysis to the continuous exposure regime and compared data from high concentration samples to control samples.
differentially expressed probe sets were selected using a p < 10− <dig>  and  <dig> -fold change-filter based on an analysis of covariance  that included pa  and exposure as factors.
a total of  <dig>  probe sets was identified across the five harvest times, with increasing numbers at later time points .
the expression levels of these genes increasingly diverged from control over time, but as expected, following washout, their expression rapidly returned to control levels .
to determine the biological processes affected by exposure, we compared this list of differentially expressed genes to lists of genes or probe sets from several data sources including gene ontology terms and microarray expression clusters from wormbase, and transcription factor binding sites from the literature
<cit>  .
the major processes identified in these analyses are metabolic regulation, innate immunity, stress response, and muscular and neurological regeneration.
daf-16-mediated stress responses
a large number of the genes differentially expressed in reaction to dichlorvos exposure are transcriptionally regulated by daf- <dig>  a forkhead transcription factor implicated in the control of metabolism, ageing, innate immunity, and stress responses in c. elegans.
the microarray cluster enrichment analysis in particular revealed a large over-representation of genes known to be controlled by daf- <dig>
<cit> .
murphy et al.
<cit>  identified genes downstream of daf- <dig> in microarray experiments using adult daf-2; daf- <dig> hypomorphic mutant worms or worms in which daf- <dig> and daf- <dig> expression was knocked down with rnai.
in our work, we observed changes in expression of these downstream genes that are largely consistent with increasing daf- <dig> activity over the course of the exposure .
<cit> .
changes in gene expression  are shown for continuous and 8-h washout  conditions.
b) graph of ins- <dig> expression.
the means by which daf- <dig> is being regulated cannot be directly inferred from the gene expression data, but it appears to be complex.
the best studied mechanism of daf- <dig> regulation is repression through a phosphorylation cascade initiated by the binding of insulin-like ligands to the daf- <dig> receptor
<cit> .
however, we observed that the ins- <dig> gene, which encodes a known agonist of daf- <dig>  is up-regulated in response to dichlorvos exposure  which should result in the suppression of daf- <dig> activity rather than the inferred activation.
this observation suggests that there is a signal activating daf-16-mediated transcription that is independent of daf- <dig> and possibly countering the normal repression elicited by ins-7/daf- <dig> binding.
a likely mechanism for activation of daf- <dig> is evident from the elevated adp/atp ratios that have been documented as a result of organophosphate exposure
<cit> .
it is possible that one of the amp-kinases—which are sensitive to adp/atp balance—is transducing this energy metabolism signal to daf- <dig>  a known ampk target
<cit> .
consistent with this possibility, we found that genes known to be up-regulated in crh- <dig> mutants, an ortholog of the camp-inducible creb transcription factor and downstream target of ampk signaling
<cit> , were over-represented in the microarray cluster analysis  and primarily up-regulated in response to dichlorvos.
aak- <dig> is a strong candidate for the ampk that could be activating daf- <dig>  as genes downstream of aak- <dig> are over-represented in the dichlorvos data set ; furthermore, the gene encoding superoxide dismutase, sod- <dig>  a known daf- <dig> target, is down-regulated in aak- <dig> mutants
<cit> .
the amp dependent kinase, kin- <dig> is a less likely candidate for activating daf- <dig>  in response to dichlorvos exposure, the message levels for kin- <dig> are up-regulated.
further, it has been shown to inactivate torc2
<cit> , an essential transducer in the daf- <dig> signaling pathway.
although, kin- <dig> may play a role in metabolic signaling since it is involved in food sensing, it appears to be chiefly involved in regulating chemoreceptor signaling in sensory neurons
<cit> , with limited expressions elsewhere in the nervous system and in body wall muscles and hypodermal cells
<cit> .
these factors argue against kin- <dig> playing an important role in daf- <dig> regulation.
in our analysis, we have inferred that daf- <dig> is activated based on the behavior of the genes proposed as downstream targets of daf- <dig> by murphy et al.
<cit> .
however, approximately 27% of the putative daf- <dig> target genes in our study are not regulated in the direction that would be predicted from murphy and coworkers’ experiments
<cit> .
functional analysis did not reveal any specific functions for the discrepant set of genes that were distinct from the functions observed for the complete set of differentially expressed genes.
we did, however, observe that this set was enriched with genes containing pqm- <dig> binding sites , suggesting an alternative transcriptional regulator event that might be responsive to reactive oxygen stress .
in our view, the differences in results in the studies most likely arise from differences in the experimental designs, particularly because it is not possible to distinguish between direct and indirect effects, or among effects due to differential activity of the daf- <dig> isoforms
<cit> .
murphy and coworkers used adult worms fed on a bacterial lawn in a purely genetic experiment, while we studied larval worms in axenic medium exposed to a chemical toxicant.
it is known that regulatory processes involving daf- <dig> are complex, and the differences in life stage of the worm, culture conditions, and experimental perturbation could easily account for the differences we see.
we conclude that the effects we have described are largely daf- <dig> mediated in response to dichlorvos exposure.
based on our observations, we believe that the daf- <dig> mediated gene response is a result of reduced energy reserves within the cells and is not a daf- <dig> mediated event.
two known mechanisms of op toxicity can contribute to this effect: 1) elevated consumption of energy because of increased muscular activity due to continual stimulation when acetylcholine is not degraded and 2) cholinergic receptor-independent mitochondrial disruption
<cit> .
furthermore, in other work, worms exposed to similar dichlorvos concentrations show reduced bacterial feeding
<cit> .
these effects, which are likely additive, appear to be driving forces in daf- <dig> mediated gene expression changes affecting metabolic processes.
innate immunity related responses
we found a number of similarities between responses elicited by dichlorvos and infection in c. elegans.
dichlorvos exposure appears to activate some innate immunity signaling cascades , and microarray cluster enrichment analysis  revealed that many genes affected by dichlorvos exposure are in data sets derived from worms infected with a diverse set of pathogens and show similar regulation.
however, dichlorvos exposure does not fully mimic the responses to infection.
microbacterium nematophilum, enterococcus faecalis, photorhabdus luminescens, yersinia pestis, and pseudomonas aeruginosa infections activate many of the same genes as dichlorvos, and show similar regulation.
many of these genes are known to be involved in innate immunity and were up-regulated in response to the exposure, including genes encoding lysozymes , invertebrate lysozymes , antibacterial factors , caenacin antimicrobial peptides  and  <dig> c-type lectins.
the response appears to overlap but to be distinct from innate immunity because other genes considered to part of the innate immune response are down-regulated, including lys- <dig>  lys- <dig>  lys- <dig>  lys- <dig>  ilys- <dig>  and  <dig> c-type lectins genes.
dichlorvos exposure also impacts three of the seven signaling cascades known to be involved in immune defense in c. elegans <cit> : the daf- <dig> insulin-like receptor pathway , the sma/tgb-β pathway, and the p38/pmk- <dig> mapk pathway.
the components of the sma/tgb-β pathway—the tgf-β receptor , sma − <dig>  and sma- <dig> —are all up-regulated.
the effects of dichlorvos on the p38/pmk- <dig> mapk and related pathways are complicated.
we observed up-regulation of the toll-domain containing receptor  upstream of p38/pmk- <dig> along with sek- <dig>  the map kinase kinase regulator of pmk- <dig>  and pmk- <dig>  we also found enrichment for genes with a binding site for the gata-family transcription factor, elt- <dig>  which has been shown to affect pmk- <dig> mediated expression of antimicrobial peptides
<cit> .
the expression of pmk- <dig> itself is regulated by npr- <dig>  a homolog of the neuropeptide y receptor, and we find that  <dig> of the  <dig> npr- <dig> regulated genes
<cit>  present in our array data  are differentially expressed.
however, the direction of gene change induced by dichlorvos for pmk- <dig> responsive genes and for npr- <dig> regulated genes is primarily in the opposite direction from those predicted for pmk- <dig> mediated innate immunity
<cit> .
a potential reason for this discrepancy is the increased expression of the gene encoding the phosphatase vhp- <dig>  which is a known antagonist of pmk-1
<cit> .
vhp- <dig> is a dual-specificity phosphatase that interacts with multiple map kinases including kgb- <dig>  pmk- <dig>  and pmk-1
<cit>  and plays a role in modulating the respective signaling pathways during response to stress.
it has been shown to be regulated by rnt- <dig> during oxidative stress
<cit> .
it is likely that the increased expression of vhp- <dig> and the concurrent repression of the pmk- <dig> pathway are the results of dichlorvos-induced stress.
the role of the daf- <dig> insulin-like pathway in the immune response triggered by dichlorvos is somewhat complex.
as mentioned above, the daf- <dig> pathway is presumably being activated by ins- <dig>  however, a secondary signal, likely via an ampk, is preventing daf- <dig> mediated suppression of daf- <dig> activity.
what makes this more intriguing is that the induction of ins- <dig> is required for p. aeruginosa virulence
<cit> , presumably due to suppression of the immune response.
three other genes encoding insulin-like peptides  are also up-regulated and one  is down-regulated, although all are affected to a much lesser extent than ins- <dig>  with the variety of signaling molecules and control mechanisms acting on the daf- <dig> pathway, it is difficult to determine whether the observed effects are part of an immune response or represent another role of the daf- <dig> signaling pathway.
it is clear that the pattern of expression for some of these genes overlaps with the immune response.
there has been some discussion in the community as to whether the c. elegans’ immune response is specifically triggered by pathogenic signaling molecules or if it is a generalized stress response
<cit> , since no homologs for several key components of the mammalian toll-like receptor pathway have been identified in the worm.
under the axenic culture conditions utilized in this work, the differentially expressed genes that are involved in the immune response are not being induced by a pathogen.
however, it is not clear whether the dichlorvos stimulation of the innate immunity response acts through a “normal” stress sensing process or whether it is affected by an aberrant “toxic” mechanism.
there are four primary mechanisms by which dichlorvos could be triggering the innate immune response: generalized stress, altered feeding pathways, autogenic stimulation, or neuronal signaling.
the simplest explanation for stimulation of the immune response is through stress pathways induced by the dichlorvos exposure, but others must also be considered.
since bacteria are the primary food source for c. elegans under natural conditions, there are overlaps between the immune response to bacteria and adaptation to different bacterial food sources
<cit> .
therefore, the adaptive response triggered by the low energy state within the worms might induce genes necessary for feeding on different bacteria, and these genes may contribute to the observed immune response.
another possibility has recently been described in mammals.
molecules present in organelles derived from ancient prokaryote symbiotes have characteristics that cause them to trigger immune receptors in the same manner as infectious bacteria
<cit> .
in light of the potential damage to muscle and nerves , the release of these molecules is likely to be occurring and might be the cause of the immune response induced by dichlorvos exposure; though it is still unclear whether c. elegans responds to these molecules in this way
<cit> .
finally, the primary mechanism of dichlorvos toxicity, the inhibition of ache, may disrupt the neuronal signaling that helps to control the innate immune response
<cit> .
it is reasonable to assume that with the variety of genes and processes that comprise the immune response that the stimulation of this response seen in exposed worms includes a combination, if not all, of these processes.
neural and musculature repair responses
the gene expression changes seen at the later time points  highlight the known myopathic and neurodegenerative effects of organophosphate poisoning
<cit>  and are consistent with increased muscular and neuronal growth and repair.
there was an increase in expression levels for a large number of muscle-specific genes that are required for myogenesis.
genes encoding proteins localized to the m-bands  are over represented in a gene ontology analysis;  <dig> of the  <dig> genes present in the data set that encode m band proteins are up-regulated.
other muscle genes up-regulated include those encoding titin , two paralogs of titin , homologs of troponin t and c , and the homolog of tropomyosin .
the evidence for neuronal growth and repair is most apparent through effects that are related to axon guidance and regeneration.
two genes in the dlk- <dig> mapk pathway, dlk- <dig>  a mapkkk, and pmk- <dig>  a p <dig> mapk, were up-regulated.
the dlk- <dig> pathway is essential for axon regeneration
<cit> .
there was also increased expression of multiple genes required for axonal guidance .
in addition, two of the pathways that are involved in innate immunity play a dual role in the worm and are also involved in neuronal cell fate determination.
the dbl-1/sma tgf-β pathway has been implicated in controlling neuronal cell fate in the male sensory array
<cit> , and the proteins encoded by two of the genes up-regulated in the pmk- <dig> mapk pathway, unc- <dig> and sek- <dig>  function to determine asymmetric cell fate decision during sensory neuron development
<cit> .
it is not possible to conclude from the data in-hand whether or not the changes in these pathways are the result of innate immunity or cell fate determination, or if they may be involved in both.
taken together, these gene expression changes suggest that at the later time points in the exposure, the worms have significant cellular damage and have initiated multiple repair processes.
the mechanism by which dichlorvos provokes myopathic and neurodegenerative effects has not been defined, but some work suggests that it is related to mitochondrial dysfunction
<cit> .
the changes in metabolic processes and potential mitochondrial dysfunction  observed in response to dichlorvos exposure make this a plausible mechanism for the damage induced in the worms.
early responses
in order to focus on the primary effects of dichlorvos intoxication, we closely examined the  <dig> genes that were differentially expressed within the first  <dig> h of the exposure .
the most notable biological processes affected at this early time point are immune response and energy metabolism.
in the microarray cluster enrichment analysis, this was highlighted by over representation of genes regulated by daf- <dig> or differentially expressed during infection by p. aeruginosa .
the down-regulation of the genes encoding a fatty acid elongase  and an acetoacetyl-coa synthetase  suggest a reduction in fatty acid elongation, while the up-regulation of the gene encoding an acyl-coa synthetase  implies an increase in β-oxidation.
<cit> , ox_stress
<cit>  and also for time after infection for pa <dig> vs op50
<cit> .
pqm- <dig> and pha-4_starved are list of genes with binding sites for the respective transcription factors
<cit> .
cells marked with yellow indicate dichlorvos-responsive genes that are present in the gene set indicated at the head of the table.
the numbers at the foot indicate the number of genes shared between the  <dig> h dichlorvos exposure gene set and the indicated experimental data.
the direction of change in gene expression is shown in the left-most column, minus sign/green—down-regulated, plus sign/red—up-regulated.
over half of the genes differentially expressed at this time have a binding site for pqm- <dig> , a transcription factor strongly up-regulated by paraquat exposure, and seven of the genes were identified as being up-regulated by oxidative stress .
the primary mechanism of paraquat toxicity is believed to be reactive oxygen species  generation mediated by mitochondria
<cit> .
taken together, these observations suggest that the early effects of dichlorvos exposure include a reduction of available energy and increased ros.
while increases in β-oxidation contribute to ros generation, the gene response suggests a higher level of ros than would be anticipated by β-oxidation alone.
a lower energy state and higher levels of ros are consistent with the idea that dichlorvos causes mitochondrial dysfunction.
the reason for overlap with genes differentially regulated during p. aeruginosa infection is not as directly evident.
there are no bacteria present, and it is unlikely to be autogenic due to damaged cells.
the two most plausible reasons are 1) a stress-mediated immune response, or 2) overlap with non-immune mediated daf- <dig> insulin signaling.
the latter is of particular relevance considering one of the pathogenic mechanisms of p. aeruginosa is the up-regulation of ins- <dig> <cit>  and while ins- <dig> was not differentially expressed at  <dig> h based on our strict criteria, it is in fact modestly up-regulated  and strongly up-regulated at later time points .
therefore, the genes regulated by p. aeruginosa infection are likely to include genes regulated by the daf- <dig> insulin signaling pathway both in an immune-specific and non-specific manner.
the majority of the genes differentially expressed as an early response to dichlorvos is directly related to energy metabolism or is regulated by the daf- <dig> insulin-signaling pathway.
these effects highlight a change in metabolism within the worms that is indicative of reduction in available energy and the utilization of energy reserves.
the inclusion of genes regulated by oxidative stress and pqm- <dig> hints at mitochondrial dysfunction as the mechanism for the reduction in energy levels.
recovery
an important goal of this study was to identify genes and biological processes that might be involved in recovery from dichlorvos intoxication.
to identify genes that participate in the recovery process and not ongoing toxicity, data were examined from unexposed and high concentration samples collected after the  <dig> h washout .
we used the same criteria for selecting differentially expressed genes in the washout conditions  that we used for the continuous exposure analysis.
a total of  <dig> probe sets meet these criteria for differential expression in at least one of the harvest times.
the effect of the exposure diminishes rapidly with time as there are  <dig> probe sets differentially expressed in the  <dig> h washout samples and only  <dig> at the  <dig> h time point.
there were no significant differences at the  <dig> h time point.
as expected, many of these probes are also differentially expressed in the continuous exposure regime .
an interesting observation is that  <dig> probe sets that were not identified as being differentially expressed at the  <dig> h time point are differentially expressed at later times by the worms in both the continuous exposure and washout protocols, although the magnitude of differential expression is lower for nearly every probe set in the washout samples than the continuously exposed samples.
while trending suggests that some of these genes are affected at  <dig> h , there is a set of genes whose expression is clearly affected only after the washout has occurred.
this observation suggests that there are processes affected by dichlorvos which produce delayed effects even after removal of the toxicant.
we were unable to find any statistically enriched processes related to these genes, but many appear to be involved in immune response and energy homeostasis.
the probe sets differentially expressed in the washout samples, but not the continuous exposure regime, are the most attractive targets for identifying recovery-specific effects.
there are seven probe sets that were identified only in the washout samples.
however, because the selection criteria were designed to stringently minimize false positives, it is possible that some of them actually could be differentially expressed during the continuous exposure regime.
we used less stringent criteria  to allow us to discern probe sets that were weakly differentially expressed in the continuous exposure regime and found that only two of the seven probe sets were uniquely differentially expressed after the washout procedure .
unfortunately, the proteins encoded by the genes targeted by these probe sets, a widely expressed heat shock protein  and an iron-binding protein of unknown function , offer little insight into recovery specific mechanisms.
in addition to looking at the genes that were uniquely identified in the recovery period, we examined the entire group of probe sets that were differentially expressed during the recovery phase to determine biological mechanisms potentially involved in recovery.
this list was subdivided to include only the  <dig> probe sets that were differentially expressed using strict criteria after the dichlorvos was removed .
in both the complete list and the “after exposure” list, the biological processes that are represented appear to exactly mimic those that are seen during the exposure, albeit to a lesser extent .
no recovery-specific mechanisms were identified.
after removal of dichlorvos from the medium, the gene expression changes seen in the worms appear to return to baseline levels rather rapidly, and those observed at the first harvest after the washout simply appear to be residual effects of ongoing biological processes resulting from the exposure.
by the  <dig> h harvest time, the study lacks sufficient power to identify any differentially expressed genes.
with the exception of expression changes that appear to be the result of delayed development, the treated worms appear to have fully recovered from the exposure.
mechanism of toxicity
many if not all of the biological effects identified in this work that were caused by dichlorvos exposure can be plausibly tied to inhibition of acetylcholinesterase.
developmental timing in c. elegans requires active nicotinic acetylcholine receptors
<cit> .
persistent cholinergic signaling causes a depletion of muscular atp that can mediate an increase in certain metabolic processes
<cit> .
acetylcholine signaling during starvation activates a mapk  in pharyngeal muscles
<cit> .
innate immunity in c. elegans is known to be at least partially regulated by neuronal cells and signaling
<cit> , and this signaling leads to the activation of pmk-1
<cit> .
from these examples, it is clear that neuronal signaling and acetylcholine levels can affect metabolic processes, innate immunity, and developmental time, which are the primary effects observed in response to dichlorvos exposure.
though ache inhibition is likely the dominant cause of toxicity in c. elegans, there are two other non-cholinergic mechanisms that have been implicated in dichlorvos toxicity that are consistent with the observations in this work; they are mitochondrial dysfunction
<cit>  and perturbation of calcium homeostasis
<cit> .
the metabolic changes highlighted by alterations in expression levels of daf- <dig> target genes underscore the lower energy reserves present in exposed worms.
some of this deficiency may be caused by increased muscle stimulation or reduced feeding, but the high levels of induction suggest a more extreme reason for this deficiency, such as a reduction in energy production triggered by the effects of dichlorvos on mitochondrial function.
the down-regulation of the genes encoding multiple subunits of the mitochondrial electron transport chain including ones from complex i , complex ii , and complex v  suggests a reduction in mitochondrial energy production.
further work will be required to delineate what roles mitochondrial dysfunction is playing in dichlorvos toxicity.
the role of calcium homeostasis in dichlorvos toxicity is another area where the data do not provide a clear answer but are very suggestive.
there are numerous calcium dependent processes that are perturbed by the exposure.
for example, the gene encoding unc- <dig>  a calmodulin-dependent protein kinase that can activate the p38/pmk- <dig> mapk pathway, is up-regulated and could be the cause of the activation of that pathway seen in this work.
there is also evidence for activation of a calcineurin pathway.
the gene encoding cnb- <dig>  an ortholog of calcineurin b and rcn- <dig>  a gene whose transcription is dependent on a calcineurin a ortholog  and free calcium, are both up-regulated.
furthermore, genes known to be regulated by tax- <dig> are over represented in the data set .
the exact role of calcium homeostasis in dichlorvos toxicity is unknown, and there is no defined mechanism by which dichlorvos effects calcium homeostasis.
this is another promising area where more work must be performed to fully understand the effects of dichlorvos toxicity.
CONCLUSIONS
upon exposure to dichlorvos, worms show a time-dependent response with progressive delays in development, alterations in metabolism, increased cellular repair, and alterations in innate immunity functions.
early in the course of the exposure, the effects appear to be primarily related to energy metabolism.
as the exposure progresses, the changes in the metabolic pathways become more pronounced, and there is increased perturbation of immune responses.
at the latest times, repair processes have been activated that are presumably a response to damaged nerve and muscles cells.
upon removal of the toxicant, gene expression quickly returns to baseline, with the only lingering effect being a synchronized delay in developmental gene expression.
this delay in development is the single largest effect on the worms and is evident as early as  <dig> h after treatment.
there are several potential mechanisms by which dichlorvos is mediating toxicity in c. elegans.
inhibition of ache is likely a key player as many of the perturbed processes involve neuronal signaling, and the observed behavior changes are consistent with increased accumulation of acetylcholine.
the observed changes in the expression of genes involved in energy metabolism suggest that dichlorvos is disrupting mitochondrial function, which has been shown to be a mechanism for neuronal degeneration
<cit> .
calcium homeostasis is potentially another process perturbed by dichlorvos exposure, but it is unclear whether dichlorvos directly or indirectly perturbs calcium homeostasis.
through this work, we have identified key genes and pathways responding to dichlorvos toxicity .
the primary mechanisms appear to be inhibition of ache and depletion of energy reserves possibly through mitochondrial dysfunction.
we have also shown the resilience of the worms by their recovery after removal of the toxicant.
this work provides insight into the detailed process of dichlorvos intoxication and sets the groundwork for future efforts.

 because organophosphate pesticides are widely used for the control of agricultural pests and of arthropod disease vectors commercially, residentially, and institutionally, large numbers of people throughout the world are routinely exposed to ops, many at toxic levels. it has been estimated that there may be as many as three million poisonings and deaths per year from op exposure worldwide. all ops exert their primary effects as acetylcholinesterase inhibitors. both short term and long-term adverse effects of ache inhibitor exposure have been described, ranging from acutely life-threatening conditions to subtle, long term behavioral deficits. however, the consequences of low level exposure and the nature and mechanism of persistent effects are poorly understood and, in many cases, controversial. without intervention, high-dose, acute exposures to ops and other ache inhibitors result in death from respiratory failure; less severe exposures may cause excessive salivation, lacrimation, and urination, as well as diaphoresis, gastrointestinal motility, and emesis followed by paralysis. people and animals that recover from acute exposures to many of these compounds may present with a delayed syndrome; organophosphate induced delayed polyneuropathy, characterized by the appearance of numbness, weakness, and parathesia in the limbs, and degeneration of peripheral nerve and central nervous system myelin sheaths 7– days after exposure. a so-called intermediate syndrome, which presents 24– h after exposure, has been identified and is characterized by weakness of the neck, proximal limb, and respiratory musculature. this syndrome is believed to result from acetylcholine receptor desensitization. the occurrence of persistent neurological and neuropsychiatric effects after low level exposure and of developmental neuro-behavioral effects has also been described. these effects remain controversial, however, because of the difficulty in obtaining persuasive evidence of low-level and transient exposures to these compounds. in an effort to resolve some of the issues surrounding acute, transient, and low-level exposures to ops, we have undertaken studies tracking both development of and recovery from op intoxication at the global gene expression level using the genomic model organism caenorhabditis elegans and affymetrix whole genome c. elegans genechip microarrays. c. elegans shows substantial similarity to mammals in the relevant biochemistry and genomics. the acute toxicity of the commercial ops, including dichlorvos, results from inhibition of acetylcholinesterase in vertebrates, and this mechanism appears to be the same in worms as well. in contrast to mammals, however, the nematode possesses four ache genes rather than a single alternatively spliced gene. mutations in at least c. elegans genes have been shown to confer resistance to the op-like carbamate pesticide aldicarb, and all these resistance genes are known—or plausibly believed—to be involved in acetylcholine metabolism, secretion, or recycling. the c. elegans genome also contains a close homolog of the vertebrate secondary op target, neuropathy target esterase. impairment of nte function is thought to underlie opidp. in previous work, we found that c. elegans is relatively resistant to op lethality but incurred developmental abnormalities after h in high concentrations of the ache inhibitors aldicarb, fenamiphos, and dichlorvos. the abnormalities included malformed cuticles, severe anatomical disorganization, and protruding vulvae, suggesting that the ops also target molecules other than ache. a large number of non-ache and non-nte targets have been proposed in humans and other vertebrates, including nicotinic, muscarinic, and cannabinoid receptors, kinases, and carboxylesterases in addition to ache and nte. the potent neurotoxicity of the ops, however, has made it difficult to study “off-target” effects. we chose to perform the present studies with dichlorvos in order to maintain comparability with previous work and also because dichlorvos is known to reversibly inactivate ache. we reasoned that it might be possible to wash dichlorvos out of the worm culture so that the worms’ recovery from acute intoxication could be followed. while it is unlikely that c. elegans experiences a syndrome like opidp because of the brevity of its life-span and the shortness of its neurons, opidp has been described in humans following dichlorvos exposure. this raises the possibility that the consequences of nte inhibition by dichlorvos might nevertheless be observable in c. elegans at the gene expression level. as noted above, we have observed developmental abnormalities in worms exposed to dichlorvos suggesting that there could be undiscovered targets for ops in c. elegans. because c. elegans is relatively resistant to op lethality, yet shows substantial similarity to mammals in the principal biochemical systems targeted by ops, we conjectured that it might be possible to highlight non-cholinergic effects from op exposures in the worm that are difficult to discern using animal models less resistant to ops. we either exposed synchronized l3/l stage c. elegans cultures to dichlorvos or washed the dichlorvos out of the culture to allow the worms to recover. we harvested worms at several time points and performed an analysis of global gene expression. in these studies, we have elucidated perturbations in key pathways related to energy metabolism, innate immunity, and cellular damage resulting from op exposure. c. elegans following dichlorvos exposure worms were exposed to various concentrations of dichlorvos either continuously for h or for h followed by removal of the dichlorvos and recovery in fresh medium for or h. worms were examined microscopically, and the amount of movement compared with controls was recorded.  organophosphate pesticides are commonly used and may account for as many as three million poisonings and deaths per year worldwide. the principal mechanism of acute toxicity, inhibition of acetylcholinesterase, has been well studied, yet little is known about the cause of “off-target” effects, including persistent and delayed neurological effects. to investigate off-target effects, we chose the nematode, c. elegans, because of its resistance to lethality from ache inhibition. range-finding experiments were performed to determine suitable exposure concentrations, and worms were exposed to two concentrations dichlorvos. in a subset of the cultures, the toxicant was removed to examine recovery specific mechanisms. gene expression changes were monitored throughout the course of the experiment using whole genome microarrays. a striking consequence of the exposure was delayed development, which complicated analysis because of the interplay between developmentally-regulated events and toxicant-induced expression changes. in general, the worms showed time-dependent responses to the toxicant, which progressively intensified. the major effects of the exposure included metabolic disruption, and innate immunity-like and cellular repair responses. effects of dichlorvos on worm motility we performed range-finding experiments by observing the effects of varying concentrations of dichlorvos on motility since impairment of worm movement is a readily observable consequence of op exposure in worms. for the definitive experiments, we chose concentrations of dichlorvos that produced observable effects on motility yet allowed the worms to recover following washout. commencing at the early l stage, worms were exposed for h to a series of 5-fold dilutions of dichlorvos beginning with μm, a concentration with effects on gene expression at h of exposure, as we have previously described. at the end of the exposure, the worms cultured with and μm dichlorvos were immobile. immobility was associated with the worms assuming a static form with multiple sharp bends, consistent with hypercontracted paralysis resulting from excess acetylcholine at the neuromuscular junction and resembling the sharp bends observed by mahoney et al. in worms exposed to aldicarb, a carbamate inhibitor of acetyl cholinesterase. although the rate at which worms writhed when exposed to μm dichlorvos was visibly slowed at h post-exposure, few, if any, worms became sharply bent or completely immobile even after h of exposure. there was no obvious difference between control cultures and cultures exposed to μm even after h of exposure. these observations are consistent with those of cole and coworkers, who developed a software package to quantify worm mobility which demonstrated that a h exposure to μm dichlorvos reduced worm motion by 50%. to verify the reversibility of the effects of dichlorvos exposure, time-course experiments were performed in which early l worms were exposed to varying concentrations of dichlorvos for h, and then the culture medium was changed to remove the dichlorvos. by visual inspection, worms exposed to μm or μm dichlorvos showed an incomplete recovery over h. in contrast, worms exposed to μm dichlorvos appeared to be fully recovered by h. by h, we could not readily distinguish the behavior of the worms in any of the exposures from control. there was no discernible effect of μm dichlorvos on the worms under any condition. although the worms resumed normal behavior based on visual observations, ache activity may still have been inhibited. in studies of the recovery of adult c. elegans from a h exposure to the ache inhibitory pesticides carbofuran and fenamiphos, the recovery of ache activity was incomplete even when behavioral evidence indicated complete recovery. gene expression analysis of dichlorvos exposure for the definitive study, worms were exposed to two concentrations of dichlorvos; μm—the loel determined in the motility assays—and μm—a concentration previously shown to elicit changes in gene expression. the exposures began at the l3/l stage cuticular molt and continued for h, which permitted a full h recovery after an h exposure but did not extend into the egg-laying period. three exposure regimes were utilized. in one, there were no changes in the medium prior to harvesting the samples. in the other two, the medium was washed out and replaced with fresh medium after or h to remove the dichlorvos, allowing the worms to recover; control flasks were also subjected to the washout procedure. samples were harvested at the zero time point, h, and then every h until the conclusion of the experiment at h. the entire exposure protocol was completed four times on separate days to be used as biological replicates for the microarray analysis. levels of expression were determined for approximately genes represented on the whole genome microarray for each of the samples; two samples failed to provide data that passed quality control standards. a total of probe sets were detected above background levels in all the replicates of at least one condition. we considered only these probe sets in subsequent analyses. in the initial evaluation of the data, it was clear that the h washout regime and the low dose exposures produced few detectable effects; therefore, we focused our investigation on the other conditions. dichlorvos-dependent delayed development previous work and preliminary data exploration clearly indicated that exposure to dichlorvos causes a delay of development in the l larvae as judged by gene expression. a principal components analysis of the samples from the continuous exposure protocol shows that the largest effect in this experiment is associated with the time of harvest, presumably resulting from changes in gene expression throughout normal development. the spread of the clusters of worms harvested at different times in the high concentration exposure regime is compressed along the first principal component toward the early harvest times for the unexposed samples, illustrating the dichlorvos-dependent developmental delay that complicates the analysis of these data. the predicted age of worms cultured under the continuous or the h washout conditions described in the text was calculated as described in methods; the age is calculated from the initiation of the exposure. the number of probe sets used for predicting each age is shown, and predicted ages are reported ± sem. statistically significant differences from control are indicated by asterisks. the data presented in this table are also plotted in figure  in an attempt to better understand the effects of exposure on development and to predict the amount of developmental delay, we developed a linear regression-based method for computing a predicted age for comparison purposes based on the expression level of developmentally-regulated genes and the time of harvest. for the pa, the beginning of the exposure is set as the zero time. as expected, the unexposed samples in both the continuous regime and washout regime have pas that are similar to their harvest times. the washout procedure itself somewhat delayed development, but the effect is only significant in the h-harvest condition. however, the exposed samples show significant differences from the harvest time for all times 8– h post-exposure. in the h- and h-harvest samples, there is approximately a 5– h delay in development in the exposed worms. surprisingly, the worms in the continuous exposure and h washout protocols showed similar developmental delays, suggesting that some key developmental process in early l larvae are susceptible to interference by dichlorvos or that the cause of the delay persists at least for a time after the dichlorvos is removed. while the profiles of the majority of the developmentally regulated genes are consistent with a well-regulated but delayed developmental process, we identified developmentally-regulated genes that do not appear to be affected by the developmental delay. these genes were identified by filtering the list of developmentally regulated genes using the following inclusion criteria applied to data from the continuous exposure protocol: 1) a ≤ -fold difference in expression level between exposed and control conditions at each time point; 2) a ≥ -fold difference between at least one pair of consecutive harvest times in the 8– h time range. the only significantly enriched biological classification identified for this set was the presence of two genes that encode proteins involved in amino acid transport and metabolism. based on the limited annotation available for these genes, we cannot provide a clear biological interpretation of the role this small set of genes plays in development. it is also possible that the genes are regulated by both dichlorvos exposure and development in such a way that some interaction between the two makes it appear that the expression profile of these genes is not delayed. dichlorvos exposure significantly delayed worm development—by more than h at the h-harvest point as judged by the pa value. overall, the coordinated progression of development does not appear to be disrupted, but we are unable to rule out desynchronization of a small number of developmental processes. biological responses to dichlorvos exposure in order to understand the mechanism of dichlorvos toxicity, we attempted to identify gene expression changes that were associated with dichlorvos exposure but were not the consequence of developmental delay. we initially limited our analysis to the continuous exposure regime and compared data from high concentration samples to control samples. differentially expressed probe sets were selected using a p < 10− and -fold change-filter based on an analysis of covariance that included pa and exposure as factors. a total of probe sets was identified across the five harvest times, with increasing numbers at later time points. the expression levels of these genes increasingly diverged from control over time, but as expected, following washout, their expression rapidly returned to control levels. to determine the biological processes affected by exposure, we compared this list of differentially expressed genes to lists of genes or probe sets from several data sources including gene ontology terms and microarray expression clusters from wormbase, and transcription factor binding sites from the literature. the major processes identified in these analyses are metabolic regulation, innate immunity, stress response, and muscular and neurological regeneration. daf-16-mediated stress responses a large number of the genes differentially expressed in reaction to dichlorvos exposure are transcriptionally regulated by daf- a forkhead transcription factor implicated in the control of metabolism, ageing, innate immunity, and stress responses in c. elegans. the microarray cluster enrichment analysis in particular revealed a large over-representation of genes known to be controlled by daf-. murphy et al. identified genes downstream of daf- in microarray experiments using adult daf-2; daf- hypomorphic mutant worms or worms in which daf- and daf- expression was knocked down with rnai. in our work, we observed changes in expression of these downstream genes that are largely consistent with increasing daf- activity over the course of the exposure.. changes in gene expression are shown for continuous and 8-h washout conditions. b) graph of ins- expression. the means by which daf- is being regulated cannot be directly inferred from the gene expression data, but it appears to be complex. the best studied mechanism of daf- regulation is repression through a phosphorylation cascade initiated by the binding of insulin-like ligands to the daf- receptor. however, we observed that the ins- gene, which encodes a known agonist of daf- is up-regulated in response to dichlorvos exposure which should result in the suppression of daf- activity rather than the inferred activation. this observation suggests that there is a signal activating daf-16-mediated transcription that is independent of daf- and possibly countering the normal repression elicited by ins-7/daf- binding. a likely mechanism for activation of daf- is evident from the elevated adp/atp ratios that have been documented as a result of organophosphate exposure. it is possible that one of the amp-kinases—which are sensitive to adp/atp balance—is transducing this energy metabolism signal to daf- a known ampk target. consistent with this possibility, we found that genes known to be up-regulated in crh- mutants, an ortholog of the camp-inducible creb transcription factor and downstream target of ampk signaling, were over-represented in the microarray cluster analysis and primarily up-regulated in response to dichlorvos. aak- is a strong candidate for the ampk that could be activating daf- as genes downstream of aak- are over-represented in the dichlorvos data set; furthermore, the gene encoding superoxide dismutase, sod- a known daf- target, is down-regulated in aak- mutants. the amp dependent kinase, kin- is a less likely candidate for activating daf- in response to dichlorvos exposure, the message levels for kin- are up-regulated. further, it has been shown to inactivate torc2, an essential transducer in the daf- signaling pathway. although, kin- may play a role in metabolic signaling since it is involved in food sensing, it appears to be chiefly involved in regulating chemoreceptor signaling in sensory neurons, with limited expressions elsewhere in the nervous system and in body wall muscles and hypodermal cells. these factors argue against kin- playing an important role in daf- regulation. in our analysis, we have inferred that daf- is activated based on the behavior of the genes proposed as downstream targets of daf- by murphy et al.. however, approximately 27% of the putative daf- target genes in our study are not regulated in the direction that would be predicted from murphy and coworkers’ experiments. functional analysis did not reveal any specific functions for the discrepant set of genes that were distinct from the functions observed for the complete set of differentially expressed genes. we did, however, observe that this set was enriched with genes containing pqm- binding sites, suggesting an alternative transcriptional regulator event that might be responsive to reactive oxygen stress. in our view, the differences in results in the studies most likely arise from differences in the experimental designs, particularly because it is not possible to distinguish between direct and indirect effects, or among effects due to differential activity of the daf- isoforms. murphy and coworkers used adult worms fed on a bacterial lawn in a purely genetic experiment, while we studied larval worms in axenic medium exposed to a chemical toxicant. it is known that regulatory processes involving daf- are complex, and the differences in life stage of the worm, culture conditions, and experimental perturbation could easily account for the differences we see. we conclude that the effects we have described are largely daf- mediated in response to dichlorvos exposure. based on our observations, we believe that the daf- mediated gene response is a result of reduced energy reserves within the cells and is not a daf- mediated event. two known mechanisms of op toxicity can contribute to this effect: 1) elevated consumption of energy because of increased muscular activity due to continual stimulation when acetylcholine is not degraded and 2) cholinergic receptor-independent mitochondrial disruption. furthermore, in other work, worms exposed to similar dichlorvos concentrations show reduced bacterial feeding. these effects, which are likely additive, appear to be driving forces in daf- mediated gene expression changes affecting metabolic processes. innate immunity related responses we found a number of similarities between responses elicited by dichlorvos and infection in c. elegans. dichlorvos exposure appears to activate some innate immunity signaling cascades, and microarray cluster enrichment analysis revealed that many genes affected by dichlorvos exposure are in data sets derived from worms infected with a diverse set of pathogens and show similar regulation. however, dichlorvos exposure does not fully mimic the responses to infection. microbacterium nematophilum, enterococcus faecalis, photorhabdus luminescens, yersinia pestis, and pseudomonas aeruginosa infections activate many of the same genes as dichlorvos, and show similar regulation. many of these genes are known to be involved in innate immunity and were up-regulated in response to the exposure, including genes encoding lysozymes, invertebrate lysozymes, antibacterial factors, caenacin antimicrobial peptides and c-type lectins. the response appears to overlap but to be distinct from innate immunity because other genes considered to part of the innate immune response are down-regulated, including lys- lys- lys- lys- ilys- and c-type lectins genes. dichlorvos exposure also impacts three of the seven signaling cascades known to be involved in immune defense in c. elegans: the daf- insulin-like receptor pathway, the sma/tgb-β pathway, and the p38/pmk- mapk pathway. the components of the sma/tgb-β pathway—the tgf-β receptor, sma − and sma- —are all up-regulated. the effects of dichlorvos on the p38/pmk- mapk and related pathways are complicated. we observed up-regulation of the toll-domain containing receptor upstream of p38/pmk- along with sek- the map kinase kinase regulator of pmk- and pmk- we also found enrichment for genes with a binding site for the gata-family transcription factor, elt- which has been shown to affect pmk- mediated expression of antimicrobial peptides. the expression of pmk- itself is regulated by npr- a homolog of the neuropeptide y receptor, and we find that of the npr- regulated genes present in our array data are differentially expressed. however, the direction of gene change induced by dichlorvos for pmk- responsive genes and for npr- regulated genes is primarily in the opposite direction from those predicted for pmk- mediated innate immunity. a potential reason for this discrepancy is the increased expression of the gene encoding the phosphatase vhp- which is a known antagonist of pmk-1. vhp- is a dual-specificity phosphatase that interacts with multiple map kinases including kgb- pmk- and pmk-1 and plays a role in modulating the respective signaling pathways during response to stress. it has been shown to be regulated by rnt- during oxidative stress. it is likely that the increased expression of vhp- and the concurrent repression of the pmk- pathway are the results of dichlorvos-induced stress. the role of the daf- insulin-like pathway in the immune response triggered by dichlorvos is somewhat complex. as mentioned above, the daf- pathway is presumably being activated by ins- however, a secondary signal, likely via an ampk, is preventing daf- mediated suppression of daf- activity. what makes this more intriguing is that the induction of ins- is required for p. aeruginosa virulence, presumably due to suppression of the immune response. three other genes encoding insulin-like peptides are also up-regulated and one is down-regulated, although all are affected to a much lesser extent than ins- with the variety of signaling molecules and control mechanisms acting on the daf- pathway, it is difficult to determine whether the observed effects are part of an immune response or represent another role of the daf- signaling pathway. it is clear that the pattern of expression for some of these genes overlaps with the immune response. there has been some discussion in the community as to whether the c. elegans’ immune response is specifically triggered by pathogenic signaling molecules or if it is a generalized stress response, since no homologs for several key components of the mammalian toll-like receptor pathway have been identified in the worm. under the axenic culture conditions utilized in this work, the differentially expressed genes that are involved in the immune response are not being induced by a pathogen. however, it is not clear whether the dichlorvos stimulation of the innate immunity response acts through a “normal” stress sensing process or whether it is affected by an aberrant “toxic” mechanism. there are four primary mechanisms by which dichlorvos could be triggering the innate immune response: generalized stress, altered feeding pathways, autogenic stimulation, or neuronal signaling. the simplest explanation for stimulation of the immune response is through stress pathways induced by the dichlorvos exposure, but others must also be considered. since bacteria are the primary food source for c. elegans under natural conditions, there are overlaps between the immune response to bacteria and adaptation to different bacterial food sources. therefore, the adaptive response triggered by the low energy state within the worms might induce genes necessary for feeding on different bacteria, and these genes may contribute to the observed immune response. another possibility has recently been described in mammals. molecules present in organelles derived from ancient prokaryote symbiotes have characteristics that cause them to trigger immune receptors in the same manner as infectious bacteria. in light of the potential damage to muscle and nerves, the release of these molecules is likely to be occurring and might be the cause of the immune response induced by dichlorvos exposure; though it is still unclear whether c. elegans responds to these molecules in this way. finally, the primary mechanism of dichlorvos toxicity, the inhibition of ache, may disrupt the neuronal signaling that helps to control the innate immune response. it is reasonable to assume that with the variety of genes and processes that comprise the immune response that the stimulation of this response seen in exposed worms includes a combination, if not all, of these processes. neural and musculature repair responses the gene expression changes seen at the later time points highlight the known myopathic and neurodegenerative effects of organophosphate poisoning and are consistent with increased muscular and neuronal growth and repair. there was an increase in expression levels for a large number of muscle-specific genes that are required for myogenesis. genes encoding proteins localized to the m-bands are over represented in a gene ontology analysis; of the genes present in the data set that encode m band proteins are up-regulated. other muscle genes up-regulated include those encoding titin, two paralogs of titin, homologs of troponin t and c, and the homolog of tropomyosin. the evidence for neuronal growth and repair is most apparent through effects that are related to axon guidance and regeneration. two genes in the dlk- mapk pathway, dlk- a mapkkk, and pmk- a p mapk, were up-regulated. the dlk- pathway is essential for axon regeneration. there was also increased expression of multiple genes required for axonal guidance. in addition, two of the pathways that are involved in innate immunity play a dual role in the worm and are also involved in neuronal cell fate determination. the dbl-1/sma tgf-β pathway has been implicated in controlling neuronal cell fate in the male sensory array, and the proteins encoded by two of the genes up-regulated in the pmk- mapk pathway, unc- and sek- function to determine asymmetric cell fate decision during sensory neuron development. it is not possible to conclude from the data in-hand whether or not the changes in these pathways are the result of innate immunity or cell fate determination, or if they may be involved in both. taken together, these gene expression changes suggest that at the later time points in the exposure, the worms have significant cellular damage and have initiated multiple repair processes. the mechanism by which dichlorvos provokes myopathic and neurodegenerative effects has not been defined, but some work suggests that it is related to mitochondrial dysfunction. the changes in metabolic processes and potential mitochondrial dysfunction observed in response to dichlorvos exposure make this a plausible mechanism for the damage induced in the worms. early responses in order to focus on the primary effects of dichlorvos intoxication, we closely examined the genes that were differentially expressed within the first h of the exposure. the most notable biological processes affected at this early time point are immune response and energy metabolism. in the microarray cluster enrichment analysis, this was highlighted by over representation of genes regulated by daf- or differentially expressed during infection by p. aeruginosa. the down-regulation of the genes encoding a fatty acid elongase and an acetoacetyl-coa synthetase suggest a reduction in fatty acid elongation, while the up-regulation of the gene encoding an acyl-coa synthetase implies an increase in β-oxidation., ox_stress and also for time after infection for pa vs op50. pqm- and pha-4_starved are list of genes with binding sites for the respective transcription factors. cells marked with yellow indicate dichlorvos-responsive genes that are present in the gene set indicated at the head of the table. the numbers at the foot indicate the number of genes shared between the h dichlorvos exposure gene set and the indicated experimental data. the direction of change in gene expression is shown in the left-most column, minus sign/green—down-regulated, plus sign/red—up-regulated. over half of the genes differentially expressed at this time have a binding site for pqm-, a transcription factor strongly up-regulated by paraquat exposure, and seven of the genes were identified as being up-regulated by oxidative stress. the primary mechanism of paraquat toxicity is believed to be reactive oxygen species generation mediated by mitochondria. taken together, these observations suggest that the early effects of dichlorvos exposure include a reduction of available energy and increased ros. while increases in β-oxidation contribute to ros generation, the gene response suggests a higher level of ros than would be anticipated by β-oxidation alone. a lower energy state and higher levels of ros are consistent with the idea that dichlorvos causes mitochondrial dysfunction. the reason for overlap with genes differentially regulated during p. aeruginosa infection is not as directly evident. there are no bacteria present, and it is unlikely to be autogenic due to damaged cells. the two most plausible reasons are 1) a stress-mediated immune response, or 2) overlap with non-immune mediated daf- insulin signaling. the latter is of particular relevance considering one of the pathogenic mechanisms of p. aeruginosa is the up-regulation of ins- and while ins- was not differentially expressed at h based on our strict criteria, it is in fact modestly up-regulated and strongly up-regulated at later time points. therefore, the genes regulated by p. aeruginosa infection are likely to include genes regulated by the daf- insulin signaling pathway both in an immune-specific and non-specific manner. the majority of the genes differentially expressed as an early response to dichlorvos is directly related to energy metabolism or is regulated by the daf- insulin-signaling pathway. these effects highlight a change in metabolism within the worms that is indicative of reduction in available energy and the utilization of energy reserves. the inclusion of genes regulated by oxidative stress and pqm- hints at mitochondrial dysfunction as the mechanism for the reduction in energy levels. recovery an important goal of this study was to identify genes and biological processes that might be involved in recovery from dichlorvos intoxication. to identify genes that participate in the recovery process and not ongoing toxicity, data were examined from unexposed and high concentration samples collected after the h washout. we used the same criteria for selecting differentially expressed genes in the washout conditions that we used for the continuous exposure analysis. a total of probe sets meet these criteria for differential expression in at least one of the harvest times. the effect of the exposure diminishes rapidly with time as there are probe sets differentially expressed in the h washout samples and only at the h time point. there were no significant differences at the h time point. as expected, many of these probes are also differentially expressed in the continuous exposure regime. an interesting observation is that probe sets that were not identified as being differentially expressed at the h time point are differentially expressed at later times by the worms in both the continuous exposure and washout protocols, although the magnitude of differential expression is lower for nearly every probe set in the washout samples than the continuously exposed samples. while trending suggests that some of these genes are affected at h, there is a set of genes whose expression is clearly affected only after the washout has occurred. this observation suggests that there are processes affected by dichlorvos which produce delayed effects even after removal of the toxicant. we were unable to find any statistically enriched processes related to these genes, but many appear to be involved in immune response and energy homeostasis. the probe sets differentially expressed in the washout samples, but not the continuous exposure regime, are the most attractive targets for identifying recovery-specific effects. there are seven probe sets that were identified only in the washout samples. however, because the selection criteria were designed to stringently minimize false positives, it is possible that some of them actually could be differentially expressed during the continuous exposure regime. we used less stringent criteria to allow us to discern probe sets that were weakly differentially expressed in the continuous exposure regime and found that only two of the seven probe sets were uniquely differentially expressed after the washout procedure. unfortunately, the proteins encoded by the genes targeted by these probe sets, a widely expressed heat shock protein and an iron-binding protein of unknown function, offer little insight into recovery specific mechanisms. in addition to looking at the genes that were uniquely identified in the recovery period, we examined the entire group of probe sets that were differentially expressed during the recovery phase to determine biological mechanisms potentially involved in recovery. this list was subdivided to include only the probe sets that were differentially expressed using strict criteria after the dichlorvos was removed. in both the complete list and the “after exposure” list, the biological processes that are represented appear to exactly mimic those that are seen during the exposure, albeit to a lesser extent. no recovery-specific mechanisms were identified. after removal of dichlorvos from the medium, the gene expression changes seen in the worms appear to return to baseline levels rather rapidly, and those observed at the first harvest after the washout simply appear to be residual effects of ongoing biological processes resulting from the exposure. by the h harvest time, the study lacks sufficient power to identify any differentially expressed genes. with the exception of expression changes that appear to be the result of delayed development, the treated worms appear to have fully recovered from the exposure. mechanism of toxicity many if not all of the biological effects identified in this work that were caused by dichlorvos exposure can be plausibly tied to inhibition of acetylcholinesterase. developmental timing in c. elegans requires active nicotinic acetylcholine receptors. persistent cholinergic signaling causes a depletion of muscular atp that can mediate an increase in certain metabolic processes. acetylcholine signaling during starvation activates a mapk in pharyngeal muscles. innate immunity in c. elegans is known to be at least partially regulated by neuronal cells and signaling, and this signaling leads to the activation of pmk-1. from these examples, it is clear that neuronal signaling and acetylcholine levels can affect metabolic processes, innate immunity, and developmental time, which are the primary effects observed in response to dichlorvos exposure. though ache inhibition is likely the dominant cause of toxicity in c. elegans, there are two other non-cholinergic mechanisms that have been implicated in dichlorvos toxicity that are consistent with the observations in this work; they are mitochondrial dysfunction and perturbation of calcium homeostasis. the metabolic changes highlighted by alterations in expression levels of daf- target genes underscore the lower energy reserves present in exposed worms. some of this deficiency may be caused by increased muscle stimulation or reduced feeding, but the high levels of induction suggest a more extreme reason for this deficiency, such as a reduction in energy production triggered by the effects of dichlorvos on mitochondrial function. the down-regulation of the genes encoding multiple subunits of the mitochondrial electron transport chain including ones from complex i, complex ii, and complex v suggests a reduction in mitochondrial energy production. further work will be required to delineate what roles mitochondrial dysfunction is playing in dichlorvos toxicity. the role of calcium homeostasis in dichlorvos toxicity is another area where the data do not provide a clear answer but are very suggestive. there are numerous calcium dependent processes that are perturbed by the exposure. for example, the gene encoding unc- a calmodulin-dependent protein kinase that can activate the p38/pmk- mapk pathway, is up-regulated and could be the cause of the activation of that pathway seen in this work. there is also evidence for activation of a calcineurin pathway. the gene encoding cnb- an ortholog of calcineurin b and rcn- a gene whose transcription is dependent on a calcineurin a ortholog and free calcium, are both up-regulated. furthermore, genes known to be regulated by tax- are over represented in the data set. the exact role of calcium homeostasis in dichlorvos toxicity is unknown, and there is no defined mechanism by which dichlorvos effects calcium homeostasis. this is another promising area where more work must be performed to fully understand the effects of dichlorvos toxicity.  upon exposure to dichlorvos, worms show a time-dependent response with progressive delays in development, alterations in metabolism, increased cellular repair, and alterations in innate immunity functions. early in the course of the exposure, the effects appear to be primarily related to energy metabolism. as the exposure progresses, the changes in the metabolic pathways become more pronounced, and there is increased perturbation of immune responses. at the latest times, repair processes have been activated that are presumably a response to damaged nerve and muscles cells. upon removal of the toxicant, gene expression quickly returns to baseline, with the only lingering effect being a synchronized delay in developmental gene expression. this delay in development is the single largest effect on the worms and is evident as early as h after treatment. there are several potential mechanisms by which dichlorvos is mediating toxicity in c. elegans. inhibition of ache is likely a key player as many of the perturbed processes involve neuronal signaling, and the observed behavior changes are consistent with increased accumulation of acetylcholine. the observed changes in the expression of genes involved in energy metabolism suggest that dichlorvos is disrupting mitochondrial function, which has been shown to be a mechanism for neuronal degeneration. calcium homeostasis is potentially another process perturbed by dichlorvos exposure, but it is unclear whether dichlorvos directly or indirectly perturbs calcium homeostasis. through this work, we have identified key genes and pathways responding to dichlorvos toxicity. the primary mechanisms appear to be inhibition of ache and depletion of energy reserves possibly through mitochondrial dysfunction. we have also shown the resilience of the worms by their recovery after removal of the toxicant. this work provides insight into the detailed process of dichlorvos intoxication and sets the groundwork for future efforts. 
because caenorhabditis elegans is relatively resistant to op lethality—particularly through the inhibition of ache—studies in this nematode provide an opportunity to observe alterations in global gene expression following op exposure that cannot be readily observed in less resistant organisms. the worms showed a time-dependent increase in the expression of genes involved in stress responses. we then analyzed gene expression using whole genome microarrays from rna obtained from worms sampled at multiple time points throughout the exposure. however, other mechanisms leading to persistent effects and neurodegeneration remain controversial and difficult to detect. the changes in gene expression observed in the worms following exposure to dichlorvos point towards two potential mechanisms of toxicity: inhibition of ache and mitochondrial disruption. early in the exposure, the predominant effect was on metabolic processes, while at later times, an immune-like response and cellular repair mechanisms dominated the expression pattern. caenorhabditis elegansorganophosphate pesticide intoxicationgene expressiondichlorvosacetylcholinesterase inhibitionmitochondrial disruptiondichlorvos-induced developmental delay following removal of dichlorvos, the gene expression in the worms appeared to relatively rapidly return to steady-state levels.
 the principal toxicity of acute organophosphate pesticide poisoning is the disruption of neurotransmission through inhibition of acetylcholinesterase. however, other mechanisms leading to persistent effects and neurodegeneration remain controversial and difficult to detect. because caenorhabditis elegans is relatively resistant to op lethality—particularly through the inhibition of ache—studies in this nematode provide an opportunity to observe alterations in global gene expression following op exposure that cannot be readily observed in less resistant organisms.  we exposed cultures of worms in axenic, defined medium to dichlorvos under three exposure protocols. in the first, worms were exposed continuously throughout the experiment. in the second and third, the worms were exposed for either or h, the dichlorvos was washed out of the culture, and the worms were allowed to recover. we then analyzed gene expression using whole genome microarrays from rna obtained from worms sampled at multiple time points throughout the exposure. the worms showed a time-dependent increase in the expression of genes involved in stress responses. early in the exposure, the predominant effect was on metabolic processes, while at later times, an immune-like response and cellular repair mechanisms dominated the expression pattern. following removal of dichlorvos, the gene expression in the worms appeared to relatively rapidly return to steady-state levels.  the changes in gene expression observed in the worms following exposure to dichlorvos point towards two potential mechanisms of toxicity: inhibition of ache and mitochondrial disruption. caenorhabditis elegansorganophosphate pesticide intoxicationgene expressiondichlorvosacetylcholinesterase inhibitionmitochondrial disruptiondichlorvos-induced developmental delay 