Assessment of mixture toxicity of (tri)azoles and their hepatotoxic effects in vitro by means of omics technologies

A pragmatic workflow for human relevance assessment of toxicological pathways and associated new approach methodologies

Currently, safety assessments of chemical substances are predominantly based on animal data. Multiple considerations call for the use of alternative testing strategies that are based on new approach methodologies (NAMs). However, the human relevance of these testing strategies is usually uncertain. This necessitates a harmonized and accepted workflow for assessing their applicability for regulatory purposes. This report proposes such a workflow, applicable for assessing the human relevance of a toxicological pathway and the relevance of NAMs related to the different components of the pathway. The workflow starts with an established toxicological pathway, of which the adverse outcome is relevant for human health risk assessment and that has sufficient weight of evidence. Human relevance is assessed through three main questions, related to the different components (steps) of the pathway, the pathology of human syndromes that have a similar adverse outcome, and quantitative aspects. The latter comprise both interspecies differences and in vitro - in vivo differences. The combined evidence is scored as 'strong', 'moderate' or 'weak' support of human relevance, based on expert judgement. The workflow developed was tested in a case study, through application to an AOP describing craniofacial malformations after in utero exposure to triazoles. Based on evidence collected for two of the three main questions, the case study provided moderate to strong support for human relevance of both the various components of the AOP and its associated NAMs. Furthermore, it demonstrated that the workflow is a promising approach that allows for a more transparent scientific evaluation of human relevance of toxicological pathways and associated NAMs. Therefore, despite some areas for improvement, we consider the workflow an important step forward for application of AOPs and related NAMs in human health risk assessment.

Mechanisms of Azole Potentiation: Insights from Drug Repurposing Approaches

The emergence of azole resistance and tolerance in pathogenic fungi has emerged as a significant public health concern, emphasizing the urgency for innovative strategies to bolster the efficacy of azole-based treatments. Drug repurposing stands as a promising and practical avenue for advancing antifungal therapy, with the potential for swift clinical translation. This review offers a comprehensive overview of azole synergistic agents uncovered through drug repurposing strategies, alongside an in-depth exploration of the mechanisms by which these agents augment azole potency. Drawing from these mechanisms, we delineate strategies aimed at enhancing azole effectiveness, such as inhibiting efflux pumps to elevate azole concentrations within fungal cells, intensifying ergosterol synthesis inhibition, mitigating fungal cell resistance to azoles, and disrupting biological processes extending beyond ergosterol synthesis. This review is beneficial for the development of these potentiators, as it meticulously examines instances and provides nuanced discussions on the mechanisms underlying the progression of azole potentiators through drug repurposing strategies.

Hydramethylnon induces mitochondria-mediated apoptosis in BEAS-2B human bronchial epithelial cells

Typically used household chemicals comprise numerous compounds. Determining mixture toxicity, as observed when using household chemicals containing multiple substances, is of considerable importance from a regulatory perspective. Upon examining the toxic effects of household chemical mixtures, we observed that hydramethylnon combined with tetramethrin resulted in synergistic toxicity. To determine the unknown toxicity mechanism of hydramethylnon, which carries the risk of inhalation exposure when using household chemicals, we conducted a further investigation using BEAS-2B cells, a human bronchial epithelial cell line. Hydramethylnon-induced cytotoxicity was determined following 24 and 48 h of exposure using the water-soluble tetrazolium 1 and lactate dehydrogenase assays. To elucidate the toxicity mechanism, we utilized flow cytometry and measured the levels of apoptosis-related proteins and caspase activities. Given that hydramethylnon, as an insecticide, disrupts the mitochondrial electron transfer chain, we analyzed the relevant mechanisms, including mitochondrial superoxide levels as well as the mitochondrial membrane potential (MMP). Hydramethylnon dose-dependently induced BEAS-2B cell apoptosis via the intrinsic pathway. Furthermore, it significantly increased mitochondrial superoxide levels and disrupted the MMP. Pre-treatment with a caspase inhibitor (Z-DEVD-FMK) confirmed that hydramethylnon induced caspase-dependent apoptosis. Apoptosis, a key event in the toxicological process of chemicals, can lead to lung diseases, including fibrosis and cancer. The results of the present study suggest a mechanism of toxicity of hydramethrylnon, an organofluorine biocide whose toxicity has been little studied, to the lung epithelium. Considering the potential risks associated with inhalation exposure, these results highlight the need for careful management and regulation of hydramethylnon.

Induction of drug metabolizing enzyme and drug transporter expression by antifungal triazole pesticides in human HepaSH hepatocytes

