A gene expression biomarker provides early prediction and mechanistic assessment of hepatic tumor induction by nongenotoxic chemicals

Urinary MicroRNA biomarkers of nephrotoxicity in Macaca fascicularis

Drug-induced kidney injury (DIKI) refers to kidney damage resulting from the administration of medications. The aim of this project was to identify reliable urinary microRNA (miRNAs) biomarkers that can be used as potential predictors of DIKI before disease diagnosis. This study quantified a panel of six miRNAs (miRs-210-3p, 423-5p, 143-3p, 130b-3p, 486-5p, 193a-3p) across multiple time points using urinary samples from a previous investigation evaluating effects of a nephrotoxicant in cynomolgus monkeys. Exosome-associated miRNA exhibited distinctive trends when compared to miRNAs quantified in whole urine, which may reflect a different urinary excretion mechanism of miRNAs than those released passively into the urine. Although further research and mechanistic studies are required to elucidate how these miRNAs regulate signaling in disease pathways, we present, for the first time, data that several miRNAs displayed strong correlations with histopathology scores, thus indicating their potential use as biomarkers to predict the development of DIKI in preclinical studies and clinical trials. Also, these findings can potentially be translated into other non-clinical species or human for the detection of DIKI.

Systematic transcriptome-wide meta-analysis across endocrine disrupting chemicals reveals shared and unique liver pathways, gene networks, and disease associations

Cardiometabolic disorders (CMD) are a growing public health problem across the world. Among the known cardiometabolic risk factors are compounds that induce endocrine and metabolic dysfunctions, such as endocrine disrupting chemicals (EDCs). To date, how EDCs influence molecular programs and cardiometabolic risks has yet to be fully elucidated, especially considering the complexity contributed by species-, chemical-, and dose-specific effects. Moreover, different experimental and analytical methodologies employed by different studies pose challenges when comparing findings across studies. To explore the molecular mechanisms of EDCs in a systematic manner, we established a data-driven computational approach to meta-analyze 30 human, mouse, and rat liver transcriptomic datasets for 4 EDCs, namely bisphenol A (BPA), bis(2-ethylhexyl) phthalate (DEHP), tributyltin (TBT), and perfluorooctanoic acid (PFOA). Our computational pipeline uniformly re-analyzed pre-processed quality-controlled microarray data and raw RNAseq data, derived differentially expressed genes (DEGs) and biological pathways, modeled gene regulatory networks and regulators, and determined CMD associations based on gene overlap analysis. Our approach revealed that DEHP and PFOA shared stable transcriptomic signatures that are enriched for genes associated with CMDs, suggesting similar mechanisms of action such as perturbations of peroxisome proliferator-activated receptor gamma (PPARγ) signaling and liver gene network regulators VNN1 and ACOT2. In contrast, TBT exhibited highly divergent gene signatures, pathways, network regulators, and disease associations from the other EDCs. In addition, we found that the rat, mouse, and human BPA studies showed highly variable transcriptomic patterns, providing molecular support for the variability in BPA responses. Our work offers insights into the commonality and differences in the molecular mechanisms of various EDCs and establishes a streamlined data-driven workflow to compare molecular mechanisms of environmental substances to elucidate the underlying connections between chemical exposure and disease risks.

Advancing the science of a read-across framework for evaluation of data-poor chemicals incorporating systematic and new approach methods

The assessment of human health hazards posed by chemicals traditionally relies on toxicity studies in experimental animals. However, most chemicals currently in commerce do not meet the minimum data requirements for hazard identification and dose-response analysis in human health risk assessment. Previously, we introduced a read-across framework designed to address data gaps for screening-level assessment of chemicals with insufficient in vivo toxicity information (Wang et al., 2012). It relies on inference by analogy from suitably tested source analogues to a target chemical, based on structural, toxicokinetic, and toxicodynamic similarity. This approach has been used for dose-response assessment of data-poor chemicals relevant to the U.S. EPA's Superfund program. We present herein, case studies of the application of this framework, highlighting specific examples of the use of biological similarity for chemical grouping and quantitative read-across. Based on practical knowledge and technological advances in the fields of read-across and predictive toxicology, we propose a revised framework. It includes important considerations for problem formulation, systematic review, target chemical analysis, analogue identification, analogue evaluation, and incorporation of new approach methods. This work emphasizes the integration of systematic methods and alternative toxicity testing data and tools in chemical risk assessment to inform regulatory decision-making.

A multi-label learning model for predicting drug-induced pathology in multi-organ based on toxicogenomics data

Drug-induced toxicity damages the health and is one of the key factors causing drug withdrawal from the market. It is of great significance to identify drug-induced target-organ toxicity, especially the detailed pathological findings, which are crucial for toxicity assessment, in the early stage of drug development process. A large variety of studies have devoted to identify drug toxicity. However, most of them are limited to single organ or only binary toxicity. Here we proposed a novel multi-label learning model named Att-RethinkNet, for predicting drug-induced pathological findings targeted on liver and kidney based on toxicogenomics data. The Att-RethinkNet is equipped with a memory structure and can effectively use the label association information. Besides, attention mechanism is embedded to focus on the important features and obtain better feature presentation. Our Att-RethinkNet is applicable in multiple organs and takes account the compound type, dose, and administration time, so it is more comprehensive and generalized. And more importantly, it predicts multiple pathological findings at the same time, instead of predicting each pathology separately as the previous model did. To demonstrate the effectiveness of the proposed model, we compared the proposed method with a series of state-of-the-arts methods. Our model shows competitive performance and can predict potential hepatotoxicity and nephrotoxicity in a more accurate and reliable way. The implementation of the proposed method is available at https://github.com/RanSuLab/Drug-Toxicity-Prediction-MultiLabel.

A Collaborative Initiative to Establish Genomic Biomarkers for Assessing Tumorigenic Potential to Reduce Reliance on Conventional Rodent Carcinogenicity Studies

There is growing recognition across broad sectors of the scientific community that use of genomic biomarkers has the potential to reduce the need for conventional rodent carcinogenicity studies of industrial chemicals, agrochemicals, and pharmaceuticals through a weight-of-evidence approach. These biomarkers fall into 2 major categories: (1) sets of gene transcripts that can identify distinct tumorigenic mechanisms of action; and (2) cancer driver gene mutations indicative of rapidly expanding growth-advantaged clonal cell populations. This call-to-action article describes a collaborative approach launched to develop and qualify biomarker gene expression panels that measure widely accepted molecular pathways linked to tumorigenesis and their activation levels to predict tumorigenic doses of chemicals from short-term exposures. Growing evidence suggests that application of such biomarker panels in short-term exposure rodent studies can identify both tumorigenic hazard and tumorigenic activation levels for chemical-induced carcinogenicity. In the future, this approach will be expanded to include methodologies examining mutations in key cancer driver gene mutation hotspots as biomarkers of both genotoxic and nongenotoxic chemical tumor risk. Analytical, technical, and biological validation studies of these complementary genomic tools are being undertaken by multisector and multidisciplinary collaborative teams within the Health and Environmental Sciences Institute. Success from these efforts will facilitate the transition from current heavy reliance on conventional 2-year rodent carcinogenicity studies to more rapid animal- and resource-sparing approaches for mechanism-based carcinogenicity evaluation supporting internal and regulatory decision-making.

Clearly weighing the evidence in read-across can improve assessments of data-poor chemicals

