An integrative omics approach to unravel toxicity mechanisms of environmental chemicals: effects of a formulated herbicide

Transcriptome analysis of a springtail, Folsomia candida, reveals energy constraint and oxidative stress during petroleum hydrocarbon exposure

Petroleum hydrocarbon (PHC) contamination in soil is ubiquitous and poses harmful consequences to many organisms. The toxicity of PHC-impacted soil is difficult to predict due to variations in mixture composition and the impacts of natural weathering processes. Hence, high-throughput methods to assess PHC-impacted soils is required to expedite land management decisions. Next-generation sequencing is a robust tool that allows researchers to investigate the effects of contaminants on the transcriptome of organisms and identify molecular biomarkers. In this study, the effects of PHCs on conventional endpoints (i.e., survival and reproduction) and gene expression rates of a model springtail species, Folsomia candida were investigated. Age-synchronized F. candida were exposed to ecologically-relevant concentrations of soils spiked with fresh crude oil to calculate the reproductive EC25 and EC50 values using conventional toxicity testing. Soils spiked to these concentrations were then used to evaluate effects on the F. candida transcriptome over a 7-day exposure period. RNA-seq analysis found 98 and 132 differentially expressed genes when compared to the control for the EC25 and EC50 treatment groups, respectively. The majority of up-regulated genes were related to xenobiotic biotransformation reactions and oxidative stress response, while down-regulated genes coded for carbohydrate and peptide metabolic processes. Promotion of the pentose phosphate pathway was also found. Results suggest that the decreased reproduction rates of F. candida exposed to PHCs is due to energy constraints caused by inhibition of carbohydrate metabolic processes and allocation of remaining energy to detoxify xenobiotics. These findings provide insights into the molecular effects in F. candida following exposure to crude oil for seven days and highlight their potential to be used as a high-throughput screening test for PHC-contaminated sites. Adverse molecular effects can be measured as early as 24 h following exposure, whereas conventional toxicity tests may require a minimum of four weeks.

Hotspots of soil pollution: Possible glyphosate and aminomethylphosphonic acid risks on terrestrial ecosystems and human health

The study presents a literature review of glyphosate (GLY) occurrence and its breakdown product, aminomethylphosphonic acid (AMPA), in soils worldwide, but with a specific focus on South America. In addition, an ecological risk approach based on the ecotoxicological endpoints for key soil biota (e.g., collembolans, and earthworms) assessed the impact of GLY and AMPA on these organisms. A generic probabilistic model for human health risk was also calculated for the different world regions. For what reports the risk for edaphic species and the level of pollution under the worst-case scenario, the South American continent was identified as the region of most concern. Nonetheless, other areas may also be in danger, but no risk could be calculated due to the lack of data. Since tropical countries are the top food exporters worldwide, the results obtained in this study must be carefully examined for their implications on a global scale. Some of the factors behind the high levels of these two chemicals in soils are debated (e.g., permissive protection policies, the extensive use of genetically modified crops), and some possible guidelines are presented that include, for example, further environmental characterisation and management of pesticide residues. The present review integrates data that can be used as a base by policymakers and decision-makers to develop and implement environmental policies.

PROTEOMAS: a workflow enabling harmonized proteomic meta-analysis and proteomic signature mapping

Toxicological evaluation of substances in regulation still often relies on animal experiments. Understanding the substances' mode-of-action is crucial to develop alternative test strategies. Omics methods are promising tools to achieve this goal. Until now, most attention was focused on transcriptomics, while proteomics is not yet routinely applied in toxicology despite the large number of datasets available in public repositories. Exploiting the full potential of these datasets is hampered by differences in measurement procedures and follow-up data processing. Here we present the tool PROTEOMAS, which allows meta-analysis of proteomic data from public origin. The workflow was designed for analyzing proteomic studies in a harmonized way and to ensure transparency in the analysis of proteomic data for regulatory purposes. It agrees with the Omics Reporting Framework guidelines of the OECD with the intention to integrate proteomics to other omic methods in regulatory toxicology. The overarching aim is to contribute to the development of AOPs and to understand the mode of action of substances. To demonstrate the robustness and reliability of our workflow we compared our results to those of the original studies. As a case study, we performed a meta-analysis of 25 proteomic datasets to investigate the toxicological effects of nanomaterials at the lung level. PROTEOMAS is an important contribution to the development of alternative test strategies enabling robust meta-analysis of proteomic data. This workflow commits to the FAIR principles (Findable, Accessible, Interoperable and Reusable) of computational protocols.

