The liver pharmacological and xenobiotic gene response repertoire

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.

Kinetic Modeling of Hepatic Metabolism and Simulation of Treatment Effects

Mathematical modeling is a promising strategy to fill the experimentally unapproachable knowledge gaps about the relative contribution of various molecular processes to cellular metabolic function. To this end, we developed detailed kinetic models of the central metabolism of different cell types, comprising multiple metabolic functionalities. We used the model to simulate metabolic changes in several cell types under different experimental settings in health and disease. In this way, we show that it is possible to decipher and characterize the relative influence of various metabolic pathways and enzymes to overall metabolic performance and phenotype.Quantitative Systems Metabolism (QSM™) allows quantitative assessment of metabolic functionality and metabolic profiling based on proteomic data. Here, we describe the technique, namely, molecular resolved kinetic modeling, underlying QSM™. We explain the necessary steps for the generation of cell-specific models to functionally interpret proteomic data and point out some unresolved challenges and open questions.

Toxicogenomics Approaches to Address Toxicity and Carcinogenicity in the Liver

Toxicogenomic technologies query the genome, transcriptome, proteome, and the epigenome in a variety of toxicological conditions. Due to practical considerations related to the dynamic range of the assays, sensitivity, cost, and technological limitations, transcriptomic approaches are predominantly used in toxicogenomics. Toxicogenomics is being used to understand the mechanisms of toxicity and carcinogenicity, evaluate the translational relevance of toxicological responses from in vivo and in vitro models, and identify predictive biomarkers of disease and exposure. In this session, a brief overview of various transcriptomic technologies and practical considerations related to experimental design was provided. The advantages of gene network analyses to define mechanisms were also discussed. An assessment of the utility of toxicogenomic technologies in the environmental and pharmaceutical space showed that these technologies are being increasingly used to gain mechanistic insights and determining the translational relevance of adverse findings. Within the environmental toxicology area, there is a broader regulatory consideration of benchmark doses derived from toxicogenomics data. In contrast, these approaches are mainly used for internal decision-making in pharmaceutical development. Finally, the development and application of toxicogenomic signatures for prediction of apical endpoints of regulatory concern continues to be area of intense research.

It is time to apply biclustering: a comprehensive review of biclustering applications in biological and biomedical data

Biclustering is a powerful data mining technique that allows clustering of rows and columns, simultaneously, in a matrix-format data set. It was first applied to gene expression data in 2000, aiming to identify co-expressed genes under a subset of all the conditions/samples. During the past 17 years, tens of biclustering algorithms and tools have been developed to enhance the ability to make sense out of large data sets generated in the wake of high-throughput omics technologies. These algorithms and tools have been applied to a wide variety of data types, including but not limited to, genomes, transcriptomes, exomes, epigenomes, phenomes and pharmacogenomes. However, there is still a considerable gap between biclustering methodology development and comprehensive data interpretation, mainly because of the lack of knowledge for the selection of appropriate biclustering tools and further supporting computational techniques in specific studies. Here, we first deliver a brief introduction to the existing biclustering algorithms and tools in public domain, and then systematically summarize the basic applications of biclustering for biological data and more advanced applications of biclustering for biomedical data. This review will assist researchers to effectively analyze their big data and generate valuable biological knowledge and novel insights with higher efficiency.

Unsupervised identification of disease states from high-dimensional physiological and histopathological profiles

The liver and kidney in mammals play central roles in protecting the organism from xenobiotics and are at high risk of xenobiotic-induced injury. Xenobiotic-induced tissue injury has been extensively studied from both classical histopathological and biochemical perspectives. Here, we introduce a machine-learning approach to analyze toxicological response. Unsupervised characterization of physiological and histological changes in a large toxicogenomic dataset revealed nine discrete toxin-induced disease states, some of which correspond to known pathology, but others were novel. Analysis of dynamics revealed transitions between disease states at constant toxin exposure, mostly toward decreased pathology, implying induction of tolerance. Tolerance correlated with induction of known xenobiotic defense genes and decrease of novel ferroptosis sensitivity biomarkers, suggesting ferroptosis as a druggable driver of tissue pathophysiology. Lastly, mechanism of body weight decrease, a known primary marker for xenobiotic toxicity, was investigated. Combined analysis of food consumption, body weight, and molecular biomarkers indicated that organ injury promotes cachexia by whole-body signaling through Gdf15 and Igf1, suggesting strategies for therapeutic intervention that may be broadly relevant to human disease.

Second exposure to acetaminophen overdose is associated with liver fibrosis in mice

Repeated administration of hepatotoxicants is usually accompanied by liver fibrosis. However, the difference in response as a result of repeated exposures of acetaminophen (APAP) compared to a single dose is not well-studied. Therefore, in the current study, the liver response after a second dose of APAP was investigated. Adult fasted Balb/C mice were exposed to two toxic doses of 300 mg/kg APAP, which were administered 72 h apart from each other. Subsequently, blood and liver from the treated mice were collected 24 h and 72 h after both APAP administrations. Liver transaminase, i.e. alanine amino transferase (ALT) and aspartate amino transferase (AST) levels revealed that the fulminant liver damage was reduced after the second APAP administration compared to that observed at the same time point after the first treatment. These results correlated with the necrotic areas as indicated by histological analyses. Surprisingly, Picro Sirius Red (PSR) staining showed that the accumulation of extracellular matrix after the second dose coincides with the upregulation of some fibrogenic signatures, e.g., alpha smooth muscle actin. Non-targeted liver tissue metabolic profiling indicates that most alterations occur 24 h after the first dose of APAP. However, the levels of most metabolites recover to basal values over time. This organ adaptation process is also confirmed by the upregulation of antioxidative systems like e.g. superoxide dismutase and catalase. From the results, it can be concluded that there is a different response of the liver to APAP toxic doses, if the liver has already been exposed to APAP. A necroinflammatory process followed by a liver regeneration was observed after the first APAP exposure. However, fibrogenesis through the accumulation of extracellular matrix is observed after a second challenge. Therefore, further studies are required to mechanistically understand the so called "liver memory".

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.

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.