This paper provides a systematic weight-of-evidence method for read-across analyses of data-poor chemicals. The read-across technique extrapolates toxicity from analogous chemicals for which suitable test data are available to a target chemical. To determine that a candidate analogue is the 'best' and is sufficiently similar, the evidence for similarity of each candidate analogue to the target is weighed. We present a systematic weight of evidence method that provides transparency and imposes a consistent and rigorous inferential process. The method assembles relevant information concerning structure, physicochemical attributes, toxicokinetics, and toxicodynamics of the target and analogues. The information is then organized by evidence types and subtypes and weighted in terms of properties: relevance, strength, and reliability into weight levels, expressed as symbols. After evidence types are weighted, the bodies of evidence are weighted for collective properties: number, diversity, and coherence. Finally, the weights for the types and bodies of evidence are weighed for each analogue, and, if the overall weight of evidence is sufficient for one or more analogues, the analogue with the greatest weight is used to estimate the endpoint effect. We illustrate this WoE approach with a read-across analysis for screening the organochlorine contaminant, p,p'-dichlorodiphenyldichloroethane (DDD), for noncancer oral toxicity.

Antioxidant and Antiproliferative Activity of Finasteride against Glioblastoma Cells

Glioblastoma is an actively growing and aggressive brain tumor with a high propensity of recurrence. Although the surgical removal of tumor mass is the primary therapeutic option against glioblastoma, supportive pharmacotherapy is highly essential due to incredibly infiltrative characteristic of glioblastoma. Temozolomide, an FDA-approved alkylating agent, has been used as a first-line standard pharmacological approach, but several evident limitations were repeatedly reported. Despite additional therapeutic options suggested, there are no medications that successfully prevent a recurrence of glioblastoma and increase the five-year survival rate. In this study, we tested the possibility that finasteride has the potential to be developed as an anti-glioblastoma drug. Finasteride, an FDA-approved medication for the treatment of benign prostate hyperplasia and androgenic alopecia, is already known to pass through the blood-brain barrier and possess antiproliferative activity of prostate epithelial cells. We showed that finasteride inhibited the maintenance of glioma stem-like cells and repressed the proliferation of glioblastoma. Mechanistically, finasteride lowered intracellular ROS level by upregulating antioxidant genes, which contributed to inefficient β-catenin accumulation. Downregulated β-catenin resulted in the reduction in stemness and cell growth in glioblastoma.

Phosphatidylserine-specific phospholipase A1: A friend or the devil in disguise

Various human tissues and cells express phospholipase A1 member A (PLA1A), including the liver, lung, prostate gland, and immune cells. The enzyme belongs to the pancreatic lipase family. PLA1A specifically hydrolyzes sn-1 fatty acid of phosphatidylserine (PS) or 1-acyl-lysophosphatidylserine (1-acyl-lysoPS). PS externalized by activated cells or apoptotic cells or extracellular vesicles is a potential source of substrate for the production of unsaturated lysoPS species by PLA1A. Maturation and functions of many immune cells, such as T cells, dendritic cells, macrophages, and mast cells, can be regulated by PLA1A and lysoPS. Several lysoPS receptors, including GPR34, GPR174 and P2Y10, have been identified. High serum levels and high PLA1A expression are associated with autoimmune disorders such as Graves' disease and systemic lupus erythematosus. Increased expression of PLA1A is associated with metastatic melanomas. PLA1A may contribute to cardiometabolic disorders through mediating cholesterol transportation and producing lysoPS. Furthermore, PLA1A is necessary for hepatitis C virus assembly and can play a role in the antivirus innate immune response. This review summarizes recent findings on PLA1A expression, lysoPS and lysoPS receptors in autoimmune disorders, cancers, cardiometabolic disorders, antivirus immune responses, as well as regulations of immune cells.

RNA-Seq analysis of the protection by Dendrobium nobile alkaloids against carbon tetrachloride hepatotoxicity in mice

OBJECTIVE

Dendrobium nobile is a genuine Chinese medicine. Dendrobium nobile Lindl. alkaloids (DNLA) protects against CCl₄-induced acute liver injury. This study used RNA-Seq to explore the mechanisms.

METHODS

Mice were pretreated with DNLA (10 and 20 mg/kg, po) for 7 days, and subsequently intoxicated with CCl₄ (20 μL/kg, ip for 24 h). Liver RNA was extracted and subjected to RNA-Seq. The bioinformatics, including PCA, GO, KEGG, two-dimensional clustering, Ingenuity Pathways Analysis (IPA), and Illumina BaseSpace Correlation Engine (BSCE) were used to analyze the data. qPCR was performed on selected genes to verify RNA-Seq results.

RESULTS

DNLA protection against CCl₄ hepatotoxicity was confirmed by histopathology. PCA revealed the distinct gene expression patterns between the different treatment groups. GO showed that CCl₄ induced the activation, adhesion and proliferation of immune cells. KEGG showed CCl₄ induced oxidative stress, diseases and compromised adaptive responses. CCl₄ induced differentially expressed genes (DEGs) were identified by DESeq2 with Padj < 0.05 and 2D-clustered with other groups. DNLA reverted CCl₄-induced DEGs in a dose-dependent manner. qPCR analysis of S100 g, Sprr1, CCL3/7, Saa2/3, IL1rn, Cox7a2 and Rad15 confirmed RNA-Seq results. IPA showed that CCl₄ treatment altered some signaling and metabolic pathways, which were ameliorated or returned to normal following DNLA treatment. The CCl₄-activated mitochondrial oxidative phosphorylation was illustrated as an example. IPA Upstream Regulator Analysis further revealed the activated or inhibited molecules and chemicals that are responsible for CCl₄-induced DEGs, and DNLA attenuated these changes. BSCE analysis verified that CCl₄-induced DEGs were highly correlated with the GEO database of CCl₄ hepatotoxicity in rodents, and DNLA dose-dependently attenuated such correlation.

CONCLUSION

RNA-Seq revealed CCl₄-induced DEGs, disruption of canonical pathways, activation or inhibition of upstream regulators, which are highly correlated with database for CCl₄ hepatotoxicity. All these changes were attenuated or returned to normal by DNLA, demonstrating the mechanisms for DNLA to protect against CCl₄ hepatotoxicity.

Investigation of MicroRNA and transcription factor mediated regulatory network for silicosis using systems biology approach

Silicosis is a major health issue among workers exposed to crystalline silica. Genetic susceptibility has been implicated in silicosis. The present research demonstrates key regulatory targets and propagated network of gene/miRNA/transcription factor (TF) with interactions responsible for silicosis by integrating publicly available microarray data using a systems biology approach. Array quality is assessed with the Quality Metrics package of Bioconductor, limma package, and the network is constructed using Cytoscape. We observed and enlist 235 differentially expressed genes (DEGs) having up-regulation expression (85 nos) and down-regulation expression (150 nos.) in silicosis; and 24 TFs for the regulation of these DEGs entangled with thousands of miRNAs. Functional enrichment analysis of the DEGs enlighten that, the maximum number of DEGs are responsible for biological process viz, Rab proteins signal transduction (11 nos.) and Cellular Senescence (20 nos.), whereas IL-17 signaling pathway (16 nos.) and Signalling by Nuclear Receptors (14 nos.) etc. are Biological Pathway involving more DEGs. From the identified 1100 high target microRNA (miRNA)s involved in silicosis, 1055 miRNAs are found to relate with down-regulated genes and 847 miRNAs with up-regulated genes. The CDK19 gene (Up-regulated) is associated with 617 miRNAs whereas down-regulated gene ARID5B is regulated by as high as 747 high target miRNAs. In Prediction of Small-molecule signatures, maximum scoring small-molecule combinations for the DEGs have shown that CGP-60774 (with 20 combinations), alvocidib (with 15 combinations) and with AZD-7762 (24 combinations) with few other drugs having the high probability of success.

Transcriptomic analyses of livers from mice exposed to 1,4-dioxane for up to 90 days to assess potential mode(s) of action underlying liver tumor development

