Benchmarking Water Quality from Wastewater to Drinking Waters Using Reduced Transcriptome of Human Cells

Recent advances in bioassays for assessing the toxicity of environmental contaminants in effect-directed analysis

Chemical cocktails in the environment can cause adverse impacts on ecosystems and human health even at low concentrations. Effect-directed analysis (EDA) has proven to be very valuable in identifying key toxic substances in environmental mixtures. For this, it is important to carefully select accurate bioassays from a wide range of tests for EDA when applying it to actual environmental samples. This article reviews studies published from 2014 to 2023 that have applied EDA and summarizes the bioassays and their corresponding biological effects. A total of 127 studies were selected from 591 publications evaluating the toxic effects of environmental samples, including wastewater, surface water, and sediments. Here, bioassays used in EDA are summarized, including the assays that measure specific receptor-mediated modes of action (MOA), induction of xenobiotic metabolism pathways, and induction of adaptive stress response pathways using either in vitro or in vivo bioassays. Also, the identified substances using EDA are discussed based on their MOA. The importance of EDA in establishing a comprehensive approach for the detection of environmental contaminants using bioanalytical methods is emphasized. The current limitations and benefits of using EDA in practical applications are outlined and strategies for moving forward are proposed.

FT-GNN Tool for Bridging HRMS Features and Bioactivity: Uncovering Unidentified Estrogen Receptor Agonists in Sewage

Identifying primary estrogen receptor (ER) agonists in municipal sewage is essential for ensuring the health of aquatic environments. Given the complex and variable chemical composition of sewage, the predominant ER agonists remain unclear. High-resolution mass spectrometry (HRMS)-based models have been developed to predict compound bioactivity in complex matrices, but further optimization is needed to effectively bridge HRMS features with ER agonists. To address this challenge, an FT-GNN (fragmentation tree-based graph neural network) model was proposed. Given limited data and class imbalance, data augmentation was performed using model predictions within the applicability domain (AD) and oversampling technique (OTE). Model development results demonstrated that integrating the FT-GNN with data augmentation improved the balanced accuracy (bACC) value by 6%-31%. The developed model, with a high bACC to identify more true ER agonists, efficiently classified tens of thousands of unidentified HRMS features in sewage, reducing postprocessing workload in nontargeted screening. Analysis of ER agonist transformation during sewage treatment revealed the anaerobic stage as key to both their removal and formation. Estrogenic effect balance analysis suggests that α-E2 and 9,11-didehydroestriol may be two previously overlooked key ER agonists. Collectively, the development and application of the FT-GNN model are crucial advancements toward credible tracking and efficient control of estrogenic risks in water.

Novel Approach to Screen Endocrine-Disrupting Chemicals via Endocrine-Enhanced Reduced Human Transcriptome

Endocrine-disrupting chemicals (EDCs) can interfere with multiple pathways and trigger different modes of action. Thus, the traditional EDC in vitro screening processes often require a battery of bioassays to cover multiple target pathways. Here we developed an endocrine-enhanced reduced human transcriptome (ERHT) focused on hormone receptor signaling induced by the EDCs regulating specific genes. ERHT was developed based on 1200 prioritized genes covering 110 endocrine-related biological pathways across eight potential adverse outcomes. The ability of this approach to identify EDCs was derived from machine learning of 1068 dose-dependent transcriptome profiles and enhanced by quantifying chemical-induced critical pathway responses, and thus, it demonstrated excellent classification performance (AUC = 0.84 ± 0.03) in internal cross-validation. We ultimately applied this approach to known EDCs and inactive substances to validate the reliability of this approach. Through external validation on 210 chemicals, the extrapolation accuracy exceeded 80%, demonstrating the outstanding practical performance of this approach. Meanwhile, the pathway responses induced by the same chemical were consistent with the experimental results reported by multiple sequencing platforms, highlighting the robustness of this approach. The above results demonstrate that this approach can provide novel insights for EDCs' high-throughput screening and comprehensive toxic mechanisms through biological pathways.

Critical Structures of Bisphenol Analogues on Embryonic Toxicity Identified by a Computational Approach

Safer chemical alternatives to bisphenol (BP) have been a major pursuit of modern green chemistry and toxicology. Using a chemical similarity-based approach, it is difficult to identify minor structural differences that contribute to the significant changes of toxicity. Here, we used omics and computational toxicology to identify chemical features associated with BP analogue-induced embryonic toxicity, offering valuable insights to inform the design of safer chemical alternatives. The zebrafish embryonic acute toxicity, behavioral effects, and concentration-dependent transcriptome analysis of 17 BP analogues were tested, and the chemical structure characteristics and key biological activities-induced embryonic toxicity were explored. BPE, BPF, BPP, BPBP, and BPS induced lower embryonic lethality than BPA. And, 8 BP analogues triggered hyperactive behavior at environmentally and human relevant concentrations. BP analogues with phenol rings linked via hydrophobic segments ("chain:alkaneBranch_neopentyl_C5") disturbed stress response, leading to embryonic lethality, and introducing hydrophobic groups on the meta position of bisphenol structure augmented their embryonic lethality effects. "3DACorr_TotChg_3" of BP analogues is a key physicochemical feature for behavioral disorders, and BP analogues with 3DACorr_TotChg_3 value < 0.11 could induce hyperactive behavior by perturbing neurodevelopment relevant biological pathways. This study provides an integrated strategy, combining data-driven profiling and mechanism-based analysis for safer chemical alternatives.

A tiered toxicity testing strategy for assessing early life stage toxicity in estuarine fish (Mugilogobius chulae): A case study on tris (1-chloro-2-propyl) phosphate ester

The estuarine ecological environment faces significant threats from contaminants of emerging concern (CECs); yet, the risk posed by CECs to resident organisms remains poorly understood. Here, we employed tiered toxicity testing to investigate the adverse effects and potential mechanisms of tris (1-chloro-2-propyl) phosphate (TCPP) on the early life stages of an estuarine fish, Mugilogobius chulae. TCPP affected the development of M. chulae embryos, including survival, morphology, hatching, and behavior. A concentration-dependent transcriptomic analysis showed that TCPP disrupted 12 neurodevelopment-related KEGG pathways in M. chulae embryos, with five of the 30 % top-ranked pathways related to neurotransmitter signaling. Besides the cholinergic synapse signaling pathway, the glutamatergic signaling pathway (including NMDAR and AMPAR subtypes) may also mediate TCPP-induced neurodevelopmental toxicity. The NMDAR subtype GRIN2B was downregulated at high concentrations. Molecular dynamics simulations revealed a strong interaction between TCPP and GRIN2B, with TCPP binding to the residues Ile153 and Ile188. The results suggest that NMDARs play a crucial role in TCPP-induced neurodevelopmental toxicity toward M. chulae. AOP network analysis predicted that TCPP may impact cognitive functions and memory. Our study provides a novel testing strategy for identifying the mechanisms of toxicity of CECs, a crucial component of ecological risk assessment.

The impact of wastewater treatment plants on the composition and toxicity of pollutants in urban rivers in Nanjing, China

Widespread wastewater pollution is one of the biggest challenges threatening the ecological health of rivers. It is crucial to identify the toxic changes of effluents after entering urban rivers as well as the toxic substances in the complex chemical mixtures found in these urban rivers. This study used HepG2 cell line for cytotoxicity test to evaluate the ecological impact of effluents on urban rivers. Water samples were collected from the Xingwu River and Yunliang River in Nanjing, China. The bacterial communities in the lower reaches of urban rivers were altered due to the differences in total nitrogen and nitrate nitrogen. The complex chemical mixtures collected in the urban rivers were divided into 10 fractions, >100 chemicals were screened in each fraction. The substances with LC50 < 1000 mg/L were listed as toxic substances, and the number of toxic substances dominated the toxicity of urban rivers. Our study highlights toxicity as a comprehensive indicator for assessing river pollutants and reveals relationship between the number of toxic substance and river toxicity. These findings have direct implications for the monitoring and management of environmental stressors and the protection of aquatic organisms and human health.

Advancing the Effect-Directed Identification in Combined Pollution: Using Pathways to Link Effects and Toxicants

The difficulty in associating diverse pollutants with mixture effects has led to significant challenges in identifying toxicants in combined pollution. In this study, pathways were used to link effects and toxicants. By pathways evaluated by the concentration-dependent transcriptome, individual effects were extended to molecular mechanisms encompassing 135 pathways corresponding to 6 biological processes. Accordingly, mechanism-based identification of toxicants was achieved by constructing a pathway toxicant database containing 2413 chemical-pathway interactions and identifying pathway active fragments of 72 pathways. The developed method was applied to two different wastewaters, industrial wastewater OB and municipal wastewater HL. Although lethality and teratogenesis were both observed at the individual level, different molecular mechanisms were revealed by pathways, with cardiotoxicity- and genotoxicity-related pathways significantly enriched in OB, and neurotoxicity- and environmental information processing-related pathways significantly enriched in HL. Further suspect and nontargeted screening generated 59 and 86 causative toxicants in OB and HL, respectively, among which 29 toxicants were confirmed, that interacted with over 90% of enriched pathways and contributed over 50% of individual effects. After upgrading treatments based on causative toxicants, consistent removal of toxicants, pathway effects, and individual effects were observed. Mediation by pathways enables mechanism-based identification, supporting the assessment and management of combined pollution.

