Circulating microRNA Profiles Identify a Patient Subgroup with High Inflammation and Severe Symptoms in Schizophrenia Experiencing Acute Psychosis.

Schizophrenia is a complex and heterogenous psychiatric disorder. This study aimed to demonstrate the potential of circulating microRNAs (miRNAs) as a clinical biomarker to stratify schizophrenia patients and to enhance understandings of their heterogenous pathophysiology. We measured levels of 179 miRNA and 378 proteins in plasma samples of schizophrenia patients experiencing acute psychosis and obtained their Positive and Negative Syndrome Scale (PANSS) scores. The plasma miRNA profile revealed three subgroups of schizophrenia patients, where one subgroup tended to have higher scores of all the PANSS subscales compared to the other subgroups. The subgroup with high PANSS scores had four distinctively downregulated miRNAs, which enriched 'Immune Response' according to miRNA set enrichment analysis and were reported to negatively regulate IL-1ß, IL-6, and TNFα. The same subgroup had 22 distinctively upregulated proteins, which enriched 'Cytokine-cytokine receptor interaction' according to protein set enrichment analysis, and all the mapped proteins were pro-inflammatory cytokines. Hence, the subgroup is inferred to have comparatively high inflammation within schizophrenia. In conclusion, miRNAs are a potential biomarker that reflects both disease symptoms and molecular pathophysiology, and identify a patient subgroup with high inflammation. These findings provide insights for the precision medicinal strategies for anti-inflammatory treatments in the high-inflammation subgroup of schizophrenia.

Activation of Naturally Occurring Lecithin:Cholesterol Acyltransferase Mutants by a Novel Activator Compound.

Lecithin:cholesterol acyltransferase (LCAT) is a unique plasma enzyme able to esterify cholesterol, and it plays an important role in HDL maturation and promotion of reverse cholesterol transport. Familial LCAT deficiency (FLD; OMIM number 245900) is a rare recessive disease that results from loss-of-function mutations in the LCAT gene and has no cure. In this study, we assessed the in vitro efficacy of a novel small-molecule LCAT activator. Cholesterol esterification rate (CER) and LCAT activity were tested in plasma from six controls and five FLD homozygous carriers of various LCAT mutations at different doses of the compound (0.1, 1, and 10 µg/ml). In control plasma, the compound significantly increased both CER (P < 0.001) and LCAT activity (P = 0.007) in a dose-dependent manner. Both CER and LCAT activity increased by 4- to 5-fold, reaching maximum activation at the dose of 1 µg/ml. Interestingly, Daiichi Sankyo compound produced an increase in CER in two of the five tested LCAT mutants (Leu372--Arg and Val309--Met), while LCAT activity increased in three LCAT mutants (Arg147--Trp, Thr274--Ile and Leu372--Arg); mutant Pro254--Ser was not activated at any of the tested doses. The present findings form the basis for personalized therapeutic interventions in FLD carriers and support the potential LCAT activation in secondary LCAT defects. SIGNIFICANCE STATEMENT: We characterized the pharmacology of a novel small-molecule LCAT activator in vitro on a subset of naturally occurring LCAT mutants. Our findings form the basis for personalized therapeutic interventions for familial LCAT deficiency carriers, who can face severe complications and for whom no cure exists.

Exploratory Analysis of Circulating miRNA Signatures in Atrial Fibrillation Patients Determining Potential Biomarkers to Support Decision-Making in Anticoagulation and Catheter Ablation.

Novel biomarkers are desired to improve risk management for patients with atrial fibrillation (AF). We measured 179 plasma miRNAs in 83 AF patients using multiplex qRT-PCR. Plasma levels of eight (i.e., hsa-miR-22-3p, hsa-miR-128-3p, hsa-miR-130a-3p, hsa-miR-140-5p, hsa-miR-143-3p, hsa-miR-148b-3p, hsa-miR-497-5p, hsa-miR-652-3p) and three (i.e., hsa-miR-144-5p, hsa-miR-192-5p, hsa-miR-194-5p) miRNAs showed positive and negative correlations with CHA2DS2-VASc scores, respectively, which also showed negative and positive correlations with catheter ablation (CA) procedure, respectively, within the follow-up observation period up to 6-month after enrollment. These 11 miRNAs were functionally associated with TGF-ß signaling and androgen signaling based on pathway enrichment analysis. Seven of possible target genes of these miRNAs, namely TGFBR1, PDGFRA, ZEB1, IGFR1, BCL2, MAPK1 and DICER1 were found to be modulated by more than four miRNAs of the eleven. Of them, TGFBR1, PDGFRA, ZEB1 and BCL2 are reported to exert pro-fibrotic functions, suggesting that dysregulations of these eleven miRNAs may reflect pro-fibrotic condition in the high-risk patients. Although highly speculative, these miRNAs may potentially serve as potential biomarkers, providing mechanistic and quantitative information for pathophysiology in daily clinical practice with AF such as possible pro-fibrotic state in left atrium, which would enhance the risk of stroke and reduce the preference for performing CA.

