Transcriptional Profiling of Dibenzo[def,p]chrysene-induced Spleen Atrophy Provides Mechanistic Insights into its Immunotoxicity in MutaMouse

Potential Toxicity Risk Assessment and Priority Control Strategy for PAHs Metabolism and Transformation Behaviors in the Environment

In this study, 16 PAHs were selected as the priority control pollutants to summarize their environmental metabolism and transformation processes, including photolysis, plant degradation, bacterial degradation, fungal degradation, microalgae degradation, and human metabolic transformation. Meanwhile, a total of 473 PAHs by-products generated during their transformation and degradation in different environmental media were considered. Then, a comprehensive system was established for evaluating the PAHs by-products' neurotoxicity, immunotoxicity, phytotoxicity, developmental toxicity, genotoxicity, carcinogenicity, and endocrine-disrupting effect through molecular docking, molecular dynamics simulation, 3D-QSAR model, TOPKAT method, and VEGA platform. Finally, the potential environmental risk (phytotoxicity) and human health risks (neurotoxicity, immunotoxicity, genotoxicity, carcinogenicity, developmental toxicity, and endocrine-disrupting toxicity) during PAHs metabolism and transformation were comprehensively evaluated. Among the 473 PAH's metabolized and transformed products, all PAHs by-products excluding ACY, CHR, and DahA had higher neurotoxicity, 152 PAHs by-products had higher immunotoxicity, and 222 PAHs by-products had higher phytotoxicity than their precursors during biological metabolism and environmental transformation. Based on the TOPKAT model, 152 PAH by-products possessed potential developmental toxicity, and 138 PAH by-products had higher genotoxicity than their precursors. VEGA predicted that 247 kinds of PAH derivatives had carcinogenic activity, and only the natural transformation products of ACY did not have carcinogenicity. In addition to ACY, 15 PAHs produced 123 endocrine-disrupting substances during metabolism and transformation. Finally, the potential environmental and human health risks of PAHs metabolism and transformation products were evaluated using metabolic and transformation pathway probability and degree of toxic risk as indicators. Accordingly, the priority control strategy for PAHs was constructed based on the risk entropy method by screening the priority control pathways. This paper assesses the potential human health and environmental risks of PAHs in different environmental media with the help of models and toxicological modules for the toxicity prediction of PAHs by-products, and thus designs a risk priority control evaluation system for PAHs.

Soil exposed to silver nanoparticles reveals significant changes in community structure and altered microbial transcriptional profiles

Anthropogenic activities can disrupt soil ecosystems, normally resulting in reduced soil microbial health. Regulatory agencies need to determine the effects of uncharacterized substances on soil microbial health to establish the safety of these chemicals if they end up in the environment. Previous work has focused on measuring traditional ecotoxicologial endpoints within the categories of microbial biomass, activity, and community structure/diversity. Because these tests can be labor intensive, lengthy to conduct, and cannot measure changes in individual gene functions, we wanted to establish whether metatranscriptomics could be used as a more sensitive endpoint and provide a perspective on community function that is more informative than taxonomic identification of microbes alone. We spiked a freshly collected sandy loam soil (Vulcan, Alberta, Canada) with 0, 60, 145, 347, 833, and 2000 mg kg^-1 of silver nanoparticles (AgNPs), a known antagonist of microorganisms due to its propensity for dissolution of toxic silver ions. Assessments performed in our previous work using traditional tests demonstrated the toxicity of AgNPs on soil microbial processes. We expanded this analysis with genomics-based tests by measuring changes in community taxonomic structure and function using 16S rDNA profiling and metatranscriptomics. In addition to identifying bacterial taxa affected by AgNPs, we found that genes involved in heavy metal resistance (e.g., the CzcA efflux pump) and other toxicity response pathways were highly upregulated in the presence of silver. Dose-response analysis using BMDExpress2 software successfully modeled many physiologically relevant genes responding to low concentrations of AgNPs. We found that the transcriptomic point of departure (BMD₅₀) was lower than the IC50s calculated using the traditional tests in our previous work. These results suggest that dose-response modeling of metatranscriptomic gene expression is a useful tool in soil microbial health assessment. SUMMARY: Genomics-based endpoints for the assessment of soil microbial health can be used to perform quantitative dose-response modeling, and soil-based RNAseq adds functional insights.