Development of a multi-omics extraction method for ecotoxicology: investigation of the reproductive cycle of Gammarus fossarum

Omics study exemplified by proteomics, lipidomics or metabolomics, provides the opportunity to get insight of the molecular modifications occurring in living organisms in response to contaminants or in different physiological conditions. However, individual omics discloses only a single layer of information leading to a partial image of the biological complexity. Multiplication of samples preparation and processing can generate analytical variations resulting from several extractions and instrumental runs. To get all the -omics information at the proteins, metabolites and lipids level coming from a unique sample, a specific sample preparation must be optimized. In this study, we streamlined a biphasic extraction procedure based on a MTBE/Methanol mixture to provide the simultaneous extraction of polar (proteins, metabolites) and apolar compounds (lipids) for multi-omics analyses from a unique biological sample by a liquid chromatography (LC)/mass spectrometry (MS)/MS-based targeted approach. We applied the methodology for the study of female amphipod Gammarus fossarum during the reproductive cycle. Multivariate data analyses including Partial Least Squares Discriminant Analysis and multiple factor analysis were applied for the integration of the multi-omics data sets and highlighted molecular signatures, specific to the different stages.

Ecotoxicological relevance of glyphosate and flazasulfuron to soil habitat and retention functions - Single vs combined exposures

Glyphosate (GLY) and flazasulfuron (FLA) are two non-selective herbicides commonly applied together. However, research focused on their single and combined ecotoxicological impacts towards non-target organisms is still inconclusive. Therefore, this study aimed to test their single effects on soil's habitat and retention functions, and to unravel their combined impacts to earthworms and terrestrial plants. For this, ecotoxicological assays were performed with plants (Medicago sativa), oligochaetes (Eisenia fetida) and collembola (Folsomia candida). Soil elutriates were also prepared and tested in macrophytes (Lemna minor) and microalgae (Raphidocelis subcapitata). FLA (82-413 µg kg^-1) reduced earthworms' and collembola's reproduction and severely impaired M. sativa growth, being much more toxic than GLY (up to 30 mg kg^-1). In fact, the latter only affected plant growth (≥ 9 mg kg^-1) and earthworms (≥ 13 mg kg^-1), especially at high concentrations, with no effects on collembola. Moreover, only elutriates from FLA-contaminated soils significantly impacted L. minor and R. sucapitata. The experiments revealed that the co-exposure to GLY and FLA enhanced the toxic effects of contaminated soils not only on plants but also on earthworms'. However, such increase in toxicity was dependent on GLY residual concentrations in soils. Overall, this work underpins that herbicides risk assessment should consider herbicides co-exposures, since the evaluation of single exposures is not representative of current phytosanitary practices and of the potential effects under field conditions, where residues of different compounds may persist in soils.

Disturbance of cellular calcium homeostasis plays a pivotal role in glyphosate-based herbicide-induced oxidative stress

Glyphosate-based herbicides (GBHs) are the most worldwide used pesticides. The wide application of GBHs contaminates the soil and, consequently, water and food resources reaching human consumption. GBHs induce oxidative stress in non-target organisms, leading to a pro-inflammatory and pro-apoptotic cellular status, promoting tissue dysfunction and, thus, metabolic and neurobehavioral changes. This review presents evidence of oxidative damage induced by GBHs and the mechanism of cell damage and health consequences. To summarize, exposure to GBHs may induce disorders in calcium homeostasis related to the activation of ion channels. Also, alterations in pathways related to redox state regulation must have a primordial role in oxidative stress caused by GBHs.

Glyphosate effects on the female reproductive systems: a systematic review

