Mapping of mouse obesity genes: A generic approach to a complex trait

Effect of an obesogenic high-fat and high-sucrose diet on hepatic gene expression signatures in male Collaborative Cross mice

Nonalcoholic fatty liver disease (NAFLD), the most prevalent chronic liver disease, is characterized by substantial variations in case-level severity. In this study, we used a genetically diverse Collaborative Cross (CC) mouse population model to analyze the global transcriptome and clarify the molecular mechanisms involved in hepatic fat accumulation that determine the level and severity of NAFLD. Twenty-four strains of male CC mice were maintained on a high-fat/high-sucrose (HF/HS) diet for 12 wk, and their hepatic gene expression profiles were determined by next-generation RNA sequencing. We found that the development of the nonalcoholic fatty liver (NAFL) phenotype in CC mice coincided with significant changes in the expression of hepatic genes at the population level, evidenced by the presence of 724 differentially expressed genes involved in lipid and carbohydrate metabolism, cell morphology, vitamin and mineral metabolism, energy production, and DNA replication, recombination, and repair. Importantly, expression of 68 of these genes strongly correlated with the extent of hepatic lipid accumulation in the overall population of HF/HS diet-fed male CC mice. Results of partial least squares (PLS) modeling showed that these derived hepatic gene expression signatures help to identify the individual mouse strains that are highly susceptible to the development of NAFLD induced by an HF/HS diet. These findings imply that gene expression profiling, combined with a PLS modeling approach, may be a useful tool to predict NAFLD severity in genetically diverse patient populations.NEW & NOTEWORTHY Feeding male Collaborative Cross mice an obesogenic diet allows modeling NAFLD at the population level. The development of NAFLD coincided with significant hepatic transcriptomic changes in this model. Genes (724) were differentially expressed and expression of 68 genes strongly correlated with the extent of hepatic lipid accumulation. Partial least squares modeling showed that derived hepatic gene expression signatures may help to identify individual mouse strains that are highly susceptible to the development of NAFLD.

Chow fed UC Davis strain female Lepr fatty Zucker rats exhibit mild glucose intolerance, hypertriglyceridemia, and increased urine volume, all reduced by a Brown Norway strain chromosome 1 congenic donor region

Our objective is to identify genes that influence the development of any phenotypes of type 2 diabetes (T2D) or kidney disease in obese animals. We use the reproductively isolated UC Davis fatty Zucker strain rat model in which the defective chromosome 4 leptin receptor (Lepr^faSte/faSte) results in fatty obesity. We previously produced a congenic strain with the distal half of chromosome 1 from the Brown Norway strain (BN) on a Zucker (ZUC) background (BN.ZUC-D1Rat183–D1Rat90). Previously published studies in males showed that the BN congenic donor region protects from some phenotypes of renal dysfunction and T2D. We now expand our studies to include females and expand phenotyping to gene expression. We performed diabetes and kidney disease phenotyping in chow-fed females of the BN.ZUC-D1Rat183-D1Rat90 congenic strain to determine the specific characteristics of the UC Davis model. Fatty Lepr^faSte/faSte animals of both BN and ZUC genotype in the congenic donor region had prediabetic levels of fasting blood glucose and blood glucose 2 hours after a glucose tolerance test. We observed significant congenic strain chromosome 1 genotype effects of the BN donor region in fatty females that resulted in decreased food intake, urine volume, glucose area under the curve during glucose tolerance test, plasma triglyceride levels, and urine glucose excretion per day. In fatty females, there were significant congenic strain BN genotype effects on non-fasted plasma urea nitrogen, triglyceride, and creatinine. Congenic region genotype effects were observed by quantitative PCR of mRNA from the kidney for six genes, all located in the chromosome 1 BN donor region, with potential effects on T2D or kidney function. The results are consistent with the hypothesis that the BN genotype chromosome 1 congenic region influences traits of both type 2 diabetes and kidney function in fatty UC Davis ZUC females and that there are many positional candidate genes.

Analysis of candidate colitis genes in the Gdac1 locus of mice deficient in glutathione peroxidase-1 and -2

Background

