Omics Advances in Ecotoxicology

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.

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.

A ubiquitous tire rubber additive induced serious eye injury in zebrafish (Danio rerio)

Previous studies have indicated that tire wear particles (TWPs) leachate exposure induced serious eye injury in fish through inhibiting the thyroid peroxidase (TPO) enzyme activity. However, the main TPO inhibitors in the leachate were still unknown. In this study, we identified 2-Mercaptobenzothiazole (MBT) as the potential TPO inhibitor in the TWPs leachate through references search, model prediction based on Danish QSAR and ToxCast database, molecular docking, and in vivo assay. We further explored the toxic mechanism of MBT under environmentally relevant concentrations. The decreased eye size of zebrafish larvae was mainly caused by the decreased lens diameter and cell density in the inner nuclear layer (INL) and outer nuclear layer (ONL) of the retina. Transcriptomics analysis demonstrated that the eye phototransduction function was significantly suppressed by inhibiting the photoreceptor cell proliferation process after MBT exposure. The altered opsin gene expression and decreased opsin protein levels were induced by weakening thyroid hormone signaling after MBT treatment. These results were comparable to those obtained from a known TPO inhibitor, methimazole. This study has identified MBT as the primary TPO inhibitor responsible for inducing eye impairment in zebrafish larvae exposed to TWPs leachate. It is crucial for reducing the toxicity of TWPs leachate in fish.

Characterization of the transcriptional effects of the plastic additive dibutyl phthalate alone and in combination with microplastic on the green-lipped mussel Perna canaliculus

The presence and persistence of microplastics (MPs) in diverse aquatic environments are of global concern. Microplastics can impact marine organisms via direct physical interaction and the release of potentially harmful chemical additives incorporated into the plastic. These chemicals are physically bound to the plastic matrix and can leach out. The hazards associated with chemical additives to exposed organisms is not well characterized. We investigated the hazards of plastic additives leaching from plastic. We used the common plasticizer dibutyl phthalate (DBP) as a chemical additive proxy and the New Zealand green-lipped mussel (Perna canaliculus) as a model. We used early-adult P. canaliculus exposed to combinations of virgin and DBP-spiked polyvinyl chloride (PVC), MPs, and DBP alone for 7 days. Whole transcriptome sequencing (RNA-seq) was conducted to assess whether leaching of DBP from MPs poses a hazard. The differences between groups were evaluated using pairwise permutational multivariate analysis of variance (PERMANOVA), and all treatments were significantly different from controls. In addition, a significant difference was seen between DBP and PVC MP treatment. Transcriptome analysis revealed that mussels exposed to DBP alone had the most differentially expressed genes (914), followed by PVC MP + DBP (448), and PVC MP (250). Gene ontology functional analysis revealed that the most enriched pathway types were in cellular metabolism, immune response, and endocrine disruption. Microplastic treatments enriched numerous pathways related to cellular metabolism and immune response. The combined exposure of PVC MP + DBP appears to cause combined effects, suggesting that DBP is bioavailable to the exposed mussels in the PVC MP + DBP treatment. Our results support the hypothesis that chemical additives are potentially an important driver of MP toxicity. Environ Toxicol Chem 2024;43:1604-1614. © 2024 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

A multigenerational study can detect the evolutionary response to BaP exposure in the non-biting freshwater midge Chironomus riparius

The release of polycyclic aromatic hydrocarbons (PAHs) into the environment is posing a threat to ecosystems and human health. Benzo(a)pyrene (BaP) is considered a biomarker of PAH exposure and is classified as a Group 1 carcinogen. However, it was not known whether BaP is mutagenic, i.e. induces inherited germline mutations. In this study, we used a recently established method, which combines short-term mutation accumulation lines (MAL) with whole genome sequencing (WGS) to assess mutagenicity in the non-biting midge Chironomus riparius. The mutagenicity analysis was supplemented by an evaluation of the development of population fitness in three successive generations in the case of chronic exposure to BaP at a high concentration (100 μg/L). In addition, the level of ROS-induced oxidative stress was examined in vivo. Exposure to the higher BaP concentration led to an increase in germline mutations relative to the control, while the lower concentration showed no mentionable effect. Against expectations, BaP exposure decreased ROS-level compared to the control and is thus probably not responsible for the increased mutation rate. Likewise, the higher BaP concentration decreased fitness measured as population growth rate per day (PGR) significantly over all generations, without signs of rapid evolutionary adaptations. Our results thus highlighted that high BaP exposure may influence the evolutionary trajectory of organisms.

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.

Exploring the Relationship Among Lipid Profile Changes, Growth, and Reproduction in Folsomia candida Exposed to Teflubenzuron Over Time

The integration of untargeted lipidomics approaches in ecotoxicology has emerged as a strategy to enhance the comprehensiveness of environmental risk assessment. Although current toxicity tests with soil microarthropods focus on species performance, that is, growth, reproduction, and survival, understanding the mechanisms of toxicity across all levels of biological organization, from molecule to community is essential for informed decision-making. Our study focused on the impacts of sublethal concentrations of the insecticide teflubenzuron on the springtail Folsomia candida. Untargeted lipidomics was applied to link changes in growth, reproduction, and the overall stress response with lipid profile changes over various exposure durations. The accumulation of teflubenzuron in organisms exposed to the highest test concentration (0.035 mg a.s. kg-1 soil dry wt) significantly impacted reproductive output without compromising growth. The results suggested a resource allocation shift from reproduction to size maintenance. This hypothesis was supported by lipid shifts on day 7, at which point reductions in triacylglycerol and diacylglycerol content corresponded with decreased offspring production on day 21. The hypermetabolism of fatty acids and N-acylethanolamines on days 2 and 7 of exposure indicated oxidative stress and inflammation in the animals in response to teflubenzuron bioaccumulation, as measured using high-performance liquid chromatography-tandem mass spectrometry. Overall, the changes in lipid profiles in comparison with phenotypic adverse outcomes highlight the potential of lipid analysis as an early-warning tool for reproductive disturbances caused by pesticides in F. candida. Environ Toxicol Chem 2024;00:1-12. © 2024 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Effects of ambient UVB light on Pacific oyster Crassostrea gigas mantle tissue based on multivariate data

Ambient ultraviolet radiation (UVB) from solar and artificial light presents serious environmental risks to aquatic ecosystems. The Pacific oyster, Crassostrea gigas, perceives changes in the external environment primarily through its mantle tissue, which contains many nerve fibers and tentacles. Changes within the mantles can typically illustrate the injury of ambient UVB. In this study, a comprehensive analysis of phenotypic, behavioral, and physiological changes demonstrated that extreme UVB radiation (10 W/m²) directly suppressed the behavioral activities of C. gigas. Conversely, under ambient UVB radiation (5 W/m²), various physiological processes exhibited significant alterations in C. gigas, despite the behavior remaining relatively unaffected. Using mathematical model analysis, the integrated analysis of the full-length transcriptome, proteome, and metabolome showed that ambient UVB significantly affected the metabolic processes (saccharide, lipid, and protein metabolism) and cellular biology processes (autophagy, apoptosis, oxidative stress) of the C. gigas mantle. Subsequently, using Procrustes analysis and Pearson correlation analysis, the association between multi-omics data and physiological changes, as well as their biomarkers, revealed the effect of UVB on three crucial biological processes: activation of autophagy signaling (key factors: Ca2+, LC3B, BECN1, caspase-7), response to oxidative stress (reactive oxygen species, heat shock 70, cytochrome c oxidase), and recalibration of energy metabolism (saccharide, succinic acid, translation initiation factor IF-2). These findings offer a fresh perspective on the integration of multi-data from non-model animals in ambient UVB risk assessment.

