Skip Navigation


ToxSci Advance Access originally published online on August 17, 2006
Toxicological Sciences 2006 94(1):206-216; doi:10.1093/toxsci/kfl078
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
94/1/206    most recent
kfl078v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (5)
Right arrowRequest Permissions
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Yin, H.-Q.
Right arrow Articles by Lee, B.-H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Yin, H.-Q.
Right arrow Articles by Lee, B.-H.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2006. Published by Oxford University Press on behalf of the Society of Toxicology. All rights reserved. For Permissions, please email: journals.permissions@oxfordjournals.org

Hepatic Gene Expression Profiling and Lipid Homeostasis in Mice Exposed to Steatogenic Drug, Tetracycline

Hu-Quan Yin*, Mingoo Kim{dagger}, Ju-Han Kim{dagger},{ddagger}, Gu Kong§, Mi-Ock Lee*, Kyung-Sun Kang, Byung-IL Yoon||, Hyung-Lae Kim||| and Byung-Hoon Lee*,1

* College of Pharmacy and Research Institute of Pharmaceutical Sciences, Seoul National University, Seoul 151-742, Republic of Korea {dagger} Seoul National University Biomedical Informatics, College of Medicine, Seoul National University, Seoul 110-799, Republic of Korea {ddagger} Human Genome Research Institute, College of Medicine, Seoul National University, Seoul 110-799, Republic of Korea § College of Medicine, Hanyang University, Seoul133-791, Republic of Korea College of Veterinary Medicine, Seoul National University, Seoul 151-742, Republic of Korea || College of Veterinary Medicine, Kangwon National University, Chuncheon 200-791, Republic of Korea ||| College of Medicine, Ewha Womans University, Seoul 158-710, Republic of Korea

1To whom correspondence should be addressed at College of Pharmacy, Seoul National University, San 56-1, Sillim-dong, Gwanak-gu, Seoul 151-742, Republic of Korea. Fax: +82-2-874-7843. E-mail: lee{at}snu.ac.kr.

Received May 30, 2006; accepted August 3, 2006


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Tetracycline is one of a group of drugs known to induce microvesicular steatosis. In the present study, we investigated the effects of tetracycline on gene expression in mouse liver, using Applied Biosystems Mouse Genome Survey Microarrays. A single oral dose of 0.1 or 1 g/kg tetracycline was administered to male ICR mice, and liver samples were obtained after 6, 24, or 72 h. Histopathological evaluation showed microvesicular steatosis in the high-dose group at 24 h. In total, 96 genes were identified as tetracycline responsive. Their level of expression differed significantly from controls (two-way analysis of variance; p < 0.05), after adjustment by the Benjamini-Hochberg multiple testing correction, and displayed a twofold or greater induction or repression. The largest groups of gene products affected by tetracycline exposure were those involved in signal transduction, nucleic acid metabolism, developmental processes, and protein metabolism. The expression of genes known to be involved in lipid metabolism was examined, using two-sample Student's t-test for each treatment group versus a corresponding control group. The overall net effects on expression of lipid metabolism genes indicated an increase in cholesterol and triglyceride biosynthesis and a decrease in ß-oxidation of fatty acids. Our data support a proposed mechanism for tetracycline-induced steatogenic hepatotoxicity that involves these processes. Moreover, we demonstrated global changes in hepatic gene expression following tetracycline exposure; many of these genes have the potential to be used as biomarkers of exposure to steatogenic hepatotoxic agents.

Key Words: tetracycline; microvesicular steatosis; toxicogenomics; lipid metabolism; microarray.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The observation of toxicity in animals or of adverse reactions in humans are compelling reasons for terminating a drug's development, as well as for postmarketing withdrawal. Among the toxic side effects seen in animal testing and in clinical trials, hepatotoxicity is the principal cause of termination for safety reasons (Ballett, 1997Go; Lumley and Walder, 1987aGo,bGo). Animals that are frequently used in toxicity tests are not fully predictive of the response in humans due to species variation in physiology, anatomy, and metabolism. The number of animals, the dose of the drug, and the duration of the experiment are also critical limiting factors. Therefore, practical needs exist for establishing new predictive toxicology methods. Toxicogenomic or in silico technologies are currently popular as fast and cost-efficient alternatives or supplements to bioassays for the identification of toxic effects at an early stage of drug development (Helma, 2005Go; Storck et al., 2002Go; Suter et al., 2004Go).

