Ultrafine mapping of SNPs from mouse strains C57BL/6J, DBA/2J, and C57BLKS/J for loci contributing to diabetes and atherosclerosis susceptibility

Increased glycolysis affects β-cell function and identity in aging and diabetes

Objective

Age is a risk factor for type 2 diabetes (T2D). We aimed to elucidate whether β-cell glucose metabolism is altered with aging and contributes to T2D.

Methods

We used senescence-accelerated mice (SAM), C57BL/6J (B6) mice, and ob/ob mice as aging models. As a diabetes model, we used db/db mice. The glucose responsiveness of insulin secretion and the [U-¹³C]-glucose metabolic flux were examined in isolated islets. We analyzed the expression of β-cell-specific genes in isolated islets and pancreatic sections as molecular signatures of β-cell identity. β cells defective in the malate-aspartate (MA) shuttle were previously generated from MIN6-K8 cells by the knockout of Got1, a component of the shuttle. We analyzed Got1 KO β cells as a model of increased glycolysis.

Results

We identified hyperresponsiveness to glucose and compromised cellular identity as dysfunctional phenotypes shared in common between aged and diabetic mouse β cells. We also observed a metabolic commonality between aged and diabetic β cells: hyperactive glycolysis through the increased expression of nicotinamide mononucleotide adenylyl transferase 2 (Nmnat2), a cytosolic nicotinamide adenine dinucleotide (NAD)-synthesizing enzyme. Got1 KO β cells showed increased glycolysis, β-cell dysfunction, and impaired cellular identity, phenocopying aging and diabetes. Using Got1 KO β cells, we show that attenuation of glycolysis or Nmnat2 activity can restore β-cell function and identity.

Conclusions

Our study demonstrates that hyperactive glycolysis is a metabolic signature of aged and diabetic β cells, which may underlie age-related β-cell dysfunction and loss of cellular identity. We suggest Nmnat2 suppression as an approach to counteract age-related T2D.

•

Glucose hypersensitivity and impaired identity are common features of aged and diabetic β cells.

•

Metabolic tracing reveals increased glycolysis and altered NAD production in aged and diabetic β cells.

•

Increased glycolysis induces β-cell dysfunction and loss of identity.

•

NAD production by Nmnat2 can be targeted to restore β-cell phenotypes.

Leptin Receptor-Deficient db/db Mice Show Significant Heterogeneity in Response to High Non-heme Iron Diet

Recent studies have shown an association between iron homeostasis, obesity and diabetes. In this work, we investigated the differences in the metabolic status and inflammation in liver, pancreas and visceral adipose tissue of leptin receptor-deficient db/db mice dependent on high iron concentration diet. 3-month-old male BKS-Leprdb/db/JOrlRj (db/db) mice were divided into two groups, which were fed with different diets containing high iron (29 g/kg, n = 57) or standard iron (0.178 g/kg; n = 42) concentrations for 4 months. As anticipated, standard iron-fed db/db mice developed obesity and diabetes. However, high iron-fed mice exhibited a wide heterogeneity. By dividing into two subgroups at the diabetes level, non-diabetic subgroup 1 (<13.5 mmol/l, n = 30) significantly differed from diabetic subgroup two (>13.5 mmol/l, n = 27). Blood glucose concentration, HbA1c value, inflammation markers interleukin six and tumor necrosis factor α and heme oxygenase one in visceral adipose tissue were reduced in subgroup one compared to subgroup two. In contrast, body weight, C-peptide, serum insulin and serum iron concentrations, pancreatic islet and signal ratio as well as cholesterol, LDL and HDL levels were enhanced in subgroup one. While these significant differences require further studies and explanation, our results might also explain the often-contradictory results of the metabolic studies with db/db mice.

DBA2J db/db mice are susceptible to early albuminuria and glomerulosclerosis that correlate with systemic insulin resistance

