Acetaminophen dosing of humans results in blood transcriptome and metabolome changes consistent with impaired oxidative phosphorylation.

UNLABELLED: The diagnosis and management of drug-induced liver injury (DILI) is hindered by the limited utility of traditional clinical chemistries. It has recently been shown that hepatotoxicants can produce compound-specific changes in the peripheral blood (PB) transcriptome in rodents, suggesting that the blood transcriptome might provide new biomarkers of DILI. To investigate in humans, we used DNA microarrays as well as serum metabolomic methods to characterize changes in the transcriptome and metabolome in serial PB samples obtained from six healthy adults treated with a 4-g bolus dose of acetaminophen (APAP) and from three receiving placebo. Treatment did not cause liver injury as assessed by traditional liver chemistries. However, 48 hours after exposure, treated subjects showed marked down-regulation of genes involved in oxidative phosphorylation/mitochondrial function that was not observed in the placebos (P < 1.66E-19). The magnitude of down-regulation was positively correlated with the percent of APAP converted to the reactive metabolite N-acetyl-p-benzoquinone-imide (NAPQI) (r= 0.739;P= 0.058). In addition, unbiased analysis of the serum metabolome revealed an increase in serum lactate from 24 to 72 hours postdosing in the treated subjects alone (P< 0.005). Similar PB transcriptome changes were observed in human overdose patients and rats receiving toxic doses. CONCLUSION: The single 4-g APAP dose produced a transcriptome signature in PB cells characterized by down-regulation of oxidative phosphorylation genes accompanied by increased serum lactate. Similar gene expression changes were observed in rats and several patients after consuming hepatotoxic doses of APAP. The timing of the changes and the correlation with NAPQI production are consistent with mechanisms known to underlie APAP hepatoxicity. These studies support the further exploration of the blood transcriptome for biomarkers of DILI.

Gene expression response in target organ and whole blood varies as a function of target organ injury phenotype.

This report details the standardized experimental design and the different data streams that were collected (histopathology, clinical chemistry, hematology and gene expression from the target tissue (liver) and a bio-available tissue (blood)) after treatment with eight known hepatotoxicants (at multiple time points and doses with multiple biological replicates). The results of the study demonstrate the classification of histopathological differences, likely reflecting differences in mechanisms of cell-specific toxicity, using either liver tissue or blood transcriptomic data.

Multicenter study of acetaminophen hepatotoxicity reveals the importance of biological endpoints in genomic analyses.

Gene expression profiling is a widely used technique with data from the majority of published microarray studies being publicly available. These data are being used for meta-analyses and in silico discovery; however, the comparability of toxicogenomic data generated in multiple laboratories has not been critically evaluated. Using the power of prospective multilaboratory investigations, seven centers individually conducted a common toxicogenomics experiment designed to advance understanding of molecular pathways perturbed in liver by an acute toxic dose of N-acetyl-p-aminophenol (APAP) and to uncover reproducible genomic signatures of APAP-induced toxicity. The nonhepatotoxic APAP isomer N-acetyl-m-aminophenol was used to identify gene expression changes unique to APAP. Our data show that c-Myc is induced by APAP and that c-Myc-centered interactomes are the most significant networks of proteins associated with liver injury. Furthermore, sources of error and data variability among Centers and methods to accommodate this variability were identified by coupling gene expression with extensive toxicological evaluation of the toxic responses. We show that phenotypic anchoring of gene expression data is required for biologically meaningful analysis of toxicogenomic experiments.

ATM requirement in gene expression responses to ionizing radiation in human lymphoblasts and fibroblasts.

The heritable disorder ataxia telangiectasia (AT) is caused by mutations in the AT-mutated (ATM) gene with manifestations that include predisposition to lymphoproliferative cancers and hypersensitivity to ionizing radiation (IR). We investigated gene expression changes in response to IR in human lymphoblasts and fibroblasts from seven normal and seven AT-affected individuals. Both cell types displayed ATM-dependent gene expression changes after IR, with some responses shared and some responses varying with cell type and dose. Interestingly, after 5 Gy IR, lymphoblasts displayed ATM-independent responses not seen in the fibroblasts at this dose, which likely reflect signaling through ATM-related kinases, e.g., ATR, in the absence of ATM function.

