Biomolecules & Therapeutics 2017; 25(2): 112-121
Expression Levels of GABA-A Receptor Subunit Alpha 3, Gabra3 and Lipoprotein Lipase, Lpl Are Associated with the Susceptibility to Acetaminophen-Induced Hepatotoxicity
Minjeong Kim1,†, Jun-Won Yun2,†, Kyeho Shin3, Yejin Cho4, Mijeong Yang4, Ki Taek Nam4,*, and Kyung-Min Lim1,*
1College of Pharmacy, Ewha Womans University, Seoul 03760, Republic of Korea, 2Department of Experimental Animal Research, Biomedical Research Institute, Seoul National University Hospital, Seoul 03080, 3Department of Beauty Coordination, Suwon Science College, Suwon 18516, Republic of Korea, 4Severance Biomedical Science Institute, Brain Korea 21 PLUS Project for Medical Science, Yonsei University College of Medicine, Seoul 03722, Republic of Korea
*Corresponding Authors E-mail: (Lim KM), (Nam KT) Tel: +82-2-3277-3055 (Lim KM), +82-2-2228-0754 (Nam KT) Fax: +82-2-3277-3760 (Lim KM), +82-2-2227-8129 (Nam KT)
Received: April 5, 2016; Revised: May 9, 2016; Accepted: May 11, 2016; Published online: August 19, 2016.
© The Korean Society of Applied Pharmacology. All rights reserved.

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Drug-induced liver injury (DILI) is the serious and fatal drug-associated adverse effect, but its incidence is very low and individual variation in severity is substantial. Acetaminophen (APAP)-induced liver injury accounts for >50% of reported DILI cases but little is known for the cause of individual variations in the severity. Intrinsic genetic variation is considered a key element but the identity of the genes was not well-established. Here, pre-biopsy method and microarray technique was applied to uncover the key genes for APAP-induced liver injury in mice, and a cause and effect experiment employing quantitative real-time PCR was conducted to confirm the correlation between the uncovered genes and APAP-induced hepatotoxicity. We identified the innately and differentially expressed genes of mice susceptible to APAP-induced hepatotoxicity in the pre-biopsied liver tissue before APAP treatment through microarray analysis of the global gene expression profiles (Affymetrix GeneChip® Mouse Gene 1.0 ST for 28,853 genes). Expression of 16 genes including Gdap10, Lpl, Gabra3 and Ccrn4l were significantly different (t-test: FDR <10%) more than 1.5 fold in the susceptible animals than resistant. To confirm the association with the susceptibility to APAP-induced hepatotoxicity, another set of animals were measured for the expression level of selected 4 genes (higher two and lower two genes) in the liver pre-biopsy and their sensitivity to APAP-induced hepatotoxicity was evaluated by post hoc. Notably, the expressions of Gabra3 and Lpl were significantly correlated with the severity of liver injury (p<0.05) demonstrating that these genes may be linked to the susceptibility to APAP-induced hepatotoxicity.

Keywords: Acetaminophen, Hepatotoxicity, Toxicogenomics, GABA-A receptor subunit alpha 3, Lipoprotein lipase

Acetaminophen (APAP) is one of the most widely used analgesics and antipyretics, and it is considered safe to use within the maximum daily dose of 3,250 mg and the maximum single dose of 650 mg (Krenzelok, 2009). However, its overuse or co-existing risk factors that can exacerbate liver function such as narcotic, alcohol drinking or polypharmacy, often cause drug-induced liver injury (DILI) (Watson et al., 2004; Larson et al., 2005). DILI, namely hepatotoxicity induced by drugs is one of the most frequent and severest forms of adverse drug reaction (Lazarou et al., 1998). The incidence of DILI is low varying from 1/100,000 to 1/10,000 but, its impact is huge since acute liver failure with high mortality can occur (Chang and Schiano, 2007; Stephens et al., 2012). More than 1,000 drugs were known to cause DILI in human (Lewis, 2000; Reuben et al., 2010). Of these APAP-induced DILI accounts for 46% of cases of fulminant liver failure in the United States (Ostapowicz et al., 2002). More than 200 million persons take APAP annually and among them about 200 persons die of fulminant hepatic failure from over-dosage (Park, 2006).

APAP induces hepatotoxicity via the generation of reactive and electrophilic metabolites like N-acetyl-p-benzoquinone imine (NAPQI) by cytochrome P450 that can result in glutathione (GSH)-adduct formation, covalent binding to vital endogenous macromolecules, and depletion of cellular antioxidant reserve (Bessems and Vermeulen, 2001). N-acetylcysteine given (NAC) within 12 h seems to be effective in protecting most patients against severe liver damage induced by APAP (Prescott et al., 1977) but frequently APAP-induced DILI could be recognized only after symptoms have developed by when NAC is no longer effective, reflecting that APAP-induced hepatotoxicity progresses in multiple stages and diverse factors may be involved that include inflammation, metabolism and lipid synthesis (Hinson et al., 2010). Moreover, APAP-induced DILI appears with variable degrees of severity ranging from minimal increases in the levels of serum alanine transaminase (ALT)/aspartate transaminase (AST) to severe hepatic necrosis and fatal hepatic failure (Larson et al., 2005; Watkins et al., 2006; Yun et al., 2014), reflecting the existence of “individual factors” in the manifestation of APAP-induced hepatotoxicity.