Mice that are deficient for glutathione peroxidases 1 and 2 (GPX) show large variations in the penetrance and severity of colitis in C57BL/6J and 129S1/SvImJ backgrounds. We mapped a locus contributing to this difference to distal chromosome 2 (∼119–133 mbp) and named it glutathione peroxidase-deficiency-associated colitis 1 (Gdac1). The aim of this study was to identify the best gene candidates within the Gdac1 locus contributing to the murine colitis phenotype.

Method/Principal Findings

We refined the boundaries of Gdac1 to 118–125 mbp (95% confidence interval) by increasing sample size and marker density across the interval. The narrowed region contains 128 well-annotated protein coding genes but it excludes Fermt1, a human inflammatory bowel disease candidate that was within the original boundaries of Gdac1. The locus we identified may be the Cdcs3 locus mapped by others studying IL10-knockout mice. Using in silico analysis of the 128 genes, based on published colon expression data, the relevance of pathways to colitis, gene mutations, presence of non-synonymous-single-nucleotide polymorphisms (nsSNPs) and whether the nsSNPs are predicted to have an impact on protein function or expression, we excluded 42 genes. Based on a similar analysis, twenty-five genes from the remaining 86 genes were analyzed for expression-quantitative-trait loci, and another 15 genes were excluded.

Conclusion/Significance

Among the remaining 10 genes, we identified Pla2g4f and Duox2 as the most likely colitis gene candidates, because GPX metabolizes PLA2G4F and DUOX2 products. Pla2g4f is a phospholipase A2 that has three potentially significant nsSNP variants and showed expression differences across mouse strains. PLA2G4F produces arachidonic acid, which is a substrate for lipoxygenases and, in turn, for GPXs. DUOX2 produces H₂O₂ and may control microbial populations. DUOX-1 and -2 control microbial populations in mammalian lung and in the gut of several insects and zebrafish. Dysbiosis is a phenotype that differentiates 129S1/SvImJ from C57BL/6J and may be due to strain differences in DUOX2 activity.

Current understanding of the genetic basis for physical activity

Although it is well known that physical activity prevents and ameliorates a large number of conditions and chronic diseases, it is also incontrovertible that physical inactivity is becoming more prevalent. This paradox has led some to suggest that genetic/biological factors influence activity levels as opposed to the classical notion that voluntary activity is solely regulated by environmental factors. There is a plethora of recent data showing that there is considerable genetic influence on activity levels in both humans and animals and emerging evidence suggesting potential genomic locations for those genetic factors. Several independent lines of evidence suggest that dopamine receptor 1 (Drd1) and nescient helix loop helix (Nhlh2) are excellent candidate genes for the regulation of physical activity, with several other potential candidate genes only partially supported. This foundation provides the basis for continuing work to identify additional candidate genes, to identify other genetic factors that are involved in the regulation of physical activity, and to investigate the mechanisms by which these genes and genetic factors regulate activity.

The contribution of animal models to the study of obesity

Obesity results from prolonged imbalance of energy intake and energy expenditure. Animal models have provided a fundamental contribution to the historical development of understanding the basic parameters that regulate the components of our energy balance. Five different types of animal model have been employed in the study of the physiological and genetic basis of obesity. The first models reflect single gene mutations that have arisen spontaneously in rodent colonies and have subsequently been characterized. The second approach is to speed up the random mutation rate artificially by treating rodents with mutagens or exposing them to radiation. The third type of models are mice and rats where a specific gene has been disrupted or over-expressed as a deliberate act. Such genetically-engineered disruptions may be generated through the entire body for the entire life (global transgenic manipulations) or restricted in both time and to certain tissue or cell types. In all these genetically-engineered scenarios, there are two types of situation that lead to insights: where a specific gene hypothesized to play a role in the regulation of energy balance is targeted, and where a gene is disrupted for a different purpose, but the consequence is an unexpected obese or lean phenotype. A fourth group of animal models concern experiments where selective breeding has been utilized to derive strains of rodents that differ in their degree of fatness. Finally, studies have been made of other species including non-human primates and dogs. In addition to studies of the physiological and genetic basis of obesity, studies of animal models have also informed us about the environmental aspects of the condition. Studies in this context include exploring the responses of animals to high fat or high fat/high sugar (Cafeteria) diets, investigations of the effects of dietary restriction on body mass and fat loss, and studies of the impact of candidate pharmaceuticals on components of energy balance. Despite all this work, there are many gaps in our understanding of how body composition and energy storage are regulated, and a continuing need for the development of pharmaceuticals to treat obesity. Accordingly, reductions in the use of animal models, while ethically desirable, will not be feasible in the short to medium term, and indeed an expansion in activity using animal models is anticipated as the epidemic continues and spreads geographically.

