- Research article
- Open Access
In vitro gene regulatory networks predict in vivo function of liver
© Deng et al; licensee BioMed Central Ltd. 2010
Received: 29 June 2010
Accepted: 12 November 2010
Published: 12 November 2010
Evolution of toxicity testing is predicated upon using in vitro cell based systems to rapidly screen and predict how a chemical might cause toxicity to an organ in vivo. However, the degree to which we can extend in vitro results to in vivo activity and possible mechanisms of action remains to be fully addressed.
Here we use the nitroaromatic 2,4,6-trinitrotoluene (TNT) as a model chemical to compare and determine how we might extrapolate from in vitro data to in vivo effects. We found 341 transcripts differentially expressed in common among in vitro and in vivo assays in response to TNT. The major functional term corresponding to these transcripts was cell cycle. Similarly modulated common pathways were identified between in vitro and in vivo. Furthermore, we uncovered the conserved common transcriptional gene regulatory networks between in vitro and in vivo cellular liver systems that responded to TNT exposure, which mainly contain 2 subnetwork modules: PTTG1 and PIR centered networks. Interestingly, all 7 genes in the PTTG1 module were involved in cell cycle and downregulated by TNT both in vitro and in vivo.
The results of our investigation of TNT effects on gene expression in liver suggest that gene regulatory networks obtained from an in vitro system can predict in vivo function and mechanisms. Inhibiting PTTG1 and its targeted cell cyle related genes could be key machanism for TNT induced liver toxicity.
High-throughput toxicity testing is predicated upon using in vitro cell based systems to rapidly screen and predict how a chemical might cause toxicity to an organ in vivo. Recent microarray studies have shown that the in vitro gene expression profiles in liver slices treated with various compounds could predict the toxicity and pathology observed in vivo. However, the degree to which in vitro results can be extended to in vivo activity and possible mechanisms of action remains to be fully addressed.
Compared with animal models, primary cell cultures have advantages for investigating mechanisms of chemical toxicity. Primary cultured cells, such as hepatocytes, can offer a convenient system that is easily genetically manipulated and can be used to test various throughput screens using different molecular and biochemical methods. Use of primary cell cultures can also reduce cost and mitigate animal welfare concerns inherent in in vivo studies . In vitro systems have a long history of use in screening new drugs for human diseases such as cancer and in studying cellular and molecular events of different molecules (e.g., pharmaceuticals and xenobiotics) [4, 5].
In this study, we used the nitroaromatic 2,4,6 -rinitrotoluene (TNT) as a model chemical to compare and determine how we might extrapolate in vitro data to in vivo systems. We compared the gene expression profiles of in vitro primary liver cells with the gene expression profiles of in vivo liver tissue of rats exposed to TNT.
TNT is a munitions compound which is released to the environment as a result of military training activities and manufacturing processes [6, 7]. The major toxic effects to rats in response to TNT exposure are methemoglobinemia, anemia, hypercholesterolemia, and hepatomegaly, splenomegaly, and testicular atrophy with their accompanying histologic lesions . Cytotoxic and genotoxic effects are also caused by TNT [9–11]. Our recent toxicity studies in rats found that TNT decreased body weight, increased liver weight, and induced erythrocytosis (Deng et al., unpublished data).
Although a variety of studies have focused on toxicity aspects of TNT, the underlying mechanisms of toxicity induced by this compound are largely unknown. Several laboratories have started to use microarray analysis to understand the molecular mechanisms of invertebrate ecotoxicity. The exposure of the earthworm Eisenia fetida with TNT regulates the expression of genes involved in multiple biological processes including muscle contraction, neuronal signaling and growth, ubiquitinylation, fibrinolysis and coagulation, iron and calcium homeostasis, oxygen transport, and immunity . Microarray analysis of TNT exposed Arabidopsis thaliana (Arabidopsis) reveals the induced expression of oxophytodienoate reductases (OPRs), the protein responsible for TNT detoxification in Arabidopsis. The same team also found that bifunctional O- and C-glucosyltransferases is involved in combating TNT induced toxicity in Arabidopsis [2, 13]. Eisentraeger et al.  investigated the mechanism of human liver cell line HepG2 treated with TNT using the toxicology cDNA microarray. They found that the detoxification metabolism in the HepG2 cells by TNT induced genes encoded phase I and phase II enzymes.
In the current study, we examined gene transcriptional responses to short term (24 h and 48 h) TNT exposures in rat in both in vitro hepatocytes and in vivo liver. We observed that a number of genes were commonly regulated by both in vitro and in vivo TNT treatments. Moreover, we found commonly affected pathways in both in vivo and in vitro systems exposed with TNT. Functional analysis indicates that both TNT treatments in vivo and in vitro impact genes involved in cell cycle, cell growth and cell death signaling, detoxification response, lipid metabolism and immune response, which can reasonably account for the physiological dysfunctions induced by TNT. In addition, we identified conserved gene networks inferred from in vivo and in vitro networks. Our results suggest that in vitro system can predict in vivo functions and mechanisms based on gene expression profiles, and gene network is a valuable approach for predicting in vivo function using in vitro data.
Commonly regulated genes by in vivo and in vitro TNT treatments
To determine if in vitro TNT exposure could be used to predict in vivo TNT exposure, the first aim of this experiment was to determine if there were common genes differentially regulated by both in vivo and in vitro TNT treatments. For the in vivo experiment, rats were gavaged with TNT at: 0 (control), 4.8, 48, 96 or 192 mg/kg and sacrificed after 24 h or 48 h. RNA was isolated for microarray hybridization. Agilent rat whole genome array (Agilent Technologies Inc., Palo Alto, CA) was employed in the study. Four biological replicates were used for each unique condition as described in the Material and Methods. For the in vitro experiments, primary rat hepatocytes (Lonza, Walkersville, MD) were cultured and treated with TNT at 10 ppm or vehicle control for 24 h, after which RNA was isolated for microarray experiments. The in vivo and in vitro micorarray data have been deposited in the GEO databases with assigned numbers GSE19628 and GSE19662, respectively.
To identify genes differentially expressed in vivo, a one-way ANOVA was applied to identify differentiated transcripts across different doses of TNT treatments at each of the 24 h or 48 h time points. A p value cut-off of 0.05 and a 1.5 fold gene transcript change by comparing at least one pair of conditions between different doses including control were used as filters to identify differentially expressed genes. A total of 4875 gene transcripts were significantly and differentially expressed (i.e., upregulated and downregulated) in in vivo TNT treatment for all doses and time points.
Similar gene expression pattern shared by commonly regulated in vivo and in vitro genes
Functional analysis of commonly regulated genes
Significant functional categories based on genes commonly regulated by TNT in vivo and in vitro
ANGPTL2, LGALS1, CYP26B1, COL1A1, CCND2, CKAP2, CDKN3, KIF20A, CXCL12, ECT2, CCNB1, RARB, PTTG1, LZTS1, TFDP2, BCL2A1, CDC2, INPP5D
MYC, CCNC, DNM1L, NRG2(includesEG:9542), MAPK1, NRG1, TFRC, PPARG, CREG1, RIOK3, NTRK1, DUSP13, ADM
GCK, MMP2, INPP5 D, PTTG1, CXCL12, LGALS1
NQO1, ADM, PLA2G7, UGDH, GCLC, PPP1R3C, NTRK1, PPARG, JMJD7, PLA2G4B, ABCC3, UGT1A6, CPT1A, PARD3, NRG1, H6PD, SLC5A3, MYC
GCK, PTTG1, CXCL12, LGALS1
MYC, SLC5A3, ABCG5, EIF2S1, GMFB, H6PD, MAPK1, NRG1, PARD3, CPT1A, TFRC, ABCC3, PPARG, AQP8, ABCC4, PPP1R3C, GCLC, ADM, NQO1
Cellular Growth and Proliferation
LGALS1, COL1A1, CCND2, KIF20A, CXCL12, CCNB1, RARB, PTTG1, LZTS1, HSD11B2, DLC1, COL1A2, BCL2A1, CDC2, INPP5 D, MMP2
MYC, DNAJB6, NRG2 (includes EG:9542), CES2 (includes EG:8824), MAFF, CDH4 (includesEG:1002), UGT2B17, MAPK1, NRG1, TRIM35, TFRC, PPARG, CXCL2, ALDH1L1, CREG1, GSTP1, HMGCR, NTRK1, CDA, PHLDA1, PFN2, GRIN2C, ADM, CYP1A1
LGALS1, CYP26B1, CCND2, CKAP2, CXCL12, CCNB1, RARB, PTTG1, HSD11B2, DLC1, BCL2A1, CDC2, INPP5 D, MMP2
MYC, DNAJB6, EIF2S1, DNM1L, NRG2(includes EG:9542), GMFB, HTATIP2, MAPK1, NRG1, TRIM35, TFRC, EPHX1, PPARG, CXCL2, SQSTM1, CYP2F1, GSTP1, GSR, NTRK1, ABCC4, PHLDA1, GCLC, PLA2G7, TXNRD1, NR1I3, GRIN2C, ADM, NQO1, NCF2
DNA Replication, Recombination, and Repair
LGALS1, CCND2, CXCL12, ECT2, CCNB1, PTTG1, MMP2
MYC, NRG2 (includes EG:9542), PDE5A, NRG1, PPARG, RIOK3, GSTP1, GCLC, NR1I3, ADM, AMPD3, NQO1
GCK, INPP5 D, HSD11B2, CXCL12, PNPLA3, LGALS1
CYP1A1, ADM, PLA2G7, PPP1R3C, UGT2B7, NTRK1, ACOT4, GSTP1, AQP8, PPARG, JMJD7-PLA2G4B, ABCC3, RDH16, CPT1A, MAPK1, H6PD, UGT2B17, CYP3A43, CYP2C18, ABCG5, MYC
Cellular Assembly and Organization
LGALS1, COL1A1, ECT2, CCNB1, PTTG1, COL1A2, CDC2, KRT20, INPP5D
MYC, DNAJB6, DNM1L, EPB41, NRG1, RIOK3, PFN2, ADM
Immune Cell Trafficking
COL1A1, CYTIP, CXCL12
Humoral Immune Response
CXCL12, BCL2A1, NPP5D
Other commonly significantly affected functional terms were cell growth and proliferation (40 genes), cell death (19 genes), DNA replication (19 genes), and recombination and repair (19 genes). Interestingly, there were 10 genes that were repeatedly present in all four functional terms: ADM, CCNB1, CCND2, CXCL12, LGALS, MYC, NRG1, NRG2, PPARG, and PTTG1. In addition, these genes were consistently regulated by both in vivo TNT and in vitro TNT in the same directions. The expression of CCNB1, CCND2, PTTG1, CXC12, and LGALS1 was decreased by TNT, and the expression of ADM, MYC, NGR1, NGR2, and PPARG was elevated by TNT.
There were 23 commonly regulated genes that were involved in molecular transport (Table 1). Nineteen of the 23 genes (83%) were upregulated, and only 4 of them were downregulated by both in vitro and in vivo TNT treatments. Some significant upregulated genes included NQO1, ADM, ABCC3 and ABCC4. Our results suggest that TNT can enhance molecular transport both in vivo and in vitro liver systems.
Twenty-three genes contributing to carbohydrate metabolism were commonly regulated by both in vivo and in vitro TNT treatments. Among them, only 6 genes were repressed and 17 genes were induced. The expression of 27 genes participating in lipid metabolism was significantly changed in response to both in vivo and in vitro TNT exposures. Similar to carbohydrate metabolism, we saw more upregulated genes (21) than downregulated genes (6) in this category. Several cytochrome P450 family members such as CYP1A1, CYP3A43 and CYP2C18 were in the upregulated gene list. Interestingly, 10 genes (MYC, PPARG, PPP1R3C, CXCL12, GCK, ABCC3, ADM, CPT1A, H6PD, and LGALS1) belong to the three functional categories: molecular transport, carbohydrate metabolism and lipid metabolism. Three of these 10 genes (LGALS1, GCK, CXCL12) were downregulated and the other 7 genes were upregulated by both in vivo and in vitro TNT treatments.
There were 17 commonly regulated genes that were involved in cellular assembly and organization. The expression of 10 genes was decreased, while expression of 7 genes were increased. Eight genes played a role in immune response, 5 of which were downregulated and 3 genes upregulated. More immune response genes were downregulated, indicating that TNT may interfere with normal immune function to induce its toxicity.
Commonly regulated canonical pathways based on in vivo, in vitro and common gene lists regulated by TNT
-Log(P-value) in vivo
-Log(P-value) in vitro
-Log(P-value) in vivo and in vitro
Metabolism of Xenobiotics by Cytochrome P450
NRF2-mediated Oxidative Stress Response
LPS/IL-1 Mediated Inhibition of RXR Function
Xenobiotic Metabolism Signaling
Pentose and Glucuronate Interconversions
Aryl Hydrocarbon Receptor Signaling
Fructose and Mannose Metabolism
Biosynthesis of Steroids
Fatty Acid Metabolism
Starch and Sucrose Metabolism
Androgen and Estrogen Metabolism
Top common pathways regulated by TNT in vitro and in vivo
Common regulated genes*
Metabolism of Xenobiotics by Cytochrome P450
CYP3A43, CYP2F1, GSTP1, CYP1A1, CYP2C18, UGT1A6, UGT2B7, GSTA5, ALDH1L1, UGT2B17, EPHX1
NRF2-mediated Oxidative Stress Response
GSR, AKR7A3, GSTP1, MAPK1, GSTA5, NQO1, GCLC, SQSTM1, DNAJB6, MAFF, TXNRD1, EPHX1
Xenobiotic Metabolism Signaling
GSTP1, CYP1A1, UGT1A6, MAPK1, GSTA5, NQO1, ALDH1L1, GCLC, CES2 (includes EG:8824), UGT2B7, NR1I3, UGT2B17, ABCC3
Aryl Hydrocarbon Receptor Signaling
MYC, GSTP1, CYP1A1, CCND2, MAPK1, GSTA5, NQO1, RARB, ALDH1L1
Pentose and Glucuronate Interconversions
TCAG7.1260, AKR7A3, UGDH, UGT1A6, UGT2B7, UGT2B17
LPS/IL-1 Mediated Inhibition of RXR Function
GSTP1, ABCG5, CPT1A, GSTA5, NR1I3, ALDH1L1, CES2 (includes EG:8824), ABCC3, ABCC4
GSR, GSTP1, GSTA5, H6PD, GCLC
CPT1A, NR1I3, CES2 (includes EG:8824), ABCC3
Interestingly, all the genes in Table 3, except for CCND2 and retinoid acid receptor beta (RARB), were commonly upregulated by both in vivo and in vitro TNT treatments. MAPK1, a critical intracellular signaling protein involved in multiple cellular functions (Table 1) was also presented in multiple pathways such as NRF-2 mediated oxidative stress response, xenobiotic metabolism signaling, and aryl hydrocarbon signaling pathways.
Inferring in vitro and in vivo gene networks
Since transcriptional regulated gene networks are trigged by transcription factors, we then counted the connections of the transcription factors (TFs) in the 341 transcripts for the both in vivo and in vitro gene networks. Top TFs with most connections in the in vitro gene network were cellular repressor of E1A-stimulated genes 1 (CREG1), pirin (iron-binding nuclear protein) (PIR), BTB (POZ) domain containing 11(BTBD11), SAM pointed domain containing ets transcription factor (SPDEF) and basic helix-loop-helix family, member e41 (BHLHE41). TFs with most connections in the in vivo gene network included HIV-1 Tat interactive protein 2, 30kDa (HTATIP2), transcription factor Dp-2 (E2F dimerization partner 2) (TFDP2), leucine zipper, putative tumor suppressor 1 (LZTS1), zinc finger protein 217 (ZNF217), PIR, and MYC. The transcription factors that had over 15 connections in the both gene networks included CREG1, PIR, BTBD11, LZTS1, ASB2, HES6, and SQSTM1(Figure 6C). Two TFs, SPDEF and BHLHE41, had at least twice more connections in the in vitro gene network than that in the in vivo gene network. Seven TFs including HTATIP2, MYC, NEQ1, ZNF394, E2F3, DNAJB6 and TFDP2 possessed over twice as many connections in the in vivo gene network than that in the in vitro gene network.
Conserved in vitro and in vivo gene networks
The TFs PIR, BTBD11, HTATIP2, and CREG1 were colocalized in the biggest conserved sub-network (herein called the PBHC network). Eight genes were connected with the TF PIR centered sub-gene network module (Figure 7), and contained highly enriched genes associated with oxidoreductase activity (ACOT4, AKR7A3, RDH16, TXNRD1, and UGDH). AKR7A3 and TXNRD1 also participate in NRF2-mediated oxidative stress pathway. All the genes in this network module were upregulated by both in vivo and in vitro TNT exposures. This result suggests that oxidative stress is activated by TNT .
The TF CREG1 centered sub-network module included 5 genes (Figure 7), 3 of which were also cell cycle related genes (CREG1, NTRK1 and DUSP13). In contrast to PTTG1 centered sub-network, all the genes in the CREG sub-network module were upregulated by both in vivo and in vitro TNT treatments. Glutathione reductase (GSR), a gene related to oxidoreductase activity, and TXNRD1 provided the connections between the CREG1, HTATIP2, and PIR sub-networks (Figure 7), creating a larger 13 gene sub-network. In addition, GSR also participated in the NRF2-mediated oxidative stress pathway. Indeed, we found that all the genes in the biggest conserved sub-network formed by CREG1, HTATIP2, BTBD11, and PIR mediated modules were consistently upregulated by both in vivo and in vitro TNT additions.
Verification of microarray responses using real time QRT-PCR
For instance, the gene ABCC3 was shown in microarray (Additional file 2, Figure S 1A) to be upregulated at both 24 h and 48 h with a clear dose response to in vivo TNT treatment, and was also induced by TNT in vitro, and the same pattern was exhibited by QRT-PCR (Additional file 2, Supplementary Figure 1A). In both microarray and QRT-PCR results, the gene AKR7A3 was upregulated at both time points in vivo as well as upregulated in vitro. Interestingly, AKR7A3 was significantly upregulated at the highest dose (192 mg/kg) in vivo at 48 h but not at other lower doses, which was shown consistently by both microarray and QRT-PCR (Additional file 2: Figure S1 A). For the gene CYP1A1, we saw it was upregulated more in response to in vitro TNT than in vivo TNT treatment, and was induced much more at 24 h than 48 h in vivo, which was consistently revealed by both microarray and QRT-PCR results (Additional file 2, Figure S 1C). Three genes PTTG1 (Figure 8), CCNB1 and CCND2 that were in PTTG1 centered subnetwork were consistently shown to be downregulated by TNT in vivo and in vitro, with both approaches (Additional file 2, Figures S 1B and F). The QRT-PCR results could exactly confirm our two observations obtained from the microarray data: the genes up or downregulated in the same directions in vitro and in vivo, and a clear dose response in vivo.
In this study, we compared gene expression profiles regulated by TNT exposure in vivo and in vitro from shared gene lists, functional terms, common pathways, and conserved networks. Overall, good consistency existed between in vitro and in vivo exposures. The Venn diagram (Figure 1) shows that a only small percentage of the in vivo transcripts were also regulated in the in vitro experiments. Because the regulated gene number in vivo is the total number from all TNT treatments including multiple doses and two time points (24 h and 48h), the gene number is much larger than that of in vitro exposures. However, the in vitro gene number comes from only one dose. The reason we used one dose is that the in vitro experiment is a large one with a total of over 105 chemicals and 531 arrays. But if we compare the overlapped transcripts with the in vitro transcripts (940), the number of overlapped transcripts is not small, with over one third of the total in vitro transcripts. This observation is found more than in overlaps examined in other publications (. Also we would like to emphasize that the gene expression patterns for the 192 mg/kg, 24 h in vivo TNT treatment and the 24 h in vitro TNT treatment are similar. If we only count the regulated transcript number for 192 mg/kg, 24 h in vivo TNT treatment (Figure 2), the number is comparable to that of the 24 h in vitro TNT treatment. In addition, in vivo and in vitro are two different systems, we cannot expect that they are exactly the same. Our purpose is to identify common transcripts between these two systems, thereby to understand the in vivo function through the in vitro system. From the perspective of functional analysis, both TNT treatments in vivo and in vitro influence genes involved in cell cycle, cell growth and cell death signaling, detoxification response, lipid metabolism and immune reponse, which can reasonably explain the physiological dysfunctions induced by TNT. We also found conserved sub-networks between inferred networks from in vivo and in vitro TNT regulated gene expression profiles.
Cell cycle, cell growth, and cell death signaling
A large number of genes involved in cell cycle, cell growth and cell death were commonly regulated in both in vitro and in vivo systems perturbed by TNT. The highly represented downregulated genes were CCNB1, CCND2, PTTG1, CXC12, and LGALS1, and the highly represented upregulated genes were ADM, MYC, NRG1, NRG2, and PPARG. The cyclin B1 (CCNB1) gene product complexes with p34(cdc2) to form the maturation-promoting factor (MPF). CCNB1 is expressed predominantly during G2/M phase because of its critical role in cell mitosis [17, 18]. The protein encoded by CCND2 belongs to the highly conserved cyclin family, whose members are characterized by a dramatic periodicity in protein abundance through the cell cycle. Cyclins function as regulators of CDK kinases. Different cyclins exhibit distinct expression and degradation patterns which contribute to the temporal coordination of each mitotic event. The CCND2-translated cyclin forms a complex with and functions as a regulatory subunit of CDK4 or CDK6, whose activity is essential for cell cycle G1/S transition. PTTG1 encoded protein is a homolog of yeast securin proteins, which prevent separins from promoting sister chromatid separation. The gene product contains 2 PXXP motifs, which are essential for its transforming and tumorigenic activities, as well as for its stimulation of basic fibroblast growth factor expression [19–21]. Besides the functions of cell cycle and growth, the major role of CXC12 is involved in immune response. The protein encoded by LGALS1 is a member of beta-galactoside-binding protein family implicated in modulating cell-cell and cell-matrix interactions, which may act as an autocrine negative growth factor that regulates cell proliferation .
The protein product of the gene Adrenomedullin (ADM) is a multifunctional peptide vasodilator that carries out its functions through calcitonin receptor-like receptor/receptor activity modifying protein-2 and -3 (CLR/RAMP2 and CLR/RAMP3). It can positively enhance cell proliferation . MYC is a well known multifunctional, nuclear phosphoprotein that contributes to cell cycle progression, apoptosis and cellular transformation . It acts as a transcription factor that regulates transcription of specific target genes. It has been implicated as an oncogene to facilitate cell growth and survival. The protein encoded by NRG1 was initially characterized as a 44-kD glycoprotein that interacts with the NEU/ERBB2 receptor tyrosine kinase to increase its phosphorylation on tyrosine residues. This protein is a signaling protein that regulates cell-cell interactions and is involved in the promotion of growth and development of multiple organ systems [24, 25]. Neuregulin 2 (NRG2) is a new member of the neuregulin family of growth and differentiation factors. By interacting with the Erbb family of receptors, NRG2 induces the growth and differentiation of epithelial, neuronal, glial, and other cell types . The gene PPARG encodes a member of the peroxisome proliferator-activated receptor (PPAR) subfamily of nuclear receptors. PPARs form heterodimers with retinoid X receptors (RXRs) and these heterodimers regulate transcription of various genes to execute various functions including lipid metabolism and cell growth [27, 28].
The gene BCL2A1, encodes a member of the BCL-2 protein family, a well known protein family functioning as anti- and pro-apoptotic regulators. The protein encoded by this BCL2A1 is able to reduce the release of pro-apoptotic cytochrome c from mitochondria and block caspase activation, thereby enhance cell survival. The reduction of the expression of this gene could partially account for the cell damage induced by TNT.
Interestingly, the downregulated genes involved in cell growth and cell death also play a role in the cell cycle process. Therefore, the cell growth and other possible cell death related phenotypes may occur primarily through interrupting cell cycle progression of cells. Recently we found that there was weight loss in rats treated with TNT (data not shown). The possible reason for the body weight loss is that cell growth is reduced by TNT. TNT could inhibit the growth of V79 and TK6 human lymphoblastic cells . It could explain cytotoxic effects caused by TNT [9–11].
Our results revealed a couple of pathways that were commonly regulated by both in vivo and in vitro TNT treatments. From the pathways listed in Table 3, we identified xenobiotic metabolism signaling as the key of all the pathways because all the other pathways are involved in xenobiotic metabolism signaling. Certainly, several phase I and II metabolizing enzymes with their family members that take part in multiple pathways play a pivotal role in the detoxification process. These enzymes mainly cover CYPs, GSTs, UGTs, AKR7A3, ALDHL1, NQO1, and EPHX1. As expected, all these enzymes were consistently elevated by both in vitro and in vivo TNT treatments.
CYP1A1 has been a hallmark for the treatments of many toxins and carcinogens [31–34]. CYP1A1 was strongly induced in both in vivo and in vitro liver cells by TNT. Due to its catalytic function, it could play dual roles in the cells in the presence of TNT. Its activation leads to the production of reactive oxygen species (ROS)  which induces oxidative stress, lipid peroxidation (LPO), protein modification and denaturation, and DNA damage, and thus induces toxicity. Meanwhile, it may also directly metabolize TNT or its metabolites, resulting in cell survival. Glutathione transferases (GSTs) catalyze the conjunction of reduced glutathione (GSH) to electrophiles and oxidative stress products, thus facilitating their removal .
UGTs catalyze the generation of glucuronide conjugates of dihydrodiols to execute their detoxification function . AKR7A3 catalyzes aflatoxin B1 (AFB 1) and suppresses AFB1 dialdehyde metabolite to its corresponding mono and dialcohols. Hence, the activation of AKR7A3 attenuates the toxicity of the AFB1-dialdehyde that reacts with proteins, and thus decreases AFB1-induced toxicity .
ALDH1L1 is well recognized as protectors against ROS caused oxidative damage. It metabolizes reactive products of toxic LPO, which include 4-hydroxy-2-nonenal (4-HNE) and malondialdehyde (MDA) . This protective mechanism could explain how ALDH responds to the cellular recovery to TNT exposure.
As a cytosolic flavoenzyme, NQO1 catalyzes the two-electron reduction of diversified substrates . NQO1 is characterized as a detoxification enzyme mainly because of its capability of degenerating quinone substrates to their less toxic hydroquinone metabolites, bypassing the redox-cycling semiquinone radical [41, 42]. NQO1 can also transform ubiquinone and vitamin E quinone to their antioxidant forms .
EPHX1 is one member of epoxide hydrolases (EPHs) that catalyze the hydrolysis of electrophiles such as epoxides to the less reactive vicinal diols, which can explain the mechanism of epoxide hydrolases as classical detoxifying enzymes. EPH is able to inactivate highly diversified reactive epoxides with different structures, and therefore plays a critical role in the enzymatic defense against adverse effects of xenobiotic compounds .
Many upregulated genes involved in xenobiotic metabolism signaling also participate in cell cycle, cell growth, and cell death. Our results suggest that, upon exposure to TNT, cells arrest the cell cycle and increase detoxification and oxidative stress enzymes to remove the xenobiotic and eliminate any cell damage caused by the xenobiotic..
Several genes participating in lipid metabolism were commonly regulated by both in vivo and in vitro TNT exposures. Most of these lipid metabolism genes showed elevated expression, while only a small number of genes were repressed. The downregulated genes included LGALS1, PNPLA3, CXCL12, HSD11B2 and INPP5 D. One well known gene, inositol polyphosphate-5-phosphatase, 145kDa (INPP5 D, or SHIP-1), encodes a protein that is largely confined to hematopoietic cells. It is a well characterized inhibitory molecule that is recruited by engagement of the inhibitory Fcγ type IIB receptor in B cells and mast cells or by engagement of Fcε type I or Fcγ type III, cytokine, and growth factor receptors in myeloid cells . Once recruited to the plasma membrane by signaling complexes, its enzymatic activity depletes PtdIns(3,4,5)P3 and prevents membrane localization of some PH domain-containing effectors, eventually leading to impaired PI3K-dependent signaling events.
The protein encoded by the gene PNPLA3 is a member of the adiponutrin family complement the hormone sensitive lipase (HSL) as responsible for adipocyte triacylglycerol lipase activity. Mice lacking HSL reveal a lean phenotype and accumulate diglycerides suggesting that HSL is the main enzyme for the second step of lipolysis . LGALS1 has been shown to be involved in many other functions in Table 1, also plays a role in lipid metabolism, by increasing the induction of levels of phosphatidylserine as well as the mobilization of phosphatidylserine . Hydroxysteroid (11-beta) dehydrogenase 2 (HSD11B2) has been reported to contribute to the metabolism of aldosterone, cortisone, glucorticoid and hydrocortisone [47, 48]. The downregulated genes involved in lipid metabolism may indicate that TNT may interrupt lipid metabolism to some degree. The lipid metabolism was also impaired in the liver of quail exposed to 2,6-DNT . Meanwhile, we saw many other genes involved in lipid metabolism upregulated by TNT, many of which overlap with xenobiotic signaling genes (Table 3), suggesting that they also function for detoxification by trying to recover normal lipid metabolism.
As described earlier in the results, 5 of 8 genes associated with immune response were downregulated by TNT both in vivo and in vitro. These 5 genes are CXCL12, COL1A1, BCL2A1, INPP5 D, and CYTIP.
Chemokine (C-X-C motif) ligand 12 (CXCL12 or SDF-1) was strongly downregulated by TNT. It is a ligand for the G-coupled receptor protein chemokine (C-X-C motif) receptor CXCR4. Activation of CXCR4 by CXCL12 is involved in many biological functions such as cell migration, growth, and survival .
Besides enhancing cell survival, BCL2A1 is a direct transcription target of NF-κB in response to inflammatory mediators, and is upregulated by different extracellular signals such as granulocyte-macrophage colony-stimulating factor (GM-CSF), CD40, phorbol ester and inflammatory cytokines TNF and IL-1, which suggests a cytoprotective function that is essential for lymphocyte activation, and plays a role in immune response .
COL1A1 encodes the pro-alpha 1 chains of type I collagen whose triple helix is composed of two alpha 1 chains and one alpha 2 chain. Type I collagen is a fibril-forming collagen found in most connective tissues and is abundant in bone, cornea, dermis and tendon. Thus this gene is critical for cell aggregation, migration, proliferation, binding, adhesion and growth . This gene carries out its immune response by primarily affecting transmigration of T lymphocytes. As a gene primarily functioning in lipid metabolism, INPP5 D participating in immune response by mainly maintaining the quantity of B lymphocytes and participating in phospholipid metabolism . Cytohesin 1 interacting protein (CYTIP) encodes a protein containing 2 leucine zipper domains and a putative C-terminal nuclear targeting signal. Its major role is involved in the quantity maintenance, migration of leukocytes and lymphocytes .
The decreased above gene expression could explain the interference of normal immune response by TNT both in vitro and in vivo, which could lead to inflammation and other immune related toxicities caused by TNT.
Network as a valuable approach for predicting in vivo function using in vitro data
Finally, not only could we find commonly regulated genes, functional terms, and pathways, but we were also able to identify conserved gene networks between in vitro and in vivo liver systems perturbed by TNT. There were two big conserved sub-networks, the TF PTTG1 centered gene network and the TFs PIR, BTBD11, HTATIP2 and CREG1-connected gene network (PBHC network). Interestingly, all the genes in the PTTG1 sub-network were associated with cell cycle function. For example, CDC2 and cyclin B (CCNB1) form a complex which is responsible for the onset of metaphase . PTTG1, CDC2, and CCNB1 all participated in the mitotic roles of Polo-Like kinase pathway [56, 57]. Furthermore, all the genes in the PTTG sub-network were downregulated by TNT in vitro and in vivo. Since cell cycle is the most significant functional term affected by TNT both in vitro and in vivo (Table 1), this conserved cell cycle relating network could well explain the common mechanism of TNT treated liver systems both in vivo and in vitro, which is that cell cycle progression is interfered by TNT, thereby leading to cell growth inhibition and cell death, ultimately cause toxicity [58–60]. We did observe liver cell death in rats exposed to high dose TNT (Deng et al., unpublished data).
In the PBHC sub-network, some genes such as CREG1, NTRK1, DUSP13, GSR, and TXNRD1 play a role in cell cycle and/or cell death. Many genes possess oxidoreductase activity , and are involved in oxidative stress process and play a crucial role in metabolizing toxic compounds and their metabolites. The TF PIR encodes a TF which is a member of the cupin superfamily. The encoded protein is a Fe(II)-containing nuclear protein expressed in all tissues of the body and concentrated within dot-like subnuclear structures . Interactions with nuclear factor I/CCAAT box transcription factor as well as B cell lymphoma 3-encoded oncoprotein suggest the encoded protein may function as a transcriptional cofactor and be involved in the regulation of DNA transcription and replication. Interestingly, the TF HTATIP2 also has oxidoreductase activity , and connects two genes GSR and TXNRD1 that both possess oxidoreductase activity as well as both take part in NRF2-mediated oxidative stress pathway. HTATIP was reported to be regulated by NRF2, too .
Overall, we can see that the genes in the PBHC sub-network are involved in cell cycle growth, and/or oxidoreductase activity. Surprisingly, 17 genes in the PBHC sub-network (Figure 7) were all consistently upregulated by TNT both in vivo and in vitro. The induction of these gene expression lead to cell survival and removal of toxic products, therefore, this conserved sub-network may well explain the common detoxification mechanism both in vivo and in vitro liver systems treated with TNT.
We did not use only TNT-exposed samples but included all samples in the pool to construct static gene networks, which aim to uncover any gene interactions. In so doing, the more samples and conditions included make the network analyses more robust. For instance, Faith et al. (2007) used 445 arrays with over 180 various conditions to build an E. Coli static gene network. Because the genes we used are TNT-regulated genes, the constructed gene networks are specifically TNT-regulated networks.
We built separate gene networks using in vitro and in vivo samples and found commonly conserved networks. These common networks are ideal because they validate both in vivo and in vitro conditions.
Our results indicate that TNT perturbed similar gene regulated networks when in vitro effects on hepatocyte cells are compared with in vitro effects on liver. Thus, gene regulatory networks obtained from an in vitro system can be predictive of in vivo function in the liver when perturbed by a chemical stressor such as TNT. Furthermore, PTTG1 regulated cell cycle may be a key targeted process indicator for TNT induced toxicity.
This study is the first report to use an in vitro transcriptional regulatory gene network to predict in vivo toxicity and mechanism induced by a chemical. Knowledge achieved from this innovative study can provide an efficient way to assess whether a soldier or civilian has been exposed to TNT and to find possible ways to prevent, treat, and reduce TNT induced adverse effects.
Microarray experimental design
Changes in gene expression were tested using Agilent commercial whole rat genome microarrays (4 X 44K). For in vivo experiment, one of five compounds (2,4-DNT, 2,6-DNT, TNT, 2-ADNT and 4-ADNT) were used to treat rats. Female Sprague-Dawley rats (175-225 grams) were from the in-house breeding colony (College of Pharmacy, University of Louisiana at Monroe [ULM] and treated in accordance with the Guide for Use and Care of Animals. Breeders were from Harlan-Sprague Dawley in Madison, WI. Housing consisted of a 12 h light/dark cycle with ad libitum access to tap water and rodent chow (Harlan/Teklad 7012, Madison, WI). Rats were housed individually in polycarbonate cages on hardwood bedding (Sani-chips, Harlan/Tekland, Madison, WI) one week prior to treatment. Food was withdrawn the night before treatments, which were administered by gavage between 8 and 10 AM. Study protocols were preapproved by the ULM Animal Care and Use Committee.
Groups of rats were weighed and randomly assigned to treatment. Treatments were vehicle (5% v/v DMSO in corn oil), 2,4-DNT (4.8, 48, 96, and 192 mg/kg), 2A-DNT (4.4, 44, 87 and 174 mg/kg) and 4A-DNT (4.7, 47, 94 and 187 mg/kg), TNT (4.8, 48, 96 and 192 mg/kg), and 2,6-DNT (5.0, 25, 50 and 99 mg/kg). Rats were observed continuously for the first hour after dosing, hourly for 8 h and daily thereafter. Moribund rats were euthanized with CO2. Livers were excised and weighed. A portion of the liver was removed and placed in RNA Later (Ambion) following manufacturer's instruction and later used for genomic analyses. Remaining liver was flash frozen in liquid N2 and stored at -70°C for further analyses.
For in vitro experiment, primary cultured rat hepatocytes (Lonza, Walkersville, MD) were kept in an incubator at 37°C and 5% CO2. After seeding in flasks (3x106 cells/flask), the cells were treated with 10 mg/l TNT and vehicle control (DMSO). Three experimental replicates were used. After 24 h, the cells were isolated and lysed for RNA analysis. Cells were also treated with one of 106 compounds, including TNT (10 mg/l), and relative controls. At least three biological replicates were used for each unique condition (Additional file 3, Table S 3).
Total RNA extraction
Total RNA was extracted from about 30 mg of liver tissue or cell pellet. Tissues or cells were homogenized in the lysis buffer with FAST Prep-24 from MP at speed 6.0/s twice, each last 30 s before using RNeasy kits (Qiagen). Total RNA concentrations were measured using NanoDrop® ND-1000 Spectrophotometer (NanoDrop technologies, Wilmington, DE, USA). The integrity and quality of total RNA was checked on an Agilent 2100 Bioanalyzer (Palo Alto, CA). The gel-like images generated by the Bioanalyzer show that total RNAs have two bands, represent 18 S and 26 S RNA of mammalian RNA. Nuclease-free water (Ambion) was used to elute total RNA.
Rat whole genome oligo arrays in the format of 4X44K were purchased from Agilent. Sample cRNA synthesis, labeling, hybridization and microarray processing were performed according to manufacturer's protocol "One-Color Microarray-Based Gene Expression Analysis" (version 1.0). The Agilent One-Color Spike-Mix (part number 5188-5282) was diluted 5000-fold and 5 μl of the diluted spike-in mix was added to 1000 ng of each of the total RNA samples prior to labeling reactions. The labeling reactions were performed using the Agilent Low RNA Input Linear Amplification Kit in the presence of cyanine 3-CTP. The labeled cRNA from each labeling reaction was hybridized to individual arrays at 65°C for 17 h using the Agilent Gene Expression Hybridization Kit. After washing, the arrays were scanned at PMT levels 350 using GenePix 4200AL scanner (Molecular Device Inc.), the Feature extraction software (V. 9.5.1) from Agilent was used to automatically find and place microarray grids, reject outlier pixels, accurately determine feature intensities and ratios, flag outlier pixels, and calculate statistical confidences.
Microarray data analysis
Microarray data analyses were processed with GeneSpring version 7.0 and 10.0. The sample quality control was based on the Pearson correlation of a sample with other samples in the whole experiment. If the average Pearson correlation with other samples was less than 80%, the sample was excluded for further analysis. If the scanned intensity was less than 5.0 for a probe, it was transformed to 5. A perchip (within) array normalization was performed using 50 percentile values of all the probe values in the array. Per gene (between) array normalization was also applied using the median value of a gene across all samples in the experiment. Probe features were first filtered using flags. A "present" or "absent" flag was defined using the Agilent Feature Extraction 9.5.1 software. Only a probe that had present flags in at least 50% samples of all the arrays was kept for further analyses. Data were subsequently log (base 2) transformed for statistical analyses.
To identify differentiated genes that responded to TNT treatment with various doses in vivo, One-Way ANOVA was performed across 5 doses at each time point (24 h or 48h). A cut off p value 0.05 was used to find statistically changed genes, and a 1.5 fold change showing at least one pair of doses at each time point was further applied to filter less significantly regulated genes. To idenitfy differential genes after TNT treatment in vitro, an un-paired t-test with cut off p value 0.05 was applied to compare control samples and TNT exposed samples. In additon, 1.5 fold change was applied to identify more signifcantly regulated genes.
Gene functional analysis and pathway analysis
Significantly regulated probes were employed for two-way hierarchical clustering (clustering both genes and samples) using GeneSpring 7.0. A Pearson correlation with average linkage was applied for the clustering. Gene functional categories were classified according to the Ingenuity knowledge base tool. A Gene functional term enrichment p value less than 0.05 was considered significant. Pathway analysis was performed using the Ingenuity canonical pathways analysis tool. Similar to functional term analysis, a pathway with an enrichment p value less than 0.05 was considered to be a significantly regulated pathway.
Reverse-transcription quantitative PCR (QRT-PCR)
Two-stage RT-QPCR were performed, 1000 ng of total RNA were first reverse transcribed into cDNA in a 20 μl reaction containing 250 ng random primers and SuperScript™ III reverse transcriptase (Invitrogen) following the manufacturer's instruction. The synthesized cDNA was diluted to 10 ng/μl as cDNA template. QPCR was performed on ABI Sequence Detector 7900. Each 20-μl reaction was run in duplicate and contained 6 μl (10ng/μl) of synthesized cDNA templates and 3 μl of nuclease-free water along with 1 μl of TaqMan gene specific assay and 10ul of 2× TaqMan universal PCR Master Mix (ABI). Cycling parameters were 95°C for 15 min to activate the DNA polymerase, then 40 cycles of 95°C for 15 s and 60°C for 1 min.
Inferring gene networks using a reverse engineering algorithm
Az and Bz are the z-score based on of A's and B's MI score distribution respectively between gene A and gene B . Ten bins were used for the binning parameters of CLR, with a spline degree of 3. A cut-off z-score was set to 2.0 so that only the most significant interaction edges could be counted. Networks were visualized by Cytoscape [68, 69].
We thank Dr. Sharon Meyer from University of Louisiana at Monroe for supervising animal experiments. The use of trade, product, or firm names in this report is for descriptive purposes only and does not imply endorsement by the U.S. Government. The findings of this report are not to be construed as an official Department of the Army position unless so designated by other authorized documents. The tests described and the resulting data presented herein, unless otherwise noted, were supported by research under the Environmental Quality Technology Program (contract #W912HZ-05-P-0145) of the United States Army Engineer Research and Development Center. Permission was granted by the Chief of Engineers to publish this information.
- Kirkland D, Aardema M, Henderson L, Muller L: Evaluation of the ability of a battery of three in vitro genotoxicity tests to discriminate rodent carcinogens and non-carcinogens I. Sensitivity, specificity and relative predictivity. Mutat Res. 2005, 584: 1-256.View ArticlePubMedGoogle Scholar
- Elferink MG, Olinga P, Draaisma AL, Merema MT, Bauerschmidt S, Polman J, et al.: Microarray analysis in rat liver slices correctly predicts in vivo hepatotoxicity. Toxicol Appl Pharmacol. 2008, 229: 300-309. 10.1016/j.taap.2008.01.037View ArticlePubMedGoogle Scholar
- Bakand S, Winder C, Khalil C, Hayes A: Toxicity assessment of industrial chemicals and airborne contaminants: transition from in vivo to in vitro test methods: a review. Inhal Toxicol. 2005, 17: 775-787. 10.1080/08958370500225240View ArticlePubMedGoogle Scholar
- Blumenthal RD, Goldenberg DM: Methods and goals for the use of in vitro and in vivo chemosensitivity testing. Mol Biotechnol. 2007, 35: 185-197. 10.1007/BF02686104View ArticlePubMedGoogle Scholar
- Navas JM, Segner H: Vitellogenin synthesis in primary cultures of fish liver cells as endpoint for in vitro screening of the (anti)estrogenic activity of chemical substances. Aquat Toxicol. 2006, 80: 1-22. 10.1016/j.aquatox.2006.07.013View ArticlePubMedGoogle Scholar
- Whong WZ, Edwards GS: Genotoxic activity of nitroaromatic explosives and related compounds in Salmonella typhimurium. Mutat Res. 1984, 136: 209-215. 10.1016/0165-1218(84)90054-5View ArticlePubMedGoogle Scholar
- Wintz H, Yoo LJ, Loguinov A, Wu YY, Steevens JA, Holland RD, et al.: Gene expression profiles in fathead minnow exposed to 2, 4-DNT: correlation with toxicity in mammals. Toxicol Sci. 2006, 94: 71-82. 10.1093/toxsci/kfl080View ArticlePubMedGoogle Scholar
- Levine BS, Furedi EM, Gordon DE, Barkley JJ, Lish PM: Toxic interactions of the munitions compounds TNT and RDX in F344 rats. Fundam Appl Toxicol. 1990, 15: 373-380. 10.1016/0272-0590(90)90062-OView ArticlePubMedGoogle Scholar
- Sabbioni G, Sepai O, Norppa H, Yan H, Hirvonen A, Zheng Y, et al.: Comparison of biomarkers in workers exposed to 2, 4, 6-trinitrotoluene. Biomarkers. 2007, 12: 21-37. 10.1080/13547500600807012View ArticlePubMedGoogle Scholar
- Sabbioni G, Liu YY, Yan H, Sepai O: Hemoglobin adducts, urinary metabolites and health effects in 2, 4, 6-trinitrotoluene exposed workers. Carcinogenesis. 2005, 26: 1272-1279. 10.1093/carcin/bgi078View ArticlePubMedGoogle Scholar
- Sabbioni G, Wei J, Liu YY: Determination of hemoglobin adducts in workers exposed to 2, 4, 6-trinitrotoluene. J Chromatogr B Biomed Appl. 1996, 682: 243-248. 10.1016/0378-4347(96)00083-7View ArticlePubMedGoogle Scholar
- Gong P, Guan X, Inouye LS, Pirooznia M, Indest KJ, Athow RS, et al.: Toxicogenomic analysis provides new insights into molecular mechanisms of the sublethal toxicity of 2, 4, 6-trinitrotoluene in Eisenia fetida. Environ Sci Technol. 2007, 41: 8195-8202. 10.1021/es0716352View ArticlePubMedGoogle Scholar
- Gandia-Herrero F, Lorenz A, Larson T, Graham IA, Bowles DJ, Rylott EL, et al.: Detoxification of the explosive 2, 4, 6-trinitrotoluene in Arabidopsis: discovery of bifunctional O- and C-glucosyltransferases. Plant J. 2008, 56: 963-974. 10.1111/j.1365-313X.2008.03653.xView ArticlePubMedGoogle Scholar
- Eisentraeger A, Reifferscheid G, Dardenne F, Blust R, Schofer A: Hazard characterization and identification of a former ammunition site using microarrays, bioassays, and chemical analysis. Environ Toxicol Chem. 2007, 26: 634-646. 10.1897/06-285R.1View ArticlePubMedGoogle Scholar
- Kumagai Y, Wakayama T, Lib S, Shinohara A, Iwamatsu A, Sun G, et al.: Zeta-crystallin catalyzes the reductive activation of 2, 4, 6-trinitrotoluene to generate reactive oxygen species: a proposed mechanism for the induction of cataracts. FEBS Lett. 2000, 478: 295-298. 10.1016/S0014-5793(00)01841-XView ArticlePubMedGoogle Scholar
- Perkins EJ, Bao W, Guan X, Ang CY, Wolfinger RD, Chu TM, et al.: Comparison of transcriptional responses in liver tissue and primary hepatocyte cell cultures after exposure to hexahydro-1, 3, 5-trinitro-1, 3, 5-triazine. BMC Bioinformatics. 2006, 7 (Suppl 4): S22- 10.1186/1471-2105-7-S4-S22PubMed CentralView ArticlePubMedGoogle Scholar
- Chow JP, Siu WY, Fung TK, Chan WM, Lau A, Arooz T, et al.: DNA damage during the spindle-assembly checkpoint degrades CDC25A, inhibits cyclin-CDC2 complexes, and reverses cells to interphase. Mol Biol Cell. 2003, 14: 3989-4002. 10.1091/mbc.E03-03-0168PubMed CentralView ArticlePubMedGoogle Scholar
- Chow JP, Siu WY, Ho HT, Ma KH, Ho CC, Poon RY: Differential contribution of inhibitory phosphorylation of CDC2 and CDK2 for unperturbed cell cycle control and DNA integrity checkpoints. J Biol Chem. 2003, 278: 40815-40828. 10.1074/jbc.M306683200View ArticlePubMedGoogle Scholar
- Pei L, Melmed S: Isolation and characterization of a pituitary tumor-transforming gene (PTTG). Mol Endocrinol. 1997, 11: 433-441. 10.1210/me.11.4.433View ArticlePubMedGoogle Scholar
- Zou H, McGarry TJ, Bernal T, Kirschner MW: Identification of a vertebrate sister-chromatid separation inhibitor involved in transformation and tumorigenesis. Science. 1999, 285: 418-422. 10.1126/science.285.5426.418View ArticlePubMedGoogle Scholar
- Kaafarani I, Fernandez-Sauze S, Berenguer C, Chinot O, Delfino C, Dussert C, et al.: Targeting adrenomedullin receptors with systemic delivery of neutralizing antibodies inhibits tumor angiogenesis and suppresses growth of human tumor xenografts in mice. FASEB J. 2009, 23: 3424-3435. 10.1096/fj.08-127852View ArticlePubMedGoogle Scholar
- Park JW, Voss PG, Grabski S, Wang JL, Patterson RJ: Association of galectin-1 and galectin-3 with Gemin4 in complexes containing the SMN protein. Nucleic Acids Res. 2001, 29: 3595-3602. 10.1093/nar/29.17.3595PubMed CentralView ArticlePubMedGoogle Scholar
- Danial NN, Korsmeyer SJ: Cell death: critical control points. Cell. 2004, 116: 205-219. 10.1016/S0092-8674(04)00046-7View ArticlePubMedGoogle Scholar
- Bersell K, Arab S, Haring B, Kuhn B: Neuregulin1/ErbB4 signaling induces cardiomyocyte proliferation and repair of heart injury. Cell. 2009, 138: 257-270. 10.1016/j.cell.2009.04.060View ArticlePubMedGoogle Scholar
- Britsch S: The neuregulin-I/ErbB signaling system in development and disease. Adv Anat Embryol Cell Biol. 2007, 190: 1-65. full_textView ArticlePubMedGoogle Scholar
- Chang H, Riese DJ, Gilbert W, Stern DF, McMahan UJ: Ligands for ErbB-family receptors encoded by a neuregulin-like gene. Nature. 1997, 387: 509-512. 10.1038/387509a0View ArticlePubMedGoogle Scholar
- Kersten S, Desvergne B, Wahli W: Roles of PPARs in health and disease. Nature. 2000, 405: 421-424. 10.1038/35013000View ArticlePubMedGoogle Scholar
- Auwerx J: PPARgamma, the ultimate thrifty gene. Diabetologia. 1999, 42: 1033-1049. 10.1007/s001250051268View ArticlePubMedGoogle Scholar
- Rasooly R, Schuster GU, Gregg JP, Xiao JH, Chandraratna RA, Stephensen CB: Retinoid × receptor agonists increase bcl2a1 expression and decrease apoptosis of naive T lymphocytes. J Immunol. 2005, 175: 7916-7929.View ArticlePubMedGoogle Scholar
- Lachance B, Robidoux PY, Hawari J, Ampleman G, Thiboutot S, Sunahara GI: Cytotoxic and genotoxic effects of energetic compounds on bacterial and mammalian cells in vitro. Mutat Res. 1999, 444: 25-39.View ArticlePubMedGoogle Scholar
- Padgham CR, Paine AJ: Altered expression of cytochrome P-450 mRNAs, and potentially of other transcripts encoding key hepatic functions, are triggered during the isolation of rat hepatocytes. Biochem J. 1993, 289 (Pt 3): 621-624.PubMed CentralView ArticlePubMedGoogle Scholar
- Hahn ME: Biomarkers and bioassays for detecting dioxin-like compounds in the marine environment. Sci Total Environ. 2002, 289: 49-69. 10.1016/S0048-9697(01)01016-6View ArticlePubMedGoogle Scholar
- Gonzalez FJ: Transgenic models in xenobiotic metabolism and toxicology. Toxicology. 2002, 181-182: 237-239. 10.1016/S0300-483X(02)00288-3View ArticlePubMedGoogle Scholar
- Hawkins SA, Billiard SM, Tabash SP, Brown RS, Hodson PV: Altering cytochrome P4501A activity affects polycyclic aromatic hydrocarbon metabolism and toxicity in rainbow trout (Oncorhynchus mykiss). Environ Toxicol Chem. 2002, 21: 1845-1853.View ArticlePubMedGoogle Scholar
- Morel Y, de W, Barouki R: A repressive cross-regulation between catalytic and promoter activities of the CYP1A1 and CYP2E1 genes: role of H(2)O(2). Mol Pharmacol. 2000, 57: 1158-1164.PubMedGoogle Scholar
- Boverhof DR, Burgoon LD, Tashiro C, Chittim B, Harkema JR, Jump DB, et al.: Temporal and dose-dependent hepatic gene expression patterns in mice provide new insights into TCDD-Mediated hepatotoxicity. Toxicol Sci. 2005, 85: 1048-1063. 10.1093/toxsci/kfi162View ArticlePubMedGoogle Scholar
- Naspinski C, Gu X, Zhou GD, Mertens-Talcott SU, Donnelly KC, Tian Y: Pregnane × receptor protects HepG2 cells from BaP-induced DNA damage. Toxicol Sci. 2008, 104: 67-73. 10.1093/toxsci/kfn058View ArticlePubMedGoogle Scholar
- Bodreddigari S, Jones LK, Egner PA, Groopman JD, Sutter CH, Roebuck BD, et al.: Protection against aflatoxin B1-induced cytotoxicity by expression of the cloned aflatoxin B1-aldehyde reductases rat AKR7A1 and human AKR7A3. Chem Res Toxicol. 2008, 21: 1134-1142. 10.1021/tx7004458View ArticlePubMedGoogle Scholar
- Lassen N, Bateman JB, Estey T, Kuszak JR, Nees DW, Piatigorsky J, et al.: Multiple and additive functions of ALDH3A1 and ALDH1A1: cataract phenotype and ocular oxidative damage in Aldh3a1(-/-)/Aldh1a1(-/-) knock-out mice. J Biol Chem. 2007, 282: 25668-25676. 10.1074/jbc.M702076200PubMed CentralView ArticlePubMedGoogle Scholar
- Lind C, Hochstein P, Ernster L: DT-diaphorase as a quinone reductase: a cellular control device against semiquinone and superoxide radical formation. Arch Biochem Biophys. 1982, 216: 178-185. 10.1016/0003-9861(82)90202-8View ArticlePubMedGoogle Scholar
- Thor H, Smith MT, Hartzell P, Bellomo G, Jewell SA, Orrenius S: The metabolism of menadione (2-methyl-1, 4-naphthoquinone) by isolated hepatocytes. A study of the implications of oxidative stress in intact cells. J Biol Chem. 1982, 257: 12419-12425.PubMedGoogle Scholar
- Anwar A, Dehn D, Siegel D, Kepa JK, Tang LJ, Pietenpol JA, et al.: Interaction of human NAD(P)H:quinone oxidoreductase 1 (NQO1) with the tumor suppressor protein p53 in cells and cell-free systems. J Biol Chem. 2003, 278: 10368-10373. 10.1074/jbc.M211981200View ArticlePubMedGoogle Scholar
- Keating AF, Rajapaksa KS, Sipes IG, Hoyer PB: Effect of CYP2E1 gene deletion in mice on expression of microsomal epoxide hydrolase in response to VCD exposure. Toxicol Sci. 2008, 105: 351-359. 10.1093/toxsci/kfn136PubMed CentralView ArticlePubMedGoogle Scholar
- Harris SJ, Parry RV, Westwick J, Ward SG: Phosphoinositide lipid phosphatases: natural regulators of phosphoinositide 3-kinase signaling in T lymphocytes. J Biol Chem. 2008, 283: 2465-2469. 10.1074/jbc.R700044200View ArticlePubMedGoogle Scholar
- Johansson LE, Johansson LM, Danielsson P, Norgren S, Johansson S, Marcus C, et al.: Genetic variance in the adiponutrin gene family and childhood obesity. PLoS One. 2009, 4: e5327- 10.1371/journal.pone.0005327PubMed CentralView ArticlePubMedGoogle Scholar
- Karmakar S, Cummings RD, McEver RP: Contributions of Ca2+ to galectin-1-induced exposure of phosphatidylserine on activated neutrophils. J Biol Chem. 2005, 280: 28623-28631. 10.1074/jbc.M414140200View ArticlePubMedGoogle Scholar
- Tagawa N, Yuda R, Kubota S, Wakabayashi M, Yamaguchi Y, Kiyonaga D, et al.: 17Beta-estradiol inhibits 11beta-hydroxysteroid dehydrogenase type 1 activity in rodent adipocytes. J Endocrinol. 2009, 202: 131-139. 10.1677/JOE-09-0021View ArticlePubMedGoogle Scholar
- Bujalska IJ, Gathercole LL, Tomlinson JW, Darimont C, Ermolieff J, Fanjul AN, et al.: A novel selective 11beta-hydroxysteroid dehydrogenase type 1 inhibitor prevents human adipogenesis. J Endocrinol. 2008, 197: 297-307. 10.1677/JOE-08-0050PubMed CentralView ArticlePubMedGoogle Scholar
- Rawat A, Gust KA, Deng Y, Garcia-Reyero N, Quinn MJ, Johnson MS, et al.: From raw materials to validated system: The construction of a genomic library and microarray to interpret systemic perturbations in Northern bobwhite. Physiol Genomics. 2010, 7;42 (2): 219-35.View ArticleGoogle Scholar
- Grundler R, Brault L, Gasser C, Bullock AN, Dechow T, Woetzel S, et al.: Dissection of PIM serine/threonine kinases in FLT3-ITD-induced leukemogenesis reveals PIM1 as regulator of CXCL12-CXCR4-mediated homing and migration. J Exp Med. 2009, 206: 1957-1970. 10.1084/jem.20082074PubMed CentralView ArticlePubMedGoogle Scholar
- Mandal M, Borowski C, Palomero T, Ferrando AA, Oberdoerffer P, Meng F, et al.: The BCL2A1 gene as a pre-T cell receptor-induced regulator of thymocyte survival. J Exp Med. 2005, 201: 603-614. 10.1084/jem.20041924PubMed CentralView ArticlePubMedGoogle Scholar
- Kitamoto S, Sukhova GK, Sun J, Yang M, Libby P, Love V, et al.: Cathepsin L deficiency reduces diet-induced atherosclerosis in low-density lipoprotein receptor-knockout mice. Circulation. 2007, 115: 2065-2075. 10.1161/CIRCULATIONAHA.107.688523View ArticlePubMedGoogle Scholar
- Kurosaki T: Genetic analysis of B cell antigen receptor signaling. Annu Rev Immunol. 1999, 17: 555-592. 10.1146/annurev.immunol.17.1.555View ArticlePubMedGoogle Scholar
- Coppola V, Barrick CA, Bobisse S, Rodriguez-Galan MC, Pivetta M, Reynolds D, et al.: The scaffold protein Cybr is required for cytokine-modulated trafficking of leukocytes in vivo. Mol Cell Biol. 2006, 26: 5249-5258. 10.1128/MCB.02473-05PubMed CentralView ArticlePubMedGoogle Scholar
- Hoffmann I, Clarke PR, Marcote MJ, Karsenti E, Draetta G: Phosphorylation and activation of human cdc25-C by cdc2--cyclin B and its involvement in the self-amplification of MPF at mitosis. EMBO J. 1993, 12: 53-63.PubMed CentralPubMedGoogle Scholar
- Donaldson MM, Tavares AA, Hagan IM, Nigg EA, Glover DM: The mitotic roles of Polo-like kinase. J Cell Sci. 2001, 114: 2357-2358.PubMedGoogle Scholar
- Donaldson MM, Tavares AA, Ohkura H, Deak P, Glover DM: Metaphase arrest with centromere separation in polo mutants of Drosophila. J Cell Biol. 2001, 153: 663-676. 10.1083/jcb.153.4.663PubMed CentralView ArticlePubMedGoogle Scholar
- Robidoux PY, Svendsen C, Sarrazin M, Hawari J, Thiboutot S, Ampleman G, et al.: Evaluation of tissue and cellular biomarkers to assess 2, 4, 6-trinitrotoluene (TNT) exposure in earthworms: effects-based assessment in laboratory studies using Eisenia andrei. Biomarkers. 2002, 7: 306-321. 10.1080/13547500210127327View ArticlePubMedGoogle Scholar
- Robidoux PY, Svendsen C, Sarrazin M, Thiboutot S, Ampleman G, Hawari J, et al.: Assessment of a 2, 4, 6-trinitrotoluene-contaminated site using Aporrectodea rosea and Eisenia andrei in mesocosms. Arch Environ Contam Toxicol. 2005, 48: 56-67. 10.1007/s00244-004-0217-7View ArticlePubMedGoogle Scholar
- Gong P, Guan X, Inouye LS, Deng Y, Pirooznia M, Perkins EJ: Transcriptomic analysis of RDX and TNT interactive sublethal effects in the earthworm Eisenia fetida. BMC Genomics. 2008, 9 (Suppl 1): S15- 10.1186/1471-2164-9-S1-S15PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Q, Pi J, Woods CG, Andersen ME: Phase I to II cross-induction of xenobiotic metabolizing enzymes: a feedforward control mechanism for potential hormetic responses. Toxicol Appl Pharmacol. 2009, 237: 345-356. 10.1016/j.taap.2009.04.005PubMed CentralView ArticlePubMedGoogle Scholar
- Adams M, Jia Z: Structural and biochemical analysis reveal pirins to possess quercetinase activity. J Biol Chem. 2005, 280: 28675-28682. 10.1074/jbc.M501034200View ArticlePubMedGoogle Scholar
- Persson B, Kallberg Y, Bray JE, Bruford E, Dellaporta SL, Favia AD, et al.: The SDR (short-chain dehydrogenase/reductase and related enzymes) nomenclature initiative. Chem Biol Interact. 2009, 178: 94-98. 10.1016/j.cbi.2008.10.040PubMed CentralView ArticlePubMedGoogle Scholar
- Kwak MK, Wakabayashi N, Itoh K, Motohashi H, Yamamoto M, Kensler TW: Modulation of gene expression by cancer chemopreventive dithiolethiones through the Keap1-Nrf2 pathway. Identification of novel gene clusters for cell survival. J Biol Chem. 2003, 278: 8135-8145. 10.1074/jbc.M211898200View ArticlePubMedGoogle Scholar
- National Research Council: Guide for the care and use of laboratory animals. 1996, National Academy Press, Washington D.CGoogle Scholar
- Faith JJ, Hayete B, Thaden JT, Mogno I, Wierzbowski J, Cottarel G, et al.: Large-scale mapping and validation of Escherichia coli transcriptional regulation from a compendium of expression profiles. PLoS Biol. 2007, 5: e8- 10.1371/journal.pbio.0050008PubMed CentralView ArticlePubMedGoogle Scholar
- Taylor RC, cquaah-Mensah G, Singhal M, Malhotra D, Biswal S: Network inference algorithms elucidate Nrf2 regulation of mouse lung oxidative stress. PLoS Comput Biol. 2008, 4: e1000166-PubMed CentralView ArticlePubMedGoogle Scholar
- Emig D, Cline MS, Lengauer T, Albrecht M: Integrating expression data with domain interaction networks. Bioinformatics. 2008, 24: 2546-2548. 10.1093/bioinformatics/btn437PubMed CentralView ArticlePubMedGoogle Scholar
- Cline MS, Smoot M, Cerami E, Kuchinsky A, Landys N, Workman C, et al.: Integration of biological networks and gene expression data using Cytoscape. Nat Protoc. 2007, 2: 2366-2382. 10.1038/nprot.2007.324PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.