A novel molecular pathway of lipid accumulation in human hepatocytes caused by PFOA and PFOS

Exposed to ubiquitously perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS) has been associated with non-alcoholic fatty liver disease (NAFLD), yet the underlying molecular mechanism remains elusive. The extrapolation of empirical studies correlating per- and polyfluoroalkyl substance (PFAS) exposure with NAFLD occurrence to real-life exposure was hindered by the limited availability of mechanistic data at environmentally relevant concentrations. Herein, a novel pathway mediating hepatocyte lipid accumulation by PFOA and PFOS at human-relevant dose (<10 μM) was identified by integrating CRISPR-Cas9 genome screening, concentration-dependent transcriptional assay in HepG2 cell and epidemiological data mining. 1) At genetic level, nudt7 showed the highest enriched potency among 569 NAFLD-related genes, and the transcription of nudt7 was significantly downregulated by PFOA and PFOS exposure (<7 μM). 2) At molecular pathway, upon exposure to ≤10-4 μM PFOA and PFOS, the downregulation of nudt7 transcriptional expression triggered the reduction of Ace-CoA hydrolase activity. 3) At cellular level, increased lipids were measured in HepG2 cells with PFOA and PFOS (<2 μM). Overall, we identified a novel mechanism mediated by transcriptional downregulation of nudt7 gene in hepatocellular lipid increase treated with PFOA and PFOS, which could potentially explain the NAFLD occurrence associated with exposure to PFASs in humans.

Identification of Key Toxic Substances Considering Metabolic Activation: A Combination of Transcriptome and Nontarget Analysis

There have been numerous studies using effect-directed analysis (EDA) to identify key toxic substances present in source and drinking water, but none of these studies have considered the effects of metabolic activation. This study developed a comprehensive method including a pretreatment process based on an in vitro metabolic activation system, a comprehensive biological effect evaluation based on concentration-dependent transcriptome (CDT), and a chemical feature identification based on nontarget chemical analysis (NTA), to evaluate the changes in the toxic effects and differences in the chemical composition after metabolism. Models for matching metabolites and precursors as well as data-driven identification methods were further constructed to identify toxic metabolites and key toxic precursor substances in drinking water samples from the Yangtze River. After metabolism, the metabolic samples showed a general trend of reduced toxicity in terms of overall biological potency (mean: 3.2-fold). However, metabolic activation led to an increase in some types of toxic effects, including pathways such as excision repair, mismatch repair, protein processing in endoplasmic reticulum, nucleotide excision repair, and DNA replication. Meanwhile, metabolic samples showed a decrease (17.8%) in the number of peaks and average peak area after metabolism, while overall polarity, hydrophilicity, and average molecular weight increased slightly (10.3%). Based on the models for matching of metabolites and precursors and the data-driven identification methods, 32 chemicals were efficiently identified as key toxic substances as main contributors to explain the different transcriptome biological effects such as cellular component, development, and DNA damage related, including 15 industrial compounds, 7 PPCPs, 6 pesticides, and 4 natural products. This study avoids the process of structure elucidation of toxic metabolites and can trace them directly to the precursors based on MS spectra, providing a new idea for the identification of key toxic pollutants of metabolites.

Insights into the synergistic toxicity mechanisms caused by nano- and microplastics with triclosan using a dose-dependent functional genomics approach in Saccharomyces cerevisiae

The emergence of polystyrene (PS) nano- and microplastics (NMPs) and triclosan (TCS) as environmental contaminants has raised concerns about their combined toxicities to organisms, but the complex toxicity arising from their interactions and the underlying molecular mechanisms remain obscure to us. In this study, we comprehensively detected the combined toxicity of PS-NMPs and TCS via the dose-dependent yeast functional genomics profiling. Firstly, our findings demonstrated that the combined exposure to PS-NMPs and TCS elicited a synergistic toxic effect in which the toxicity depended on the size of the PS-NMPs. Secondly, we found that TCS exposure, either alone or in combination with PS-NMPs, influenced lipid biosynthetic processes and ATP export pathways, while the unique responsive genes triggered by combined exposure to TCS and PS-NMPs are significantly enriched in mitochondrial translation, ribosomal small subunit assembly, and tRNA wobble uridine modification. Thirdly, our results demonstrated that point of departure (POD) at the pathway level was positively correlated with IC50, and POD was a more sensitive predictor of toxicity than the apical toxicity endpoints. More importantly, our findings suggested that the combined exposure of PS-NMPs in a size-dependent manner not only alleviated the harmful effects of TCS on glycerophospholipid metabolism, but also exacerbated its negative impact on oxidative phosphorylation. Collectively, our study not only provides new insights into the intricate molecular mechanisms that control the combined toxicity of PS-NMPs and TCS, but also confirms the effectiveness of the dose-dependent functional genomics approach in elucidating the molecular mechanisms of the combined toxicity of pollutants.

High-efficiency effect-directed analysis (EDA) advancing toxicant identification in aquatic environments: Latest progress and application status

Facing the great threats to ecosystems and human health posed by the continuous release of chemicals into aquatic environments, effect-directed analysis (EDA) has emerged as a powerful tool for identifying causative toxicants. However, traditional EDA shows problems of low-coverage, labor-intensive and low-efficiency. Currently, a number of high-efficiency techniques have been integrated into EDA to improve toxicant identification. In this review, the latest progress and current limitations of high-efficiency EDA, comprising high-coverage effect evaluation, high-resolution fractionation, high-coverage chemical analysis, high-automation causative peak extraction and high-efficiency structure elucidation, are summarized. Specifically, high-resolution fractionation, high-automation data processing algorithms and in silico structure elucidation techniques have been well developed to enhance EDA. While high-coverage effect evaluation and chemical analysis should be further emphasized, especially omics tools and data-independent mass acquisition. For the application status in aquatic environments, high-efficiency EDA is widely applied in surface water and wastewater. Estrogenic, androgenic and aryl hydrocarbon receptor-mediated activities are the most concerning, with causative toxicants showing the typical structural features of steroids and benzenoids. A better understanding of the latest progress and application status of EDA would be beneficial to further advance in the field and greatly support aquatic environment monitoring.

AOP-Anchored Transcriptome Analysis Catalogue Accelerates the Discovery of Environmental Toxicants in Zebrafish

