Transcriptomics and metabolomics analyses provide insights into the difference in toxicity of benzo[a]pyrene and 6-chlorobenzo[a]pyrene to human hepatic cells

Multi-omics integration analysis: Tools and applications in environmental toxicology

Nowadays, traditional single-omics study is not enough to explain the causality between molecular alterations and toxicity endpoints for environmental pollutants. With the development of high-throughput sequencing technology and high-resolution mass spectrometry technology, the integrative analysis of multi-omics has become an efficient strategy to understand holistic biological mechanisms and to uncover the regulation network in specific biological processes. This review summarized sample preparation methods, integration analysis tools and the application of multi-omics integration analyses in environmental toxicology field. Currently, omics methods have been widely applied being as the sensitivity of early biological response, especially for low-dose and long-term exposure to environmental pollutants. Integrative omics can reveal the overall changes of genes, proteins, and/or metabolites in the cells, tissues or organisms, which provide new insights into revealing the overall toxicity effects, screening the toxic targets, and exploring the underlying molecular mechanism of pollutants.

Mechanism of cytochrome P450s mediated interference with glutathione and amino acid metabolisms from halogenated PAHs exposure

Epidemiological evidence indicates that exposure to halogenated polycyclic aromatic hydrocarbons (HPAHs) is associated with many adverse effects. However, the mechanisms of metabolic disorder of HPAHs remains limited. Herein, effects of pyrene (Pyr), and its halogenated derivatives (1-chloropyrene (1-Cl-Pyr), 1-bromopyrene (1-Br-Pyr)) on endogenous metabolic pathways were investigated, in human hepatoma (HepG2) and HepG2-derived cell lines expressing various human cytochrome P450s (CYPs). Non-targeted metabolomics results suggested that 1-Br-Pyr and Pyr exposure (625 nM) induced disruption in glutathione and riboflavin metabolism which associated with redox imbalance, through abnormal accumulation of oxidized glutathione, mediated by bioactivation of CYP2E1. Conversely, CYP2C9-mediated 1-Cl-Pyr significantly interfered with glutathione metabolism intermediates, including glycine, L-glutamic acid and pyroglutamic acid. Notably, CYP1A1-mediated Pyr-induced perturbation of amino acid metabolism which associated with nutrition and glycolipid metabolism, resulting in significant upregulation of most amino acids, whereas halogenated derivatives mediated by CYP1A2 substantially downregulated amino acids. In conclusion, this study suggested that Pyr and its halogenated derivatives exert potent effects on endogenous metabolism disruption under the action of various exogenous metabolic enzymes (CYPs). Thus, new evidence was provided to toxicological mechanisms of HPAHs, and reveals potential health risks of HPAHs in inducing diseases caused by redox and amino acid imbalances.

Insights into the effect of benzotriazoles in liver using integrated metabolomic and transcriptomic analysis

Benzotriazoles (BTRs) are a class of benzoheterocyclic chemicals that are frequently used as metal-corrosive inhibitors, both in industry and daily use. However, the exposure effect information on BTRs remains relatively limited. In this study, an integrated metabolomic and transcriptomic approach was utilized to evaluate the effect of three BTRs, benzotriazole, 6-chloro-1-hydroxi-benzotriazole, and 1-hydroxy-benzotriazole, in the mouse liver with results showing disrupted basal metabolic processes and vitamin and cofactor metabolism after 28 days. The expression of several genes that are related to the inflammatory response and aryl hydrocarbon receptor pathways, such as Gstt2 and Arntl, was altered by the exposure to BTRs. Exposure to BTRs also affected metabolites and genes that are involved in the immune system and xenobiotic responses. The altered expression of several cytochrome P450 family genes reveal a potential detoxification mechanism in the mouse liver. Taken together, our findings provide new insights into the multilayer response of the mouse liver to BTRs exposure as well as a resource for further exploration of the molecular mechanisms by which the response occurs.

PI3K/Akt/FoxO Pathway Mediates Antagonistic Toxicity in HepG2 Cells Coexposed to Deoxynivalenol and Enniatins

The emerging mycotoxins enniatins (ENNs) and the traditional mycotoxin deoxynivalenol (DON) often co-contaminate various grain raw materials and foods. While the liver is their common target organ, the mechanism of their combined effect remains unclear. In this study, the combined cytotoxic effects of four ENNs (ENA, ENA1, ENB, and ENB1) with DON and their mechanisms were investigated using the HepG2 cell line. Additionally, a population exposure risk assessment of these mycotoxins was performed by using in vitro experiments and computer simulations. The results showed that only ENA at 1/4 IC50 and ENB1 at 1/8 IC50 coexposed with DON showed an additive effect, while ENB showed the strongest antagonism at IC50 (CI = 3.890). Co-incubation of ENNs regulated the signaling molecule levels which were disrupted by DON. Transcriptome analysis showed that ENB (IC50) up-regulated the PI3K/Akt/FoxO signaling pathway and inhibited the expression of apoptotic genes (Bax, P53, Caspase 3, etc.) via phosphorylation of FoxO, thereby reducing the cytotoxic effects caused by DON. Both types of mycotoxins posed serious health risks, and the cumulative risk of coexposure was particularly important for emerging mycotoxins.