To identify factors that determine the severity of APAP-induced hepatotoxicity or to develop a biomarker to screen susceptible individuals, various approaches have been attempted. Liu et al. (2010) have analyzed hepatic gene expressions related to APAP toxicity using resistant strain SJL/J and three sensitive strains C57BL/6J, DBA/2J, and SM/J. Stamper et al. (2010) have found several genes with significant differences in the expression between non-toxic APAP regio-isomer 3-hydroxyacetanilide and APAP in TGF-α transgenic mouse hepatocytes. Also, Umbright et al. (2010) have reported that many blood genes associated with inflammation, immune, stress responses and energy metabolism were statistically different in their expression levels following APAP treatment.

While these approaches are effective in the identification of target molecules or markers that can diagnose the extent of APAP-induced hepatotoxicity, their applicability as biomarkers to screen-out susceptible individuals is limited since the exposure to APAP might have extensively altered the genetic landscape due to direct or collateral tissue damages and subsequent inflammatory responses. Recently, to draw an intact genetic picture without the interference from APAP-induced toxicity, we examined through the comparison of pre-dose blood genes and the severity of APAP-induced hepatotoxicity posterior in rats in vivo, which demonstrated that protein kinase A (PKA) inhibitor alpha (Pkia) expression level in predose blood can be employed to predict susceptible individuals without the interference from the exposure to APAP (Yun et al., 2014). More importantly, this study has provided an insight into the role of PKA in the manifestation of APAP-induced hepatotoxicity. This approach has been further corroborated by Lu et al. (2015) who demonstrated that four genes, including Incenp, Rpgrip1, Sbf1, and Mmp12, which are associated with cell proliferation and tissue repair functions, in the blood collected from individuals prior and posterior to APAP administration can be used for identifying susceptible population to DILI. However, these studies offered surrogate genetic biomarkers in blood rather than that of direct target, liver, therefore, its utility and implication might be limited. In this study, we further investigated the transcriptome in the biopsy of liver pre-dose to uncover the genetic factors for individual susceptibility of APAP-induced hepatotoxicity in mice in vivo in an effort to provide a clue to understand the progression of, and defense mechanisms against APAP-induced DILI.



Outbred ICR mice was selected as the experimental species for this study since rats do not manifest a clear profile of blood chemistry to APAP-induced hepatotoxicity when compared with the responses to carbon tetrachloride or D-galactosamine (Shin et al., 2014). Male ICR mice aged 7 weeks were purchased from Jung-Ang Lab Animals (Seoul, Korea) and housed in a specific pathogen-free (SPF) facility of Ewha Womans University. We used ICR mice, one of the outbred strains, which retain a certain level of genetic diversity in the test population in present study. The mice were kept under controlled environmental conditions (23 ± 3°C, 40–60% relative humidity, 12/12 h dark-light cycle) with ad libitum access to laboratory normal diets (Purina Co., Seoul, Korea) and tab water. All protocol involving animals were carried out in accordance with the guidelines of the Institutional Animal Care and Use Committee of Ewha Womans University.

Experimental protocol

To identify the innate genes for individual variation of liver injury, we analyzed the gene expression profiles in the liver pre-biopsied from individual mice prior to the oral administration of APAP (Sigma, St. Louis, MO, USA) and compared them with the individual severity of liver injury after APAP administration according to the method described previously (Yun et al., 2009, 2010) (Fig. 1). In brief, the mice were randomly assigned into three groups consisting of negative control group (no biopsy group) (N=5), biopsy control group (N=5), and APAP group (N=32). The operation was performed for minimum liver biopsies (about 10 mg) from left lobe of the liver of anesthetized mice in biopsy groups (biopsy control group and APAP group) using an isoflurane vaporizer (Midmark, Orchard Park, OH, USA) with an isoflurane of 1.5% to 3% and an oxygen flow of 0.5 L/min. After conducting the sutures to close the subcutis tissue and the skin using Surgifit 6-0 (AILEE, Busan, Korea) and Black silk 6-0 (AILEE). For the assessment of postoperative recovery, ALT, AST, albumin, albumin/globulin ratio (A/G ratio), total bilirubin and bile acid were measured using serum samples from blood collected via retro-orbital plexus of anesthetized animals into a gel serum separator plain blood tubes (MiniCollect 0.8 ml Z Serum Sep, Greiner Bio-One, Frickenhausen, Germany). After recovery for 3 weeks, the mice were administered with APAP through oral gavage at a dose of 300 mg/kg (dissolved in deionized water) according to a previous method (Saha and Nandi, 2009) with minor modifications. At 24 h after administrations of APAP, we performed the biochemical analysis using the blood collected via the postcaval vein from anesthetized animals, and then conducted microarray analysis using the pre-collected liver samples of the 10 animals (5 susceptible and 5 resistant) selected based on the results of the biochemical analysis and compared the gene expressions with the individual severity of liver injury after APAP administration to identify the innate genes for individual variation of liver injury. Additionally, to identify whether the innate genes selected by microarray experiment can predict the susceptibility of liver injury, we conducted the real-time PCR analysis for selected genes using liver samples biopsied from another set of animals (N=32) and then compared the gene expression with individual severity of liver injury after APAP administration (300 mg/kg, p.o.). The severity of liver injury were analyzed with the biochemical indicators, including ALT, total bilirubin, AST, lactate dehydrogenase (LDH) and bile acid at 24 h after administrations of APAP.

RNA isolation and microarray analysis

The liver tissue biopsied before APAP administration was processed with TRIzol reagent (Invitrogen, Carlsbad, CA, USA) for isolation of total RNA. RNA precipitates were dissolved in RNase free DEPC treated water (USB, Cleveland, OH, USA). The concentration of RNA was determined Nano drop 1000 spectrophotometer (NanoDrop Technologies, Inc., Wilmington, DE, USA). Affymetrix (Santa Clara, CA, USA) GeneChip1 Mouse Gene 1.0 ST arrays were used to analyze the differential gene expressions, as described previously (Yun et al., 2009). The normalized scanned probe array data were compared between the groups to generate p-value and signal log ratio (fold change). Unpaired t-test was applied to determine statistically reliable probe sets.