Computational drug repositioning for rare diseases in the era of precision medicine

There are tremendous unmet needs in drug development for rare diseases. Computational drug repositioning is a promising approach and has been successfully applied to the development of treatments for diseases. However, how to utilize this knowledge and effectively conduct and implement computational drug repositioning approaches for rare disease therapies is still an open issue. Here, we focus on the means of utilizing accumulated genomic data for accelerating and facilitating drug repositioning for rare diseases. First, we summarize the current genome landscape of rare diseases. Second, we propose several promising bioinformatics approaches and pipelines for computational drug repositioning for rare diseases. Finally, we discuss recent regulatory incentives and other enablers in rare disease drug development and outline the remaining challenges.

A transcriptomics data-driven gene space accurately predicts liver cytopathology and drug-induced liver injury

Predicting unanticipated harmful effects of chemicals and drug molecules is a difficult and costly task. Here we utilize a 'big data compacting and data fusion'-concept to capture diverse adverse outcomes on cellular and organismal levels. The approach generates from transcriptomics data set a 'predictive toxicogenomics space' (PTGS) tool composed of 1,331 genes distributed over 14 overlapping cytotoxicity-related gene space components. Involving ∼2.5 × 10⁸ data points and 1,300 compounds to construct and validate the PTGS, the tool serves to: explain dose-dependent cytotoxicity effects, provide a virtual cytotoxicity probability estimate intrinsic to omics data, predict chemically-induced pathological states in liver resulting from repeated dosing of rats, and furthermore, predict human drug-induced liver injury (DILI) from hepatocyte experiments. Analysing 68 DILI-annotated drugs, the PTGS tool outperforms and complements existing tests, leading to a hereto-unseen level of DILI prediction accuracy.

Mining kidney toxicogenomic data by using gene co-expression modules

Background

Acute kidney injury (AKI) caused by drug and toxicant ingestion is a serious clinical condition associated with high mortality rates. We currently lack detailed knowledge of the underlying molecular mechanisms and biological networks associated with AKI. In this study, we carried out gene co-expression analyses using DrugMatrix—a large toxicogenomics database with gene expression data from rats exposed to diverse chemicals—and identified gene modules associated with kidney injury to probe the molecular-level details of this disease.

Results

We generated a comprehensive set of gene co-expression modules by using the Iterative Signature Algorithm and found distinct clusters of modules that shared genes and were associated with similar chemical exposure conditions. We identified two module clusters that showed specificity for kidney injury in that they 1) were activated by chemical exposures causing kidney injury, 2) were not activated by other chemical exposures, and 3) contained known AKI-relevant genes such as Havcr1, Clu, and Tff3. We used the genes in these AKI-relevant module clusters to develop a signature of 30 genes that could assess the potential of a chemical to cause kidney injury well before injury actually occurs. We integrated AKI-relevant module cluster genes with protein-protein interaction networks and identified the involvement of immunoproteasomes in AKI. To identify biological networks and processes linked to Havcr1, we determined genes within the modules that frequently co-express with Havcr1, including Cd44, Plk2, Mdm2, Hnmt, Macrod1, and Gtpbp4. We verified this procedure by showing that randomized data did not identify Havcr1 co-expression genes and that excluding up to 10 % of the data caused only minimal degradation of the gene set. Finally, by using an external dataset from a rat kidney ischemic study, we showed that the frequently co-expressed genes of Havcr1 behaved similarly in a model of non-chemically induced kidney injury.

Conclusions

Our study demonstrated that co-expression modules and co-expressed genes contain rich information for generating novel biomarker hypotheses and constructing mechanism-based molecular networks associated with kidney injury.

Electronic supplementary material

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

Characterization of Conserved Toxicogenomic Responses in Chemically Exposed Hepatocytes across Species and Platforms

BACKGROUND

Genome-wide expression profiling is increasingly being used to identify transcriptional changes induced by drugs and environmental stressors. In this context, the Toxicogenomics Project-Genomics Assisted Toxicity Evaluation system (TG-GATEs) project generated transcriptional profiles from rat liver samples and human/rat cultured primary hepatocytes exposed to more than 100 different chemicals.

OBJECTIVES

To assess the capacity of the cell culture models to recapitulate pathways induced by chemicals in vivo, we leveraged the TG-GATEs data set to compare the early transcriptional responses observed in the liver of rats treated with a large set of chemicals with those of cultured rat and human primary hepatocytes challenged with the same compounds in vitro.

METHODS

We developed a new pathway-based computational pipeline that efficiently combines gene set enrichment analysis (GSEA) using pathways from the Reactome database with biclustering to identify common modules of pathways that are modulated by several chemicals in vivo and in vitro across species.

RESULTS

We found that some chemicals induced conserved patterns of early transcriptional responses in in vitro and in vivo settings, and across human and rat genomes. These responses involved pathways of cell survival, inflammation, xenobiotic metabolism, oxidative stress, and apoptosis. Moreover, our results support the transforming growth factor beta receptor (TGF-βR) signaling pathway as a candidate biomarker associated with exposure to environmental toxicants in primary human hepatocytes.

CONCLUSIONS

Our integrative analysis of toxicogenomics data provides a comprehensive overview of biochemical perturbations affected by a large panel of chemicals. Furthermore, we show that the early toxicological response occurring in animals is recapitulated in human and rat primary hepatocyte cultures at the molecular level, indicating that these models reproduce key pathways in response to chemical stress. These findings expand our understanding and interpretation of toxicogenomics data from human hepatocytes exposed to environmental toxicants.

Toward the Replacement of Animal Experiments through the Bioinformatics-driven Analysis of 'Omics' Data from Human Cell Cultures