Current toxicity screening approaches to evaluate the vast number of environmental chemicals that require assessment are hampered due to their significant costs, time requirements, and reliance on live animal testing. The aim of the present study was to develop an adverse outcome pathway (AOP)-anchored transcriptome analysis (AATA) catalogue to expedite the discovery of environmental toxicants. 437 AOPs from the AOPwiki (https://aopwiki.org/) and 2280 transcriptomics data sets from NCBI Gene Expression Omnibus (GEO) and EMBL-EBI ArrayExpress (AE) repositories were comprehensively reviewed and analyzed. By using the differentially expressed molecular key event (mKE) genes as connection nodes, we created a large-scale environmental substance─target gene (mKE)─predicted adverse outcomes (SGAs) network that included 78 substances, 1099 genes, and 354 adverse outcomes (AOs). To validate the reliability of the network, comprehensive literature verification was conducted. We demonstrated that 164 of the 354 AOs identified have been previously characterized in the literature. The results for 136 of these AOs were consistent with the predictions of the AATA catalogue, representing an accuracy rate of 82.9%. Besides, distinct patterns in molecular KEs and AOs among categories of substances, such as biocides and metals, were demonstrated. Some representative substances, including atrazine and copper, pose significant risks to fish at various levels of biological organization. Moreover, experimental verification of the AATA predictions was conducted, including exposures of zebrafish to perfluorooctanesulfonate, cresyl diphenyl phosphate, and lanthanum. Results demonstrated consistency with predictions of the AATA catalogue, with an accuracy rate of 92.3%. Collectively, the present findings support the AATA catalogue as an efficient and promising platform for identifying environmental toxicants to fish and thereby provide novel insights into the understanding of potential risks of environmental contaminants.

High-Efficiency Effect-Directed Analysis Leveraging Five High Level Advancements: A Critical Review

Organic contaminants are ubiquitous in the environment, with mounting evidence unequivocally connecting them to aquatic toxicity, illness, and increased mortality, underscoring their substantial impacts on ecological security and environmental health. The intricate composition of sample mixtures and uncertain physicochemical features of potential toxic substances pose challenges to identify key toxicants in environmental samples. Effect-directed analysis (EDA), establishing a connection between key toxicants found in environmental samples and associated hazards, enables the identification of toxicants that can streamline research efforts and inform management action. Nevertheless, the advancement of EDA is constrained by the following factors: inadequate extraction and fractionation of environmental samples, limited bioassay endpoints and unknown linkage to higher order impacts, limited coverage of chemical analysis (i.e., high-resolution mass spectrometry, HRMS), and lacking effective linkage between bioassays and chemical analysis. This review proposes five key advancements to enhance the efficiency of EDA in addressing these challenges: (1) multiple adsorbents for comprehensive coverage of chemical extraction, (2) high-resolution microfractionation and multidimensional fractionation for refined fractionation, (3) robust in vivo/vitro bioassays and omics, (4) high-performance configurations for HRMS analysis, and (5) chemical-, data-, and knowledge-driven approaches for streamlined toxicant identification and validation. We envision that future EDA will integrate big data and artificial intelligence based on the development of quantitative omics, cutting-edge multidimensional microfractionation, and ultraperformance MS to identify environmental hazard factors, serving for broader environmental governance.

Bioactivity assessment of organophosphate flame retardants via a dose-dependent yeast functional genomics approach

Organophosphate flame retardants (OPFRs) have been widely detected in multiple environment media and have many adverse effects with complex toxicity mechanisms. However, the early molecular responses to OPFRs have not been fully elucidated, thereby making it difficult to assess their risks accurately. In this work, we systematically explored the point of departure (POD) of biological pathways at genome-wide level perturbed by 14 OPFRs with three substituents (alkyl, halogen, and aryl) using a dose-dependent functional genomics approach in Saccharomyces cerevisiae at 24 h exposure. Firstly, our results demonstrated that the overall biological potency at gene level (PODDRG20) ranged from 0.013 to 35.079 μM for 14 OPFRs, especially the tributyl phosphate (TnBP) exhibited the strongest biological potency with the least PODDRG20. Secondly, we found that structural characteristics of carbon number and logKow were significantly negatively correlated with POD, and carbon number and logKow also significantly affected lipid metabolism associated processes. Thirdly, these early biological pathways of OPFRs toxification were found to be involved in lipid metabolism, oxidative stress, DNA damage, MAPK signaling pathway, and amino acid and carbohydrate metabolism, among which the lipid metabolism was the most sensitive molecular response perturbed by most OPFRs. More importantly, we identified one resistant mutant strain with knockout of ERG2 (YMR202W) gene participated in steroid biosynthesis pathway, which can serve as a key yeast strain of OPFRs toxification. Overall, our study demonstrated an effective platform for accurately assessing OPFRs risks and provided a basis for further green OPFRs development.

Identify organic contaminants of high-concern based on non-targeted toxicity testing and non-targeted LC-HRMS analysis in tap water and source water along the Yangtze River

Many organic pollutants were detected in tap water (TW) and source water (SW) along the Yangtze River. However, the potential toxic effects and the high-concern organics (HCOs) which drive the effect are still unknown. Here, a non-targeted toxicity testing method based on the concentration-dependent transcriptome and non-targeted LC-HRMS analysis combining tiered filtering were used to reveal the overall biological effects and chemical information. Subsequently, we developed a qualitative pathway-structure relationship (QPSR) model to effectively match the biological and chemical information and successfully identified HCOs in TW and SW along the Yangtze River by potential substructures of HCOs. Non-targeted toxicity testing found that the biological potency of both TW and SW was stronger in the downstream of the Yangtze River, and disruption of the endocrine system and cancer were the main drivers of the effect. In addition, non-targeted LC-HRMS analysis combined with retention time prediction results identified 3220 and 631 high-confidence compound structures in positive and negative ion modes, respectively. Then, QPSR model was further implied and identified a total of 103 HCOs, containing 35 industrial chemicals, 30 PPCPs, 26 pesticides, and 12 hormones in TW and SW, respectively. Among them, the neuroactive and hormonal compounds oxoamide, 8-iso-16-cyclohexyl-tetranor prostaglandin E2, E Keppra, and Tocris-0788 showed the highest frequency of detection, which were identified in more than 1/3 of the samples. The strategy of combining non-targeted toxicity testing and non-targeted LC-HRMS analysis will support comprehensive biological effect assessment, identification of HCOs, and risk control of mixtures.

High-Throughput Transcriptomics of Water Extracts Detects Reductions in Biological Activity with Water Treatment Processes

The presence of numerous chemical contaminants from industrial, agricultural, and pharmaceutical sources in water supplies poses a potential risk to human and ecological health. Current chemical analyses suffer from limitations, including chemical coverage and high cost, and broad-coverage in vitro assays such as transcriptomics may further improve water quality monitoring by assessing a large range of possible effects. Here, we used high-throughput transcriptomics to assess the activity induced by field-derived water extracts in MCF7 breast carcinoma cells. Wastewater and surface water extracts induced the largest changes in expression among cell proliferation-related genes and neurological, estrogenic, and antibiotic pathways, whereas drinking and reclaimed water extracts that underwent advanced treatment showed substantially reduced bioactivity on both gene and pathway levels. Importantly, reclaimed water extracts induced fewer changes in gene expression than laboratory blanks, which reinforces previous conclusions based on targeted assays and improves confidence in bioassay-based monitoring of water quality.

Congener Variation of Genetic Dependent-Developmental Toxicology in Two Emerging Classes of Dioxin-like Compounds

Emerging classes of dioxin-like compounds (DLCs) like hydroxylated/methoxylated polybrominated diphenyl ethers (HO-/MeO-PBDEs) and polychlorinated diphenyl sulfides (PCDPSs) could lead to diverse adverse outcomes in humans and wildlife, yet knowledge gaps exist in their molecular mechanisms associated with different structures following early life environmental exposure. This study integrated a genetic knockout technique and concentration-dependent reduced zebrafish transcriptome approach (CRZT) to unravel the toxicological pathways underpinning developmental toxicity of four HO-/MeO-PBDEs and five PCDPSs at environmentally relevant doses. Generally, the dependence of aryl hydrocarbon receptor (AhR) on the embryotoxicity and transcriptomic potencies induced by the HO-PBDEs and PCDPSs varied across different congeners. The knockout of the ahr2 gene led to 1.02- to 76.48-fold decreases of DLC-induced embryotoxicities and reduced the transcriptome-based potencies ranging from 1.38 to 2124.74 folds in the CRZT test. The fold changes denoting AhR-mediated potentials significantly increased with the increasing chlorination degrees of MeO-PBDEs and PCDPSs (p < 0.05). Moreover, ahr2 knockout primarily affected the DLC-induced early molecular responses relevant to DNA damage, enzyme activation, and organ development. Our integrated approach revealed the differential role of AhR in mediating the developmental toxicity of emerging DLCs possessing varied structures at environmentally relevant doses.

Time-dependent hormetic dose responses of skin care product mixtures to Vibrio qinghaiensis sp.-Q67: Appearance and quantification

Hormesis is a widely recognized and extensively studied phenomenon. However, few studies have described the quantitative characteristics of hormesis required for appropriate risk assessment. Although skin care product (SCP) mixtures and their active ingredients can induce the hormesis of Vibrio qinghaiensis sp.-Q67 (Q67), the quantitative characteristics of time-dependent hormetic dose responses in SCPs have not yet been investigated. In this study, 28 SCP mixtures were tested for luminescence toxicity against Q67 after five exposure durations (0.25, 3, 6, 9, and 12 h). With increasing exposure duration, the concentration response curves (CRCs) were classified as constant monotonic nonlinear (S-shaped) for four SCPs, S- to hormetic (J-shaped) for 13 SCPs, and constant J-shaped for 11 SCPs. Of 140 CRCs, 98 were J-shaped. An increased frequency of SCPs inducing hormesis was observed. The toxicity (pEC50) of the SCPs was independent of the exposure duration and product type. The maximum stimulatory effect (Emin) of the 12 SCPs increased with exposure duration. We proposed a modified parameter, the width of inhibition dose zone (WIDZ; EC50/EC10), to depict the width of inhibition dose zone. The WIDZ of S-shaped CRCs were significantly larger than that of J-shaped CRCs. In addition, the characteristic parameters reported in the general literature were analyzed. The good linear relationship between EC50 and the maximum stimulatory effective concentration (ECmin) indicated that toxicity may be transformed into stimulatory effects over exposure durations. The width of stimulation dose zone (WSDZ) and Emin of the seven SCPs had the same increasing trends with increasing exposure duration. The combination of WIDZ with other characteristic parameters (e.g., zero effective concentration point, ECmin, etc.) could better depict hormesis with low-dose stimulation and high-dose inhibition. The quantitative characteristics of the dose-responses of hormesis-inducing SCPs could provide reference basis for the risk assessment of SCP mixtures.

Characterizing temporal variability and repeatability of dose-dependent functional genomics approach for evaluating triclosan toxification

Dose-dependent functional genomics approach has shown great advantage in identifying the molecular initiating event (MIE) of chemical toxification and yielding point of departure (POD) at genome-wide scale. However, POD variability and repeatability derived from experimental design (settings of dose, replicate number, and exposure time) has not been fully determined. In this work, we evaluated POD profiles perturbed by triclosan (TCS) using dose-dependent functional genomics approach in Saccharomyces cerevisiae at multiple time points (9 h, 24 h and 48 h). The full dataset (total 9 concentrations with 6 replicates per treatment) at 9 h was subsampled 484 times to generate subsets of 4 dose groups (Dose A - Dose D with varied concentration range and spacing) and 5 replicate numbers (2 reps - 6 reps). Firstly, given the accuracy of POD and the experimental cost, the POD profiles from 484 subsampled datasets demonstrated that the Dose C group (space narrow at high concentrations and wide dose range) with three replicates was best choice at both gene and pathway levels. Secondly, the variability of POD was found to be relatively robustness and stability across different experimental designs, but POD was more dependent on the dose range and interval than the number of replicates. Thirdly, MIE of TCS toxification was identified to be the glycerophospholipid metabolism pathway at all-time points, supporting the ability of our approach to accurately recognize MIE of chemical toxification at both short- and long-term exposure. Finally, we identified and validated 13 key mutant strains involved in MIE of TCS toxification, which could serve as biomarkers for TCS exposure. Taken together, our work evaluated the repeatability of dose-dependent functional genomics approach and the variability of POD and MIE of TCS toxification, which will benefit the experimental design for future dose-dependent functional genomics study.

Risk assessment based on dose-responsive and time-responsive genes to build PLS-DA models for exogenously induced lung injury

Xenobiotics can easily harm human lungs owing to the openness of the respiratory system. Identifying pulmonary toxicity remains challenging owing to several reasons: 1) no biomarkers for pulmonary toxicity are available that might help to detect lung injury; 2) traditional animal experiments are time-consuming; 3) traditional detection methods solely focus on poisoning accidents; 4) analytical chemistry methods hardly achieve universal detection. An in vitro testing system able to identify the pulmonary toxicity of contaminants from food, the environment, and drugs is urgently needed. Compounds are virtually infinite, whereas toxicological mechanisms are countable. Therefore, universal methods to identify and predict the risks of contaminants can be designed based on these well-known toxicity mechanisms. In this study, we established a dataset based on transcriptome sequencing of A549 cells upon treatment with different compounds. The representativeness of our dataset was analyzed using bioinformatics methods. Artificial intelligence methods, namely partial least squares discriminant analysis (PLS-DA) models, were employed for toxicity prediction and toxicant identification. The developed model predicted the pulmonary toxicity of compounds with a 92 % accuracy. These models were submitted to an external validation using highly heterogeneous compounds, which supported the accuracy and robustness of our developed methodology. This assay exhibits universal potential applications for water quality monitoring, crop pollution detection, food and drug safety evaluation, as well as chemical warfare agent detection.