Exploring ecological effects of arsenic and cadmium combined exposure on cropland soil: from multilevel organisms to soil functioning by multi-omics coupled with high-throughput quantitative PCR

Arsenic (As) and cadmium (Cd) pose potential ecological threats to cropland soils; however, few studies have investigated their combined effects on multilevel organisms and soil functioning. Here, we used collembolans and soil microbiota as test organisms to examine their responses to soil As and Cd co-contamination at the gene, individual, and community levels, respectively, and further uncovered ecological relationships between pollutants, multilevel organisms, and soil functioning. At the gene level, collembolan transcriptome revealed that elevated As concentrations stimulated As-detoxifying genes AS3MT and GST, whereas the concurrent Cd restrained GST gene expression. At the individual level, collembolan reproduction was sensitive to pollutants while collembolan survival wasn't. At the community level, significant but inconsistent correlations were observed between the biodiversity of different soil keystone microbial clusters and soil As levels. Moreover, soil functioning related to nutrient (e.g., carbon, nitrogen, phosphorus, and sulfur) cycles was inhibited under As and Cd co-exposure only through the mediation of plant pathogens. Overall, these findings suggested multilevel bioindicators (i.e., AS3MT gene expression in collembolans, collembolan reproduction, and biodiversity of soil keystone microbial clusters) in cropland soils co-contaminated with As and Cd, thus improving the understanding of the ecotoxicological impact of heavy metal co-contamination on soil ecosystems.

Considerations for the inclusion of metabarcoding data in the plant protection product risk assessment of the soil microbiome

There is increasing interest in further developing the plant protection product (PPP) environmental risk assessment, particularly within the European Union, to include the assessment of soil microbial community composition, as measured by metabarcoding approaches. However, to date, there has been little discussion as to how this could be implemented in a standardized, reliable, and robust manner suitable for regulatory decision-making. Introduction of metabarcoding-based assessments of the soil microbiome into the PPP risk assessment would represent a significant increase in the degree of complexity of the data that needs to be processed and analyzed in comparison to the existing risk assessment on in-soil organisms. The bioinformatics procedures to process DNA sequences into community compositional data sets currently lack standardization, while little information exists on how these data should be used to generate regulatory endpoints and the ways in which these endpoints should be interpreted. Through a thorough and critical review, we explore these challenges. We conclude that currently, we do not have a sufficient degree of standardization or understanding of the required bioinformatics and data analysis procedures to consider their use in an environmental risk assessment context. However, we highlight critical knowledge gaps and the further research required to understand whether metabarcoding-based assessments of the soil microbiome can be utilized in a statistically and ecologically relevant manner within a PPP risk assessment. Only once these challenges are addressed can we consider if and how we should use metabarcoding as a tool for regulatory decision-making to assess and monitor ecotoxicological effects on soil microorganisms within an environmental risk assessment of PPPs. Integr Environ Assess Manag 2024;20:337-358. © 2023 SETAC.

Reporting and reproducibility: Proteomics of fish models in environmental toxicology and ecotoxicology

Environmental toxicology and ecotoxicology research efforts are employing proteomics with fish models as New Approach Methodologies, along with in silico, in vitro and other omics techniques to elucidate hazards of toxicants and toxins. We performed a critical review of toxicology studies with fish models using proteomics and reported fundamental parameters across experimental design, sample preparation, mass spectrometry, and bioinformatics of fish, which represent alternative vertebrate models in environmental toxicology, and routinely studied animals in ecotoxicology. We observed inconsistencies in reporting and methodologies among experimental designs, sample preparations, data acquisitions and bioinformatics, which can affect reproducibility of experimental results. We identified a distinct need to develop reporting guidelines for proteomics use in environmental toxicology and ecotoxicology, increased QA/QC throughout studies, and method optimization with an emphasis on reducing inconsistencies among studies. Several recommendations are offered as logical steps to advance development and application of this emerging research area to understand chemical hazards to public health and the environment.

The posterity of Zebrafish in paradigm of in vivo molecular toxicological profiling

The aggrandised advancement in utility of advanced day-to-day materials and nanomaterials has raised serious concern on their biocompatibility with human and other biotic members. In last few decades, understanding of toxicity of these materials has been given the centre stage of research using many in vitro and in vivo models. Zebrafish (Danio rerio), a freshwater fish and a member of the minnow family has garnered much attention due to its distinct features, which make it an important and frequently used animal model in various fields of embryology and toxicological studies. Given that fertilization and development of zebrafish eggs take place externally, they serve as an excellent model organism for studying early developmental stages. Moreover, zebrafish possess a comparable genetic composition to humans and share almost 70% of their genes with mammals. This particular model organism has become increasingly popular, especially for developmental research. Moreover, it serves as a link between in vitro studies and in vivo analysis in mammals. It is an appealing choice for vertebrate research, when employing high-throughput methods, due to their small size, swift development, and relatively affordable laboratory setup. This small vertebrate has enhanced comprehension of pathobiology and drug toxicity. This review emphasizes on the recent developments in toxicity screening and assays, and the new insights gained about the toxicity of drugs through these assays. Specifically, the cardio, neural, and, hepatic toxicology studies inferred by applications of nanoparticles have been highlighted.

A novel multi-omics data analysis of dose-dependent and temporal changes in regulatory pathways due to chemical perturbation: a case study on caffeine

Comprehensive analysis of multi-omics data can reveal alterations in regulatory pathways induced by cellular exposure to chemicals by characterizing biological processes at the molecular level. Data-driven omics analysis, conducted in a dose-dependent or dynamic manner, can facilitate comprehending toxicity mechanisms. This study introduces a novel multi-omics data analysis designed to concurrently examine dose-dependent and temporal patterns of cellular responses to chemical perturbations. This analysis, encompassing preliminary exploration, pattern deconstruction, and network reconstruction of multi-omics data, provides a comprehensive perspective on the dynamic behaviors of cells exposed to varying levels of chemical stimuli. Importantly, this analysis is adaptable to any number of omics layers, including site-specific phosphoproteomics. We implemented this analysis on multi-omics data obtained from HepG2 cells exposed to a range of caffeine doses over varying durations and identified six response patterns, along with their associated biomolecules and pathways. Our study demonstrates the effectiveness of the proposed multi-omics data analysis in capturing multidimensional patterns of cellular response to chemical perturbation, enhancing understanding of pathway regulation for chemical risk assessment.