Triazole pesticides are widely used fungicides, to which humans are rather highly exposed. They are known to activate drug-sensing receptors regulating expression of hepatic drug metabolizing enzymes and drug transporters, thus suggesting that the hepatic drug detoxification system is modified by these agrochemicals. To investigate this hypothesis, the effects of 9 triazole fungicides towards expression of drug metabolizing enzymes and transporters were characterized in cultured human HepaSH cells, that are human hepatocytes deriving from chimeric humanized liver TK-NOG mice. Most of triazoles used at 10 μM were found to act as inducers of cytochrome P-450 (CYP) 1A1, CYP1A2, CYP2B6, CYP3A4 and UDP-glucuronosyltransferase 1A1 mRNA levels and of CYP3A4 protein; some triazoles also enhanced mRNA expression of the canalicular transporters P-glycoprotein/MDR1, multidrug resistance-associated protein 2 and breast cancer resistance protein. Triazoles however concomitantly inhibited CYP2B6 and CYP3A4 activities and thus appeared as dual regulators of these CYPs, being both inducers of their expression and inhibitors of their activity. The inducing effect however predominated, at least for bromuconazole, propiconazole and tebuconazole. Bromuconazole was moreover predicted to enhance CYP2B6 and CYP3A4 expression in humans exposed to this fungicide in a chronic, acute or occupational context. These data demonstrate that key-actors of the human hepatic detoxification system are impacted by triazole pesticides, which may have to be considered for the risk assessment of these agrochemicals. They additionally highlight that the use of human HepaSH cells as surrogates to primary human hepatocytes represents an attractive and promising way for studying hepatic effects of environmental chemicals.

Effects of pesticide exposure on the expression of selected genes in normal and cancer samples: Identification of predictive biomarkers for risk assessment

Pesticides pivotal in controlling pests, can represent a threat for human health. Regulatory agencies constantly monitor their harmful effects, regulating their use. Several studies support a positive association between long-term exposure to pesticides and chronic pathologies, such as cancer. Geno-toxicological biomonitoring has proven to be valuable to assess genetic risks associated with exposure to pesticides, representing a promising tool to improve preventive measures and identify workers at higher risk. In this study, a differential gene expression analysis of 70 candidate genes deregulated upon pesticide exposure, was performed in 10 GEO human gene expression DataSets. It was found that six genes (PMAIP1, GCLM, CD36, SQSTM1, ABCC3, NR4A2) had significant AUC predictive values. Also, CD36 was upregulated in non-transformed cell samples and healthy workers, but downregulated in cancer cells. Further validation in larger groups of workers will corroborate the importance of the identified candidates as biomarkers of exposure/effect.

Toxic interactions at the physiological and biochemical levels of green algae under stress of mixtures of three azole fungicides

Assessing the interactions between environmental pollutants and these mixtures is of paramount significance in understanding their negative effects on aquatic ecosystems. However, existing research often lacks comprehensive investigations into the physiological and biochemical mechanisms underlying these interactions. This study aimed to reveal the toxic mechanisms of cyproconazole (CYP), imazalil (IMA), and prochloraz (PRO) and corresponding these mixtures on Auxenochlorella pyrenoidosa by analyzing the interactions at physiological and biochemical levels. Higher concentrations of CYP, IMA, and PRO and these mixtures resulted in a reduction in chlorophyll (Chl) content and increased total protein (TP) suppression, and malondialdehyde (MDA) content exhibited a negative correlation with algal growth. The activity of catalase (CAT) and superoxide dismutase (SOD) decreased with increasing azole fungicides and their mixture concentrations, correlating positively with growth inhibition. Azole fungicides induced dose-dependent apoptosis in A. pyrenoidosa, with higher apoptosis rates indicative of greater pollutant toxicity. The results revealed concentration-dependent toxicity effects, with antagonistic interactions at low concentrations and synergistic effects at high concentrations within the CYP-IMA mixtures. These interactions were closely linked to the interactions observed in Chl-a, carotenoid (Car), CAT, and cellular apoptosis. The antagonistic effects of CYP-PRO mixtures on A. pyrenoidosa growth inhibition can be attributed to the antagonism observed in Chl-a, Chl-b, Car, TP, CAT, SOD, and cellular apoptosis. This study emphasized the importance of gaining a comprehensive understanding of the physiological and biochemical interactions within algal cells, which may help understand the potential mechanism of toxic interaction.

In vitro and in vivo investigation of a thyroid hormone system-specific interaction with triazoles

Alterations in thyroid hormones (TH) and thyroid-stimulating hormone levels are frequently found following exposure to chemicals of concern. Dysregulation of TH levels can severely perturb physiological growth, metabolism, differentiation, homeostasis in the adult and developmental processes in utero. A frequently identified mode of action for this interaction is the induction of hepatic detoxification mechanisms (e.g. SULTs and UGTs), which lead to TH conjugation and elimination and therefore interfere with hormonal homeostasis, fulfilling the endocrine disruptors (EDs) definition. A short-term study in rats with dietary exposure to cyproconazole, epoxiconazole and prochloraz was conducted and hepatocyte hypertrophy, hepatic UGT activity and Phase 1/2 gene expression inductions were observed together with changes in TH levels and thyroid follicular hypertrophy and hyperplasia. To test for specific interaction with the thyroid hormone system, in vitro assays were conducted covering thyroidal I-uptake (NIS), TH transmembranal transport via MCT8 and thyroid peroxidase (TPO) function. Assays for iodothyronine deiodinases (DIO1-DIO3) and iodotyrosine deiodinase (DEHAL1) were included, and from the animal experiment, Dio1 and Dehal1 activities were measured in kidney and liver as relevant local indicators and endpoints. The fungicides did not affect any TH-specific KEs, in vitro and in vivo, thereby suggesting hepatic conjugation as the dominant MoA.

