Genetic architecture of adiposity in the cross of LG/J and SM/J inbred mice

Potential Mechanisms of Precision Nutrition-Based Interventions for Managing Obesity

Precision nutrition (PN) considers multiple individual-level and environmental characteristics or variables to better inform dietary strategies and interventions for optimizing health, including managing obesity and metabolic disorders. Here, we review the evidence on potential mechanisms-including ones to identify individuals most likely to respond-that can be leveraged in the development of PN interventions addressing obesity. We conducted a review of the literature and included laboratory, animal, and human studies evaluating biochemical and genetic data, completed and ongoing clinical trials, and public programs in this review. Our analysis describes the potential mechanisms related to 6 domains including genetic predisposition, circadian rhythms, physical activity and sedentary behavior, metabolomics, the gut microbiome, and behavioral and socioeconomic characteristics, i.e., the factors that can be leveraged to design PN-based interventions to prevent and treat obesity-related outcomes such as weight loss or metabolic health as laid out by the NIH 2030 Strategic Plan for Nutrition Research. For example, single nucleotide polymorphisms can modify responses to certain dietary interventions, and epigenetic modulation of obesity risk via physical activity patterns and macronutrient intake have also been demonstrated. Additionally, we identified limitations including questions of equitable implementation across a limited number of clinical trials. These include the limited ability of current PN interventions to address systemic influences such as supply chains and food distribution, healthcare systems, racial or cultural inequities, and economic disparities, particularly when designing and implementing PN interventions in low- and middle-income communities. PN has the potential to help manage obesity by addressing intra- and inter-individual variation as well as context, as opposed to "one-size fits all" approaches though there is limited clinical trial evidence to date.

Genetic variants of PKLR are associated with acute pain in sickle cell disease

Acute pain, the most prominent complication of sickle cell disease (SCD), results from vaso-occlusion triggered by sickling of deoxygenated red blood cells (RBCs). Concentration of 2,3-diphosphoglycerate (2,3-DPG) in RBCs promotes deoxygenation by preferentially binding to the low-affinity T conformation of HbS. 2,3-DPG is an intermediate substrate in the glycolytic pathway in which pyruvate kinase (gene PKLR, protein PKR) is a rate-limiting enzyme; variants in PKLR may affect PKR activity, 2,3-DPG levels in RBCs, RBC sickling, and acute pain episodes (APEs). We performed a candidate gene association study using 2 cohorts: 242 adult SCD-HbSS patients and 977 children with SCD-HbSS or SCD-HbSβ0 thalassemia. Seven of 47 PKLR variants evaluated in the adult cohort were associated with hospitalization: intron 4, rs2071053; intron 2, rs8177970, rs116244351, rs114455416, rs12741350, rs3020781, and rs8177964. All 7 variants showed consistent effect directions in both cohorts and remained significant in weighted Fisher's meta-analyses of the adult and pediatric cohorts using P < .0071 as threshold to correct for multiple testing. Allele-specific expression analyses in an independent cohort of 52 SCD adults showed that the intronic variants are likely to influence APE by affecting expression of PKLR, although the causal variant and mechanism are not defined.

Giant Island Mice Exhibit Widespread Gene Expression Changes in Key Metabolic Organs

Island populations repeatedly evolve extreme body sizes, but the genomic basis of this pattern remains largely unknown. To understand how organisms on islands evolve gigantism, we compared genome-wide patterns of gene expression in Gough Island mice, the largest wild house mice in the world, and mainland mice from the WSB/EiJ wild-derived inbred strain. We used RNA-seq to quantify differential gene expression in three key metabolic organs: gonadal adipose depot, hypothalamus, and liver. Between 4,000 and 8,800 genes were significantly differentially expressed across the evaluated organs, representing between 20% and 50% of detected transcripts, with 20% or more of differentially expressed transcripts in each organ exhibiting expression fold changes of at least 2×. A minimum of 73 candidate genes for extreme size evolution, including Irs1 and Lrp1, were identified by considering differential expression jointly with other data sets: 1) genomic positions of published quantitative trait loci for body weight and growth rate, 2) whole-genome sequencing of 16 wild-caught Gough Island mice that revealed fixed single-nucleotide differences between the strains, and 3) publicly available tissue-specific regulatory elements. Additionally, patterns of differential expression across three time points in the liver revealed that Arid5b potentially regulates hundreds of genes. Functional enrichment analyses pointed to cell cycling, mitochondrial function, signaling pathways, inflammatory response, and nutrient metabolism as potential causes of weight accumulation in Gough Island mice. Collectively, our results indicate that extensive gene regulatory evolution in metabolic organs accompanied the rapid evolution of gigantism during the short time house mice have inhabited Gough Island.

Genetic background and diet affect brown adipose gene coexpression networks associated with metabolic phenotypes

Adipose is a dynamic endocrine organ that is critical for regulating metabolism and is highly responsive to nutritional environment. Brown adipose tissue is an exciting potential therapeutic target; however, there are no systematic studies of gene-by-environment interactions affecting function of this organ. We leveraged a weighted gene coexpression network analysis to identify transcriptional networks in brown adipose tissue from LG/J and SM/J inbred mice fed high- or low-fat diets and correlate these networks with metabolic phenotypes. We identified eight primary gene network modules associated with variation in obesity and diabetes-related traits. Four modules were enriched for metabolically relevant processes such as immune and cytokine response, cell division, peroxisome functions, and organic molecule metabolic processes. The relative expression of genes in these modules is highly dependent on both genetic background and dietary environment. Genes in the immune/cytokine response and cell division modules are particularly highly expressed in high fat-fed SM/J mice, which show unique brown adipose-dependent remission of diabetes. The interconnectivity of genes in these modules is also heavily dependent on diet and strain, with most genes showing both higher expression and coexpression under the same context. We highlight several genes of interest, Col28a1, Cyp26b1, Bmp8b, and Ngef, that have distinct expression patterns among strain-by-diet contexts and fall under metabolic quantitative trait loci previously mapped in an F16 generation of an advanced intercross between LG/J and SM/J. Each of these genes have some connection to obesity and diabetes-related traits, but have not been studied in brown adipose tissue. Our results provide important insights into the relationship between brown adipose and systemic metabolism by being the first gene-by-environment study of brown adipose transcriptional networks.

Genome-wide Associations Reveal Human-Mouse Genetic Convergence and Modifiers of Myogenesis, CPNE1 and STC2

Muscle bulk in adult healthy humans is highly variable even after height, age, and sex are accounted for. Low muscle mass, due to fewer and/or smaller constituent muscle fibers, would exacerbate the impact of muscle loss occurring in aging or disease. Genetic variability substantially influences muscle mass differences, but causative genes remain largely unknown. In a genome-wide association study (GWAS) on appendicular lean mass (ALM) in a population of 85,750 middle-aged (aged 38-49 years) individuals from the UK Biobank (UKB), we found 182 loci associated with ALM (p < 5 × 10-8). We replicated associations for 78% of these loci (p < 5 × 10-8) with ALM in a population of 181,862 elderly (aged 60-74 years) individuals from UKB. We also conducted a GWAS on hindlimb skeletal muscle mass of 1,867 mice from an advanced intercross between two inbred strains (LG/J and SM/J); this GWAS identified 23 quantitative trait loci. Thirty-eight positional candidates distributed across five loci overlapped between the two species. In vitro studies of positional candidates confirmed CPNE1 and STC2 as modifiers of myogenesis. Collectively, these findings shed light on the genetics of muscle mass variability in humans and identify targets for the development of interventions for treatment of muscle loss. The overlapping results between humans and the mouse model GWAS point to shared genetic mechanisms across species.

Proteomic analysis of plasma from children with sickle cell anemia and silent cerebral infarction

