3' RNA-seq is superior to standard RNA-seq in cases of sparse data but inferior at identifying toxicity pathways in a model organism

Introduction: The application of RNA-sequencing has led to numerous breakthroughs related to investigating gene expression levels in complex biological systems. Among these are knowledge of how organisms, such as the vertebrate model organism zebrafish (Danio rerio), respond to toxicant exposure. Recently, the development of 3' RNA-seq has allowed for the determination of gene expression levels with a fraction of the required reads compared to standard RNA-seq. While 3' RNA-seq has many advantages, a comparison to standard RNA-seq has not been performed in the context of whole organism toxicity and sparse data. Methods and results: Here, we examined samples from zebrafish exposed to perfluorobutane sulfonamide (FBSA) with either 3' or standard RNA-seq to determine the advantages of each with regards to the identification of functionally enriched pathways. We found that 3' and standard RNA-seq showed specific advantages when focusing on annotated or unannotated regions of the genome. We also found that standard RNA-seq identified more differentially expressed genes (DEGs), but that this advantage disappeared under conditions of sparse data. We also found that standard RNA-seq had a significant advantage in identifying functionally enriched pathways via analysis of DEG lists but that this advantage was minimal when identifying pathways via gene set enrichment analysis of all genes. Conclusions: These results show that each approach has experimental conditions where they may be advantageous. Our observations can help guide others in the choice of 3' RNA-seq vs standard RNA sequencing to query gene expression levels in a range of biological systems.

Do Carbon Nanotubes and Asbestos Fibers Exhibit Common Toxicity Mechanisms?

During the last two decades several nanoscale materials were engineered for industrial and medical applications. Among them carbon nanotubes (CNTs) are the most exploited nanomaterials with global production of around 1000 tons/year. Besides several commercial benefits of CNTs, the fiber-like structures and their bio-persistency in lung tissues raise serious concerns about the possible adverse human health effects resembling those of asbestos fibers. In this review, we present a comparative analysis between CNTs and asbestos fibers using the following four parameters: (1) fibrous needle-like shape, (2) bio-persistent nature, (3) high surface to volume ratio and (4) capacity to adsorb toxicants/pollutants on the surface. We also compare mechanisms underlying the toxicity caused by certain diameters and lengths of CNTs and asbestos fibers using downstream pathways associated with altered gene expression data from both asbestos and CNT exposure. Our results suggest that indeed certain types of CNTs are emulating asbestos fiber as far as associated toxicity is concerned.

Toxicometabolomics: Small Molecules to Answer Big Toxicological Questions

Given the high biological impact of classical and emerging toxicants, a sensitive and comprehensive assessment of the hazards and risks of these substances to organisms is urgently needed. In this sense, toxicometabolomics emerged as a new and growing field in life sciences, which use metabolomics to provide new sets of susceptibility, exposure, and/or effects biomarkers; and to characterize in detail the metabolic responses and altered biological pathways that various stressful stimuli cause in many organisms. The present review focuses on the analytical platforms and the typical workflow employed in toxicometabolomic studies, and gives an overview of recent exploratory research that applied metabolomics in various areas of toxicology.

2018 Lush Science Prize

The Lush Prize supports animal-free testing by awarding money prizes of up to £350,000 per year to the most effective projects and individuals who have been working towards the goal of replacing animals in product or ingredient safety testing. Since its inception in 2012, the Lush Prize has distributed almost £2 million. Prizes are awarded for developments in five strategic areas: Science; Lobbying; Training; Public Awareness; and Young Researchers. In 2015, the judges also awarded a Black Box prize for the development of the skin sensitisation Adverse Outcome Pathway and its associated in vitro assays. The Science Prize is awarded to researchers whose work the judging panel believe to have made the most significant contribution, in the preceding year, to the replacement of animal testing. This 2018 Science Background paper outlines the research projects that were presented to the Prize judges as potential candidates for the 2018 Lush Science Prize award. To obtain an overview of developments in the field of animal replacement in toxicity research, recent work by the relevant scientific institutions and projects in this area, including the OECD, CAAT, ECVAM, UK NC3Rs, US Tox21 Programme, the ToxCast programme and EU-ToxRisk, was reviewed. Recent developments in toxicity testing research were investigated by searching the relevant literature. Abstracts from conferences focusing on animal replacement in toxicity testing that were held in the preceding 12 months, were also analysed, including those from the 2017 10th World Congress on Alternatives and Animals in the Life Sciences and the 2018 Society of Toxicology annual conference.