Modeling mixtures interactions in environmental toxicology

In the environment, organisms are exposed to mixtures of different toxicants, which may interact in ways that are difficult to predict when only considering each component individually. Adapting and expanding tools from pharmacology, the toxicology field uses analytical, graphical, and computational methods to identify and quantify interactions in multi-component mixtures. The two general frameworks are concentration addition, where components have similar modes of action and their effects sum together, or independent action, where components have dissimilar modes of action and do not interact. Other interaction behaviors include synergism and antagonism, where the combined effects are more or less than the additive sum of individual effects. This review covers foundational theory, methods, an in-depth survey of original research from the past 20 years, current trends, and future directions. As humans and ecosystems are exposed to increasingly complex mixtures of environmental contaminants, analyzing mixtures interactions will continue to become a more critical aspect of toxicological research.

Combined cytotoxicity of phthalate esters on HepG2 cells: A comprehensive analysis of transcriptomics and metabolomics

Phthalate esters (PAEs), widely used as plasticizers, may pose a potential environmental and human hazard. The aim of this study was to compare the cytotoxicity of di(2-ethylhexyl) phthalates (DEHP) and dibutyl phthalate (DBP)) after their exposure to HepG2 cells alone or in combination. HepG2 cells treated with individual/combined DEHP and DBP at a dose of 10-2 M for 24 h were selected for metabolome and transcriptome analysis. The results demonstrated that exposure to the mixtures of DEHP and DBP caused enhanced or reduced toxic effects regarding 8 pathways with 1065 downregulated genes and 643 upregulated genes, in comparison with those of single chemicals. The combined toxicity of mixture revealed both synergistic and antagonistic interactions between DEHP and DBP. Besides, combined exposure to DEHP and DBP promoted TCA cycle, pyrimidine, and purine metabolism, while an antagonistic effect on fatty acid derangement should require further investigation. To summarize, our results suggest that DEHP exposed alone or combined with DBP caused a variety of metabolic disorders, and the type of combination effects varied among metabolic pathways.

The assessment of the potential hepatotoxicity of new drugs by in vitro metabolomics

Drug hepatotoxicity assessment is a relevant issue both in the course of drug development as well as in the post marketing phase. The use of human relevant in vitro models in combination with powerful analytical methods (metabolomic analysis) is a promising approach to anticipate, as well as to understand and investigate the effects and mechanisms of drug hepatotoxicity in man. The metabolic profile analysis of biological liver models treated with hepatotoxins, as compared to that of those treated with non-hepatotoxic compounds, provides useful information for identifying disturbed cellular metabolic reactions, pathways, and networks. This can later be used to anticipate, as well to assess, the potential hepatotoxicity of new compounds. However, the applicability of the metabolomic analysis to assess the hepatotoxicity of drugs is complex and requires careful and systematic work, precise controls, wise data preprocessing and appropriate biological interpretation to make meaningful interpretations and/or predictions of drug hepatotoxicity. This review provides an updated look at recent in vitro studies which used principally mass spectrometry-based metabolomics to evaluate the hepatotoxicity of drugs. It also analyzes the principal drawbacks that still limit its general applicability in safety assessment screenings. We discuss the analytical workflow, essential factors that need to be considered and suggestions to overcome these drawbacks, as well as recent advancements made in this rapidly growing field of research.

New Approach Methods for Hazard Identification: A Case Study with Azole Fungicides Affecting Molecular Targets Associated with the Adverse Outcome Pathway for Cholestasis

Triazole fungicides such as propiconazole (Pi) or tebuconazole (Te) show hepatotoxicity in vivo, e.g., hypertrophy and vacuolization of liver cells following interaction with nuclear receptors such as PXR (pregnane-X-receptor) and CAR (constitutive androstane receptor). Accordingly, azoles affect gene expression associated with these adverse outcomes in vivo but also in human liver cells in vitro. Additionally, genes indicative of liver cholestasis are affected in vivo and in vitro. We therefore analyzed the capability of Pi and Te to cause cholestasis in an adverse outcome pathway (AOP)-driven approach in hepatic cells of human origin in vitro, considering also previous in vivo studies. Bile salt export pump (BSEP) activity assays confirmed that both azoles are weak inhibitors of BSEP. They alternate the expression of various cholestasis-associated target genes and proteins as well as the mitochondrial membrane function. Published in vivo data, however, demonstrate that neither Pi nor Te cause cholestasis in rodent bioassays. This discrepancy can be explained by the in vivo concentrations of both azoles being well below their EC50 for BSEP inhibition. From a regulatory perspective, this illustrates that toxicogenomics and human in vitro models are valuable tools to detect the potential of a substance to cause a specific type of toxicity. To come to a sound regulatory conclusion on the in vivo relevance of such a finding, results will have to be considered in a broader context also including toxicokinetics in a weight-of-evidence approach.

Triticonazole enantiomers induced enantioselective metabolic phenotypes in Fusarium graminearum and HepG2 cells