Diabetic nephropathy (DN) is the leading cause of kidney failure in the world. To understand important mechanisms underlying this condition, and to develop new therapies, good animal models are required. In mouse models of type 1 diabetes, the DBA/2J strain has been shown to be more susceptible to develop kidney disease than other common strains. We hypothesized this would also be the case in type 2 diabetes. We studied db/db and wild-type (wt) DBA/2J mice and compared these with the db/db BLKS/J mouse, which is currently the most widely used type 2 DN model. Mice were analyzed from age 6 to 12 wk for systemic insulin resistance, albuminuria, and glomerular histopathological and ultrastructural changes. Body weight and nonfasted blood glucose were increased by 8 wk in both genders, while systemic insulin resistance commenced by 6 wk in female and 8 wk in male db/db DBA/2J mice. The urinary albumin-to-creatinine ratio (ACR) was closely linked to systemic insulin resistance in both sexes and was increased ~50-fold by 12 wk of age in the db/db DBA/2J cohort. Glomerulosclerosis, foot process effacement, and glomerular basement membrane thickening were observed at 12 wk of age in db/db DBA/2J mice. Compared with db/db BLKS/J mice, db/db DBA/2J mice had significantly increased levels of urinary ACR, but similar glomerular histopathological and ultrastructural changes. The db/db DBA/2J mouse is a robust model of early-stage albuminuric DN, and its levels of albuminuria correlate closely with systemic insulin resistance. This mouse model will be helpful in defining early mechanisms of DN and ultimately the development of novel therapies.

Rapamycin Increases Mortality in db/db Mice, a Mouse Model of Type 2 Diabetes

We examined the effect of rapamycin on the life span of a mouse model of type 2 diabetes, db/db mice. At 4 months of age, male and female C57BLKSJ-lepr (db/db) mice (db/db) were placed on either a control diet, lacking rapamycin or a diet containing rapamycin and maintained on these diets over their life span. Rapamycin was found to reduce the life span of the db/db mice. The median survival of male db/db mice fed the control and rapamycin diets was 349 and 302 days, respectively, and the median survival of female db/db mice fed the control and rapamycin diets was 487 and 411 days, respectively. Adjusting for gender differences, rapamycin increased the mortality risk 1.7-fold in both male and female db/db mice. End-of-life pathological data showed that suppurative inflammation was the main cause of death in the db/db mice, which is enhanced slightly by rapamycin treatment.

Divergent compensatory responses to high-fat diet between C57BL6/J and C57BLKS/J inbred mouse strains

Impaired glucose tolerance (IGT) and type 2 diabetes (T2DM) are polygenic disorders with complex pathophysiologies; recapitulating them with mouse models is challenging. Despite 70% genetic homology, C57BL/6J (BL6) and C57BLKS/J (BLKS) inbred mouse strains differ in response to diet- and genetic-induced obesity. We hypothesized these differences would yield insight into IGT and T2DM susceptibility and response to pharmacological therapies. To this end, male 8-wk-old BL6 and BLKS mice were fed normal chow (18% kcal from fat), high-fat diet (HFD; 42% kcal from fat), or HFD supplemented with the PPARγ agonist pioglitazone (PIO; 140 mg PIO/kg diet) for 16 wk. Assessments of body composition, glucose homeostasis, insulin production, and energy metabolism, as well as histological analyses of pancreata were undertaken. BL6 mice gained weight and adiposity in response to HFD, leading to peripheral insulin resistance that was met with increased β-cell proliferation and insulin production. By contrast, BLKS mice responded to HFD by restricting food intake and increasing activity. These behavioral responses limited weight gain and protected against HFD-induced glucose intolerance, which in this strain was primarily due to β-cell dysfunction. PIO treatment did not affect HFD-induced weight gain in BL6 mice, and decreased visceral fat mass, whereas in BLKS mice PIO increased total fat mass without improving visceral fat mass. Differences in these responses to HFD and effects of PIO reflect divergent human responses to a Western lifestyle and underscore the careful consideration needed when choosing mouse models of diet-induced obesity and diabetes treatment.

QTL analysis of dietary obesity in C57BL/6byj X 129P3/J F2 mice: diet- and sex-dependent effects