Standardizing global gene expression analysis between laboratories and across platforms.

To facilitate collaborative research efforts between multi-investigator teams using DNA microarrays, we identified sources of error and data variability between laboratories and across microarray platforms, and methods to accommodate this variability. RNA expression data were generated in seven laboratories, which compared two standard RNA samples using 12 microarray platforms. At least two standard microarray types (one spotted, one commercial) were used by all laboratories. Reproducibility for most platforms within any laboratory was typically good, but reproducibility between platforms and across laboratories was generally poor. Reproducibility between laboratories increased markedly when standardized protocols were implemented for RNA labeling, hybridization, microarray processing, data acquisition and data normalization. Reproducibility was highest when analysis was based on biological themes defined by enriched Gene Ontology (GO) categories. These findings indicate that microarray results can be comparable across multiple laboratories, especially when a common platform and set of procedures are used.

Differential gene expression profiling in whole blood during acute systemic inflammation in lipopolysaccharide-treated rats.

Microarrays have been used to evaluate the expression of thousands of genes in various tissues. However, few studies have investigated the change in gene expression profiles in one of the most easily accessible tissues, whole blood. We utilized an acute inflammation model to investigate the possibility of using a cDNA microarray to measure the gene expression profile in the cells of whole blood. Blood was collected from male Sprague-Dawley rats at 2 and 6 h after treatment with 5 mg/kg (ip) LPS. Hematology showed marked neutrophilia accompanied by lymphopenia at both time points. TNF-alpha and IL-6 levels were markedly elevated at 2 h, indicating acute inflammation, but by 6 h the levels had declined. Total RNA was isolated from whole blood and hybridized to the National Institute of Environmental Health Sciences Rat Chip v.3.0. LPS treatment caused 226 and 180 genes to be differentially expressed at 2 and 6 h, respectively. Many of the differentially expressed genes are involved in inflammation and the acute phase response, but differential expression was also noted in genes involved in the cytoskeleton, cell adhesion, oxidative respiration, and transcription. Real-time RT-PCR confirmed the differential regulation of a representative subset of genes. Principal component analysis of gene expression discriminated between the acute inflammatory response apparent at 2 h and the observed recovery underway at 6 h. These studies indicate that, in whole blood, changes in gene expression profiles can be detected that are reflective of inflammation, despite the adaptive shifts in leukocyte populations that accompany such inflammatory processes.

Gene expression profiling of rat livers reveals indicators of potential adverse effects.

This study tested the hypothesis that gene expression profiling can reveal indicators of subtle injury to the liver induced by a low dose of a substance that does not cause overt toxicity as defined by conventional criteria of toxicology (e.g., abnormal clinical chemistry and histopathology). For the purpose of this study we defined this low dose as subtoxic, i.e., a dose that elicits effects which are below the detection of conventional toxicological parameters. Acetaminophen (APAP) was selected as a model hepatotoxicant because (1) considerable information exists concerning the mechanism of APAP hepatotoxicity that can occur following high doses, (2) intoxication with APAP is the leading cause of emergency room visits involving acute liver failure within the United States, and (3) conventional clinical markers have poor predictive value. Rats treated with a single dose of 0, 50, 150, or 1500 mg/kg APAP were examined at 6, 24, or 48 h after exposure for conventional toxicological parameters and for gene expression alterations. Patterns of gene expression were found which indicated cellular energy loss as a consequence of APAP toxicity. Elements of these patterns were apparent even after exposure to subtoxic doses. With increasing dose, the magnitude of changes increased and additional members of the same biological pathways were differentially expressed. The energy loss suggested by gene expression changes was confirmed at the 1500 mg/kg dose exposure by measuring ATP levels. Only by ultrastructural examination could any indication of toxicity be identified after exposure to a subtoxic dose of APAP and that was occasional mitochondrial damage. In conclusion, this study provides evidence that supports the hypothesis that gene expression profiling may be a sensitive means of identifying indicators of potential adverse effects in the absence of the occurrence of overt toxicity.