Fine-mapping gene-by-diet interactions on chromosome 13 in a LG/J x SM/J murine model of obesity

Obesity is one of the most serious threats to human health today. Although there is general agreement that environmental factors such as diet have largely caused the current obesity pandemic, the environmental changes have not affected all individuals equally. To model gene-by-environment interactions in a mouse model system, our group has generated an F(16) advanced intercross line (AIL) from the SM/J and LG/J inbred strains. Half of our sample was fed a low-fat (15% energy from fat) diet while the other half was fed a high-fat (43% energy from fat) diet. The sample was assayed for a variety of obesity- and diabetes-related phenotypes such as growth rate, response to glucose challenge, organ and fat pad weights, and serum lipids and insulin. An examination in the F(16) sample of eight adiposity quantitative trait loci previously identified in an F(2) intercross of SM/J and LG/J mouse strains reveals locus-by-diet interactions for all previously mapped loci. Adip7, located on proximal chromosome 13, demonstrated the most interactions and therefore was selected for fine mapping with microsatellite markers. Three phenotypic traits, liver weight in male animals, serum insulin in male animals, and reproductive fat pad weight, show locus-by-diet interactions in the 127-kb region between markers D13Mit1 and D13Mit302. The phosphofructokinase (PFK) C (Pfkp) and the pitrilysin metalloprotease 1 (Pitrm1) genes are compelling positional candidate genes in this region that show coding sequence differences between the parental strains in functional domains.

Genetic evidence for discordance between obesity- and diabetes-related traits in the LGXSM recombinant inbred mouse strains

Obesity and its comorbidities, particularly type 2 diabetes, have become serious public health problems over the past few decades. Although the current pandemic is largely caused by societal environmental changes in diet, variation in response to these changes have, in part, a genetic basis. Here we address the genetic basis for both obesity- and diabetes-related traits themselves and dietary fat responses for these traits in a set of recombinant inbred mouse strains formed from the cross of LG/J with SM/J (LGXSM lines) fed a standard low-fat (15% calories from fat) or high-fat (42% calories from fat) diet. We found substantial genetic variation for most of the traits studied. Weight at time of death, liver weight, and weight of the reproductive fat pad had especially high heritabilities, whereas heart weight and serum levels of free fatty acids and triglycerides had low heritabilities. Genetic correlations were very high among fat pad weights and serum leptin, indicating shared genetic variation between fat levels and hormonal appetite control. These obesity traits were moderately correlated with adult growth, liver weight, and serum insulin and cholesterol levels. A majority of traits also displayed genetic variation in response to a high-fat diet, especially the weight of the reproductive and renal fat pads as well as the liver. Genetic correlations in dietary response followed a pattern similar to that found for the traits themselves. Several strains manifested discordant responses for obesity, glucose, and insulin, consistent with the presence of genotypes protective for diabetes in the presence of obesity. These recombinant inbred strains represent potentially valuable new models for dissecting the complex physiological relationships among obesity and diabetes.

Immunity to tuberculosis

Only 5 to 10% of immunocompetent humans are susceptible to tuberculosis, and over 85% of them develop the disease exclusively in the lungs. Human immunodeficiency virus (HIV)-infected humans, in contrast, can develop systemic disease that is more quickly lethal. This is in keeping with other evidence showing that susceptible humans generate some level of Th1 immunity to Mycobacterium tuberculosis (Mtb) infection. Tuberculosis in mice is also exclusively a lung disease that is progressive and lethal, in spite of the generation of Th1-mediated immunity. Thus mouse tuberculosis is a model of tuberculosis in susceptible humans, as is tuberculosis in guinea pigs and rabbits. Inability to resolve infection and prevent disease may not be a consequence of the generation of an inadequate number of Th1 cells but of an intrinsic deficiency in macrophage function that prevents these cells from expressing immunity. If this proves to be true, vaccinating susceptible humans against tuberculosis will be a difficult task.

