Identification of quantitative trait loci affecting body composition in a mouse intercross

A GWA study reveals genetic loci for body conformation traits in Chinese Laiwu pigs and its implications for human BMI

Pigs share numerous physiological and phenotypic similarities with human and thus have been considered as a good model in nonrodent mammals for the study of genetic basis of human obesity. Researches on candidate genes for obesity traits have successfully identified some common genes between humans and pigs. However, few studies have assessed how many similarities exist between the genetic architecture of obesity in pigs and humans by large-scale comparative genomics. Here, we performed a genome-wide association study (GWAS) using the porcine 60 K SNP Beadchip for BMI and other four conformation traits at three different ages in a Chinese Laiwu pig population, which shows a large variability in fat deposition. In total, 35 SNPs were found to be significant at Bonferroni-corrected 5 % chromosome-wise level (P = 2.13 × 10^-5) and 88 SNPs had suggestive (P < 10^-4) association with the conformation traits. Some SNPs showed age-dependent association. Intriguingly, out of 32 regions associated with BMI in pigs, 18 were homologous with the loci for BMI in humans. Furthermore, five closest genes to GWAS peaks including HIF1AN, SMYD3, COX10, SLMAP, and GBE1 have been already associated with BMI in humans, which makes them very promising candidates for these QTLs. The result of GO analysis provided strong support to the fact that mitochondria and synapse play important roles in obesity susceptibility, which is consistent with previous findings on human obesity, and it also implicated new gene sets related to chromatin modification and Ig-like C2-type 5 domain. Therefore, these results not only provide new insights into the genetic architecture of BMI in pigs but also highlight that humans and pigs share the significant overlap of obesity-related genes.

Evaluation of growth patterns and body composition in C57Bl/6J mice using dual energy X-ray absorptiometry

The normal growth pattern of female C57BL/6J mice, from 5 to 30 weeks of age, has been investigated in a longitudinal study. Weight, body surface area (BS), and body mass index (BMI) were evaluated in forty mice. Lean mass and fat mass, bone mineral content (BMC), and bone mineral density (BMD) were monitored by dual energy X-ray absorptiometry (DEXA). Weight and BS increased linearly (16.15 ± 0.64-27.64 ± 1.42 g; 51.13 ± 0.74-79.57 ± 2.15 cm(2), P < 0.01), more markedly from 5 to 9 weeks of age (P < 0.001). BMD showed a peak at 17 weeks (0.0548 ± 0.0011 g/cm(2) ∗ m, P < 0.01). Lean mass showed an evident gain at 9 (15.8 ± 0.8 g, P < 0.001) and 25 weeks (20.5 ± 0.3 g, P < 0.01), like fat mass from 13 to 17 weeks (2.0 ± 0.4-3.6 ± 0.7 g, P < 0.01). BMI and lean mass index (LMI) reached the highest value at 21 weeks (3.57 ± 0.02-0.284 ± 0.010 g/cm(2), resp.), like fat mass index (FMI) at 17 weeks (0.057 ± 0.009 g/cm(2)) (P < 0.01). BMI, weight, and BS showed a moderate positive correlation (0.45-0.85) with lean mass from 5 to 21 weeks. Mixed linear models provided a good prediction for lean mass, fat mass, and BMD. This study may represent a baseline reference for a future comparison of wild-type C57BL/6J mice with models of altered growth.

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

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

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.

A proximal genomic region of mouse chromosome 10 contains quantitative trait loci affecting fatness

Body weight and fatness are quantitative traits of agricultural and medical importance. In previous genome-wide quantitative trait locus (QTL) analyses, two QTLs for body weight and weight gain at an early postnatal growth period were discovered on mouse chromosome 10 from a gene pool of wild subspecies mice, Mus musculus castaneus. In this study, we developed a congenic strain with an approximately 63-Mb wild-derived genomic region on which the two growth QTLs could be located, by recurrent backcrossing to the common inbred strain C57BL/6J. We compared body weights at 1-10 weeks of age, body weight gains at 1-3, 3-6 and 6-10 weeks, internal organ weights and body lengths between the congenic strain developed and C57BL/6J. Unfortunately, no effects of the two growth QTLs on body weights and weight gains were confirmed. However, at least two new QTLs affecting fatness traits were discovered within the introgressed congenic region. The wild-derived allele at one QTL increased body mass index, whereas at another one it decreased white fat pad weight and adiposity index. Thus, the congenic mouse strain developed here is a useful model animal for understanding the genetic and molecular basis of fat deposition in livestock as well as humans.

