Identification of quantitative trait loci affecting murine long bone length in a two-generation intercross of LG/J and SM/J Mice

Genetic architecture of trait variance in craniofacial morphology

The genetic architecture of trait variance has long been of interest in genetics and evolution. One of the earliest attempts to understand this architecture was presented in Lerner's Genetic Homeostasis (1954). Lerner proposed that heterozygotes should be better able to tolerate environmental perturbations because of functional differences between the alleles at a given locus, with each allele optimal for slightly different environments. This greater robustness to environmental variance, he argued, would result in smaller trait variance for heterozygotes. The evidence for Lerner's hypothesis has been inconclusive. To address this question using modern genomic methods, we mapped loci associated with differences in trait variance (vQTL) on 1,101 individuals from the F34 of an advanced intercross between LG/J and SM/J mice. We also mapped epistatic interactions for these vQTL in order to understand the influence of epistasis for the architecture of trait variance. We did not find evidence supporting Lerner's hypothesis, that heterozygotes tend to have smaller trait variances than homozygotes. We further show that the effects of most mapped loci on trait variance are produced by epistasis affecting trait means and that those epistatic effects account for about a half of the differences in genotypic-specific trait variances. Finally, we propose a model where the different interactions between the additive and dominance effects of the vQTL and their epistatic partners can explain Lerner's original observations but can also be extended to include other conditions where heterozygotes are not the least variable genotype.

Independent evolution toward larger body size in the distinctive Faroe Island mice

Most phenotypic traits in nature involve the collective action of many genes. Traits that evolve repeatedly are particularly useful for understanding how selection may act on changing trait values. In mice, large body size has evolved repeatedly on islands and under artificial selection in the laboratory. Identifying the loci and genes involved in this process may shed light on the evolution of complex, polygenic traits. Here, we have mapped the genetic basis of body size variation by making a genetic cross between mice from the Faroe Islands, which are among the largest and most distinctive natural populations of mice in the world, and a laboratory mouse strain selected for small body size, SM/J. Using this F2 intercross of 841 animals, we have identified 111 loci controlling various aspects of body size, weight and growth hormone levels. By comparing against other studies, including the use of a joint meta-analysis, we found that the loci involved in the evolution of large size in the Faroese mice were largely independent from those of a different island population or other laboratory strains. We hypothesize that colonization bottleneck, historical hybridization, or the redundancy between multiple loci have resulted in the Faroese mice achieving an outwardly similar phenotype through a distinct evolutionary path.

Endochondral ossification and the evolution of limb proportions

Mammals have remarkably diverse limb proportions hypothesized to have evolved adaptively in the context of locomotion and other behaviors. Mechanistically, evolutionary diversity in limb proportions is the result of differential limb bone growth. Longitudinal limb bone growth is driven by the process of endochondral ossification, under the control of the growth plates. In growth plates, chondrocytes undergo a tightly orchestrated life cycle of proliferation, matrix production, hypertrophy, and cell death/transdifferentiation. This life cycle is highly conserved, both among the long bones of an individual, and among homologous bones of distantly related taxa, leading to a finite number of complementary cell mechanisms that can generate heritable phenotype variation in limb bone size and shape. The most important of these mechanisms are chondrocyte population size in chondrogenesis and in individual growth plates, proliferation rates, and hypertrophic chondrocyte size. Comparative evidence in mammals and birds suggests the existence of developmental biases that favor evolutionary changes in some of these cellular mechanisms over others in driving limb allometry. Specifically, chondrocyte population size may evolve more readily in response to selection than hypertrophic chondrocyte size, and extreme hypertrophy may be a rarer evolutionary phenomenon associated with highly specialized modes of locomotion in mammals (e.g., powered flight, ricochetal bipedal hopping). Physical and physiological constraints at multiple levels of biological organization may also have influenced the cell developmental mechanisms that have evolved to produce the highly diverse limb proportions in extant mammals. This article is categorized under: Establishment of Spatial and Temporal Patterns > Regulation of Size, Proportion, and Timing Comparative Development and Evolution > Regulation of Organ Diversity Comparative Development and Evolution > Organ System Comparisons Between Species.