Delving into the molecular initiating event of cadmium toxification via the dose-dependent functional genomics approach in Saccharomyces cerevisiae

Determining dose-response relationship is essential for comprehensively revealing chemical-caused effects on organisms. However, uncertainty and complexity of gene/protein interactions cause the inability of traditional toxicogenomic methods (e.g., transcriptomics, proteomics and metabolomics) to effectively establish the direct relationship between chemical exposure and genes. In this work, we built an effective dose-dependent yeast functional genomics approach, which can clearly identify the direct gene-chemical link in the process of cadmium (Cd) toxification from a genome-wide scale with wide range concentrations (0.83, 2.49, 7.48, 22.45, 67.34, 202.03 and 606.1 μM). Firstly, we identified 220 responsive strains, and found that 142, 110, 91, 34, 8, 0 and 0 responsive strains can be respectively modulated by seven different Cd exposure concentrations ranging from high to low. Secondly, our results demonstrated that these genes induced by the high Cd exposure were mainly enriched in the process of cell autophagy, but ones caused by the low Cd exposure were primarily involved in oxidative stress. Thirdly, we found that the top-ranked GO biological processes with the lowest point of departure (POD) were transmembrane transporter complex and mitochondrial respiratory chain complex III, suggesting that mitochondrion might be the toxicity target of Cd. Similarly, nucleotide excision repair was ranked first in KEGG pathway with the least POD, indicating that this dose-dependent functional genomics approach can effectively detect the molecular initiating event (MIE) of cadmium toxification. Fourthly, we identified four key mutant strains (RIP1, QCR8, CYT1 and QCR2) as biomarkers for Cd exposure. Finally, the dose-dependent functional genomics approach also performed well in identifying MIE for additional genotoxicity chemical 4-nitroquinoline-1-oxide (4-NQO) data. Overall, our study developed a dose-dependent functional genomics approach, which is powerful to delve into the MIE of chemical toxification and is beneficial for guiding further chemical risk assessment.

Prediction of zebrafish embryonic developmental toxicity by integrating omics with adverse outcome pathway

New approach methodologies (NAMs), especially omics-based high-throughput bioassays have been developed rapidly, providing rich mechanistic information such as molecular initiation events (MIEs) and (sub)cellular key events (KEs) in adverse outcome pathways (AOPs). However, how to apply the knowledge of MIEs/KEs to predict adverse outcomes (AOs) induced by chemicals represents a new challenge for computational toxicology. Here, an integrated method named ScoreAOP was developed and evaluated to predict chemicals' developmental toxicity for zebrafish embryos by integrating four related AOPs and dose-dependent reduced zebrafish transcriptome (RZT). The rules of ScoreAOP included 1) sensitivity of responsive KEs demonstrated by point of departure of KEs (POD_KE), 2) evidence reliability and 3) distance between KEs and AOs. Moreover, eleven chemicals with different modes of action (MoAs) were tested to evaluate ScoreAOP. Results showed that eight of the eleven chemicals caused developmental toxicity at tested concentration in apical tests. All the tested chemicals' developmental defects were predicted using ScoreAOP, whereas eight out of the eleven chemicals predicted by ScoreMIE which was developed to score MIEs disturbed by chemicals based on in vitro bioassays data. Finally, in terms of mechanism explanation, ScoreAOP clustered chemicals with different MoAs while ScoreMIE failed, and ScoreAOP revealed the activation of aryl hydrocarbon receptor (AhR) plays a significant role in dysfunction of cardiovascular system, resulting in zebrafish developmental defects and mortality. In conclusion, ScoreAOP represents a promising approach to apply mechanism information obtained from omics to predict AOs induced by chemicals.

Reduced Transcriptome Analysis of Zebrafish Embryos Prioritizes Environmental Compounds with Adverse Cardiovascular Activities

Cardiovascular diseases are the leading cause of premature death in humans and remain a global public health challenge. While age, sex, family history, and false nutrition make a contribution, our understanding of compounds acting as cardiovascular disruptors is far from complete. Here, we aim to identify cardiovascular disruptors via a reduced transcriptome atlas (RTA) approach, which integrates large-scale transcriptome data sets of zebrafish and compiles a specific gene panel related to cardiovascular diseases. Among 767 gene expression profiles covering 81 environmental compounds, 11 priority compounds are identified with the greatest effects on the cardiovascular system at the transcriptional level. Among them, metals (AgNO₃, Ag nanoparticles, arsenic) and pesticides/biocides (linuron, methylparaben, triclosan, and trimethylchlorotin) are identified with the most significant effects. Distinct transcriptional signatures are further identified by the percentage values, indicating that different physiological endpoints exist among prioritized compounds. In addition, cardiovascular dysregulations are experimentally confirmed for the prioritized compounds via alterations of cardiovascular physiology and lipid profiles of zebrafish. The accuracy rate of experimental verification reaches up to 62.9%. The web-based RTA analysis tool, Cardionet, for rapid cardiovascular disruptor discovery was further provided at http://www.envh.sjtu.edu.cn/cardionet.jsp. Our integrative approach yields an efficient platform to discover novel cardiovascular-disrupting chemicals in the environment.

Liquid Crystal Monomer: A Potential PPARγ Antagonist

Liquid crystal monomers (LCMs) are a large family of artificial ingredients that have been widely used in global liquid crystal display (LCD) industries. As a major constituent in LCDs as well as the end products of e-waste dismantling, LCMs are of growing research interest with regard to their environmental occurrences and biochemical consequences. Many studies have analyzed LCMs in multiple environmental matrices, yet limited research has investigated the toxic effects upon exposure to them. In this study, we combined in silico simulation and in vitro assay validation along with omics integration analysis to achieve a comprehensive toxicity elucidation as well as a systematic mechanism interpretation of LCMs for the first time. Briefly, the high-throughput virtual screen and reporter gene assay revealed that peroxisome proliferator-activated receptor gamma (PPARγ) was significantly antagonized by certain LCMs. Besides, LCMs induced global metabolome and transcriptome dysregulation in HK2 cells. Notably, fatty acid β-oxidation was conspicuously dysregulated, which might be mediated through multiple pathways (IL-17, TNF, and NF-kB), whereas the activation of AMPK and ligand-dependent PPARγ antagonism may play particularly important parts. This study illustrated LCMs as a potential PPARγ antagonist and explored their toxicological mode of action on the trans-omics level, which provided an insightful overview in future chemical risk assessment.

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 path via pathway-based approaches towards safety assessment: A concise review

For decades, chemical safety assessment has been proposed to shift from animal testing to in vitro testing systems in response to the call for the 3R. In Europe, the answer was to combine various information sources in integrated testing strategies (ITS); In the US, it was in 2007 when the landmark report by the National Research Council put forward a vision of in vitro toxicity testing paradigm. Since then, efforts to develop pathway-based assessment framework have been on the track. In 2010, systems biology brought out a conceptual framework called adverse outcome pathway (AOP), which took one step further from toxicity pathway to regulatory toxicology. Computational modeling, high-throughput screening, high-content omics have all been approached to facilitate this progress. This paper briefly reviewed the achievement of pathway-based chemical assessment since 2007, discussed potential pitfalls and challenges that mechanism-driven chemical assessment may undergo, and presented future perspectives of safety assessment that is to be based on computational system biology.

