A multi-omics concentration-response framework uncovers novel understanding of triclosan effects in the chlorophyte Scenedesmus vacuolatus

Modeling HepaRG metabolome responses to pyrrolizidine alkaloid exposure for insight into points of departure and modes of action

The new challenges in toxicology demand novel and innovative in vitro approaches for deriving points of departure (PODs) and determining the mode of action (MOA) of chemicals. Therefore, the aim of this original study was to couple in vitro studies with untargeted metabolomics to model the concentration-response of extra- and intracellular metabolome data on human HepaRG cells treated for 48 h with three pyrrolizidine alkaloids (PAs): heliotrine, retrorsine and lasiocarpine. Modeling revealed that the three PAs induced various monotonic and, importantly, biphasic curves of metabolite content. Based on unannotated metabolites, the endometabolome was more sensitive than the exometabolome in terms of metabolomic effects, and benchmark concentrations (BMCs) confirmed that lasiocarpine was the most hepatotoxic PA. Regarding its MOA, impairment of lipid metabolism was highlighted at a very low BMC (first quartile, 0.003 µM). Moreover, results confirmed that lasiocarpine targets bile acids, as well as amino acid and steroid metabolisms. Analysis of the endometabolome, based on coupling concentration-response and PODs, gave encouraging results for ranking toxins according to their hepatotoxic effects. Therefore, this novel approach is a promising tool for next-generation risk assessment, readily applicable to a broad range of compounds and toxic endpoints.

Meta-metabolomic responses of river biofilms to cobalt exposure and use of dose-response model trends as an indicator of effects

The response of the meta-metabolome is rarely used to characterize the effects of contaminants on a whole community. Here, the meta-metabolomic fingerprints of biofilms were examined after 1, 3 and 7 days of exposure to five concentrations of cobalt (from background concentration to 1 × 10-5 M) in aquatic microcosms. The untargeted metabolomic data were processed using the DRomics tool to build dose-response models and to calculate benchmark-doses. This approach made it possible to use 100% of the chemical signal instead of being limited to the very few annotated metabolites (7%). These benchmark-doses were further aggregated into an empirical cumulative density function. A trend analysis of the untargeted meta-metabolomic feature dose-response curves after 7 days of exposure suggested the presence of a concentration range inducing defense responses between 1.7 × 10-9 and 2.7 × 10-6 M, and of a concentration range inducing damage responses from 2.7 × 10-6 M and above. This distinction was in good agreement with changes in the other biological parameters studied (biomass and chlorophyll content). This study demonstrated that the molecular defense and damage responses can be related to contaminant concentrations and represents a promising approach for environmental risk assessment of metals.

Evolutionary approach for pollution study: The case of ionizing radiation

Estimating the consequences of environmental changes, specifically in a global change context, is essential for conservation issues. In the case of pollutants, the interest in using an evolutionary approach to investigate their consequences has been emphasized since the 2000s, but these studies remain rare compared to the characterization of direct effects on individual features. We focused on the study case of anthropogenic ionizing radiation because, despite its potential strong impact on evolution, the scarcity of evolutionary approaches to study the biological consequences of this stressor is particularly true. In this study, by investigating some particular features of the biological effects of this stressor, and by reviewing existing studies on evolution under ionizing radiation, we suggest that evolutionary approach may help provide an integrative view on the biological consequences of ionizing radiation. We focused on three topics: (i) the mutagenic properties of ionizing radiation and its disruption of evolutionary processes, (ii) exposures at different time scales, leading to an interaction between past and contemporary evolution, and (iii) the special features of contaminated areas called exclusion zones and how evolution could match field and laboratory observed effects. This approach can contribute to answering several key issues in radioecology: to explain species differences in the sensitivity to ionizing radiation, to improve our estimation of the impacts of ionizing radiation on populations, and to help identify the environmental features impacting organisms (e.g., interaction with other pollution, migration of populations, anthropogenic environmental changes). Evolutionary approach would benefit from being integrated to the ecological risk assessment process.