Obesity is a heritable trait caused by complex interactions between genes and environment, including diet. Gene-by-diet interactions are difficult to study in humans because the human diet is hard to control. Here, we used mice to study dietary obesity genes, by four methods. First, we bred 213 F₂ mice from strains that are susceptible [C57BL/6ByJ (B6)] or resistant [129P3/J (129)] to dietary obesity. Percent body fat was assessed after mice ate low-energy diet and again after the same mice ate high-energy diet for 8 weeks. Linkage analyses identified QTLs associated with dietary obesity. Three methods were used to filter candidate genes within the QTL regions: (a) association mapping was conducted using >40 strains; (b) differential gene expression and (c) comparison of genomic DNA sequence, using two strains closely related to the progenitor strains from Experiment 1. The QTL effects depended on whether the mice were male or female or which diet they were recently fed. After feeding a low-energy diet, percent body fat was linked to chr 7 (LOD = 3.42). After feeding a high-energy diet, percent body fat was linked to chr 9 (Obq5; LOD = 3.88), chr 12 (Obq34; LOD = 3.88), and chr 17 (LOD = 4.56). The Chr 7 and 12 QTLs were sex dependent and all QTL were diet-dependent. The combination of filtering methods highlighted seven candidate genes within the QTL locus boundaries: Crx, Dmpk, Ahr, Mrpl28, Glo1, Tubb5, and Mut. However, these filtering methods have limitations so gene identification will require alternative strategies, such as the construction of congenics with very small donor regions.

Genetic basis of atherosclerosis: insights from mice and humans

Atherosclerosis is a complex and heritable disease involving multiple cell types and the interactions of many different molecular pathways. The genetic and molecular mechanisms of atherosclerosis have, in part, been elucidated by mouse models; at least 100 different genes have been shown to influence atherosclerosis in mice. Importantly, unbiased genome-wide association studies have recently identified a number of novel loci robustly associated with atherosclerotic coronary artery disease. Here, we review the genetic data elucidated from mouse models of atherosclerosis, as well as significant associations for human coronary artery disease. Furthermore, we discuss in greater detail some of these novel human coronary artery disease loci. The combination of mouse and human genetics has the potential to identify and validate novel genes that influence atherosclerosis, some of which may be candidates for new therapeutic approaches.

Systems genetics of susceptibility to obesity-induced diabetes in mice

Inbred strains of mice are strikingly different in susceptibility to obesity-driven diabetes. For instance, deficiency in leptin receptor (db/db) leads to hyperphagia and obesity in both C57BL/6 and DBA/2 mice, but only on the DBA/2 background do the mice develop beta-cell loss leading to severe diabetes, while C57BL/6 mice are relatively resistant. To further investigate the genetic factors predisposing to diabetes, we have studied leptin receptor-deficient offspring of an F2 cross between C57BL/6J (db/+) males and DBA/2J females. The results show that the genetics of diabetes susceptibility are enormously complex and a number of quantitative trait loci (QTL) contributing to diabetes-related traits were identified, notably on chromosomes 4, 6, 7, 9, 10, 11, 12, and 19. The Chr. 4 locus is likely due to a disruption of the Zfp69 gene in C57BL/6J mice. To identify candidate genes and to model coexpression networks, we performed global expression array analysis in livers of the F2 mice. Expression QTL (eQTL) were identified and used to prioritize candidate genes at clinical trait QTL. In several cases, clusters of eQTLs colocalized with clinical trait QTLs, suggesting a common genetic basis. We constructed coexpression networks for both 5 and 12 wk old mice and identified several modules significantly associated with clinical traits. One module in 12 wk old mice was associated with several measures of hepatic fat content as well as with other lipid- and diabetes-related traits. These results add to the understanding of the complex genetic interactions contributing to obesity-induced diabetes.

Quantitative trait loci for electrical seizure threshold mapped in C57BLKS/J and C57BL/10SnJ mice

We mapped the quantitative trait loci (QTL) that contribute to the robust difference in maximal electroshock seizure threshold (MEST) between C57BLKS/J (BKS) and C57BL10S/J (B10S) mice. BKS, B10S, BKS × B10S F1 and BKS × B10S F2 intercross mice were tested for MEST at 8-9 weeks of age. Results of F2 testing showed that, in this cross, MEST is a continuously distributed trait determined by polygenic inheritance. Mice from the extremes of the trait distribution were genotyped using microarray technology. MEST correlated significantly with body weight and sex; however, because of the high correlation between these factors, the QTL mapping was conditioned on sex alone. A sequential series of statistical analyses was used to map QTLs including single-point, multipoint and multilocus methods. Two QTLs reached genome-wide levels of significance based upon an empirically determined permutation threshold: chromosome 6 (LOD = 6.0 at ∼69 cM) and chromosome 8 (LOD = 5.7 at ∼27 cM). Two additional QTLs were retained in a multilocus regression model: chromosome 3 (LOD = 2.1 at ∼68 cM) and chromosome 5 (LOD = 2.7 at ∼73 cM). Together the four QTLs explain one third of the total phenotypic variance in the mapping population. Lack of overlap between the major MEST QTLs mapped here in BKS and B10S mice and those mapped previously in C57BL/6J and DBA/2J mice (strains that are closely related to BKS and B10S) suggest that BKS and B10S represent a new polygenic mouse model for investigating susceptibility to seizures.