Effects of sublethal fipronil exposure on cross-generational functional responses and gene expression in Binodoxys communis

The effective systemic insecticide fipronil is widely used on a variety of crops and in public spaces to control insect pests. Binodoxys communis (Gahan) (Hymenoptera: Braconidae) is the dominant natural enemy of Aphis gossypii Glover (Homoptera: Aphididae), an important cotton pest, and has good efficiency in inhibiting aphid populations. The direct effects of environmental residues of sublethal fipronil doses on adult B. communis have not previously been reported. This study therefore aimed to evaluate the side effects and transcriptomic impacts of sublethal fipronil doses on B. communis. The results showed that exposure to the LC10 dose of fipronil significantly reduced the survival rate and parasitism rate of the F0 generation, but did not affect these indicators in the F1 generation. The LC25 dose did not affect the survival or parasitic rates of the F0 generation, but did significantly reduce the survival rate of F1 generation parasitoids. These results indicated that sublethal doses of fipronil affected B. communis population growth. Transcriptome analysis showed that differentially expressed genes (DEGs) in B. communis at 1 h after treatment were primarily enriched in pathways associated with fatty acid elongation, biosynthesis of fatty acids, and fatty acid metabolism. DEGs at 3 days after treatment were mainly enriched in ribosomal functions, glycolysis/gluconeogenesis, and tyrosine metabolism. Six DEGs (PY, ELOVL, VLCOAR, MRJP1, ELOVL AAEL008004-like, and RPL13) were selected for validation with real-time fluorescent quantitative PCR. This is the first report of sublethal, trans-generational, and transcriptomic side effects of fipronil on the dominant parasitoid of A. gossypii. The results of this study show that adaptation of parasitoids to high concentrations of pesticides may be at the expense of their offspring. These findings broaden our overall understanding of the intergenerational adjustments used by insects to respond to pesticide stress and call for risk assessments of the long-term impacts and intergenerational effects of other pesticides.

Insights in Pharmaceutical Pollution: The Prospective Role of eDNA Metabarcoding

Environmental pollution is a growing threat to natural ecosystems and one of the world's most pressing concerns. The increasing worldwide use of pharmaceuticals has elevated their status as significant emerging contaminants. Pharmaceuticals enter aquatic environments through multiple pathways related to anthropogenic activity. Their high consumption, insufficient waste treatment, and the incapacity of organisms to completely metabolize them contribute to their accumulation in aquatic environments, posing a threat to all life forms. Various analytical methods have been used to quantify pharmaceuticals. Biotechnology advancements based on next-generation sequencing (NGS) techniques, like eDNA metabarcoding, have enabled the development of new methods for assessing and monitoring the ecotoxicological effects of pharmaceuticals. eDNA metabarcoding is a valuable biomonitoring tool for pharmaceutical pollution because it (a) provides an efficient method to assess and predict pollution status, (b) identifies pollution sources, (c) tracks changes in pharmaceutical pollution levels over time, (d) assesses the ecological impact of pharmaceutical pollution, (e) helps prioritize cleanup and mitigation efforts, and (f) offers insights into the diversity and composition of microbial and other bioindicator communities. This review highlights the issue of aquatic pharmaceutical pollution while emphasizing the importance of using modern NGS-based biomonitoring actions to assess its environmental effects more consistently and effectively.

Application of CRISPR/Cas9-based genome editing in ecotoxicology

Chemicals are widely used and released into the environment, and their degradation, accumulation, migration, and transformation processes in the environment can pose a threat to the ecosystem. The advancement in analytical methods with high-throughput screening of biomolecules has revolutionized the way toxicologists used to explore the effects of chemicals on organisms. CRISPR/Cas is a newly developed tool, widely used in the exploration of basic science and biologically engineered products given its high efficiency and low cost. For example, it can edit target genes efficiently, and save loss of the crop yield caused by environmental pollution as well as gain a better understanding of the toxicity mechanisms from various chemicals. This review briefly introduces the development history of CRISPR/Cas and summarizes the current application of CRISPR/Cas in ecotoxicology, including its application on improving crop yield and drug resistance towards agricultural pollution, antibiotic pollution and other threats. The benefits by applying the CRISPR/Cas9 system in conventional toxicity mechanism studies are fully demonstrated here together with its foreseeable expansions in other area of ecotoxicology. Finally, the prospects and disadvantages of CRISPR/Cas system in the field of ecotoxicology are also discussed.

Toxicity effects of pesticides based on zebrafish (Danio rerio) models: Advances and perspectives

Pesticides inevitably enter aquatic environments, posing potential risks to organisms. The common aquatic model organism, zebrafish (Danio rerio), are widely used to evaluate the toxicity of pesticides. In this review, we searched the Web of Science database for articles published between 2012 and 2022, using the keywords "pesticide", "zebrafish", and "toxicity", retrieving 618 publications. Furthermore, we described the main pathways by which pesticides enter aquatic environments and the fate of their residues in these environments. We systematically reviewed the toxicity effects of pesticides on zebrafish, including developmental toxicity, endocrine-disrupting effects, reproductive toxicity, neurotoxicity, immunotoxicity, and genotoxicity. Importantly, we summarized the latest research progress on the toxicity mechanism of pesticides to zebrafish based on omics technologies, including transcriptomics, metabolomics, and microbiomics. Finally, we discussed future research prospects, focusing on the combined exposure of multiple pollutants including pesticides, the risk of multigenerational exposure to pesticides, and the chronic toxicity of aquatic nanopesticides. This review provides essential data support for ecological risk assessments of pesticides in aquatic environments, and has implications for water management in the context of pesticide pollution.

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.

Microbially mediated fossil concretions and their characterization by the latest methodologies: a review

The study of well-preserved organic matter (OM) within mineral concretions has provided key insights into depositional and environmental conditions in deep time. Concretions of varied compositions, including carbonate, phosphate, and iron-based minerals, have been found to host exceptionally preserved fossils. Organic geochemical characterization of concretion-encapsulated OM promises valuable new information of fossil preservation, paleoenvironments, and even direct taxonomic information to further illuminate the evolutionary dynamics of our planet and its biota. Full exploitation of this largely untapped geochemical archive, however, requires a sophisticated understanding of the prevalence, formation controls and OM sequestration properties of mineral concretions. Past research has led to the proposal of different models of concretion formation and OM preservation. Nevertheless, the formation mechanisms and controls on OM preservation in concretions remain poorly understood. Here we provide a detailed review of the main types of concretions and formation pathways with a focus on the role of microbes and their metabolic activities. In addition, we provide a comprehensive account of organic geochemical, and complimentary inorganic geochemical, morphological, microbial and paleontological, analytical methods, including recent advancements, relevant to the characterization of concretions and sequestered OM. The application and outcome of several early organic geochemical studies of concretion-impregnated OM are included to demonstrate how this underexploited geo-biological record can provide new insights into the Earth's evolutionary record. This paper also attempts to shed light on the current status of this research and major challenges that lie ahead in the further application of geo-paleo-microbial and organic geochemical research of concretions and their host fossils. Recent efforts to bridge the knowledge and communication gaps in this multidisciplinary research area are also discussed, with particular emphasis on research with significance for interpreting the molecular record in extraordinarily preserved fossils.