Accumulation of polyethylene microplastics in river biofilms and effect on the uptake, biotransformation and toxicity of the antimicrobial triclosan

The interaction of multiple stressors in freshwater ecosystems may lead to adverse effects on aquatic communities and their ecological functions. Microplastics (MPs) are a class of contaminants of emerging concern that can exert both direct and indirect ecotoxicological effects. A growing number of studies have investigated MPs-attached microbial communities, but the interaction between MPs and substrate-associated biofilm (i.e., on natural river substrates, such as stones and sediments) remains poorly studied. In this work, the combined effects of polyethylene MPs (PE-MPs) with a particle size of 10-45 μm (2 mg/L) and the antimicrobial triclosan (TCS) (20 μg/L) were investigated on river biofilms through a short-term exposure experiment (72 h). To the best of authors' knowledge, this is the first time that the combined effects of MPs and chemical contaminants in substrate-associated river biofilms were assessed. Different response parameters were evaluated, including (i) exposure assessment and ii) contaminants effects at different levels: bacterial community composition, antibiotic resistance, extracellular polymeric substances (EPS), photosynthetic efficiency (Yeff), and leucine aminopeptidase activity (LAPA). Triclosan was accumulated in river biofilms (1189-1513 ng/g dw) alongside its biotransformation product methyl-triclosan (20-29 ng/g dw). Also, PE-MPs were detected on biofilms (168-292 MP/cm2), but they had no significant influence on the bioaccumulation and biotransformation of TCS. A moderate shift in bacterial community composition was driven by TCS, regardless of PE-MPs co-exposure (e.g., increased relative abundance of Sphingomonadaceae family). Additionally, Yeff and EPS content were significantly disrupted in TCS-exposed biofilms. Therefore, the most remarkable effects on river biofilms were related to the antimicrobial TCS, whereas single PE-MPs exposure did not alter any of the evaluated parameters. These results demonstrate that biofilms might act as environmental sink of MPs. Although no interaction between PE-MPs and TCS was observed, the possible indirect impact of other MPs-adsorbed contaminants on biofilms should be further assessed.

Editorial trend: adverse outcome pathway (AOP) and computational strategy - towards new perspectives in ecotoxicology

The adverse outcome pathway (AOP) has been conceptualized in 2010 as an analytical construct to describe a sequential chain of causal links between key events, from a molecular initiating event leading to an adverse outcome (AO), considering several levels of biological organization. An AOP aims to identify and organize available knowledge about toxic effects of chemicals and drugs, either in ecotoxicology or toxicology, and it can be helpful in both basic and applied research and serve as a decision-making tool in support of regulatory risk assessment. The AOP concept has evolved since its introduction, and recent research in toxicology, based on integrative systems biology and artificial intelligence, gave it a new dimension. This innovative in silico strategy can help to decipher mechanisms of action and AOP and offers new perspectives in AOP development. However, to date, this strategy has not yet been applied to ecotoxicology. In this context, the main objective of this short article is to discuss the relevance and feasibility of transferring this strategy to ecotoxicology. One of the challenges to be discussed is the level of organisation that is relevant to address for the AO (population/community). This strategy also offers many advantages that could be fruitful in ecotoxicology and overcome the lack of time, such as the rapid identification of data available at a time t, or the identification of "data gaps". Finally, this article proposes a step forward with suggested priority topics in ecotoxicology that could benefit from this strategy.

Multi-omics reveals the increased biofilm formation of Salmonella Typhimurium M3 by the induction of tetracycline at sub-inhibitory concentrations

