Targeted disruption of glycerol kinase gene in mice: expression analysis in liver shows alterations in network partners related to glycerol kinase activity.

Glycerol kinase deficiency (GKD) is an X-linked inborn error of metabolism with metabolic and neurological crises. Liver shows the highest level of glycerol kinase (GK) activity in humans and mice. Absence of genotype-phenotype correlations in patients with GKD indicates the involvement of modifier genes, including other network partners. To understand the molecular pathogenesis of GKD, we performed microarray analysis on liver mRNA from neonatal glycerol kinase (Gyk) knockout (KO) and wild-type (WT) mice. Unsupervised learning revealed that the overall gene expression profile of the KO mice was different from that of WT. Real-time PCR confirmed the differences for selected genes. Functional gene enrichment analysis was used to find 56 increased and 37 decreased gene functional categories. PathwayAssist analysis identified changes in gene expression levels of genes involved in organic acid metabolism indicating that GK was part of the same metabolic network which correlates well with the patients with GKD having metabolic acidemia during their episodic crises. Network component analysis (NCA) showed that transcription factors sterol regulatory element-binding protein (SREBP)-1c, carbohydrate response element-binding protein (ChREBP), hepatocyte nuclear factor-4 alpha (HNF-4alpha) and peroxisome proliferative-activated receptor-alpha (PPARalpha) had increased activity in the Gyk KO mice compared with WT mice, whereas SREBP-2 was less active in the Gyk KO mice. These studies show that Gyk deletion causes alterations in expression of genes in several regulatory networks and is the first time NCA has been used to expand on microarray data from a mouse KO model of a human disease.