Reverse transcription-polymerase chain reaction (RT-PCR)

Relative expression levels of mRNAs were measured by quantitative real-time PCR. Total RNA, extracted from liver prior to APAP treatment, was used to synthesize cDNA using pre-master mix with oligo dT (Bioepis, Seoul, Korea). Each reaction was performed using power SYBR Green PCR master mix in a 7300 real-time PCR machine (Applied Biosystems, Warrington, UK). The sequence of primers of mice liver were as follows: forward Gdap10, 5′-GCC TGG CAC GTA GAA CAA AC-3′; reverse Gdap10, 5′-ACT AGT TTG CTG CCC TCG TT-3′; forward Lpl, 5′-TCA GAG CCA AGA GAA GCA GC-3′; reverse Lpl, 5′-ATC TCG AAG GCC TGG TTG TG-3′; forward Klf10, 5′-AGC AAG GGT CAC TCC TCA GA-3′; reverse Klf10, 5′-AGG TTT TTC CCC TGT GTG TG-3′; forward Malat1, 5′-GAC AAA AGG CTA AAG TGG ATG-3′; reverse Malat1, 5′-TTG TGC TGG CTC TAC CAT T-3′; forward Ccrn4l, 5′-CCT GGA GTG CAA GGA GTC TG-3′; reverse Ccrn4l, 5′-CTG GGT GAT GTT CTG CAG GT-3′; forward Gabra3, 5′-CCA ACA GCG ATT GCT TCA CC-3′; reverse Gabra3, 5′-GGC AAA GAG CAC AGG GAA GA-3′; forward GAPDH, 5′-CTA CCC CCA ATG TGT CCG TC-3′; reverse GAPDH, 5′-AAG TCG CAG GAG ACA ACC TG-3′. Annealing temperature were 51°C. For real-time PCR data, fold change results were calculated using 2−ΔΔCt relative to the internal reference gene (GAPDH) and the mean of all samples. The change in the Ct (ΔCt) of the target genes was calculated as ΔCt=(Ct of target genes)-(Ct of GAPDH). The ratio of the target gene to the housekeeping gene was calculated and expressed as 2−ΔCt. This ratio was then used to evaluate the expression level of the target gene within each animal. To determine the fold changes in gene expression among animals, the normalized gene expression of the target genes was divided by the normalized expression of the same target gene in the sample with the lowest level of the normalized gene expression of the target genes, expressed as 2−ΔΔCt.

Statistical analysis

All data were analyzed by student t-test or one-way analysis of variance (ANOVA) followed by post hoc Tukey’s multiple range tests to determine treatment effect and to compare differences between group means. Differences were considered to be significant at p<0.05.


Assessment of the influences of liver biopsy on liver function of ICR mice

Pre-biopsy scheme needs the liver biopsy at the naïve condition prior to APAP administration to identify the innate genetic factors associated with the individual variation of APAP-induced liver injury afterwards (Fig. 1). Firstly, to check if pre-biopsy procedure does not affect normal liver function of mice, small amount of liver tissue (about 10 mg) was pre-biopsied from the identical region of left lobe and liver function was evaluated after sufficient postoperative recovery (3 weeks) by measuring ALT, AST, albumin, A/G ratio, total bilirubin, and bile acid in the serum. As a result, liver function of biopsy control group was not statistically different from those in Negative control group (Fig. 2), suggesting that the effects of pre-biopsy scheme on the liver function are minimal for mice as shown in rats previously (Yun et al., 2010).

Hepatotoxicity induced by APAP administration after liver biopsy

After liver pre-biopsy and postoperative recovery for 3 weeks, APAP was orally administered to the ICR mice (N= 32), an outbred strain at a dose of 300 mg/kg. The APAP-treated mice exhibited clear signs of liver injury, as indicated by increased serum ALT and total bilirubin (Fig. 3). Importantly, there was a substantial inter-animal variation in the level of these biochemical indicators, suggesting that mice can display different susceptibility to APAP. On the basis of ALT and total bilirubin levels, animals could be grouped into 2 groups, that is, susceptible and resistant animals (top 5, susceptible and bottom 5, resistant, Fig. 3A). Serum ALT and total bilirubin in susceptible group were substantially higher when compared with both Negative control group and resistant group (Fig. 3B, 3C). Further histological examination (Fig. 4A) and blood biochemistry (Fig. 4B) also revealed that the animals grouped as susceptible manifested higher toxic response to APAP than resistant which could be readily determined by significantly higher liver toxicity markers and massive hepatic injury including numerous apoptotic cells and inflammatory cell infiltration in the liver tissue.

Gene expression analysis with microarray to identify genes associated with susceptibility to APAP-induced hepatotoxicity

To determine the genetic factors associated with these inter-individual variations to APAP-induced hepatotoxicity, microarray analysis was performed with the liver samples of 5 susceptible and 5 resistant animals pre-obtained before APAP treatment as described in scheme (Fig. 1). The reliability of the transformed and normalized data was statistically analyzed using one-way ANOVA. This was visualized by hierarchical clustering of the calculated data from the experiment (Fig. 5A). From this analysis, 16 genes (excluding unknown sequences and noncoding genes) were found to be different with statistical reliability at p<0.05 with more than 1.5 fold difference between two groups in their expression levels (Table 1). Among them, uppermost two genes expressed higher innately in susceptible group were found to be Gdap10 and Lpl, and lowermost two genes expressed lower innately in susceptible group were Gabra3 and Ccrn4l (Fig. 5B, 5C).