Strain- and tissue-dependent induction of plasminogen activator inhibitor-1 gene expression in fasted mice

Plasminogen activator inhibitor-1 (PAI-1), the primary physiological inhibitor of plasminogen activators, is an important contributor to hypofibrinolysis in the presence of metabolic disorders such as diabetes and obesity. The C57BLKS/J (BKS) inbred mouse strain is a popular animal model of type 2 diabetes. We previously described that food deprivation (FD) induces adipose PAI-1 expression in both lean BKS mice and BKS-db/db mice carrying a mutation in the leptin receptor gene. To evaluate the effects of the background of mouse strains, we examined FD-induced PAI-1 expression in the liver, heart and epididymal adipose tissues of BKS, C57BL/6J (B6), C3H/HeN and ICR mice. We found that PAI-1 expression is significantly induced in the heart and liver of fasted mice, although levels of expression in adipose tissues are strain-dependent. The effect of FD on plasma PAI-1 levels is also strain-dependent. Genetic background seems to be an important factor that should be considered when investigating thrombosis and fibrinolysis relative to metabolic changes in mice.

Early hepatic insulin resistance precedes the onset of diabetes in obese C57BLKS-db/db mice

OBJECTIVE

To identify metabolic derangements contributing to diabetes susceptibility in the leptin receptor-deficient obese C57BLKS/J-db/db (BKS-db) mouse strain.

RESEARCH DESIGN AND METHODS

Young BKS-db mice were used to identify metabolic pathways contributing to the development of diabetes. Using the diabetes-resistant B6-db strain as a comparison, in vivo and in vitro approaches were applied to identify metabolic and molecular differences between the two strains.

RESULTS

Despite higher plasma insulin levels, BKS-db mice exhibit lower lipogenic gene expression, rate of lipogenesis, hepatic triglyceride and glycogen content, and impaired insulin suppression of gluconeogenic genes. Hepatic insulin receptor substrate (IRS)-1 and IRS-2 expression and insulin-stimulated Akt-phosphorylation are decreased in BKS-db primary hepatocytes. Hyperinsulinemic-euglycemic clamp studies indicate that in contrast to hepatic insulin resistance, skeletal muscle is more insulin sensitive in BKS-db than in B6-db mice. We also demonstrate that elevated plasma triglyceride levels in BKS-db mice are associated with reduced triglyceride clearance due to lower lipase activities.

CONCLUSIONS

Our study demonstrates the presence of metabolic derangements in BKS-db before the onset of beta-cell failure and identifies early hepatic insulin resistance as a component of the BKS-db phenotype. We propose that defects in hepatic insulin signaling contribute to the development of diabetes in the BKS-db mouse strain.

PP2A regulatory subunit PP2A-B' counteracts S6K phosphorylation

The insulin/TOR signaling pathway plays a crucial role in animal homeostasis, sensing nutrient status to regulate organismal growth and metabolism. We identify here the Drosophila B' regulatory subunit of PP2A (PP2A-B') as a novel, conserved component of the insulin pathway that specifically targets the PP2A holoenzyme to dephosphorylate S6K. PP2A-B' knockout flies have elevated S6K phosphorylation and exhibit phenotypes typical of elevated insulin signaling such as reduced total body triglycerides and reduced longevity. We show that PP2A-B' interacts with S6K both physically and genetically. The human homolog of PP2A-B', PPP2R5C, also counteracts S6K1 phosphorylation, indicating a conserved mechanism in mammals. Since S6K affects development of cancer and metabolic disease, our data identify PPP2R5C as a novel factor of potential medical relevance.

Identification of quantitative trait loci underlying proteome variation in human lymphoblastoid cells