Silent cerebral infarction is the most common neurological abnormality in children with sickle cell anemia, affecting 30-40% of 14 year olds. There are no known biomarkers to identify children with silent cerebral infarcts, and the pathological basis is also unknown. We used an unbiased proteomic discovery approach to identify plasma proteins differing in concentration between children with and without silent cerebral infarcts. Clinical parameters and plasma samples were analysed from 51 children (mean age 11.8 years, range 6-18) with sickle cell anemia (HbSS). A total of 19 children had silent cerebral infarcts and 32 normal MRI; the children with silent infarcts had lower HbF levels (8.6 vs. 16.1%, P=0.049) and higher systolic blood pressures (115 vs. 108.6, P=0.027). Plasma proteomic analysis showed 13 proteins increased more than 1.3 fold in the SCI patients, including proteins involved in hypercoagulability (α2-antiplasmin, fibrinogen−γ chain, thrombospondin-4), inflammation (α2-macroglobulin, complement C1s and C3), and atherosclerosis (apolipoprotein B-100). Higher levels of gelsolin and retinol-binding protein 4 were also found in the population with silent infarcts, both of which have been linked to stroke. We investigated the genetic basis of these differences by studying 359 adults with sickle cell disease (199 with silent cerebral infarcts, 160 normal MRIs), who had previously undergone a genome-wide genotyping array. None of the genes coding for the differentially expressed proteins were significantly associated with silent infarction. Our study suggests that silent cerebral infarcts in sickle cell anemia may be associated with higher systolic blood pressure, lower HbF levels, hypercoagulability, inflammation and atherosclerotic lipoproteins.

Epistatic effects between pairs of the growth hormone secretagogue receptor 1a, growth hormone, growth hormone receptor, non-SMC condensin I complex, subunit G and stearoyl-CoA desaturase genes on carcass, price-related and fatty acid composition traits in Japanese Black cattle

Growth hormone secretagogue receptor 1a (GHSR1a), growth hormone (GH), growth hormone receptor (GHR), non-SMC condensin I complex, subunit G (NCAPG) and stearoyl-CoA desaturase (SCD), are known to play important roles in growth and lipid metabolisms. Single and epistatic effects of the five genes on carcass, price-related and fatty acid (FA) composition traits were analyzed in a commercial Japanese Black cattle population of Ibaraki Prefecture. A total of 650 steers and 116 heifers for carcass and price-related traits, and 158 steers for FA composition traits were used in this study. Epistatic effects between pairs of the five genes were found in several traits. Alleles showing strain-specific differences in the five genes had significant single and epistatic effects in some traits. The data suggest that a TG-repeat polymorphism of the GHSR1a.5'UTR-(TG)_n locus plays a central role in gene-gene epistatic interaction of FA composition traits in the adipose tissue of Japanese Black cattle.

Kidney function is associated with an altered protein composition of high-density lipoprotein

Patients with chronic kidney disease (CKD) exhibit a myriad of metabolic derangements, including dyslipidemia characterized by low plasma concentrations of high-density lipoprotein (HDL)-associated cholesterol. However, the effects of kidney disease on HDL composition have not been comprehensively determined. Here we used a targeted mass spectrometric approach to quantify 38 proteins contained in the HDL particles within a CKD cohort of 509 participants with a broad range of estimated glomerular filtration rates (eGFRs) (CKD stages I-V, and a mean eGFR of 45.5 mL/min/1.73m²). After adjusting for multiple testing, demographics, comorbidities, medications, and other characteristics, eGFR was significantly associated with differences in four HDL proteins. Compared to participants with an eGFR of 60 mL/min/1.73m² or more, those with an eGFR under 15 mL/min/1.73m² exhibited 1.89-fold higher retinol-binding protein 4 (95% confidence interval 1.34-2.67), 1.52-fold higher apolipoprotein C-III (1.25-1.84), 0.70-fold lower apolipoprotein L1 (0.55-0.92), and 0.64-fold lower vitronectin (0.48-0.85). Although the HDL apolipoprotein L1 was slightly lower among African Americans than among Caucasian individuals, the relationship to eGFR did not differ by race. After adjustment, no HDL-associated proteins associated with albuminuria. Thus, modest changes in the HDL proteome provide preliminary evidence for an association between HDL proteins and declining kidney function, but this needs to be replicated. Future analyses will determine if HDL proteomics is indeed a clinical predictor of declining kidney function or cardiovascular outcomes.

do Monte BG, Watanabe IM, Cunha J, Peripato AC. A link between thrifty phenotype and maternal care across two generations of intercrossed mice

Maternal effects are causal influences from mother to offspring beyond genetic information, and have lifelong consequences for multiple traits. Previously, we reported that mice whose mothers did not nurse properly had low birth weight followed by rapid fat accumulation and disturbed development of some organs. That pattern resembles metabolic syndromes known collectively as the thrifty phenotype, which is believed to be an adaptation to a stressful environment which prepares offspring for reduced nutrient supply. The potential link between maternal care, stress reactivity, and the thrifty phenotype, however, has been poorly explored in the human and animal literature: only a couple of studies even mention (much less, test) these concepts under a cohesive framework. Here, we explored this link using mice of the parental inbred strains SM/J and LG/J–who differ dramatically in their maternal care–and the intercrossed generations F1 and F2. We measured individual differences in 15 phenotypes and used structural equation modeling to test our hypotheses. We found a remarkable relationship between thrifty phenotype and lower quality of maternal behaviors, including nest building, pup retrieval, grooming/licking, and nursing. To our knowledge, this is the first study to show, in any mammal, a clear connection between the natural variation in thrifty phenotype and maternal care. Both traits in the mother also had a substantial effect on survival rate in the F3 offspring. To our surprise, however, stress reactivity seemed to play no role in our models. Furthermore, the strain of maternal grandmother, but not of paternal grandmother, affected the variation of maternal care in F2 mice, and this effect was mediated by thrifty phenotype in F2. Since F1 animals were all genetically identical, this finding suggests that maternal effects pass down both maternal care and thrifty phenotype in these mice across generations via epigenetic transmission.

Modeling of a negative feedback mechanism explains antagonistic pleiotropy in reproduction in domesticated Caenorhabditis elegans strains

Most biological traits and common diseases have a strong but complex genetic basis, controlled by large numbers of genetic variants with small contributions to a trait or disease risk. The effect-size of most genetic variants is not absolute and is instead dependent upon multiple factors such as the age and genetic background of an organism. In order to understand the mechanistic basis of these changes, we characterized heritable trait differences between two domesticated strains of C. elegans. We previously identified a major effect locus, caused in part by a mutation in a component of the NURF chromatin remodeling complex, that regulates reproductive output in an age-dependent manner. The effect-size of this locus changes from positive to negative over the course of an animal’s reproductive lifespan. Here, we use a previously published macroscale model of the egg-laying rate in C. elegans to show that time-dependent effect-size is explained by an unequal use of sperm combined with negative feedback between sperm and ovulation rate. We validate key predictions of this model with controlled mating experiments and quantification of oogenesis and sperm use. Incorporation of this model into QTL mapping allows us to identify and partition new QTLs into specific aspects of the egg-laying process. Finally, we show how epistasis between two genetic variants is predicted by this modeling as a consequence of the unequal use of sperm. This work demonstrates how modeling of multicellular communication systems can improve our ability to predict and understand the role of genetic variation on a complex phenotype. Negative autoregulatory feedback loops, common in transcriptional regulation, could play an important role in modifying genetic architecture in other traits.

Complex traits are influenced by the individual effects of genetic variants in addition to the interactions of the variants with the environment, age, and each other. While complex genetic architectures are ubiquitous in natural traits, little is known about the causal mechanisms that create their complex genetic architectures. Here we identify an example of age-dependent genetic architecture controlling the rate and timing of reproduction in the hermaphroditic nematode C. elegans. We use computational modeling to demonstrate how age-dependent genetic architecture can arise as a consequence of two factors: hormonal feedback on oocytes mediated by major sperm protein (MSP) released by sperm stored in the spermatheca and life history differences in sperm use caused by genetic variants. Our work also suggests how antagonistic pleiotropy can emerge from multicellular feedback systems.

Armour GE, Mackay TFC, Anholt RRH. Epistatic partners of neurogenic genes modulate Drosophila olfactory behavior