1,4-Dioxane is a volatile organic compound with industrial and commercial applications as a solvent and in the manufacture of other chemicals. 1,4-Dioxane has been demonstrated to induce liver tumors in chronic rodent bioassays conducted at very high doses. The available evidence for 1,4-dioxane-induced liver tumors in rodents aligns with a threshold-dependent mode of action (MOA), with the underlying mechanism being less clear in the mouse than in rats. To gain a better understanding of the underlying molecular mechanisms related to liver tumor development in mice orally exposed to 1,4-dioxane, transcriptomics analysis was conducted on liver tissue collected from a 90-day drinking water study in female B6D2F₁/Crl mice (Lafranconi et al., 2020). Using tissue samples from female mice exposed to 1,4-dioxane in the drinking water at concentrations of 0, 40, 200, 600, 2,000 or 6,000 ppm for 7, 28, and 90 days, transcriptomic analyses demonstrate minimal treatment effects on global gene expression at concentrations below 600 ppm. At higher concentrations, genes involved in phase II metabolism and mitotic cell cycle checkpoints were significantly upregulated. There was an overall lack of enrichment of genes related to DNA damage response. The increase in mitotic signaling is most prevalent in the livers of mice exposed to 1,4-dioxane at the highest concentrations for 90 days. This finding aligns with phenotypic changes reported by Lafranconi et al. (2020) after 90-days of exposure to 6,000 ppm 1,4-dioxane in the same tissues. The transcriptomics analysis further supports overarching study findings demonstrating a non-mutagenic, threshold-based, mitogenic MOA for 1,4-dioxane-induced liver tumors.

Thresholds Derived From Common Measures in Rat Studies Are Predictive of Liver Tumorigenic Chemicals

We hypothesized that typical tissue and clinical chemistry (ClinChem) end points measured in rat toxicity studies exhibit chemical-independent biological thresholds beyond which cancer occurs. Using the rat in vivo TG-GATES study, 75 chemicals were examined across chemical-dose-time comparisons that could be linked to liver tumor outcomes. Thresholds for liver weight to body weight (LW/BW) and 21 serum ClinChem end points were defined as the maximum and minimum values for those exposures that did not lead to liver tumors in rats. Upper thresholds were identified for LW/BW (117%), aspartate aminotransferase (195%), alanine aminotransferase (141%), alkaline phosphatase (152%), and total bilirubin (115%), and lower thresholds were identified for phospholipids (82%), relative albumin (93%), total cholesterol (82%), and total protein (94%). Thresholds derived from the TG-GATES data set were consistent across other acute and subchronic rat studies. A training set of ClinChem and LW/BW thresholds derived from a 38 chemical training set from TG-GATES was predictive of liver tumor outcomes for a test set of 37 independent TG-GATES chemicals (91%). The thresholds were most predictive when applied to 7d treatments (98%). These findings provide support that biological thresholds for common end points in rodent studies can be used to predict chemical tumorigenic potential.

A set of six Gene expression biomarkers and their thresholds identify rat liver tumorigens in short-term assays

Traditional methods for cancer risk assessment are retrospective, resource-intensive, and not feasible for the vast majority of environmental chemicals. In earlier studies, we used a set of six biomarkers to accurately identify liver tumorigens in transcript profiles derived from chemically-treated rats using either a Toxicological Priority Index (ToxPi) approach or using derived biomarker thresholds for cancer. The biomarkers consisting of 7-113 genes are used to predict the most common liver cancer molecular initiating events: genotoxicity, cytotoxicity and activation of the xenobiotic receptors AhR, CAR, ER, and PPARα. In the present study, we apply and evaluate the performance of these methods for cancer prediction in an independent rat liver study of 44 chemicals (6 h-7d exposures) examined by Affymetrix arrays. In the first approach, ToxPi ranking of biomarker scores consistently gave the highest scores to tumorigenic chemical-dose pairs; balanced accuracies for identification of liver tumorigenic chemicals were up to 89 %. The second approach used tumorigenic thresholds derived in the present study or from our earlier study that were set at the maximum value for chemical-dose exposures without detectable liver tumor outcomes. Using these thresholds, balanced accuracies were up to 90 %. Both approaches identified all tumorigenic chemicals. Almost all of the tumorigenic chemicals activated more than one MIE. We also compared biomarker responses between two types of profiling platforms (Affymetrix full-genome array, TempO-Seq 1500+ array containing ∼2600 genes) and found that the lack of the full set of biomarker genes on the 1500+ array resulted in decreased ability to identify chemicals that activate the MIEs. Overall, these results demonstrate that predictive approaches based on the 6 biomarkers could be used in short-term assays to identify chemicals and their doses that induce liver tumors, the most common endpoint in rodent bioassays.

A Set of Six Gene Expression Biomarkers Identify Rat Liver Tumorigens in Short-term Assays

Chemical-induced liver cancer occurs in rodents through well-characterized adverse outcome pathways. We hypothesized that measurement of the 6 most common molecular initiating events (MIEs) in liver cancer adverse outcome pathways in short-term assays using only gene expression will allow early identification of chemicals and their associated doses that are likely to be tumorigenic in the liver in 2-year bioassays. We tested this hypothesis using transcript data from a rat liver microarray compendium consisting of 2013 comparisons of 146 chemicals administered at doses with previously established effects on rat liver tumor induction. Five MIEs were measured using previously characterized gene expression biomarkers composed of gene sets predictive for genotoxicity and activation of 1 or more xenobiotic receptors (aryl hydrocarbon receptor, constitutive activated receptor, estrogen receptor, and peroxisome proliferator-activated receptor α). Because chronic injury can be important in tumorigenesis, we also developed a biomarker for cytotoxicity that had a 96% balanced accuracy. Characterization of the genes in each biomarker set using the unsupervised TXG-MAP network model demonstrated that the genes were associated with distinct functional coexpression modules. Using the Toxicological Priority Index to rank chemicals based on their ability to activate the MIEs showed that chemicals administered at tumorigenic doses clearly gave the highest ranked scores. Balanced accuracies using thresholds derived from either TG-GATES or DrugMatrix data sets to predict tumorigenicity in independent sets of chemicals were up to 93%. These results show that a MIE-directed approach using only gene expression biomarkers could be used in short-term assays to identify chemicals and their doses that cause tumors.

Pharmacodynamic and safety considerations for influenza vaccine and adjuvant design

INTRODUCTION

A novel adjuvant evaluation system for safety and immunogenicity is needed. Vaccination is important for infection prevention, for example, from influenza viruses. Adjuvants are considered critical for improving the effectiveness of influenza vaccines. Adjuvant development is an important issue in influenza vaccine design.

AREAS COVERED

A conventional in vivo evaluation method for vaccine safety has been limited in analyzing phenotypic and pathological changes. Therefore, it is difficult to obtain information on the changes at the molecular level. This review aims to explain the recently developed genomics analysis-based vaccine adjuvant safety evaluation tools verified by AddaVax^TM and polyinosinic-polycytidylic acid (poly I:C) using 18 biomarker genes and whole-virion inactivated influenza vaccine as a toxicity control. Genomics analyzes would help provide safety and efficacy information regarding influenza vaccine design by facilitating appropriate adjuvant selection.

EXPERT OPINION

The efficacy and safety profiles of influenza vaccines and adjuvants using genomics technologies provide useful information regarding immunogenicity, which is related to safety and efficacy. This approach provides important information to select appropriate inoculation routes, combinations of vaccine antigens and adjuvants, and dosing amounts. The efficacy of vaccine adjuvant evaluation by genomics analysis should be verified by various studies using various vaccines in the future.

Identification of Time-Invariant Biomarkers for Non-Genotoxic Hepatocarcinogen Assessment