Prediction of the susceptibility of animals to APAP-induced liver injury by real-time PCR analysis of the selected genes

To further confirm whether the individual expression levels of Gdap10, Lpl, Ccrn4l and Gabra3 are related to the interindividual variation in the susceptibility of APAP-induced liver injury indeed, the expression of these genes was analyzed with liver biopsy samples pre-collected from 32 animals before APAP-treatment using quantitative real-time PCR analysis. While a meaningful relationship could not be found with the expression levels of Gdap10 and Ccrn4l with the severity of APAP-induced hepatotoxicity, the innate expression level of Gabra3 matched well with the severity of APAP-induced liver injury as could be determined by the strong correlation with APAP-induced AST (Spearman Correlation Analysis, p<0.05; Correlation coefficients, −0.4209) (Fig. 6A) and ALT (Spearman Correlation Analysis, p<0.05; Correlation coefficient, −0.4012) decrease (Fig. 6B). Moreover, the innate expression level of Lpl matched well with the severity of APAP-induced liver injury as could be determined by the strong correlation with APAP-induced total bilirubin (Spearman Correlation Analysis, p<0.05; Correlation coefficient, 0.4052) increases (Fig. 6C), of which relationships were further corroborated by the comparison of AST, ALT and total bilirubin of grouping based on the genes (Fig. 6D–6F).


Here, we explored the genetic markers in the liver that are associated with the individual susceptibility to APAP-induced hepatotoxicity in ICR mice employing pre-biopsy scheme. Through this approach, we could identify that low Gabra3 and high Lpl expression may be related to the susceptibility to APAP-induced hepatotoxicity. More importantly, we could confirm that susceptible animals could roughly be predicted through determining the expression levels of Gabra3 and Lpl in liver. In addition, this study suggests that Gabra3 may play a protective role against APAP-induced hepatotoxicity in the liver while Lpl may contribute to the development or aggravation of liver damages, which warrant further studies.

In the present study, a substantial inter-individual variation on APAP-induced hepatotoxicity was found among individual mouse, a well-known species susceptible to APAP, as evidenced by the histopathological examination and biochemistry analysis. Many previous transcriptomic approaches have examined the genetic landscape after APAP treatment and hence, can have limited applicability for the prediction of the individual susceptibility to APAP-induced hepatotoxicity. We previously demonstrated (Yun et al., 2009, 2010) that the predose scheme using pre-biopsied liver, which has a unique ability to fully regenerate after injury for the maintenance of its functions on metabolism and detoxification (Fausto et al., 2006; Michalopoulos, 2007) may provide a useful tool for prescreening of the susceptible individual and for the discovery of key molecules in the manifestation of toxicity. In the current microarray experiment using the pre-biopsy method, we identified 16 genes including Gdap10, Klf10, Malat1, and Ccrn41 (unpaired t-test: FDR <10% and fold change >1.5) as candidate genes predictive of severity on APAP-induced liver injury.

Malat1, which has been known to be associated with a prognostic marker for metastasis in early stages of lung adenocarcinoma (Ji et al., 2003; Gutschner et al., 2013) and ganglioside-induced differentiation-associated-protein 10 (Gdap10), which is one of Gdap genes involved in different signal transduction pathway (Liu et al., 1999). Klf10 has been known to regulate TGF-β signaling by blocking expression of the negative regulator, Smad7 (Johnsen et al., 2002), and activating expression of the positive effector, Smad2 (Johnsen et al., 2002). Ccrn4l is a gene that encodes a circadian deadenylase, and its disruption in mice showed lower body weight and reduced visceral fat, reflecting resistance to fatty liver and diet-induced obesity (Green et al., 2007). Genes encoding cellular components like Ankrd55, Cspp1, Dennd4a, Chd7, Snord58b, Taf1d and Pdilt, are also shown to be associated. Other genes, Gas5, Dclre1c, and Meig1 also need further studies to examine their roles in the manifestation of APAP-induced hepatotoxicity.

Although DNA microarray can simultaneously quantitate the expression of thousands of genes, a second methodology including quantitative real-time RT-PCR is required to assess the accuracy of the candidate genes discovered by microarray measurements (Draghici et al., 2006). Furthermore, the margin of difference in the expression is largely small, which necessitates a verification step. To confirm if the gene expression profiles detected with the microarray is reproducible and can be used to predict the animal to be susceptible or resistant to APAP-induced hepatotoxicity indeed, the liver samples pre-biopsied from another set of animals were undergone realtime PCR analysis for the expression levels of 4 genes including higher two and lower two genes selected from microarray analysis.

Among four candidate genes, it is notable that the innate expression of Lpl in the liver was significantly correlated well with total bilirubin which was well-correlated with the severity of APAP-induced hepatotoxicity (Spearman correlation coefficients 0.4052, p<0.05). The expression of Lpl in liver was expressed higher innately in susceptible group to APAP-induced hepatotoxicity when compared to resistant group. Lpl, which is a member of the lipase superfamily that includes pancreatic, hepatic and endothelial lipase, is widely expressed in many tissues such as liver, brain, heart and adipose tissues (Wang and Eckel, 2009). Kim et al. (2001) have demonstrated that the overexpression of Lpl, which has been known to be the rate-limiting enzyme involved with triglyceride hydrolysis (Goldberg, 1996), causes profound insulin resistance in liver (Baron et al., 1988). Insulin resistance can lead to selective accumulation of fatty acid-derived metabolites (i.e., fatty acyl CoA, ceramide, diacylglycerol) in liver and may be implicated in the development of acute liver failure (ALF) (Clark et al., 2001) through impaired peripheral glucose utilization and a failure to fully suppress endogenous glucose production, contributing to the catabolic state that occurs in ALF. These suggest potential involvement of Lpl in the susceptibility to hepatotoxicity although little is known currently for the exact mechanism on link with the APAP-induced hepatotoxicity.