Glyphosate-based herbicides (GBHs) are organophosphate pesticides, which interrupt the chemicals involved in the endocrine system and cause lifelong disorders in women's reproductive system. The current study was designed to systematically evaluate the association between GBH exposure and the female reproductive tract. According to PRISMA Guidelines, the systematic review was performed, searching online databases, including Google Scholar, Web of Science, PubMed, and Scopus, throughout April 2020. Studies with Rodent, lamb, and fish or exposed to GBH to affect the female reproductive system were selected. All studies were in the English language. Two investigators independently assessed the articles. The first author's name, publication date, animal model, age, sample size, gender, dose, duration, and route of exposure and outcomes were extracted from each publication. The present review summarizes 14 publications on uterus alterations and oocytes, histological changes ovary, and assessed mRNA expression, protein expression, serum levels progesterone, and estrogen and intracellular Reaction Oxygen Species (ROS) in rodents, fish, and lamb exposed to GHB exposure. Most of the studies reported histological changes in ovarian and uterus tissue, alterations in serum levels, and increased oxidative stress level following exposure to GBH. Additionally, due to alterations in the reproductive systems (e.g., histomorphological changes, reduction of the mature follicles, higher atretic follicles, and interstitial fibrosis), it seems the GBH-induced female these alterations are both dose- and time-dependent. The present findings support an association between GBH exposure and female reproductive system diseases. However, more studies are needed to identify the mechanisms disrupting the effects of GBH and their underlying mechanisms. Considering the current literature, it is recommended that further investigations be focused on the possible effects of various pesticides on the human reproductive system.

Soil Ecotoxicology Needs Robust Biomarkers: A Meta-Analysis Approach to Test the Robustness of Gene Expression-Based Biomarkers for Measuring Chemical Exposure Effects in Soil Invertebrates

Gene expression-based biomarkers are regularly proposed as rapid, sensitive, and mechanistically informative tools to identify whether soil invertebrates experience adverse effects due to chemical exposure. However, before biomarkers could be deployed within diagnostic studies, systematic evidence of the robustness of such biomarkers to detect effects is needed. In our study, we present an approach for conducting a meta-analysis of the robustness of gene expression-based biomarkers in soil invertebrates. The approach was developed and trialed for two measurements of gene expression commonly proposed as biomarkers in soil ecotoxicology: earthworm metallothionein (MT) gene expression for metals and earthworm heat shock protein 70 (HSP70) gene expression for organic chemicals. We collected 294 unique gene expression data points from the literature and used linear mixed-effect models to assess concentration, exposure duration, and species effects on the quantified response. The meta-analysis showed that the expression of earthworm MT was strongly metal concentration dependent, stable over time and species independent. The metal concentration-dependent response was strongest for cadmium, indicating that this gene is a suitable biomarker for this metal. For copper, no clear concentration-dependent response of MT gene expression in earthworms was found, indicating MT is not a reliable biomarker for this metal. For HSP70, overall marginal up-regulation and lack of a concentration-dependent response indicated that this gene is not suitable as a biomarker for organic pollutant effects in earthworms. The present study demonstrates how meta-analysis can be used to assess the status of biomarkers. We encourage colleagues to apply this open-access approach to other biomarkers, as such quantitative assessment is a prerequisite to ensuring that the suitability and limitations of proposed biomarkers are known and stated. Environ Toxicol Chem 2022;41:2124-2138. © 2022 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Omics-based ecosurveillance for the assessment of ecosystem function, health, and resilience

Current environmental monitoring efforts often focus on known, regulated contaminants ignoring the potential effects of unmeasured compounds and/or environmental factors. These specific, targeted approaches lack broader environmental information and understanding, hindering effective environmental management and policy. Switching to comprehensive, untargeted monitoring of contaminants, organism health, and environmental factors, such as nutrients, temperature, and pH, would provide more effective monitoring with a likely concomitant increase in environmental health. However, even this method would not capture subtle biochemical changes in organisms induced by chronic toxicant exposure. Ecosurveillance is the systematic collection, analysis, and interpretation of ecosystem health-related data that can address this knowledge gap and provide much-needed additional lines of evidence to environmental monitoring programs. Its use would therefore be of great benefit to environmental management and assessment. Unfortunately, the science of ‘ecosurveillance’, especially omics-based ecosurveillance is not well known. Here, we give an overview of this emerging area and show how it has been beneficially applied in a range of systems. We anticipate this review to be a starting point for further efforts to improve environmental monitoring via the integration of comprehensive chemical assessments and molecular biology-based approaches. Bringing multiple levels of omics technology-based assessment together into a systems-wide ecosurveillance approach will bring a greater understanding of the environment, particularly the microbial communities upon which we ultimately rely to remediate perturbed ecosystems.

Roundup causes embryonic development failure and alters metabolic pathways and gut microbiota functionality in non-target species

BACKGROUND