Non-genotoxic hepatocarcinogens (NGHCs) can only be confirmed by 2-year rodent studies. Toxicogenomics (TGx) approaches using gene expression profiles from short-term animal studies could enable early assessment of NGHCs. However, high variance in the modulation of the genes had been noted among exposure styles and datasets. Expanding from our previous strategy in identifying consensus biomarkers in multiple experiments, we aimed to identify time-invariant biomarkers for NGHCs in short-term exposure styles and validate their applicability to long-term exposure styles. In this study, nine time-invariant biomarkers, namely A2m, Akr7a3, Aqp7, Ca3, Cdc2a, Cdkn3, Cyp2c11, Ntf3, and Sds, were identified from four large-scale microarray datasets. Machine learning techniques were subsequently employed to assess the prediction performance of the biomarkers. The biomarker set along with the Random Forest models gave the highest median area under the receiver operating characteristic curve (AUC) of 0.824 and a low interquartile range (IQR) variance of 0.036 based on a leave-one-out cross-validation. The application of the models to the external validation datasets achieved high AUC values of greater than or equal to 0.857. Enrichment analysis of the biomarkers inferred the involvement of chronic inflammatory diseases such as liver cirrhosis, fibrosis, and hepatocellular carcinoma in NGHCs. The time-invariant biomarkers provided a robust alternative for NGHC prediction.

Application of transcriptomic and microRNA profiling in the evaluation of potential liver carcinogens

Hepatocarcinogens are agents that increase the incidence of liver cancer in exposed animals or humans. It is now established that carcinogenic exposures have a widespread impact on the transcriptome, inducing both adaptive and adverse changes in the activities of genes and pathways. Chemical hepatocarcinogens have also been shown to affect expression of microRNA (miRNA), the evolutionarily conserved noncoding RNA that regulates gene expression posttranscriptionally. Considerable effort has been invested into examining the involvement of mRNA in chemical hepatocarcinogenesis and their potential usage for the classification and prediction of new chemical entities. For miRNA, there has been an increasing number of studies reported over the past decade, although not to the same degree as for transcriptomic studies. Current data suggest that it is unlikely that any gene or miRNA signature associated with short-term carcinogen exposure can replace the rodent bioassay. In this review, we discuss the application of transcriptomic and miRNA profiles to increase mechanistic understanding of chemical carcinogens and to aid in their classification.

Prospects for incorporation of epigenetic biomarkers in human health and environmental risk assessment of chemicals

Epigenetic mechanisms have gained relevance in human health and environmental studies, due to their pivotal role in disease, gene × environment interactions and adaptation to environmental change and/or contamination. Epigenetic mechanisms are highly responsive to external stimuli and a wide range of chemicals has been shown to determine specific epigenetic patterns in several organisms. Furthermore, the mitotic/meiotic inheritance of such epigenetic marks as well as the resulting changes in gene expression and cell/organismal phenotypes has now been demonstrated. Therefore, epigenetic signatures are interesting candidates for linking environmental exposures to disease as well as informing on past exposures to stressors. Accordingly, epigenetic biomarkers could be useful tools in both prospective and retrospective risk assessment but epigenetic endpoints are currently not yet incorporated into risk assessments. Achieving a better understanding on this apparent impasse, as well as identifying routes to promote the application of epigenetic biomarkers within environmental risk assessment frameworks are the objectives of this review. We first compile evidence from human health studies supporting the use of epigenetic exposure-associated changes as reliable biomarkers of exposure. Then, specifically focusing on environmental science, we examine the potential and challenges of developing epigenetic biomarkers for environmental fields, and discuss useful organisms and appropriate sequencing techniques to foster their development in this context. Finally, we discuss the practical incorporation of epigenetic biomarkers in the environmental risk assessment of chemicals, highlighting critical data gaps and making key recommendations for future research within a regulatory context.

The Carcinogenome Project: In Vitro Gene Expression Profiling of Chemical Perturbations to Predict Long-Term Carcinogenicity

BACKGROUND

Most chemicals in commerce have not been evaluated for their carcinogenic potential. The de facto gold-standard approach to carcinogen testing adopts the 2-y rodent bioassay, a time-consuming and costly procedure. High-throughput in vitro assays are a promising alternative for addressing the limitations in carcinogen screening.

OBJECTIVES

We developed a screening process for predicting chemical carcinogenicity and genotoxicity and characterizing modes of actions (MoAs) using in vitro gene expression assays.

METHODS

We generated a large toxicogenomics resource comprising [Formula: see text] expression profiles corresponding to 330 chemicals profiled in HepG2 (human hepatocellular carcinoma cell line) at multiple doses and replicates. Predictive models of carcinogenicity and genotoxicity were built using a random forest classifier. Differential pathway enrichment analysis was performed to identify pathways associated with carcinogen exposure. Signatures of carcinogenicity and genotoxicity were compared with external sources, including Drugmatrix and the Connectivity Map.

RESULTS

Among profiles with sufficient bioactivity, our classifiers achieved 72.2% Area Under the ROC Curve (AUC) for predicting carcinogenicity and 82.3% AUC for predicting genotoxicity. Chemical bioactivity, as measured by the strength and reproducibility of the transcriptional response, was not significantly associated with long-term carcinogenicity in doses up to [Formula: see text]. However, sufficient bioactivity was necessary for a chemical to be used for prediction of carcinogenicity. Pathway enrichment analysis revealed pathways consistent with known pathways that drive cancer, including DNA damage and repair. The data is available at https://clue.io/CRCGN_ABC , and a portal for query and visualization of the results is accessible at https://carcinogenome.org .

DISCUSSION

We demonstrated an in vitro screening approach using gene expression profiling to predict carcinogenicity and infer MoAs of chemical perturbations. https://doi.org/10.1289/EHP3986.

Resveratrol Improves Recovery and Survival of Diet-Induced Obese Mice Undergoing Extended Major (80%) Hepatectomy

INTRODUCTION

Loss of hepatic epidermal growth factor receptor (EGFR) expression is a cause for the increased perioperative risk for complications and death in patients with obesity and fatty liver undergoing liver resection. Herein, we set out to identify agents that might increase EGFR expression and improve recovery for patients with fatty liver undergoing resection. Using the diet-induced obese (DIO) mouse model of fatty liver, we examined resveratrol as a therapy to induce EGFR expression and improve outcomes following 80% partial hepatectomy (PH) in a murine model.

METHODS

DIO mice were fed resveratrol or carrier control by gavage. EGFR expression and the response to major (80%) PH were examined.

RESULTS

Based on an Illumina analysis, resveratrol was identified as increasing EGFR gene expression in A549 cells. Resveratrol was observed to also increase EGFR protein expression in A549 cells. DIO mice fed resveratrol by gavage (75 mg/kg) demonstrated an increased EGFR expression without the identified hepatic toxicity. Resveratrol and control mice subjected to 80% PH, a model of high mortality hepatectomy in DIO mice, demonstrated macroscopically decreased fatty liver and fewer liver hemorrhagic petechiae. Resveratrol pretreatment ameliorated liver injury and accelerated regeneration of the hepatic remnant after 80% PH including decreasing serum ALT and bilirubin, while increasing hepatic PCNA expression. Resveratrol increased induction of p-STAT3 and p-AKT after 80% hepatectomy. Resveratrol pretreatment significantly improved survival rates in DIO mice undergoing extended 80% PH.

CONCLUSIONS

Oral resveratrol restores EGFR expression in fatty liver. Resveratrol may be a promising protective agent in instances where extensive hepatic resection of fatty liver is required.

Robust Co-clustering to Discover Toxicogenomic Biomarkers and Their Regulatory Doses of Chemical Compounds Using Logistic Probabilistic Hidden Variable Model