Circulating miRNA Signature as a Potential Biomarker for the Prediction of Analgesic Efficacy of Hydromorphone.

No practical biomarkers currently exist for the prediction of the analgesic efficacy of opioids. Previously, we reported circulating miRNA signatures differentially regulated by µ-opioid receptor (MOR) agonists in healthy subjects. We hypothesized that these miRNAs could be potential pharmacodynamic biomarkers to estimate MOR stimulation, and predict the efficacy of opioids; i.e., patients with low MOR stimulation may be more vulnerable to strengthening of the MOR signal upon hydromorphone treatment. To test this hypothesis, plasma samples were obtained from 25 patients with cancer pain prior to the initiation of hydromorphone treatment and the circulating miRNA levels were evaluated, focusing on four miRNAs (i.e., hsa-miR-423-3p, hsa-let-7a-5p, hsa-miR-26a-5p, and hsa-let-7f-5p) and four miRNAs (i.e., hsa-miR-144-3p, hsa-miR-451a, hsa-miR-215, and hsa-miR-363-3p) that were most clearly up and downregulated by hydromorphone and oxycodone. The patients were classified into two classes with putative high and low MOR signal, estimated based on the plasma miRNA signature. A significant correlation was observed between the analgesic efficacy and the putative MOR signal level, and patients with low MOR signal achieved better pain control (i.e., ΔVAS < 0) through hydromorphone. These results suggested that plasma miRNA signatures could serve as clinical biomarkers for the prediction of the analgesic efficacy of hydromorphone.

Identification of Circulating miRNAs Differentially Regulated by Opioid Treatment.

Emerging evidence demonstrates functional contributions of microRNAs (miRNAs) to µ-opioid receptor (MOR) signaling, but the information so far has been mostly limited to their intracellular regulatory mechanisms. The present study aimed to investigate changes in plasma miRNA profiles elicited by opioid treatment in blood samples collected from clinical studies. Healthy male subjects were orally administered with hydromorphone or oxycodone and blood samples were collected at a specified time after the drug treatment. A total of 179 plasma miRNAs were measured using multiplex qRT-PCR. Nine and seventeen miRNAs were commonly upregulated (let-7a-5p, miR-423-3p, miR-199a-3p, miR-146a-5p, miR-23b-3p, miR-24-3p, miR-221-3p, miR-223-3p, and miR-146b-5p) and downregulated (miR-144-3p, miR-215, miR-363-3p, etc.), respectively, following opioid treatment. The MOR signaling-associated miRNAs, namely let-7 family miRNAs (i.e., let-7d-5p, let-7f-5p, let-7c, let-7e-5p), miR-103a-3p, miR-339-3p, miR-146a-5p, miR-23b-3p, miR-23a-3p, and miR-181a-5p, were differentially expressed following drug treatment. These differentially expressed miRNAs are circulating biomarker candidates that can be used to evaluate MOR stimulation and serve as novel clinical diagnostic tools for improving clinical outcomes.

Data-intensive drug development in the information age: applications of Systems Biology/Pharmacology/Toxicology.