Tetracyclines are a well-known family of antibiotics that are active against a wide range of gram-positive and gram-negative bacteria. In addition, anti-inflammatory properties of tetracycline, which are independent of its antibacterial activity, have been described (Gabler et al., 1992Go). Tetracycline induces microvesicular steatosis in the liver, which is severe and even fatal in some vulnerable patients (Bhagavan et al., 1982Go; Westphal et al., 1994Go). Although formerly considered a benign and fully reversible condition, hepatic fatty infiltration is an important pathogenic factor in the development of many drug-induced liver diseases. Tetracyclines downregulate the expression of peroxisome proliferator activated receptors, which are widely used in the well-established Tet-off systems for studying inducible gene expression (Tachibana et al., 2005Go; Wu et al., 1996Go). Freneaux et al. (1988)Go demonstrated the inhibition of mitochondrial ß-oxidation of fatty acids by tetracycline in mice. Similarly, Amacher and Martin (1997)Go showed that the canine hepatocyte is susceptible to tetracycline-induced steatosis, when triglyceride (TG) accumulation was concomitant with the inhibition of mitochondrial lipid metabolism.

Large-scale gene expression analysis provides a logical approach to studying the detailed mechanisms of chemical-induced organ toxicity as well as to identifying potential biomarkers of toxicity. To identify global gene expression changes that are associated with hepatotoxicity and to better understand the molecular mechanisms underlying tetracycline-induced steatosis, we investigated the effects of tetracycline on gene expression in mouse liver across multiple doses and time points, with special emphasis on genes associated with lipid metabolism.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Animal treatments.
Specific pathogen-free male ICR mice (20–25 g) were obtained from Jung-Ang Laboratory Animal Co. (Seoul, Korea) and allowed free access to standard chow and tap water. They were kept in temperature-controlled and filter-sterilized animal quarters under a 12-h light:12-h dark cycle. The use of animals was in compliance with the guidelines established by the Animal Care Committee of this institute. A total of 27 mice (nine per group) were given a single oral dose of 0.1 g/kg (low-dose group, L), 1 g/kg (high-dose group, H), or distilled water (control group, C), and three from each group were subsequently killed at 6 (C6, L6, and H6), 24 (C24, L24, and H24), or 72 h (C72, L72, and H72). The dose was selected by preliminary experiments that induce mild hepatic steatosis at 24 h. A cross section of the left lateral lobe of the liver was collected in 10% neutral-buffered formalin for histopathology. The remaining liver tissue was collected in RNase-free tubes and snap frozen in liquid nitrogen. Frozen tissues were stored at – 70°C until processed for RNA extraction.

Clinical chemistry and histological analysis.
Serum samples were prepared from blood withdrawn by heart puncture, using standard methods. Serum levels of alanine aminotransferase (ALT), alkaline phosphatase, blood urea nitrogen, creatinine, total bilirubin, TG, and albumin were monitored by standard clinical chemistry assays on an Automated Chemistry Analyzer (Prestige 24I; Tokyo Boeki Medical System, Tokyo, Japan). Liver tissues collected in formalin were dehydrated, embedded in paraffin, sectioned at 4 µm, and stained with hematoxylin and eosin. Histopathologic examinations of the liver sections were conducted by a pathologist and were peer reviewed.

RNA isolation.
Total RNA was extracted using an Easy-Blue total RNA extraction kit (Intron Biotech, Sungnam, Korea), purified using Qiagen RNeasy Mini Kits (Qiagen, Basel, Switzerland), and examined for integrity using an Agilent 2100 Bioanalyzer (Ambion, Austin, TX).