This paper outlines the work for which Roland Grafström and Pekka Kohonen were awarded the 2014 Lush Science Prize. The research activities of the Grafström laboratory have, for many years, covered cancer biology studies, as well as the development and application of toxicity-predictive in vitro models to determine chemical safety. Through the integration of in silico analyses of diverse types of genomics data (transcriptomic and proteomic), their efforts have proved to fit well into the recently-developed Adverse Outcome Pathway paradigm. Genomics analysis within state-of-the-art cancer biology research and Toxicology in the 21st Century concepts share many technological tools. A key category within the Three Rs paradigm is the Replacement of animals in toxicity testing with alternative methods, such as bioinformatics-driven analyses of data obtained from human cell cultures exposed to diverse toxicants. This work was recently expanded within the pan-European SEURAT-1 project (Safety Evaluation Ultimately Replacing Animal Testing), to replace repeat-dose toxicity testing with data-rich analyses of sophisticated cell culture models. The aims and objectives of the SEURAT project have been to guide the application, analysis, interpretation and storage of 'omics' technology-derived data within the service-oriented sub-project, ToxBank. Particularly addressing the Lush Science Prize focus on the relevance of toxicity pathways, a 'data warehouse' that is under continuous expansion, coupled with the development of novel data storage and management methods for toxicology, serve to address data integration across multiple 'omics' technologies. The prize winners' guiding principles and concepts for modern knowledge management of toxicological data are summarised. The translation of basic discovery results ranged from chemical-testing and material-testing data, to information relevant to human health and environmental safety.

Temporal endogenous gene expression profiles in response to lipid-mediated transfection

BACKGROUND

Design of efficient nonviral gene delivery systems is limited as a result of the rudimentary understanding of the specific molecules and processes that facilitate DNA transfer.

METHODS

Lipoplexes formed with Lipofectamine 2000 (LF2000) and plasmid-encoding green fluorescent protein (GFP) were delivered to the HEK 293T cell line. After treating cells with lipoplexes, HG-U133 Affymetrix microarrays were used to identify endogenous genes differentially expressed between treated and untreated cells (2 h exposure) or between flow-separated transfected cells (GFP+) and treated, untransfected cells (GFP-) at 8, 16 and 24 h after lipoplex treatment. Cell priming studies were conducted using pharmacologic agents to alter endogenous levels of the identified differentially expressed genes to determine effect on transfection levels.

RESULTS

Relative to untreated cells 2 h after lipoplex treatment, only downregulated genes were identified ≥ 30-fold: ALMS1, ITGB1, FCGR3A, DOCK10 and ZDDHC13. Subsequently, relative to GFP- cells, the GFP+ cell population showed at least a five-fold upregulation of RAP1A and PACSIN3 (8 h) or HSPA6 and RAP1A (16 and 24 h). Pharmacologic studies altering endogenous levels for ALMS1, FCGR3A, and DOCK10 (involved in filopodia protrusions), ITGB1 (integrin signaling), ZDDHC13 (membrane trafficking) and PACSIN3 (proteolytic shedding of membrane receptors) were able to increase or decrease transgene production.

CONCLUSIONS

RAP1A, PACSIN3 and HSPA6 may help lipoplex-treated cells overcome a transcriptional shutdown due to treatment with lipoplexes and provide new targets for investigating molecular mechanisms of transfection or for enhancing transfection through cell priming or engineering of the nonviral gene delivery system.

Temporal endogenous gene expression profiles in response to polymer-mediated transfection and profile comparison to lipid-mediated transfection

BACKGROUND

Design of efficient nonviral gene delivery systems is limited by the rudimentary understanding of specific molecules that facilitate transfection.

METHODS

Polyplexes using 25-kDa polyethylenimine (PEI) and plasmid-encoding green fluorescent protein (GFP) were delivered to HEK 293T cells. After treating cells with polyplexes, microarrays were used to identify endogenous genes differentially expressed between treated and untreated cells (2 h of exposure) or between flow-separated transfected cells (GFP+) and treated, untransfected cells (GFP-) at 8, 16 and 24 h after lipoplex treatment. Cell priming studies were conducted using pharmacologic agents to alter endogenous levels of the identified differentially expressed genes to determine effect on transfection levels. Differentially expressed genes in polyplex-mediated transfection were compared with those differentially expressed in lipoplex transfection to identify DNA carrier-dependent molecular factors.

RESULTS

Differentially expressed genes were RGS1, ARHGAP24, PDZD2, SNX24, GSN and IGF2BP1 after 2 h; RAP1A and ACTA1 after 8 h; RAP1A, WDR78 and ACTA1 after 16 h; and RAP1A, SCG5, ATF3, IREB2 and ACTA1 after 24 h. Pharmacologic studies altering endogenous levels for ARHGAP24, GSN, IGF2BP1, PDZD2 and RGS1 were able to increase or decrease transgene production. Comparing differentially expressed genes for polyplexes and lipoplexes, no common genes were identified at the 2-h time point, whereas, after the 8-h time point, RAP1A, ATF3 and HSPA6 were similarly expressed. SCG5 and PGAP1 were only upregulated in polyplex-transfected cells.

CONCLUSIONS

The identified genes and pharmacologic agents provide targets for improving transfection systems, although polyplex or lipoplex dependencies must be considered.

Meta-Analysis of Large-Scale Toxicogenomic Data Finds Neuronal Regeneration Related Protein and Cathepsin D to Be Novel Biomarkers of Drug-Induced Toxicity