The management of Fusarium head blight relies heavily on triazole fungicides. Most of triazole fungicides are chiral, and their enantioselective effects on metabolic phenotypes are poorly understood. Herein, we analyzed the bioactivity of triticonazole against Fusarium graminearum, and ¹H-nuclear magnetic resonance-based metabolomics was used to assess the metabolic disturbances of triticonazole enantiomers in Fusarium graminearum and human hepatocarcinoma cells. Results indicated that the bioactivity of R-triticonazole was 4.28-fold higher than its antipode since it bound stronger with fungal CYP51B and induced more abnormal metabolic processes of Fusarium graminearum, including lipid metabolism, glycolysis, and amino acid metabolism. In human hepatocarcinoma cells, pathways of "alanine, aspartic acid and glutamate metabolism" and "pyruvate metabolism" were disturbed significantly by R-triticonazole; "phenylalanine metabolism" and "taurine-hypotaurine metabolism" were abnormal in the exposure of S-triticonazole. These results suggested that R- and S-triticonazole could affect different metabolic pathways of human hepatocarcinoma cells, and the massively use of inefficient S-triticonazole should be avoided. Our data will help to better understand the enantioselectivity of chiral pesticides and provide a reference for the development of green pesticides.

Use of transcriptomics in hazard identification and next generation risk assessment: A case study with clothianidin

Toxicological risk assessment is essential in the evaluation and authorization of different classes of chemical substances. Genotoxicity and mutagenicity testing are of highest priority and rely on established in vitro systems with bacterial and mammalian cells, sometimes followed by in vivo testing using rodent animal models. Transcriptomic approaches have recently also shown their value to determine transcript signatures specific for genotoxicity. Here, we studied how transcriptomic data, in combination with in vitro tests with human cells, can be used for the identification of genotoxic properties of test compounds. To this end, we used liver samples from a 28-day oral toxicity study in rats with the pesticidal active substances imazalil, thiacloprid, and clothianidin, a neonicotinoid-type insecticide with, amongst others, known hepatotoxic properties. Transcriptomic results were bioinformatically evaluated and pointed towards a genotoxic potential of clothianidin. In vitro Comet and γH2AX assays in human HepaRG hepatoma cells, complemented by in silico analyses of mutagenicity, were conducted as follow-up experiments to check if the genotoxicity alert from the transcriptomic study is in line with results from a battery of guideline genotoxicity studies. Our results illustrate the combined use of toxicogenomics, classic toxicological data and new approach methods in risk assessment. By means of a weight-of-evidence decision, we conclude that clothianidin does most likely not pose genotoxic risks to humans.

•

Analysis of clothianidin genotoxicity in silico, in vitro and in vivo.

•

Application of a toxicogenomics approach to analyze genotoxicity.

•

Weight-of-evidence decision supports classification as “non-genotoxic”.

Bird exposure to fungicides through the consumption of treated seeds: A study of wild red-legged partridges in central Spain

Sown seeds are a key component of many farmland birds' diets due to natural food shortages in autumn and winter. Because these seeds are often treated with pesticides, their ingestion by birds can result in toxic effects. For risk assessment, data on treated seed toxicity should be combined with information about exposure risk for wild birds and the factors that modulate it. We characterized the exposure of red-legged partridges to pesticide-treated seeds through the analysis of digestive contents of birds shot by hunters (n = 194) in an agricultural region in central Spain. We measured the contribution of sown seeds to the partridges' diet and how it related to pesticide exposure. Moreover, we evaluated the influence of landscape composition on the intake of sown seeds and pesticides by partridges. During peak sowing time, seeds constituted half (50.7%) of the fresh biomass ingested by partridges, which consumed mostly winter cereal seeds (42.3% of biomass). Residues of seven fungicides and one insecticide (active ingredients) were detected in 33.0% of birds. The presence of pesticides in digestive contents was linked to the ingestion of cereal sown seeds. Moreover, dietary exposure of birds to pesticides was modulated by landscape characteristics, being lower in areas with heterogeneous landscapes, greater habitat mosaic and more natural vegetation. The estimated dietary intake of pesticides resulting from our field observations, in combination with experimental data on pesticide toxicity, raise concerns about the risks that pesticide-treated cereal seeds pose to granivorous bird populations. Our results highlight the importance of farming landscape composition and diversification, which should be considered as a priority in the agricultural policy to mitigate pesticide risks to farmland birds through the consumption of treated seeds.

A targeted transcriptomics approach for the determination of mixture effects of pesticides

While real-life exposure occurs to complex chemical mixtures, toxicological risk assessment mostly focuses on individual compounds. There is an increasing demand for in vitro tools and strategies for mixture toxicity analysis. Based on a previously established set of hepatotoxicity marker genes, we analyzed mixture effects of non-cytotoxic concentrations of different pesticides in exposure-relevant binary mixtures in human HepaRG hepatocarcinoma cells using targeted transcriptomics. An approach for mixture analysis at the level of a complex endpoint such as a transcript pattern is presented, including mixture design based on relative transcriptomic potencies and similarities. From a mechanistic point of view, goal of the study was to evaluate combinations of chemicals with varying degrees of similarity in order to determine whether differences in mechanisms of action lead to different mixtures effects. Using a model deviation ratio-based approach for assessing mixture effects, it was revealed that most data points are consistent with the assumption of dose addition. A tendency for synergistic effects was only observed at high concentrations of some combinations of the test compounds azoxystrobin, cyproconazole, difenoconazole, propiconazole and thiacloprid, which may not be representative of human real-life exposure. In summary, the findings of our study suggest that, for the pesticide mixtures investigated, risk assessment based on the general assumption of dose addition can be considered sufficiently protective for consumers. The way of data analysis presented in this paper can pave the way for a more comprehensive use of multi-gene expression data in experimental studies related to mixture toxicity.