Population-based variability in protein expression patterns, especially in humans, is often observed but poorly understood. Moreover, very little is known about how interindividual genetic variation contributes to protein expression patterns. To begin to address this, we describe elements of technical and biological variations contributing to expression of 544 proteins in a population of 24 individual human lymphoblastoid cell lines that have been extensively genotyped as part of the International HapMap Project. We determined that expression levels of 10% of the proteins were tightly correlated to cell doubling rates. Using the publicly available genotypes for these lymphoblastoid cell lines, we applied a genetic association approach to identify quantitative trait loci associated with protein expression variation. Results identified 24 protein forms corresponding to 15 proteins for which genetic elements were responsible for >50% of the expression variation. The genetic variation associated with protein expression levels were located in cis with the gene coding for the transcript of the protein for 19 of these protein forms. Four of the genetic elements identified were coding non-synonymous single nucleotide polymorphisms that resulted in migration pattern changes in the two-dimensional gel. This is the first description of large scale proteomics analysis demonstrating the direct relationship between genome and proteome variations in human cells.

Pancreatic islet expression profiling in diabetes-prone C57BLKS/J mice reveals transcriptional differences contributed by DBA loci, including Plagl1 and Nnt

Background

C57BLKS/J (BLKS) mice are susceptible to islet exhaustion in insulin-resistant states as compared with C57BL6/J (B6) mice, as observed by the presence of the leptin receptor (Lepr) allele, Lepr^db/db. Furthermore, DBA2/J (DBA) mice are also susceptible to β-cell failure and share 25% of their genome with BLKS; thus the DBA genome may contribute to β-cell dysfunction in BLKS mice.

Results

Here we show that BLKS mice exhibit elevated insulin secretion, as evidenced by improved glucose tolerance and increased islet insulin secretion compared with B6 mice, and describe interstrain transcriptional differences in glucose response. Transcriptional differences between BLKS and B6 mice were identified by expression profiling of isolated islets from both strains. Genomic mapping of gene expression differences demonstrated a significant association of expression differences with DBA loci in BLKS mice (P = 4×10^-27).

Conclusion

Two genes, Nicotinamide nucleotide transhydrogenase (Nnt) and Pleiomorphic adenoma gene like 1 (Plagl1), were 4 and 7.2-fold higher respectively in BLKS islets, and may be major contributors to increased insulin secretion by BLKS islets. Contrary to reports for B6 mice, BLKS mice do not harbor a mutant Nnt gene. We detected 16 synonymous polymorphisms and a two-amino acid deletion in the Plagl1 gene in BLKS mice. Several inflammatory glucose-responsive genes are expressed at a higher level in BLKS, suggesting an inflammatory component to BLKS islet dysfunction. This study describes physiological differences between BLKS and B6 mice, and provides evidence for a causative role of the DBA genome in β-cell dysfunction in BLKS mice.

Exon and junction microarrays detect widespread mouse strain- and sex-bias expression differences

Background

Studies have shown that genetic and sex differences strongly influence gene expression in mice. Given the diversity and complexity of transcripts produced by alternative splicing, we sought to use microarrays to establish the extent of variation found in mouse strains and genders. Here, we surveyed the effect of strain and sex on liver gene and exon expression using male and female mice from three different inbred strains.

Results

71 liver RNA samples from three mouse strains – DBA/2J, C57BL/6J and C3H/HeJ – were profiled using a custom-designed microarray monitoring exon and exon-junction expression of 1,020 genes representing 9,406 exons. Gene expression was calculated via two different methods, using the 3'-most exon probe ("3' gene expression profiling") and using all probes associated with the gene ("whole-transcript gene expression profiling"), while exon expression was determined using exon probes and flanking junction probes that spanned across the neighboring exons ("exon expression profiling"). Widespread strain and sex influences were detected using a two-way Analysis of Variance (ANOVA) regardless of the profiling method used. However, over 90% of the genes identified in 3' gene expression profiling or whole transcript profiling were identified in exon profiling, along with 75% and 38% more genes, respectively, showing evidence of differential isoform expression. Overall, 55% and 32% of genes, respectively, exhibited strain- and sex-bias differential gene or exon expression.

Conclusion

Exon expression profiling identifies significantly more variation than both 3' gene expression profiling and whole-transcript gene expression profiling. A large percentage of genes that are not differentially expressed at the gene level demonstrate exon expression variation suggesting an influence of strain and sex on alternative splicing and a need to profile expression changes at sub-gene resolution.

Hepatic lipin 1beta expression is diminished in insulin-resistant obese subjects and is reactivated by marked weight loss

OBJECTIVE