Undesirable toxicity is one of the main reasons for withdrawing drugs from the market or eliminating them as candidates in clinical trials. Although numerous studies have attempted to identify biomarkers capable of predicting pharmacotoxicity, few have attempted to discover robust biomarkers that are coherent across various species and experimental settings. To identify such biomarkers, we conducted meta-analyses of massive gene expression profiles for 6,567 in vivo rat samples and 453 compounds. After applying rigorous feature reduction procedures, our analyses identified 18 genes to be related with toxicity upon comparisons of untreated versus treated and innocuous versus toxic specimens of kidney, liver and heart tissue. We then independently validated these genes in human cell lines. In doing so, we found several of these genes to be coherently regulated in both in vivo rat specimens and in human cell lines. Specifically, mRNA expression of neuronal regeneration-related protein was robustly down-regulated in both liver and kidney cells, while mRNA expression of cathepsin D was commonly up-regulated in liver cells after exposure to toxic concentrations of chemical compounds. Use of these novel toxicity biomarkers may enhance the efficiency of screening for safe lead compounds in early-phase drug development prior to animal testing.

Systematic toxicity mechanism analysis of proton pump inhibitors: an in silico study

Proton pump inhibitors (PPIs) are extensively used for the treatment of gastric acid-related disorders. PPIs appear to be well tolerated and almost have no short-term side effects. However, the clinical adverse reactions of long-term PPI usage are increasingly reported in recent years. So far, there is no study that elucidates the side effect mechanisms of PPIs comprehensively and systematically. In this study, a well-defined small molecule perturbed microarray data set of 344 compounds and 1695 samples was analyzed. With this high-throughput data set, a new index (Identity, I) was designed to identify PPI-specific differentially expressed genes. Results indicated that (1) up-regulated genes, such as RETSAT, CYP1A1, CYP1A2, and UGT, enhanced vitamin A's metabolism processes in the cellular retinol metabolism pathway; and that (2) down-regulated genes, such as C1QA, C1QC, C4BPA, C4BPB, CFI, and SERPING1, enriched in the complement and coagulation cascades pathway. In addition, strong association was observed between these PPI-specific differentially expressed genes and the reported side effects of PPIs by the gene-disease association network analysis. One potential toxicity mechanism of PPIs as suggested from this systematic PPI-specific gene expression analysis is that PPIs are enriched in acidic organelles where they are activated and inhibit V-ATPases and acid hydrolases, and consequently block the pathways of antigen presentation, the synthesis and secretion of cytokines, and complement component proteins and coagulation factors. The strategies developed in this work could be extended to studies on other drugs.

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.

Characterization of chemically induced liver injuries using gene co-expression modules

Liver injuries due to ingestion or exposure to chemicals and industrial toxicants pose a serious health risk that may be hard to assess due to a lack of non-invasive diagnostic tests. Mapping chemical injuries to organ-specific damage and clinical outcomes via biomarkers or biomarker panels will provide the foundation for highly specific and robust diagnostic tests. Here, we have used DrugMatrix, a toxicogenomics database containing organ-specific gene expression data matched to dose-dependent chemical exposures and adverse clinical pathology assessments in Sprague Dawley rats, to identify groups of co-expressed genes (modules) specific to injury endpoints in the liver. We identified 78 such gene co-expression modules associated with 25 diverse injury endpoints categorized from clinical pathology, organ weight changes, and histopathology. Using gene expression data associated with an injury condition, we showed that these modules exhibited different patterns of activation characteristic of each injury. We further showed that specific module genes mapped to 1) known biochemical pathways associated with liver injuries and 2) clinically used diagnostic tests for liver fibrosis. As such, the gene modules have characteristics of both generalized and specific toxic response pathways. Using these results, we proposed three gene signature sets characteristic of liver fibrosis, steatosis, and general liver injury based on genes from the co-expression modules. Out of all 92 identified genes, 18 (20%) genes have well-documented relationships with liver disease, whereas the rest are novel and have not previously been associated with liver disease. In conclusion, identifying gene co-expression modules associated with chemically induced liver injuries aids in generating testable hypotheses and has the potential to identify putative biomarkers of adverse health effects.

Prediction of the number of activated genes in multiple independent Cd(+2)- and As(+3)-induced malignant transformations of human urothelial cells (UROtsa)

BACKGROUND

Many toxic environmental agents such as cadmium and arsenic are found to be epidemiologically linked to increasing rates of cancers. In vitro studies show that these toxic agents induced malignant transformation in human cells. It is not clear whether such malignant transformation induced by a single toxic agent is driven by a common set of genes. Although the advancement of high-throughput technology has facilitated the profiling of global gene expression, it remains a question whether the sample size is sufficient to identify this common driver gene set.

RESULTS

We have developed a statistical method, SOFLR, to predict the number of common activated genes using a limited number of microarray samples. We conducted two case studies, cadmium and arsenic transformed human urothelial cells. Our method is able to precisely predict the number of stably induced and repressed genes and the number of samples to identify the common activated genes. The number of independent transformed isolates required for identifying the common activated genes is also estimated. The simulation study further validated our method and identified the important parameters in this analysis.

CONCLUSIONS

Our method predicts the degree of similarity and diversity during cell malignant transformation by a single toxic agent. The results of our case studies imply the existence of common driver and passenger gene sets in toxin-induced transformation. Using a pilot study with small sample size, this method also helps microarray experimental design by determining the number of samples required to identify the common activated gene set.

Beyond modules and hubs: the potential of gene coexpression networks for investigating molecular mechanisms of complex brain disorders

In a research environment dominated by reductionist approaches to brain disease mechanisms, gene network analysis provides a complementary framework in which to tackle the complex dysregulations that occur in neuropsychiatric and other neurological disorders. Gene-gene expression correlations are a common source of molecular networks because they can be extracted from high-dimensional disease data and encapsulate the activity of multiple regulatory systems. However, the analysis of gene coexpression patterns is often treated as a mechanistic black box, in which looming 'hub genes' direct cellular networks, and where other features are obscured. By examining the biophysical bases of coexpression and gene regulatory changes that occur in disease, recent studies suggest it is possible to use coexpression networks as a multi-omic screening procedure to generate novel hypotheses for disease mechanisms. Because technical processing steps can affect the outcome and interpretation of coexpression networks, we examine the assumptions and alternatives to common patterns of coexpression analysis and discuss additional topics such as acceptable datasets for coexpression analysis, the robust identification of modules, disease-related prioritization of genes and molecular systems and network meta-analysis. To accelerate coexpression research beyond modules and hubs, we highlight some emerging directions for coexpression network research that are especially relevant to complex brain disease, including the centrality-lethality relationship, integration with machine learning approaches and network pharmacology.