Artificial selection sheds light on developmental mechanisms of limb elongation

Species diversity in limb lengths and proportions is thought to have evolved adaptively in the context of locomotor and habitat specialization, but the heritable cellular processes that drove this evolution within species are poorly understood. In this study, we take a novel "micro-evo-devo" approach, using artificial selection on relative limb length to amplify phenotypic variation in a population of mice, known as Longshanks, to examine the cellular mechanisms of postnatal limb development that contribute to intraspecific limb length variation. Cross-sectional growth data indicate that differences in bone length between Longshanks and random-bred controls are not due to prolonged growth, but to accelerated growth rates. Histomorphometric and cell proliferation assays on proximal tibial growth plates show that Longshanks' increased limb bone length is associated with an increased number of proliferative chondrocytes. In contrast, we find no differences in other growth plate cellular features known to underlie interspecific differences in limb bone size and shape, such as the rates of chondrocyte proliferation or the size and number of hypertrophic cells in the growth plate. These data suggest that small differences among individuals in the number of proliferating chondrocytes are a potentially important determinant of selectable intraspecific variation in individual limb bone lengths, independent of body size.

Genetic Dissection of a QTL Affecting Bone Geometry

Parameters of bone geometry such as width, length, and cross-sectional area are major determinants of bone strength. Although these traits are highly heritable, few genes influencing bone geometry have been identified. Here, we dissect a major quantitative trait locus (QTL) influencing femur size. This QTL was originally identified in an F2 cross between the C57BL/6J-hg/hg (HG) and CAST/EiJ strains and was referred to as femur length in high growth mice 2 (Feml2). Feml2 was located on chromosome (Chr.) 9 at ∼20 cM. Here, we show that the HG.CAST-(D9Mit249-D9Mit133)/Ucd congenic strain captures Feml2. In an F2 congenic cross, we fine-mapped the location of Feml2 to an ∼6 Mbp region extending from 57.3 to 63.3 Mbp on Chr. 9. We have identified candidates by mining the complete genome sequence of CAST/EiJ and through allele-specific expression (ASE) analysis of growth plates in C57BL/6J × CAST/EiJ F1 hybrids. Interestingly, we also find that the refined location of Feml2 overlaps a cluster of six independent genome-wide associations for human height. This work provides the foundation for the identification of novel genes affecting bone geometry.

Genetics of Skeletal Evolution in Unusually Large Mice from Gough Island

Organisms on islands often undergo rapid morphological evolution, providing a platform for understanding mechanisms of phenotypic change. Many examples of evolution on islands involve the vertebrate skeleton. Although the genetic basis of skeletal variation has been studied in laboratory strains, especially in the house mouse Mus musculus domesticus, the genetic determinants of skeletal evolution in natural populations remain poorly understood. We used house mice living on the remote Gough Island-the largest wild house mice on record-to understand the genetics of rapid skeletal evolution in nature. Compared to a mainland reference strain from the same subspecies (WSB/EiJ), the skeleton of Gough Island mice is considerably larger, with notable expansions of the pelvis and limbs. The Gough Island mouse skeleton also displays changes in shape, including elongations of the skull and the proximal vs. distal elements in the limbs. Quantitative trait locus (QTL) mapping in a large F₂ intercross between Gough Island mice and WSB/EiJ reveals hundreds of QTL that control skeletal dimensions measured at 5, 10, and/or 16 weeks of age. QTL exhibit modest, mostly additive effects, and Gough Island alleles are associated with larger skeletal size at most QTL. The QTL with the largest effects are found on a few chromosomes and affect suites of skeletal traits. Many of these loci also colocalize with QTL for body weight. The high degree of QTL colocalization is consistent with an important contribution of pleiotropy to skeletal evolution. Our results provide a rare portrait of the genetic basis of skeletal evolution in an island population and position the Gough Island mouse as a model system for understanding mechanisms of rapid evolution in nature.