Evaluation of dioxin induced transcriptomic responses in a 3D human liver microtissue model

Three-dimensional human liver microtissue model provides a promising method for predicting the human hepatotoxicity of environmental chemicals. However, the dynamics of transcriptional responses of 3D human liver microtissue model to dioxins exposure remain unclear. Herein, time-series transcriptomic analysis was used to characterize modulation of gene expression over 14 days in 3D human liver microtissues exposed to 2,3,7,8-tetra-chlorodibenzo-p-dioxin (TCDD, 31 nM, 10 ng/ml). Changes in gene expression and modulation of biological pathways were evaluated at several time points. The results showed that microtissues stably expressed genes related to toxicological pathways (e.g. highly of genes involved in external stimuli and maintenance of cell homeostasis pathways) during the 14-day culture period. Furthermore, a weekly phenomenon pattern was observed for the number of the differentially expressed genes in microtissues exposed to TCDD at each time point. TCDD led to an induction of genes involved in cell cycle regulation at day three. Metabolic pathways were the main significantly induced pathways during the subsequent days, with the immune/inflammatory response enriched on the fifth day, and the cellular response to DNA damage was identified at the end of the exposure. Finally, relevant transcription patterns identified in microtissues were compared with published data on rodent and human cell-line studies to elucidate potential species-specific responses to TCDD over time. Cell development and cytochrome P450 pathway were mainly affected after a 3-day exposure, with the DNA damage response identified at the end of exposure in the human microtissue system but not in mouse/rat primary hepatocytes models. Overall, the 3D human liver microtissue model is a valuable tool to predict the toxic effects of environmental chemicals with a relatively long exposure.

Effect-Directed Analysis Based on the Reduced Human Transcriptome (RHT) to Identify Organic Contaminants in Source and Tap Waters along the Yangtze River

Since a large number of contaminants are detected in source waters (SWs) and tap waters (TWs), it is important to perform a comprehensive effect evaluation and key contributor identification. A reduced human transcriptome (RHT)-based effect-directed analysis, which consisted of a concentration-dependent RHT to reveal the comprehensive effects and noteworthy pathways and systematic identification of key contributors based on the interactions between compounds and pathway effects, was developed and applied to typical SWs and TWs along the Yangtze River. By RHT, 42% more differentially expressed genes and 33% more pathways were identified in the middle and lower reaches, indicating heavier pollution. Hormone and immune pathways were prioritized based on the detection frequency, sensitivity, and removal efficiency, among which the estrogen receptor pathway was the most noteworthy. Consistent with RHT, estrogenic effects were widespread along the Yangtze River based on in vitro evaluations. Furthermore, 38 of 100 targets, 39 pathway-related suspects, and 16 estrogenic nontargets were systematically identified. Among them, diethylstilbestrol was the dominant contributor, with the estradiol equivalent (EEQ) significantly correlated with EEQ_water. In addition, zearalenone and niclosamide explained up to 54% of the EEQ_water. The RHT-based EDA method could support the effect evaluation, contributor identification, and risk management of micropolluted waters.

Assessment of genotoxic chemicals using chemogenomic profiling based on gene-knockout library in Saccharomyces cerevisiae

Understanding the adverse effects of genotoxic chemicals and identifying them effectively from non-genotoxic chemicals are of great worldwide concerns. Here, Saccharomyces cerevisiae (yeast) genome-wide single-gene knockout screening approach was conducted to assess two genotoxic chemicals (4-nitroquinoline-1-oxide (4-NQO) and formaldehyde (FA)) and environmental pollutant dichloroacetic acid (DCA, genotoxicity is controversial). DNA repair was significant enriched in the gene ontology (GO) biology process (BP) terms and KEGG pathways when exposed to low concentrations of 4-NQO and FA. Higher concentrations of 4-NQO and FA influenced some RNA metabolic and biosynthesis pathways. Moreover, replication and repair associated pathways were top ranked KEGG pathways with high fold-change for low concentrations of 4-NQO and FA. The similar gene profiles perturbed by DCA with three test concentrations identified, the common GO BP terms associated with aromatic amino acid family biosynthetic process and ubiquitin-dependent protein catabolic process via the multivesicular body sorting pathway. DCA has no obvious genotoxicity as there was no enriched DNA damage and repair pathways and fold-change of replication and repair KEGG pathways were very low. Five genes (RAD18, RAD59, MUS81, MMS4, and BEM4) could serve as candidate genes for genotoxic chemicals. Overall, the yeast functional genomic profiling showed great performance for assessing the signatures and potential molecular mechanisms of genotoxic chemicals.

Cross-Model Comparison of Transcriptomic Dose-Response of Short-Chain Chlorinated Paraffins

Short-chain chlorinated paraffins (SCCPs) have attracted attention because of their toxicological potential in humans and wildlife at environmentally relevant doses. However, limited information is available regarding mechanistic differences across species in terms of the biological pathways that are impacted by SCCP exposure. Here, a concentration-dependent reduced human transcriptome (RHT) approach was conducted to evaluate 15 SCCPs in HepG2 cells and compared with our previous results using a reduced zebrafish transcriptome (RZT) approach in zebrafish embryos (ZFEs). Generally, SCCPs induced a broader suite of biological pathways in ZFEs than HepG2 cells, and all of the 15 SCCPs were more potent in HepG2 cells compared to ZFEs. Despite these general differences, the transcriptional potency of SCCPs in both model systems showed a significant linear relationship (p = 0.0017, r² = 0.57), and the average ratios of transcriptional potency for each SCCP in RZT to that in RHT were ∼100,000. C₁₀H₁₄Cl₈ was the most potent SCCP, while C₁₀H₁₇Cl₅ was the least potent in both ZFEs and HepG2 cells. An adverse outcome pathway network-based analysis demonstrated model-specific responses, such as xenobiotic metabolism that may be mediated by different nuclear receptor-mediated pathways between HepG2 cells (e.g., CAR and AhR activation) and ZFEs (e.g., PXR activation). Moreover, induced transcriptional changes in ZFEs associated with pathways and molecular initiating events (e.g., activation of nicotinic acetylcholine receptor) suggest that SCCPs may disrupt neural development processes. The cross-model comparison of concentration-dependent transcriptomics represents a promising approach to assess and prioritize SCCPs.

Assessment of fibrotic pathways induced by environmental chemicals using 3D-human liver microtissue model

Exposure to environmental chemicals, particularly those with persistent and bioaccumulative properties have been linked to liver diseases. Induction of fibrotic pathways is considered as a pre-requirement of chemical induced liver fibrosis. Here, we applied 3D in vitro human liver microtissues (MTs) composed of HepaRG, THP-1 and hTERT-HSC that express relevant hepatic pathways (bile acid, sterol, and xenobiotic metabolism) and can recapitulate key events of liver fibrosis (e.g. extracellular matrix-deposition). The liver MTs were exposed to a known profibrotic chemical, thioacetamide (TAA) and three representative environmental chemicals (TCDD, benzo [a] pyrene (BaP) and PCB126). Both TAA and BaP triggered fibrotic pathway related events such as hepatocellular damage (cytotoxicity and decreased albumin release), hepatic stellate cell activation (transcriptional upregulation of α-SMA and Col1α1) and extracellular matrix remodelling. TCDD or PCB126 at measured concentrations did not elicit these responses in the 3D liver MTs system, though they caused cytotoxicity in HepaRG monoculture at high concentrations. Reduced human transcriptome (RHT) analysis captured molecular responses involved in liver fibrosis when MTs were treated with TAA and BaP. The results suggest that 3D, multicellular, human liver microtissues represent an alternative, human-relevant, in vitro liver model for assessing fibrotic pathways induced by environmental chemicals.

Photodegradation of carbon dots cause cytotoxicity

Carbon dots (CDs) are photoluminescent nanomaterials with wide-ranging applications. Despite their photoactivity, it remains unknown whether CDs degrade under illumination and whether such photodegradation poses any cytotoxic effects. Here, we show laboratory-synthesized CDs irradiated with light degrade into molecules that are toxic to both normal (HEK-293) and cancerous (HeLa and HepG2) human cells. Eight days of irradiation photolyzes 28.6-59.8% of the CDs to <3 kilo Dalton molecules, 1431 of which are detected by high-throughput, non-target high-performance liquid chromatography-quadrupole time-of-flight mass spectrometry. Molecular network and community analysis further reveal 499 cytotoxicity-related molecules, 212 of which contain polyethylene glycol, glucose, or benzene-related structures. Photo-induced production of hydroxyl and alkyl radicals play important roles in CD degradation as affected by temperature, pH, light intensity and wavelength. Commercial CDs show similar photodegraded products and cytotoxicity profiles, demonstrating that photodegradation-induced cytotoxicity is likely common to CDs regardless of their chemical composition. Our results highlight the importance of light in cytocompatibility studies of CDs.

Carbon dots have attracted much attention for biomedical applications but potential degradation and associated toxicity are still poorly understood. Here, the authors report on a study into the photo-degradation of carbon dots, the products produced and associated cytotoxicity.