Exposure to sub-inhibitory concentrations (sub-MICs) of antibiotics could induce the biofilm formation of microorganisms, but its underlying mechanisms still remain elusive. In the present work, biofilm formation by Salmonella Typhimurium M3 was increased when in the presence of tetracycline at sub-MIC, and the highest induction was observed with tetracycline at 1/8 MIC. The integration of RNA-sequencing and untargeted metabolomics was applied in order to further decipher the potential mechanisms for this observation. In total, 439 genes and 144 metabolites of S. Typhimurium M3 were significantly expressed after its exposure to 1/8 MIC of tetracycline. In addition, the co-expression analysis revealed that 6 genes and 8 metabolites play a key role in response to 1/8 MIC of tetracycline. The differential genes and metabolites were represented in 12 KEGG pathways, including five pathways of amino acid metabolism (beta-alanine metabolism, tryptophan metabolism, arginine and proline metabolism, phenylalanine, tyrosine and tryptophan biosynthesis, and glutathione metabolism), three lipid metabolism pathways (biosynthesis of unsaturated fatty acids, fatty acid degradation, and fatty acid biosynthesis), two nucleotide metabolism pathways (purine metabolism, and pyrimidine metabolism), pantothenate and CoA biosynthesis, and ABC transporters. Metabolites (anthranilate, indole, and putrescine) from amino acid metabolism may act as signaling molecules to promote the biofilm formation of S. Typhimurium M3. The results of this work highlight the importance of low antimicrobial concentrations on foodborne pathogens of environmental origin.

Triclosan toxicity in a model cyanobacterium (Anabaena flos-aquae): Growth, photosynthesis and transcriptomic response

Exposure to triclosan (TCS) has been reported to reduce photosynthetic pigments, suppress photosynthesis, and inhibit growth in both prokaryotic and eukaryotic algae including Anabaena flos-aquae (a model cyanobacterium). In particular, cyanobacteria are more sensitive to TCS toxicity compared to eukaryotic algae possibly due to the structural similarity to bacteria (target organisms); however, whether TCS exerts its toxicity to cyanobacteria by targeting signaling pathways of fatty acid biosynthesis as in bacteria remains virtually unknown, particularly at environmental exposure levels. With the complete genome sequence of A. flos-aquae presented in this study, the transcriptomic alterations and potential toxic mechanisms in A. flos-aquae under TCS stress were revealed. The growth, pigments and photosynthetic activity of A. flos-aquae were markedly suppressed following a 7-day TCS exposure at 0.5 µg/L but not 0.1 µg/L (both concentrations applied are environmentally relevant). The transcriptomic sequencing analysis showed that signaling pathways, such as biofilm formation - Pseudomonas aeruginosa, two-component system, starch and sucrose metabolism, and photosynthesis were closely related to the TCS-induced growth inhibition in the 0.5 µg/L TCS treatment. Photosynthesis systems and potentially two-component system were identified to be sensitive targets of TCS toxicity in A. flos-aquae. The present study provides novel insights on TCS toxicity at the transcriptomic level in A. flos-aquae.

Effects of Chemical Compounds on the Activity of the N-acetyl-β-D-Glucosaminidase of the Marine Prawn, Palaemon serratus: Screening In Vitro

N-acetyl-β-D-glucosaminidase (NAGase) is important for crustaceans because the enzyme activity is necessary for the molting process. The present study aimed to assess the sensitivity of Palaemon serratus NAGase activity to a set of compounds of diverse chemical families in the context of in vitro exposures. Compounds representing different chemical families were selected according to their abundance, impact in the environment, and relevance as disruptors of the molting process. In a first step, four solvents (dimethylsulfoxide [DMSO], methanol, acetone, and ethanol) were tested to determine their suitability to dissolve hydrophobic compounds without affecting NAGase activity. Exclusively, ethanol had no effect on enzyme activity and on the integrity of the proteins present in the enzyme extract. The 18 other compounds were tested and four of these compounds, pentoxifylline, fenoxycarb, dithiocarbamate, and RH5849, showed a specific alteration on the activity of NAGase, without affecting the protein content. However, cadmium, zinc, and glyphosate showed a nonspecific alteration, affecting both the enzyme activity and the proteins, whereas ibuprofen exclusively altered the protein content. Finally, 10 of the 22 tested compounds (including DMSO, acetone, and methanol) showed a direct alteration of NAGase activity. Environ Toxicol Chem 2023;42:846-858. © 2023 SETAC.