A prospective whole-mixture approach to assess risk of the food and chemical exposome

Many widely used chemicals result in ubiquitous human exposure from multiple sources, including diet. Legislation mainly deals with the toxicological evaluation of single substances owing to a methodological and conceptual lack of alternatives, and does so within defined silos subject to over 40 distinct regulations in the EU alone. Furthermore, much of the research and many of the initiatives concerned with the assessment and evaluation of chemical mixtures and their potential effects on human health rely on retrospective analysis. Here we propose an approach for the prospective identification, assessment and regulation of mixtures relevant to human health. We address two distinct aspects of toxicology-which chemicals actually do occur together, and how potential mixture-related health hazards can be predicted-with an adapted concept of the exposome and large-scale hazard screens. The proactive use of the likelihood of co-exposure, together with the new approach of methods-based testing, may be a timely and feasible way of identifying those substances and mixtures where hazards may have been overlooked and regulatory action is needed. Ideally, we would generate co-exposure patterns for specific consumer groups, depending on lifestyle and dietary habits, to assess the specific risk of identified mixtures.

Fungicide bromuconazole has the potential to induce hepatotoxicity at the physiological, metabolomic and transcriptomic levels in rats

Bromuconazole (BROMU), a representative triazole fungicide, has been widely used in agriculture for its low cost and highly efficiency against various fungi. BROMU residue was often detected in the environment and food chain, even though there is indication of health risk to animals, and in humans. However, the data related to the toxicity of BROMU in animals remains unclear, and the mechanism is still not fully elucidated. Here, male adult rats were exposed to 0, 13.8, 32.8 and 65.6 mg/kg/d of BROMU for 10 days by oral gavage. It was observed that short time BROMU exposure not only caused liver histological damage, including vacuolar degeneration of hepatocytes with pyknotic nuclei, but also changed the levels of some hepatic physiological parameters, including aspartate transaminase (AST), triglyceride (TG), pyruvate and total cholesterol (TC), indicating that BROMU causes hepatotoxicity in rats. In addition, according to the transcriptomics and metabolomics analysis, a total of 58 metabolites and 259 genes significantly changed in the high-dose BROMU treated group. Although several different pathways are involved, lipid metabolism- and bile acids metabolism-related pathways were highlighted in both metabolomics and transcriptomics analysis. More importantly, further validation had proven that BROMU could not only interact with peroxisome proliferator-activated receptor γ (PPAR-γ), but also significantly decrease its protein and gene expression in the liver, supporting that BROMU decreased the TG synthesis via inhibiting the PPAR-γ pathway. These results clearly showed that BROMU exposure could result in hepatotoxicity at metabolomic and transcriptomic level in rats. These observations could provide some important steps toward understanding the mechanism underlying BROMU-induced mammalian toxicity.

An adverse outcome pathway on the disruption of retinoic acid metabolism leading to developmental craniofacial defects

Adverse outcome pathway (AOP) is a conceptual framework that links a molecular initiating event (MIE) via intermediate key events (KEs) with adverse effects (adverse outcomes, AO) relevant for risk assessment, through defined KE relationships (KERs). The aim of the present work is to describe a linear AOP, supported by experimental data, for skeletal craniofacial defects as the AO. This AO was selected in view of its relative high incidence in humans and the suspected relation to chemical exposure. We focused on inhibition of CYP26, a retinoic acid (RA) metabolizing enzyme, as MIE, based on robust previously published data. Conazoles were selected as representative stressors. Intermediate KEs are RA disbalance, aberrant HOX gene expression, disrupted specification, migration, and differentiation of neural crest cells, and branchial arch dysmorphology. We described the biological basis of the postulated events and conducted weight of evidence (WoE) assessments. The biological plausibility and the overall empirical evidence were assessed as high and moderate, respectively, the latter taking into consideration the moderate evidence for concordance of dose-response and temporal relationships. Finally, the essentiality assessment of the KEs, considered as high, supported the robustness of the presented AOP. This AOP, which appears of relevance to humans, thus contributes to mechanistic underpinning of selected test methods, thereby supporting their application in integrated new approach test methodologies and strategies and application in a regulatory context.

Transcriptomics analysis of hepatotoxicity induced by the pesticides imazalil, thiacloprid and clothianidin alone or in binary mixtures in a 28-day study in female Wistar rats