Enhanced cartilage repair in 'healer' mice-New leads in the search for better clinical options for cartilage repair

Adult articular cartilage has a poor capacity to undergo intrinsic repair. Current strategies for the repair of large cartilage defects are generally unsatisfactory because the restored cartilage does not have the same resistance to biomechanical loading as authentic articular cartilage and degrades over time. Recently, an exciting new research direction, focused on intrinsic cartilage regeneration rather than fibrous repair by external means, has emerged. This review explores the new findings in this rapidly moving field as they relate to the clinical goal of restoration of structurally robust, stable and non-fibrous articular cartilage following injury.

Multivariate Analysis of Genotype-Phenotype Association

With the advent of modern imaging and measurement technology, complex phenotypes are increasingly represented by large numbers of measurements, which may not bear biological meaning one by one. For such multivariate phenotypes, studying the pairwise associations between all measurements and all alleles is highly inefficient and prevents insight into the genetic pattern underlying the observed phenotypes. We present a new method for identifying patterns of allelic variation (genetic latent variables) that are maximally associated-in terms of effect size-with patterns of phenotypic variation (phenotypic latent variables). This multivariate genotype-phenotype mapping (MGP) separates phenotypic features under strong genetic control from less genetically determined features and thus permits an analysis of the multivariate structure of genotype-phenotype association, including its dimensionality and the clustering of genetic and phenotypic variables within this association. Different variants of MGP maximize different measures of genotype-phenotype association: genetic effect, genetic variance, or heritability. In an application to a mouse sample, scored for 353 SNPs and 11 phenotypic traits, the first dimension of genetic and phenotypic latent variables accounted for >70% of genetic variation present in all 11 measurements; 43% of variation in this phenotypic pattern was explained by the corresponding genetic latent variable. The first three dimensions together sufficed to account for almost 90% of genetic variation in the measurements and for all the interpretable genotype-phenotype association. Each dimension can be tested as a whole against the hypothesis of no association, thereby reducing the number of statistical tests from 7766 to 3-the maximal number of meaningful independent tests. Important alleles can be selected based on their effect size (additive or nonadditive effect on the phenotypic latent variable). This low dimensionality of the genotype-phenotype map has important consequences for gene identification and may shed light on the evolvability of organisms.

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.

Impacts of genetic correlation on the independent evolution of body mass and skeletal size in mammals

Background

Mammals show a predictable scaling relationship between limb bone size and body mass. This relationship has a genetic basis which likely evolved via natural selection, but it is unclear how much the genetic correlation between these traits in turn impacts their capacity to evolve independently. We selectively bred laboratory mice for increases in tibia length independent of body mass, to test the hypothesis that a genetic correlation with body mass constrains evolutionary change in tibia length.

Results

Over 14 generations, we produced mean tibia length increases of 9-13%, while mean body mass was unchanged, in selectively bred mice and random-bred controls. Using evolutionary scenarios with different selection and quantitative genetic parameters, we also found that this genetic correlation impedes the rate of evolutionary change in both traits, slowing increases in tibia length while preventing decreases in body mass, despite the latter’s negative effect on fitness.

Conclusions

Overall, results from this ongoing selection experiment suggest that parallel evolution of relatively longer hind limbs among rodents, for example in the context of strong competition for resources and niche partitioning in heterogeneous environments, may have occurred very rapidly on geological timescales, in spite of a moderately strong genetic correlation between tibia length and body mass.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0258-0) contains supplementary material, which is available to authorized users.

Environmental temperature impact on bone and cartilage growth