The extent to which epistasis affects the genetic architecture of complex traits is difficult to quantify, and identifying variants in natural populations with epistatic interactions is challenging. Previous studies in Drosophila implicated extensive epistasis between variants in genes that affect neural connectivity and contribute to natural variation in olfactory response to benzaldehyde. In this study, we implemented a powerful screen to quantify the extent of epistasis as well as identify candidate interacting variants using 203 inbred wild-derived lines with sequenced genomes of the Drosophila melanogaster Genetic Reference Panel (DGRP). We crossed the DGRP lines to P[GT1]-element insertion mutants in Sema-5c and neuralized (neur), two neurodevelopmental loci which affect olfactory behavior, and to their coisogenic wild-type control. We observed significant variation in olfactory responses to benzaldehyde among F1 genotypes and for the DGRP line by mutant genotype interactions for both loci, showing extensive nonadditive genetic variation. We performed genome-wide association analyses to identify the candidate modifier loci. None of these polymorphisms were in or near the focal genes; therefore, epistasis is the cause of the nonadditive genetic variance. Candidate genes could be placed in interaction networks. Several candidate modifiers are associated with neural development. Analyses of mutants of candidate epistatic partners with neur (merry-go-round (mgr), prospero (pros), CG10098, Alhambra (Alh) and CG12535) and Sema-5c (CG42540 and bruchpilot (brp)) showed aberrant olfactory responses compared with coisogenic controls. Thus, integrating genome-wide analyses of natural variants with mutations at defined genomic locations in a common coisogenic background can unmask specific epistatic modifiers of behavioral phenotypes.

Body Composition QTLs Identified in Intercross Populations Are Reproducible in Consomic Mouse Strains

Genetic variation contributes to individual differences in obesity, but defining the exact relationships between naturally occurring genotypes and their effects on fatness remains elusive. As a step toward positional cloning of previously identified body composition quantitative trait loci (QTLs) from F2 crosses of mice from the C57BL/6ByJ and 129P3/J inbred strains, we sought to recapture them on a homogenous genetic background of consomic (chromosome substitution) strains. Male and female mice from reciprocal consomic strains originating from the C57BL/6ByJ and 129P3/J strains were bred and measured for body weight, length, and adiposity. Chromosomes 2, 7, and 9 were selected for substitution because previous F2 intercross studies revealed body composition QTLs on these chromosomes. We considered a QTL confirmed if one or both sexes of one or both reciprocal consomic strains differed significantly from the host strain in the expected direction after correction for multiple testing. Using these criteria, we confirmed two of two QTLs for body weight (Bwq5-6), three of three QTLs for body length (Bdln3-5), and three of three QTLs for adiposity (Adip20, Adip26 and Adip27). Overall, this study shows that despite the biological complexity of body size and composition, most QTLs for these traits are preserved when transferred to consomic strains; in addition, studying reciprocal consomic strains of both sexes is useful in assessing the robustness of a particular QTL.

Using whole-genome sequences of the LG/J and SM/J inbred mouse strains to prioritize quantitative trait genes and nucleotides

Background

The laboratory mouse is the most commonly used model for studying variation in complex traits relevant to human disease. Here we present the whole-genome sequences of two inbred strains, LG/J and SM/J, which are frequently used to study variation in complex traits as diverse as aging, bone-growth, adiposity, maternal behavior, and methamphetamine sensitivity.

Results

We identified small nucleotide variants (SNVs) and structural variants (SVs) in the LG/J and SM/J strains relative to the reference genome and discovered novel variants in these two strains by comparing their sequences to other mouse genomes. We find that 39% of the LG/J and SM/J genomes are identical-by-descent (IBD). We characterized amino-acid changing mutations using three algorithms: LRT, PolyPhen-2 and SIFT. We also identified polymorphisms between LG/J and SM/J that fall in regulatory regions and highly informative transcription factor binding sites (TFBS). We intersected these functional predictions with quantitative trait loci (QTL) mapped in advanced intercrosses of these two strains. We find that QTL are both over-represented in non-IBD regions and highly enriched for variants predicted to have a functional impact. Variants in QTL associated with metabolic (231 QTL identified in an F₁₆ generation) and developmental (41 QTL identified in an F₃₄ generation) traits were interrogated and we highlight candidate quantitative trait genes (QTG) and nucleotides (QTN) in a QTL on chr13 associated with variation in basal glucose levels and in a QTL on chr6 associated with variation in tibia length.

Conclusions

We show how integrating genomic sequence with QTL reduces the QTL search space and helps researchers prioritize candidate genes and nucleotides for experimental follow-up. Additionally, given the LG/J and SM/J phylogenetic context among inbred strains, these data contribute important information to the genomic landscape of the laboratory mouse.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1592-3) contains supplementary material, which is available to authorized users.

High-resolution mapping of a genetic locus regulating preferential carbohydrate intake, total kilocalories, and food volume on mouse chromosome 17

The specific genes regulating the quantitative variation in macronutrient preference and food intake are virtually unknown. We fine mapped a previously identified mouse chromosome 17 region harboring quantitative trait loci (QTL) with large effects on preferential macronutrient intake-carbohydrate (Mnic1), total kilcalories (Kcal2), and total food volume (Tfv1) using interval-specific strains. These loci were isolated in the [C57BL/6J.CAST/EiJ-17.1-(D17Mit19-D17Mit50); B6.CAST-17.1] strain, possessing a ∼ 40.1 Mb region of CAST DNA on the B6 genome. In a macronutrient selection paradigm, the B6.CAST-17.1 subcongenic mice eat 30% more calories from the carbohydrate-rich diet, ∼ 10% more total calories, and ∼ 9% more total food volume per body weight. In the current study, a cross between carbohydrate-preferring B6.CAST-17.1 and fat-preferring, inbred B6 mice was used to generate a subcongenic-derived F2 mapping population; genotypes were determined using a high-density, custom SNP panel. Genetic linkage analysis substantially reduced the 95% confidence interval for Mnic1 (encompassing Kcal2 and Tfv1) from 40.1 to 29.5 Mb and more precisely established its boundaries. Notably, no genetic linkage for self-selected fat intake was detected, underscoring the carbohydrate-specific effect of this locus. A second key finding was the separation of two energy balance QTLs: Mnic1/Kcal2/Tfv1 for food intake and a newly discovered locus regulating short term body weight gain. The Mnic1/Kcal2/Tfv1 QTL was further de-limited to 19.0 Mb, based on the absence of nutrient intake phenotypes in subcongenic HQ17IIa mice. Analyses of available sequence data and gene ontologies, along with comprehensive expression profiling in the hypothalamus of non-recombinant, cast/cast and b6/b6 F2 controls, focused our attention on candidates within the QTL interval. Zfp811, Zfp870, and Btnl6 showed differential expression and also contain stop codons, but have no known biology related to food intake regulation. The genes Decr2, Ppard and Agapt1 are more appealing candidates because of their involvement in lipid metabolism and down-regulation in carbohydrate-preferring animals.

Estimating directional epistasis

Epistasis, i.e., the fact that gene effects depend on the genetic background, is a direct consequence of the complexity of genetic architectures. Despite this, most of the models used in evolutionary and quantitative genetics pay scant attention to genetic interactions. For instance, the traditional decomposition of genetic effects models epistasis as noise around the evolutionarily-relevant additive effects. Such an approach is only valid if it is assumed that there is no general pattern among interactions—a highly speculative scenario. Systematic interactions generate directional epistasis, which has major evolutionary consequences. In spite of its importance, directional epistasis is rarely measured or reported by quantitative geneticists, not only because its relevance is generally ignored, but also due to the lack of simple, operational, and accessible methods for its estimation. This paper describes conceptual and statistical tools that can be used to estimate directional epistasis from various kinds of data, including QTL mapping results, phenotype measurements in mutants, and artificial selection responses. As an illustration, I measured directional epistasis from a real-life example. I then discuss the interpretation of the estimates, showing how they can be used to draw meaningful biological inferences.

Quantitative trait loci for energy balance traits in an advanced intercross line derived from mice divergently selected for heat loss

Obesity in human populations, currently a serious health concern, is considered to be the consequence of an energy imbalance in which more energy in calories is consumed than is expended. We used interval mapping techniques to investigate the genetic basis of a number of energy balance traits in an F₁₁ advanced intercross population of mice created from an original intercross of lines selected for increased and decreased heat loss. We uncovered a total of 137 quantitative trait loci (QTLs) for these traits at 41 unique sites on 18 of the 20 chromosomes in the mouse genome, with X-linked QTLs being most prevalent. Two QTLs were found for the selection target of heat loss, one on distal chromosome 1 and another on proximal chromosome 2. The number of QTLs affecting the various traits generally was consistent with previous estimates of heritabilities in the same population, with the most found for two bone mineral traits and the least for feed intake and several body composition traits. QTLs were generally additive in their effects, and some, especially those affecting the body weight traits, were sex-specific. Pleiotropy was extensive within trait groups (body weights, adiposity and organ weight traits, bone traits) and especially between body composition traits adjusted and not adjusted for body weight at sacrifice. Nine QTLs were found for one or more of the adiposity traits, five of which appeared to be unique. The confidence intervals among all QTLs averaged 13.3 Mb, much smaller than usually observed in an F₂ cross, and in some cases this allowed us to make reasonable inferences about candidate genes underlying these QTLs. This study combined QTL mapping with genetic parameter analysis in a large segregating population, and has advanced our understanding of the genetic architecture of complex traits related to obesity.