Multiomics Point of Departure (moPOD) Modeling Supports an Adverse Outcome Pathway Network for Ionizing Radiation

While adverse biological effects of acute high-dose ionizing radiation have been extensively investigated, knowledge on chronic low-dose effects is scarce. The aims of the present study were to identify hazards of low-dose ionizing radiation to Daphnia magna using multiomics dose-response modeling and to demonstrate the use of omics data to support an adverse outcome pathway (AOP) network development for ionizing radiation. Neonatal D. magna were exposed to γ radiation for 8 days. Transcriptomic analysis was performed after 4 and 8 days of exposure, whereas metabolomics and confirmative bioassays to support the omics analyses were conducted after 8 days of exposure. Benchmark doses (BMDs, 10% benchmark response) as points of departure (PODs) were estimated for both dose-responsive genes/metabolites and the enriched KEGG pathways. Relevant pathways derived using the BMD modeling and additional functional end points measured by the bioassays were overlaid with a previously published AOP network. The results showed that several molecular pathways were highly relevant to the known modes of action of γ radiation, including oxidative stress, DNA damage, mitochondrial dysfunction, protein degradation, and apoptosis. The functional assays showed increased oxidative stress and decreased mitochondrial membrane potential and ATP pool. Ranking of PODs at the pathway and functional levels showed that oxidative damage related functions had relatively low PODs, followed by DNA damage, energy metabolism, and apoptosis. These were supportive of causal events in the proposed AOP network. This approach yielded promising results and can potentially provide additional empirical evidence to support further AOP development for ionizing radiation.

The effect of disinfectants and antiseptics on co- and cross-selection of resistance to antibiotics in aquatic environments and wastewater treatment plants

The outbreak of the SARS-CoV-2 pandemic led to increased use of disinfectants and antiseptics (DAs), resulting in higher concentrations of these compounds in wastewaters, wastewater treatment plant (WWTP) effluents and receiving water bodies. Their constant presence in water bodies may lead to development and acquisition of resistance against the DAs. In addition, they may also promote antibiotic resistance (AR) due to cross- and co-selection of AR among bacteria that are exposed to the DAs, which is a highly important issue with regards to human and environmental health. This review addresses this issue and provides an overview of DAs structure together with their modes of action against microorganisms. Relevant examples of the most effective treatment techniques to increase the DAs removal efficiency from wastewater are discussed. Moreover, insight on the resistance mechanisms to DAs and the mechanism of DAs enhancement of cross- and co-selection of ARs are presented. Furthermore, this review discusses the impact of DAs on resistance against antibiotics, the occurrence of DAs in aquatic systems, and DA removal mechanisms in WWTPs, which in principle serve as the final barrier before releasing these compounds into the receiving environment. By recognition of important research gaps, research needs to determine the impact of the majority of DAs in WWTPs and the consequences of their presence and spread of antibiotic resistance were identified.

High-throughput screening paradigms in ecotoxicity testing: Emerging prospects and ongoing challenges