Are New Phthalate Ester Substitutes Safer than Traditional DBP and DiBP? Comparative Endocrine-Disrupting Analyses on Zebrafish Using In Vivo, Transcriptome, and In Silico Approaches

Although previous studies have confirmed the association between phthalate esters (PAEs) exposure and endocrine disorders in humans, few studies to date have systematically assessed the threats of new PAE alternatives to endocrine disruptions. Herein, zebrafish embryos were continuously exposed to two PAEs [di-n-butyl phthalate (DBP) and diisobutyl phthalate (DiBP)], two structurally related alternatives [diiononyl phthalate (DINP) and diisononyl hexahydrophthalate (DINCH)], and two non-PAE substitutes [dipropylene glycol dibenzoate (DGD) and glyceryl triacetate (GTA)], and the endocrine-disrupting effects were investigated during the early stages (8-48 hpf). For five endogenous hormones, including progesterone, testosterone, 17β-estradiol, triiodothyronine (T3), and cortisol, the tested chemicals disturbed the contents of at least one hormone at environmentally relevant concentrations (≤3.9 μM), except DINCH and GTA. Then, the concentration-dependent reduced zebrafish transcriptome analysis was performed. Thyroid hormone (TH)- and androgen/estrogen-regulated adverse outcome pathways (AOPs) were the two types of biological pathways most sensitive to PAE exposure. Notably, six compounds disrupted four TH-mediated AOPs, from the inhibition of deiodinases (molecular initiating event, MIE), a decrease in T3 levels (key event, KE), to mortality (adverse outcome, AO) with the quantitatively linear relationships between MIE-KE (|r| = 0.96, p = 0.002), KE-AO (|r| = 0.88, p = 0.02), and MIE-AO (|r| = 0.89, p = 0.02). Multiple structural analyses showed that benzoic acid is the critical toxicogenic fragment. Our data will facilitate the screening and development of green alternatives.

Mechanistic insight into the adverse outcome of tire wear and road particle leachate exposure in zebrafish (Danio rerio) larvae

Tire wear particles (TWP) have become the major microplastic pollution in China. Road runoff containing TWP leachate can decrease the eye size and even induced mortality in the aquatic organisms. However, the toxic mechanism of TWP and road particles (RP) leachate on aquatic organisms is still unclear. In this study, the zebrafish embryos were exposed to TWP or RP leachate for 5 days at both environmental relevant and high concentrations. The adverse outcome pathways (AOPs) were screened from individual to molecular levels. The morphological and behavioral analysis demonstrated that the leachate exposure mainly impaired the eye development of zebrafish larvae and inhibited the larval swim behavior and phototactic response, which are the adverse outcomes. The phototransduction modulated by zebrafish retina was significantly down-regulated through transcriptomics and metabolomics analysis. The eye histopathological analysis showed that the decreased thickness of the retinal outer nuclear layer (ONL) and retinal pigmented epithelium (RPE) after leachate exposure were caused by the decreased photoreceptor cells. Moreover, the expression of NR2E3 and TPO genes showed concentration-dependent down-regulation after leachate exposure. The inhibition of photoreceptor cell proliferation was identified as the main reason for photoreceptor cell decrease in zebrafish larval eye. This study, for the first time, uncovered the underlying toxic mechanism of TWP and RP on zebrafish larval eyes.

Ecogenotoxicity assessment with land snails: A mini-review

In the context of the increasing environmental and sanitary crisis, it is accepted that soil pollution can cause health alterations and disturb natural population dynamics. Consequently, the assessment of the genotoxic potential of compounds found in contaminated soils is important. Indeed, the alteration of genomic integrity may increase the risk of cancer development and may impair reproduction and long-term population dynamics. Among the methodologies to assess terrestrial genotoxic potential, there has been growing interest during the last decade in monitoring alterations of the genome in bioindicators of soil quality. As some land snail species are recognized bioindicators of soil quality, especially to assess the environmental and toxicological bioavailability of compounds, this review focuses on current knowledge regarding the genotoxicology of land snails. Classical biomarkers to assess genotoxic effects have been used (e.g., DNA breakage, micronuclei, random amplification polymorphic DNA) at various stages of the life cycle, including embryos. The studies were performed in vitro, in vivo, in situ and ex situ and covered a diverse set of contaminants (nanoparticles, metal(loid)s, pesticides, polycyclic aromatic hydrocarbons) and snail species (Cantareus aspersus, Eobania vermiculata, Theba pisana, Helix lucorum). Based on recent studies reviewed here, the use of land snails to map soil genotoxic potential is promising due to their ability to reveal pollution and subsequent environmental risks. Moreover, the position of snails in the trophic chain and the existing bridges between contaminant bioavailability to snails and bioaccessibility to humans reinforce the value of land snail-based ecotoxicological assessment.

Current Status of Omics in Biological Quality Elements for Freshwater Biomonitoring

Freshwater ecosystems have been experiencing various forms of threats, mainly since the last century. The severity of this adverse scenario presents unprecedented challenges to human health, water supply, agriculture, forestry, ecological systems, and biodiversity, among other areas. Despite the progress made in various biomonitoring techniques tailored to specific countries and biotic communities, significant constraints exist, particularly in assessing and quantifying biodiversity and its interplay with detrimental factors. Incorporating modern techniques into biomonitoring methodologies presents a challenging topic with multiple perspectives and assertions. This review aims to present a comprehensive overview of the contemporary advancements in freshwater biomonitoring, specifically by utilizing omics methodologies such as genomics, metagenomics, transcriptomics, proteomics, metabolomics, and multi-omics. The present study aims to elucidate the rationale behind the imperative need for modernization in this field. This will be achieved by presenting case studies, examining the diverse range of organisms that have been studied, and evaluating the potential benefits and drawbacks associated with the utilization of these methodologies. The utilization of advanced high-throughput bioinformatics techniques represents a sophisticated approach that necessitates a significant departure from the conventional practices of contemporary freshwater biomonitoring. The significant contributions of omics techniques in the context of biological quality elements (BQEs) and their interpretations in ecological problems are crucial for biomonitoring programs. Such contributions are primarily attributed to the previously overlooked identification of interactions between different levels of biological organization and their responses, isolated and combined, to specific critical conditions.

The influence of temperature rise on the metabolic response of Ruditapes philippinarum clams to 17-α-ethinylestradiol