Research around the weedkiller Roundup is among the most contentious of the twenty-first century. Scientists have provided inconclusive evidence that the weedkiller causes cancer and other life-threatening diseases, while industry-paid research reports that the weedkiller has no adverse effect on humans or animals. Much of the controversial evidence on Roundup is rooted in the approach used to determine safe use of chemicals, defined by outdated toxicity tests. We apply a system biology approach to the biomedical and ecological model species Daphnia to quantify the impact of glyphosate and of its commercial formula, Roundup, on fitness, genome-wide transcription and gut microbiota, taking full advantage of clonal reproduction in Daphnia. We then apply machine learning-based statistical analysis to identify and prioritize correlations between genome-wide transcriptional and microbiota changes.

RESULTS

We demonstrate that chronic exposure to ecologically relevant concentrations of glyphosate and Roundup at the approved regulatory threshold for drinking water in the US induce embryonic developmental failure, induce significant DNA damage (genotoxicity), and interfere with signaling. Furthermore, chronic exposure to the weedkiller alters the gut microbiota functionality and composition interfering with carbon and fat metabolism, as well as homeostasis. Using the "Reactome," we identify conserved pathways across the Tree of Life, which are potential targets for Roundup in other species, including liver metabolism, inflammation pathways, and collagen degradation, responsible for the repair of wounds and tissue remodeling.

CONCLUSIONS

Our results show that chronic exposure to concentrations of Roundup and glyphosate at the approved regulatory threshold for drinking water causes embryonic development failure and alteration of key metabolic functions via direct effect on the host molecular processes and indirect effect on the gut microbiota. The ecological model species Daphnia occupies a central position in the food web of aquatic ecosystems, being the preferred food of small vertebrates and invertebrates as well as a grazer of algae and bacteria. The impact of the weedkiller on this keystone species has cascading effects on aquatic food webs, affecting their ability to deliver critical ecosystem services. Video Abstract.

Commercial glyphosate-based herbicides effects on springtails (Collembola) differ from those of their respective active ingredients and vary with soil organic matter content

Glyphosate-based herbicides (GBH) are currently the most widely used agrochemicals for weed control. Environmental risk assessments (ERA) on nontarget organisms mostly consider the active ingredients (AIs) of these herbicides, while much less is known on effects of commercial GBH formulations that are actually applied in the field. Moreover, it is largely unknown to what extent different soil characteristics alter potential side effects of herbicides. We conducted a greenhouse experiment growing a model weed population of Amaranthus retroflexus in arable field soil with either 3.0 or 4.1% soil organic matter (SOM) content and treated these weeds either with GBHs (Roundup LB Plus, Touchdown Quattro, Roundup PowerFlex) or their respective AIs (isopropylammonium, diammonium or potassium salts of glyphosate) at recommended dosages. Control pots were mechanically weeded. Nontarget effects were assessed on the surface activity of the springtail species Sminthurinus niger (pitfall trapping) and litter decomposition in the soil (teabag approach). Both GBHs and AIs increased the surface activity of springtails compared to control pots; springtail activity was higher under GBHs than under corresponding AIs. Stimulation of springtail activity was much higher in soil with higher SOM content than with low SOM content (significant treatment x SOM interaction). Litter decomposition was unaffected by GBHs, AIs or SOM levels. We suggest that ERAs for pesticides should be performed with actually applied herbicides rather than only on AIs and should also consider influences of different soil properties.

An integrated transcriptomic- and proteomic-based approach to evaluate the human skin sensitization potential of glyphosate and its commercial agrochemical formulations