The rapidly increasing number of new production chemicals coupled with stringent implementation of global chemical management programs necessities a paradigm shift towards boarder uses of low-cost and high-throughput ecotoxicity testing strategies as well as deeper understanding of cellular and sub-cellular mechanisms of ecotoxicity that can be used in effective risk assessment. The latter will require automated acquisition of biological data, new capabilities for big data analysis as well as computational simulations capable of translating new data into in vivo relevance. However, very few efforts have been so far devoted into the development of automated bioanalytical systems in ecotoxicology. This is in stark contrast to standardized and high-throughput chemical screening and prioritization routines found in modern drug discovery pipelines. As a result, the high-throughput and high-content data acquisition in ecotoxicology is still in its infancy with limited examples focused on cell-free and cell-based assays. In this work we outline recent developments and emerging prospects of high-throughput bioanalytical approaches in ecotoxicology that reach beyond in vitro biotests. We discuss future importance of automated quantitative data acquisition for cell-free, cell-based as well as developments in phytotoxicity and in vivo biotests utilizing small aquatic model organisms. We also discuss recent innovations such as organs-on-a-chip technologies and existing challenges for emerging high-throughput ecotoxicity testing strategies. Lastly, we provide seminal examples of the small number of successful high-throughput implementations that have been employed in prioritization of chemicals and accelerated environmental risk assessment.

Microalgae-based wastewater treatment - Microalgae-bacteria consortia, multi-omics approaches and algal stress response

Sustainable environmental management is one of the important aspects of sustainable development goals. Increasing amounts of wastewaters (WW) from exponential economic growth is a major challenge, and conventional treatment methods entail a huge carbon footprint in terms of energy use and GHG emissions. Microalgae-based WW treatment is a potential candidate for sustainable WW treatment. The nutrients which are otherwise unutilized in the conventional processes are recovered in the beneficial microalgal biomass. This review presents comprehensive information regarding the potential of microalgae as sustainable bioremediation agents. Microalgae-bacterial consortia play a critical role in synergistic nutrient removal, supported by the complex nutritional and metabolite exchange between microalgae and the associated bacteria. Design of effective microalgae-bacteria consortia either by screening or by recent technologies such as synthetic biology approaches are highly required for efficient WW treatment. Furthermore, this review discusses the crucial research gap in microalgal WW treatment - the application of a multi-omics platform for understanding microalgal response towards WW conditions and the design of effective microalgal or microalgae-bacteria consortia based on genetic information. While metagenomics helps in the identification and monitoring of the microbial community throughout the treatment process, transcriptomics, proteomics and metabolomics aid in studying the algal cellular response towards the nutrients and pollutants in WW. It has been established that the integration of microalgal processes into conventional WW treatment systems is feasible. In this direction, future research directions for microalgal WW treatment emphasize the need for identifying the niche in WW treatment, while highlighting the pilot sale plants in existence. Microalgae-based WW treatment could be a potential phase in the waste hierarchy of circular economy and sustainability, considering WWs are a rich secondary source of finite resources such as nitrogen and phosphorus.

Daphnia as a Sentinel Species for Environmental Health Protection: A Perspective on Biomonitoring and Bioremediation of Chemical Pollution

Despite available technology and the knowledge that chemical pollution damages human and ecosystem health, chemical pollution remains rampant, ineffectively monitored, rarely prevented, and only occasionally mitigated. We present a framework that helps address current major challenges in the monitoring and assessment of chemical pollution by broadening the use of the sentinel species Daphnia as a diagnostic agent of water pollution. And where prevention has failed, we propose the application of Daphnia as a bioremediation agent to help reduce hazards from chemical mixtures in the environment. By applying "omics" technologies to Daphnia exposed to real-world ambient chemical mixtures, we show improvements at detecting bioactive components of chemical mixtures, determining the potential effects of untested chemicals within mixtures, and identifying targets of toxicity. We also show that using Daphnia strains that naturally adapted to chemical pollution as removal agents of ambient chemical mixtures can sustainably improve environmental health protection. Expanding the use of Daphnia beyond its current applications in regulatory toxicology has the potential to improve both the assessment and the remediation of environmental pollution.

Community metabolomics provides insights into mechanisms of pollution-induced community tolerance of periphyton