Environmental temperature can have a surprising impact on extremity growth in homeotherms, but the underlying mechanisms have remained elusive for over a century. Limbs of animals raised at warm ambient temperature are significantly and permanently longer than those of littermates housed at cooler temperature. These remarkably consistent lab results closely resemble the ecogeographical tenet described by Allen's "extremity size rule," that appendage length correlates with temperature and latitude. This phenotypic growth plasticity could have adaptive significance for thermal physiology. Shortened extremities help retain body heat in cold environments by decreasing surface area for potential heat loss. Homeotherms have evolved complex mechanisms to maintain tightly regulated internal temperatures in challenging environments, including "facultative extremity heterothermy" in which limb temperatures can parallel ambient. Environmental modulation of tissue temperature can have direct and immediate consequences on cell proliferation, metabolism, matrix production, and mineralization in cartilage. Temperature can also indirectly influence cartilage growth by modulating circulating levels and delivery routes of essential hormones and paracrine regulators. Using an integrated approach, this article synthesizes classic studies with new data that shed light on the basis and significance of this enigmatic growth phenomenon and its relevance for treating human bone elongation disorders. Discussion centers on the vasculature as a gateway to understanding the complex interconnection between direct (local) and indirect (systemic) mechanisms of temperature-enhanced bone lengthening. Recent advances in imaging modalities that enable the dynamic study of cartilage growth plates in vivo will be key to elucidating fundamental physiological mechanisms of long bone growth regulation.

Genomic correlates of relationship QTL involved in fore- versus hind limb divergence in mice

Divergence of serially homologous elements of organisms is a common evolutionary pattern contributing to increased phenotypic complexity. Here, we study the genomic intervals affecting the variational independence of fore- and hind limb traits within an experimental mouse population. We use an advanced intercross of inbred mouse strains to map the loci associated with the degree of autonomy between fore- and hind limb long bone lengths (loci affecting the relationship between traits, relationship quantitative trait loci [rQTL]). These loci have been proposed to interact locally with the products of pleiotropic genes, thereby freeing the local trait from the variational constraint due to pleiotropic mutations. Using the known polymorphisms (single nucleotide polymorphisms [SNPs]) between the parental strains, we characterized and compared the genomic regions in which the rQTL, as well as their interaction partners (intQTL), reside. We find that these two classes of QTL intervals harbor different kinds of molecular variation. SNPs in rQTL intervals more frequently reside in limb-specific cis-regulatory regions than SNPs in intQTL intervals. The intQTL loci modified by the rQTL, in contrast, show the signature of protein-coding variation. This result is consistent with the widely accepted view that protein-coding mutations have broader pleiotropic effects than cis-regulatory polymorphisms. For both types of QTL intervals, the underlying candidate genes are enriched for genes involved in protein binding. This finding suggests that rQTL effects are caused by local interactions among the products of the causal genes harbored in rQTL and intQTL intervals. This is the first study to systematically document the population-level molecular variation underlying the evolution of character individuation.

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.

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.

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.

Synchronizing allelic effects of opposing quantitative trait loci confirmed a major epistatic interaction affecting acute lung injury survival in mice

Increased oxygen (O₂) levels help manage severely injured patients, but too much for too long can cause acute lung injury (ALI), acute respiratory distress syndrome (ARDS) and even death. In fact, continuous hyperoxia has become a prototype in rodents to mimic salient clinical and pathological characteristics of ALI/ARDS. To identify genes affecting hyperoxia-induced ALI (HALI), we previously established a mouse model of differential susceptibility. Genetic analysis of backcross and F₂ populations derived from sensitive (C57BL/6J; B) and resistant (129X1/SvJ; X1) inbred strains identified five quantitative trait loci (QTLs; Shali1-5) linked to HALI survival time. Interestingly, analysis of these recombinant populations supported opposite within-strain effects on survival for the two major-effect QTLs. Whereas Shali1 alleles imparted the expected survival time effects (i.e., X1 alleles increased HALI resistance and B alleles increased sensitivity), the allelic effects of Shali2 were reversed (i.e., X1 alleles increased HALI sensitivity and B alleles increased resistance). For in vivo validation of these inverse allelic effects, we constructed reciprocal congenic lines to synchronize the sensitivity or resistance alleles of Shali1 and Shali2 within the same strain. Specifically, B-derived Shali1 or Shali2 QTL regions were transferred to X1 mice and X1-derived QTL segments were transferred to B mice. Our previous QTL results predicted that substituting Shali1 B alleles onto the resistant X1 background would add sensitivity. Surprisingly, not only were these mice more sensitive than the resistant X1 strain, they were more sensitive than the sensitive B strain. In stark contrast, substituting the Shali2 interval from the sensitive B strain onto the X1 background markedly increased the survival time. Reciprocal congenic lines confirmed the opposing allelic effects of Shali1 and Shali2 on HALI survival time and provide unique models to identify their respective quantitative trait genes and to critically assess the apparent bidirectional epistatic interactions between these major-effect loci.

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.