We also observed that Gabra3 expression level was found to be correlated well with the APAP-induced hepatotoxicity as shown by significant Spearman correlation coefficient of −0.4209 (p<0.05) and −0.4012 (p<0.05) for AST and ALT increase, respectively. Gabra3 has been originally found in the CNS but its expression in the liver tissue has been reported (Moe et al., 2008; Oh et al., 2009). Moreover, Biju et al. (2001) have showed that hepatic GABA-A receptor can give inhibitory signal for hepatic cell proliferation. Bozogluer et al. (2012) also showed that flumazenil, a GABA-A receptor antagonist, attenuated the APAP-induced hepatotoxicity. Conversely, the expression of Gabra3 was down-regulated following the exposure to a hepatotoxicant, 4,4′-methylenedianiline in mice (Oh et al., 2009). Despite many contradicting results on the role of GABA-A receptor in the liver, we could speculate that the susceptibility to APAP-induced hepatotoxicity in the animals with innately lower expression of Gabra3 can be attributable to the altered signaling for hepatocyte proliferation which is important for the recuperation processes.

Previously, we demonstrated that high expression level of Pkia in pre-dose blood may be related to the susceptibility to APAP-induced hepatotoxicity (Yun et al., 2014). In our microarray analysis, however, the difference of Pkia expression between susceptible and resistant was only marginal (1.12 fold higher in susceptible with p-value of 0.055, t-test, data not shown) which may be considered confirmatory but unapparent. This discrepancy may be from different species employed for the studies (SD rats for blood genes, ICR mice for liver genes). Here, we employed ICR mice since this species manifests clear blood chemistry profiles for APAP-induced hepatotoxicity in contrast to the rat which needs histology to monitor APAP-induced hepatotoxicity additional to blood chemistry. These results suggest that other experimental animal species or strains may produce different genetic markers owing to distinct physiology and translational clinical research in human is ultimately necessary to draw the solid answers.

In conclusion, we demonstrated that the two gene biomarkers in the liver pre-biopsied prior to administration were related to the inter-individual variation in the severity of APAP-induced hepatotoxicity although other important factors like different gastrointestinal absorption of APAP (Sanaka et al., 1998) or alteration of metabolic capacity like CYP2E1 (Lee et al., 1996) could not be examined due to the limitation of study design and animal species employed. However, the pre-dose expressions of Lpl and Gabra3 in the liver were all correlated well with the post-dose changes of ALT, AST, and total bilirubin. Accordingly, the data presented in this study suggest a novel role of Lpl and Gabra3 in the liver in the manifestation of APAP-induced hepatotoxicity. In addition, we could speculate that these genes can be employed to screen out individuals susceptible to APAP-induced hepatotoxicity although further studies should be conducted.


This study was supported by the Korea Mouse Phenotyping Center, and the Brain Korea 21 PLUS Project for Medical Science, Yonsei University, the Bio and Medical Technology Development Program of the National Research Foundation (NRF) funded by the Ministry of Science, ICT and Future Planning (NRF-2013M3A9D5072551), and Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2014R1A1A2058389) to K.T.N.

Fig. 1. Experimental scheme for pre-biopsy. To identify the innate genes for individual variation of liver injury, we analyzed the gene expression profiles in the liver pre-biopsied from individual mice prior to APAP administration and compared them with the individual severity of liver injury after APAP administration. ALT: alanine aminotransferase, T.bilirubin: Total Bilirubin, Sus: susceptible, Res: resistant.
Fig. 2. Effects of the liver biopsy on serum biochemistry after recovery for 3 weeks. Levels of ALT, AST, albumin, A/G ratio, total bilirubin and bile acid in the serum of mice recovered from or not undergone pre-biopsy. The data are presented as means ± SE. NC, Negative control group; Biopsy, Biopsy control group.
Fig. 3. Selection of the groups based on the serum ALT and total bilirubin after APAP administration. (A) Selection of the samples for microarray analysis (N=5/group). (B) ALT and total bilirubin levels of selected groups. The data are presented as means ± SE (*p<0.05 and **p<0.01). Con, Negative control group; Res, resistant group; Sus, susceptible group.
Fig. 4. Biochemical and histological analysis between groups after APAP treatment. The mice were orally administrated with APAP at a dose of 300 mg/kg. At 24 h after APAP administration, histological examination (A) and blood biochemistry for lactate dehydrogenase (LDH, U/L), aspartate aminotransferase (AST, U/L)) and bile acid (μmol/L) (B) were conducted. Representative H&E stained liver tissue (bar = 200 μm). The data are presented as means ± SE. *Significantly different from Con group (p<0.05). Con, Negative control group; Res, resistant group; Sus, susceptible group.
Fig. 5. Hierarchical clustering that exhibited altered expression (p<0.05, fold change >1.5) of the microarray assay. (A) Comparison between the 5 most susceptible mice and the 5 resistant mice. The x-axis of the dendrogram represents the gene symbol. The y-axis represents the severity of liver injury. Colors range from bright green (log2 value −1) to bright red (log2 value +1). (B, C) The difference of gene expressions (intensity) in microarray analysis between susceptible and resistance group.
Fig. 6. Verification test with a new set of animals (Correlation between expression of candidate genes and serum chemistry, N=32). The expression levels of Gabra3 and Lpl in the pre-collected liver biopsy were measured with quantitative real-time PCR analysis and compared with the blood chemistry data after APAP treatment. Correlation of expression level of Gabra3 with serum ALT (A) or AST (B) and that between Lpl and total bilirubin (C). *Represents statistical significance analyzed by Spearman correlation test and r notes Spearman’s rank correlation coefficient (p<0.05). Comparison of serum ALT (D), AST (E) or total bilirubin (F) levels of animals in 1st quartile of Gabra3 and 4th of Lpl expression levels (Predicted Res, N=8) with those of 4th of Gabra3 or 1stLpl (Predicted Sus, N=8). *p<0.05 and **p<0.01.