Microarray analysis of differential gene expression.
Applied Biosystems (Foster City, CA) Mouse Genome Survey Microarrays were used to analyze differential gene expression profiles. These arrays contain 32,996 probes that cover 44,498 individual transcripts and target a complete, annotated, and fully curated set of 32,381 mouse genes from the public and Celera databases. RNA extracted from the livers of three mice was analyzed for each time point (6, 24, and 72 h) and for each treatment. Digoxigenin-UTP–labeled cRNA was generated and linearly amplified from 2 µg of total RNA using an Applied Biosystems Chemiluminescent RT-IVT Labeling Kit v2.0. Array hybridization, chemiluminescence detection, image acquisition, and analysis were performed using an Applied Biosystems Chemiluminescence Detection Kit and an Applied Biosystems 1700 Chemiluminescent Microarray Analyzer, following the manufacturer's protocol. Briefly, each microarray was first prehybridized at 55°C for 1 h in hybridization buffer with blocking reagent. Ten micrograms of labeled cRNA targets were incubated with fragmentation buffer at 60°C for 30 min and hybridized to each microarray in a 0.5-ml volume at 55°C for 16 h. Images were autogridded, and spots were spatially normalized. Chemiluminescent signals were quantified, corrected for background, and the final images and feature data were processed using the Applied Biosystems 1700 Chemiluminescent Microarray Analyzer v1.1.

Data analysis.
The quality of hybridization and overall chip performance were monitored by visual inspection of both internal quality control checks and the raw scanned data. Assay-normalized signals with quality flags lower than 100 were used for the analysis. Filtered values were imputed by a k-nearest neighbors imputation algorithm (Troyanskaya et al., 2001Go), transformed by the variance-stabilizing normalization method (Huber et al., 2002Go), and normalized by quantile normalization across the arrays (Bolstad et al., 2003Go).