Untargeted Nuclear Magnetic Resonance metabolomics was employed to study the effects of warming conditions (17-21 °C) and exposure to 17-α-ethinylestradiol (EE2) on the polar metabolome of Ruditapes philippinarum clams, to identify metabolic markers for monitoring/prediction of deviant environmental conditions. Warming alone triggered changes in alanine/aspartate/glutamate, aromatic amino acids, taurine/hypotaurine and homarine/trigonelline pathways, as well as in energy metabolism, suggesting osmoregulatory adaptations and glycolytic/tricarboxylic acid (TCA) cycle activation, possibly accompanied to some extent by gluconeogenesis to preserve glycogen reserves. At 17 °C, the lowest EE2 concentration (5 ng/L) specifically engaged branched-chain and aromatic amino acids to activate the glycolysis/TCA cycle. Notably, a partial metabolic recovery was observed at 25 ng/L, whereas higher EE2 concentrations (125 and 625 ng/L) again induced significant metabolic disturbances. These included enhanced glycogen biosynthesis and increased lipid reserves, sustained by low-level glutathione-based antioxidative mechanisms that seemed active. At 21 °C, response to EE2 was notably weak at low/intermediate concentrations, becoming particularly significant at the highest EE2 concentration (625 ng/L), suggesting higher protection capacity of Ruditapes philippinarum clams under warming conditions. At 625 ng/L, disturbances in alanine/aspartate/glutamate and taurine/hypotaurine metabolisms were observed, with no evidence of enhanced carbohydrate/protein catabolism. This low energy function profile was accompanied by marked antioxidative mechanisms and choline compounds modulation for cell membrane protection/repair. These results help monitor clams´ response to temperature rise and EE2 exposure, paving the way for future effective guidance and prediction of environmental damaging effects.

Comparison between droplet digital PCR and reverse transcription-quantitative PCR methods to measure ecotoxicology biomarkers

Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is currently the gold-standard technique for detecting and quantifying messenger RNA. However, without proper validation, the method may produce artefactual and non-reproducible cycle threshold values generating poor-quality data. The newer droplet digital PCR (ddPCR) method allows for the absolute quantification of targeted nucleic acids providing more sensitive and accurate measurements without requiring external standards. This study compared these two PCR-based methods to measure the expression of well-documented genes used in ecotoxicology studies. We exposed Mediterranean mussels (Mytilus galloprovincialis) to copper and analyzed gene expression in gills and digestive glands using RT-qPCR and ddPCR assays. A step-by-step methodology to optimize and compare the two technologies is described. After ten-fold serial complementary DNA dilution, both RT-qPCR and ddPCR exhibited comparable linearity and efficiency and produced statistically similar results. We conclude that ddPCR is a suitable method to assess gene expression in an ecotoxicological context. However, RT-qPCR has a shorter processing time and remains more cost-effective.

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.

Freshwater crustacean exposed to active pharmaceutical ingredients: ecotoxicological effects and mechanisms

Concerns over the ecotoxicological effects of active pharmaceutical ingredients (APIs) on aquatic invertebrates have been raised in the last decade. While numerous studies have reported the toxicity of APIs in invertebrates, no attempt has been made to synthesize and interpret this dataset in terms of different exposure scenarios (acute, chronic, multigenerational), multiple crustacean species, and the toxic mechanisms. In this study, a thorough literature review was performed to summarize the ecotoxicological data of APIs tested on a range of invertebrates. Therapeutic classes including antidepressants, anti-infectives, antineoplastic agents, hormonal contraceptives, immunosuppressants, and neuro-active drugs exhibited higher toxicity to crustaceans than other API groups. The species sensitivity towards APIs exposure is compared in D. magna and other crustacean species. In the case of acute and chronic bioassays, ecotoxicological studies mainly focus on the apical endpoints including growth and reproduction, whereas sex ratio and molting frequency are commonly used for evaluating the substances with endocrine-disrupting properties. The multigenerational and "Omics" studies, primarily transcriptomics and metabolomics, were confined to a few API groups including beta-blocking agents, blood lipid-lowing agents, neuroactive agents, anticancer drugs, and synthetic hormones. We emphasize that in-depth studies on the multigenerational effects and the toxic mechanisms of APIs on the endocrine systems of freshwater crustacean are warranted.

Transcriptomic hallmarks of in vitro TiO2 nanotubes toxicity in Chlamydomonas reinhardtii

Recently, more accessible transcriptomic approaches have provided a new and deeper understanding of environmental toxicity. The present study focuses on the transcriptomic profiles of green microalgae Chlamydomonas reinhardtii exposed to new industrially promising material, TiO₂ nanotubes (NTs), as an example of a widely used one-dimensional nanomaterial. The first algal in vitro assay included 2.5 and 7.5 mg/L TiO₂ NTs, resulting in a dose-dependent negative effect on biological endpoints. At a working concentration of 7.5 mg/L, RNA-sequencing showed a mainly negative effect on the cells. In summary, the results indicated metabolic disruption, ^{such as ATP loss, damage to mitochondria and chloroplasts, loss of solutes due to permeated membranes,} and cell wall damage. Moreover, apoptosis-induced transcripts were detected. Interestingly, reactivation of transposons was observed. In signalling and transcription pathways, including chromatin remodelling and locking, the annotated genes were downregulated.

Daphnia as a model organism to probe biological responses to nanomaterials-from individual to population effects via adverse outcome pathways

The importance of the cladoceran Daphnia as a model organism for ecotoxicity testing has been well-established since the 1980s. Daphnia have been increasingly used in standardised testing of chemicals as they are well characterised and show sensitivity to pollutants, making them an essential indicator species for environmental stress. The mapping of the genomes of D. pulex in 2012 and D. magna in 2017 further consolidated their utility for ecotoxicity testing, including demonstrating the responsiveness of the Daphnia genome to environmental stressors. The short lifecycle and parthenogenetic reproduction make Daphnia useful for assessment of developmental toxicity and adaption to stress. The emergence of nanomaterials (NMs) and their safety assessment has introduced some challenges to the use of standard toxicity tests which were developed for soluble chemicals. NMs have enormous reactive surface areas resulting in dynamic interactions with dissolved organic carbon, proteins and other biomolecules in their surroundings leading to a myriad of physical, chemical, biological, and macromolecular transformations of the NMs and thus changes in their bioavailability to, and impacts on, daphnids. However, NM safety assessments are also driving innovations in our approaches to toxicity testing, for both chemicals and other emerging contaminants such as microplastics (MPs). These advances include establishing more realistic environmental exposures via medium composition tuning including pre-conditioning by the organisms to provide relevant biomolecules as background, development of microfluidics approaches to mimic environmental flow conditions typical in streams, utilisation of field daphnids cultured in the lab to assess adaption and impacts of pre-exposure to pollution gradients, and of course development of mechanistic insights to connect the first encounter with NMs or MPs to an adverse outcome, via the key events in an adverse outcome pathway. Insights into these developments are presented below to inspire further advances and utilisation of these important organisms as part of an overall environmental risk assessment of NMs and MPs impacts, including in mixture exposure scenarios.

Increasing Temperatures Potentiate the Damage of Rare Earth Element Yttrium to the Crop Plant Triticum aestivum L