The effect of dietary fat intake on hepatic gene expression in LG/J AND SM/J mice

Background

The liver plays a major role in regulating metabolic homeostasis and is vital for nutrient metabolism. Identifying the genetic factors regulating these processes could lead to a greater understanding of how liver function responds to a high-fat diet and how that response may influence susceptibilities to obesity and metabolic syndrome. In this study we examine differences in hepatic gene expression between the LG/J and SM/J inbred mouse strains and how gene expression in these strains is affected by high-fat diet. LG/J and SM/J are known to differ in their responses to a high-fat diet for a variety of obesity- and diabetes-related traits, with the SM/J strain exhibiting a stronger phenotypic response to diet.

Results

Dietary intake had a significant effect on gene expression in both inbred lines. Genes up-regulated by a high-fat diet were involved in biological processes such as lipid and carbohydrate metabolism; protein and amino acid metabolic processes were down regulated on a high-fat diet. A total of 259 unique transcripts exhibited a significant diet-by-strain interaction. These genes tended to be associated with immune function. In addition, genes involved in biochemical processes related to non-alcoholic fatty liver disease (NAFLD) manifested different responses to diet between the two strains. For most of these genes, SM/J had a stronger response to the high-fat diet than LG/J.

Conclusions

These data show that dietary fat impacts gene expression levels in SM/J relative to LG/J, with SM/J exhibiting a stronger response. This supports previous data showing that SM/J has a stronger phenotypic response to high-fat diet. Based upon these findings, we suggest that SM/J and its cross with the LG/J strain provide a good model for examining non-alcoholic fatty liver disease and its role in metabolic syndrome.

Epistasis and quantitative traits: using model organisms to study gene-gene interactions

The role of epistasis in the genetic architecture of quantitative traits is controversial, despite the biological plausibility that nonlinear molecular interactions underpin the genotype-phenotype map. This controversy arises because most genetic variation for quantitative traits is additive. However, additive variance is consistent with pervasive epistasis. In this Review, I discuss experimental designs to detect the contribution of epistasis to quantitative trait phenotypes in model organisms. These studies indicate that epistasis is common, and that additivity can be an emergent property of underlying genetic interaction networks. Epistasis causes hidden quantitative genetic variation in natural populations and could be responsible for the small additive effects, missing heritability and the lack of replication that are typically observed for human complex traits.

Quantitative trait loci for bone mineral density and femoral morphology in an advanced intercross population of mice

Osteoporosis, characterized by low levels of bone mineral density (BMD), is a prevalent medical condition in humans. We investigated its genetic and environmental basis by searching for quantitative trait loci (QTLs) affecting six skeletal (including three BMD) traits in a G10 advanced intercross population produced from crosses of mice from the inbred strain C57BL/6J with mice from a strain selected for high voluntary wheel running. The mice in this population were fed either a high-fat or a matched control diet throughout the study, allowing us to test for QTL by diet interactions for the skeletal traits. Our genome scan uncovered a number of QTLs, the great majority of which were different from QTLs previously found for these same traits in an earlier (G4) generation of the same intercross. Further, the confidence intervals for the skeletal trait QTLs were reduced from an average of 18.5 Mb in the G4 population to an equivalent of about 9 Mb in the G10 population. We uncovered a total of 50 QTLs representing 32 separate genomic sites affecting these traits, with a distal region on chromosome 1 harboring several QTLs with large effects on the BMD traits. One QTL was located on chromosome 5 at 4.0 Mb with a confidence interval spanning from 4.0 to 4.6 Mb. Only three protein coding genes reside in this interval, and one of these, Cyp51, is an attractive candidate as others have shown that developing Cyp51 knockout embryos exhibit shortened and bowed limbs and synotosis of the femur and tibia. Several QTLs showed significant interactions with sex, although only two QTLs interacted with diet, both affecting only mice fed the high-fat diet.

Deconstructing Mus gemischus: advances in understanding ancestry, structure, and variation in the genome of the laboratory mouse

The laboratory mouse is an artificial construct with a complex relationship to its natural ancestors. In 2002, the mouse became the first mammalian model organism with a reference genome. Importantly, the mouse genome sequence was assembled from data on a single inbred laboratory strain, C57BL/6. Several large-scale genetic variant discovery efforts have been conducted, resulting in a catalog of tens of millions of SNPs and structural variants. High-density genotyping arrays covering a subset of those variants have been used to produce hundreds of millions of genotypes in laboratory stocks and a small number of wild mice. These landmark resources now enable us to determine relationships among laboratory mice, assign local ancestry at fine scale, resolve important controversies, and identify a new set of challenges-most importantly, the troubling scarcity of genetic data on the very natural populations from which the laboratory mouse was derived. Our aim with this review is to provide the reader with an historical context for the mouse as a model organism and to explain how practical decisions made in the past have influenced both the architecture of the laboratory mouse genome and the design and execution of current large-scale resources. We also provide examples on how the accomplishments of the past decade can be used by researchers to streamline the use of mice in their experiments and correctly interpret results. Finally, we propose future steps that will enable the mouse community to extend its successes in the decade to come.

Milk ejection in mice LG/J x SM/J

In mammals, milk provision is crucial to offspring survival and growth from birth to weaning. Milk deficiency early in life may cause death or changes in the progeny metabolism that later may lead to obesity and metabolic disorders. This study investigates milk ejection (ME) the first day after birth (D1) in F(2) females from the intercross of LG/J and SM/J inbred mice strains. The absence of milk in F(3) pups' stomach at D1 is directly associated with their survival (p < 0.001) and growth pattern (p < 0.001) in the early stages of life. Furthermore, late growth pattern is also affected by this lack of nutrients at D1 because pups that survive this absence, mostly males, are heavier at weaning (p < 0.001) which, after necropsy, is shown to be due to significant higher total fat deposition (p < 0.01). We performed QTL analysis for ME at D1 in these F(2) females. Maternal performance of ME revealed a complex genetic architecture which even though it contains only a single QTL (accounting for 8 % of the variation in ME), it is totally context-dependent on the genetic background. We discovered many regions involved in epistatic interactions that together with the single QTL explain 19 % of the genetic variation for this trait. Milk ejection is an important component of maternal care, and understanding the mechanisms modulating its variation, along with other maternal features, may help to disentangle the complexity that is the mother/offspring relationship.

Subcongenic analysis of tabw2 obesity QTL on mouse chromosome 6

Background

We previously established a congenic mouse strain with TALLYHO/Jng (TH) donor segment on chromosome 6 in a C57BL/6 (B6) background that harbors an obesity quantitative trait locus, tabw2. The B6.TH-tabw2 congenic mice developed increased adiposity that became exacerbated upon feeding a high fat-high sucrose (HFS) diet. To fine map the tabw2, in this study we generated and characterized subcongenic lines with smaller TH donor segments.

Results

We fixed four subcongenic lines, with maximum size of donor segment retained in the lines ranging from 10.8 – 92.5 Mb. For mapping, all the subcongenic mice, along with B6.TH-tabw2 congenic and B6-homozygous control mice were fed either chow or HFS diets, and their post-mortem fat pads were weighed. Mice were also characterized for energy expenditure, respiratory exchange ratio, locomotor activity, and food intake. As previously reported, B6.TH-tabw2 congenic mice showed a significantly larger fat mass than controls on both diets. On chow, a subcongenic line retaining the distal region of the TH donor congenic interval exhibited significantly larger fat mass than B6-homozygous controls, and comparable that to B6.TH-tabw2 congenic mice. Two nested subcongenic lines within that region suggested that the effect of tabw2 on obesity could be attributed to at least two subloci. On HFS diets, on the other hand, all the subcongenic mice had significantly larger fat mass than controls without genotype differences, but none of them had fat mass as large as the original congenic mice. This possibly implicates that further genetic complexity involves in the effect of tabw2 on diet-induced obesity. Significantly reduced locomotor activity was exhibited in B6.TH-tabw2 congenic and subcongenic mice compared to controls when animals were fed HFS diets. B6.TH-tabw2 congenic mice, but not subcongenic mice, also had significantly increased food intake on HFS diets.