Haplotype probabilities in advanced intercross populations

Advanced intercross populations, in which multiple inbred strains are mated at random for many generations, have the advantage of greater precision of genetic mapping because of the accumulation of recombination events across the multiple generations. Related designs include heterogeneous stock and the diversity outcross population. In this article, I derive the two-locus haplotype probabilities on the autosome and X chromosome with these designs. These haplotype probabilities provide the key quantities for developing hidden Markov models for the treatment of missing genotype information. I further derive the map expansion in these populations, which is the frequency of recombination breakpoints on a random chromosome.

Forward genetic approaches to understanding complex behaviors

Assigning function to genes has long been a focus of biomedical research.Even with complete knowledge of the genomic sequences of humans, mice and other experimental organisms, there is still much to be learned about gene function and control. Ablation or overexpression of single genes using knockout or transgenic technologies has provided functional annotation for many genes, but these technologies do not capture the extensive genetic variation present in existing experimental mouse populations. Researchers have only recently begun to truly appreciate naturally occurring genetic variation resulting from single nucleotide substitutions,insertions, deletions, copy number variation, epigenetic changes (DNA methylation,histone modifications, etc.) and gene expression differences and how this variation contributes to complex phenotypes. In this chapter, we will discuss the benefits and limitations of different forward genetic approaches that capture the genetic variation present in inbred mouse strains and present the utility of these approaches for mapping QTL that influence complex behavioral phenotypes.

Identification of quantitative trait loci influencing skeletal architecture in mice: emergence of Cdh11 as a primary candidate gene regulating femoral morphology

Bone strength is influenced by many properties intrinsic to bone, including its mass, geometry, and mineralization. To further advance our understanding of the genetic basis of bone-strength-related traits, we used a large (n = 815), moderately (G(4) ) advanced intercross line (AIL) of mice derived from a high-runner selection line (HR) and the C57BL/6J inbred strain. In total, 16 quantitative trait loci (QTLs) were identified that affected areal bone mineral density (aBMD) and femoral length and width. Four significant (p < .05) and one suggestive (p < .10) QTLs were identified for three aBMD measurements: total body, vertebral, and femoral. A QTL on chromosome (Chr.) 3 influenced all three aBMD measures, whereas the other four QTLs were unique to a single measure. A total of 10 significant and one suggestive QTLs were identified for femoral length (FL) and two measures of femoral width, anteroposterior (AP) and mediolateral (ML). FL QTLs were distinct from loci affecting AP and ML width, and of the 7 AP QTLs, only three affected ML. A QTL on Chr. 8 that explained 7.1% and 4.0% of the variance in AP and ML, respectively, was mapped to a 6-Mb region harboring 12 protein-coding genes. The pattern of haplotype diversity across the QTL region and expression profiles of QTL genes suggested that of the 12, cadherin 11 (Cdh11) was most likely the causal gene. These findings, when combined with existing data from gene knockouts, identify Cdh11 as a strong candidate gene within which genetic variation may affect bone morphology.

Genetic factors and diet affect long-bone length in the F34 LG,SM advanced intercross

