In functional genomics, DNA microarrays for gene expression profiling are increasingly being used to provide insights into biological function or pathology. To better understand the significance of the multiple transcriptional changes across a time period, the temporal changes in phenotype must be described. Orotic acid-induced fatty liver disease was investigated at the transcriptional and metabolic levels using microarrays and metabolic profiling in two strains of rats. High-resolution 1H-NMR spectroscopic analysis of liver tissue indicated that Kyoto rats compared with Wistar rats are predisposed to the insult. Metabolite analysis and gene expression profiling following orotic acid treatment identified perturbed metabolic pathways, including those involved in fatty acid, triglyceride, and phospholipid synthesis, beta-oxidation, altered nucleotide, methyl donor, and carbohydrate metabolism, and stress responses. Multivariate analysis and statistical bootstrapping were used to investigate co-responses with transcripts involved in metabolism and stress responses. This reverse functional genomic strategy highlighted the relationship between changes in the transcription of stearoyl-CoA desaturase 1 and those of other lipid-related transcripts with changes in NMR-derived lipid profiles. The results suggest that the integration of 1H-NMR and gene expression data sets represents a robust method for identifying a focused line of research in a complex system.