Lipin 1 plays critical roles in controlling energy metabolism. We sought to determine the expression of lipin 1 isoforms (lipin 1alpha and -beta) in liver and adipose tissue of obese subjects and to evaluate cellular mechanisms involved in the regulation of lipin 1 expression by physiologic stimuli.

RESEARCH DESIGN AND METHODS

The expression of lipin 1alpha and -beta was quantified in liver and adipose tissue of extremely obese (average BMI 60.8 kg/m(2)) human subjects undergoing gastric bypass surgery (GBS). Second, the expression of lipin 1 was evaluated in HepG2 cells in response to overexpression of peroxisome proliferator-activated receptor-gamma coactivator (PGC)-1alpha under normal or hyperinsulinemic conditions.

RESULTS

The expression of lipin 1beta in liver and adipose tissue was inversely related to BMI, fasting plasma insulin concentration, and the homeostasis model assessment of insulin resistance but was significantly increased by marked weight loss and insulin sensitization following GBS. Hepatic lipin 1beta mRNA levels were strongly correlated with the expression of PGC-1alpha, and overexpression of PGC-1alpha in HepG2 cells increased lipin 1 expression. Conversely, hyperinsulinemic culture conditions downregulated the expression of lipin 1beta, PGC-1alpha, and their known target genes involved in mitochondrial metabolism in HepG2 cells. Finally, overexpression of lipin 1beta or PGC-1alpha reversed the effect of hyperinsulinemia on the expression of their target genes.

CONCLUSIONS

These studies suggest that hepatic lipin 1beta and PGC-1alpha expression are downregulated by obesity and obesity-related metabolic perturbations in human subjects, likely due to alterations in insulin concentration or sensitivity.

Genome-level analysis of genetic regulation of liver gene expression networks

The liver is the primary site for the metabolism of nutrients, drugs, and chemical agents. Although metabolic pathways are complex and tightly regulated, genetic variation among individuals, reflected in variations in gene expression levels, introduces complexity into research on liver disease. This study dissected genetic networks that control liver gene expression through the combination of large-scale quantitative mRNA expression analysis with genetic mapping in a reference population of BXD recombinant inbred mouse strains for which extensive single-nucleotide polymorphism, haplotype, and phenotypic data are publicly available. We profiled gene expression in livers of naive mice of both sexes from C57BL/6J, DBA/2J, B6D2F1, and 37 BXD strains using Agilent oligonucleotide microarrays. These data were used to map quantitative trait loci (QTLs) responsible for variations in the expression of about 19,000 transcripts. We identified polymorphic local and distant QTLs, including several loci that control the expression of large numbers of genes in liver, by comparing the physical transcript position with the location of the controlling QTL.

CONCLUSION

The data are available through a public web-based resource (www.genenetwork.org) that allows custom data mining, identification of coregulated transcripts and correlated phenotypes, cross-tissue, and cross-species comparisons, as well as testing of a broad array of hypotheses.

Mouse models of diabetic neuropathy

Diabetic neuropathy (DN) is a debilitating complication of type 1 and type 2 diabetes. Rodent models of DN do not fully replicate the pathology observed in human patients. We examined DN in streptozotocin (STZ)-induced [B6] and spontaneous type 1 diabetes [B6Ins2(Akita)] and spontaneous type 2 diabetes [B6-db/db, BKS-db/db]. Despite persistent hyperglycemia, the STZ-treated B6 and B6Ins2(Akita) mice were resistant to the development of DN. In contrast, DN developed in both type 2 diabetes models: the B6-db/db and BKS-db/db mice. The persistence of hyperglycemia and development of DN in the B6-db/db mice required an increased fat diet while the BKS-db/db mice developed severe DN and remained hyperglycemic on standard mouse chow. Our data support the hypothesis that genetic background and diet influence the development of DN and should be considered when developing new models of DN.

The genetic landscape of type 2 diabetes in mice

Inbred mouse strains provide genetic diversity comparable to that of the human population. Like humans, mice have a wide range of diabetes-related phenotypes. The inbred mouse strains differ in the response of their critical physiological functions, such as insulin sensitivity, insulin secretion, beta-cell proliferation and survival, and fuel partitioning, to diet and obesity. Most of the critical genes underlying these differences have not been identified, although many loci have been mapped. The dramatic improvements in genomic and bioinformatics resources are accelerating the pace of gene discovery. This review describes how mouse genetics can be used to discover diabetes-related genes, summarizes how the mouse strains differ in their diabetes-related phenotypes, and describes several examples of how loci identified in the mouse may directly relate to human diabetes.