Whole Genome Sequencing of the Mutamouse Model Reveals Strain- and Colony-Level Variation, and Genomic Features of the Transgene Integration Site

The MutaMouse transgenic rodent model is widely used for assessing in vivo mutagenicity. Here, we report the characterization of MutaMouse's whole genome sequence and its genetic variants compared to the C57BL/6 reference genome. High coverage (>50X) next-generation sequencing (NGS) of whole genomes from multiple MutaMouse animals from the Health Canada (HC) colony showed ~5 million SNVs per genome, ~20% of which are putatively novel. Sequencing of two animals from a geographically separated colony at Covance indicated that, over the course of 23 years, each colony accumulated 47,847 (HC) and 17,677 (Covance) non-parental homozygous single nucleotide variants. We found no novel nonsense or missense mutations that impair the MutaMouse response to genotoxic agents. Pairing sequencing data with array comparative genomic hybridization (aCGH) improved the accuracy and resolution of copy number variants (CNVs) calls and identified 300 genomic regions with CNVs. We also used long-read sequence technology (PacBio) to show that the transgene integration site involved a large deletion event with multiple inversions and rearrangements near a retrotransposon. The MutaMouse genome gives important genetic context to studies using this model, offers insight on the mechanisms of structural variant formation, and contributes a framework to analyze aCGH results alongside NGS data.

Comparative mechanisms of PAH toxicity by benzo[a]pyrene and dibenzo[def,p]chrysene in primary human bronchial epithelial cells cultured at air-liquid interface

Current assumption for assessing carcinogenic risk of polycyclic aromatic hydrocarbons (PAHs) is that they function through a common mechanism of action; however, recent studies demonstrate that PAHs can act through unique mechanisms potentially contributing to cancer outcomes in a non-additive manner. Using a primary human 3D bronchial epithelial culture (HBEC) model, we assessed potential differences in mechanism of toxicity for two PAHs, benzo[a]pyrene (BAP) and dibenzo[def,p]chrysene (DBC), compared to a complex PAH mixture based on short-term biosignatures identified from transcriptional profiling. Differentiated bronchial epithelial cells were treated with BAP (100-500 μg/ml), DBC (10 μg/ml), and coal tar extract (CTE 500-1500 μg/ml, SRM1597a) for 48 h and gene expression was measured by RNA sequencing or quantitative PCR. Comparison of BAP and DBC gene signatures showed that the majority of genes (~60%) were uniquely regulated by treatment, including signaling pathways for inflammation and DNA damage by DBC and processes for cell cycle, hypoxia and oxidative stress by BAP. Specifically, BAP upregulated targets of AhR, NRF2, and KLF4, while DBC downregulated these same targets, suggesting a chemical-specific pattern in transcriptional regulation involved in antioxidant response, potentially contributing to differences in PAH potency. Other processes were regulated in common by all PAH treatments, BAP, DBC and CTE, including downregulation of genes involved in cell adhesion and reduced functional measurements of barrier integrity. This work supports prior in vivo studies and demonstrates the utility of profiling short-term biosignatures in an organotypic 3D model to identify mechanisms linked to carcinogenic risk of PAHs in humans.

Hydroxyfasudil alleviates demyelination through the inhibition of MOG antibody and microglia activation in cuprizone mouse model