Chemical pollution is a major concern for freshwater ecosystems, but the impact and mechanisms of chemical stressors on communities are barely understood. Pollution stress beyond natural homeostatic capacities can trigger succession of tolerant species within a community, enhancing the overall community tolerance. This process was operationalized in the Pollution-Induced Community Tolerance (PICT) concept and applied in many case studies, however, the molecular mechanisms of community tolerance and implications for ecological functions remain largely unexplored. Our study aimed to demonstrate that 1) community metabolomics can unravel potential mechanisms of PICT in periphyton and 2) induced tolerance helps to maintain primary production under re-occuring pollution. To this end, we grew periphyton for 5 weeks with and without the model herbicide diuron in microcosms, quantified PICT, and determined the related metabolic fingerprint of periphyton by GC-MS-based untargeted metabolomics. Further, we explored the autotrophic community based on pigment composition and functional parameters including photosynthesis and gross primary production. Chronic diuron exposure resulted in a shift in pigment composition, higher community tolerance and an individual metabolic fingerprint in the contaminated communities. Opposing responses of selected metabolites during a short-term exposure indicated differences in diuron pre-adaptation in the different communities. Metabolites (threonic acid and two sugar acid lactones) were found to be related to tolerance development, suggesting that ascorbate metabolism was induced in contaminated communities. Despite these compensating mechanism, contaminated communities were compromised in production-to-respiration ratio and biomass. A ranking of sensitivity thresholds of different biological endpoints revealed that metabolites were less sensitive than photosynthetic parameters, which reflects the mode-of-action of the herbicide. In conclusion, we could demonstrate that community metabolomics is able to unravel complex biochemical changes and allows mechanistic insights into community tolerance. Moreover, we were able to show that induced community tolerance was insufficient to safeguard functions like primary production.

Metabolomics insight into the influence of environmental factors in responses of freshwater biofilms to the model herbicide diuron

Freshwater biofilms have been increasingly used during the last decade in ecotoxicology due to their ecological relevance to assess the effect(s) of environmental stress at the community level. Despite growing knowledge about the effect of various stressors on the structure and the function of these microbial communities, a strong research effort is still required to better understand their response to chemical stress and the influence of environmental stressors in this response. To tackle this challenge, untargeted metabolomics is an approach of choice because of its capacity to give an integrative picture of the exposure to multiple stress and associated effect as well as identifying the molecular pathways involved in these responses. In this context, the present study aimed to explore the use of an untargeted metabolomics approach to unravel at the molecular/biochemical level the response of the whole biofilm to chemical stress and the influence of various environmental factors in this response. To this end, archived high-resolution mass spectrometry data from previous experiments at our laboratory on the effect of the model photosynthesis inhibitor diuron on freshwater biofilm were investigated by using innovative solutions for OMICs data (e.g., DRomics) and more usual chemometric approaches (multivariate and univariate statistical analyses). The results showed a faster (1 min) and more sensitive response of the metabolome to diuron than usual functional descriptors, including photosynthesis. Also, the metabolomics response to diuron resulted from metabolites following various trends (increasing, decreasing, U/bell shape) along increasing concentration and time. This metabolomics response was influenced by the temperature, photoperiod, and flow. A focus on a plant-specific omega-3 (eicosapentaenoic acid) playing a key role in the trophic chain highlighted the potential relevance of metabolomics approach to establish the link between molecular alteration and ecosystem structure/functioning impairment but also how complex is the response and the influence of all the tested factors on this response at the metabolomics level. Altogether, our results underline that more fundamental researches are needed to decipher the metabolomics response of freshwater biofilm to chemical stress and its link with physiological, structural, and functional responses toward the unraveling of adverse outcome pathways (AOP) for key ecosystem functions (e.g., primary production).

Study on toxicity effects of environmental pollutants based on metabolomics: A review