Pharmaceutical companies continuously face challenges to deliver new drugs with true medical value. R&D productivity of drug development projects depends on 1) the value of the drug concept and 2) data and in-depth knowledge that are used rationally to evaluate the drug concept's validity. A model-based data-intensive drug development approach is a key competitive factor used by innovative pharmaceutical companies to reduce information bias and rationally demonstrate the value of drug concepts. Owing to the accumulation of publicly available biomedical information, our understanding of the pathophysiological mechanisms of diseases has developed considerably; it is the basis for identifying the right drug target and creating a drug concept with true medical value. Our understanding of the pathophysiological mechanisms of disease animal models can also be improved; it can thus support rational extrapolation of animal experiment results to clinical settings. The Systems Biology approach, which leverages publicly available transcriptome data, is useful for these purposes. Furthermore, applying Systems Pharmacology enables dynamic simulation of drug responses, from which key research questions to be addressed in the subsequent studies can be adequately informed. Application of Systems Biology/Pharmacology to toxicology research, namely Systems Toxicology, should considerably improve the predictability of drug-induced toxicities in clinical situations that are difficult to predict from conventional preclinical toxicology studies. Systems Biology/Pharmacology/Toxicology models can be continuously improved using iterative learn-confirm processes throughout preclinical and clinical drug discovery and development processes. Successful implementation of data-intensive drug development approaches requires cultivation of an adequate R&D culture to appreciate this approach.

(+)-Usnic acid-induced myocardial toxicity in rats.

(+)-Usnic acid (UA) has been known to be a strong uncoupler, and mitochondrial and endoplasmic reticulum (ER)-related stresses are suggested to be involved in the mechanism of hepatotoxicity. However, it has not been clarified whether UA causes toxicity in other mitochondria-rich organs such as the heart. We elucidated whether UA induces cardiotoxicity and its mechanism. UA was orally administered to rats for 14 days, and laboratory and histopathological examinations were performed in conjunction with toxicogenomic analysis. As a result, there was no alteration in blood chemistry, whereas cytoplasmic rarefaction of myocardium was observed microscopically. This finding corresponded to the swollen mitochondria observed ultrastructurally. Immunohistochemically, expression of prohibitin, indicating mitochondrial imbalance, increased in the sarcoplasmic area. Toxicogenomic analysis highlighted the upregulation of gene groups consisting of oxidative stress, ER stress, and amino acid limitation. Interestingly, the number of upregulated genes was larger in the amino acid limitation-related gene group than that in other groups, implying that amino acid limitation might be one of the sources of oxidative stress, not only mitochondria and ER-originated stresses. In conclusion, the heart was manifested to be one of the target organs of UA. Mitochondrial imbalance with complex stresses may be involved in the toxic mechanism.

Developing a practical toxicogenomics data analysis system utilizing open-source software.

Comprehensive gene expression analysis has been applied to investigate the molecular mechanism of toxicity, which is generally known as toxicogenomics (TGx). When analyzing large-scale gene expression data obtained by microarray analysis, typical multivariate data analysis methods performed with commercial software such as hierarchical clustering or principal component analysis usually do not provide conclusive outputs by themselves. To best utilize the TGx data for toxicity evaluation in the drug development process, fit-for-purpose customization of the analytical algorithm with user-friendly interface and intuitive outputs are required to practically address the toxicologists' demands. However, commercial software is usually not very flexible in the customization of their functions or outputs. Owing to the recent advancement and accumulation of open-source software contributed by bioinformaticians all over the world, it becomes easier for us to develop practical and fit-for-purpose analytical software by ourselves with fairly low cost and efforts. The aim of this article is to present an example of developing an automated TGx data processing system (ATP system), which implements gene set-level analysis toxicogenomic profiling by D-score method and generates straightforward output that makes it easy to interpret the biological and toxicological significance of the TGx data. Our example will provide basic clues for readers to develop and customize their own TGx data analysis system which complements the function of existing commercial software.

Toxicogenomic investigation on rat testicular toxicity elicited by 1,3-dinitrobenzene.

Rats were treated with a single oral dose of 10, 25 and 50mg/kg of 1,3-dinitrobenzene (DNB), and the testis was subjected to a GeneChip microarray analysis. A total of 186 and 304 gene probe sets were up- and down-regulated, respectively, by the DNB treatment, where spermatocyte death and Sertoli cell vacuolation in testis and increased debris of spermatogenic cell in epididymis were noted. The expression profile for four sets of genes were investigated, whose expressions are reported to localize in specific cell types in the seminiferous epithelium, namely Sertoli cells, spermatogonia plus early spermtocytes, pachytene spermatocytes and round spermatids. The data demonstrated that pachytene spermatocyte-specific genes elicited explicit down-regulation in parallel with the progression of spermatocyte death, while other gene sets did not show characteristic expression changes. In addition, Gene Ontology analysis indicated that genes associated with cell adhesion-related genes were significantly enriched in the up-regulated genes following DNB treatment. Cell adhesion-related genes, namely Cdh2, Ctnna1, Vcl, Zyx, Itgb1, Testin, Lamc3, Pvrl2 and Gsn, showed an increase in microarray and the up-regulation of Cdh2 and Testin were confirmed by real time RT-PCR. The gene expression changes of pachytene spermatocyte-specific genes and cell adhesion-related genes were thought to reflect a decrease in the number of spermatocytes and dysfunction of Sertoli-germ cells adhesion junction, and therefore these genes would be potential genomic biomarkers for assessing DNB-type testicular toxicity.