Conclusions

It appears that at least two subloci explaining the tabw2 effect under chow feeding map to the distal region of the congenic interval, whereas the diet-induced obesity mediated by tabw2 is attributed to more complex genetic mechanism.

The Human Obesity Gene Map: The 2003 Update

This is the tenth update of the human obesity gene map, incorporating published results up to the end of October 2003 and continuing the previous format. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, quantitative trait loci (QTLs) from human genome-wide scans and animal crossbreeding experiments, and association and linkage studies with candidate genes and other markers is reviewed. Transgenic and knockout murine models relevant to obesity are also incorporated (N = 55). As of October 2003, 41 Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and causal genes or strong candidates have been identified for most of these syndromes. QTLs reported from animal models currently number 183. There are 208 human QTLs for obesity phenotypes from genome-wide scans and candidate regions in targeted studies. A total of 35 genomic regions harbor QTLs replicated among two to five studies. Attempts to relate DNA sequence variation in specific genes to obesity phenotypes continue to grow, with 272 studies reporting positive associations with 90 candidate genes. Fifteen such candidate genes are supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. Overall, more than 430 genes, markers, and chromosomal regions have been associated or linked with human obesity phenotypes. The electronic version of the map with links to useful sites can be found at http://obesitygene.pbrc.edu.

The Human Obesity Gene Map: The 2004 Update

This paper presents the eleventh update of the human obesity gene map, which incorporates published results up to the end of October 2004. Evidence from single-gene mutation obesity cases, Mendelian disorders exhibiting obesity as a clinical feature, transgenic and knockout murine models relevant to obesity, quantitative trait loci (QTLs) from animal cross-breeding experiments, association studies with candidate genes, and linkages from genome scans is reviewed. As of October 2004, 173 human obesity cases due to single-gene mutations in 10 different genes have been reported, and 49 loci related to Mendelian syndromes relevant to human obesity have been mapped to a genomic region, and causal genes or strong candidates have been identified for most of these syndromes. There are 166 genes which, when mutated or expressed as transgenes in the mouse, result in phenotypes that affect body weight and adiposity. The number of QTLs reported from animal models currently reaches 221. The number of human obesity QTLs derived from genome scans continues to grow, and we have now 204 QTLs for obesity-related phenotypes from 50 genome-wide scans. A total of 38 genomic regions harbor QTLs replicated among two to four studies. The number of studies reporting associations between DNA sequence variation in specific genes and obesity phenotypes has also increased considerably with 358 findings of positive associations with 113 candidate genes. Among them, 18 genes are supported by at least five positive studies. The obesity gene map shows putative loci on all chromosomes except Y. Overall, >600 genes, markers, and chromosomal regions have been associated or linked with human obesity phenotypes. The electronic version of the map with links to useful publications and genomic and other relevant sites can be found at http://obesitygene.pbrc.edu.

Identification of exercise capacity QTL using association mapping in inbred mice

There are large interindividual differences in exercise capacity. It is well established that there is a genetic basis for these differences. However, the genetic factors underlying this variation are undefined. Therefore, the purpose of this study was to identify novel putative quantitative trait loci (QTL) for exercise capacity by measuring exercise capacity in inbred mice and performing genome-wide association mapping. Exercise capacity, defined as run time and work, was assessed in male mice (n = 6) from 34 strains of classical and wild-derived inbred mice performing a graded treadmill test. Genome-wide association mapping was performed with an efficient mixed-model association (EMMA) algorithm to identify QTL. Exercise capacity was significantly different across strains. Run time varied by 2.7-fold between the highest running strain (C58/J) and the lowest running strain (A/J). These same strains showed a 16.5-fold difference in work. Significant associations were identified for exercise time on chromosomes 1, 2, 7, 11, and 13. The QTL interval on chromosome 2 (~168 Mb) contains one gene, Nfatc2, and overlaps with a suggestive QTL for training responsiveness in humans. These results provide phenotype data on the widest range of inbred strains tested thus far and indicate that genetic background significantly influences exercise capacity. Furthermore, the novel QTLs identified in the current study provide new targets for investigating the underlying mechanisms for variation in exercise capacity.

Hypothalamic expression of Peg3 gene is associated with maternal care differences between SM/J and LG/J mouse strains

Maternal care is essential in mammals, and variations in the environment provided by mothers may directly influence the viability of newborns and emotional behavior later in life. A previous study investigated genetic variations associated with maternal care in an intercross of LG/J and SM/J inbred mouse strains and identified two single-locus QTLs (quantitative trait loci). Here, we selected three candidate genes located within these QTLs intervals; Oxt on chromosome 2, and FosB and Peg3 on chromosome 7 and tested their association with maternal care. LG/J females showed impaired postpartum nest building and pup retrieval, a one-day delay in milk ejection, reduced exploratory activity, and higher anxiety-like behavior when compared to SM/J females. The nucleotide sequences of Oxt and FosB were similar between strains, as were their hypothalamic expression levels. Conversely, Peg3 nucleotide sequences showed four nonsynonymous replacement substitutions on LG/J dams, T11062G, G13744A, A13808G, and G13813A, and a 30 base pair (10 aa) in tandem repeat in the coding region with three copies in SM/J and five copies in LG/J. Maternal care impaired LG/J mothers express 37% lower Peg3 mRNA levels in the hypothalamus on the second postpartum day. We also found an association of the Peg3 repeat-variant and poor maternal care in F(2) heterozygote females derived from a LG/J × SM/J intercross. These results may suggest that the maternally imprinted Peg3 gene is responsible for the single-locus QTL on chromosome 7 that has been shown to influence maternal care in these strains. Furthermore, these data provide additional support for an epigenetic regulation of maternal behavior.

Weak genetic relationship between trabecular bone morphology and obesity in mice

Obesity, in addition to being associated with metabolic diseases, such as diabetes, has also been found to lower the risk of osteoporotic fractures. The relationship between obesity and bone trabecular structure is complex, involving responses to mechanical loading and the effects of adipocyte-derived hormones, both directly interacting with bone tissue and indirectly through central nervous system signaling. Here we examine the effects of sex, a high fat diet, and genetics on the trabecular density and structure of the lumbar and caudal vertebra and the proximal tibia along with body weight, fat pad weight, and serum leptin levels in a murine obesity model, the LGXSM recombinant inbred (RI) mouse strains. The sample included 481 mice from 16 RI strains. We found that vertebral trabecular density was higher in males while the females had higher tibial trabecular density. The high fat diet led to only slightly higher trabecular density in both sexes despite its extreme effects on obesity and serum leptin levels. Trait heritabilities are moderate to strong and genetic correlations among trabecular features are high. Most genetic variation contrasts strains with large numbers of thick, closely-spaced, highly interconnected, plate-like trabeculae with a high bone volume to total volume ratio against strains displaying small numbers of thin, widely-spaced, sparsely connected, rod-like trabeculae with a low bone volume to total volume ratio. Genetic correlations between trabecular and obesity-related traits were low and not statistically significant. We mapped trabecular properties to 20 genomic locations. Only one-quarter of these locations also had effects on obesity. In this population obesity has a relatively minor effect on trabecular bone morphology.

Heritability of articular cartilage regeneration and its association with ear wound healing in mice

OBJECTIVE

Emerging evidence suggests that genetic components contribute significantly to cartilage degeneration in osteoarthritis pathophysiology, but little information is available on the genetics of cartilage regeneration. Therefore, this study was undertaken to investigate cartilage regeneration in genetic murine models using common inbred strains and a set of recombinant inbred (RI) lines generated from LG/J (healer of ear wounds) and SM/J (nonhealer) inbred mouse strains.

METHODS