Multiple sclerosis (MS) is an immune-mediated demyelinating disease of the central nervous system characterized by oligodendrocyte loss and progressive neurodegeneration. The cuprizone (CPZ)-induced demyelination is widely used to investigate the demyelination/remyelination. Here, we explored the therapeutic effects of Hydroxyfasudil (HF), an active metabolite of Fasudil, in CPZ model. HF improved behavioral abnormality and reduced myelin damage in the corpus callosum. Splenic atrophy and myelin oligodendrocyte glycoprotein (MOG) antibody were observed in CPZ model, which were partially restored and obviously inhibited by HF, therefore reducing pathogenic binding of MOG antibody to oligodendrocytes. HF inhibited the percentages of CD4⁺IL-17⁺ T cells from splenocytes and infiltration of CD4⁺ T cells and CD68⁺ macrophages in the brain. HF also declined microglia-mediated neuroinflammation, and promoted the production of astrocyte-derived brain derived neurotrophic factor (BDNF) and regeneration of NG2⁺ oligodendrocyte precursor cells. These results provide potent evidence for the therapeutic effects of HF in CPZ-induced demyelination.

Integration of the TGx-28.65 genomic biomarker with the flow cytometry micronucleus test to assess the genotoxicity of disperse orange and 1,2,4-benzenetriol in human TK6 cells

In vitro gene expression signatures to predict toxicological responses can provide mechanistic context for regulatory testing. We previously developed the TGx-28.65 genomic biomarker from a database of gene expression profiles derived from human TK6 cells exposed to 28 well-known compounds. The biomarker comprises 65 genes that can classify chemicals as DNA damaging or non-DNA damaging. In this study, we applied the TGx-28.65 genomic biomarker in parallel with the in vitro micronucleus (MN) assay to determine if two chemicals of regulatory interest at Health Canada, disperse orange (DO: the orange azo dye 3-[[4-[(4-Nitrophenyl)azo]phenyl] benzylamino]propanenitrile) and 1,2,4-benzenetriol (BT: a metabolite of benzene) are genotoxic or non-genotoxic. Both chemicals caused dose-dependent declines in relative survival and increases in apoptosis. A strong significant increase in MN induction was observed for all concentrations of BT; the top two concentrations of DO also caused a statistically significant increase in MN, but these increases were <2-fold above controls. TGx-28.65 analysis classified BT as genotoxic at all three concentrations and DO as genotoxic at the mid and high concentrations. Thus, although DO only caused a small increase in MN, this response was sufficient to induce a cellular DNA damage response. Benchmark dose modeling confirmed that BT is much more potent than DO. The results strongly suggest that follow-up work is required to assess whether DO and BT are also genotoxic in vivo. This is particularly important for DO, which may require metabolic activation by bacterial gut flora to fully induce its genotoxic potential. Our previously published data and this proof of concept study suggest that the TGx-28.65 genomic biomarker has the potential to add significant value to existing approaches used to assess genotoxicity.

Recommended approaches in the application of toxicogenomics to derive points of departure for chemical risk assessment

There is increasing interest in the use of quantitative transcriptomic data to determine benchmark dose (BMD) and estimate a point of departure (POD) for human health risk assessment. Although studies have shown that transcriptional PODs correlate with those derived from apical endpoint changes, there is no consensus on the process used to derive a transcriptional POD. Specifically, the subsets of informative genes that produce BMDs that best approximate the doses at which adverse apical effects occur have not been defined. To determine the best way to select predictive groups of genes, we used published microarray data from dose–response studies on six chemicals in rats exposed orally for 5, 14, 28, and 90 days. We evaluated eight approaches for selecting genes for POD derivation and three previously proposed approaches (the lowest pathway BMD, and the mean and median BMD of all genes). The relationship between transcriptional BMDs derived using these 11 approaches and PODs derived from apical data that might be used in chemical risk assessment was examined. Transcriptional BMD values for all 11 approaches were remarkably aligned with corresponding apical PODs, with the vast majority of toxicogenomics PODs being within tenfold of those derived from apical endpoints. We identified at least four approaches that produce BMDs that are effective estimates of apical PODs across multiple sampling time points. Our results support that a variety of approaches can be used to derive reproducible transcriptional PODs that are consistent with PODs produced from traditional methods for chemical risk assessment.

Iron mediated toxicity and programmed cell death: A review and a re-examination of existing paradigms