ATM-dependent and -independent gene expression changes in response to oxidative stress, gamma irradiation, and UV irradiation.

Ataxia telangiectasia (AT) is an autosomal recessive disorder characterized by progressive cerebellar degeneration, immunodeficiencies, telangiectasias, sensitivity to ionizing radiation, and high predisposition for malignancies. The ataxia telangiectasia mutated (ATM) gene encodes a protein (ATM) with serine/threonine kinase activity. DNA-double strand breaks are known to increase its kinase activity. While cells from individuals with AT are attenuated in their G(1)-, S- and G(2)-phase cell cycle checkpoint functions in response to gamma irradiation and oxidative stress, their response to UV irradiation appears to be equivalent to that of wild-type cells. In this study, we investigated changes in gene expression in response to gamma irradiation, oxidative stress, and UV irradiation, focusing on the dependence on ATM. Doses for all three treatments were selected that resulted in roughly an equivalent induction of a G(1) checkpoint response and inhibition of progression through S phase. To investigate gene expression changes, logarithmically growing wild-type and AT dermal diploid fibroblasts were exposed to either gamma radiation (5 Gy), oxidative stress (75 micro M t-butyl-hydroperoxide), or UV radiation (7.5 J/m(2)), and RNA was harvested 6 h after treatment. Gene expression analysis was performed using the NIEHS Human ToxChip 2.0 with approximately 1900 cDNA clones representing known genes and ESTs. All three treatments resulted in distinct patterns of gene expression changes, as shown previously. ATM-dependent and ATM-independent components were detected within these patterns, as were novel indications of involvement of ATM in regulation of transcription factors such as SP1, AP1 and MTF1.

Identification of distinct and common gene expression changes after oxidative stress and gamma and ultraviolet radiation.

The human genome is exposed to many different kinds of DNA-damaging agents. While most damage is detected and repaired through complex damage recognition and repair machineries, some damage has the potential to escape these mechanisms. Unrepaired DNA damage can give rise to alterations and mutations in the genome in an individual cell, which can result in malignant transformation, especially when critical genes are deregulated. In this study, we investigated gene expression changes in response to oxidative stress, gamma (gamma) radiation, and ultraviolet (UV) radiation and their potential implications in cancer development. Doses were selected for each of the three treatments, based on their ability to cause a similar G(1) checkpoint induction and slow down in early S-phase progression, as reflected by a comparable reduction in cyclin E-associated kinase activity of at least 75% in logarithmically growing human dermal diploid fibroblasts. To investigate gene expression changes, logarithmically growing dermal diploid fibroblasts were exposed to either gamma radiation (5 Gy), oxidative stress (75 microM of tert-butyl hydroperoxide (t-butyl-OOH)), or UV radiation (UVC) (7.5 J/m(2)) and RNA was harvested 6 h after treatment. Gene expression was analyzed using the NIEHS Human ToxChip 2.0 with approximately 1901 cDNA clones representing known genes and expressed sequence tags (ESTs). We were able to identify common and distinct responses in dermal diploid fibroblasts to the three different stimuli used. Within our analysis, gene expression profiles in response to gamma radiation and oxidative stress appeared to be more similar than profiles expressed after UV radiation. Interestingly, equivalent cyclin E-associated kinase activity reduction with all the three treatments was associated with greater transcriptional changes after UV radiation than after gamma radiation and oxidative stress. While samples treated with UV radiation displayed modulations of their mitogen activated protein kinase (MAPK) pathway, gamma radiation had its major influence on cell-cycle progression in S-phase and mitosis. In addition, cell cultures from different individuals displayed significant differences in their gene expression responses to DNA damage.