We investigated the skin sensitization hazard of glyphosate, the surfactant polyethylated tallow amine (POEA) and two commercial glyphosate-containing formulations using different omics-technologies based on a human dendritic cell (DC)-like cell line. First, the GARD™skin assay, investigating changes in the expression of 200 transcripts upon cell exposure to xenobiotics, was used for skin sensitization prediction. POEA and the formulations were classified as skin sensitizers while glyphosate alone was classified as a non-sensitizer. Interestingly, the mixture of POEA together with glyphosate displayed a similar sensitizing prediction as POEA alone, indicating that glyphosate likely does not increase the sensitizing capacity when associated with POEA. Moreover, mass spectrometry analysis identified differentially regulated protein groups and predicted molecular pathways based on a proteomic approach in response to cell exposures with glyphosate, POEA and the glyphosate-containing formulations. Based on the protein expression data, predicted pathways were linked to immunologically relevant events and regulated proteins further to cholesterol biosynthesis and homeostasis as well as to autophagy, identifying novel aspects of DC responses after exposure to xenobiotics. In summary, we here present an integrative analysis involving advanced technologies to elucidate the molecular mechanisms behind DC activation in the skin sensitization process triggered by the investigated agrochemical materials. SIGNIFICANCE: The use of glyphosate has increased worldwide, and much effort has been made to improve risk assessments and to further elucidate the mechanisms behind any potential human health hazard of this chemical and its agrochemical formulations. In this context, omics-based techniques can provide a multiparametric approach, including several biomarkers, to expand the mechanistic knowledge of xenobiotics-induced toxicity. Based on this, we performed the integration of GARD™skin and proteomic data to elucidate the skin sensitization hazard of POEA, glyphosate and its two commercial mixtures, and to investigate cellular responses more in detail on protein level. The proteomic data indicate the regulation of immune response-related pathways and proteins associated with cholesterol biosynthesis and homeostasis as well as to autophagy, identifying novel aspects of DC responses after exposure to xenobiotics. Therefore, our data show the applicability of a multiparametric integrated approach for the mechanism-based hazard evaluation of xenobiotics, eventually complementing decision making in the holistic risk assessment of chemicals regarding their allergenic potential in humans.

Oppia nitens C.L. Koch, 1836 (Acari: Oribatida): Current Status of Its Bionomics and Relevance as a Model Invertebrate in Soil Ecotoxicology

The oribatid soil mite Oppia nitens C.L. Koch, 1836, is a model microarthropod in soil ecotoxicity testing. This species has a significant role in supporting soil functions and as a suitable indicator of soil contamination. Despite its significance to the environment and to ecotoxicology, however, very little is known of its biology, ecology, and suborganismal responses to contaminants in the soil. In the present review, we present detailed and critical insights into the biology and ecology of O. nitens in relation to traits that are crucial to its adaptive responses to contaminants in soil. We used a species sensitivity distribution model to rank the species sensitivity to heavy metals (cadmium and zinc) and neonicotinoids (imidacloprid and thiacloprid) compared with other standardized soil invertebrates. Although the International Organization for Standardization and Environment and Climate Change Canada are currently standardizing a protocol for the use of O. nitens in soil toxicity testing, we believe that O. nitens is limited as a model soil invertebrate until the molecular pathways associated with its response to contaminants are better understood. These pathways can only be elucidated with information from the mites' genome or transcriptome, which is currently lacking. Despite this limitation, we propose a possible molecular pathway to metal tolerance and a putative adverse outcome pathway to heavy metal toxicity in O. nitens. Environ Toxicol Chem 2019;38:2593-2613. © 2019 SETAC.

Assessment of fipronil toxicity to the freshwater midge Chironomus riparius: Molecular, biochemical, and organismal responses

Fipronil is a phenylpyrazole insecticide that entered the market to replace organochlorides and organophosphates. Fipronil impairs the regular inhibition of nerve impulses that ultimately result in paralysis and death of insects. Because of its use as a pest control, and due to runoff events, fipronil has been detected in freshwater systems near agricultural areas, and therefore might represent a threat to non-target aquatic organisms. In this study, the toxicity of fipronil to the freshwater midge Chironomus riparius was investigated at biochemical, molecular, and whole organism (e.g. growth, emergence, and behavior) levels. At the individual level, chronic (28 days) exposure to fipronil resulted in reduced larval growth and emergence with a lowest observed effect concentration (LOEC) of 0.081 μg L^-1. Adult weight, which is directly linked to the flying performance and fecundity of midges, was also affected (LOEC = 0.040 μg L^-1). Additionally, behavioral changes such as irregular burrowing behavior of C. riparius larvae (EC₅₀ = 0.084 μg L^-1) and impairment of adult flying performance were observed. At a biochemical level, acute (48 h) exposure to fipronil increased cellular oxygen consumption (as indicated by the increase of electron transport system (ETS) activity) and decreased antioxidant and detoxification defenses (as suggested by the decrease in catalase (CAT) and glutathione S-transferase (GST) activities). Exposure to fipronil also caused alterations in the fatty acid profile of C. riparius, since high levels of stearidonic acid (SDA) were observed. A comparison between exposed and non-exposed larvae also revealed alterations in the expression of globins, cytoskeleton and motor proteins, and proteins involved in protein biosynthesis. These alterations may aid in the interpretation of potential mechanisms of action that lead to the effects observed at the organism level. Present results show that environmentally relevant concentrations of fipronil are toxic to chironomid populations which call for monitoring of phenylpyrazole insecticides and of their ecological effects in freshwaters. Present results also emphasize the importance of complementing ecotoxicological data with molecular approaches such as proteomics, for a better interpretation of the mode of action of insecticides in aquatic invertebrates.