Metabolomic and Transcriptomic Analysis of MCF-7 Cells Exposed to 23 Chemicals at Human-Relevant Levels: Estimation of Individual Chemical Contribution to Effects

BACKGROUND

Humans are constantly being exposed to various xenobiotics at relatively low concentrations. To date, limited evidence is available to ascertain whether a complex xenobiotic mixture at human-relevant levels causes any health effect. Moreover, there is no effective method to pinpoint the contribution of each chemical toward such an effect.

OBJECTIVES

This study aims to understand the responses of cells to a mixture containing 23 xenobiotics at human-relevant levels and develop a feasible method to decipher the chemical(s) that contribute significantly to the observed effect.

METHODS

We characterized the metabolome and transcriptome of breast cancer cells (MCF-7) before and after exposure to the mixture at human-relevant levels; preexposure levels were derived from existing large-scale biomonitoring data. A high-throughput metabolomics-based "leave-one-out" method was proposed to understand the relative contribution of each component by comparing the metabolome with and without the particular chemical in the mixture.

RESULTS

The metabolomic analysis suggested that the mixture altered metabolites associated with cell proliferation and oxidative stress. For the transcriptomes, gene ontology terms and pathways including "cell cycle," "cell proliferation," and "cell division" were significantly altered after mixture exposure. The mixture altered genes associated with pathways such as "genotoxicity" and "nuclear factor erythroid 2-related factor 2 (Nrf2)." Through joint pathways analysis, metabolites and genes were observed to be well-aligned in pyrimidine and purine metabolisms. The leave-one-out results showed that many chemicals made their contributions to specific metabolic pathways. The overall metabolome pattern of the absence of 2,4-dihyroxybenzophenone (DHB) or bisphenol A (BPA) showed great resemblance to controls, suggesting their higher relative contribution to the observed effect.

DISCUSSION

The omics results showed that exposure to the mixture at human-relevant levels can induce significant in vitro cellular changes. Also, the leave one out method offers an effective approach for deconvoluting the effects of the mixture. https://doi.org/10.1289/EHP6641.

Molecular fingerprints of conazoles via functional genomic profiling of Saccharomyces cerevisiae

Conazoles were designed to inhibit ergosterol biosynthesis. Conazoles have been widely used as agricultural fungicides and are frequently detected in the environment. Although conazoles have been reported to have adverse effects, such as potential carcinogenic effects, the underlying molecular mechanisms of toxicity remain unclear. Here, the molecular fingerprints of five conazoles (propiconazole (Pro), penconazole (Pen), tebuconazole (Teb), flusilazole (Flu) and epoxiconazole (Epo)) were assessed in Saccharomyces cerevisiae (yeast) via functional genome-wide knockout mutant profiling. A total of 169 (4.49%), 176 (4.67%), 198 (5.26%), 218 (5.79%) and 173 (4.59%) responsive genes were identified at three concentrations (IC₅₀, IC₂₀ and IC₁₀) of Pro, Pen, Teb, Flu and Epo, respectively. The five conazoles tended to have similar gene mutant fingerprints and toxicity mechanisms. "Ribosome" (sce03010) and "cytoplasmic translation" (GO: 0002181) were the common KEGG pathway and GO biological process term by gene set enrichment analysis of the responsive genes, which suggested that conazoles influenced protein synthesis. Conazoles also affected fatty acids synthesis because "biosynthesis of unsaturated fatty acids" pathway was among the top-ranked KEGG pathways. Moreover, two genes, YGR037C (acyl-CoA-binding protein) and YCR034W (fatty acid elongase), were key fingerprints of conazoles because they played vital roles in conazole-induced toxicity. Overall, the fingerprints derived from the yeast functional genomic screening provide an alternative approach to elucidate the molecular mechanisms of environmental pollutant conazoles.

Impact of Mixture Effects between Emerging Organic Contaminants on Cytotoxicity: A Systems Biological Understanding of Synergism between Tris(1,3-dichloro-2-propyl)phosphate and Triphenyl Phosphate

Humans are exposed to many xenobiotics simultaneously, but little is known about the toxic effects based on chemical-chemical interactions. This study aims at evaluating the binary interactions between 13 common environmental organic compounds (resulting in 78 pairs) by observing their cytotoxicity on HepG2 cells. Among all of the tested pairs, the combination of flame-retardant triphenyl phosphate (TPP) and tris(1,3-dichloro-2-propyl)phosphate (TDCPP) exhibited one of the most significant synergistic effects. We further characterized the transcriptome and metabolome after combined exposure to TPP and TDCPP and individual exposure. The results suggested that the coexposure caused many more changes in gene expressions and cellular activities. The transcriptome data showed that the coexposure triggered significant pathway changes including "cholesterol biosynthesis" and "ATF6-Alpha activated chaperone genes", together with distinct gene ontology (GO) terms such as the "negative regulation of the ERK1 and ERK2 cascade". Additionally, coexposure enhanced the biological activity of liver X receptors and nuclear factor erythroid 2-related factor 2 (Nrf2). The metabolome data showed that coexposure significantly elevated oxidative stress and affected the purine and pyrimidine metabolism. Overall, this study showed that interactions, which may enhance or suppress the biological processes, are common among environmental chemicals, although their environmental relevance should be studied in the future.

Toward Sustainable Environmental Quality: Priority Research Questions for Asia

Environmental and human health challenges are pronounced in Asia, an exceptionally diverse and complex region where influences of global megatrends are extensive and numerous stresses to environmental quality exist. Identifying priorities necessary to engage grand challenges can be facilitated through horizon scanning exercises, and to this end we identified and examined 23 priority research questions needed to advance toward more sustainable environmental quality in Asia, as part of the Global Horizon Scanning Project. Advances in environmental toxicology, environmental chemistry, biological monitoring, and risk-assessment methodologies are necessary to address the adverse impacts of environmental stressors on ecosystem services and biodiversity, with Asia being home to numerous biodiversity hotspots. Intersections of the food-energy-water nexus are profound in Asia; innovative and aggressive technologies are necessary to provide clean water, ensure food safety, and stimulate energy efficiency, while improving ecological integrity and addressing legacy and emerging threats to public health and the environment, particularly with increased aquaculture production. Asia is the largest chemical-producing continent globally. Accordingly, sustainable and green chemistry and engineering present decided opportunities to stimulate innovation and realize a number of the United Nations Sustainable Development Goals. Engaging the priority research questions identified herein will require transdisciplinary coordination through existing and nontraditional partnerships within and among countries and sectors. Answering these questions will not be easy but is necessary to achieve more sustainable environmental quality in Asia. Environ Toxicol Chem 2020;39:1485-1505. © 2020 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

A Tiered Approach for Screening and Assessment of Environmental Mixtures by Omics and In Vitro Assays

New methodology approaches with a broad coverage of the biological effects are urgently needed to evaluate the safety of the universe of environmentally relevant chemicals. Here, we propose a tiered approach incorporating transcriptomics and in vitro bioassays to assess environmental mixtures. The mixture samples and the perturbed biological pathways are prioritized by concentration-dependent transcriptome (CDT) and then used to guide the selection of in vitro bioassays for toxicant identification. To evaluate omics' screening capability, we first applied a CDT technique to test mixture samples by HepG2 and MCF7 cells. The effect recoveries of large-volume solid-phase extraction on the overall bioactivity of the mixture were 48.9% in HepG2 and 58.3% in MCF7. The overall bioactivity potencies obtained by transcriptomics were positively correlated with the panel of 8 bioassays among 14 mixture samples combined with the previous data. Transcriptomics could predict their activation status (AUC = 0.783) and the relative potency (p < 0.05) of bioassays for four of the eight receptors (AhR, ER, AR, and Nrf2). Furthermore, the CDT identified other biological pathways perturbated by mixture samples, such as the pathway related to TP53, CAR, FXR, HIF, THRA, etc. Overall, this study demonstrates the potential of concentration-dependent omics for effect-based water quality assessment.

Toxicogenomics provides insights to toxicity pathways of neonicotinoids to aquatic insect, Chironomus dilutus