Cerivastatin induces type-I fiber-, not type-II fiber-, predominant muscular toxicity in the young male F344 rats.

3-hydroxy-3-methylglutaryl coenzyme A reductase inhibitors (statins) are associated with adverse skeletal muscle toxicity, but the underlying mechanism remains unclear. To investigate the pathological mechanism of statin-induced myotoxicity, cerivastatin (20 ppm; corresponding to 2 mg/kg/day) was dietarily administered to young male F344 rats for 10 days, and time-course clinical observations, measurement of plasma creatine kinase activity, and light and electron microscopy of type I fiber-predominant skeletal muscle (soleus) or type II fiber-predominant skeletal muscles (extensor digitorum longus and tibialis anterior), were performed. Clinical symptoms including weakness of hind limbs, staggering gait and body weight loss, accompanied by marked plasma creatinine kinase elevation in rats fed cerivastatin at around Day 6 to 8. Interestingly, microscopic examination revealed that cerivastatin-induced muscle damages characterized by hypercontraction (opaque) and necrosis of the fibers were of particular abundance in the soleus muscle at Day 8, whereas these histological lesions in the extensor digitorum longus and tibialis anterior were negligible, even at Day 9. Prior to manifestation of muscle damage, swollen mitochondria and autophagic vacuoles in the soleus were observed as the earliest ultra structural changes at Day 6; then activated lysosomes, disarray of myofibril and dilated sarcoplasmic reticulum vesicles became ubiquitous at Day 8. These results demonstrate that cerivastatin induces type I fiber-predominant muscles injury, which is associated with mitochondrial damage, in young male F344 rats. Since the rat exhibiting type I fiber-targeted injury is a unique animal model for statin-induced myotoxicity, it will be useful for gaining insight into mechanisms of statin-induced myotoxicity.

Hepatic transcriptome and proteome responses against diethyl maleate-induced glutathione depletion in the rat.

Hepatic transcriptome and proteome responses against glutathione depletion were investigated by Affymetrix GeneChip Microarray and 2-dimensional gel electrophoresis (2D-DIGE), followed by MALDI-TOF-MS analysis and utilizing a glutathione-depleted rat model treated with diethyl maleate (DEM). Hepatic glutathione content decreased to 1.29 µmol/g liver (25.5% compared to control) after DEM treatment, and there were no apparent hepatotoxic signs estimated by blood chemistry examinations. A total of 247 and 213 annotated gene probe sets exhibited greater than twofold up- and down-regulation compared with controls, respectively. The up-regulated gene list contained a number of glutathione depletion-responsive genes reported previously, such as Trib3, Srxn1, Myc, Asns, Igfbp1, Txnrd1, or Hmox1, suggesting that these genes are robust mRNA biomarkers for evaluating hepatic glutathione depletion. In the 2D-DIGE analysis, proteins for a total of 361 spots were identified by MALDI-TOF-MS analysis. Of the identified proteins, 5 and 14 proteins showed up- and down-regulation, respectively. Some proteins exhibited differential expression in the protein level but not in the mRNA level, including L-FABP, MAWDBP, aldo-keto reductase family 1 member A1, catalase and ATP synthase subunit beta, suggesting that these proteins would be potential protein biomarkers for evaluating glutathione depletion. Moreover, up-regulation of FABP1 protein along with up-regulation of PPARα-regulated gene transcripts (i.e., Acot2 and Acot4) is indicative of PPARα activation, which may contribute to hepatocellular protection against glutathione depletion-induced oxidative stress. The up-regulation of L-FABP1 was detected by proteome data but not by transcriptome data, demonstrating the advantage of utilizing transcriptomics and proteomics combination to investigate glutathione depletion-induced molecular dynamics.

Practical application of toxicogenomics for profiling toxicant-induced biological perturbations.