Previous studies on the LG,SM advanced intercross line have identified approximately 40 quantitative trait loci (QTL) for long -bone (humerus, ulna, femur, and tibia) lengths. In this study, long-bone-length QTL were fine-mapped in the F(34) generation (n = 1424) of the LG,SM advanced intercross. Environmental effects were assessed by dividing the population by sex between high-fat and low-fat diets, producing eight sex/diet cohorts. We identified 145 individual bone-length QTL comprising 45 pleiotropic QTL; 69 replicated QTL from previous studies, 35 were new traits significant at previously identified loci, and 41 were novel QTL. Many QTL affected only a subset of the population based on sex and/or diet. Eight of ten known skeletal growth genes were upregulated in 3-week-old LG/J male proximal tibial growth plates relative to SM/J. The sequences of parental strains LG/J and SM/J indicated the presence of over half a million polymorphisms in the confidence intervals of these 45 QTL. We examined 526 polymorphisms and found that 97 represented radical changes to amino acid composition while 40 were predicted to be deleterious to protein function. Additional experimentation is required to understand how changes in gene regulation or protein function can alter the genetic architecture and interact with the environment to produce phenotypic variation.

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.

Developmental and genetic origins of murine long bone length variation

If we wish to understand whether development influences the rate or direction of morphological evolution, we must first understand the developmental bases of morphological variation within species. However, quantitative variation in adult morphology is the product of molecular and cellular processes unfolding from embryonic development through juvenile growth to maturity. The Atchley-Hall model provides a useful framework for dissecting complex morphologies into their component parts as a way of determining which developmental processes contribute to variation in adult form. We have examined differences in postnatal allometry and the patterns of genetic correlation between age-specific traits for ten recombinant inbred strains of mice generated from an intercross of LG/J and SM/J. Long bone length is closely tied to body size, but variation in adult morphology is more closely tied to differences in growth rate between 3 and 5 weeks of age. These analyses show that variation generated during early development is overridden by variation generated later in life. To more precisely determine the cellular processes generating this variation we then examined the cellular dynamics of long bone growth plates at the time of maximum elongation rate differences in the parent strains. Our analyses revealed that variation in long bone length is the result of faster elongation rates of the LG/J stain. The developmental bases for these differences in growth rate involve the rate of cell division and chondrocyte hypertrophy in the growth plate.

Calpain-10 is a component of the obesity-related quantitative trait locus Adip1

We previously mapped Adip1, an obesity quantitative trait locus (QTL), to the central portion of murine chromosome 1 containing the calpain-10 (Capn10) gene. Human studies have associated calpain-10 (CAPN10) variants with type 2 diabetes and various metabolic traits. We performed a quantitative hybrid complementation test (QHCT) to determine whether differences attributed to Adip1 are the result of variant Capn10 alleles in LG/J and SM/J mice. We crossed LG/J and SM/J to wild-type (C57BL/6J) and Capn10 knockout (Capn10(-/-)) mice to form four F(1) hybrid groups: LG/J by wild-type, LG/J by Capn10(-/-), SM/J by wild-type, and SM/J by Capn10(-/-). We performed a two-way ANOVA with the experimental strain, tester strain, and their interaction as the factors. Significant interaction indicates a quantitative failure to complement. We found failure to complement for fat, organ, and body weights, and leptin, female free fatty acid, and triglyceride levels. Capn10(-/-) resulted in heavier weights and higher serum levels in LG/J crosses but not in SM/J crosses. For glucose tolerance and insulin response tests, the Capn10(-/-) allele resulted in lower glucose levels in crosses with SM/J but had no effect in the LG/J crosses. Differences between the LG/J and SM/J Capn10 alleles are the likely source of some of the QTL effects mapped to Adip1 in the LG/J-by-SM/J cross. Capn10 plays an important role in regulating obesity and diabetes in mice.

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.

Genetically determined phenotype covariation networks control bone strength