Physiological homology between Drosophila melanogaster and vertebrate cardiovascular systems

The physiology of the Drosophila melanogaster cardiovascular system remains poorly characterized compared with its vertebrate counterparts. Basic measures of physiological performance remain unknown. It also is unclear whether subtle physiological defects observed in the human cardiovascular system can be reproduced in D. melanogaster. Here we characterize the cardiovascular physiology of D. melanogaster in its pre-pupal stage by using high-speed dye angiography and optical coherence tomography. The heart has vigorous pulsatile contractions that drive intracardiac, aortic and extracellular-extravascular hemolymph flow. Several physiological measures, including weight-adjusted cardiac output, body-length-adjusted aortic velocities and intracardiac shear forces, are similar to those in the closed vertebrate cardiovascular systems, including that of humans. Extracellular-extravascular flow in the pre-pupal D. melanogaster circulation drives convection-limited fluid transport. To demonstrate homology in heart dysfunction, we showed that, at the pre-pupal stage, a troponin I mutant, held-up2 (hdp2), has impaired systolic and diastolic heart wall velocities. Impaired heart wall velocities occur in the context of a non-dilated phenotype with a mildly depressed fractional shortening. We additionally derive receiver operating characteristic curves showing that heart wall velocity is a potentially powerful discriminator of systolic heart dysfunction. Our results demonstrate physiological homology and support the use of D. melanogaster as an animal model of complex cardiovascular disease.

Mouse BMD quantitative trait loci show improved concordance with human genome-wide association loci when recalculated on a new, common mouse genetic map

Bone mineral density (BMD) is a heritable trait, and in mice, over 100 quantitative trait loci (QTLs) have been reported, but candidate genes have been identified for only a small percentage. Persistent errors in the mouse genetic map have negatively affected QTL localization, spurring the development of a new, corrected map. In this study, QTLs for BMD were remapped in 11 archival mouse data sets using this new genetic map. Since these QTLs all were mapped in a comparable way, direct comparisons of QTLs for concordance would be valid. We then compared human genome-wide association study (GWAS) BMD loci with the mouse QTLs. We found that 26 of the 28 human GWAS loci examined were located within the confidence interval of a mouse QTL. Furthermore, 14 of the GWAS loci mapped to within 3 cM of a mouse QTL peak. Lastly, we demonstrated that these newly remapped mouse QTLs can substantiate a candidate gene for a human GWAS locus, for which the peak single-nucleotide polymorphism (SNP) fell in an intergenic region. Specifically, we suggest that MEF2C (human chromosome 5, mouse chromosome 13) should be considered a candidate gene for the genetic regulation of BMD. In conclusion, use of the new mouse genetic map has improved the localization of mouse BMD QTLs, and these remapped QTLs show high concordance with human GWAS loci. We believe that this is an opportune time for a renewed effort by the genetics community to identify the causal variants regulating BMD using a synergistic mouse-human approach.

Bivariate genome-wide linkage analysis of femoral bone traits and leg lean mass: Framingham study

The risk of osteoporotic fracture is a function of both applied muscle mass and bone tissue distribution. Leg lean mass (LLM) and femoral bone geometry are both known to have substantial genetic components. Therefore, we estimated shared heritability (h(2)) and performed linkage analysis to identify chromosomal regions governing both LLM and bone geometry. A genome-wide scan (using 636 microsatellite markers) for linkage analyses was performed on 1346 adults from 327 extended families of the Framingham study. DXA measures were LLM, femoral neck length, neck-shaft angle (NSA), subperiosteal width, cross-sectional area (CSA), and section modulus (Z) at the femoral narrow neck and shaft (S) regions. Variance component linkage analysis was performed on normalized residuals (adjusted for age, height, BMI, and estrogen status in women). The results indicated substantial h(2) for LLM (0.42 +/- 0.07) that was comparable to bone geometry traits. Phenotypic correlations between LLM and bone geometry phenotypes ranged from 0.033 with NSA (p > 0.05) to 0.251 with S_Z (p < 0.001); genetic correlations ranged from 0.087 (NSA, p > 0.05) to 0.454 (S_Z, p < 0.001). Univariate linkage analysis of covariate-adjusted LLM identified no chromosomal regions with LOD scores >or=2.0; however, bivariate analysis identified two loci with LOD scores >3.0, shared by LLM with S_CSA on chromosome 12p12.3-12p13.2, and with NSA, on 14q21.3-22.1. In conclusion, we identified chromosomal regions potentially linked to both LLM and femoral bone geometry. Identification and subsequent characterization of these shared loci may further elucidate the genetic contributions to both osteoporosis and sarcopenia.