Detection of biomarker genes and their regulatory doses of chemical compounds (DCCs) is one of the most important tasks in toxicogenomic studies as well as in drug design and development. There is an online computational platform "Toxygates" to identify biomarker genes and their regulatory DCCs by co-clustering approach. Nevertheless, the algorithm of that platform based on hierarchical clustering (HC) does not share gene-DCC two-way information simultaneously during co-clustering between genes and DCCs. Also it is sensitive to outlying observations. Thus, this platform may produce misleading results in some cases. The probabilistic hidden variable model (PHVM) is a more effective co-clustering approach that share two-way information simultaneously, but it is also sensitive to outlying observations. Therefore, in this paper we have proposed logistic probabilistic hidden variable model (LPHVM) for robust co-clustering between genes and DCCs, since gene expression data are often contaminated by outlying observations. We have investigated the performance of the proposed LPHVM co-clustering approach in a comparison with the conventional PHVM and Toxygates co-clustering approaches using simulated and real life TGP gene expression datasets, respectively. Simulation results show that the proposed method improved the performance over the conventional PHVM in presence of outliers; otherwise, it keeps equal performance. In the case of real life TGP data analysis, three DCCs (glibenclamide-low, perhexilline-low, and hexachlorobenzene-medium) for glutathione metabolism pathway dataset as well as two DCCs (acetaminophen-medium and methapyrilene-low) for PPAR signaling pathway dataset were incorrectly co-clustered by the Toxygates online platform, while only one DCC (hexachlorobenzene-low) for glutathione metabolism pathway was incorrectly co-clustered by the proposed LPHVM approach. Our findings from the real data analysis are also supported by the other findings in the literature.

Establishment of a novel safety assessment method for vaccine adjuvant development

Vaccines effectively prevent infectious diseases. Many types of vaccines against various pathogens that threaten humans are currently in widespread use. Recently, adjuvant adaptation has been attempted to activate innate immunity to enhance the effectiveness of vaccines. The effectiveness of adjuvants for vaccinations has been demonstrated in many animal models and clinical trials. Although a highly potent adjuvant tends to have high effectiveness, it also has the potential to increase the risk of side effects such as pain, edema, and fever. Indeed, highly effective adjuvants, such as poly(I:C), have not been clinically applied due to their high risks of toxicity in humans. Therefore, the task in the field of adjuvant development is to clinically apply highly effective and non- or low-toxic adjuvant-containing vaccines. To resolve this issue, it is essential to ensure a low risk of side effects and the high efficacy of an adjuvant in the early developmental phases. This review summarizes the theory and history of the current safety assessment methods for adjuvants, using the inactivated influenza vaccine as a model. Our novel method was developed as a system to judge the safety of a candidate compound using biomarkers identified by genomic technology and statistical tools. A systematic safety assessment tool for adjuvants would be of great use for predicting toxicity during novel adjuvant development, screening, and quality control.

Transcriptional Responses Reveal Similarities Between Preclinical Rat Liver Testing Systems

Toxicogenomics (TGx) is an important tool to gain an enhanced understanding of toxicity at the molecular level. Previously, we developed a pair ranking (PRank) method to assess in vitro to in vivo extrapolation (IVIVE) using toxicogenomic datasets from the Open Toxicogenomics Project-Genomics Assisted Toxicity Evaluation System (TG-GATEs) database. With this method, we investiagted three important questions that were not addressed in our previous study: (1) is a 1-day in vivo short-term assay able to replace the 28-day standard and expensive toxicological assay? (2) are some biological processes more conservative across different preclinical testing systems than others? and (3) do these preclinical testing systems have the similar resolution in differentiating drugs by their therapeutic uses? For question 1, a high similarity was noted (PRank score = 0.90), indicating the potential utility of shorter term in vivo studies to predict outcome in longer term and more expensive in vivo model systems. There was a moderate similarity between rat primary hepatocytes and in vivo repeat-dose studies (PRank score = 0.71) but a low similarity (PRank score = 0.56) between rat primary hepatocytes and in vivo single dose studies. To address question 2, we limited the analysis to gene sets relevant to specific toxicogenomic pathways and we found that pathways such as lipid metabolism were consistently over-represented in all three assay systems. For question 3, all three preclinical assay systems could distinguish compounds from different therapeutic categories. This suggests that any noted differences in assay systems was biological process-dependent and furthermore that all three systems have utility in assessing drug responses within a certain drug class. In conclusion, this comparison of three commonly used rat TGx systems provides useful information in utility and application of TGx assays.

Meloxicam increases epidermal growth factor receptor expression improving survival after hepatic resection in diet-induced obese mice

BACKGROUND

Patients with fatty liver have delayed regenerative responses, increased hepatocellular injury, and increased risk for perioperative mortality. Currently, no clinical therapy exists to prevent liver failure or improve regeneration in patients with fatty liver. Previously we demonstrated that obese mice have markedly reduced levels of epidermal growth factor receptor in liver. We sought to identify pharmacologic agents to increase epidermal growth factor receptor expression to improve hepatic regeneration in the setting of fatty liver resection.

METHODS

Lean (20% calories from fat) and diet-induced obese mice (60% calories from fat) were subjected to 70% or 80% hepatectomy.

RESULTS

Using the BaseSpace Correlation Engine of deposited gene arrays we identified agents that increased hepatic epidermal growth factor receptor. Meloxicam was identified as inducing epidermal growth factor receptor expression across species. Meloxicam improved hepatic steatosis in diet-induced obese mice both grossly and histologically. Immunohistochemistry and Western blot analysis demonstrated that meloxicam pretreatment of diet-induced obese mice dramatically increased epidermal growth factor receptor protein expression in hepatocytes. After 70% hepatectomy, meloxicam pretreatment ameliorated liver injury and significantly accelerated mitotic rates of hepatocytes in obese mice. Recovery of liver mass was accelerated in obese mice pretreated with meloxicam (by 26% at 24 hours and 38% at 48 hours, respectively). After 80% hepatectomy, survival was dramatically increased with meloxicam treatment.

CONCLUSION

Low epidermal growth factor receptor expression is a common feature of fatty liver disease. Meloxicam restores epidermal growth factor receptor expression in steatotic hepatocytes. Meloxicam pretreatment may be applied to improve outcome after fatty liver resection or transplantation with steatotic graft.

Identification of consensus biomarkers for predicting non-genotoxic hepatocarcinogens

The assessment of non-genotoxic hepatocarcinogens (NGHCs) is currently relying on two-year rodent bioassays. Toxicogenomics biomarkers provide a potential alternative method for the prioritization of NGHCs that could be useful for risk assessment. However, previous studies using inconsistently classified chemicals as the training set and a single microarray dataset concluded no consensus biomarkers. In this study, 4 consensus biomarkers of A2m, Ca3, Cxcl1, and Cyp8b1 were identified from four large-scale microarray datasets of the one-day single maximum tolerated dose and a large set of chemicals without inconsistent classifications. Machine learning techniques were subsequently applied to develop prediction models for NGHCs. The final bagging decision tree models were constructed with an average AUC performance of 0.803 for an independent test. A set of 16 chemicals with controversial classifications were reclassified according to the consensus biomarkers. The developed prediction models and identified consensus biomarkers are expected to be potential alternative methods for prioritization of NGHCs for further experimental validation.

Detection of non-genotoxic hepatocarcinogens and prediction of their mechanism of action in rats using gene marker sets