Co-occurrence of pesticide residues in food commodities raises a potential safety issue as their mixture effects on human health are largely unknown. In a previous study, we reported the toxicological effects (pathology and histopathology) of imazalil (IMZ), thiacloprid (THI), and clothianidin (CTD) alone and in binary mixtures in a 28-day oral gavage study in female Wistar rats. Five dose levels (up to 350 mg/kg body weight/day) ranging from a typical toxicological reference value to a clear effect dose were applied. In the present study, we undertook a transcriptomics analysis of rat livers by means of total RNA sequencing (RNA-Seq). Bioinformatic data analysis involving Ingenuity Pathway Analysis (IPA) was used to gain mechanistic information on hepatotoxicity-related pathways affected after treatment with the pesticides, alone and in mixtures. Our data show that 2986 genes were differentially regulated by CTD while IMZ and THI had effects on 194 and 225 genes, respectively. All three individual compounds shared a common subset of genes whose network is associated with xenobiotic metabolism and nuclear receptor activation. Similar networks were retrieved for the mixtures. Alterations in the expression of individual genes were in line with the assumption of dose addition. Our results bring new insight into the hepatotoxicity mechanisms of IMZ, THI, and CTD and their mixtures.

Risk Characterization and Probabilistic Concentration-Response Modeling of Complex Environmental Mixtures Using New Approach Methodologies (NAMs) Data from Organotypic in Vitro Human Stem Cell Assays

BACKGROUND

Risk assessment of chemical mixtures or complex substances remains a major methodological challenge due to lack of available hazard or exposure data. Therefore, risk assessors usually infer hazard or risk from data on the subset of constituents with available toxicity values.

OBJECTIVES

We evaluated the validity of the widely used traditional mixtures risk assessment paradigms, Independent Action (IA) and Concentration Addition (CA), with new approach methodologies (NAMs) data from human cell-based in vitro assays.

METHODS

A diverse set of 42 chemicals was tested both individually and as mixtures for functional and cytotoxic effects in vitro. A panel of induced pluripotent stem cell (iPSCs)-derived models (hepatocytes, cardiomyocytes, endothelial, and neurons) and one primary cell type (HUVEC) were used. Bayesian concentration-response modeling of individual chemicals or their mixtures was performed for a total of 47 phenotypes to derive point-of-departure (POD) values. Probabilistic IA or CA was conducted to estimate the mixture effects based on the bioactivity profiles from the individual chemicals and compared with mixture bioactivity.

RESULTS

All mixtures showed significant bioactivity, even though some were constructed using individual chemical concentrations considered "low" or "safe." Even though CA is much more accurate as a predictor of mixture effects in comparison with IA, with CA-based POD typically within an order of magnitude of the actual mixture, in some cases, the bioactivity of the mixtures appeared to be much greater than that of their components under either additivity assumption.

DISCUSSION

These results suggest that CA is a preferred first approximation for predicting mixture toxicity when data for all constituents are available. However, because the accuracy of additivity assumptions varies greatly across phenotypes, we posit that mixtures and complex substances need to be directly tested for their hazard potential. NAMs provide a practical solution that rapidly yields highly informative data for mixtures risk assessment. https://doi.org/10.1289/EHP7600.

Cross-species analysis of hepatic cytochrome P450 and transport protein expression

Most drugs and xenobiotics are metabolized in the liver. Amongst others, different cytochrome P450 (CYP) enzymes catalyze the metabolic conversion of foreign compounds, and various transport proteins are engaged in the excretion of metabolites from the hepatocytes. Inter-species and inter-individual differences in the hepatic levels and activities of drug-metabolizing enzymes and transporters result from genetic as well as from environmental factors, and play a decisive role in determining the pharmacokinetic properties of a compound in a given test system. To allow for a meaningful comparison of results from metabolism studies, it is, therefore, of utmost importance to know about the specific metabolic properties of the test systems, especially about the levels of metabolic enzymes such as the CYPs. Using a targeted proteomics approach, we, therefore, compared the hepatic levels of important CYP enzymes and transporters in different experimental systems in vivo and in vitro, namely Wistar rats, C57/Bl6 mice, mice humanized for the two xeno-sensing receptors PXR (pregnane-X-receptor) and CAR (constitutive androstane receptor), mice with human hepatocyte-repopulated livers, human HepaRG hepatocarcinoma cells, primary human hepatocytes, and human liver biopsies. In addition, the effects of xenobiotic inducers of drug metabolism on CYP enzymes and transporters were analyzed in selected systems. This study for the first time presents a comprehensive overview of similarities and differences in important drug metabolism-related proteins among the different experimental models.

Triazoles and aromatase: The impact of copper cocktails

Triazoles are used as antifungal agents, they mostly inhibit two enzymes: 14α-demethylase and aromatase. These enzymes are utilised also in other species and therefore the affection in non-target species in the environment is expected as well. Besides, triazoles are often being applied in a mixture and they can also interact with other substances present. This study clarifies how three selected representative triazoles (tebuconazole, penconazole and cyproconazole) interact with each other (group effect) and in mixtures (cocktail effect) with copper, essential/toxic for all organisms. Within the experiments on electrospray and collision-induced dissociations (both ESI-MS), it has been found that the fragments correspond to typical triazole metabolites. For their formation, the presence of copper ions is crucial. The inhibitory effect of Cu cocktails on aromatase enzymatic activity has been studied. The presence of Cu ions together with triazole(s) significantly increases the inhibitory effect on aromatase activity. The highest inhibitory effect (more than 60%) on aromatase activity is produced by cocktails containing penconazole and Cu ions, namely by penconazole/Cu and penconazole/tebuconazole/Cu. The reactivity of triazoles in groups is not significantly affected by the interactions among them. Additionally, the role of triazoles in copper Fenton reaction regulation has been observed and described. These changes may be attributed to the formation and stabilization of the complexes with the central Cu ion, with usually one, two or three triazolic ligands, depending on the mixture. The study demonstrates that the interaction of triazoles and Cu ions is a complex process; their impact on metabolism seems to be rather extensive and must be evaluated in the context of biochemical reactions.