Given that increasing temperature may aggravate the toxicity of pollutants, it is a daunting challenge to evaluate the realistic risks of rare earth elements (REEs) under global warming. Here, we studied how elevated temperatures (27 and 32 °C) impact the effect of yttrium (Y) on wheat plants (Triticum aestivum L.) at concentrations not causing effects (0, 0.5, and 1 μM) at the control temperature (22 °C) in a hydroponic system. After 14 days of exposure, significant inhibition (p < 0.05, 29.5%) of root elongation was observed only at 1 μM of Y at 32 °C. Exposure to Y at 27 °C showed no visible effects on root length, but induced significant (p < 0.05) metabolic disorders of a range of carbohydrates and amino acids related to galactose, phenylalanine, and glutamate metabolisms. Such cases were even shifted to substantial perturbation of the nucleotide pool reallocation involved in the disruption of purine and pyrimidine metabolism at 32 °C. These observations were regulated by sets of genes involved in these perturbed pathways. Using weighted gene co-expression network analysis, the disorder of nucleotide metabolism was shown to be responsible for the aggravated Y phytotoxicity at the extreme high temperature. Although the temperature fluctuation considered seems to be in an extreme range, unexpected implications driven by high temperature cannot be neglected. Our findings thus reduce the gaps of knowledge in REE toxicity to plants under future climate warming scenarios and highlight the importance of incorporating environmental temperature into the framework of the risk assessment of REEs.

Molecular fingerprints of polar narcotic chemicals based on heterozygous essential gene knockout library in Saccharomyces cerevisiae

Cytotoxicity of non-polar narcotic chemicals can be predicted by quantitative structure activity relationship (QSAR) models, but the polar narcotic chemicals' actual cytotoxicity exceeds the predicted values by their chemical structures. This discrepancy indicates that the molecular mechanism by which polar narcotic chemicals exert their toxicity is unclear. Taking advantage of Saccharomyces cerevisiae (yeast) functional genome-wide heterozygous essential gene knockout mutants, we here have identified the specific molecular fingerprints of two main chemical structure groups (phenols and anilines) of polar narcotic chemicals (dichlorophen (DCP), 4-chlorophenol (4-CP), 2, 4, 6-trichlorophenol (TCP), 3, 4-dichloroaniline (DCA) and N-methylaniline (NMA)) and one non-polar narcotic chemical 2, 2, 2-trichloroethanol (TCE). Especially, we identify 33, 57, 54, 46, 59 and 53 responsive strains through exposure to TCE, DCP, 4-CP, TCP, DCA and NMA with three test concentrations, respectively, revealing that these polar narcotic chemicals have more responsive strains than the non-polar narcotic chemical. Remarkably, we find that the molecular fingerprints of polar narcotic chemicals in different chemical structure groups are obviously varied, particularly phenols and anilines have their own specific molecular fingerprints. Interestingly, our results demonstrate that the molecular toxicity mechanisms of anilines are associated with DNA replication, but phenols are related with pathway of RNA degradation. Additionally, we find that the two knockout strains (SME1 and DIS3) and the three knockout strains (TSC11, RSP5 and HSF1) can specifically respond to exposure to phenols and anilines, respectively. Thus, they may be served as potential biomarkers to distinguish phenols from anilines. Collectively, our works demonstrate that the functional genomic platform of yeast essential gene mutants can not only act as an effective tool to identify key specific molecular fingerprints for polar narcotic chemicals, but also help to understand the molecular mechanisms of polar narcotic chemicals.

Application of proteomics and metabolomics to assess ammonia stress response and tolerance mechanisms of juvenile ornate rock lobster Panulirus ornatus

Ammonia is a common pollutant in the aquatic ecosystem and closed aquaculture systems. It may pose a threat to the lobster growth, reproduction and survival. However, there is lack of research of the mechanisms on the toxic effects ammonia at molecular levels. In this work, proteomics and metabolomics were integrated to analyze the proteome and metabolome responses in the ornate spiny lobster Panulirus ornatus treated with ammonia (20 mg L^-1) for 48 h. A total of 199 proteins and 176 metabolites were significantly altered in P. ornatus following ammonia treatment. The responsive proteins and metabolites were predominantly involved in immune response, phase I and phase II biotransformation, carbohydrate metabolism, amino acid metabolism, and lipid metabolism. Furthermore, an increase in urea levels was observed, and amino acid metabolism was induced, indicating that the urea cycle was utilized to biotransform ammonia so as to reduce endogenous ammonia content. Ammonia exposure also affected the antioxidant system and induced cellular apoptosis. Overall, our results provide comprehensive insights into the molecular mechanisms underlying the response of P. ornatus to ammonia stress. We believe that the data reported herein should contribute to the development of novel, efficient methods for P. ornatus aquaculture.

Advances in Single-Cell Toxicogenomics in Environmental Toxicology

The toxicity evaluation system of environmental pollutants has undergone numerous changes due to the application of new technologies. Single-cell toxicogenomics is rapidly changing our view on environmental toxicology by increasing the resolution of our analysis to the level of a single cell. Applications of this technology in environmental toxicology have begun to emerge and are rapidly expanding the portfolio of existing technologies and applications. Here, we first summarized different methods involved in single-cell isolation and amplification in single-cell sequencing process, compared the advantages and disadvantages of different methods, and analyzed their development trends. Then, we reviewed the main advances of single-cell toxicogenomics in environmental toxicology, emphatically analyzed the application prospects of this technology in identifying the target cells of pollutants in early embryos, clarifying the heterogeneous response of cell subtypes to pollutants, and finding pathogenic bacteria in unknown microbes, and highlighted the unique characteristics of this approach with high resolution, high throughput, and high specificity by examples. We also offered a prediction of the further application of this technology and the revolution it brings in environmental toxicology. Overall, these advances will provide practical solutions for controlling or mitigating exogenous toxicological effects that threaten human and ecosystem health, contribute to improving our understanding of the physiological processes affected by pollutants, and lead to the emergence of new methods of pollution control.

Critical features identification for chemical chronic toxicity based on mechanistic forecast models

Facing billions of tons of pollutants entering the ocean each year, aquatic toxicity is becoming a crucial endpoint for evaluating chemical adverse effects on ecosystems. Notably, huge amount of toxic chemicals at environmental relevant doses can cause potential adverse effects. However, chronic aquatic toxicity effects of chemicals are much scarcer, especially at population level. Rotifers are highly sensitive to toxicants even at chronic low-doses and their communities are usually considered as effective indicators for assessing the status of aquatic ecosystems. Therefore, the no observed effect concentration (NOEC) for population abundance of rotifers were selected as endpoints to develop machine learning models for the prediction of chemical aquatic chronic toxicity. In this study, forty-eight binary models were built by eight types of chemical descriptors combined with six machine learning algorithms. The best binary model was 1D & 2D molecular descriptors - random trees model (RT) with high balanced accuracy (BA) (0.83 for training and 0.83 for validation set), and Matthews correlation coefficient (MCC) (0.72 for training set and 0.67 for validation set). Moreover, the optimal model identified the primary factors (SpMAD_Dzp, AMW, MATS2v) and filtered out three high alerting substructures [c1cc(Cl)cc1, CNCO, CCOP(=S)(OCC)O] influencing the chronic aquatic toxicity. These results showed that the compounds with low molecular volume, high polarity and molecular weight could contribute to adverse effects on rotifers, facilitating the deeper understanding of chronic toxicity mechanisms. In addition, forecast models had better performances than the common models embedded into ECOSAR software. This study provided insights into structural features responsible for the toxicity of different groups of chemicals and thereby allowed for the rational design of green and safer alternatives.