Significantly different genes in the liver pre-biopsy of Susceptible group vs. Resistant

GeneNCBI RefSeq No.Descriptionlog2 ratioAbsolute fold changep-valueGO Process
Gdap10BC052902Ganglioside-induced differentiation-associated-protein
LplNM_008509Lipoprotein lipase0.911.880.0004Lipid metabolic process // positive regulation of macrophage derived foam cell differentiation // lipid catabolic process // triglyceride biosynthetic process // triglyceride catabolic process
Klf10NM_013692Kruppel-like factor 100.821.770.0011Transcription // induction of apoptosis // bone mineralization // regulation of cell proliferation // regulation of transcription // positive regulation of osteoclast differentiation // response to protein stimulus
Malat1NR_002847Metastasis associated lung adenocarcinoma transcript 10.811.750.0040Biological_process
Gas5NR_002840Growth arrest specific0.801.750.0054-
Ankrd55NM_029898Ankyrin repeat domain 550.791.730.0066Molecular_function // cellular_component // biological_process
Cspp1NM_026493Centrosome and spindle pole associated protein 10.731.660.0024Biological_process // cellular_component // molecular_function
Dennd4aNM_001162917DENN/MADD domain containing 4A0.701.620.0048Biological_process // cellular_component // molecular_function
Chd7NM_001081417chromodomain helicase DNA binding protein 70.641.550.0014Blood vessel development // in utero embryonic development // heart morphogenesis // chromatin assembly or disassembly // transcription // adult heart development // sensory perception of sound // locomotory behavior // adult walking behavior // blood circulation // chromatin modification // female genitalia development // embryonic hindlimb morphogenesis // positive regulation of multicellular organism growth // ear morphogenesis // inner ear morphogenesis // camera-type eye development // nose development // regulation of transcription // palate development chromatin // cellular_component // nucleus nucleotide binding // molecular_function // nucleic acid binding // DNA binding // chromatin binding // helicase activity // ATP binding // hydrolase activity
Snord58bNR_028552Small nucleolar RNA, C/D box 58B0.621.540.0011Biological_process //cellular_component // molecular_function //
Dclre1cNM_146114DNA cross-link repair 1C, PSO2 homolog (S. cerevisiae)0.621.530.0008Telomere maintenance // DNA repair // double-strand break repair // DNA recombination // response to DNA damage stimulus // response to ionizing radiation // B cell differentiation // chromosome organization nucleus single-stranded DNA specific endodeoxyribonuclease activity // nuclease activity // endonuclease activity // exonuclease activity // 5′-3′ exonuclease activity // hydrolase activity
Taf1dBC056964TATA box binding protein (Tbp)-associated factor, RNA polymerase I, D0.601.520.0036Transcription // biological_process // regulation of transcription cellular_component // nucleus molecular_function // DNA binding
PdiltNM_027943Protein disulfide isomerase-like, testis expressed0.601.510.0044Multicellular organismal development // spermatogenesis // biological_process // metabolic process // cell differentiation // cell redox homeostasis // cellular_component // endoplasmic reticulum molecular_function // isomerase activity
Meig1NM_008579Meiosis expressed gene 10.601.510.0066Meiosis nucleus protein binding
Ccrn4lNM_009834CCR4 carbon catabolite repression 4-like (S. cerevisiae)−0.631.550.0008Rhythmic process
Gabra3NM_008067Gamma-aminobutyric acid (Saeed et al.) A receptor, subunit alpha 3−1.562.940.0034Transport // ion transport // chloride transport // signal transduction // signal transduction // gamma-aminobutyric acid signaling pathway // synaptic transmission // ion transmembrane transport
  1. Baron, AD, Brechtel, G, Wallace, P, and Edelman, SV (1988). Rates and tissue sites of non-insulin- and insulin-mediated glucose uptake in humans. Am J Physiol. 255, E769-E774.
  2. Bessems, JG, and Vermeulen, NP (2001). Paracetamol (acetaminophen)-induced toxicity: molecular and biochemical mechanisms, analogues and protective approaches. Crit Rev Toxicol. 31, 55-138.
    Pubmed CrossRef
  3. Biju, MP, Pyroja, S, Rajeshkumar, NV, and Paulose, CS (2001). Hepatic GABA(A) receptor functional regulation during rat liver cell proliferation. Hepatol. Res. 21, 136-146.
    Pubmed CrossRef
  4. Bozogluer, E, Madenoglu, H, Aksu, R, Bicer, C, Yazici, C, and Boyaci, A (2012). The effect of different doses of flumazenil on acetaminophen toxicity in rats. Bratisl. Lek. Listy. 113, 525-528.
  5. Chang, CY, and Schiano, TD (2007). Review article: drug hepatotoxicity. Aliment Pharmacol Ther. 25, 1135-1151.
    Pubmed CrossRef
  6. Clark, SJ, Shojaee-Moradie, F, Croos, P, Seed, PT, Umpleby, AM, Wendon, JA, and Miell, J (2001). Temporal changes in insulin sensitivity following the development of acute liver failure secondary to acetaminophen. Hepatology. 34, 109-115.
    Pubmed CrossRef
  7. Draghici, S, Khatri, P, Eklund, AC, and Szallasi, Z (2006). Reliability and reproducibility issues in DNA microarray measurements. Trends Genet. 22, 101-109.
  8. Fausto, N, Campbell, JS, and Riehle, KJ (2006). Liver regeneration. Hepatology. 43, S45-S53.
    Pubmed CrossRef
  9. Goldberg, IJ (1996). Lipoprotein lipase and lipolysis: central roles in lipoprotein metabolism and atherogenesis. J Lipid Res. 37, 693-707.
  10. Green, CB, Douris, N, Kojima, S, Strayer, CA, Fogerty, J, Lourim, D, Keller, SR, and Besharse, JC (2007). Loss of Nocturnin, a circadian deadenylase, confers resistance to hepatic steatosis and diet-induced obesity. Proc Natl Acad Sci USA. 104, 9888-9893.
    Pubmed KoreaMed CrossRef
  11. Gutschner, T, Hammerle, M, Eissmann, M, Hsu, J, Kim, Y, Hung, G, Revenko, A, Arun, G, Stentrup, M, Gross, M, Zornig, M, MacLeod, AR, Spector, DL, and Diederichs, S (2013). The Noncoding RNA MALAT1 Is a Critical Regulator of the Metastasis Phenotype of Lung Cancer Cells. Cancer Res. 73, 1180-1189.
  12. Hinson, JA, Roberts, DW, and James, LP (2010). Mechanisms of acetaminophen-induced liver necrosis. Handb Exp Pharmacol. 196, 369-405.
  13. Ji, P, Diederichs, S, Wang, WB, Boing, S, Metzger, R, Schneider, PM, Tidow, N, Brandt, B, Buerger, H, Bulk, E, Thomas, M, Berdel, WE, Serve, H, and Muller-Tidow, C (2003). MALAT-1, a novel noncoding RNA, and thymosin beta 4 predict metastasis and survival in early-stage non-small cell lung cancer. Oncogene. 22, 8031-8041.
    Pubmed CrossRef
  14. Johnsen, SA, Subramaniam, M, Janknecht, R, and Spelsberg, TC (2002). TGF beta inducible early gene enhances TGF beta/Smad-dependent transcriptional responses. Oncogene. 21, 5783-5790.
    Pubmed CrossRef
  15. Kim, JK, Fillmore, JJ, Chen, Y, Yu, C, Moore, IK, Pypaert, M, Lutz, EP, Kako, Y, Velez-Carrasco, W, Goldberg, IJ, Breslow, JL, and Shulman, GI (2001). Tissue-specific overexpression of lipoprotein lipase causes tissue-specific insulin resistance. Proc Natl Acad Sci USA. 98, 7522-7527.
    Pubmed KoreaMed CrossRef
  16. Krenzelok, EP (2009). The FDA Acetaminophen Advisory Committee Meeting - what is the future of acetaminophen in the United States? The perspective of a committee member. Clin. Toxicol. (Phila.). 47, 784-789.
  17. Larson, AM, Polson, J, Fontana, RJ, Davern, TJ, Lalani, E, Hynan, LS, Reisch, JS, Schiødt, FV, Ostapowicz, G, Shakil, AO, and Lee, WM (2005). Acetaminophen-induced acute liver failure: results of a United States multicenter, prospective study. Hepatology. 42, 1364-1372.
    Pubmed CrossRef
  18. Lazarou, J, Pomeranz, BH, and Corey, PN (1998). Incidence of adverse drug reactions in hospitalized patients: a meta-analysis of prospective studies. JAMA. 279, 1200-1205.
    Pubmed CrossRef
  19. Lee, SS, Buters, JT, Pineau, T, Fernandez-Salguero, P, and Gonzalez, FJ (1996). Role of CYP2E1 in the hepatotoxicity of acetaminophen. J Biol Chem. 271, 12063-12067.
    Pubmed CrossRef
  20. Lewis, JH (2000). Drug-induced liver disease. Med Clin North Am. 84, 1275-1311.
    Pubmed CrossRef
  21. Liu, H, Nakagawa, T, Kanematsu, T, Uchida, T, and Tsuji, S (1999). Isolation of 10 differentially expressed cDNAs in differentiated Neuro2a cells induced through controlled expression of the GD3 synthase gene. J Neurochem. 72, 1781-1790.
    Pubmed CrossRef
  22. Liu, HH, Lu, P, Guo, YY, Farrell, E, Zhang, X, Zheng, M, Bosano, B, Zhang, Z, Allard, J, Liao, G, Fu, S, Chen, J, Dolim, K, Kuroda, A, Usuka, J, Cheng, J, Tao, W, Welch, K, Liu, Y, Pease, J, de Keczer, SA, Masjedizadeh, M, Hu, JS, Weller, P, Garrow, T, and Peltz, G (2010). An integrative genomic analysis identifies Bhmt2 as a diet-dependent genetic factor protecting against acetaminophen-induced liver toxicity. Genome Res. 20, 28-35.
    KoreaMed CrossRef
  23. Lu, XY, Hu, B, Zheng, J, Ji, C, Fan, XH, and Gao, Y (2015). Predose and postdose blood gene expression profiles identify the individuals susceptible to acetaminophen-induced liver injury in rats. PLoS ONE. 10, e0141750.
    Pubmed KoreaMed CrossRef
  24. Michalopoulos, GK (2007). Liver regeneration. J Cell Physiol. 213, 286-300.
    Pubmed KoreaMed CrossRef
  25. Moe, M, Lien, S, Bendixen, C, Hedegaard, J, Hornshøj, H, Berget, I, Meuwissen, TH, and Grindflek, E (2008). Gene expression profiles in liver of pigs with extreme high and low levels of androstenone. BMC Vet Res. 4, 29.
    Pubmed KoreaMed CrossRef
  26. Oh, J-H, Yoon, H-J, Lim, J-S, Park, H-J, Cho, J-W, Kwon, M-S, and Yoon, S-J (2009). Analysis of Gene Expression in 4, 4′-Methylenedianiline-induced Acute Hepatotoxicity. Toxicol Res. 25, 85-92.
  27. Ostapowicz, G, Fontana, RJ, Schiodt, FV, Larson, A, Davern, TJ, Han, SH, McCashland, TM, Shakil, AO, Hay, JE, Hynan, L, Crippin, JS, Blei, AT, Samuel, G, Reisch, J, and Lee, WM (2002). Results of a prospective study of acute liver failure at 17 tertiary care centers in the United States. Ann Intern Med. 137, 947-954.
    Pubmed CrossRef
  28. Park, A (2006). The FDA’s painkiller warning: how to avoid taking too much [Internet]. TIME.Available from:,8599,1572012,00.html.
  29. Prescott, LF, Park, J, Ballantyne, A, Adriaenssens, P, and Proudfoot, AT (1977). Treatment of paracetamol (acetaminophen) poisoning with N-acetylcysteine. Lancet. 2, 432-434.
    Pubmed CrossRef
  30. Reuben, A, Koch, DG, and Lee, WM (2010). Drug-induced acute liver failure: results of a U.S. multicenter, prospective study. Hepatology. 52, 2065-2076.
    Pubmed KoreaMed CrossRef
  31. Saha, B, and Nandi, D (2009). Farnesyltransferase inhibitors reduce ras activation and ameliorate acetaminophen-induced liver injury in mice. Hepatology. 50, 1547-1557.
    Pubmed CrossRef
  32. Sanaka, M, Kuyama, Y, and Yamanaka, M (1998). Guide for judicious use of the paracetamol absorption technique in a study of gastric emptying rate of liquids. J Gastroenterol. 33, 785-791.
    Pubmed CrossRef
  33. Shin, JH, Lee, CW, Oh, SJ, Yun, J, Kang, MR, Han, SB, Park, H, Jung, JC, Chung, YH, and Kang, JS (2014). Hepato-protective Effect of Aged Black Garlic Extract in Rodents. Toxicol Res. 30, 49-54.
    Pubmed KoreaMed CrossRef
  34. Stamper, BD, Bammler, TK, Beyer, RP, Farin, FM, and Nelson, SD (2010). Differential regulation of mitogen-activated protein kinase pathways by acetaminophen and its nonhepatotoxic regioisomer 3′-hydroxyacetanilide in TAMH cells. Toxicol Sci. 116, 164-173.
    Pubmed KoreaMed CrossRef
  35. Stephens, C, Lucena, MI, and Andrade, RJ (2012). Genetic variations in drug-induced liver injury (DILI): resolving the puzzle. Front Genet. 3, 253.
    Pubmed KoreaMed CrossRef
  36. Umbright, C, Sellamuthu, R, Li, SQ, Kashon, M, Luster, M, and Joseph, P (2010). Blood gene expression markers to detect and distinguish target organ toxicity. Mol Cell Biochem. 335, 223-234.
  37. Wang, H, and Eckel, RH (2009). Lipoprotein lipase: from gene to obesity. Am J Physiol Endocrinol Metab. 297, E271-E288.
    Pubmed CrossRef
  38. Watkins, PB, Kaplowitz, N, Slattery, JT, Colonese, CR, Colucci, SV, Stewart, PW, and Harris, SC (2006). Aminotransferase elevations in healthy adults receiving 4 grams of acetaminophen daily: a randomized controlled trial. JAMA. 296, 87-93.
    Pubmed CrossRef
  39. Watson, WA, Litovitz, TL, Klein-Schwartz, W, Rodgers, GC, Youniss, J, Reid, N, Rouse, WG, Rembert, RS, and Borys, D (2004). 2003 Annual Report of the American Association of Poison Control Centers Toxic Exposure Surveillance System. Am J Emerg Med. 22, 335-404.
    Pubmed CrossRef
  40. Yun, JW, Kim, CW, Bae, IH, Park, YH, Chung, JH, Lim, KM, and Kang, KS (2009). Determination of the key innate genes related to individual variation in carbon tetrachloride-induced hepatotoxicity using a pre-biopsy procedure. Toxicol Appl Pharmacol. 239, 55-63.
    Pubmed CrossRef
  41. Yun, JW, Kim, CW, Bae, IH, Park, YH, Chung, JH, Lim, KM, and Kang, KS (2010). Expression levels of pituitary tumor transforming 1 and glutathione-S-transferase theta 3 are associated with the individual susceptibility to D-galactosamine-induced hepatotoxicity. Toxicol Appl Pharmacol. 242, 91-99.
  42. Yun, JW, Kim, M, Cho, SD, Lee, JY, Bae, ON, and Lim, KM (2014). Highly expressed protein kinase A inhibitor alpha and suppression of protein kinase A may potentiate acetaminophen-induced hepatotoxicity. Toxicol Lett. 229, 59-65.
    Pubmed CrossRef

This Article

Cited By Articles
  • CrossRef (0)

Funding Information

Social Network Service