Several studies have successfully detected hepatocarcinogenicity in rats based on gene expression data. However, prediction of hepatocarcinogens with certain mechanisms of action (MOAs), such as enzyme inducers and peroxisome proliferator-activated receptor α (PPARα) agonists, can prove difficult using a single model and requires a highly toxic dose. Here, we constructed a model for detecting non-genotoxic (NGTX) hepatocarcinogens and predicted their MOAs in rats. Gene expression data deposited in the Open Toxicogenomics Project-Genomics Assisted Toxicity Evaluation System (TG-GATEs) was used to investigate gene marker sets. Principal component analysis (PCA) was applied to discriminate different MOAs, and a support vector machine algorithm was applied to construct the prediction model. This approach identified 106 probe sets as gene marker sets for PCA and enabled the prediction model to be constructed. In PCA, NGTX hepatocarcinogens were classified as follows based on their MOAs: cytotoxicants, PPARα agonists, or enzyme inducers. The prediction model detected hepatocarcinogenicity with an accuracy of more than 90% in 14- and 28-day repeated-dose studies. In addition, the doses capable of predicting NGTX hepatocarcinogenicity were close to those required in rat carcinogenicity assays. In conclusion, our PCA and prediction model using gene marker sets will help assess the risk of hepatocarcinogenicity in humans based on MOAs and reduce the number of two-year rodent bioassays.

Critical role of toxicologic pathology in a short-term screen for carcinogenicity

Carcinogenic potential of chemicals is currently evaluated using a two year bioassay in rodents. Numerous difficulties are known for this assay, most notably, the lack of information regarding detailed dose response and human relevance of any positive findings. A screen for carcinogenic activity has been proposed based on a 90 day screening assay. Chemicals are first evaluated for proliferative activity in various tissues. If negative, lack of carcinogenic activity can be concluded. If positive, additional evaluation for DNA reactivity, immunosuppression, and estrogenic activity are evaluated. If these are negative, additional efforts are made to determine specific modes of action in the animal model, with a detailed evaluation of the potential relevance to humans. Applications of this approach are presented for liver and urinary bladder. Toxicologic pathology is critical for all of these evaluations, including a detailed histopathologic evaluation of the 90 day assay, immunohistochemical analyses for labeling index, and involvement in a detailed mode of action analysis. Additionally, the toxicologic pathologist needs to be involved with molecular evaluations and evaluations of new molecularly developed animal models. The toxicologic pathologist is uniquely qualified to provide the expertise needed for these evaluations.

Prediction of Non-Genotoxic Carcinogenicity Based on Genetic Profiles of Short Term Exposure Assays

Non-genotoxic carcinogens are substances that induce tumorigenesis by non-mutagenic mechanisms and long term rodent bioassays are required to identify them. Recent studies have shown that transcription profiling can be applied to develop early identifiers for long term phenotypes. In this study, we used rat liver expression profiles from the NTP (National Toxicology Program, Research Triangle Park, USA) DrugMatrix Database to construct a gene classifier that can distinguish between non-genotoxic carcinogens and other chemicals. The model was based on short term exposure assays (3 days) and the training was limited to oxidative stressors, peroxisome proliferators and hormone modulators. Validation of the predictor was performed on independent toxicogenomic data (TG-GATEs, Toxicogenomics Project-Genomics Assisted Toxicity Evaluation System, Osaka, Japan). To build our model we performed Random Forests together with a recursive elimination algorithm (VarSelRF). Gene set enrichment analysis was employed for functional interpretation. A total of 770 microarrays comprising 96 different compounds were analyzed and a predictor of 54 genes was built. Prediction accuracy was 0.85 in the training set, 0.87 in the test set and increased with increasing concentration in the validation set: 0.6 at low dose, 0.7 at medium doses and 0.81 at high doses. Pathway analysis revealed gene prominence of cellular respiration, energy production and lipoprotein metabolism. The biggest target of toxicogenomics is accurately predict the toxicity of unknown drugs. In this analysis, we presented a classifier that can predict non-genotoxic carcinogenicity by using short term exposure assays. In this approach, dose level is critical when evaluating chemicals at early time points.

An integrative test strategy for cancer hazard identification

Assessment of genotoxic and carcinogenic potential is considered one of the basic requirements when evaluating possible human health risks associated with exposure to chemicals. Test strategies currently in place focus primarily on identifying genotoxic potential due to the strong association between the accumulation of genetic damage and cancer. Using genotoxicity assays to predict carcinogenic potential has the significant drawback that risks from non-genotoxic carcinogens remain largely undetected unless carcinogenicity studies are performed. Furthermore, test systems already developed to reduce animal use are not easily accepted and implemented by either industries or regulators. This manuscript reviews the test methods for cancer hazard identification that have been adopted by the regulatory authorities, and discusses the most promising alternative methods that have been developed to date. Based on these findings, a generally applicable tiered test strategy is proposed that can be considered capable of detecting both genotoxic as well as non-genotoxic carcinogens and will improve understanding of the underlying mode of action. Finally, strengths and weaknesses of this new integrative test strategy for cancer hazard identification are presented.

Mechanism Profiling of Hepatotoxicity Caused by Oxidative Stress Using Antioxidant Response Element Reporter Gene Assay Models and Big Data

BACKGROUND

Hepatotoxicity accounts for a substantial number of drugs being withdrawn from the market. Using traditional animal models to detect hepatotoxicity is expensive and time-consuming. Alternative in vitro methods, in particular cell-based high-throughput screening (HTS) studies, have provided the research community with a large amount of data from toxicity assays. Among the various assays used to screen potential toxicants is the antioxidant response element beta lactamase reporter gene assay (ARE-bla), which identifies chemicals that have the potential to induce oxidative stress and was used to test > 10,000 compounds from the Tox21 program.

OBJECTIVE

The ARE-bla computational model and HTS data from a big data source (PubChem) were used to profile environmental and pharmaceutical compounds with hepatotoxicity data.

METHODS

Quantitative structure-activity relationship (QSAR) models were developed based on ARE-bla data. The models predicted the potential oxidative stress response for known liver toxicants when no ARE-bla data were available. Liver toxicants were used as probe compounds to search PubChem Bioassay and generate a response profile, which contained thousands of bioassays (> 10 million data points). By ranking the in vitro-in vivo correlations (IVIVCs), the most relevant bioassay(s) related to hepatotoxicity were identified.

RESULTS

The liver toxicants profile contained the ARE-bla and relevant PubChem assays. Potential toxicophores for well-known toxicants were created by identifying chemical features that existed only in compounds with high IVIVCs.

CONCLUSION

Profiling chemical IVIVCs created an opportunity to fully explore the source-to-outcome continuum of modern experimental toxicology using cheminformatics approaches and big data sources.

CITATION

Kim MT, Huang R, Sedykh A, Wang W, Xia M, Zhu H. 2016. Mechanism profiling of hepatotoxicity caused by oxidative stress using antioxidant response element reporter gene assay models and big data. Environ Health Perspect 124:634-641; http://dx.doi.org/10.1289/ehp.1509763.

Risk assessment of Soulatrolide and Mammea (A/BA+A/BB) coumarins from Calophyllum brasiliense by a toxicogenomic and toxicological approach

Calophyllum brasiliense (Calophyllaceae) is a tropical rain forest tree distributed in Central and South America. It is an important source of tetracyclic dipyrano coumarins (Soulatrolide) and Mammea type coumarins. Soulatrolide is a potent inhibitor of HIV-1 reverse transcriptase and displays activity against Mycobacterium tuberculosis. Meanwhile, Mammea A/BA and A/BB, pure or as a mixture, are highly active against several human leukemia cell lines, Trypanosoma cruzi and Leishmania amazonensis. Nevertheless, there are few studies evaluating their safety profile. In the present work we performed toxicogenomic and toxicological analysis for both type of compounds. Soulatrolide, and the Mammea A/BA + A/BB mixture (2.1) were slightly toxic accordingly to Lorke assay classification (DL50 > 3000 mg/kg). After a short-term administration (100 mg/kg/daily, orally, 1 week) liver toxicogenomic analysis revealed 46 up and 72 downregulated genes for Mammea coumarins, and 665 up and 1077 downregulated genes for Soulatrolide. Gene enrichment analysis identified transcripts involved in drug metabolism for both compounds. In addition, network analysis through protein-protein interactions, tissue evaluation by TUNEL assay, and histological examination revealed no tissue damage on liver, kidney and spleen after treatments. Our results indicate that both type of coumarins displayed a safety profile, supporting their use in further preclinical studies to determine its therapeutic potential.