Metabolomic analysis predicted changes in growth rate in Daphnia magna exposed to acetaminophen

As urbanization and the global population increases, pollutants associated with municipal wastewater such as pharmaceuticals are becoming more prevalent in aquatic environments. Acetaminophen (paracetamol) is a widely used drug worldwide and one of the most frequently detected pharmaceuticals in freshwater ecosystems. This study investigated the impact of acetaminophen on the metabolite profile of Daphnia magna at two life stages; and used these metabolomic findings to hypothesize a potential impact at a higher organismal level which was subsequently tested experimentally. Targeted polar metabolite analysis by liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to measure changes in the concentration of 51 metabolites in the neonate (> 24 h old) and adult (8 day-old) daphnids following a 48-h exposure to sub-lethal concentrations of acetaminophen. The impact of acetaminophen on the metabolic profile of neonates was widely different from adults. Also, acetaminophen exposure perturbed the abundance of nucleotides more extensively than other metabolites. The acute metabolomic experimental results led to the hypotheses that exposure to sub-lethal concentrations of acetaminophen upregulates protein synthesis in D. magna and subsequently increases growth during early life stages and has an opposite impact on adults. Accordingly, a 10 day growth rate experiment indicated that exposure to acetaminophen elevated biomass production in neonates but not in adults. These novel findings demonstrate that a targeted analysis and interpretation of the changes in the polar metabolic profile of organisms in response to environmental stressors could be used as a tool to predict changes at higher biological levels. As such, this study further emphasizes the incorporation of molecular-level platforms as critical and robust tools in environmental assessment frameworks and biomonitoring of aquatic ecosystems.

Integration of Omics Approaches Enhances the Impact of Scientific Research in Environmental Applications

In the original article [...]

Comparative genome analysis of the monogonont marine rotifer Brachionus manjavacas Australian strain: Potential application for ecotoxicology and environmental genomics

This is the first study to analyze the whole-genome sequence of B. manjavacas Australian (Aus.) strain through combination of Oxford Nanopore long-read seq, resulting in a total length of 108.1 Mb and 75 contigs. Genome-wide detoxification related gene families in B. manjavacas Aus. strain were comparatively analyzed with those previously identified in other Brachionus spp., including B. manjavacas German (Ger.) strain. Most of the subfamilies in detoxification related families (CYPs, GSTs, and ABCs) were highly conserved and confirmed orthologous relationship with Brachionus spp., and with accumulation of genome data, clear differences between genomic repertoires were demonstrated the marine and the freshwater species. Furthermore, strain-specific genetic variations were present between the Aus. and Ger. strains of B. manjavacas. This whole-genome analysis provides in-depth review on the genomic structural differences for detoxification-related gene families and further provides useful information for comparative ecotoxicological studies and evolution of detoxification mechanisms in Brachionus spp.

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.

From single-omics to interactomics: How can ligand-induced perturbations modulate single-cell phenotypes?

Cells suffer from perturbations by different stimuli, which, consequently, rise to individual alterations in their profile and function that may end up affecting the tissue as a whole. This is no different if we consider the effect of a therapeutic agent on a biological system. As cells are exposed to external ligands their profile can change at different single-omics levels. Detecting how these changes take place through different sequencing technologies is key to a better understanding of the effects of therapeutic agents. Single-cell RNA-sequencing stands out as one of the most common approaches for cell profiling and perturbation analysis. As a result, single-cell transcriptomics data can be integrated with other omics data sources, such as proteomics and epigenomics data, to clarify the perturbation effects and mechanism at the cell level. Appropriate computational tools are key to process and integrate the available information. This chapter focuses on the recent advances on ligand-induced perturbation and single-cell omics computational tools and algorithms, their current limitations, and how the deluge of data can be used to improve the current process of drug research and development.

The genome of the freshwater monogonont rotifer Brachionus rubens: Identification of phase I, II, and III detoxification genes

Monogonont rotifers are common species in aquatic environments and make model species for ecotoxicology studies. Whole genomes of several species of the genus Brachionus have been assembled, but no information on the freshwater rotifer Brachionus rubens has been reported. In this study, the whole-genome sequence of B. rubens was successfully assembled using NextDenovo. The total length of the genome was 132.7 Mb (N50 = 2.51 Mb), including 122 contigs. The GC contents accounted for 29.96% of the genome. Aquatic organisms are always exposed to various external stresses, and a comprehensive genomic analysis is needed to better understand the adverse effects on organisms. This paper focuses on the ecotoxicological aspect and conducted genome analysis of representative gene families involved in detoxification mechanisms against environmental stressors. Specifically, we identified cytochrome P450 genes (CYPs) of phase I, glutathione S-transferase genes (GSTs) of phase II, and ATP-binding cassette transporter genes (ABCs) of phase III in the genome of B. rubens. Gene duplications were found in CYP, GST, and ABC genes, as is the case for other Brachionus rotifers. Our results suggest that these detoxification-related gene families have evolved in a species-specific and/or lineage-specific manner. This paper improves our understanding of how the freshwater Brachionus rotifers respond to environmental stressors in a molecular ecotoxicology context.

Evaluation of the influences of low dose polybrominated diphenyl ethers exposure on human early retinal development

Increasing evidence in animal models has suggested that polybrominated diphenyl ethers (PBDEs), a class of brominated flame retardants, can cause retinotoxicity. However, data on the influence of PBDE treatment on human retinal development are scarce due to the lack of appropriate models. In the present study, we report the utilization of human embryonic stem cell-derived retinal organoids (hESC-ROs) for toxicity assessment of the most common PBDE congener (BDE-47) during the early stages of retinal development. Exposure to BDE-47 decreased the thickness and area of the neural retina (NR) of hESC-ROs in a dose- and time-dependent manner. Abnormal retinal cell distributions, disordered NR structures, and neural rosette-like structures were found on hESC-ROs after low-level BDE-47 exposure. Moreover, BDE-47 exposure decreased cell proliferation, promoted cell apoptosis, and caused abnormal differentiation. Transcriptomic analysis demonstrated that differentially expressed genes, caused by BDE-47, were enriched in extracellular matrix organization. Metabolomic studies of hESC-ROs revealed significant changes in the metabolism of purine and glutathione after BDE-47 exposure for five weeks. This study clarifies the retinotoxicity of low-level BDE-47 treatment and highlights the powerfulness of the hESC-RO model, deepening our understanding of BDE-47-driven human early retina developmental toxicity.

Metabolic response of earthworms (Pheretima guillemi) to silver nanoparticles in sludge-amended soil