The reproductive toxicity and potential mechanisms of combined exposure to dibutyl phthalate and diisobutyl phthalate in male zebrafish (Danio rerio)

Dibutyl phthalate (DBP) and diisobutyl phthalate (DiBP) are phthalate compounds frequently detected in the environment. Despite increasing awareness of their toxicity in human and animals, the male reproductive toxicity of their combined exposure remains elusive. The purposes of this study were to investigate whether combined exposure to DBP and DiBP could induce male reproductive toxicity, and to explore the potential toxicological mechanisms. Adult male zebrafish were exposed to DBP (11, 113 and 1133 μg L^-1), DiBP (10, 103 and 1038 μg L^-1) and their mixtures (Mix) (11 + 10, 113 + 103, 1133 + 1038 μg L^-1) for 30 days, and their effects on plasma hormone secretion, testis histology and transcriptomics were examined. Highest concentrations of Mix exposure caused greater imbalance ratio of T/E2 and more severe structural damage to testis than single exposure. These effects were consistent with the testis transcriptome analysis for which 4570 genes were differentially expressed in Mix exposure, while 2795 and 1613 genes were differentially expressed in DBP and DiBP, respectively. KEGG pathway analysis showed that both single and combined exposure of DBP and DiBP could affect cytokine-cytokine receptor interaction. The difference was that combined exposure could also affect steroid hormone synthesis, extracellular matrix receptor interaction, retinol metabolism, and PPAR signaling pathways. These results demonstrated that combined exposure to DBP and DiBP could disrupt spermatogenesis and elicit male reproductive toxicity in zebrafish.

Multi-level analysis of exposure to triazole fungicides through treated seed ingestion in the red-legged partridge

Triazole fungicides are the most widely used products to treat cereal seeds. Granivorous birds, such as red-legged partridges (Alectoris rufa), which consume seeds left on the surface of fields after sowing, have a high risk of exposure. As triazole fungicides can affect sterol synthesis, we tested the hypothesis that treated seed consumption could alter the synthesis of sex hormones and reduce the reproductive capacity of partridges. We exposed adult partridges to seeds treated with four different formulations containing triazoles as active ingredients (flutriafol, prothioconazole, tebuconazole, and a mixture of the latter two) simulating a field exposure during the late autumn sowing season. All treatments produced biochemical changes and an overexpression of genes encoding for enzymes involved in the biosynthesis of sterols and steroid hormones, such as PMVK, ABCA1, MVD, PSCK9, DHCR7 and HSD17B7. Plasma levels of oestradiol were reduced in partridges exposed to tebuconazole. We also monitored reproduction 3 months after exposure (laying date, egg fertilization and hatching rates). We observed a 14-day delay in the laying onset of partridges that had been exposed to flutriafol as compared to controls. These results show that the consumption of seeds treated with triazole fungicides has the potential to affect granivorous bird reproduction. We recommend the evaluation of lagged reproductive effects as part of the protocols of environmental risk assessment of pesticides in wild birds in light of the effects resulting from the exposure to triazole-treated seeds.

Induction and repression effects on CYP and transporter protein abundance by azole mixture uptake in rat liver

Detection of mixture effects is a major challenge in current experimental and regulatory toxicology. Robust markers are needed that are easy to quantify and responsive to chemical stressors in a broad dose range. Several hepatic enzymes and proteins related to drug metabolism like cytochrome-P-450 (CYP) enzymes and transporters have been shown to be responsive to pesticide active substances in a broad dose range and are therefore good candidates to be used as markers for mixture toxicity. Even though they can be well quantified at the mRNA level, quantification on the protein level is challenging because most of these proteins are membrane bound. Here we report the development of mass spectrometry-based assays using triple-x-proteomics (TXP) antibodies in combination with targeted selected ion monitoring (tSIM) to quantify changes of protein levels due to exposure to mixtures of pesticide active substances. Our results indicate that changes on the protein level of CYP1A1, ABCB2, ABCC3 are in line with observations on the mRNA and enzyme activity level and are indicative of mixture effects. Therefore, the tests are promising to reveal effects by chemical mixture effects in toxicological studies in rats.

The Connection of Azole Fungicides with Xeno-Sensing Nuclear Receptors, Drug Metabolism and Hepatotoxicity