Increased lipids in non-lipogenic tissues are indicators of the severity of type 2 diabetes in mice

We hypothesised that the molecular changes triggered in type 2 diabetes might cause phenotypic changes in the lipid fraction of tissues. We compared tissue lipid profiles of inbred lean B6-Bom with those of the obese B6-ob/ob and diabetic BKS-db/db mice and found that genetically diabetic mice significantly accumulate fat (especially monounsaturated fatty acids, MUFA) in non-lipogenic tissues such as the eye (MUFA, 2-fold), skeletal muscle (MUFA, 13-fold) and pancreas (MUFA, 16-fold). In contrast, the B6-ob/ob mice which manifest a milder form of type 2 diabetes use the liver as their predominant lipid depot (MUFA 91-fold increase, as compared to lean mice values). The lipids in the BKS-db/db skeletal muscle and pancreas were also significantly enriched with linoleic acid (LA, (9-fold and 6-fold, respectively); and alpha-linolenic acid (ALA, 8.5-fold and 8-fold, respectively). MUFA, LA and ALA accumulation in the non-lipogenic tissues of BKS-db/db mice was associated with reduced liver stearoyl-CoA desaturase-1 expression.

The idd4 locus displays sex-specific epistatic effects on type 1 diabetes susceptibility in nonobese diabetic mice

The nonobese diabetic (NOD) mouse recapitulates many aspects of the pathogenesis of type 1 diabetes in humans, including inheritance as a complex trait. More than 20 Idd loci have been linked to type 1 diabetes susceptibility in NOD mice. Previously, we used linkage analysis of NOD crossed to the nonobese diabetes-resistant (NOR) strain and NOD congenic strains to map susceptibility to both spontaneous and cyclophosphamide-accelerated type 1 diabetes to the Idd4 locus on chromosome 11 that displayed a sex-specific effect on diabetes susceptibility. Here, we elucidate the complex genetic architecture of Idd4 by analysis of congenic strains on the NOD and NOR backgrounds. We previously refined Idd4.1 to 1.4 Mb and demonstrated an impact of this interval on type 1 interferon pathways in antigen-presenting cells. Here, we identify a second subregion, the 0.92 Mb Idd4.2 locus located telomeric to Idd4.1. Strikingly, Idd4.2 displayed a sex-specific, epistatic interaction with Idd4.1 in NOR.NOD congenic females that was not observed in syngenic males. Idd4.2 contains 29 genes, and promising candidates for the Idd4.2 effect on type 1 diabetes are described. These data demonstrate sex-dependent interaction effects on type 1 diabetes susceptibility and provide a framework for functional analysis of Idd4.2 candidate genes.

Impact of chromosome 2 obesity loci on cardiovascular complications of insulin resistance in LDL receptor-deficient C57BL/6 mice

Previous characterization of mouse chromosome 2 identified genomic intervals that influence obesity, insulin resistance, and dyslipidemia. For this, resistant CAST/Ei (CAST) alleles were introgressed onto a susceptible C57BL/6J background to generate congenic strains with CAST alleles encompassing 67-162 Mb (multigenic obesity 6 [MOB6]) and 84-180 Mb (MOB5) from mouse chromosome 2. To examine the effects of each congenic locus on atherosclerosis and glucose disposal, we bred each strain onto a sensitizing LDL receptor-null (LDLR(-/-)) C57BL/6J background to predispose them to hypercholesterolemia and insulin resistance. LDLR(-/-) congenics and controls were characterized for measures of atherogenesis, insulin sensitivity, and obesity. We identified a genomic interval unique to the MOB6 congenic (72-84 Mb) that dramatically decreased atherosclerosis by approximately threefold and decreased insulin resistance. This region also reduced adiposity twofold. Conversely, the congenic region unique to MOB5 (162-180 Mb) increased insulin resistance but had little effect on atherosclerosis and adiposity. The MOB congenic intervals are concordant to human and rat quantitative trait loci influencing diabetes and atherosclerosis traits. Thus, our results define a strategy for studying the poorly understood interactions between diabetes and atherosclerosis and for identifying genes underlying the cardiovascular complications of insulin resistance.

Integrating genetic and gene expression data: application to cardiovascular and metabolic traits in mice

The millions of common DNA variations that occur in the human population, or among inbred strains of mice and rats, perturb the expression (transcript levels) of a large fraction of the genes expressed in a particular tissue. The hundreds or thousands of common cis-acting variations that occur in the population may in turn affect the expression of thousands of other genes by affecting transcription factors, signaling molecules, RNA processing, and other processes that act in trans. The levels of transcripts are conveniently quantitated using expression arrays, and the cis- and trans-acting loci can be mapped using quantitative trait locus (QTL) analysis, in the same manner as loci for physiologic or clinical traits. Thousands of such expression QTL (eQTL) have been mapped in various crosses in mice, as well as other experimental organisms, and less detailed maps have been produced in studies of cells from human pedigrees. Such an integrative genetics approach (sometimes referred to as “genetical genomics”) is proving useful for identifying genes and pathways that contribute to complex clinical traits. The coincidence of clinical trait QTL and eQTL can help in the prioritization of positional candidate genes. More importantly, mathematical modeling of correlations between levels of transcripts and clinical traits in genetic crosses can allow prediction of causal interactions and the identification of “key driver” genes. An important objective of such studies will be to model biological networks in physiologic processes. When combined with high-density single nucleotide polymorphism (SNP) mapping, it should be feasible to identify genes that contribute to transcript levels using association analysis in outbred populations. In this review we discuss the basic concepts and applications of this integrative genomic approach to cardiovascular and metabolic diseases.

Molecular Genetic Analysis of the Idd4 Locus Implicates the IFN Response in Type 1 Diabetes Susceptibility in Nonobese Diabetic Mice

High-resolution mapping and identification of the genes responsible for type 1 diabetes (T1D) has proved difficult because of the multigenic etiology and low penetrance of the disease phenotype in linkage studies. Mouse congenic strains have been useful in refining Idd susceptibility loci in the NOD mouse model and providing a framework for identification of genes underlying complex autoimmune syndromes. Previously, we used NOD and a nonobese diabetes-resistant strain to map the susceptibility to T1D to the Idd4 locus on chromosome 11. Here, we report high-resolution mapping of this locus to 1.4 megabases. The NOD Idd4 locus was fully sequenced, permitting a detailed comparison with C57BL/6 and DBA/2J strains, the progenitors of T1D resistance alleles found in the nonobese diabetes-resistant strain. Gene expression arrays and quantitative real-time PCR were used to prioritize Idd4 candidate genes by comparing macrophages/dendritic cells from congenic strains where allelic variation was confined to the Idd4 interval. The differentially expressed genes either were mapped to Idd4 or were components of the IFN response pathway regulated in trans by Idd4. Reflecting central roles of Idd4 genes in Ag presentation, arachidonic acid metabolism and inflammation, phagocytosis, and lymphocyte trafficking, our combined analyses identified Alox15, Alox12e, Psmb6, Pld2, and Cxcl16 as excellent candidate genes for the effects of the Idd4 locus.

Genetic analysis of the diabetes-prone C57BLKS/J mouse strain reveals genetic contribution from multiple strains

The C57BLKS/J (BKS) inbred mouse strain is a widely used animal model of type 2 diabetes. In the presence of the diabetes (db) mutation, obese BKS-db mice develop severe diabetes. Genetic studies of diabetes-susceptibility in this strain are facilitated by the fact that BKS is a genetic composite between the diabetes-resistant C57BL/6J (B6) and susceptible DBA/2J (DBA) strains. On this basis, it has been hypothesized that diabetes-susceptibility in BKS is conferred by DBA-derived alleles. However, recent studies revealed non-B6/non-DBA genetic material in BKS. To identify the origin of this genetic component, we generated a genomic map of BKS using 537 microsatellite markers. Our results demonstrate that, in addition to B6 and DBA, BKS contains alleles from at least three other strains, including 129, C57BL/10 and an unidentified mouse strain. We also analyzed two congenic strains, B6-db and BKS-db, which are widely used for the genetic mapping of diabetes-susceptibility loci. We identified several donor-derived genomic regions introduced during the generation of these congenic strains. In summary, our study reveals novel aspects of the genetic fine-structure of BKS and related strains and facilitates the identification of diabetes-susceptibility loci in this mouse model.