Iron is an essential micronutrient that is problematic for biological systems since it is toxic as it generates free radicals by interconverting between ferrous (Fe²⁺) and ferric (Fe³⁺) forms. Additionally, even though iron is abundant, it is largely insoluble so cells must treat biologically available iron as a valuable commodity. Thus elaborate mechanisms have evolved to absorb, re-cycle and store iron while minimizing toxicity. Focusing on rarely encountered situations, most of the existing literature suggests that iron toxicity is common. A more nuanced examination clearly demonstrates that existing regulatory processes are more than adequate to limit the toxicity of iron even in response to iron overload. Only under pathological or artificially harsh situations of exposure to excess iron does it become problematic. Here we review iron metabolism and its toxicity as well as the literature demonstrating that intracellular iron is not toxic but a stress responsive programmed cell death-inducing second messenger.

A framework for the use of single-chemical transcriptomics data in predicting the hazards associated with complex mixtures of polycyclic aromatic hydrocarbons

The assumption of additivity applied in the risk assessment of environmental mixtures containing carcinogenic polycyclic aromatic hydrocarbons (PAHs) was investigated using transcriptomics. MutaTMMouse were gavaged for 28 days with three doses of eight individual PAHs, two defined mixtures of PAHs, or coal tar, an environmentally ubiquitous complex mixture of PAHs. Microarrays were used to identify differentially expressed genes (DEGs) in lung tissue collected 3 days post-exposure. Cancer-related pathways perturbed by the individual or mixtures of PAHs were identified, and dose–response modeling of the DEGs was conducted to calculate gene/pathway benchmark doses (BMDs). Individual PAH-induced pathway perturbations (the median gene expression changes for all genes in a pathway relative to controls) and pathway BMDs were applied to models of additivity [i.e., concentration addition (CA), generalized concentration addition (GCA), and independent action (IA)] to generate predicted pathway-specific dose–response curves for each PAH mixture. The predicted and observed pathway dose–response curves were compared to assess the sensitivity of different additivity models. Transcriptomics-based additivity calculation showed that IA accurately predicted the pathway perturbations induced by all mixtures of PAHs. CA did not support the additivity assumption for the defined mixtures; however, GCA improved the CA predictions. Moreover, pathway BMDs derived for coal tar were comparable to BMDs derived from previously published coal tar-induced mouse lung tumor incidence data. These results suggest that in the absence of tumor incidence data, individual chemical-induced transcriptomics changes associated with cancer can be used to investigate the assumption of additivity and to predict the carcinogenic potential of a mixture.

Tissue-specific in vivo genetic toxicity of nine polycyclic aromatic hydrocarbons assessed using the Muta™Mouse transgenic rodent assay

Test batteries to screen chemicals for mutagenic hazard include several endpoints regarded as effective for detecting genotoxic carcinogens. Traditional in vivo methods primarily examine clastogenic endpoints in haematopoietic tissues. Although this approach is effective for identifying systemically distributed clastogens, some mutagens may not induce clastogenic effects; moreover, genotoxic effects may be restricted to the site of contact and/or related tissues. An OECD test guideline for transgenic rodent (TGR) gene mutation assays was released in 2011, and the TGR assays permit assessment of mutagenicity in any tissue. This study assessed the responses of two genotoxicity endpoints following sub-chronic oral exposures of male Muta™Mouse to 9 carcinogenic polycyclic aromatic hydrocarbons (PAHs). Clastogenicity was assessed via induction of micronuclei in peripheral blood, and mutagenicity via induction of lacZ transgene mutations in bone marrow, glandular stomach, small intestine, liver, and lung. Additionally, the presence of bulky PAH-DNA adducts was examined. Five of the 9 PAHs elicited positive results across all endpoints in at least one tissue, and no PAHs were negative or equivocal across all endpoints. All PAHs were positive for lacZ mutations in at least one tissue (sensitivity=100%), and for 8 PAHs, one or more initial sites of chemical contact (i.e., glandular stomach, liver, small intestine) yielded a greater response than bone marrow. Five PAHs were positive in the micronucleus assay (sensitivity=56%). Furthermore, all PAHs produced DNA adducts in at least one tissue. The results demonstrate the utility of the TGR assay for mutagenicity assessment, especially for compounds that may not be systemically distributed.