Using molecular features of xenobiotics to predict hepatic gene expression response

Despite recent advances in molecular medicine and rational drug design, many drugs still fail because toxic effects arise at the cellular and tissue level. In order to better understand these effects, cellular assays can generate high-throughput measurements of gene expression changes induced by small molecules. However, our understanding of how the chemical features of small molecules influence gene expression is very limited. Therefore, we investigated the extent to which chemical features of small molecules can reliably be associated with significant changes in gene expression. Specifically, we analyzed the gene expression response of rat liver cells to 170 different drugs and searched for genes whose expression could be related to chemical features alone. Surprisingly, we can predict the up-regulation of 87 genes (increased expression of at least 1.5 times compared to controls). We show an average cross-validation predictive area under the receiver operating characteristic curve (AUROC) of 0.7 or greater for each of these 87 genes. We applied our method to an external data set of rat liver gene expression response to a novel drug and achieved an AUROC of 0.7. We also validated our approach by predicting up-regulation of Cytochrome P450 1A2 (CYP1A2) in three drugs known to induce CYP1A2 that were not in our data set. Finally, a detailed analysis of the CYP1A2 predictor allowed us to identify which fragments made significant contributions to the predictive scores.

Characterization of drug-induced transcriptional modules: towards drug repositioning and functional understanding

Drug-induced transcriptional modules (biclusters) were identified and annotated in three human cell lines and rat liver. These were used to assess conservation across systems and to infer and experimentally validate novel drug effects and gene functions.

Biclustering of drug-induced gene expression profiles resulted in modules of drugs and genes, which were enriched in both drug and gene annotations.

Identifying drug-induced transcriptional modules separately in three human cell lines and rat liver allows assessment of their conservation across model systems. About 70% of modules are conserved across cell lines, a lower bound of 15% was estimated for their conservation across organisms, and between the in vitro and in vivo systems.

Drug-induced transcriptional modules can predict novel gene functions. A conserved module associated with (chole)sterol metabolism revealed novel regulators of cellular cholesterol homeostasis; 10 of them were validated in functional imaging assays.

Analysis of drugs clustered into modules can give new insights into their mechanisms of action and provide leads for drug repositioning. We predicted and experimentally validated novel cell cycle inhibitors and modulators of PPARγ, estrogen and adrenergic receptors, with potential for developing new therapies against diabetes and cancer.

In pharmacology, it is crucial to understand the complex biological responses that drugs elicit in the human organism and how well they can be inferred from model organisms. We therefore identified a large set of drug-induced transcriptional modules from genome-wide microarray data of drug-treated human cell lines and rat liver, and first characterized their conservation. Over 70% of these modules were common for multiple cell lines and 15% were conserved between the human in vitro and the rat in vivo system. We then illustrate the utility of conserved and cell-type-specific drug-induced modules by predicting and experimentally validating (i) gene functions, e.g., 10 novel regulators of cellular cholesterol homeostasis and (ii) new mechanisms of action for existing drugs, thereby providing a starting point for drug repositioning, e.g., novel cell cycle inhibitors and new modulators of α-adrenergic receptor, peroxisome proliferator-activated receptor and estrogen receptor. Taken together, the identified modules reveal the conservation of transcriptional responses towards drugs across cell types and organisms, and improve our understanding of both the molecular basis of drug action and human biology.

Integrative chemical-biological read-across approach for chemical hazard classification

Traditional read-across approaches typically rely on the chemical similarity principle to predict chemical toxicity; however, the accuracy of such predictions is often inadequate due to the underlying complex mechanisms of toxicity. Here, we report on the development of a hazard classification and visualization method that draws upon both chemical structural similarity and comparisons of biological responses to chemicals measured in multiple short-term assays ("biological" similarity). The Chemical-Biological Read-Across (CBRA) approach infers each compound's toxicity from both chemical and biological analogues whose similarities are determined by the Tanimoto coefficient. Classification accuracy of CBRA was compared to that of classical RA and other methods using chemical descriptors alone or in combination with biological data. Different types of adverse effects (hepatotoxicity, hepatocarcinogenicity, mutagenicity, and acute lethality) were classified using several biological data types (gene expression profiling and cytotoxicity screening). CBRA-based hazard classification exhibited consistently high external classification accuracy and applicability to diverse chemicals. Transparency of the CBRA approach is aided by the use of radial plots that show the relative contribution of analogous chemical and biological neighbors. Identification of both chemical and biological features that give rise to the high accuracy of CBRA-based toxicity prediction facilitates mechanistic interpretation of the models.

Integrative genomics strategies to elucidate the complexity of drug response

Pharmacogenomic investigation from both genome-wide association studies and experiments focused on candidate loci involved in drug mechanism and metabolism has yielded a substantial and increasing list of robust genetic effects on drug therapy in humans. At the same time, reasonably comprehensive molecular data such as gene expression, proteomic and metabolomic data are now available for collections of hundreds to thousands of individuals. If these data are structured in a statistically robust and computationally tractable way, such as a network model, they can aid in the analysis of new pharmacogenomics studies by suggesting novel hypotheses for the regulation of genes involved in drug metabolism and response. Similarly, hypotheses taken from these same models can direct genome-wide association studies by focusing the genome-wide association studies analysis on a number of specific hypotheses informed by the relationships customarily seen between a gene's expression or protein activity and genetic variation at a particular locus. Network models based on other sorts of systematic biological data such as cell-based surveys of drug effect on gene expression and mining of literature and electronic medical records for associations between clinical and molecular phenotypes also promise similar utility. Although surely primitive in comparison with what will be developed, these model-based approaches to leveraging the increasing volume of data generated in the course of patient care and medical research nevertheless suggest a huge opportunity to improve our understanding of biological systems involved in pharmacogenomics and apply them to questions of medical relevance.