To identify genes affecting bone strength, we studied how genetic variants regulate components of a phenotypic covariation network that was previously shown to accurately characterize the compensatory trait interactions involved in functional adaptation during growth. Quantitative trait loci (QTLs) regulating femoral robustness, morphologic compensation, and mineralization (tissue quality) were mapped at three ages during growth using AXB/BXA Recombinant Inbred (RI) mouse strains and adult B6-i(A) Chromosome Substitution Strains (CSS). QTLs for robustness were identified on chromosomes 8, 12, 18, and 19 and confirmed at all three ages, indicating that genetic variants established robustness postnatally without further modification. A QTL for morphologic compensation, which was measured as the relationship between cortical area and body weight, was identified on chromosome 8. This QTL limited the amount of bone formed during growth and thus acted as a setpoint for diaphyseal bone mass. Additional QTLs were identified from the CSS analysis. QTLs for robustness and morphologic compensation regulated bone structure independently (ie, in a nonpleiotropic manner), indicating that each trait may be targeted separately to individualize treatments aiming to improve strength. Multiple regression analyses showed that variation in morphologic compensation and tissue quality, not bone size, determined femoral strength relative to body weight. Thus an individual inheriting slender bones will not necessarily inherit weak bones unless the individual also inherits a gene that impairs compensation. This systems genetic analysis showed that genetically determined phenotype covariation networks control bone strength, suggesting that incorporating functional adaptation into genetic analyses will advance our understanding of the genetic basis of bone strength.

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.

Quantitative trait loci for tibial bone strength in C57BL/6J and C3H/HeJ inbred strains of mice

Three-point bending technology has been widely used in the measurement of bone strength. Quantitative trait loci (QTLs) for bone strength have been identified using mouse femurs. In this study, we investigate the use of mouse tibiae in identification of QTLs that regulate bone strength. Mouse tibiae were from a F(2) population derived from C57BL/6J (B6) and C3H/HeJ (C3H). Three-point bending was measured using ISO 4049, with the support width adjustable to accommodate specimen sizes outside the scope of ISO 4049. The strain rate is selectable from 0.05 to 500 mm per min. All stress strain diagrams are recorded and retrieved in digital electronic form. Genome scan was performed in The Jackson Laboratory (TJL). QTL mapping was conducted using Map Manager QTX software. Data show that (i) both elastic modulus (stiffness) and maximum loading (strength) value appear as normal distributions, suggesting that multiple genetic factors control the bone strength; (ii) 11 QTLs, accounting for 90% of variation for strength, have been detected. More than half QTLs of three-point bending are located on the same locations of bone density earlier identified from mouse femurs; (iii) a major QTL of femoral and vertebral bone mineral density (BMD) was not detected for bone strength of tibiae; (iv) the QTL on chromosome 4 has extremely high LOD score of 31.8 and represents 60% of the variation of bone strength; and (v) four QTLs of stiffness (chromosomes 2, 11, 15 and 19) have been identified.

Replication of long-bone length QTL in the F9-F10 LG,SM advanced intercross

Quantitative trait locus (QTL) mapping techniques are frequently used to identify genomic regions associated with variation in phenotypes of interest. However, the F(2) intercross and congenic strain populations usually employed have limited genetic resolution resulting in relatively large confidence intervals that greatly inhibit functional confirmation of statistical results. Here we use the increased resolution of the combined F(9) and F(10) generations (n = 1455) of the LG,SM advanced intercross to fine-map previously identified QTL associated with the lengths of the humerus, ulna, femur, and tibia. We detected 81 QTL affecting long-bone lengths. Of these, 49 were previously identified in the combined F(2)-F(3) population of this intercross, while 32 represent novel contributors to trait variance. Pleiotropy analysis suggests that most QTL affect three to four long bones or serially homologous limb segments. We also identified 72 epistatic interactions involving 38 QTL and 88 novel regions. This analysis shows that using later generations of an advanced intercross greatly facilitates fine-mapping of confidence intervals, resolving three F(2)-F(3) QTL into multiple linked loci and narrowing confidence intervals of other loci, as well as allowing identification of additional QTL. Further characterization of the biological bases of these QTL will help provide a better understanding of the genetics of small variations in long-bone length.

Comparison of Mandibular Phenotypic and Genetic Integration between Baboon and Mouse

