Targeted disruption of glycerol kinase gene in mice: expression analysis in liver. Mus musculus

Glycerol kinase deficiency (GKD) is an X-linked inborn error of metabolism with metabolic and neurologic crises. Liver shows the highest level of glycerol kinase (GK) activity in humans and mice. Absence of genotype-phenotype correlations in patients with GKD indicate 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 the overall gene expression profile of the KO mice was different from that of WT. Real time PCR confirmed differences for selected genes. Functional gene enrichment analysis was used to find 56 increased and 37 decreased gene functional categories. Pathway Assist 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 SREBP-1c, ChREBP, HNF-4alpha, and PPAR-alpha, had increased activity in the Gyk KO mice compared with WT mice; while SREBP-2 was less active in the Gyk KO mice. These studies show that Gyk deletion causes alterations in gene expression of genes in several regulatory networks and is the first time NCA has been used to expand on microarray data from a mouse knockout model of a human disease. Keywords: Glycerol kinase (Gyk) knockout mouse; mouse model of human Glycerol Kinase Deficiency; Gyk KO versus WT liver expression analysis; Affymetrix mus 430 2.0 GeneChip Overall design: Male WT and KO mouse pups were sacrificed on day of life (dol) 3 and each liver was harvested. Total RNA from 4 KO and 4 WT livers was isolated individually. cDNA was synthesized from the poly(A)+ mRNA in the total RNA, Biotin-tagged and fragmented to an average strand length of 100 bases (range 35-200 bases). Ten µg of each cRNA was hybridized onto an Affymetrix mus 430 2.0 GeneChip to analyze differences in liver gene expression between KO and WT mice. Day of life three was chosen because the mice are phenotypically symptomatic with statistically different parameters for hypoglycemia, acidosis; low bicarbonate and decreased base excess. On day of life 2 they are not significantly different from wild type in all of these important clinical phenotypes.

甘油激酶缺乏症（Glycerol kinase deficiency, GKD）是一种X连锁遗传性代谢病，可引发代谢与神经系统危象。在人类与小鼠体内，肝脏的甘油激酶（glycerol kinase, GK）活性最高。GKD患者中基因型与表型无相关性的现象，提示存在修饰基因（包括其他网络互作因子）的参与。为解析GKD的分子发病机制，本研究对新生期甘油激酶（Gyk）敲除（knockout, KO）与野生型（wild type, WT）小鼠的肝脏mRNA进行了基因微阵列分析。无监督学习分析显示，KO小鼠的整体基因表达谱与WT小鼠存在显著差异。实时定量PCR（real time PCR）验证了部分筛选基因的表达差异。基因功能富集分析共筛选得到56个表达上调、37个表达下调的基因功能类别。通路辅助分析（Pathway Assist analysis）发现，参与有机酸代谢的基因表达水平发生改变，提示GK隶属于同一代谢网络，这与GKD患者在发作危象期间出现代谢性酸血症的临床表型高度吻合。网络组分分析（Network component analysis, NCA）结果显示，与WT小鼠相比，Gyk KO小鼠体内的转录因子SREBP-1c、ChREBP、HNF-4α及PPAR-α活性上调，而SREBP-2活性下调。本研究证实，Gyk敲除会导致多个调控网络中的基因表达发生改变，且本研究首次将NCA应用于人类疾病的小鼠敲除模型的微阵列数据拓展分析。
关键词：甘油激酶（Gyk）敲除小鼠；人类甘油激酶缺乏症小鼠模型；Gyk KO与WT肝脏表达谱分析；Affymetrix小鼠基因组430 2.0基因芯片（Affymetrix mus 430 2.0 GeneChip）
实验设计：选取雄性WT与KO小鼠幼崽，于出生后第3天（dol 3）处死并采集肝脏组织。分别提取4只KO小鼠与4只WT小鼠肝脏的总RNA，以总RNA中的poly(A)+ mRNA反转录合成cDNA，经生物素标记与片段化处理后，平均链长为100碱基（范围35~200碱基）。取10 μg上述cRNA分别与Affymetrix小鼠基因组430 2.0基因芯片（Affymetrix mus 430 2.0 GeneChip）杂交，以分析KO与WT小鼠肝脏的基因表达差异。选择出生后第3天作为实验时间点，是因为该时期小鼠已出现表型症状，其低血糖、酸中毒、碳酸氢盐水平降低及碱剩余减少等参数均存在统计学差异；而出生后第2天，小鼠在上述关键临床表型上与WT小鼠无显著差异。