A systems-level understanding of molecular perturbations is crucial for evaluating chemical-induced toxicity risks appropriately, and for this purpose comprehensive gene expression analysis or toxicogenomics investigation is highly advantageous. The recent accumulation of toxicity-associated gene sets (toxicogenomic biomarkers), enrichment in public or commercial large-scale microarray database and availability of open-source software resources facilitate our utilization of the toxicogenomic data. However, toxicologists, who are usually not experts in computational sciences, tend to be overwhelmed by the gigantic amount of data. In this paper we present practical applications of toxicogenomics by utilizing biomarker gene sets and a simple scoring method by which overall gene set-level expression changes can be evaluated efficiently. Results from the gene set-level analysis are not only an easy interpretation of toxicological significance compared with individual gene-level profiling, but also are thought to be suitable for cross-platform or cross-institutional toxicogenomics data analysis. Enrichment in toxicogenomics databases, refinements of biomarker gene sets and scoring algorithms and the development of user-friendly integrative software will lead to better evaluation of toxicant-elicited biological perturbations.

Ethylene glycol monomethyl ether-induced toxicity is mediated through the inhibition of flavoprotein dehydrogenase enzyme family.

Ethylene glycol monomethyl ether (EGME) is a widely used industrial solvent known to cause adverse effects to human and other mammals. Organs with high metabolism and rapid cell division, such as testes, are especially sensitive to its actions. In order to gain mechanistic understanding of EGME-induced toxicity, an untargeted metabolomic analysis was performed in rats. Male rats were administrated with EGME at 30 and 100 mg/kg/day. At days 1, 4, and 14, serum, urine, liver, and testes were collected for analysis. Testicular injury was observed at day 14 of the 100 mg/kg/day group only. Nearly 1900 metabolites across the four matrices were profiled using liquid chromatography-mass spectrometry/mass spectrometry and gas chromatography-mass spectrometry. Statistical analysis indicated that the most significant metabolic perturbations initiated from the early time points by EGME were the inhibition of choline oxidation, branched-chain amino acid catabolism, and fatty acid ß-oxidation pathways, leading to the accumulation of sarcosine, dimethylglycine, and various carnitine- and glycine-conjugated metabolites. Pathway mapping of these altered metabolites revealed that all the disrupted steps were catalyzed by enzymes in the primary flavoprotein dehydrogenase family, suggesting that inhibition of flavoprotein dehydrogenase-catalyzed reactions may represent the mode of action for EGME-induced toxicity. Similar urinary and serum metabolite signatures are known to be the hallmarks of multiple acyl-coenzyme A dehydrogenase deficiency in humans, a genetic disorder because of defects in primary flavoprotein dehydrogenase reactions. We postulate that disruption of key biochemical pathways utilizing flavoprotein dehydrogenases in conjugation with downstream metabolic perturbations collectively result in the EGME-induced tissue damage.

Methemoglobinemia induced by 1,2-dichloro-4-nitrobenzene in mice with a disrupted glutathione S-transferase Mu 1 gene.

A specific substrate to Mu class glutathione S-transferase (GST), 1,2-dichloro-4-nitrobenzene (DCNB), was administered to mice with a disrupted GST Mu 1 gene (Gstm1-null mice) to investigate the in vivo role of murine Gstm1 in toxicological responses to DCNB. A single oral administration of DCNB at doses of 500 and 1000 mg/kg demonstrated a marked increase in blood methemoglobin (MetHB) in Gstm1-null mice but not in wild-type mice. Therefore, Gstm1-null mice were considered to be more predisposed to methemoglobinemia induced by a single dosing of DCNB. In contrast, 14-day repeated-dose studies of DCNB at doses up to 600 mg/kg demonstrated a marked increase in blood MetHB in both wild-type and Gstm1-null mice. However, marked increases in the blood reticulocyte count, relative spleen weight, and extramedullary hematopoiesis in the spleen were observed in Gstm1-null mice compared with wild-type mice. In addition, microarray and quantitative reverse transcription-polymerase chain reaction analyses in the spleen showed exclusive up-regulation of hematopoiesis-related genes in Gstm1-null mice. These changes were considered to be adaptive responses to methemoglobinemia and attenuated the higher predisposition to methemoglobinemia observed in Gstm1-null mice in the single-dose study. In toxicokinetics monitoring, DCNB concentrations in plasma and blood cells were higher in Gstm1-null mice than those in wild-type mice, resulting from the Gstm1 disruption. In conclusion, it is suggested that the higher exposure to DCNB due to Gstm1 disruption was reflected in methemoglobinemia in the single-dose study and in adaptive responses in the 14-day repeated-dose study.