Genetic influences on growth and body composition in mice: multilocus interactions

BACKGROUND

The genetic architecture of body weight and body composition is complex because these traits are normally influenced by multiple genes and their interactions, even after controlling for the environment. Bayesian methodology provides an efficient way of estimating these interactions.

SUBJECTS AND MEASUREMENTS

We used Bayesian model selection techniques to simultaneously estimate the main effects, epistasis and gene-sex interactions on age-related body weight (at 3, 6 and 10 weeks, denoted as WT3wk, WT6wk and WT10wk) and body composition (organ weights and fat-related traits) in an F(2) sample obtained from a cross between high-growth (M16i) mice and low-growth (L6) mice.

RESULTS

We observed epistatic and main-effect quantitative trait loci (QTL) that controlled both body weight and body composition. Epistatic effects were generally more significant for WT6wk than WT10wk. Chromosomes 5 and 13 interacted strongly to control body weight at 3 weeks. A pleiotropic QTL on chromosome 2 was associated with body weight and some body composition phenotypes. Testis weight was regulated by a QTL on chromosome 13 with a significantly large main effect (2log(e)BF approximately 15).

CONCLUSION

By analyzing epistatic interactions, we detected QTL not found in a previous analysis of this mouse population. Hence, the detection of gene-gene interactions may provide new information about the genetic architecture of complex obesity-related traits and may lead to the detection of additional obesity genes.

Hitchhiking mapping reveals a candidate genomic region for natural selection in three-spined stickleback chromosome VIII

Identification of genes and genomic regions under directional natural selection has become one of the major goals in evolutionary genetics, but relatively little work to this end has been done by applying hitchhiking mapping to wild populations. Hitchhiking mapping starts from a genome scan using a randomly spaced set of molecular markers followed by a fine-scale analysis in the flanking regions of the candidate regions under selection. We used the hitchhiking mapping approach to narrow down a selective sweep in the genomic region flanking a candidate locus (Stn90) in chromosome VIII in the three-spined stickleback (Gasterosteus aculeatus). Twenty-four microsatellite markers were screened in an approximately 800-kb region around the candidate locus in three marine and four freshwater populations. The patterns of genetic diversity and differentiation in the candidate region were compared to those of a putatively neutral set of markers. The Bayesian FST-test indicated an elevated genetic differentiation, deviating significantly from neutral expectations, at a continuous region of approximately 20 kb upstream from the candidate locus. Furthermore, a method developed for an array of microsatellite markers rejected neutrality in a region of approximately 90 kb flanking the candidate locus supporting the selective sweep hypothesis. Likewise, the genomewide pattern of genetic diversity differed from the candidate region in a bottleneck analysis suggesting that selection, rather than demography, explains the reduced genetic diversity at the candidate interval. The neutrality tests suggest that the selective sweep had occurred mainly in the Lake Pulmanki population, but the results from bottleneck analyses indicate that selection might have operated in other populations as well. These results suggest that the narrow interval around locus Stn90 has likely been under directional selection, but the region contains several predicted genes, each of which can be the actual targets of selection. Understanding of the functional significance of this genomic region in an ecological context will require a more detailed sequence analysis.

Comparison between surrogate indexes of insulin sensitivity and resistance and hyperinsulinemic euglycemic clamp estimates in mice

Insulin resistance contributes to the pathophysiology of diabetes, obesity, and their cardiovascular complications. Mouse models of these human diseases are useful for gaining insight into pathophysiological mechanisms. The reference standard for measuring insulin sensitivity in both humans and animals is the euglycemic glucose clamp. Many studies have compared surrogate indexes of insulin sensitivity and resistance with glucose clamp estimates in humans. However, regulation of metabolic physiology in humans and rodents differs and comparisons between surrogate indexes and the glucose clamp have not been directly evaluated in rodents previously. Therefore, in the present study, we compared glucose clamp-derived measures of insulin sensitivity (GIR and SI(Clamp)) with surrogate indexes, including quantitative insulin-sensitivity check index (QUICKI), homeostasis model assessment (HOMA), 1/HOMA, log(HOMA), and 1/fasting insulin, using data from 87 mice with a wide range of insulin sensitivities. We evaluated simple linear correlations and performed calibration model analyses to evaluate the predictive accuracy of each surrogate. All surrogate indexes tested were modestly correlated with both GIR and SI(Clamp). However, a stronger correlation between body weight per se and both GIR and SI(Clamp) was noted. Calibration analyses of surrogate indexes adjusted for body weight demonstrated improved predictive accuracy for GIR [e.g., R = 0.68, for QUICKI and log(HOMA)]. We conclude that linear correlations of surrogate indexes with clamp data and predictive accuracy of surrogate indexes in mice are not as substantial as in humans. This may reflect intrinsic differences between human and rodent physiology as well as increased technical difficulties in performing glucose clamps in mice.