In the past few decades, the toxic effects of environmental pollutants on non-target organisms have received more and more attention. As a new omics technology, metabolomics can clarify the metabolic homeostasis of the organism at the overall level by studying the changes in the relative contents of endogenous metabolites in the organism. Recently, a large number of studies have used metabolomics technology to study the toxic effects of environmental pollutants on organisms. In this review, we reviewed the analysis processes and data processes of metabolomics and its application in the study of the toxic effects of environmental pollutants including heavy metals, pesticides, polychlorinated biphenyls, polycyclic aromatic hydrocarbons, polybrominated diphenyl ethers and microplastics. In addition, we emphasized that the combination of metabolomics and other omics technologies will help to explore the toxic mechanism of environmental pollutants and provide new research ideas for the toxicological evaluation of environmental pollutants.

Prediction of adverse effects of effluents containing phenolic compounds in the Ba River on the ovary of fish (Hemiculter leucisculus) using transcriptomic and metabolomic analyses

The aim of this work was to evaluate the endocrine disrupting effects on the ovarian development of sharpbelly (Hemiculter leucisculus) caused by effluents containing phenolic compounds. This was achieved using integrated transcriptomic and metabolomic analyses, along with histopathological examinations. Sharpbelly, an indigenous freshwater fish widely distributed in East Asia, were collected by pole fishing from three sampling sites in the Ba River. These sampling sites include a mid-stream site near a wastewater outfall and a reference site located upstream and a far field comparison site located downstream. In sharpbelly collected near the wastewater discharge, the oocyte development was activated, compared to the other two sites. Histopathological alterations in the fish ovaries were likely due to the upregulated steroid hormone biosynthesis process, as suggested by the differentially expressed genes (e.g., hsd3b, hsd17b1) and differentially accumulated metabolites (e.g., pregnenolone). Additionally, under the stress of effluents containing phenolic compounds, genes related to the signaling pathways for oxidative phosphorylation and leukocyte transendothelial migration were dysregulated, suggesting the potential induction of inflammation and several ovarian diseases. Overall, these findings suggest that effluents containing phenolic compounds influence ovary development and reproductive function of female sharpbelly. Whether there is any resulting dysfunction of folliculogenesis, abnormality of ovulation, production of premature eggs and/or potential induction of ovarian cancers remains to be determined by further studies, for a better evaluation on effluents containing phenolic compounds to the fish fertility and the health of their offspring, and even the stability of the wild fish population. Notably, the integration of transcriptomics and metabolomics can complement the routine chemical analysis to comprehensively monitor the effects of wastewater treatment plant effluents on the health of wild fish.

Untargeted and targeted analysis of sarin poisoning biomarkers in rat urine by liquid chromatography and tandem mass spectrometry

Chemical warfare agents continue to pose a real threat to humanity, despite their prohibition under the Chemical Weapons Convention. Sarin is one of the most toxic and lethal representatives of nerve agents. The methodology for the targeted analysis of known sarin metabolites has reached great heights, but little attention has been paid to the untargeted analysis of biological samples of victims exposed to this deadly poisonous substance. At present, the development of computational and statistical methods of analysis offers great opportunities for finding new metabolites or understanding the mechanisms of action or effect of toxic substances on the organism. This study presents the targeted LC-MS/MS determination of methylphosphonic acid and isopropyl methylphosphonic acid in the urine of rats exposed to a non-lethal dose of sarin, as well as the untarget urine analysis by LC-HRMS. Targeted analysis of polar acidic sarin metabolites was performed on a mixed-mode reversed-phase anion-exchange column, and untargeted analysis on a conventional reversed-phase C18 column. Isopropyl methylphosphonic acid was detected and quantified within 5 days after subcutaneous injection of sarin at a dose of 1/4 LD₅₀. A combination of generalized additive mixed models and dose-response analysis with database searches using accurate mass of precursor ions and corresponding MS/MS spectra enabled us to propose new six potential biomarkers of biological response to exposure. The results confirm the well-known fact that sarin poisoning has a significant impact on the victims' metabolome, with inhibition of acetylcholinesterase being just the first step and trigger of the complex toxicodynamic response.