Potential toxic effects of glyphosate and its commercial formulations below regulatory limits

Glyphosate-based herbicides (GlyBH), including Roundup, are the most widely used pesticides worldwide. Their uses have increased exponentially since their introduction on the market. Residue levels in food or water, as well as human exposures, are escalating. We have reviewed the toxic effects of GlyBH measured below regulatory limits by evaluating the published literature and regulatory reports. We reveal a coherent body of evidence indicating that GlyBH could be toxic below the regulatory lowest observed adverse effect level for chronic toxic effects. It includes teratogenic, tumorigenic and hepatorenal effects. They could be explained by endocrine disruption and oxidative stress, causing metabolic alterations, depending on dose and exposure time. Some effects were detected in the range of the recommended acceptable daily intake. Toxic effects of commercial formulations can also be explained by GlyBH adjuvants, which have their own toxicity, but also enhance glyphosate toxicity. These challenge the assumption of safety of GlyBH at the levels at which they contaminate food and the environment, albeit these levels may fall below regulatory thresholds. Neurodevelopmental, reproductive, and transgenerational effects of GlyBH must be revisited, since a growing body of knowledge suggests the predominance of endocrine disrupting mechanisms caused by environmentally relevant levels of exposure.

Bayesian Network Inference Enables Unbiased Phenotypic Anchoring of Transcriptomic Responses to Cigarette Smoke in Humans

Microarray-based transcriptomic analysis has been demonstrated to hold the opportunity to study the effects of human exposure to, e.g., chemical carcinogens at the whole genome level, thus yielding broad-ranging molecular information on possible carcinogenic effects. Since genes do not operate individually but rather through concerted interactions, analyzing and visualizing networks of genes should provide important mechanistic information, especially upon connecting them to functional parameters, such as those derived from measurements of biomarkers for exposure and carcinogenic risk. Conventional methods such as hierarchical clustering and correlation analyses are frequently used to address these complex interactions but are limited as they do not provide directional causal dependence relationships. Therefore, our aim was to apply Bayesian network inference with the purpose of phenotypic anchoring of modified gene expressions. We investigated a use case on transcriptomic responses to cigarette smoking in humans, in association with plasma cotinine levels as biomarkers of exposure and aromatic DNA-adducts in blood cells as biomarkers of carcinogenic risk. Many of the genes that appear in the Bayesian networks surrounding plasma cotinine, and to a lesser extent around aromatic DNA-adducts, hold biologically relevant functions in inducing severe adverse effects of smoking. In conclusion, this study shows that Bayesian network inference enables unbiased phenotypic anchoring of transcriptomics responses. Furthermore, in all inferred Bayesian networks several dependencies are found which point to known but also to new relationships between the expression of specific genes, cigarette smoke exposure, DNA damaging-effects, and smoking-related diseases, in particular associated with apoptosis, DNA repair, and tumor suppression, as well as with autoimmunity.

The adverse outcome pathway for rodent liver tumor promotion by sustained activation of the aryl hydrocarbon receptor

An Adverse Outcome Pathway (AOP) represents the existing knowledge of a biological pathway leading from initial molecular interactions of a toxicant and progressing through a series of key events (KEs), culminating with an apical adverse outcome (AO) that has to be of regulatory relevance. An AOP based on the mode of action (MOA) of rodent liver tumor promotion by dioxin-like compounds (DLCs) has been developed and the weight of evidence (WoE) of key event relationships (KERs) evaluated using evolved Bradford Hill considerations. Dioxins and DLCs are potent aryl hydrocarbon receptor (AHR) ligands that cause a range of species-specific adverse outcomes. The occurrence of KEs is necessary for inducing downstream biological responses and KEs may occur at the molecular, cellular, tissue and organ levels. The common convention is that an AOP begins with the toxicant interaction with a biological response element; for this AOP, this initial event is binding of a DLC ligand to the AHR. Data from mechanistic studies, lifetime bioassays and approximately thirty initiation-promotion studies have established dioxin and DLCs as rat liver tumor promoters. Such studies clearly show that sustained AHR activation, weeks or months in duration, is necessary to induce rodent liver tumor promotion--hence, sustained AHR activation is deemed the molecular initiating event (MIE). After this MIE, subsequent KEs are 1) changes in cellular growth homeostasis likely associated with expression changes in a number of genes and observed as development of hepatic foci and decreases in apoptosis within foci; 2) extensive liver toxicity observed as the constellation of effects called toxic hepatopathy; 3) cellular proliferation and hyperplasia in several hepatic cell types. This progression of KEs culminates in the AO, the development of hepatocellular adenomas and carcinomas and cholangiolar carcinomas. A rich data set provides both qualitative and quantitative knowledge of the progression of this AOP through KEs and the KERs. Thus, the WoE for this AOP is judged to be strong. Species-specific effects of dioxins and DLCs are well known--humans are less responsive than rodents and rodent species differ in sensitivity between strains. Consequently, application of this AOP to evaluate potential human health risks must take these differences into account.

Release of (and lessons learned from mining) a pioneering large toxicogenomics database

AIM

We release the Janssen Toxicogenomics database. This rat liver gene-expression database was generated using Codelink microarrays, and has been used over the past years within Janssen to derive signatures for multiple end points and to classify proprietary compounds.

MATERIALS & METHODS

The release consists of gene-expression responses to 124 compounds, selected to give a broad coverage of liver-active compounds. A selection of the compounds were also analyzed on Affymetrix microarrays.

RESULTS

The release includes results of an in-house reannotation pipeline to Entrez gene annotations, to classify probes into different confidence classes. High confidence unambiguously annotated probes were used to create gene-level data which served as starting point for cross-platform comparisons. Connectivity map-based similarity methods show excellent agreement between Codelink and Affymetrix runs of the same samples. We also compared our dataset with the Japanese Toxicogenomics Project and observed reasonable agreement, especially for compounds with stronger gene signatures. We describe an R-package containing the gene-level data and show how it can be used for expression-based similarity searches.

CONCLUSION

Comparing the same biological samples run on the Affymetrix and the Codelink platform, good correspondence is observed using connectivity mapping approaches. As expected, this correspondence is smaller when the data are compared with an independent dataset such as TG-GATE. We hope that this collection of gene-expression profiles will be incorporated in toxicogenomics pipelines of users.

Technical guide for applications of gene expression profiling in human health risk assessment of environmental chemicals

Toxicogenomics promises to be an important part of future human health risk assessment of environmental chemicals. The application of gene expression profiles (e.g., for hazard identification, chemical prioritization, chemical grouping, mode of action discovery, and quantitative analysis of response) is growing in the literature, but their use in formal risk assessment by regulatory agencies is relatively infrequent. Although additional validations for specific applications are required, gene expression data can be of immediate use for increasing confidence in chemical evaluations. We believe that a primary reason for the current lack of integration is the limited practical guidance available for risk assessment specialists with limited experience in genomics. The present manuscript provides basic information on gene expression profiling, along with guidance on evaluating the quality of genomic experiments and data, and interpretation of results presented in the form of heat maps, pathway analyses and other common approaches. Moreover, potential ways to integrate information from gene expression experiments into current risk assessment are presented using published studies as examples. The primary objective of this work is to facilitate integration of gene expression data into human health risk assessments of environmental chemicals.