Corrigendum: Investigation of Nrf2, AhR and ATF4 Activation in Toxicogenomic Databases

[This corrects the article DOI: 10.3389/fgene.2018.00429.].

Investigation of Nrf2, AhR and ATF4 Activation in Toxicogenomic Databases

Toxicological responses to chemical insult are largely regulated by transcriptionally activated pathways that may be independent, correlated and partially or fully overlapping. Investigating the dynamics of the interactions between stress responsive transcription factors from toxicogenomic data and defining the signature of each of them is an additional step toward a system level understanding of perturbation driven mechanisms. To this end, we investigated the segregation of the genes belonging to the three following transcriptionally regulated pathways: the AhR pathway, the Nrf2 pathway and the ATF4 pathway. Toxicogenomic datasets from three projects (carcinoGENOMICS, Predict-IV and TG-GATEs) obtained in various experimental conditions (in human and rat in vitro liver and kidney models and rat in vivo, with bolus administration and with repeated doses) were combined and consolidated where overlaps between datasets existed. A bioinformatic analysis was performed to refine pathways' signatures and to create chemical activation capacity scores to classify chemicals by their potency and selectivity of activation of each pathway. With some refinement such an approach may improve chemical safety classification and allow biological read across on a pathway level.

Dose and Time Dependencies in Stress Pathway Responses during Chemical Exposure: Novel Insights from Gene Regulatory Networks

Perturbation of biological networks is often observed during exposure to xenobiotics, and the identification of disturbed processes, their dynamic traits, and dose–response relationships are some of the current challenges for elucidating the mechanisms determining adverse outcomes. In this scenario, reverse engineering of gene regulatory networks (GRNs) from expression data may provide a system-level snapshot embedded within accurate molecular events. Here, we investigate the composition of GRNs inferred from groups of chemicals with two distinct outcomes, namely carcinogenicity [azathioprine (AZA) and cyclophosphamide (CYC)] and drug-induced liver injury (DILI; diclofenac, nitrofurantoin, and propylthiouracil), and a non-carcinogenic/non-DILI group (aspirin, diazepam, and omeprazole). For this, we analyzed publicly available exposed in vitro human data, taking into account dose and time dependencies. Dose–Time Network Identification (DTNI) was applied to gene sets from exposed primary human hepatocytes using four stress pathways, namely endoplasmic reticulum (ER), NF-κB, NRF2, and TP53. Inferred GRNs suggested case specificity, varying in interactions, starting nodes, and target genes across groups. DILI and carcinogenic compounds were shown to directly affect all pathway-based GRNs, while non-DILI/non-carcinogenic chemicals only affected NF-κB. NF-κB-based GRNs clearly illustrated group-specific disturbances, with the cancer-related casein kinase CSNK2A1 being a target gene only in the carcinogenic group, and opposite regulation of NF-κB subunits being observed in DILI and non-DILI/non-carcinogenic groups. Target genes in NRF2-based GRNs shared by DILI and carcinogenic compounds suggested markers of hepatotoxicity. Finally, we indicate several of these group-specific interactions as potentially novel. In summary, our reversed-engineered GRNs are capable of revealing dose dependent, chemical-specific mechanisms of action in stress-related biological networks.

Chemical-agnostic hazard prediction: statistical inference of in vitro toxicity pathways from proteomics responses to chemical mixtures

Toxicity pathways have been defined as normal cellular pathways that, when sufficiently perturbed as a consequence of chemical exposure, lead to an adverse outcome. If an exposure alters one or more normal biological pathways to an extent that leads to an adverse toxicity outcome, a significant correlation must exist between the exposure, the extent of pathway alteration, and the degree of adverse outcome. Biological pathways are regulated at multiple levels, including transcriptional, post-transcriptional, post-translational, and targeted degradation, each of which can affect the levels and extents of modification of proteins involved in the pathways. Significant alterations of toxicity pathways resulting from changes in regulation at any of these levels therefore are likely to be detectable as alterations in the proteome. We hypothesize that significant correlations between exposures, adverse outcomes, and changes in the proteome have the potential to identify putative toxicity pathways, facilitating selection of candidate targets for high throughput screening, even in the absence of a priori knowledge of either the specific pathways involved or the specific agents inducing the pathway alterations. We explored this hypothesis in vitro in BEAS-2B human airway epithelial cells exposed to different concentrations of Ni²⁺, Cd²⁺, and Cr⁶⁺, alone and in defined mixtures. Levels and phosphorylation status of a variety of signaling pathway proteins and cytokines were measured after 48 hours exposure, together with cytotoxicity. Least Absolute Shrinkage and Selection Operator (LASSO) multiple regression was used to identify a subset of these proteins that constitute a putative toxicity pathway capable of predicting cytotoxicity. The putative toxicity pathway for cytotoxicity of these metals and metal mixtures identified by LASSO is composed of phospho-RPS6KB1, phospho-p53, cleaved CASP3, phospho-MAPK8, IL-10, and Hif-1α. As this approach does not depend on knowledge of the chemical composition of the mixtures, it may be generally useful for identifying sets of proteins predictive of adverse effects for a variety of mixtures, including complex environmental mixtures of unknown composition.