Genetic loci affecting body weight and fatness in a C57BL/6J x PWK/PhJ mouse intercross

To determine the genetic variation that contributes to body composition in the mouse, we interbred a wild-derived strain (PWK/PhJ; PWK) with a common laboratory strain (C57BL/6J; B6). The parental, F(1), and F(2) mice were phenotyped at 18 weeks old for body weight and composition using dual-energy X-ray absorptiometry (DEXA). A total of 479 (244 male and 235 female) F(2) mice were genotyped for 117 polymorphic markers spanning the autosomes. Twenty-eight suggestive or significant linkages for four traits (body weight, adjusted lean and fat weight, and percent fat) were detected. Of these, three QTLs were novel: one on the proximal portion of Chr 5 for body weight (Bwq8; LOD = 4.7), one on Chr 3 for lean weight (Bwtq13; LOD = 3.6), and one on Chr 11 for percent fat (Adip19; LOD = 5.8). The remaining QTLs overlapped previously identified linkages, e.g., Adip5 on Chr 9. One QTL was sex-specific (present in males only) and seven were sex-biased (more prominent in one sex than the other). Most alleles that increased body weight were contributed by the B6 strain, and most alleles that increased percent fat were contributed by the PWK strain. Eight pairs of interacting loci were identified, none of which exactly overlapped the main-effect QTLs. Many of the QTLs found in the B6 x PWK cross map to the location of previously reported linkages, suggesting that some QTLs are common to many strains (consensus QTLs), but three new QTLs appear to be particular to the PWK strain. The location and type of QTLs detected in this new cross will assist in future efforts to identify the genetic variation that determines the ratio of lean to fat weight as well as body size in mice.

Chemical mutagenesis induced two high bone density mouse mutants map to a concordant distal chromosome 4 locus

Phenotype-driven mutagenesis approach in the mouse holds much promise as a method for revealing gene function. Earlier, we have described an N-ethyl-N-nitrosourea (ENU) mutagenesis screen to create genome-wide dominant mutations in the mouse model. Using this approach, we describe identification of two high bone density mutants in C57BL/6J (B6) background. The mutants, named as 12184 and 12137, have been bred more than five generations with wild-type B6 mice, each producing >200 backcross progeny. The average total body areal bone mineral density (aBMD) was 13-17% higher in backcrossed progeny from both mutant lines between 6 and 10 weeks of age, as compared to wild-type (WT) B6 mice (n=60-107). At 3 weeks of age the aBMD of mutant progeny was not significantly affected as compared to WT B6 mice. Data from 10- and 16-week old progeny show that increased aBMD was mainly related to a 14-20% higher bone mineral content, whereas bone size was marginally increased. In addition, the average volumetric BMD (vBMD) was 5-15% higher at the midshaft tibia or femur, as compared to WT mice. Histomorphometric analysis revealed that bone resorption was 23-34% reduced in both mutant mice. Consistent with histomorphometry data, the mRNA expression of genes that regulate osteoclast differentiation and survival were altered in the 12137 mutant mice. To determine the chromosomal location of the ENU mutation, we intercrossed both mutant lines with C3H/HeJ (C3H) mice to generate B6C3H F2 mice (n=164 for line 12137 and n=137 F2 for line 12184). Interval mapping using 60 microsatellite markers and aBMD phenotype revealed only one significant or suggestive linkage on chromosome 4. Since body weight was significantly higher in mutant lines, we also used body weight as additive and interactive covariate for interval mapping; both analyses showed higher LOD scores for both 12137 and 12184 mutants without affecting the chromosomal location. The large phenotype in the mutant mice compared to generally observed QTL effects (<5%) would increase the probability of identifying the mutant gene.