Halogenated polycyclic aromatic hydrocarbons in Chinese traditional sausages with high salt: Profiles in market samples and formation during home cooking

Halogenated polycyclic aromatic hydrocarbons (XPAHs) are likely to be generated by the reaction between polycyclic aromatic hydrocarbons (PAHs) and halide ions and therefore pose a great threat to high salt food safety. The aim is to explore the profiles of PAHs/XPAHs in market sausages and their formation during home cooking. Concentrations of PAH24 and XPAH18 in 36 market samples were 5.18-408.52 μg/kg and 0.05-0.41 μg/kg, respectively. Smoked sausages exhibited significantly higher concentrations of PAHs than non-smoked sausages. While ready-to-eat sausages presented notably higher XPAH levels than raw sausages. Furthermore, overcooking, such as baking at 220 °C, could result in an exaggerated increase in PAHs. Meanwhile, the increased chlorinated PAH levels after cooking indicated the unintentional formation of XPAHs during sausage cooking. Based on the ILCR model, the intake of 12.7 g/d for males and 10.8 g/d for females is the maximum threshold to achieve negligible risk levels (10-6).

Integration approach of transcriptomics and metabolomics reveals the toxicity of Anthracene and its chlorinated derivatives on human hepatic cells

Polycyclic aromatic hydrocarbons (PAHs) and Chlorinated PAHs (Cl-PAHs) are ubiquitous environmental contaminants. The toxicological information of anthracene (Ant) and its chlorinated derivatives is quite limited. In this study, an integrated metabolomic and transcriptomic analysis approach was adopted to assess the toxic effects triggered by Ant and its chlorinated derivatives, 2-chloroanthracene (2-ClAnt) and 9,10-dichloroanthracen (9,10-Cl2Ant), at human-relevant levels on human normal hepatocyte L02 cells. The cell viability test showed no significant effects on the viability of L02 cells exposed to Ant, 2-ClAnt and 9,10-Cl2Ant at doses of 5-500 nM for 24 h. However, based on transcriptomic analysis, Ant, 2-ClAnt and 9,10-Cl2Ant exposure at human-relevant levels obviously perturbed global gene expression in L02 cells and induced the differential expression of several genes related to cancer development. As the number of genes related to cancer development altered by 9,10-Cl2Ant is the largest, 9,10-Cl2Ant posed greater risks of tumor development than Ant and 2-ClAnt did. Metabolomics analysis demonstrated that Ant, 2-ClAnt and 9,10-Cl2Ant caused significant metabolic perturbation in L02 cells. Pathway enrichment analysis indicated that Ant, 2-ClAnt and 9,10-Cl2Ant mainly perturbed the lipid metabolism and nucleotide metabolism pathway. However, 9,10-Cl2Ant caused a wider perturbation to metabolic pathways than Ant and 2-ClAnt did. In addition, dysregulation of nucleotide metabolism perturbed by Ant, 2-ClAnt and 9,10-Cl2Ant may be associated with the genomic instability and further carcinogenesis.

Photochemical transformation of anthracene (ANT) in surface soil: Chlorination and hydroxylation

To reveal the fate of anthracene (ANT) in soil, the photodegradation behavior of ANT was systematically studied using SiO₂ to simulate a soil environment. Under xenon lamp irradiation, more than 90% of ANT loaded on SiO₂ could be removed after 240 min. Moreover, the effects of water content, chloride ions (Cl^-) and humic acid (HA) were examined. It was found that the presence of water and HA can significantly inhibit the photolysis of ANT on SiO₂, while the addition of chloride alone has no obvious effect. However, when water is present, the inhibition effect of chloride became more obvious. According to radical quenching experiments and electron paramagnetic resonance (EPR) spectra, hydroxyl radicals (•OH) and chlorine radicals (Cl•) were formed in the system. Possible reaction pathways were speculated based on products identified by mass spectrometry. ANT was attacked by •OH to form hydroxylated products, which can be further hydroxylated and oxidized with the final formation of ring-opening products. ANT directly excited by light may also react with Cl• to produce chlorinated polycyclic aromatic hydrocarbons (Cl-PAHs). Finally, the experimental results were verified on real soil. This study provides important information for understanding the photochemical transformation mechanism of ANT at the soil/air interface.

Thiram-induced hyperglycemia causes tibial dyschondroplasia by triggering aberrant ECM remodeling via the gut-pancreas axis in broiler chickens

Pesticide thiram is widely used in agriculture and has been demonstrated to cause tibial dyschondroplasia (TD) in birds. However, the underlying mechanism remains unclear. This work used multi-omics analysis to evaluate the molecular pathways of TD in broilers that were exposed to low level of thiram. Integrative analysis of transcriptomic, proteomic, and metabolomic revealed thiram activity in enhancing pathological ECM remodeling via attenuating the glycolysis pathway and activating the hexosamine and glucuronic acid pathways. Intriguingly, we found hyperglycemia as a crucial factor for ECM overproduction, which resulted in the development of TD. We further demonstrated that high glucose levels are caused by islet secretion dysfunction in thiram-treated broilers. A combination of factors, including lipid disorder, low-grade inflammation, and gut flora disturbance, might contribute to the dysregulation of insulin secretion. The current work revealed the underlying toxicological mechanisms of thiram-induced tibial dyschondroplasia through blood glucose disorder via the gut-pancreas axis in chickens for the first time, which makes it easier to figure out the health risks of pesticides for worldwide policy decisions.