Neonicotinoid insecticides have posed a great threat to non-target organisms, yet the mechanisms underlying their toxicity are not well characterized. Major modes of action (MoAs) of imidacloprid were analyzed in an aquatic insect Chironomus dilutus. Lethal and sublethal outcomes were assessed in the midges after 96-h exposure to imidacloprid. Global transcriptomic profiles were determined using de novo RNA-sequencing to more holistically identify toxicity pathways. Transcriptional 10% biological potency values derived from ranked KEGG pathways and GO terms were 0.02 (0.01-0.08) (mean (95% confidence interval) and 0.05 (0.04-0.06) μg L^-1, respectively, which were more sensitive than those from phenotypic traits (10% lethal concentration: 0.44 (0.23-0.79) μg L^-1; 10% burrowing behavior concentration: 0.30 (0.22-0.43) μg L^-1). Major MoAs of imidacloprid in aquatic species were identified as follows: the activation of nicotinic acetylcholine receptors (nAChRs) induced by imidacloprid impaired organisms' nerve system through calcium ion homeostasis imbalance and mitochondrial dysfunction, which posed oxidative stress and DNA damage and eventually caused death of organisms. The current investigation highlighted that imidacloprid affected C. dilutus at environmentally relevant concentrations, and elucidated toxicity pathways derived from gene alteration to individual outcomes, calling for more attention to toxicity of neonicotinoids to aquatic organisms.

Transcriptome analyses unravel CYP1A1 and CYP1B1 as novel biomarkers for disinfection by-products (DBPs) derived from chlorinated algal organic matter

Disinfection by-products (DBPs) are generated during chlorination of drinking water. Previous studies demonstrate that DBPs are cytotoxic, genotoxic and associated with an increased risk of human cancer. However, the molecular basis of DBPs-induced toxic effects remains unclear. Here, we chlorinated samples of algal-derived organic matter (AOM) and sediment organic matter (SOM) from a local drinking water reservoir. Chemical properties, toxicities and transcriptomic profiles of human Caco-2 cell exposed to AOM and SOM were compared before and after chlorination. We analyzed chlorination-caused distinct gene expression patterns between AOM and SOM, and identified a set of 22 differentially expressed genes under chlorination of AOM that are different from chlorinated SOM. Consequent network analysis indicates that differential CYP1A1, CYP1B1, ID1 and ID2 are common targets of the upstream regulators predicted in the AOM group, but not the SOM group. Through experimental validation and data integration from previous reports related to DBPs or environmental stressors, we found that CYP1A1 and CYP1B1 are specifically up-regulated after chlorinating AOM. Our study demonstrates that the two CYP1 genes likely act as novel biomarkers of AOM derived DBPs, and this would be helpful for testing drinking water DBPs toxicity and further monitoring drinking water safety.

Concentration-dependent transcriptome of zebrafish embryo for environmental chemical assessment

Mechanistic information is essential to screen and predict the adverse effects of a large number of chemicals during early-life exposure. Concentration-dependent omics can capture the extent of perturbations of biological pathways or processes and provide information on the mechanism of toxicity. However, the application of concentration-dependent transcriptome to assess the developmental toxicity of environmental chemicals is still limited. Here, twelve chemicals representing five different modes of action (MOAs) were tested by the concentration-dependent reduced zebrafish transcriptome approach (CRZT) in combination with a phenotype-based high content screen (PHCS). The responsiveness, sensitivity and mechanistic differentiation of CRZT were validated in comparison with PHCS. First, PHCS identified 10 chemicals with obvious embryotoxicity (LD₅₀ range: 2.11-70.68 μM), while the potencies of the biological pathways perturbed by 12 chemicals (POD_path20 range: 0.002-2.1 μM) were demonstrated by CRZT. Second, although the potency of the transcriptome perturbations was positively correlated with lethality (LD₅₀) (R² = 0.64, P-value < 0.05) for most tested chemicals, BbF was non-embryotoxic but was the most potent on the perturbance of biological pathways. Finally, the profiles of the perturbed biological processes and the transcriptome potency (POD_path20) captured by CRZT could effectively classify most chemicals corresponding to their known MOAs. In summary, CRZT could significantly improve testing the developmental toxicity of environmental chemicals.

Pathway-based assessment of single chemicals and mixtures by a high-throughput transcriptomics approach

The ever-increasing number of chemicals and complex mixtures demands a high-throughput and cost-effective approach for chemical safety assessment. High-throughput transcriptomics (HTT) is promising in investigating genome-scale perturbation of chemical exposure in concentration-dependent manner. However, the application of HTT has been limited due to lack of methodology for single chemicals and mixture assessment. This study aimed to evaluate the ability of a newly-developed human reduced transcriptomics (RHT) approach to assess pathway-based profiles of single chemicals, and to develop a biological pathway-based approach for benchmarking mixture potency using single chemical-based prediction model. First, concentration-dependent RHT were used to qualitatively and quantitatively differentiate pathway-based patterns of different chemicals, using three model toxicants, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), triclosan (TCS) and 5-Chloro-6-hydroxy-2,2',4,4'-tetrabromodiphenyl ether (5-Cl-6-OH-BDE-47). AHR-regulated genes and pathways were most sensitively induced by TCDD, while TCS and 5-Cl-6-OH-BDE-47 were much less potent in AHR-associated activation, which was concordant with known MoA of each single chemical. Second, two artificial mixtures and their components of twelve individual chemicals were performed with concentration-dependent RHT. Concentration addition (CA) and independent action (IA) models were used to predict transcriptional potency of mixtures from transcriptomics of individual chemicals. For overall bioactivity, CA and IA models can both predict potency of observed responses within 95% confidence interval. For specific biological processes, multiple biological processes such as hormone signaling and DNA damage can be predicted using CA models for mixtures. The concentration-dependent RHT can provide a powerful approach for qualitative and quantitative assessment of biological pathway perturbated by environment chemical and mixtures.

Early Life Stage Bioactivity Assessment of Short-Chain Chlorinated Paraffins at Environmentally Relevant Concentrations by Concentration-Dependent Transcriptomic Analysis of Zebrafish Embryos

Short-chain chlorinated paraffins (SCCPs), a class of ubiquitous pollutants, are considered to be embryotoxic and teratogenic. However, little is known regarding the bioactivity and mechanisms at environmentally relevant concentrations at the embryonic period. Here, a concentration-dependent reduced transcriptomic approach was used to evaluate the environmental dose (<100 ppb) effects of nine SCCP congeners and eight commercial mixtures on zebrafish embryos at 8 hpf. After 24 h of exposure, the overall biological potency of all the SCCPs, in terms of interference with 20% of the differentially expressed genes (POD_DEG20), in zebrafish embryos ranged from 0.83 to 67.61 ppb. C₁₀H₁₄Cl₈ (POD_GO20 = 3.80 ppb) and C_10-13 51.5% Cl (POD_GO20 = 3.31 ppb) exhibited the strongest interference with biological processes compared to other SCCP homologs and mixtures, respectively. The most sensitive early molecular responses induced by SCCPs were associated with pathways of genetic damage, energy metabolite interference, and metal ion binding. Furthermore, the carbon number was positively correlated with the transcriptomic potency (POD_GO20) of SCCP congeners (with chlorine content > 60%) (p = 0.038), and the chlorine content of SCCP congeners affected the bioactivity associated with genotoxic pathways. The concentration-dependent reduced transcriptomic approach significantly improved the understanding of the ecological risk of environmental contaminants at early life stages.

Surface water extracts impair gene profiles and differentiation in human mesenchymal stem cells

Low concentrations of pollutants in surface water challenge the assessment of chronic effects on human health. Human bone mesenchymal stem cells (hBMSCs) were employed as a sensitive and relevant in vitro model to evaluate the potential biological effects caused by mixtures of pollutants in surface water. Organic extracts of surface water collected from Hun River inhibited cell viability in a dose-dependent manner. Surface water extracts at noncytotoxic concentrations induced 533 to 1055 differentially expressed genes (DEGs) in hBMSCs after 48 h of exposure. Total of 370 genes were commonly affected by surface water from different sites and accounted for 35-69% of DEGs impaired by individual sample. Pathways related to human diseases, genetic information processing and organismal systems were enriched based on DEGs. Interleukins (IL1B, IL6 and IL8) were affected and involved in most human diseases related pathways. The significantly downregulation of COL1A1 and the variation of rheumatoid arthritis pathway suggested that surface water potentially inhibited osteogenic differentiation of hBMSCs. Clustering analysis and principle component analysis with DEGs distinguish the surface water from tributary and mainstream. The crossing-species comparison of transcriptomic changes identified 923 and 2715 differentially expressed orthologs in hBMSCs and zebrafish, respectively. After the exposure ceased, the followed osteogenic and adipogenic differentiation in hBMSCs for 14 days were inhibited by the treatment of surface water during undifferentiated period, whereas the non-polar fraction exhibited stronger potency in affecting differentiation than the mid to polar fractions. hBMSCs, combining unsupervised transcriptomic technique and specific endpoints test, are promising in screening the health effects of environmental mixtures in surface water.

Assessing the impact of wastewater treatment plant effluent on downstream drinking water-source quality using a zebrafish (Danio Rerio) liver cell-based metabolomics approach

Cell-based metabolomics was used in a proof-of-concept fashion to investigate the biological effects of contaminants as they traveled from a wastewater treatment plant (WWTP) discharge to a drinking water treatment plant (DWTP) intake in a surface-water usage cycle. Zebrafish liver (ZFL) cells were exposed to water samples collected along a surface-water flowpath, where a WWTP was located ∼14.5 km upstream of a DWTP. The sampling sites included: 1) upstream of the WWTP, 2) the WWTP effluent discharging point, 3) a proximal location downstream of the WWTP outfall, 4) a distal location downstream of the WWTP outfall, 5) the drinking water intake, and 6) the treated drinking water collected prior to discharge to the distribution system. After a 48-h laboratory exposure, the hydrophilic and lipophilic metabolites in ZFL cell extracts were analyzed by proton nuclear magnetic resonance (¹H NMR) spectroscopy and gas chromatography-mass spectrometry (GC-MS), respectively. Multivariate statistical analysis revealed distinct changes in metabolite profiles in response to WWTP effluent exposure. These effects on the hydrophilic metabolome gradually diminished downstream of the WWTP, becoming non-significant at the drinking water intake (comparable to upstream of the WWTP, p = 0.98). However, effects on the lipophilic metabolome increased significantly as the river flowed from the distal location downstream of the WWTP to the drinking water intake (p < 0.001), suggesting a source of bioactive compounds in this watershed other than the WWTP. ZFL cells exposed to treated drinking water did not exhibit significant changes in either the hydrophilic (p = 0.15) or lipophilic metabolome (p = 0.83) compared to the upstream site, suggesting that constituents in the WWTP effluent were efficiently removed by the drinking water treatment process. Impacts on ZFL cells from the WWTP effluent included disrupted energy metabolism, a global decrease in amino acids, and altered lipid metabolism pathways. Overall, this study demonstrated the utility of cell-based metabolomics as an effective tool for assessing the biological effects of complex pollutant mixtures, particularly when used as a complement to conventional chemical monitoring.

