A proteomics approach combined with multivariate data analysis was used to examine the hepatotoxic effect of hydrazine in 30 male Sprague Dawley rats, assigned to four treatment groups and two control groups. Liver samples from the individual animals were resolved by two-dimensional differential gel electrophoresis (2-D DIGE) and protein patterns from the 2-D gels were analyzed by principal component analysis (PCA) and partial least squares regression (PLSR). The PCA plot was able to describe the variation in the protein expression related to dose and time, by separation or clustering of different animal groups. PLSR followed by variable selection (Jack-knifing) was used to select proteins that varied significantly in relation to the dose related response of the hydrazine treatment. The 10 up-regulated and 10 down-regulated proteins with highest rank in the PLSR model were identified by mass spectrometry. Hydrazine treatment induced altered expression of proteins related to lipid metabolism, Ca(2+) homeostasis, thyroid hormone pathways and stress response. Several of the identified proteins have not previously been implicated in hydrazine toxicity and may thus be regarded as new potential biomarkers of induced liver toxicity.