A large-sample QTL study in mice: II. Body composition

Using lines of mice having undergone long-term selection for high and low growth, a large-sample (n = approximately 1,000 F2) experiment was conducted to gain further understanding of the genetic architecture of complex polygenic traits. Composite interval mapping on data from male F2 mice (n = 552) detected 50 QTL on 15 chromosomes impacting weights of various organ and adipose subcomponents of growth, including heart, liver, kidney, spleen, testis, and subcutaneous and epididymal fat depots. Nearly all aggregate growth QTL could be interpreted in terms of the organ and fat subcomponents measured. More than 25% of QTL detected map to MMU2, accentuating the relevance of this chromosome to growth and fatness in the context of this cross. Regions of MMU7, 15, and 17 also emerged as important obesity "hot-spots." Average degrees of directional dominance are close to additivity, matching expectations for body composition traits. A strong QTL congruency is evident among heart, liver, kidney, and spleen weights. Liver and testis are organs whose genetic architectures are, respectively, most and least aligned with that for aggregate body weight. In this study, growth and body weight are interpreted in terms of organ subcomponents underlying the macro aggregate traits, and anchored on the corresponding genomic locations.

Diet, obesity, and hyperglycemia in LG/J and SM/J mice

OBJECTIVE

To examine the differential response of obesity- and diabetes-related traits to a high- or low-fat diet in LG/J and SM/J mice. We also examined food consumption in these strains.

RESEARCH METHODS AND PROCEDURES

Mice were placed on a high- or low-fat diet after weaning. Animals were weighed once per week and subjected to glucose tolerance tests at 20 weeks. At sacrifice, fat pads and internal organs were removed along with serum samples. For food consumption, LG/J and SM/J mice of each sex were assigned to a high-fat or low-fat diet after reaching maturity. Mice were weighed three times per week, and food consumed was determined by subtraction.

RESULTS

LG/J animals consume more total food, but SM/J animals consume more food per gram of body weight. LG/J mice grow faster to 10 weeks but slower from 10 to 20 weeks, have higher cholesterol and free fatty acid levels, and have lower basal glucose levels and better response to a glucose challenge than SM/J mice. For most traits, SM/J mice respond more strongly to a high-fat diet than LG/J mice, including body weight and growth, basal glucose levels, organ weights, fat distribution, and circulating triglycerides and cholesterol levels.

DISCUSSION

Obesity-related phenotypes, as well as response to increased dietary fat, differ genetically between LG/J and SM/J and can, therefore, be mapped. This study indicates that the cross of SM/J and LG/J mice would be an excellent model system for the study of gene-by-diet interaction in obesity.

Susceptibility and negative epistatic loci contributing to type 2 diabetes and related phenotypes in a KK/Ta mouse model

The KK/Ta mouse strain serves as a suitable polygenic model for human type 2 diabetes. Using 93 microsatellite markers in 208 KK/Ta x (BALB/c x KK/Ta)F1 male backcross mice, we carried out a genome-wide linkage analysis of KK/Ta alleles contributing to type 2 diabetes and related phenotypes, such as obesity and dyslipidemia. We identified three major chromosomal intervals significantly contributing to impaired glucose metabolism: one quantitative trait locus for impaired glucose tolerance on chromosome 6 and two loci for fasting blood glucose levels on chromosomes 12 and 15. The latter two loci appeared to act in a complementary fashion. Two intervals showed significant linkages for serum triglyceride levels, one on chromosome 4 and the other on chromosome 8. The KK allele on chromosome 8 acts to promote serum triglyceride levels, whereas the KK allele on chromosome 4 acts to suppress this effect in a recessive fashion. In addition, it is suggested that the chromosome 4 locus also acts to downregulate body weight and that the chromosome 8 locus acts to upregulate serum insulin levels. Our data clearly showed that each disease phenotype of type 2 diabetes and related disorders in KK/Ta mice is under the control of separate genetic mechanisms. However, there appear to be common genes contributing to different disease phenotypes. There are potentially important candidate genes that may be relevant to the disease.

Genetics of body-weight regulation

