¡¡

Comprehensive analysis of differential gene expression profiles on

 D-galactosamine-induced acute mouse liver injury and regeneration

 

Heekyoung Chung,a,b Hyun-Jun Kim, a Ki-Seok Jang, a Mingoo Kim,c Jungeun Yang,d Kyung-Sun Kang, b,e

 Hyung-Lae Kim, b,f Byung-Il Yoon, b,g Mi-Ock Lee, b,h Byung-Hoon Lee, b,i Ju Han Kim,b,c Yong-Sung Lee,d and Gu Kong a,b*

Affiliation and Addresses:
 a Department of Pathology, and d Department of Biochemistry, College of Medicine, Hanyang University, Seoul 133-791, Republic of Korea;
 b Toxicogenomics Research Consortium (TGRC), Hanyang University, Seoul 133-791, Republic of Korea;
 c Seoul National University Biomedical Informatics (SNUBI), Seoul National University College of Medicine, Seoul 110-799, Republic of Korea;
 e
Department of Veterinary Public Health, College of Veterinary Medicine, Seoul National University, Seoul 151-742, Republic of Korea;
 f Department of Biochemistry, College of Medicine, Ewha Womans University, Seoul 158-710, Republic of Korea;
 g School of Veterinary Medicine, Kangwon National University, Chuncheon 200-701, Republic of Korea;
 h College of Pharmacy and Bio-MAX Institute, Seoul National University, Seoul 151-742, Republic of Korea; and
 i College of Pharmacy and Research Institute of Pharmaceutical Sciences, Seoul National University, Seoul 151-742, Republic of Korea.

 *Corresponding Author. Department of Pathology, College of Medicine, Hanyang University, 17 Haengdang-dong, Seongdong-gu, Seoul, 133-791, Republic of Korea. Tel: +82-2-2290-8251, Fax: +82-2-2295-1091,
 E-mail address: gkong@hanyang.ac.kr (G. Kong).
 Key Words: D-galactosamine (GalN); mouse; liver; toxicogenomics; Applied Biosystems Mouse Genome Survey Microarray

Abstract

Microarray analysis of RNA from D-galactosamine (GalN)-administered mouse livers was performed to establish a global gene expression profile during injury and regeneration stages at two different doses. A single dose of GalN at 266 or 26.6 mg/kg body weight was given intraperitoneally, and the liver samples were obtained after 6, 24, and 72 h. Histopathologic studies enabled the classification of the D-galactosamine effect into injury (6, 24 h) and regeneration (72 h) stages. By using the Applied Biosystems Mouse Genome Survey Microarray, a total of 7,267 out of 33,315 (21.8%) genes were found to be statistically reliable at p < 0.05 by 2-way ANOVA, and 1,469 (4.4%) probes at false discovery rate < 5% by Significance Analysis of Microarray. Among the statistically reliable clones by both analytical methods, 389 genes were differentially expressed when compared with non-treated control, with more than a 1.625-fold difference (which equals 0.7 in log2 scale) at one or more GalN treatment conditions and with less than 1.625-fold difference at all three vehicle-treated conditions. Three hundred thirty six genes and 13 genes were identified as injury- and regeneration-specific genes, respectively, showing that most of the transcriptomic changes were seen during the injury stage. Furthermore, multiple genes involved in protein synthesis and degradation, mRNA processing and binding, and cell cycle regulation showed variable transcript levels upon acute GalN administration.

¡¡

¡¡

¡¡

1. Experimental Design  

 

Array: Applied Biosystems AB1700 mouse chip (33315 probes) 

30 arrays in total with triple biological replicates

¡¡

2. Data

2.1. Raw data

Name Files
galn_High_72hr_3  
galn_High_72hr_2  
galn_High_72hr_1  
galn_High_24hr_3  
galn_High_24hr_2  
galn_High_24hr_1  
galn_High_06hr_3  
galn_High_06hr_2  
galn_High_06hr_1  
galn_Low_72hr_3  
galn_Low_72hr_2  
galn_Low_72hr_1  
galn_Low_24hr_3  
galn_Low_24hr_2  
galn_Low_24hr_1  
galn_Low_06hr_3  
galn_Low_06hr_2  
galn_Low_06hr_1  
galn_Vehicle_72hr_3  
galn_Vehicle_72hr_2  
galn_Vehicle_72hr_1  
galn_Vehicle_24hr_3  
galn_Nontreat_NA_1  
galn_Nontreat_NA_2  
galn_Nontreat_NA_3  
galn_Vehicle_06hr_1  
galn_Vehicle_06hr_2  
galn_Vehicle_06hr_3  
galn_Vehicle_24hr_1  
galn_Vehicle_24hr_2  