From laboratory to the field: Validating molecular markers of effect in Folsomia candida exposed to a fungicide-based formulation

Under controlled laboratory conditions, toxicity data tend to be less variable than in more realistic in-field studies and responses may thus differ from those in the natural environment, creating uncertainty. The validation of data under environmental conditions is therefore a major asset in environmental risk assessment of chemicals. The present study aimed to validate the mode of action of a commercial fungicide formulation in the soil invertebrate F. candida, under more realistic exposure scenarios (in-field bioassay), by targeting specific molecular biomarkers retrieved from laboratory experiments. Organisms were exposed in soil cores under minimally controlled field conditions for 4 days to a chlorothalonil fungicide dosage causing 75% reduction of reproduction in a previous laboratory experiment (127 mg a.i. kg^-1) and half this concentration (60 mg a.i. kg^-1). After exposure, organisms were retrieved and RNA was extracted from each pool of organisms. According to previous laboratorial omics results with the same formulation, ten genes were selected for gene expression analysis by qRT-PCR, corresponding to key genes of affected biological pathways including glutathione metabolism, oxidation-reduction, body morphogenesis, and reproduction. Six of these genes presented a dose-response trend with higher up- or down-regulation with increasing pesticide concentrations. Highly significant correlations between their expression patterns in laboratory and in-field experiments were observed. This work shows that effects of toxicants can be clearly demonstrated in more realistic conditions using validated biomarkers. Our work outlines a set of genes that can be used to assess the early effects of pesticides in a realistic agricultural scenario.

Glyphosate, but not its metabolite AMPA, alters the honeybee gut microbiota

The honeybee (Apis mellifera) has to cope with multiple environmental stressors, especially pesticides. Among those, the herbicide glyphosate and its main metabolite, the aminomethylphosphonic acid (AMPA), are among the most abundant and ubiquitous contaminant in the environment. Through the foraging and storing of contaminated resources, honeybees are exposed to these xenobiotics. As ingested glyphosate and AMPA are directly in contact with the honeybee gut microbiota, we used quantitative PCR to test whether they could induce significant changes in the relative abundance of the major gut bacterial taxa. Glyphosate induced a strong decrease in Snodgrassella alvi, a partial decrease of a Gilliamella apicola and an increase in Lactobacillus spp. abundances. In vitro, glyphosate reduced the growth of S. alvi and G. apicola but not Lactobacillus kunkeei. Although being no bee killer, we confirmed that glyphosate can have sublethal effects on the honeybee microbiota. To test whether such imbalanced microbiota could favor pathogen development, honeybees were exposed to glyphosate and to spores of the intestinal parasite Nosema ceranae. Glyphosate did not significantly enhance the effect of the parasite infection. Concerning AMPA, while it could reduce the growth of G. apicola in vitro, it did not induce any significant change in the honeybee microbiota, suggesting that glyphosate is the active component modifying the gut communities.

Using time-lapse omics correlations to integrate toxicological pathways of a formulated fungicide in a soil invertebrate

The use of an integrative molecular approach can actively improve the evaluation of environmental health status and impact of chemicals, providing the knowledge to develop sentinel tools that can be integrated in risk assessment studies, since gene and protein expressions represent the first response barriers to anthropogenic stress. This work aimed to determine the mechanisms of toxic action of a widely applied fungicide formulation (chlorothalonil), following a time series approach and using a soil model arthropod, Folsomia candida. To link effects at different levels of biological organization, data were collected on reproduction, gene expression and protein levels, in a time series during exposure to a natural soil. Results showed a mechanistic mode of action for chlorothalonil, affecting pathways of detoxification and excretion, immune response, cellular respiration, protein metabolism and oxidative stress defense, causing irregular cell signaling (JNK and NOD ½ pathways), DNA damage and abnormal cell proliferation, leading to impairment in developmental features such as molting cycle and reproduction. The omics datasets presented highly significant positive correlations between the gene expression levels at a certain time-point and the corresponding protein products 2-3 days later. The integrated omics in this study has provided useful insights into pesticide mechanisms of toxicity, evidencing the relevance of such analyses in toxicological studies, and highlighting the importance of considering a time-series when integrating these datasets.