Gene set-level network analysis using a toxicogenomics database.

Toxicogenomics data sets on rat livers covering 118 compounds were subjected to inference of a gene set-level, not individual gene-level, network structure. Expression changing levels for 58 gene sets was used for network inference with a Gaussian graphical model algorithm. The established network contained reasonable relationships, such as ones between the blood glucose level and glycolysis-related genes or the blood transaminase level and cellular injury-related genes, indicating that the gene set-level network inference successfully highlighted biological pathway-level interactions. In addition, the robustness of the inferred network structure was investigated using microarray data on bromobenzene-treated rat livers, where the gene set-level activation exhibited time-dependent propagation through neighbored nodes (i.e. gene sets) on the network, indicating that the network structure was robust and comparable with an external microarray data set. Accumulating such robust gene sets with toxicity-associated subnetwork structures would lead to a better understanding of the molecular mechanisms of drug-elicited toxicities.

Gene expression profiling in rat liver treated with various hepatotoxic-compounds inducing coagulopathy.

A large-scale transcriptome database of rat liver (TG-GATEs) has been established by the Toxicogenomics Project in Japan. In the present study, we focused on 8 hepatotoxic compounds within TG-GATEs, i.e., clofibrate, omeprazole, ethionine, thioacetamide, benzbromarone, propylthiouracil, Wy-14,643 and amiodarone, which induced coagulation abnormalities. Aspirin was selected as a reference compound that directly causes coagulation abnormality, but not through liver toxicity. In blood chemical examinations, for all the coagulopathic compounds there was little elevation of aspartate aminotransferase (AST) and/or alanine aminotransferase (ALT), suggesting no severe cell death by treatment with the compounds. We extracted 344 probe sets from the data for these 8 typical drugs, which induced this phenotype at any time from 3 to 28 days of repeated administration. Principal component analysis using these probe sets clearly separated dose- and time-dependent clusters of the treated groups from their controls, except aspirin and propylthiouracil, both of which were considered to cause coagulopathy not due to their hepatotoxicity but due to their direct effects on the blood coagulation system. Reviewing the extracted genes, changes in lipid metabolism were found to be dominant. Genes related to blood coagulation were generally down-regulated by these drugs except that vitamin K epoxide reductase complex subunit 1 (Vkorc1) like 1, a paralogous gene of Vkorc1, was up-regulated. As expected, expression changes of these genes were least prominent in aspirin or propylthiouracil-treated liver. We concluded that these probe sets could be a good starting point in developing mechanism-based biomarkers for diagnosis or prognosis of hepatotoxicity-related coagulation abnormalities in the early stage of drug development.

Scoring multiple toxicological endpoints using a toxicogenomic database.

As information regarding microarray data sets and toxicogenomic biomarkers grows rapidly, the process of analyzing data and interpreting the results is increasingly complicated. To facilitate data analysis, a simple expression ratio-based scoring method called the TGP1 score was previously proposed [Kiyosawa, N., Shiwaku, K., Hirode, M., Omura, K., Uehara, T., Shimizu, T., Mizukawa, Y., Miyagishima, T., Ono, A., Nagao, T., Urushidani, T., 2006. Utilization of a one-dimensional score for surveying chemical-induced changes in expression levels of multiple biomarker gene sets using a large-scale toxicogenomics database. J. Toxicol. Sci. 31, 433-448]. Although the TGP1 score has demonstrated its efficacy for rapid comprehension of large-scale toxicogenomic data sets, inclusion of low quality gene expression data in the biomarker gene set produced flaws in the calculated score. To overcome this shortcoming, we tested a new scoring method called the differentially expressed gene score (D-score), where Detection Call as well as signal log ratios generated by MAS5 algorithm on Affymetrix GeneChip data were considered for the calculation. Four prototypical toxicants, namely acetaminophen, phenobarbital, clofibrate and acetamidofluorene, were used for detailed analysis. A toxicogenomics database (TG-GATEs) was utilized as a reference data set. The D-score successfully alleviated the effects of low quality data on the score calculation, and captured the overall direction of expression changes as well as the magnitude of expression change level of a set of genes, highlighting the affected toxicological endpoints elicited by chemical treatment. The D-score will be useful for high-throughput toxicity screening using a toxicogenomic database and biomarkers.

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.