In silico models for drug-induced liver injury--current status

INTRODUCTION

Drug-induced liver injury (DILI) is one of the most important reasons for drug attrition at both pre-approval and post-approval stages. Therefore, it is crucial to develop methods that will detect potential hepatotoxicity among drug candidates as early and quickly as possible. However, the complexity of hepatotoxicity endpoint makes it very difficult to predict. In addition, there is still a lack of sensitive and specific biomarkers for DILI that consequently leads to a scarcity of reliable hepatotoxic data, which are the key to any modelling approach.

AREAS COVERED

This review explores the current status of existing in silico models predicting hepatotoxicity. Over the past decade, attempts have been made to compile hepatotoxicity data and develop in silico models, which can be used as a first-line screening of drug candidates for further testing.

EXPERT OPINION

Most of the predictive methods discussed in this review are based on the structural properties of chemicals and do not take into account genetic and environmental factors; therefore, their predictions are still uncertain. To improve the predictability of in silico models for DILI, it is essential to better understand its mechanisms as well as to develop sensitive toxicogenomics biomarkers, which show relatively good differentiation between hepatotoxins and non-hepatotoxins.

Maximizing biomarker discovery by minimizing gene signatures

Background

The use of gene signatures can potentially be of considerable value in the field of clinical diagnosis. However, gene signatures defined with different methods can be quite various even when applied the same disease and the same endpoint. Previous studies have shown that the correct selection of subsets of genes from microarray data is key for the accurate classification of disease phenotypes, and a number of methods have been proposed for the purpose. However, these methods refine the subsets by only considering each single feature, and they do not confirm the association between the genes identified in each gene signature and the phenotype of the disease. We proposed an innovative new method termed Minimize Feature's Size (MFS) based on multiple level similarity analyses and association between the genes and disease for breast cancer endpoints by comparing classifier models generated from the second phase of MicroArray Quality Control (MAQC-II), trying to develop effective meta-analysis strategies to transform the MAQC-II signatures into a robust and reliable set of biomarker for clinical applications.

Results

We analyzed the similarity of the multiple gene signatures in an endpoint and between the two endpoints of breast cancer at probe and gene levels, the results indicate that disease-related genes can be preferably selected as the components of gene signature, and that the gene signatures for the two endpoints could be interchangeable. The minimized signatures were built at probe level by using MFS for each endpoint. By applying the approach, we generated a much smaller set of gene signature with the similar predictive power compared with those gene signatures from MAQC-II.

Conclusions

Our results indicate that gene signatures of both large and small sizes could perform equally well in clinical applications. Besides, consistency and biological significances can be detected among different gene signatures, reflecting the studying endpoints. New classifiers built with MFS exhibit improved performance with both internal and external validation, suggesting that MFS method generally reduces redundancies for features within gene signatures and improves the performance of the model. Consequently, our strategy will be beneficial for the microarray-based clinical applications.

Drug discovery in the age of systems biology: the rise of computational approaches for data integration

The increased availability of large-scale open-access resources on bioactivities of small molecules has a significant impact on pharmacology facilitated mainly by computational approaches that digest the vast amounts of data. We discuss here how computational data integration enables systemic views on a drug's action and allows to tackle complex problems such as the large-scale prediction of drug targets, drug repurposing, the molecular mechanisms, cellular responses or side effects. We particularly focus on computational methods that leverage various cell-based transcriptional, proteomic and phenotypic profiles of drug response in order to gain a systemic view of drug action at the molecular, cellular and whole-organism scale.

Predictive Power Estimation Algorithm (PPEA)--a new algorithm to reduce overfitting for genomic biomarker discovery

Toxicogenomics promises to aid in predicting adverse effects, understanding the mechanisms of drug action or toxicity, and uncovering unexpected or secondary pharmacology. However, modeling adverse effects using high dimensional and high noise genomic data is prone to over-fitting. Models constructed from such data sets often consist of a large number of genes with no obvious functional relevance to the biological effect the model intends to predict that can make it challenging to interpret the modeling results. To address these issues, we developed a novel algorithm, Predictive Power Estimation Algorithm (PPEA), which estimates the predictive power of each individual transcript through an iterative two-way bootstrapping procedure. By repeatedly enforcing that the sample number is larger than the transcript number, in each iteration of modeling and testing, PPEA reduces the potential risk of overfitting. We show with three different cases studies that: (1) PPEA can quickly derive a reliable rank order of predictive power of individual transcripts in a relatively small number of iterations, (2) the top ranked transcripts tend to be functionally related to the phenotype they are intended to predict, (3) using only the most predictive top ranked transcripts greatly facilitates development of multiplex assay such as qRT-PCR as a biomarker, and (4) more importantly, we were able to demonstrate that a small number of genes identified from the top-ranked transcripts are highly predictive of phenotype as their expression changes distinguished adverse from nonadverse effects of compounds in completely independent tests. Thus, we believe that the PPEA model effectively addresses the over-fitting problem and can be used to facilitate genomic biomarker discovery for predictive toxicology and drug responses.

Comparative analysis of predictive models for nongenotoxic hepatocarcinogenicity using both toxicogenomics and quantitative structure-activity relationships