Azole fungicides, especially triazole compounds, are widely used in agriculture and as pharmaceuticals. For a considerable number of agricultural azole fungicides, the liver has been identified as the main target organ of toxicity. A number of previous studies points towards an important role of nuclear receptors such as the constitutive androstane receptor (CAR), the pregnane-X-receptor (PXR), or the aryl hydrocarbon receptor (AHR), within the molecular pathways leading to hepatotoxicity of these compounds. Nuclear receptor-mediated hepatic effects may comprise rather adaptive changes such as the induction of drug-metabolizing enzymes, to hepatocellular hypertrophy, histopathologically detectable fatty acid changes, proliferation of hepatocytes, and the promotion of liver tumors. Here, we present a comprehensive review of the current knowledge of the interaction of major agricultural azole-class fungicides with the three nuclear receptors CAR, PXR, and AHR in vivo and in vitro. Nuclear receptor activation profiles of the azoles are presented and related to histopathological findings from classic toxicity studies. Important issues such as species differences and multi-receptor agonism and the consequences for data interpretation and risk assessment are discussed.

RNA-protein correlation of liver toxicity markers in HepaRG cells

The liver is a main target organ for the toxicity of many different compounds. While in general, in vivo testing is still routinely used for assessing the hepatotoxic potential of test chemicals, the use of in vitro models offers advantages with regard to throughput, consumption of resources, and animal welfare aspects. Using the human hepatoma cell line HepaRG, we performed a comparative evaluation of a panel of hepatotoxicity marker mRNAs and proteins after exposure of the cells to 30 different pesticidal active compounds comprising herbizides, fungicides, insecticides, and others. The panel of hepatotoxicity markers included nuclear receptor target genes, key players of fatty acid and bile acid metabolism-related pathways, as well as recently identified biomarkers of drug-induced liver injury. Moreover, marker genes and proteins were identified, for example, S100P, ANXA10, CYP1A1, and CYP7A1. These markers respond with high sensitivity to stimulation with chemically diverse test compounds already at non-cytotoxic concentrations. The potency of the test compounds, determined as an overall parameter of their ability to deregulate marker expression in vitro, was very similar between the mRNA and protein levels. Thus, this study does not only characterize the response of human liver cells to 30 different pesticides but also demonstrates that hepatotoxicity testing in human HepaRG cells yields well comparable results at the mRNA and protein levels. Furthermore, robust hepatotoxicity marker genes and proteins were identified in HepaRG cells.

The State-of-the Art of Environmental Toxicogenomics: Challenges and Perspectives of "Omics" Approaches Directed to Toxicant Mixtures

The last decade witnessed extraordinary advances in “omics” methods, particularly transcriptomics, proteomics and metabolomics, enabling toxicologists to integrate toxicokinetics and toxicodynamics with mechanistic insights on the mode-of-action of noxious chemicals, single or combined. The toxicology of mixtures is, nonetheless, a most challenging enterprise, especially for environmental toxicologists and ecotoxicologists, who invariably deal with chemical mixtures, many of which contain unknowns. Despite costs and demanding computations, the systems toxicology framework, of which “omics” is a major component, endeavors extracting adverse outcome pathways for complex mixtures. Still, the interplay between the multiple components of gene expression and cell metabolism tends to be overlooked. As an example, the proteome allocates DNA methyltransferases whose altered transcription or loss of function by action of chemicals can have a global impact on gene expression in the cell. On the other hand, chemical insult can produce reactive metabolites and radicals that can intercalate or bind to DNA as well as to enzymes and structural proteins, compromising their activity. These examples illustrate the importance of exploring multiple “omes” and the purpose of “omics” and multi-“omics” for building truly predictive models of hazard and risk. Here we will review the state-of-the-art of toxicogenomics highlighting successes, shortcomings and perspectives for next-generation environmental toxicologists.

The application of omics-based human liver platforms for investigating the mechanism of drug-induced hepatotoxicity in vitro

Drug-induced liver injury (DILI) complicates safety assessment for new drugs and poses major threats to both patient health and drug development in the pharmaceutical industry. A number of human liver cell-based in vitro models combined with toxicogenomics methods have been developed as an alternative to animal testing for studying human DILI mechanisms. In this review, we discuss the in vitro human liver systems and their applications in omics-based drug-induced hepatotoxicity studies. We furthermore present bioinformatic approaches that are useful for analyzing toxicogenomic data generated from these models and discuss their current and potential contributions to the understanding of mechanisms of DILI. Human pluripotent stem cells, carrying donor-specific genetic information, hold great potential for advancing the study of individual-specific toxicological responses. When co-cultured with other liver-derived non-parenchymal cells in a microfluidic device, the resulting dynamic platform enables us to study immune-mediated drug hypersensitivity and accelerates personalized drug toxicology studies. A flexible microfluidic platform would also support the assembly of a more advanced organs-on-a-chip device, further bridging gap between in vitro and in vivo conditions. The standard transcriptomic analysis of these cell systems can be complemented with causality-inferring approaches to improve the understanding of DILI mechanisms. These approaches involve statistical techniques capable of elucidating regulatory interactions in parts of these mechanisms. The use of more elaborated human liver models, in harmony with causality-inferring bioinformatic approaches will pave the way for establishing a powerful methodology to systematically assess DILI mechanisms across a wide range of conditions.