An acute full-thickness cartilage injury was introduced in the trochlear groove of 8-week-old mice (n=265) through microsurgery. Mouse knee joints were sagittally sectioned and stained with toluidine blue to evaluate regeneration. For the ear wound phenotype, a bilateral 2-mm through-and-through puncture was created in 6-week-old mice (n=229), and healing outcomes were measured after 30 days. Broad-sense heritability and genetic correlations were calculated for both phenotypes.

RESULTS

Time-course analysis of the RI mouse lines showed no significant regeneration until 16 weeks after surgery; at that time, the strains could be segregated into 3 categories: good, intermediate, and poor healers. Analysis of heritability (H2) showed that both cartilage regeneration (H2=26%; P=0.006) and ear wound closure (H2=53%; P<0.00001) were significantly heritable. The genetic correlations between the two healing phenotypes for common inbred mouse strains (r=0.92) and RI mouse lines (r=0.86) were found to be extremely high.

CONCLUSION

Our findings indicate that articular cartilage regeneration in mice is heritable, the differences between the mouse lines are due to genetic differences, and a strong genetic correlation between the two phenotypes exists, indicating that they plausibly share a common genetic basis. We therefore surmise that LG/J by SM/J intercross mice can be used to dissect the genetic basis of variation in cartilage regeneration.

Genetic architecture of nest building in mice LG/J × SM/J

Maternal care is critical to offspring growth and survival, which is greatly improved by building an effective nest. Some suggest that genetic variation and underlying genetic effects differ between fitness-related traits and other phenotypes. We investigated the genetic architecture of a fitness-related trait, nest building, in F₂ female mice intercrossed from inbred strains SM/J and LG/J using a QTL analysis for six related nest phenotypes (Presence and Structure pre- and postpartum, prepartum Material Used and postpartum Temperature). We found 15 direct-effect QTLs explaining from 4 to 13% of the phenotypic variation in nest building, mostly with non-additive effect. Epistatic analyses revealed 71 significant epistatic interactions which together explain from 28.4 to 75.5% of the variation, indicating an important role for epistasis in the adaptive process of nest building behavior in mice. Our results suggest a genetic architecture with small direct effects and a larger number of epistatic interactions as expected for fitness-related phenotypes.

Extensive epistasis for olfactory behaviour, sleep and waking activity in Drosophila melanogaster

Epistasis is an important feature of the genetic architecture of quantitative traits, but the dynamics of epistatic interactions in natural populations and the relationship between epistasis and pleiotropy remain poorly understood. Here, we studied the effects of epistatic modifiers that segregate in a wild-derived Drosophila melanogaster population on the mutational effects of P-element insertions in Semaphorin-5C (Sema-5c) and Calreticulin (Crc), pleiotropic genes that affect olfactory behaviour and startle behaviour and, in the case of Crc, sleep phenotypes. We introduced Canton-S B (CSB) third chromosomes with or without a P-element insertion at the Crc or Sema-5c locus in multiple wild-derived inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and assessed the effects of epistasis on the olfactory response to benzaldehyde and, for Crc, also on sleep. In each case, we found substantial epistasis and significant variation in the magnitude of epistasis. The predominant direction of epistatic effects was to suppress the mutant phenotype. These observations support a previous study on startle behaviour using the same D. melanogaster chromosome substitution lines, which concluded that suppressing epistasis may buffer the effects of new mutations. However, epistatic effects are not correlated among the different phenotypes. Thus, suppressing epistasis appears to be a pervasive general feature of natural populations to protect against the effects of new mutations, but different epistatic interactions modulate different phenotypes affected by mutations at the same pleiotropic gene.

Advances in comparative genetics: influence of genetics on obesity

Obesity has reached epidemic proportions and is recognised as a significant global health problem. Increased food intake and decreased physical activity are traditionally to blame for the development of obesity; however, many variables such as behaviour, diet, environment, social structures and genetics also contribute to this multifactorial disease. Complex interactions among these variables (for example, gene-environment, gene-diet and gene-gene) contribute not only to individual differences in the development of obesity, but also in treatment response. Mouse models have historically played valuable roles in understanding the genetics of traits related to energy balance and obesity. In the present review, we survey past use and examine new advances in mouse models designed to uncover the genetic architecture of obesity and its component traits. We discuss traditional models such as inbred strains and selectively bred lines and their contributions and shortcomings. We consider the evolution of mouse models into more informative resources such as outbred crosses and the Hybrid Mouse Diversity Panel, as well as novel next-generation approaches such as the Collaborative Cross. Moreover, the genetic architecture of voluntary exercise and the interactive relationship between host genetics and the gut microbiome are presented as novel phenotypes that augment studies using body weight and body fat percentage as endpoints. Understanding the intricate network of phenotypic, genotypic and environmental variables that predispose individuals to obesity will elucidate biological networks involved in the development of obesity. Knowledge obtained from advances in mouse models will inform human health and provide insight into inter-individual variability in the aetiology of obesity-related diseases.

Sex-, diet-, and cancer-dependent epistatic effects on complex traits in mice

The genetic basis of quantitative traits such as body weight and obesity is complex, with several hundred quantitative trait loci (QTLs) known to affect these and related traits in humans and mice. It also has become increasingly evident that the single-locus effects of these QTLs vary considerably depending on factors such as the sex of the individuals and their dietary environment, and we were interested to know whether this context-dependency also applies to two-locus epistatic effects of QTLs as well. We therefore conducted a genome scan to search for epistatic effects on 13 different weight and adiposity traits in an F₂ population of mice (created from an original intercross of the FVB strain with M16i, a polygenic obesity model) that were fed either a control or a high-fat diet and half of which harbored a transgene (PyMT) that caused the development of metastatic mammary cancer. We used a conventional interval mapping approach with SNPs to scan all 19 autosomes, and found extensive epistasis affecting all of these traits. More importantly, we also discovered that the majority of these epistatic effects exhibited significant interactions with sex, diet, and/or presence of PyMT. Analysis of these interactions showed that many of them appeared to involve QTLs previously identified as affecting these traits, but whose single-locus effects were variously modified by two-locus epistatic effects of other QTLs depending on the sex, diet, or PyMT environment. It was concluded that this context-dependency of epistatic effects is an important component of the genetic architecture of complex traits such as those contributing to weight and obesity.

Progeny-testing of full-sibs IBD in a SSC2 QTL region highlights epistatic interactions for fatness traits in pigs

Background

Many QTL have been detected in pigs, but very few of them have been fine-mapped up to the causal mutation. On SSC2, the IGF2-intron3-G3072A mutation has been described as the causative polymorphism for a QTL underlying muscle mass and backfat deposition, but further studies have demonstrated that at least one additional QTL should segregate downstream of this mutation. A marker-assisted backcrossing design was set up in order to confirm the segregation of this second locus, reduce its confidence interval and better understand its mode of segregation.

Results

Five recombinant full-sibs, with genotype G/G at the IGF2 mutation, were progeny-tested. Only two of them displayed significant QTL for fatness traits although four inherited the same paternal and maternal chromosomes, thus exhibiting the same haplotypic contrast in the QTL region. The hypothesis of an interaction with another region in the genome was proposed to explain these discrepancies and after a genome scan, four different regions were retained as potential interacting regions with the SSC2 QTL. A candidate interacting region on SSC13 was confirmed by the analysis of an F2 pedigree, and in the backcross pedigree one haplotype in this region was found to mask the SSC2 QTL effect.

Conclusions

Assuming the hypothesis of interactions with other chromosomal regions, the QTL could be unambiguously mapped to a 30 cM region delimited by recombination points. The marker-assisted backcrossing design was successfully used to confirm the segregation of a QTL on SSC2 and, because full-sibs that inherited the same alleles from their two parents were analysed, the detection of epistatic interactions could be performed between alleles and not between breeds as usually done with the traditional Line-Cross model. Additional analyses of other recombinant sires should provide more information to further improve the fine-mapping of this locus, and confirm or deny the interaction identified between chromosomes 2 and 13.

Fine-mapping alleles for body weight in LG/J × SM/J F₂ and F(34) advanced intercross lines