The overall flow of data analysis is illustrated in Figure 1. Two-way analysis of variance (ANOVA) was applied to determine both time and dose effects and differentially expressed sets of genes. Statistical significances were adjusted by Benjamini-Hochberg multiple testing correction. Hierarchical clustering was performed with the Cluster/TreeView analytic package (Eisen, Stanford University, CA). Biological pathway and ontology-based analysis were performed by using Panther database (http://www.pantherdb.org) and ArrayX Path (http://www.snubi.org/software/ArrayXPath/). Post hoc two-sample Student's t-tests were performed. Default parameters from R statistical package and Avadis 4.3 (Strand Life Sciences, India) software applied when unspecified. For further analysis, gene sets with different dose and time effects according to two-way ANOVA were classified into time-dependent, dose-dependent, and combined effect (i.e., both time- and dose-dependent, MTD) groups. The three groups were then classified into with- and without-interaction groups according to the significance of the interaction terms under saturated linear models.


Figure 1
View larger version (28K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 1 Flow chart for data analysis. The microarray data were analyzed by two-way ANOVA for identifying global gene expression changes and by two-sample Student's t-test for investigating genes associated with lipid metabolism at specific treatment conditions.

 
Quantitative real-time polymerase chain reaction.
Total RNA was prepared using the Easy-Blue total RNA extraction kit (Intron Biotech), and single-strand cDNA was synthesized from the RNA in a reaction mixture containing random hexamers and Superscript II reverse transcriptase (Invitrogen, Carlsbad, CA). Gene-specific primers designed using Oligo 6.0 software (Molecular Biology Insights, Cascade, CO) were used (Table 1). Quantitative real-time polymerase chain reaction (qRT-PCR) amplifications were performed using Fast Start DNA Master SYBR Green I Mixture Kit (Roche Diagnostics, Mannheim, Germany) in a Light Cycler system (Roche Diagnostics) following manufacturer's protocol. To determine the specificity of amplification, melting curve analysis was applied to all final PCR products.


View this table:
[in this window]
[in a new window]

 
TABLE 1 Gene-Specific Primers Used in qRT-PCR

 

    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Clinical Chemistry and Histological Endpoints
Treatment of ICR mice with 0.1 or 1 g/kg tetracycline did not change the level of plasma ALT activity in any of the animals (Fig. 2A). In the histological analysis, no biologically significant liver damage was observed at 6 or at 72 h after treatment (data not shown). However, high-dose tetracycline at 24 h caused histological lesions including microvesicular steatosis, which is characterized by the presence of numerous small vesicles of fat that do not displace the nucleus (Fig. 2B). Other serum biomarkers tested in this study were not changed significantly by the drug treatment, indicating that the treatment condition used in this study did not induce any serious lesions.


Figure 2
View larger version (86K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 2 Biochemical and pathological analysis of the serum and liver obtained from ICR mice after oral exposure to tetracycline. (A) The levels of serum ALT were measured at each time point in vehicle (open bar), low-dose (filled bar), and high-dose (checked bar) treatment groups. Dashed line indicates the normal range of ALT level in ICR mice. Results are mean ± SD of three mice. (B) Histopathological analysis of hematoxylin-eosin–stained paraffin sections from (a) C24, (b) L24, and (c) H24 liver samples (x400).

 
Global Gene Expression Analysis Following Tetracycline Treatment
A total of 27 Applied Biosystems Mouse Genome Survey Microarrays were used to analyze global gene expression changes in mouse liver exposed to vehicle or to tetracycline (0.1 or 1 g/kg) at 6, 24, and 72 h posttreatment. Two-way ANOVA of the dose- and time-response data identified 4934 genes (from a total of 32,381 probe sets) that showed statistically significant changes in expression (p < 0.05, Supplementary data 1). Functional categorization of the data was performed using the annotation information in the PANTHER database for ontology. Tetracycline-induced changes in expression of genes involved in a variety of biological processes. The largest group of genes was those associated with signal transduction, while genes of nucleic acid metabolism, protein metabolism, and developmental processes were the next most abundant annotation groups. Tree analysis of the genes classified patterns of relative responses for all treatment conditions relative to controls (Fig. 3). The greatest distance between nodes was between the L6 and the H72 groups. Low- and high-dose groups at each time point were in the same node, and the distances between nodes increased in response to the time of treatment, indicating a stronger effect of time than dose in our experimental regimen. Further classification of the results into four groups under the terms of a saturated linear model demonstrated that 3229 of 4934 differentially expressed genes showed a time-dependent response, and only 667 genes were dose dependent. These results again show that the effect of time on differential gene expression was more significant than that of dose in our tetracycline-induced gene expression analysis.


Figure 3
View larger version (25K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 3 Hierarchical clustering analysis of differentially expressed genes in tetracycline-treated mice liver. Hierarchical clustering analysis was performed with 298 genes that were obtained by the adjustment of the data using Benjamini-Hochberg multiple test correction. Red color and green color in matrix indicate the relative gene induction and repression, respectively. Dendrograms of each group (above the matrix) and probe sets (left of the matrix) represent the overall similarities in gene expression profiles.

 
Adjustment of the data by the Benjamini-Hochberg multiple testing correction (Benjamini and Hochberg, 1995Go) reduced the number of affected genes to 298 (Supplementary data 2), and the expression of a total 96 genes was found to be changed by more than twofold (up or down) by at least one treatment condition. Functional categorization of the genes using the PANTHER database for ontology is presented in Table 2. Genes without clear annotations, such as expressed sequence tags, were excluded from the table.


View this table:
[in this window]
[in a new window]

 
TABLE 2 Genes Altered More than Twofold (two-way ANOVA; p < 0.05) Compared with Control in the Liver of Mice Given Tetracycline

 
Expression of Genes in Lipid Metabolism
We next investigated the expression of genes known to be involved in lipid metabolism. The whole data set obtained from each treatment group was analyzed by two-sample Student's t-test against that from the corresponding vehicle control group. In general, genes were more likely to be induced than repressed by tetracycline treatment. The greatest number of changes in gene expression occurred at 24 h, when signals from 2429 genes in the high-dose (H24) group and 1009 genes in the low-dose (L24) group were changed. Protein tyrosine phosphatase receptor type C and the RIKEN cDNA 8430438D04 were upregulated at all three time points. Lipid metabolism–related genes were present in the high- and low-dose gene groups (141 and 62, respectively; Tables 3 and 4). As shown in Figure 4, the overall net effects of altered gene expression associated with lipid metabolism were to increase cholesterol and TG biosynthesis and to decrease ß-oxidation of fatty acid.


View this table:
[in this window]
[in a new window]

 
TABLE 3 Significantly Upregulated Genes (Student's t-test; p < 0.05) in H24 Compared with C24, Which Are Involved in Lipid Metabolism (fold change ≥ 1.5)

 

View this table:
[in this window]
[in a new window]

 
TABLE 4 Significantly Downregulated Genes (Student's t-Test; p < 0.05) in H24 Compared with C24, Which Are Involved in Lipid Metabolism (fold change ≥ 1.5)

 

Figure 4
View larger version (14K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 4 Pathways for major lipid metabolism (TG [A], cholesterol biosynthesis [B], and ß-oxidation of fatty acid [C]) and the relative induction or repression of the genes, being statistically different from vehicle control based on a two-sample Student's t-test. To make data more intelligible, they are presented as fold induction in natural numbers (not in log2 scale). (**p < 0.05 and fold induction or repression ≥ 1.5; *p < 0.05 and fold induction or repression ≥ 1.25).

 
Validation of Microarray Data
Microarray results were confirmed for a subset of genes using qRT-PCR. In general, the results obtained from microarray experiments correlated well with qRT-PCR. For simple validation, 12 genes were verified by qRT-PCR, all of which displayed expression patterns comparable with the microarray data (Table 5). Because the steatosis was induced only in H24 group and was apparently reversed between 24 and 72 h, we analyzed the data by Tukey's post hoc test to find the genes that were upregulated in H24 and downregulated in H72. There was good agreement in the magnitude of the fold change when comparing microarray and qRT-PCR data, except for Myt1L that showed discrepancy between the two methods. These genes represent steatosis-associated injury and recovery specific genes (Fig. 5).


View this table:
[in this window]
[in a new window]

 
TABLE 5 qRT-PCR Validation of Selected Genes from Microarray Data in L24 and H24 Groups

 

Figure 5
View larger version (9K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 5 Real-time PCR validation of the genes that were upregulated in H24 and downregulated in H72, which represent directional changes in the steatosis-associated injury and recovery specific genes. Data represent fold change on log2 scale compared with the corresponding control group.

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Tetracycline is a widely used broad-spectrum antibiotic that is known to induce microvesicular steatosis in the liver. Although the hepatotoxic effect of tetracycline is evident, and is severe and even fatal in some vulnerable patients, a systematic approach to identifying a toxicity mechanism has not yet been applied. In the present study, we investigated the temporal and dose-dependent effects of acute tetracycline on hepatic gene expression, using high-density oligonucleotide microarrays, with special emphasis on the genes associated with lipid metabolism. The goals of this study were (1) to investigate the effects of tetracycline on global gene expression, in order to develop biomarkers and (2) to identify genes that are involved in the steatogenic action of tetracycline, in order to elucidate its mode of toxicity.

Tetracycline-Induced Global Gene Expression Changes
Accumulating evidence has demonstrated that microarray technology for global gene expression profiling is a powerful tool in toxicological research, especially in the field of predictive toxicogenomics. A critical hypothesis in this research area is that the resulting transcript profile can provide both reliable advance information on possible toxic outcomes and mechanistic insight into the toxicity itself. To use short-term toxicogenomic data in predicting results of longer duration studies, a database of compounds with known chronic toxicity must be developed. Although many technical problems limit the usefulness of the concept at present, pattern recognition experiments may hold considerable promise (Hamadeh et al., 2002Go). For this purpose, the Korean Toxicogenomic Research Consortium has embarked on a project to construct a toxicogenomic database of known hepatotoxicants, including tetracycline as a steatogenic drug.

When the microarray data from our nine experimental groups were analyzed statistically using two-way ANOVA, about 15% (4934 of 32,381) of the genes were significantly affected by tetracycline—a relatively large subset. To determine the dominant factor influencing the gene expression pattern, we next classified and adjusted the data using gene tree analysis, a saturated linear model, and the Benjamini-Hochberg multiple testing correction. When thousands of hypotheses are tested simultaneously, the chance of detecting false positives increases, and multiple testing methods should be applied to assess the statistical significance of findings. To control for the false discovery rate (FDR), defined as the expected proportion of false positives among the genes with significant expression changes, we applied the Benjamini-Hochberg multiple testing correction method. In the present study, the estimated number of false positives was presumed to be about 15 out of 298 genes (FDR < 0.05). The combined results of the three methods of analysis indicate unequivocally that the effect of time on differential gene expression is more significant than that of dose in tetracycline-induced gene expression.

Tetracycline-Induced Changes in Lipid Metabolism
The typical phenotypic change induced by tetracycline is fatty liver. We analyzed data from each treatment group against the corresponding vehicle control. Only at 24 h was the number of genes induced or suppressed directly proportional to the tetracycline dose. At that time point, of 1009 and 2429 genes that were changed by more than 1.5-fold in the low- and high-dose groups, respectively, 62 and 141 genes were involved in lipid metabolism, including mitochondrial and peroxisomal ß-oxidation of free fatty acid, TG synthesis, and cholesterol synthesis.

Mechanisms so far reported by which tetracycline induces steatosis include the inhibition both of ß-oxidation of free fatty acids and of lipoprotein secretion from the liver (Freneaux et al., 1988Go; Letteron et al., 2003Go). However, the effect of tetracycline on TG and cholesterol synthesis is not known. We observed altered patterns of expression for genes involved in these synthetic pathways. Glycerol-3-phosphate dehydrogenase (Gpd2) catalyzes the reduction of dihydroxyacetone phosphate to glycerol-3-phosphate, which is in turn acylated by glycerol-3-phosphate-acyltransferase and acylglycerol-phosphate-acyltransferase (Agpat2) to give phosphatidic acid. Phosphatidic acid is hydrolyzed to diacylglycerol, and the third acyl group is introduced by diacylglycerol acyltransferase (Dgat2) to form a TG. The reaction partner in all acylation reactions, fatty acyl-CoA, is synthesized by acyl-CoA synthase (Acsl1). In the current study, we found significant increases in expression of many genes implicated in TG synthesis (Gpd2, Agpat2, and Acsl1). A change in Dgat2 expression was also statistically significant but reflected an increase of only 27%. As a result of these changes in expression, serum TG concentrations in tetracycline-treated mice were increased significantly. Studies have demonstrated that a significant increase in acylglycerol-phosphate-acyltransferase activity by overexpression of the Agpat1 gene (Ruan and Pownall, 2001Go) increases TG synthesis, and a decrease in Agpat1 activity by mutation (Gangar et al., 2002Go) in adipocytes leads to decreased production of TGs. Thus, upregulation of genes involved in TG synthesis is a response that can lead to tetracycline-induced fat accumulation in the liver (Fig. 4A).

Cholesterol synthesis occurs in four phases: the synthesis of mevalonic acid from acetate, the conversion of mevalonic acid into isoprene units, the polymerization of isoprene units to squalene, and further processing to cholesterol. We found that tetracycline affects several genes in each steps of the cholesterol biosynthetic pathway. They include Hmgcs1 for the synthesis of HMG-CoA, Mvk for the phosphorylation of mevalonic acid, and Fdft1 for the formation of farnesyl diphosphate.

Sqle, encoding an enzyme adding one oxygen atom to the end of the squalene chain, is regulated by the cholesterol pool in the body. Cholesterol feeding reduced hepatic Sqle activity to inhibit synthesis of redundant cholesterol (Satoh et al., 1990Go). Similarly, a specific inhibitor of Sqle suppressed the secretion of cholesterol (Horie et al., 1993Go). Cell and animal models with decreased expression of 24-dehydrocholesterol reductase (Dhcr24) exhibit decreased cholesterol synthesis, and mutations in the Dhcr24 gene cause a cholesterol biosynthesis disorder (Waterham et al., 2001Go). In mice with a targeted disruption of Dhcr24, plasma and tissues contain almost no cholesterol, and desmosterol was shown to account for 99% of total sterols (Wechsler et al., 2003Go). Thus, tetracycline induction of Hmgcs1, Mvk, Idi1, Fdft1, Sqle, and Dhcr24 would facilitate increased hepatic cholesterol synthesis and a resultant fatty liver (Fig. 4B).

Mitochondrial ß-oxidation is a major process in fatty acid degradation. Although the first step is the conversion of fatty acid into fatty acyl-CoA by acyl-CoA synthase (Acsl1), the rate-limiting step is the transport of fatty acyl-CoA through the mitochondrial membrane via fatty acylcarnitine, catalyzed by carnitine acyltransferase (Crat). Tetracycline downregulated the expression of Crat1 and of 2,4-dienoyl-CoA reductase 1 (Decr1) and dodecenoyl-CoA delta isomerase (Dci), which play pivotal roles in mitochondrial ß-oxidation (Fig. 4C). Similarly, many other chemicals that induce fatty liver have also been reported to influence fatty acid metabolism (Guruge et al., 2006Go; Richards et al., 2004Go). Lipid catabolism is also mediated via microsomal {omega}-oxidation catalyzed by the Cyp4f family of enzymes (Le Quere et al., 2004Go; Sanders et al., 2005Go). Upregulation of Cyp4f14 and Cyp4f18 by tetracycline in our study suggests that the ß-oxidation pathway is inadequate to deal with the lipid accumulation induced by tetracycline, and a greater capacity for {omega}-oxidation is needed to increase the breakdown of hepatic lipids.

In summary, we documented a global gene expression profile for mouse liver following acute tetracycline treatment, as determined on oligonucleotide microarrays. We found that tetracycline affects gene expression associated with lipid metabolism in a direction to increase TG and cholesterol biosynthesis and decrease ß-oxidation of fatty acids. These results contribute not only to produce a database that can be used for prediction of toxicity through pattern recognition but also to our understanding of the mechanism of tetracycline-induced steatogenic hepatotoxicity.


    ACKNOWLEDGMENTS
 
This work was supported by a grant from Korea Food and Drug Administration (KFDA-05122-TGP-584).


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Amacher DE and Martin BA. (1997) Tetracycline-induced steatosis in primary canine hepatocyte cultures. Fundam. Appl. Toxicol. 40:256–263.[CrossRef][ISI][Medline]

Ballet F. (1997) Hepatotoxicity in drug development: Detection, significance and solutions. J. Hepatol 26:Suppl. 2, 26–36.[Medline]

Benjamini Y and Hochberg Y. (1995) Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. B 57:289–300.

Bhagavan BS, Wenk RE, McCarthy EF, Gebhardt FC, Lustgarten JA. (1982) Long-term use of tetracycline. JAMA 247:2780.[CrossRef][ISI][Medline]

Bolstad BM, Irizarry RA, Astrand M, Speed TP. (2003) A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 19:185–193.[Abstract/Free Full Text]

Freneaux E, Labbe G, Letteron P, The LD, Degott C, Geneve J, Larrey D, Pessayre D. (1988) Inhibition of the mitochondrial oxidation of fatty acids by tetracycline in mice and in man: Possible role in microvesicular steatosis induced by this antibiotic. Hepatology 8:1056–1062.[ISI][Medline]

Gabler WL, Smith J, Tsukuda N. (1992) Comparison of doxycycline and a chemically modified tetracycline inhibition of leukocyte functions. Res. Commun. Chem. Pathol. Pharmacol. 78:151–160.[ISI][Medline]

Gangar A, Raychaudhuri S, Rajasekharan R. (2002) Alteration in the cytosolic triacylglycerol biosynthetic machinery leads to decreased cell growth and triacylglycerol synthesis in oleaginous yeast. Biochem. J 365:577–589.[ISI][Medline]

Guruge KS, Yeung LWY, Yamanaka N, Miyazaki S, Lam PKS, Giesy JP, Jones PD, Yamashita N. (2006) Gene expression profiles in rat liver treated with perfluorooctanoic acid (PFOA). Toxicol. Sci. 89:93–107.[Abstract/Free Full Text]

Hamadeh HK, Bushel PR, Jayadev S, Martin K, DiSorbo O, Sieber S, Bennett L, Tennant R, Stoll R, Barrett JC, et al. (2002) Gene expression analysis reveals chemical-specific profiles. Toxicol. Sci. 67:219–231.[Abstract/Free Full Text]

Helma C. (2005) In silico predictive toxicology: The state-of-the-art and strategies to predict human health effects. Curr. Opin. Drug Discov. Dev. 8:27–31.[ISI][Medline]

Horie M, Hayashi M, Satoh T, Hotta H, Nagata Y, Ishida F, Kamei T. (1993) An inhibitor of squalene epoxidase, NB-598, suppresses the secretion of cholesterol and triacylglycerol and simultaneously reduces apolipoprotein B in HepG2 cells. Biochim. Biophys. Acta 1168:45–51.[Medline]

Huber W, von Heydebreck A, Sultmann H, Poustka A, Vingron M. (2002) Variance stabilization applied to microarray data calibration and to the quantification of differential expression. Bioinformatics 18:(Suppl. 1), S96–S104.[Abstract]

Le Quere V, Plee-Gautier E, Potin P, Madec S, Salaun JP. (2004) Human CYP4F3s are the main catalysts in the oxidation of fatty acid epoxides. J. Lipid Res. 45:1446–1458.[Abstract/Free Full Text]

Letteron P, Sutton A, Mansouri A, Fromenty B, Pessayre D. (2003) Inhibition of microsomal triglyceride transfer protein: Another mechanism for drug-induced steatosis in mice. Hepatology 38:133–140.[CrossRef][ISI][Medline]

Lumley CE and Walker SR. (1987a) Predicting the safety of medicines from animal toxicity tests. I. Rodents alone. Arch. Toxicol. Suppl. 11:295–299.[Medline]

Lumley CE and Walker SR. (1987b) Predicting the safety of medicines from animal toxicity tests. II. Rodents and non-rodents. Arch. Toxicol. Suppl. 11:300–304.[Medline]

Richards VE, Chau B, White MR, McQueen CA. (2004) Hepatic gene expression and lipid homeostasis in C57BL/6 mice exposed to hydrazine or acetylhydrazine. Toxicol. Sci. 82:318–332.[Abstract/Free Full Text]

Ruan H and Pownall HJ. (2001) Overexpression of 1-acyl-glycerol-3-phosphate acyltransferase-alpha enhances lipid storage in cellular models of adipose tissue and skeletal muscle. Diabetes 50:233–240.[Abstract/Free Full Text]

Sanders RJ, Ofman R, Valianpour F, Kemp S, Wanders RJ. (2005) Evidence for two enzymatic pathways for omega-oxidation of docosanoic acid in rat liver microsomes. J. Lipid Res. 46:1001–1008.[Abstract/Free Full Text]

Satoh T, Hidaka Y, Kamei T. (1990) Regulation of squalene epoxidase activity in rat liver. J. Lipid Res. 31:2095–2101.[Abstract]

Storck T, von Brevern MC, Behrens CK, Scheel J, Bach A. (2002) Transcriptomics in predictive toxicology. Curr. Opin. Drug Discov. Dev. 5:90–97.[ISI][Medline]

Suter L, Babiss LE, Wheeldon EB. (2004) Toxicogenomics in predictive toxicology in drug development. Chem. Biol. 11:161–171.[CrossRef][ISI][Medline]

Tachibana K, Kobayashi Y, Tanaka T, Tagami M, Sugiyama A, Katayama T, Ueda C, Yamasaki D, Ishimoto K, Sumitomo M, et al. (2005) Gene expression profiling of potential peroxisome proliferator-activated receptor (PPAR) target genes in human hepatoblastoma cell lines inducibly expressing different PPAR isoforms. Nucl. Recept. 3:3.[CrossRef][Medline]

Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, Botstein D, Altman RB. (2001) Missing value estimation methods for DNA microarrays. Bioinformatics 17:520–525.[Abstract/Free Full Text]

Waterham HR, Koster J, Romeijn GJ, Hennekam RC, Vreken P, Andersson HC, FitzPatrick DR, Kelley RI, Wanders RJ. (2001) Mutations in the 3beta-hydroxysterol Delta24-reductase gene cause desmosterolosis, an autosomal recessive disorder of cholesterol biosynthesis. Am. J. Hum. Genet. 69:685–694.[CrossRef][ISI][Medline]

Wechsler A, Brafman A, Shafir M, Heverin M, Gottlieb H, Damari G, Gozlan-Kelner S, Spivak I, Moshkin O, Fridman E, et al. (2003) Generation of viable cholesterol-free mice. Science 302:2087.[Free Full Text]

Westphal JF, Vetter D, Brogard JM. (1994) Hepatic side-effects of antibiotics. J. Antimicrob. Chemother 33:387–401.[Abstract/Free Full Text]

Wu Z, Bucher NL, Farmer SR. (1996) Induction of peroxisome proliferator-activated receptor gamma during the conversion of 3T3 fibroblasts into adipocytes is mediated by C/EBPbeta, C/EBPdelta, and glucocorticoids. Mol. Cell Biol. 16:4128–4136.[Abstract]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
94/1/206    most recent
kfl078v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (5)
Right arrowRequest Permissions
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Yin, H.-Q.
Right arrow Articles by Lee, B.-H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Yin, H.-Q.
Right arrow Articles by Lee, B.-H.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?