¡¡

2.2. Normalized data

    Download normalized data (33315 probes x 30 samples)

    Data are normalized in following procedure.

  1. Assay Normalized Signals are used.
  2. Only probes with Flag < 100 are used.
  3. Missing and filtered valued are imputed with KNN imputation algorithm.
  4. Values are VSN transformed.
  5. Arrays are scaled with quantile normalization.
  6. Normalized valued are rescaled with

¡¡

3. Diagnostic Plots

3.1. Single-chip diagnostic plots

Name Analysis
gain_Nontreat_NA_1
gain_Nontreat_NA_2
gain_Nontreat_NA_3
gain_Vehicle_06hr_1
gain_Vehicle_06hr_2
gain_Vehicle_06hr_3
gain_Vehicle_24hr_1
gain_Vehicle_24hr_2
gain_Vehicle_24hr_3
gain_Vehicle_72hr_1
gain_Vehicle_72hr_2
gain_Vehicle_72hr_3
gain_Low_06hr_1
gain_Low_06hr_2
gain_Low_06hr_3
gain_Low_24hr_1
gain_Low_24hr_2
gain_Low_24hr_3
gain_Low_72hr_1
gain_Low_72hr_2
gain_Low_72hr_3
gain_High_06hr_1
gain_High_06hr_2
gain_High_06hr_3
gain_High_24hr_1
gain_High_24hr_2
gain_High_24hr_3
gain_High_72hr_1
gain_High_72hr_2
gain_High_72hr_3

¡¡

3.2. Between-chip diagnostic plots

[View between-chip diagnostic plots]

¡¡

4. Differentially Expressed Genes (DEGs)

        4.1.  Lists of DEGs

Supplementary Table 1. Microarray results of 7,267 probes found to be statistically reliable at p < 0.05 by 2-way ANOVA
Supplementary Table 2. Microarray results of 4,240 "differentially regulated" probes found to be statistically reliable at p < 0.05 by 2-way ANOVA
Supplementary Table 3. Microarray results of 1,469 probes found to be statistically reliable by SAM at FDR < 5%
Supplementary Table 4. Microarray results of 460 "differentially regulated" probes found to be statistically reliable by SAM at FDR < 5%
Supplementary Table 5. Microarray results of 389 "differentially regulated" probes identified by SAM and 2-way ANOVA
Supplementary Table 6. Microarray results of 138 injury-specific downregulated genes identified by 2-way ANOVA and SAM
Supplementary Table 7. Microarray results of 198 injury-specific upregulated genes identified by 2-way ANOVA and SAM
Supplementary Table 8. Microarray results of 13 regeneration-specific genes identified by 2-way ANOVA and SAM
Supplementary Table 9. GenMAPP result of injuy- or regeneration-specific genes

DOWNLOAD AS EXCEL FILE

¡¡

4.2. SOM clustering and Gene ontology and pathway functional enrichment study

                    SOM clustering was performed on each ANOVA group. Then Enrichment study using Gene ontology and pathway information was performed on each cluster.
                    Broken line graphs show general tendency of expression profile for each cluster. (Installation of SVG viewer under Internet Explore is requiered.)   

P-value < 0.05
MTC MABI MABx MAxI MAxx MxBI MxBx MxxI Mxxx
un 314 C A 679 C A 201 C A 645 C A 977 C A 3697 C A 637 C A 117 C A
bf 2 C A 0 C A 0 C A 0 C A 2 C A 264 C A 0 C A 0 C A
bh 14 C A 0 C A 0 C A 0 C A 33 C A 2276 C A 0 C A 0 C A
by 1 C A 0 C A 0 C A 0 C A 2 C A 872 C A 0 C A 0 C A

D = Dose
T = Time
Click C for Cluster viewer, GO enrichement analysis
Click A for ArrayXPath, Pathway enrichment analysis
UN = Unadjusted
BF = Bonferoni multiple test correction
BH = Benjamini :: Hochberg multiple test correction
BY = Benjamini & Yekutieli multiple test correction

¡¡

¡¡

Last updated: 2006-06-14