Species-specific regulation of PXR/CAR/ER-target genes in the mouse and rat liver elicited by o, p'-DDT.

BACKGROUND: Dichlorodiphenyltrichloroethane (DDT) is a persistent estrogenic organochlorine pesticide that is a rodent hepatic tumor promoter, with inconclusive carcinogenicity in humans. We have previously reported that o, p'-DDT elicits primarily PXR/CAR-mediated activity, rather than ER-mediated hepatic responses, and suggested that CAR-mediated effects, as opposed to ER-mediated effects, may be more important in tumor promotion in the rat liver. To further characterize species-specific hepatic responses, gene expression analysis, with complementary histopathology and tissue level analyses were investigated in immature, ovariectomized C57BL/6 mice treated with 300 mg/kg o, p'-DDT, and compared to Sprague-Dawley rat data. RESULTS: Rats and mice exhibited negligible histopathology with rapid o, p'-DDT metabolism. Gene expression profiles were also similar, exhibiting PXR/CAR regulation with the characteristic induction of Cyp2b10 and Cyp3a11. However, PXR-specific target genes such as Apoa4 or Insig2 exhibited more pronounced induction compared to CAR-specific genes in the mouse. In addition, mouse Car mRNA levels decreased, possibly contributing to the preferential activation of mouse PXR. ER-regulated genes Cyp17a1 and Cyp7b1 were also induced, suggesting o, p'-DDT also elicits ER-mediated gene expression in the mouse, while ER-mediated effects were negligible in the rat, possibly due to the inhibitory effects of CAR on ER activities. In addition, o, p'-DDT induced Gadd45a, Gadd45b and Cdkn1, suggesting DNA damage may be an additional risk factor. Furthermore, elevated blood DHEA-S levels at 12 h after treatment in the mouse may also contribute to the endocrine-related effects of o, p'-DDT. CONCLUSION: Although DDT is known to cause rodent hepatic tumors, the marked species differences in PXR/CAR structure, expression patterns and ligand preference as well as significant species-specific differences in steroidogenesis, especially CYP17A1 expression and activity, confound the extrapolation of these results to humans. Nevertheless, the identification of potential modes of action as well as species-specific responses may assist in the selection and further development of more appropriate models for assessing the toxicity of DDT to humans and wildlife.

Toxicokintic and toxicodynamic analysis of clofibrate based on free drug concentrations in nagase analbuminemia rats (NAR).

Toxicokinetics (TK) is usually performed by measurement of the total drug concentrations in plasma. However, free drug concentrations in plasma are considered to correlate directly with toxicodynamics (TD). In the present study, to evaluate the applicability of TK/TD analysis based on free drug concentrations, we investigated the TK/TD of clofibrate, which binds to albumin with a higher ratio, using an albumin-deficient mutant strain, Nagase analbuminemia rats (NAR). TK, blood chemistry, histopathology, drug and fatty acid metabolizing enzymes and microarray analysis in the liver were examined after a 4-day oral administration of clofibrate. Compared to Sprague-Dawley (SD) rats, the parent strain of NAR, 4.1-fold higher AUC(0-24hr) based on free drug concentrations (3445 versus 844 microg.hr/ml) was observed in NAR when both rats showed the same level of AUC(0-24hr) based on the total drug concentrations (4436 versus 4237microg.hr/ml). Additionally, more severe hepatocellular hypertrophy, increase in aspartate transaminase (AST), alanine aminotransferase (ALT) and lactate dehydrogenase (LDH), decrease in total cholesterol (T.CHO), phospholipid (PL), triglyceride (TG), and non-esterified fatty acid (NEFA), and increase in the mRNA levels of fatty acid metabolizing enzymes (FAOS, CAT, and CPT) were observed in NAR at the same dose. These results demonstrated that NAR developed more severe toxicities and pharmacological effects than SD rats correlating with the higher AUC of the free drug concentrations. The results also suggested that TK/TD analysis based on the free drug concentration is appropriate to interpret the relationship between exposure and toxicity in cases of protein binding saturation including protein decrease or species differences on protein binding, especially when drugs showing a higher protein binding ratio are dosed.