The primary testing strategy to identify nongenotoxic carcinogens largely relies on the 2-year rodent bioassay, which is time-consuming and labor-intensive. There is an increasing effort to develop alternative approaches to prioritize the chemicals for, supplement, or even replace the cancer bioassay. In silico approaches based on quantitative structure-activity relationships (QSAR) are rapid and inexpensive and thus have been investigated for such purposes. A slightly more expensive approach based on short-term animal studies with toxicogenomics (TGx) represents another attractive option for this application. Thus, the primary questions are how much better predictive performance using short-term TGx models can be achieved compared to that of QSAR models, and what length of exposure is sufficient for high quality prediction based on TGx. In this study, we developed predictive models for rodent liver carcinogenicity using gene expression data generated from short-term animal models at different time points and QSAR. The study was focused on the prediction of nongenotoxic carcinogenicity since the genotoxic chemicals can be inexpensively removed from further development using various in vitro assays individually or in combination. We identified 62 chemicals whose hepatocarcinogenic potential was available from the National Center for Toxicological Research liver cancer database (NCTRlcdb). The gene expression profiles of liver tissue obtained from rats treated with these chemicals at different time points (1 day, 3 days, and 5 days) are available from the Gene Expression Omnibus (GEO) database. Both TGx and QSAR models were developed on the basis of the same set of chemicals using the same modeling approach, a nearest-centroid method with a minimum redundancy and maximum relevancy-based feature selection with performance assessed using compound-based 5-fold cross-validation. We found that the TGx models outperformed QSAR in every aspect of modeling. For example, the TGx models' predictive accuracy (0.77, 0.77, and 0.82 for the 1-day, 3-day, and 5-day models, respectively) was much higher for an independent validation set than that of a QSAR model (0.55). Permutation tests confirmed the statistical significance of the model's prediction performance. The study concluded that a short-term 5-day TGx animal model holds the potential to predict nongenotoxic hepatocarcinogenicity.

Functional analysis of multiple genomic signatures demonstrates that classification algorithms choose phenotype-related genes

Gene expression signatures of toxicity and clinical response benefit both safety assessment and clinical practice; however, difficulties in connecting signature genes with the predicted end points have limited their application. The Microarray Quality Control Consortium II (MAQCII) project generated 262 signatures for ten clinical and three toxicological end points from six gene expression data sets, an unprecedented collection of diverse signatures that has permitted a wide-ranging analysis on the nature of such predictive models. A comprehensive analysis of the genes of these signatures and their nonredundant unions using ontology enrichment, biological network building and interactome connectivity analyses demonstrated the link between gene signatures and the biological basis of their predictive power. Different signatures for a given end point were more similar at the level of biological properties and transcriptional control than at the gene level. Signatures tended to be enriched in function and pathway in an end point and model-specific manner, and showed a topological bias for incoming interactions. Importantly, the level of biological similarity between different signatures for a given end point correlated positively with the accuracy of the signature predictions. These findings will aid the understanding, and application of predictive genomic signatures, and support their broader application in predictive medicine.

Application of toxicogenomics in hepatic systems toxicology for risk assessment: acetaminophen as a case study

Hepatic systems toxicology is the integrative analysis of toxicogenomic technologies, e.g., transcriptomics, proteomics, and metabolomics, in combination with traditional toxicology measures to improve the understanding of mechanisms of hepatotoxic action. Hepatic toxicology studies that have employed toxicogenomic technologies to date have already provided a proof of principle for the value of hepatic systems toxicology in hazard identification. In the present review, acetaminophen is used as a model compound to discuss the application of toxicogenomics in hepatic systems toxicology for its potential role in the risk assessment process, to progress from hazard identification towards hazard characterization. The toxicogenomics-based parallelogram is used to identify current achievements and limitations of acetaminophen toxicogenomic in vivo and in vitro studies for in vitro-to-in vivo and interspecies comparisons, with the ultimate aim to extrapolate animal studies to humans in vivo. This article provides a model for comparison of more species and more in vitro models enhancing the robustness of common toxicogenomic responses and their relevance to human risk assessment. To progress to quantitative dose-response analysis needed for hazard characterization, in hepatic systems toxicology studies, generation of toxicogenomic data of multiple doses/concentrations and time points is required. Newly developed bioinformatics tools for quantitative analysis of toxicogenomic data can aid in the elucidation of dose-responsive effects. The challenge herein is to assess which toxicogenomic responses are relevant for induction of the apical effect and whether perturbations are sufficient for the induction of downstream events, eventually causing toxicity.

Biotransformation pathway maps in WikiPathways enable direct visualization of drug metabolism related expression changes

In recent decades, our knowledge of the genetics and functional genomics of drug-metabolizing enzymes has increased and a wealth of data on drug-related 'omics' has become available. Despite the availability of large amounts of biological information on xenobiotic biotransformation, the number of available biotransformation pathway maps that can easily be used for visualization of multiple omics data is limited. Here, we created integrated biotransformation pathway maps suitable for multiple omics analysis using PathVisio. The ease of visualizing data on these maps was demonstrated by using published microarray data from human hepatocyte-like cell models, exemplifying - where a sufficient capacity for metabolizing chemicals is a prerequisite for a suited model - how the biotransformation pathway maps can be used for model selection.

Tissue-level modeling of xenobiotic metabolism in liver: An emerging tool for enabling clinical translational research

This review summarizes some of the recent developments and identifies critical challenges associated with in vitro and in silico representations of the liver and assesses the translational potential of these models in the quest of rationalizing the process of evaluating drug efficacy and toxicity. It discusses a wide range of research efforts that have produced, during recent years, quantitative descriptions and conceptual as well as computational models of hepatic processes such as biotransport and biotransformation, intra- and intercellular signal transduction, detoxification, etc. The above mentioned research efforts cover multiple scales of biological organization, from molecule-molecule interactions to reaction network and cellular and histological dynamics, and have resulted in a rapidly evolving knowledge base for a "systems biology of the liver." Virtual organ/organism formulations represent integrative implementations of particular elements of this knowledge base, usually oriented toward the study of specific biological endpoints, and provide frameworks for translating the systems biology concepts into computational tools for quantitative prediction of responses to stressors and hypothesis generation for experimental design.