Silver nanoparticles (AgNPs) can enter soils via the application of sludge and pose risks to soil invertebrates. However, current knowledge regarding the toxicity of AgNPs at environmentally relevant concentration is insufficient, especially at the molecular level. Therefore, we examined the effects of low-level AgNPs (7.2 mg kg^-1, dry weight) on the bioaccumulation, pathology and metabolism of earthworms (Pheretima guillemi). After exposure for 28 d, earthworms were dissected into digestive system and the rest of the body to explore the response of different body parts to AgNPs. Ag concentration in the digestive system of exposed group (2.5 mg kg^-1, dry weight) was significantly higher than that of the control group (0.5 mg kg^-1, dry weight). AgNPs exposure had no significant effects on the survival and growth, but induced intestinal damage and metabolic interference to earthworms relative to the control. Metabolomics analysis showed that AgNPs exposure disturbed the glycerophospholipid metabolism, glutathione metabolism and energy metabolism in the digestive system and the energy metabolism in the rest of the body. AgNPs exposure also induced lipid peroxidation in the digestive system. The different metabolic responses between two body parts highlighted the importance of the uptake routes of Ag. These results provide a biochemical insight for the risk assessment of low-level AgNPs in terrestrial environment.

Acute temperature adaptation mechanisms in the native reptile species Eremias argus

Reptiles are sensitive to temperature changes as ectotherm animals. The climate warming may pose more serious threat to reptiles. Although the behavior effect and reproduction biology have been well studied, little information is available about the adaptation mechanisms of reptiles to temperature stress. In this study, the native Chinese species, Eremias argus were incubated at 15 (cold stress), 25 (control group) and 35 °C (thermal stress) for 24 h. The transcriptome and metabolome technology were applied to investigate the molecular regulation mechanisms of lizards to acute temperature changes. The CIRBP and HSPA8 were hub genes in response to temperature adaptation. The increased expression of PER gene in lizard circadian rhythm is associated with tyrosine metabolism after cold or thermal stress. The poly-unsaturated fatty acids in female lizard liver were significantly increased with up-regulation of FASN and ACACA genes after thermal stress, which proved the disruption of fatty acid biosynthesis pathway in corporation with the altered body weight. The cortisol and testosterone were important steroid hormones in response to temperature changes especially in male lizard liver. The increased CIRBP gene expression in lizard gonads suppressed the KDM6B gene, which regulates the testis development and may induce sex reversal in male lizard after thermal stress. The adaptation responses of lizards to temperature stress may threaten the health status of wild population.

Overlooked environmental risks deriving from aqueous transformation of bisphenol alternatives: Integration of chemical and toxicological insights

Owing to the widespread prevalence and ecotoxicity of bisphenol alternatives such as bisphenol S, bisphenol F, and bisphenol AF, the past decade has witnessed the publication of a remarkable number of studies related to their transformation and remediation in natural waters. However, the reactivity, removal efficiency, transformation products (TPs), and mechanisms of such emerging pollutants by different treatment processes have not been well elucidated. Particularly, the transformation-driven environmental risks have been mostly overlooked. Therefore, we present a review to address these issues from chemical and toxicological viewpoints. Four degradation systems can be largely classified as catalytic persulfate (PS) oxidation, non-catalytic oxidation, photolysis and photocatalysis, and biodegradation. It was found that bisphenol alternatives possess distinct reactivities with different oxidizing species, with the highest performance for hydroxyl radicals. All systems exhibit superior elimination efficiency for these compounds. The inadequate mineralization suggests the formation of recalcitrant TPs, from which the overall reaction pathways are proposed. The combined experimental and in silico analysis indicates that many TPs have developmental toxicity, endocrine-disrupting effects, and genotoxicity. Notably, catalytic PS systems and non-catalytic oxidation result in the formation of coupling products as well as halogenated TPs with higher acute and chronic toxicity and lower biodegradability than the parent compounds. In contrast, photolysis and photocatalysis generate hydroxylated and bond-cleavage TPs with less toxicity. Overall, this review highlights the secondary environmental risks from the transformation of bisphenol alternatives by conventional and emerging treatment processes. Finally, future perspectives are recommended to address the knowledge gaps of these contaminants in aquatic ecosystems.

Lanthanum and cerium disrupt similar biological pathways and interact synergistically in Triticum aestivum as revealed by metabolomic profiling and quantitative modeling

The industrial and agricultural applications of rare earth elements (REEs) lead to considerable REE emissions into environment. Yet, little is known about the molecular-level effects and interactions of REEs in terrestrial plants. Herein, the individual and joint effects of La and Ce in Triticum aestivum were investigated using mass spectrometry-based metabolomics. Metabolic effect level index (MELI) was utilized as a readable endpoint for quantifying mixture interactions. Exposure to single La/Ce at environmentally relevant levels induced significant dose-dependent metabolic changes. The highly overlap of differential metabolites and perturbed pathways of La and Ce suggested their similar mode of action. Exposure to La-Ce mixtures did not induce additional metabolic pathway perturbation. Specifically, metabolism of amino sugar and nucleotide sugar, starch and sucrose, fructose and mannose, glycerophospholipid and purine were disrupted for both single and binary exposures. These results, together with physiological indicators, point to REE-induced oxidative stress, energy expenditure, DNA damage and membrane disturbance. The MELI calculations showed that La and Ce interacted synergistically at the overall metabolic level, which could be causally linked to synergistic interaction at the individual level (root elongation). This work proved metabolomics could be an important and effective strategy for interpreting toxicity and interactions of REE mixtures.

The aryl hydrocarbon receptor: A predominant mediator for the toxicity of emerging dioxin-like compounds

The aryl hydrocarbon receptor (AHR) is a member of the basic helix-loop-helix/Per-ARNT-Sim (bHLH-PAS) family of transcription factors and has broad biological functions. Early after the identification of the AHR, most studies focused on its roles in regulating the expression of drug-metabolizing enzymes and mediating the toxicity of dioxins and dioxin-like compounds (DLCs). Currently, more diverse functions of AHR have been identified, indicating that AHR is not just a dioxin receptor. Dioxins and DLCs occur ubiquitously and have diverse health/ecological risks. Additional research is required to identify both shared and compound-specific mechanisms, especially for emerging DLCs such as polyhalogenated carbazoles (PHCZs), polychlorinated diphenyl sulfides (PCDPSs), and others, of which only a few investigations have been performed at present. Many of the toxic effects of emerging DLCs were observed to be predominantly mediated by the AHR because of their structural similarity as dioxins, and the in vitro TCDD-relative potencies of certain emerging DLC congeners are comparable to or even greater than the WHO-TEFs of OctaCDD, OctaCDF, and most coplanar PCBs. Due to the close relationship between AHR biology and environmental science, this review begins by providing novel insights into AHR signaling (canonical and non-canonical), AHR's biochemical properties (AHR structure, AHR-ligand interaction, AHR-DNA binding), and the variations during AHR transactivation. Then, AHR ligand classification and the corresponding mechanisms are discussed, especially the shared and compound-specific, AHR-mediated effects and mechanisms of emerging DLCs. Accordingly, a series of in vivo and in vitro toxicity evaluation methods based on the AHR signaling pathway are reviewed. In light of current advances, future research on traditional and emerging DLCs will enhance our understanding of their mechanisms, toxicity, potency, and ecological impacts.