In this study we compare patterns of mandibular integration between mice and baboons using both phenotypic and quantitative genetic data. Specifically, we test how well each species fits with the mosaic model of mandibular integration suggested by Atchley and Hall (Biol Rev Camb Philos Soc 66:101-157, 1991) based on developmental modules. We hypothesize that patterns of integration will be similar for mice and baboons and that both species will show strong integration within developmental modules and weaker integration between modules. Corresponding landmark data were collected from the hemi-mandibles of an advanced intercross mouse sample (N = 1239) and mandibles from a baboon sample of known pedigree from the Southwest Foundation for Biomedical Research (N = 430). We used four methods of analysis to quantify and compare the degree of mandibular integration between species including two methods based on a priori assumptions, and two a posteriori analyses. We found that patterns of integration are broadly similar for baboon and mouse mandibles, with both species displaying a modular pattern of integration. While there is a general trend of similarity in integration patterns between species, there were some marked differences. Mice are strongly correlated among distances within the coronoid process and the incisive alveolar region, whereas baboons are strongly integrated within the condylar process. We discuss the potential evolutionary implications of the similar patterns of integration between these species with an emphasis on the role of modularity.

A genome wide linkage scan of metacarpal size and geometry in the Framingham Study

Bone geometry is a significant component of bone strength, and has a clinical utility in predicting fractures and quantifying bone loss. Bone geometry is known to have a substantial genetic component. We performed linkage analysis to identify chromosomal regions governing metacarpal bone geometry. A genome-wide scan (with a set of 615 markers with spacing of approximately 5.7 cM) was performed on 1,702 individuals from 330 extended families of the Framingham Study. Midshaft width was measured and several indices calculated, namely Metacarpal Cortical Thickness (MCT), Cortical Index (MCI), and Section Modulus (MZ), using digitized X-rays of 1,380 participants (men, n = 666, mean age 55.2 yr, women, n = 714, 55.5 yr). Metacarpals 2, 3, and 4 were averaged. Heritability was significant for all indices, ranging from 0.51 to 0.72. Linkage analysis of indices adjusted for age, age(2), and estrogen status in women, identified chromosomal regions 6p21, 9p21, 11q21-q22, and Xq26-Xq27, with LOD scores >2.0. Additional adjustment for smoking, height, and BMI, generally reduced the LOD scores. Finally, bivariate linkage analysis confirmed that a QTL on chr. 6 (51 cM) was shared by midshaft width and MZ (LOD = 2.40, adjusted for all covariates). Neither MCT nor MCI shared linkage loci with width or MZ. In conclusion, we have identified chromosomal regions potentially linked to bone geometry. Genes in these regions may regulate bone geometry via effects on body size. Identification and subsequent characterization of loci for bone geometry can further elucidate the genetic contributions to bone's resistance to stress.

Pleiotropic patterns of quantitative trait loci for 70 murine skeletal traits

Quantitative trait locus (QTL) studies of a skeletal trait or a few related skeletal components are becoming commonplace, but as yet there has been no investigation of pleiotropic patterns throughout the skeleton. We present a comprehensive survey of pleiotropic patterns affecting mouse skeletal morphology in an intercross of LG/J and SM/J inbred strains (N = 1040), using QTL analysis on 70 skeletal traits. We identify 798 single-trait QTL, coalescing to 105 loci that affect on average 7-8 traits each. The number of traits affected per locus ranges from only 1 trait to 30 traits. Individual traits average 11 QTL each, ranging from 4 to 20. Skeletal traits are affected by many, small-effect loci. Significant additive genotypic values average 0.23 standard deviation (SD) units. Fifty percent of loci show codominance with heterozygotes having intermediate phenotypic values. When dominance does occur, the LG/J allele tends to be dominant to the SM/J allele (30% vs. 8%). Over- and underdominance are relatively rare (12%). Approximately one-fifth of QTL are sex specific, including many for pelvic traits. Evaluating the pleiotropic relationships of skeletal traits is important in understanding the role of genetic variation in the growth and development of the skeleton.