A novel approach to predict the comprehensive EROD potency: Mechanism-based curve fitting of CYP1A1 activity by PAHs

Cytochrome P450 1A1 (CYP1A1) plays critical roles in polycyclic aromatic hydrocarbon (PAH) toxicity, including DNA adduction and ROS generation. Therefore, CYP1A1 activity quantified by the 7-ethoxyresorufin-O-deethylase (EROD) assay (named EROD potency) has been considered a typical biomarker of PAH exposure and toxicity. The EROD dose-response curve always presents a biphasic style, increasing at low concentrations and decreasing at high concentrations of PAHs, but relative effect potency (REP) commonly used in PAH risk assessment is only involved in the increasing phase. In this study, a full bell-shaped EROD curve fitting formula Eq. (1) was obtained by considering both CYP1A1 mRNA induction and enzyme inhibition to completely assess the EROD potency of PAHs. Correspondingly, in silico models of QSAR and docking methods successfully predicted the full EROD curves of PAHs, and the structure-activity relationship indicated that PAHs with heavy molecular weight and large diameter showed stronger EROD potency. Further EROD potency with predicted curve parameters (EC_50,ind and area index) was confirmed by the reported REP (R² = 0.697-0.977) and experimental data from human and mouse cells (R² = 0.700-0.804). This study provides a novel curve fitting for the EROD dose-response relationship and a prediction model for PAH EROD potency.

Parent and Halogenated Polycyclic Aromatic Hydrocarbons in the Serum of Coal-Fired Power Plant Workers: Levels, Sex Differences, Accumulation Trends, and Risks

Workers in coal-fired power plants are at a high risk of exposure to polycyclic aromatic hydrocarbons (PAHs) and their halogenated derivatives (HPAHs), yet no studies have investigated such exposure of HPAHs. In this study, 12 PAHs and 8 chlorinated PAHs, but no brominated PAHs, were detected in >80% of serum samples from workers of a coal-fired power plant in eastern China. Serum HPAH concentrations were higher in plant workers (16-273 ng/g lipid) than in people without occupational exposure (12-51 ng/g lipid), and serum PAH and HPAH concentrations both in male and female workers were positively correlated with the occupational exposure duration, with an estimated doubling time of 11-17 years. Correlations were found between concentrations of ∑₈HPAHs and ∑₁₂PAHs but not between 7-chlorobenz[a]anthracene (7-ClBaA) and 1-chloropyrene (1-ClPyr) and their respective parent PAHs. In males, total concentrations of PAHs and HPAHs were positively correlated with pulmonary hypofunction and hypertension but not with abnormal electrocardiogram. The benzo[a]pyrene equivalents ratio of ∑₈HPAHs/∑₁₂PAHs was 0.3 ± 0.1. Among the HPAHs in the serum, 9-chlorophenanthrene, 7-ClBaA, and 1-ClPyr showed high health risks. This study is the first report on HPAH exposure in coal-fired power plant workers and provides new evidence on the health risks of PAHs and HPAHs in humans.

Insights into the hepatotoxicity of pyrene and 1-chloropyrene using an integrated approach of metabolomics and transcriptomics

The toxicity of pyrene (Pyr) and its chlorinated species have not be comprehensively and clearly elucidated. In this study, an integrated approach of metabolomics and transcriptomics were applied to evaluate the hepatotoxicity of Pyr and 1-chloropyrene (1-Cl-Pyr) at human exposure level, using human L02 hepatocytes. After 24 h exposure to Pyr and 1-Cl-Pyr at 5-500 nM, cell viability was not significantly changed. Transcriptomics results showed that exposure to Pyr and 1-Cl-Pyr at 5 and 50 nM obviously altered the gene expression profiles, but did not significantly induce the expression of genes strongly related to the activation of aryl hydrocarbon receptor (AhR), such as CYP1A1, CYP1B1, AHR, ARNT. Pyr and 1-Cl-Pyr both induced a notable metabolic perturbation to L02 cells. Glycerophospholipid metabolism was found to be the most significantly perturbed pathway after exposure to Pyr and 1-Cl-Pyr, indicating their potential damage to the cell membrane. The other significantly perturbed pathways were identified to be oxidative phosphorylation (OXPHOS), glycolysis, and fatty acid β oxidation, all of which are related to energy production. Exposure to Pyr at 5 and 50 nM induced the up-regulation of fatty acid β oxidation and OXPHOS. The similar result was observed after exposure to 5 nM 1-Cl-Pyr. In contrast, exposure to 50 nM 1-Cl-Pyr induced the down-regulation of OXPHOS by inhibiting the activity of complex I. The obtained results suggested that the modes of action of Pyr and 1-Cl-Pyr on energy production remarkably varied not only with molecular structure change but also with exposure concentration.