Transcriptomic analysis of zebrafish (Danio rerio) embryos to assess integrated biotoxicity of Xitiaoxi River waters

Assessing the toxicity posed by mixtures of unknown chemicals to aquatic organisms is challenging. In this study, water samples from six cross-sections along the Xitiaoxi River Basin (XRB) were monthly or bimonthly collected in 2014. The year-period physiochemical parameters as well as one-month-water sample based acute biotoxicity tests showed that the river water quality of the year was generally in a good status. High performance liquid chromatography (HPLC) screening based on one-month-water samples suggested that the organic pollutants might be non-to-moderately-polar chemicals in very low concentrations. One-month-water sample based RNA-seq was performed to measure the mRNA differential expression profile of zebrafish larvae to furtherly explore the potential bioeffect and the spatial water quality change of the river. Result indicated that the number of deferentially expressed genes (DEGs) tended to increase along the downstream direction of the river. Gene ontology (GO) enrichment analysis implied that the key pollutants might mainly be the function disruptors of biological processes. Principle components analysis (PCA) combining with transcripts and one-month-water sample based physiochemical parameters indicated that the pollution might be similar at TP, DP and CTB sites while pollution homology existed on some extent between YBQ and JW sites. Although the water quality of the river had a complex time-space alternation during the year, and the one-month-data based RNA-seq could not reflex the whole year-water quality of a watershed, the gene expression profile via RNA-seq provided an alternative way for assessing integrated biotoxicity of surface water, and it was relatively fit for early-warning of water quality of a watershed with unobservable acute toxicity. However, the identification of detail toxicants and the links between DEGs and pollution level as well as physiological-biochemical toxicity needed further investigation.

In vitro bioanalysis of drinking water from source to tap

The presence of chemical pollutants in sources of drinking water is a key environmental problem threatening public health. Efficient removal of pollutants in drinking water treatment plants (DWTPs) is needed as well as methods for assessment of the total impact of all present chemicals on water quality. In the present study we have analyzed the bioactivity of water samples from source to tap, including effects of various water treatments in a DWTP, using a battery of cell-based bioassays, covering health-relevant endpoints. Reporter gene assays were used to analyze receptor activity of the aryl hydrocarbon receptor (AhR), estrogen receptor (ER), androgen receptor (AR), peroxisome proliferator-activated receptor alpha (PPARα) and induction of oxidative stress by the nuclear factor erythroid 2-related factor 2 (Nrf2). DNA damage was determined by Comet assay. Grab water samples were concentrated by HLB or ENV solid phase extraction and the water samples assayed at a relative enrichment factor of 50. The enrichment procedure did not induce any bioactivity. No bioactivity was detected in Milli-Q water or drinking water control samples. Induction of AhR, ER and Nrf2 activities was revealed in source to tap water samples. No cytotoxicity, PPARα or AR antagonist activity, or DNA damage were observed in any of the water samples. A low AR agonist activity was detected in a few samples of surface water, but not in the samples from the DWTP. The treatment steps at the DWTP, coagulation, granulated activated carbon filtration, UV disinfection and NH₂Cl dosing had little or no effect on the AhR, Nrf2 and ER bioactivity. However, nanofiltration and passage through the distribution network drastically decreased AhR activity, while the effect on Nrf2 activity was more modest and no apparent effect was observed on ER activity. The present results suggest that bioassays are useful tools for evaluation of the efficiency of different treatment steps in DWTPs in reducing toxic activities. Bioassays of AhR and Nrf2 are useful for screening of effects of a broad range of chemicals in drinking water and ER activity can be monitored with a high sensitivity.

Reduced Zebrafish Transcriptome Atlas toward Understanding Environmental Neurotoxicants

Transcriptomic approaches monitoring gene responses at genome-scale are increasingly used in toxicological research and help to clarify the molecular mechanisms of adverse effects caused by environmental toxicants. However, their applications for chemical assessment are hampered due to high expenses required and more importantly the lack of in-depth data mining and mechanistic perspectives. Here, we described a reduced transcriptome atlas (RTA) approach which integrates transcriptomic data sets and a comprehensive panel of genes generated to represent neurogenesis and the early neuronal development of zebrafish, to determine the potential neurodevelopmental toxicities of environmental chemicals. Transcriptomic data sets of 74 chemicals and 736 related gene expression profiles were integrated resulting in 135 exposure signatures. Chemical prioritization demonstrated four sets of hits to be neurotoxic: neuro-active chemicals (representatively, Valproic acid, VPA and Carbamazepine, CAR), xenoestrogens (Bisphenol A, BPA; Genistein, GEN; 17-α ethinylestradiol, EE2), microcystins (cyanopeptolin, CP1020; microcystin-LR, MCLR) and heavy metals (AgNO₃, AgNPs). The enriched biological pathways and processes were distinct among the four sets, while the overlapping functional enrichments were observed within each set, for example, over 25% differentially expressed genes and four of top five KEGG pathways were shared between VPA and CAR. Furthermore, gene expression index (GEI) analysis demonstrated that a gene panel with 300 genes was sufficient to effectively characterize and cluster chemicals and therefore offer an efficient and cost-effective tool for the prioritization of neurotoxicants. Thus, the RTA approach provides novel insights into the understanding of the in-depth molecular mechanisms of environmental neurotoxicants and can be used as an indication for potential adverse outcomes.

Omics Advances in Ecotoxicology

Toxic substances in the environment generate adverse effects at all levels of biological organization from the molecular level to community and ecosystem. Given this complexity, it is not surprising that ecotoxicologists have struggled to address the full consequences of toxic substance release at ecosystem level, due to the limits of observational and experimental tools to reveal the changes in deep structure at different levels of organization. -Omics technologies, consisting of genomics and ecogenomics, have the power to reveal, in unprecedented detail, the cellular processes of an individual or biodiversity of a community in response to environmental change with high sample/observation throughput. This represents a historic opportunity to transform the way we study toxic substances in ecosystems, through direct linkage of ecological effects with the systems biology of organisms. Three recent examples of -omics advance in the assessment of toxic substances are explored here: (1) the use of functional genomics in the discovery of novel molecular mechanisms of toxicity of chemicals in the environment; (2) the development of laboratory pipelines of dose-dependent, reduced transcriptomics to support high-throughput chemical testing at the biological pathway level; and (3) the use of eDNA metabarcoding approaches for assessing chemical effects on biological communities in mesocosm experiments and through direct observation in field monitoring. -Omics advances in ecotoxicological studies not only generate new knowledge regarding mechanisms of toxicity and environmental effect, improving the relevance and immediacy of laboratory toxicological assessment, but can provide a wholly new paradigm for ecotoxicology by linking ecological models to mechanism-based, systems biology approaches.

A Reduced Transcriptome Approach to Assess Environmental Toxicants Using Zebrafish Embryo Test

Omics approaches can monitor responses and alterations of biological pathways at genome-scale, which are useful to predict potential adverse effects by environmental toxicants. However, high throughput application of transcriptomics in chemical assessment is limited due to the high cost and lack of "standardized" toxicogenomic methods. Here, a reduced zebrafish transcriptome (RZT) approach was developed to represent the whole transcriptome and to profile bioactivity of chemical and environmental mixtures in zebrafish embryo. RZT gene set of 1637 zebrafish Entrez genes was designed to cover a wide range of biological processes, and to faithfully capture gene-level and pathway-level changes by toxicants compared with the whole transcriptome. Concentration-response modeling was used to calculate the effect concentrations (ECs) of DEGs and corresponding molecular pathways. To validate the RZT approach, quantitative analysis of gene expression by RNA-ampliseq technology was used to identify differentially expressed genes (DEGs) at 32 hpf following exposure to seven serial dilutions of reference chemical BPA (10-10E^-5μM) or each of four water samples ranging from wastewater to drinking water (relative enrichment factors 10-6.4 × 10^-4). The RZT-ampliseq-embryo approach was both sensitive and able to identify a wide spectrum of biological activities associated with BPA exposure. Water quality was benchmarked based on the sensitivity distribution curve of biological pathways detected using RZT-ampliseq-embryo. Finally, the most sensitive biological pathways were identified, including those linked with adverse reproductive outcomes, genotoxicity and development outcomes. RZT-ampliseq-embryo approach provides an efficient and cost-effective tool to prioritize toxicants based on responsiveness of biological pathways.