Proteome mapping of epidermal growth factor induced hepatocellular carcinomas identifies novel cell metabolism targets and mitogen activated protein kinase signalling events

Background

Hepatocellular carcinoma (HCC) is on the rise and the sixth most common cancer worldwide. To combat HCC effectively research is directed towards its early detection and the development of targeted therapies. Given the fact that epidermal growth factor (EGF) is an important mitogen for hepatocytes we searched for disease regulated proteins to improve an understanding of the molecular pathogenesis of EGF induced HCC. Disease regulated proteins were studied by 2DE MALDI-TOF/TOF and a transcriptomic approach, by immunohistochemistry and advanced bioinformatics.

Results

Mapping of EGF induced liver cancer in a transgenic mouse model identified n = 96 (p < 0.05) significantly regulated proteins of which n = 54 were tumour-specific. To unravel molecular circuits linked to aberrant EGFR signalling diverse computational approaches were employed and this defined n = 7 key nodes using n = 82 disease regulated proteins for network construction. STRING analysis revealed protein-protein interactions of > 70% disease regulated proteins with individual proteins being validated by immunohistochemistry. The disease regulated network proteins were mapped to distinct pathways and bioinformatics provided novel insight into molecular circuits associated with significant changes in either glycolysis and gluconeogenesis, argine and proline metabolism, protein processing in endoplasmic reticulum, Hif- and MAPK signalling, lipoprotein metabolism, platelet activation and hemostatic control as a result of aberrant EGF signalling. The biological significance of the findings was corroborated with gene expression data derived from tumour tissues to evntually define a rationale by which tumours embark on intriguing changes in metabolism that is of utility for an understanding of tumour growth. Moreover, among the EGF tumour specific proteins n = 11 were likewise uniquely expressed in human HCC and for n = 49 proteins regulation in human HCC was confirmed using the publically available Human Protein Atlas depository, therefore demonstrating clinical significance.

Conclusion

Novel insight into the molecular pathogenesis of EGF induced liver cancer was obtained and among the 37 newly identified proteins several are likely candidates for the development of molecularly targeted therapies and include the nucleoside diphosphate kinase A, bifunctional ATP-dependent dihydroyacetone kinase and phosphatidylethanolamine-binding protein1, the latter being an inhibitor of the Raf-1 kinase.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1312-z) contains supplementary material, which is available to authorized users.

Exploiting dependencies of pairwise comparison outcomes to predict patterns of gene response

BACKGROUND

The analysis of gene expression has played an important role in medical and bioinformatics research. Although it is known that a large number of samples is needed to determine the patterns of gene expression accurately, practical designs of gene expression studies occasionally have insufficient numbers of samples, making it difficult to ascertain true response patterns of variantly expressed genes.

RESULTS

We describe an approach to cope with the challenge of predicting true orders of gene response to treatments. We show that true patterns of gene response must be orderable sets. In experiments with few samples, we modify the conventional pairwise comparison tests and increase the significance level α intelligently to deduce orderable patterns, which are most likely true orders of gene response. Additionally, motivated by the fact that a gene can be involved in multiple biological functions, our method further resamples experimental replicates and predicts multiple response patterns for each gene.

CONCLUSIONS

This method can be useful in designing cost-effective experiments with small sample sizes. Patterns of highly-variantly expressed genes can be predicted by varying α intelligently. Furthermore, clusters are labeled meaningfully with patterns that describe precisely how genes in such clusters respond to treatments.

Genomic models of short-term exposure accurately predict long-term chemical carcinogenicity and identify putative mechanisms of action

BACKGROUND

Despite an overall decrease in incidence of and mortality from cancer, about 40% of Americans will be diagnosed with the disease in their lifetime, and around 20% will die of it. Current approaches to test carcinogenic chemicals adopt the 2-year rodent bioassay, which is costly and time-consuming. As a result, fewer than 2% of the chemicals on the market have actually been tested. However, evidence accumulated to date suggests that gene expression profiles from model organisms exposed to chemical compounds reflect underlying mechanisms of action, and that these toxicogenomic models could be used in the prediction of chemical carcinogenicity.

RESULTS

In this study, we used a rat-based microarray dataset from the NTP DrugMatrix Database to test the ability of toxicogenomics to model carcinogenicity. We analyzed 1,221 gene-expression profiles obtained from rats treated with 127 well-characterized compounds, including genotoxic and non-genotoxic carcinogens. We built a classifier that predicts a chemical's carcinogenic potential with an AUC of 0.78, and validated it on an independent dataset from the Japanese Toxicogenomics Project consisting of 2,065 profiles from 72 compounds. Finally, we identified differentially expressed genes associated with chemical carcinogenesis, and developed novel data-driven approaches for the molecular characterization of the response to chemical stressors.

CONCLUSION

Here, we validate a toxicogenomic approach to predict carcinogenicity and provide strong evidence that, with a larger set of compounds, we should be able to improve the sensitivity and specificity of the predictions. We found that the prediction of carcinogenicity is tissue-dependent and that the results also confirm and expand upon previous studies implicating DNA damage, the peroxisome proliferator-activated receptor, the aryl hydrocarbon receptor, and regenerative pathology in the response to carcinogen exposure.

The Predictive Safety Testing Consortium: A synthesis of the goals, challenges and accomplishments of the Critical Path

The qualification of biomarkers of drug safety requires data on many compounds and nonclinical and clinical studies. The cost and effort associated with these qualifications cannot be easily covered by a single pharmaceutical company. Intellectual property associated with safety biomarkers is also held by many different companies. Consortia between different pharmaceutical companies can overcome cost and intellectual property hurdles to biomarker qualification. The Predictive Safety Testing Consortium (PSTC) is a collaborative effort between 16 different pharmaceutical companies to generate data supporting biomarker qualification. This Consortium is coordinated through the C-Path Institute, and currently has five biomarker qualification working groups engaged in this collaboration: nephrotoxicity, hepatotoxicity, vascular injury, myopathy, and non-genotoxic carcinogenicity. These working groups are aided by a data management team and a translational strategy team. Qualification studies of promising biomarkers are already progressing in several of the working groups, and results in the nephrotoxicity working group warranted a data submission to the FDA and EMEA for regulatory qualification of new nephrotoxicity biomarkers.:

Cross-platform toxicogenomics for the prediction of non-genotoxic hepatocarcinogenesis in rat

In the area of omics profiling in toxicology, i.e. toxicogenomics, characteristic molecular profiles have previously been incorporated into prediction models for early assessment of a carcinogenic potential and mechanism-based classification of compounds. Traditionally, the biomarker signatures used for model construction were derived from individual high-throughput techniques, such as microarrays designed for monitoring global mRNA expression. In this study, we built predictive models by integrating omics data across complementary microarray platforms and introduced new concepts for modeling of pathway alterations and molecular interactions between multiple biological layers. We trained and evaluated diverse machine learning-based models, differing in the incorporated features and learning algorithms on a cross-omics dataset encompassing mRNA, miRNA, and protein expression profiles obtained from rat liver samples treated with a heterogeneous set of substances. Most of these compounds could be unambiguously classified as genotoxic carcinogens, non-genotoxic carcinogens, or non-hepatocarcinogens based on evidence from published studies. Since mixed characteristics were reported for the compounds Cyproterone acetate, Thioacetamide, and Wy-14643, we reclassified these compounds as either genotoxic or non-genotoxic carcinogens based on their molecular profiles. Evaluating our toxicogenomics models in a repeated external cross-validation procedure, we demonstrated that the prediction accuracy of our models could be increased by joining the biomarker signatures across multiple biological layers and by adding complex features derived from cross-platform integration of the omics data. Furthermore, we found that adding these features resulted in a better separation of the compound classes and a more confident reclassification of the three undefined compounds as non-genotoxic carcinogens.