The role of genetics in obesity is twofold. Studying rare mutations in humans and model organisms provides fundamental insight into a complex physiological process, and complements population-based studies that seek to reveal primary causes. Remarkable progress has been made on both fronts, and the pace of advance is likely to accelerate as functional genomics and the human genome project expand and mature. Approaches based on mendelian and quantitative genetics may well converge, and lead ultimately to more rational and selective therapies.

Developing a new model for non-insulin dependent diabetes mellitus (NIDDM) by using the Philippine wild mouse, Mus musculus castaneus

The Philippine wild-caught castaneus mouse (Mus musculus castaneus) and laboratory mouse (C57BL/6J: B6) were used to develop a new non-insulin dependent diabetes mellitus (NIDDM) model. Offspring from the cross between a wild male and B6 female were backcrossed to the sire. One male which exhibited highest fasting hyperglycemia (190 mg/dl) among eighty-seven backcross offspring was selected at 10 weeks of age, and crossed with a B6 female to comprise the fundamental stock (F0). Thereafter, full-sib mating was performed to develop a new inbred strain named CBD (Castaneus-B6 diabetic) mouse. Mice with relatively higher fasting hyperglycemia among F0 and F1 generations were selected for breeding. From the F2 generation, mice were defined as diabetic when blood glucose levels exceeded 200 mg/dl at 120 min in intraperitoneal glucose tolerance test (IPGTT) at 10 weeks of age, and have been selectively bred. The incidence of diabetic males from the F3-F6 generation fluctuated 45-75% at 10 weeks of age and 59-72% at 20 weeks of age. Diabetic males had about two-fold higher fasting glucose and insulin levels than B6 males. Glucose-stimulated insulin secretion was impaired in diabetic CBD mice compared to B6 males at 20 weeks. Moreover, diabetic mice had slight obesity compared to B6 mice. These facts indicated that diabetic features of CBD mice resemble NIDDM in humans. The CBD strain, characterized by high incidence and early onset of diabetes with mild obesity would be of value as a new NIDDM model. The method, utilizing wild castaneus mouse of different origin from laboratory mice, maybe useful in the development of other animal models.

Genetic background determines the extent of atherosclerosis in ApoE-deficient mice

Two strains of ApoE-deficient mice were found to have markedly different plasma lipoprotein profiles and susceptibility to atherosclerosis when fed either a low-fat chow or a high-fat Western-type diet. FVB/NJ ApoE-deficient (FVB E0) mice had higher total cholesterol, HDL cholesterol, ApoA1, and ApoA2 levels when compared with C57BL/6J ApoE-deficient (C57 E0) mice. At 16 weeks of age, mean aortic root atherosclerotic lesion area was 7- to 9-fold higher in chow diet-fed C57 E0 mice and 3.5-fold higher in Western diet-fed C57 E0 mice compared with FVB E0 mice fed similar diets. Lesion area in chow diet-fed first-generation mice from a strain intercross was intermediate in size compared with parental values. The distribution of the lesion area in 150 chow diet-fed second-generation progeny spanned the range of the lesion area in both parental strains. There were no correlations between total cholesterol, non-HDL cholesterol, HDL cholesterol, ApoA1, ApoA2, ApoJ, or anti-cardiolipin antibodies and lesion area in the second-generation progeny. Thus, a genomic approach may succeed in identifying the genes responsible for the variation in atherosclerosis susceptibility in these 2 strains of ApoE-deficient mice, which could not be explained by measured plasma parameters.

Mouse models of airway responsiveness: physiological basis of observed outcomes and analysis of selected examples using these outcome indicators

The mouse is an ideal species for investigation at the interface of lung biology and lung function. As detailed in this review, there are well-developed methods for the quantitative study of lung function in mice. These methods can be applied to mice in both terminal and nonterminal experiments. Terminal experimental approaches provide more detailed physiological information, but nonterminal measurements provide adequate data for certain experiments. In this review, we provide two examples of how these models can be used to further understanding of the primary pathobiology of airway responsiveness in both the absence and the presence of induced airway inflammation. The first model is a dissection of chromosomal loci linked to the variance in airway responsiveness observed in the absence of any manipulation to induce airway inflammation. The second model explores the role of T-cell costimulatory signals in the induction of airway hyperresponsiveness. As the number of mice with targeted deletions of effector genes or insertion of informative transgenes grows, additional examples are likely to accrue.