Profiling dose-dependent activation of p53-mediated signaling pathways by chemicals with distinct mechanisms of DNA damage

As part of a larger effort to provide proof-of-concept in vitro-only risk assessments, we have developed a suite of high-throughput assays for key readouts in the p53 DNA damage response toxicity pathway: double-strand break DNA damage (p-H2AX), permanent chromosomal damage (micronuclei), p53 activation, p53 transcriptional activity, and cell fate (cell cycle arrest, apoptosis, micronuclei). Dose-response studies were performed with these protein and cell fate assays, together with whole genome transcriptomics, for three prototype chemicals: etoposide, quercetin, and methyl methanesulfonate. Data were collected in a human cell line expressing wild-type p53 (HT1080) and results were confirmed in a second p53 competent cell line (HCT 116). At chemical concentrations causing similar increases in p53 protein expression, p53-mediated protein expression and cellular processes showed substantial chemical-specific differences. These chemical-specific differences in the p53 transcriptional response appear to be determined by augmentation of the p53 response by co-regulators. More importantly, dose-response data for each of the chemicals indicate that the p53 transcriptional response does not prevent micronuclei induction at low concentrations. In fact, the no observed effect levels and benchmark doses for micronuclei induction were less than or equal to those for p53-mediated gene transcription regardless of the test chemical, indicating that p53's post-translational responses may be more important than transcriptional activation in the response to low dose DNA damage. This effort demonstrates the process of defining key assays required for a pathway-based, in vitro-only risk assessment, using the p53-mediated DNA damage response pathway as a prototype.

Modeling drug- and chemical-induced hepatotoxicity with systems biology approaches

We provide an overview of computational systems biology approaches as applied to the study of chemical- and drug-induced toxicity. The concept of “toxicity pathways” is described in the context of the 2007 US National Academies of Science report, “Toxicity testing in the 21st Century: A Vision and A Strategy.” Pathway mapping and modeling based on network biology concepts are a key component of the vision laid out in this report for a more biologically based analysis of dose-response behavior and the safety of chemicals and drugs. We focus on toxicity of the liver (hepatotoxicity) – a complex phenotypic response with contributions from a number of different cell types and biological processes. We describe three case studies of complementary multi-scale computational modeling approaches to understand perturbation of toxicity pathways in the human liver as a result of exposure to environmental contaminants and specific drugs. One approach involves development of a spatial, multicellular “virtual tissue” model of the liver lobule that combines molecular circuits in individual hepatocytes with cell–cell interactions and blood-mediated transport of toxicants through hepatic sinusoids, to enable quantitative, mechanistic prediction of hepatic dose-response for activation of the aryl hydrocarbon receptor toxicity pathway. Simultaneously, methods are being developing to extract quantitative maps of intracellular signaling and transcriptional regulatory networks perturbed by environmental contaminants, using a combination of gene expression and genome-wide protein-DNA interaction data. A predictive physiological model (DILIsym™) to understand drug-induced liver injury (DILI), the most common adverse event leading to termination of clinical development programs and regulatory actions on drugs, is also described. The model initially focuses on reactive metabolite-induced DILI in response to administration of acetaminophen, and spans multiple biological scales.

Functional toxicogenomics: mechanism-centered toxicology

Traditional toxicity testing using animal models is slow, low capacity, expensive and assesses a limited number of endpoints. Such approaches are inadequate to deal with the increasingly large number of compounds found in the environment for which there are no toxicity data. Mechanism-centered high-throughput testing represents an alternative approach to meet this pressing need but is limited by our current understanding of toxicity pathways. Functional toxicogenomics, the global study of the biological function of genes on the modulation of the toxic effect of a compound, can play an important role in identifying the essential cellular components and pathways involved in toxicity response. The combination of the identification of fundamental toxicity pathways and mechanism-centered targeted assays represents an integrated approach to advance molecular toxicology to meet the challenges of toxicity testing in the 21^st century.