The present study measured variation in body weight using a combined analysis in an F(2) intercross and an F(34) advanced intercross line (AIL). Both crosses were derived from inbred LG/J and SM/J mice, which were selected for large and small body size prior to inbreeding. Body weight was measured at 62 (± 5) days of age. Using an integrated GWAS and forward model selection approach, we identified 11 significant QTLs that affected body weight on ten different chromosomes. With these results we developed a full model that explained over 18% of the phenotypic variance. The median 1.5-LOD support interval was 5.55 Mb, which is a significant improvement over most prior body weight QTLs. We identified nonsynonymous coding SNPs between LG/J and SM/J mice in order to further narrow the list of candidate genes. Three of the genes with nonsynonymous coding SNPs (Rad23b, Stk33, and Anks1b) have been associated with adiposity, waist circumference, and body mass index in human GWAS, thus providing evidence that these genes may underlie our QTLs. Our results demonstrate that a relatively small number of loci contribute significantly to the phenotypic variance in body weight, which is in marked contrast to the situation in humans. This difference is likely to be the result of strong selective pressure and the simplified genetic architecture, both of which are important advantages of our system.

Four out of eight genes in a mouse chromosome 7 congenic donor region are candidate obesity genes

We previously identified a region of mouse chromosome 7 that influences body fat mass in F2 littermates of congenic × background intercrosses. Current analyses revealed that alleles in the donor region of the subcongenic B6.C-D7Mit318 (318) promoted a twofold increase in adiposity in homozygous lines of 318 compared with background C57BL/6ByJ (B6By) mice. Parent-of-origin effects were discounted through cross-fostering studies and an F1 reciprocal cross. Mapping of the donor region revealed that it has a maximal size of 2.8 Mb (minimum 1.8 Mb) and contains a maximum of eight protein coding genes. Quantitative PCR in whole brain, liver, and gonadal white adipose tissue (GWAT) revealed differential expression between genotypes for three genes in females and two genes in males. Alpha-2,8-sialyltransferase 8B (St8sia2) showed reduced 318 mRNA levels in brain for females and males and in GWAT for females only. Both sexes of 318 mice had reduced Repulsive guidance molecule-a (Rgma) expression in GWAT. In brain, Family with sequence similarity 174 member b (Fam174b) had increased expression in 318 females, whereas Chromodomain helicase DNA binding protein 2 (Chd2-2) had reduced expression in 318 males. No donor region genes were differentially expressed in liver. Sequence analysis of coding exons for all genes in the 318 donor region revealed only one single nucleotide polymorphism that produced a nonsynonymous missense mutation, Gln7Pro, in Fam174b. Our findings highlight the difficulty of using expression and sequence to identify quantitative trait genes underlying obesity even in small genomic regions.

The effect of a population bottleneck on the evolution of genetic variance/covariance structure

It is well known that standard population genetic theory predicts decreased additive genetic variance (V(a) ) following a population bottleneck and that theoretical models including interallelic and intergenic interactions indicate such loss may be avoided. However, few empirical data from multicellular model systems are available, especially regarding variance/covariance (V/CV) relationships. Here, we compare the V/CV structure of seventeen traits related to body size and composition between control (60 mating pairs/generation) and bottlenecked (2 mating pairs/generation; average F = 0.39) strains of mice. Although results for individual traits vary considerably, multivariate analysis indicates that V(a) in the bottlenecked populations is greater than expected. Traits with patterns and amounts of epistasis predictive of enhanced V(a) also show the largest deviations from additive expectations. Finally, the correlation structure of weekly weights is not significantly different between control and experimental lines but correlations between necropsy traits do differ, especially those involving the heart, kidney and tail length.

Dark matter: are mice the solution to missing heritability?

Genome-wide association studies (GWAS) in humans have identified hundreds of single nucleotide polymorphisms associated with complex traits, yet for most traits studied, the sum total of all these identified variants fail to explain a significant portion of the heritable variation. Reasons for this “missing heritability” are thought to include the existence of rare causative variants not captured by current genotyping arrays, structural variants that go undetected by existing technology, insufficient power to identify multi-gene interactions, small sample sizes, and the influence of environmental and epigenetic effects. As genotyping technologies have evolved it has become inexpensive and relatively straightforward to perform GWAS in mice. Mice offer a powerful tool for elucidating the genetic architecture of behavioral and physiological traits, and are complementary to human studies. Unlike F₂ crosses of inbred strains, advanced intercross lines, heterogeneous stocks, outbred, and wild-caught mice have more rapid breakdown of linkage disequilibrium which allow for increasingly high resolution mapping. Because some of these populations are created using a small number of founder chromosomes they are not expected to harbor rare alleles. We discuss the differences between these mouse populations and examine their potential to overcome some of the pitfalls that have plagued human GWAS studies.

Epistatic interactions of genes influence within-individual variation of physical activity traits in mice

A number of quantitative trait loci (QTLs) recently have been discovered that affect various activity traits in mice, but their collective impact does not appear to explain the consistently moderate to high heritabilities for these traits. We previously suggested interactions of genes, or epistasis, might account for additional genetic variability of activity, and tested this for the average distance, duration and speed run by mice during a 3 week period. We found abundant evidence for epistasis affecting these traits, although, recognized that epistatic effects may well vary within individuals over time. We therefore conducted a full genome scan for epistatic interactions affecting these traits in each of seven three-day intervals. Our intent was to assess the extent and trends in epistasis affecting these traits in each of the intervals. We discovered a number of epistatic interactions of QTLs that influenced the activity traits in the mice, the majority of which were not previously found and appeared to affect the activity traits (especially distance and speed) primarily in the early or in the late age intervals. The overall impact of epistasis was considerable, its contribution to the total phenotypic variance varying from an average of 22-35% in the three traits across all age intervals. It was concluded that epistasis is more important than single-locus effects of genes on activity traits at specific ages and it is therefore an essential component of the genetic architecture of physical activity.

Architecture of energy balance traits in emerging lines of the Collaborative Cross

The potential utility of the Collaborative Cross (CC) mouse resource was evaluated to better understand complex traits related to energy balance. A primary focus was to examine if genetic diversity in emerging CC lines (pre-CC) would translate into equivalent phenotypic diversity. Second, we mapped quantitative trait loci (QTL) for 15 metabolism- and exercise-related phenotypes in this population. We evaluated metabolic and voluntary exercise traits in 176 pre-CC lines, revealing phenotypic variation often exceeding that seen across the eight founder strains from which the pre-CC was derived. Many phenotypic correlations existing within the founder strains were no longer significant in the pre-CC population, potentially representing reduced linkage disequilibrium (LD) of regions harboring multiple genes with effects on energy balance or disruption of genetic structure of extant inbred strains with substantial shared ancestry. QTL mapping revealed five significant and eight suggestive QTL for body weight (Chr 4, 7.54 Mb; CI 3.32-10.34 Mb; Bwq14), body composition, wheel running (Chr 16, 33.2 Mb; CI 32.5-38.3 Mb), body weight change in response to exercise (1: Chr 6, 77.7Mb; CI 72.2-83.4 Mb and 2: Chr 6, 42.8 Mb; CI 39.4-48.1 Mb), and food intake during exercise (Chr 12, 85.1 Mb; CI 82.9-89.0 Mb). Some QTL overlapped with previously mapped QTL for similar traits, whereas other QTL appear to represent novel loci. These results suggest that the CC will be a powerful, high-precision tool for examining the genetic architecture of complex traits such as those involved in regulation of energy balance.

Evolution of pleiotropy: epistatic interaction pattern supports a mechanistic model underlying variation in genotype-phenotype map

The genotype-phenotype (GP) map consists of developmental and physiological mechanisms mapping genetic onto phenotypic variation. It determines the distribution of heritable phenotypic variance on which selection can act. Comparative studies of morphology as well as of gene regulatory networks show that the GP map itself evolves, yet little is known about the actual evolutionary mechanisms involved. The study of such mechanisms requires exploring the variation in GP maps at the population level, which presently is easier to quantify by statistical genetic methods rather than by regulatory network structures. We focus on the evolution of pleiotropy, a major structural aspect of the GP map. Pleiotropic genes affect multiple traits and underlie genetic covariance between traits, often causing evolutionary constraints. Previous quantitative genetic studies have demonstrated population-level variation in pleiotropy in the form of loci, at which genotypes differ in the genetic covariation between traits. This variation can potentially fuel evolution of the GP map under selection and/or drift. Here, we propose a developmental mechanism underlying population genetic variation in covariance and test its predictions. Specifically, the mechanism predicts that the loci identified as responsible for genetic variation in pleiotropy are involved in trait-specific epistatic interactions. We test this prediction for loci affecting allometric relationships between traits in an advanced intercross between inbred mouse strains. The results consistently support the prediction. We further find a high degree of sign epistasis in these interactions, which we interpret as an indication of adaptive gene complexes within the diverged parental lines.