Biomarkers for drug-induced liver injury

Of the estimated 10,000 documented human drugs, more than 1000 have been associated with drug-induced liver injury (DILI), although causality has not always been established clearly. Numerous biomarkers for DILI have been explored, but less than ten are adopted or qualified as valid by the US FDA. The biomarkers for DILI are individual or a panel of proteins, nucleic acids or metabolites from various sources, such as the liver, blood and urine. While most DILI biomarkers are drug independent, some possibly 'drug-specific' DILIs have been explored, but specificity and sensitivity of both types need to be improved for the diagnosis of DILI during drug development and in clinical practice. Novel approaches for DILI biomarkers have been actively investigated recently, but produced mainly animal-based biomarkers, which are possibly useful for drug development, but are not suitable or have not been validated for clinical applications. This review summarizes the current practice and future perspectives for DILI biomarkers.

Modelling the genesis and treatment of cancer: the potential role of physiologically based pharmacodynamics

Physiologically based modelling of pharmacodynamics/toxicodynamics requires an a priori knowledge on the underlying mechanisms causing toxicity or causing the disease. In the context of cancer, the objective of the expert meeting was to discuss the molecular understanding of the disease, modelling approaches used so far to describe the process, preclinical models of cancer treatment and to evaluate modelling approaches developed based on improved knowledge. Molecular events in cancerogenesis can be detected using 'omics' technology, a tool applied in experimental carcinogenesis, but also for diagnostics and prognosis. The molecular understanding forms the basis for new drugs, for example targeting protein kinases specifically expressed in cancer. At present, empirical preclinical models of tumour growth are in great use as the development of physiological models is cost and resource intensive. Although a major challenge in PKPD modelling in oncology patients is the complexity of the system, based in part on preclinical models, successful models have been constructed describing the mechanism of action and providing a tool to establish levels of biomarker associated with efficacy and assisting in defining biologically effective dose range selection for first dose in man. To follow the concentration in the tumour compartment enables to link kinetics and dynamics. In order to obtain a reliable model of tumour growth dynamics and drug effects, specific aspects of the modelling of the concentration-effect relationship in cancer treatment that need to be accounted for include: the physiological/circadian rhythms of the cell cycle; the treatment with combinations and the need to optimally choose appropriate combinations of the multiple agents to study; and the schedule dependence of the response in the clinical situation.

Prediction of pharmacological and xenobiotic responses to drugs based on time course gene expression profiles

More and more people are concerned by the risk of unexpected side effects observed in the later steps of the development of new drugs, either in late clinical development or after marketing approval. In order to reduce the risk of the side effects, it is important to look out for the possible xenobiotic responses at an early stage. We attempt such an effort through a prediction by assuming that similarities in microarray profiles indicate shared mechanisms of action and/or toxicological responses among the chemicals being compared. A large time course microarray database derived from livers of compound-treated rats with thirty-four distinct pharmacological and toxicological responses were studied. The mRMR (Minimum-Redundancy-Maximum-Relevance) method and IFS (Incremental Feature Selection) were used to select a compact feature set (141 features) for the reduction of feature dimension and improvement of prediction performance. With these 141 features, the Leave-one-out cross-validation prediction accuracy of first order response using NNA (Nearest Neighbor Algorithm) was 63.9%. Our method can be used for pharmacological and xenobiotic responses prediction of new compounds and accelerate drug development.

Toxicogenomic biomarkers for liver toxicity

Toxicogenomics (TGx) is a widely used technique in the preclinical stage of drug development to investigate the molecular mechanisms of toxicity. A number of candidate TGx biomarkers have now been identified and are utilized for both assessing and predicting toxicities. Further accumulation of novel TGx biomarkers will lead to more efficient, appropriate and cost effective drug risk assessment, reinforcing the paradigm of the conventional toxicology system with a more profound understanding of the molecular mechanisms of drug-induced toxicity. In this paper, we overview some practical strategies as well as obstacles for identifying and utilizing TGx biomarkers based on microarray analysis. Since clinical hepatotoxicity is one of the major causes of drug development attrition, the liver has been the best documented target organ for TGx studies to date, and we therefore focused on information from liver TGx studies. In this review, we summarize the current resources in the literature in regard to TGx studies of the liver, from which toxicologists could extract potential TGx biomarker gene sets for better hepatotoxicity risk assessment.

Data-driven methods to discover molecular determinants of serious adverse drug events

The dangers of serious adverse drug reactions (SADRs) are well known to clinicians, pharmacologists, and the lay public. Efforts to elucidate the molecular mechanisms behind SADRs have made significant progress through genetics and gene expression measurements. However, as the field of pharmacology adopts the same novel higher-density measurement modalities that have proven successful in other areas of biology, one wonders whether there can be more ways to benefit from the explosion of data created by these tools. The development of analytic tools and algorithms to interpret these biological data to create tools for medicine is central to the field of translational bioinformatics. In this review we introduce some of the types of SADR predictors that are required, and we discuss several databases that are publicly available for the study of SADRs, ranging from clinical to molecular measurements. We also describe recent examples of how bioinformatics methods coupled with data repositories can advance the science of SADRs.

The systematic production of cells for cell therapies

Stem cells have emerged as the starting material of choice for bioprocesses to produce cells and tissues to treat degenerative, genetic, and immunological disease. Translating the biological properties and potential of stem cells into therapies will require overcoming significant cell-manufacturing and regulatory challenges. Bioprocess engineering fundamentals, including bioreactor design and process control, need to be combined with cellular systems biology principles to guide the development of next-generation technologies capable of producing cell-based products in a safe, robust, and cost-effective manner. The step-wise implementation of these bioengineering strategies will enhance cell therapy product quality and safety, expediting clinical development.