Diet-dependent genetic and genomic imprinting effects on obesity in mice

Although the current obesity epidemic is of environmental origin, there is substantial genetic variation in individual response to an obesogenic environment. In this study, we perform a genome-wide scan for quantitative trait loci (QTLs) affecting obesity per se, or an obese response to a high-fat diet in mice from the LG/J by SM/J Advanced Intercross (AI) Line (Wustl:LG,SM-G16). A total of 1,002 animals from 78 F₁₆ full sibships were weaned at 3 weeks of age and half of each litter placed on high- and low-fat diets. Animals remained on the diet until 20 weeks of age when they were necropsied and the weights of the reproductive, kidney, mesenteric, and inguinal fat depots were recorded. Effects on these phenotypes, along with total fat depot weight and carcass weight at necropsy, were mapped across the genome using 1,402 autosomal single-nucleotide polymorphism (SNP) markers. Haplotypes were reconstructed and additive, dominance, and imprinting genotype scores were derived every 1 cM along the F₁₆ map. Analysis was performed using a mixed model with additive, dominance, and imprinting genotype scores, their interactions with sex, diet, and with sex-by-diet as fixed effects and with family and its interaction with sex, diet, and sex-by-diet as random effects. We discovered 95 trait-specific QTLs mapping to 40 locations. Most QTLs had additive effects with dominance and imprinting effects occurring at two-thirds of the loci. Nearly every locus interacted with sex and/or diet in important ways demonstrating that gene effects are primarily context dependent, changing depending on sex and/or diet.

Mapping the epistatic network underlying murine reproductive fatpad variation

Genome-wide mapping analyses are now commonplace in many species and several networks of interacting loci have been reported. However, relatively few details regarding epistatic interactions and their contribution to complex trait variation in multicellular organisms are available and the identification of positional candidate loci for epistatic QTL (epiQTL) is hampered, especially in mammals, by the limited genetic resolution inherent in most study designs. Here we further investigate the genetic architecture of reproductive fatpad weight in mice using the F(10) generation of the LG,SM advanced intercross (AI) line. We apply multiple mapping techniques including a single-locus model, locus-specific composite interval mapping (CIM), and tests for multiple QTL per chromosome to the 12 chromosomes known to harbor single-locus QTL (slQTL) affecting obesity in this cross. We also perform a genome-wide scan for pairwise epistasis. Using this combination of approaches we detect 199 peaks spread over all 19 autosomes, which potentially contribute to trait variation including all eight original F(2) loci (Adip1-8), novel slQTL peaks on chromosomes 7 and 9, and several novel epistatic loci. Extensive epistasis is confirmed involving both slQTL confidence intervals (C.I.) as well as regions that show no significant additive or dominance effects. These results provide important new insights into mapping complex genetic architectures and the role of epistasis in complex trait variation.

A wild derived quantitative trait locus on mouse chromosome 2 prevents obesity

Background

The genetic architecture of multifactorial traits such as obesity has been poorly understood. Quantitative trait locus (QTL) analysis is widely used to localize loci affecting multifactorial traits on chromosomal regions. However, large confidence intervals and small phenotypic effects of identified QTLs and closely linked loci are impeding the identification of causative genes that underlie the QTLs. Here we developed five subcongenic mouse strains with overlapping and non-overlapping wild-derived genomic regions from an F2 intercross of a previously developed congenic strain, B6.Cg-Pbwg1, and its genetic background strain, C57BL/6J (B6). The subcongenic strains developed were phenotyped on low-fat standard chow and a high-fat diet to fine-map a previously identified obesity QTL. Microarray analysis was performed with Affymetrix GeneChips to search for candidate genes of the QTL.

Results

The obesity QTL was physically mapped to an 8.8-Mb region of mouse chromosome 2. The wild-derived allele significantly decreased white fat pad weight, body weight and serum levels of glucose and triglyceride. It was also resistant to the high-fat diet. Among 29 genes residing within the 8.8-Mb region, Gpd2, Upp2, Acvr1c, March7 and Rbms1 showed great differential expression in livers and/or gonadal fat pads between B6.Cg-Pbwg1 and B6 mice.

Conclusions

The wild-derived QTL allele prevented obesity in both mice fed a low-fat standard diet and mice fed a high-fat diet. This finding will pave the way for identification of causative genes for obesity. A further understanding of this unique QTL effect at genetic and molecular levels may lead to the discovery of new biological and pathologic pathways associated with obesity.

Epistasis contributes to the genetic buffering of plasma HDL cholesterol in mice

Stressful environmental factors, such as a high-fat diet, can induce responses in the expression of genes that act to maintain physiological homeostasis. We observed variation in plasma concentrations of high-density lipoprotein (HDL) cholesterol across inbred mouse strains in response to high dietary fat intake. Several strains, including C57BL/6J, have stable levels of plasma HDL independent of diet, whereas other strains, including DBA2/J, show marked changes in plasma HDL. To explore this phenomenon further, we used publicly available data from a C57BL/6J × DBA/2J intercross to identify genetic factors that associate with HDL under high-fat diet conditions. Our analysis identified an epistatic interaction that plays a role in the buffering of HDL levels in C57BL/6J mice, and we have identified Arl4d as a candidate gene that mediates this effect. Structural modeling further elucidates the interaction of genetic factors that contribute to the robustness of HDL in response to high-fat diet in the C57BL/6J strain.

Fine-mapping of obesity-related quantitative trait loci in an F9/10 advanced intercross line

Obesity develops in response to a combination of environmental effects and multiple genes of small effect. Although there has been significant progress in characterizing genes in many pathways contributing to metabolic disease, knowledge about the relationships of these genes to each other and their joint effects upon obesity lags behind. The LG,SM advanced intercross line (AIL) model of obesity has been used to characterize over 70 loci involved in fatpad weight, body weight, and organ weights. Each of these quantitative trait loci (QTLs) encompasses large regions of the genome and require fine-mapping to isolate causative sequence changes and possible mechanisms of action as indicated by the genetic architecture. In this study we fine-map QTLs first identified in the F(2) and F(2/3) populations in the combined F(9/10) advanced intercross generations. We observed significantly narrowed QTL confidence regions, identified many single QTL that resolve into multiple QTL peaks, and identified new QTLs that may have been previously masked due to opposite gene effects at closely linked loci. We also present further characterization of the pleiotropic and epistatic interactions underlying these obesity-related traits.

Directionality of epistasis in a murine intercross population

Directional epistasis describes a situation in which epistasis consistently increases or decreases the effect of allele substitutions, thereby affecting the amount of additive genetic variance available for selection in a given direction. This study applies a recent parameterization of directionality of epistasis to empirical data. Data stems from a QTL mapping study on an intercross between inbred mouse (Mus musculus) strains LG/J and SM/J, originally selected for large and small body mass, respectively. Results show a negative average directionality of epistasis for body-composition traits, predicting a reduction in additive allelic effects and in the response to selection for increased size. Focusing on average modification of additive effect of single loci, we find a more complex picture, whereby the effects of some loci are enhanced consistently across backgrounds, while effects of other loci are decreased, potentially contributing to either enhancement or reduction of allelic effects when selection acts at single loci. We demonstrate and discuss how the interpretation of the overall measurement of directionality depends on the complexity of the genotype-phenotype map. The measure of directionality changes with the power of scale in a predictable way; however, its expected effect with respect to the modification of additive genetic effects remains constant.

Genotype-phenotype associations in understanding the role of corticosteroid-binding globulin in health and disease animal models

Corticosteroid-binding globulin (CBG) is a plasma glycoprotein discovered more than 60 years ago for its high-affinity for glucocorticoids. Although its molecular structure and its biochemical properties have been described, its various biological roles and its importance are not yet fully understood. This review focuses first on studies that have used no-hypothesis-driven genetic approaches in animal models to reveal the higher than expected importance of CBG in particular in glucocorticoid stress responses. Then the dissection of some CBG physiological roles in an animal model of genetic CBG deficiency is reported. Finally, studies on the role of CBG genetic variability in human obesity traits are reviewed and discussed.