Whole-genome scan for linkage to bone strength and structure in inbred Fischer 344 and Lewis rats

A review of trabecular bone functional adaptation: what have we learned from trabecular analyses in extant hominoids and what can we apply to fossils?

Many of the unresolved debates in palaeoanthropology regarding evolution of particular locomotor or manipulative behaviours are founded in differing opinions about the functional significance of the preserved external fossil morphology. However, the plasticity of internal bone morphology, and particularly trabecular bone, allowing it to respond to mechanical loading during life means that it can reveal greater insight into how a bone or joint was used during an individual's lifetime. Analyses of trabecular bone have been commonplace for several decades in a human clinical context. In contrast, the study of trabecular bone as a method for reconstructing joint position, joint loading and ultimately behaviour in extant and fossil non-human primates is comparatively new. Since the initial 2D studies in the late 1970s and 3D analyses in the 1990 s, the utility of trabecular bone to reconstruct behaviour in primates has grown to incorporate experimental studies, expanded taxonomic samples and skeletal elements, and improved methodologies. However, this work, in conjunction with research on humans and non-primate mammals, has also revealed the substantial complexity inherent in making functional inferences from variation in trabecular architecture. This review addresses the current understanding of trabecular bone functional adaptation, how it has been applied to hominoids, as well as other primates and, ultimately, how this can be used to better interpret fossil hominoid and hominin morphology. Because the fossil record constrains us to interpreting function largely from bony morphology alone, and typically from isolated bones, analyses of trabecular structure, ideally in conjunction with that of cortical structure and external morphology, can offer the best resource for reconstructing behaviour in the past.

Fine mapping of bone structure and strength QTLs in heterogeneous stock rat

We previously demonstrated that skeletal structure and strength phenotypes vary considerably in heterogeneous stock (HS) rats. These phenotypes were found to be strongly heritable, suggesting that the HS rat model represents a unique genetic resource for dissecting the complex genetic etiology underlying bone fragility. The purpose of this study was to identify and localize genes associated with bone structure and strength phenotypes using 1524 adult male and female HS rats between 17 to 20 weeks of age. Structure measures included femur length, neck width, head width; femur and lumbar spine (L3-5) areas obtained by DXA; and cross-sectional areas (CSA) at the midshaft, distal femur and femoral neck, and the 5th lumbar vertebra measured by CT. In addition, measures of strength of the whole femur and femoral neck were obtained. Approximately 70,000 polymorphic SNPs distributed throughout the rat genome were selected for genotyping, with a mean linkage disequilibrium coefficient between neighboring SNPs of 0.95. Haplotypes were estimated across the entire genome for each rat using a multipoint haplotype reconstruction method, which calculates the probability of descent at each locus from each of the 8 HS founder strains. The haplotypes were then tested for association with each structure and strength phenotype via a mixed model with covariate adjustment. We identified quantitative trait loci (QTLs) for structure phenotypes on chromosomes 3, 8, 10, 12, 17 and 20, and QTLs for strength phenotypes on chromosomes 5, 10 and 11 that met a conservative genome-wide empiric significance threshold (FDR=5%; P<3×10(-6)). Importantly, most QTLs were localized to very narrow genomic regions (as small as 0.3 Mb and up to 3 Mb), each harboring a small set of candidate genes, both novel and previously shown to have roles in skeletal development and homeostasis.

Genetic regulation of bone strength: a review of animal model studies

Population- and family-based studies have established that fragility fracture risk is heritable; yet, the genome-wide association studies published to date have only accounted for a small fraction of the known variation for fracture risk of either the femur or the lumbar spine. Much work has been carried out using animal models toward finding genetic loci that are associated with bone strength. Studies using animal models overcome some of the issues associated with using patient data, but caution is needed when interpreting the results. In this review, we examine the types of tests that have been used for forward genetics mapping in animal models to identify loci and/or genes that regulate bone strength and discuss the limitations of these test methods. In addition, we present a summary of the quantitative trait loci that have been mapped for bone strength in mice, rats and chickens. The majority of these loci co-map with loci for bone size and/or geometry and thus likely dictate strength via modulating bone size. Differences in bone matrix composition have been demonstrated when comparing inbred strains of mice, and these matrix differences may be associated with differences in bone strength. However, additional work is needed to identify loci that act on bone strength at the materials level.

High-resolution genome screen for bone mineral density in heterogeneous stock rat

We previously demonstrated that skeletal mass, structure, and biomechanical properties vary considerably in heterogeneous stock (HS) rat strains. In addition, we observed strong heritability for several of these skeletal phenotypes in the HS rat model, suggesting that it represents a unique genetic resource for dissecting the complex genetics underlying bone fragility. The purpose of this study was to identify and localize genes associated with bone mineral density in HS rats. We measured bone phenotypes from 1524 adult male and female HS rats between 17 and 20 weeks of age. Phenotypes included dual-energy X-ray absorptiometry (DXA) measurements for bone mineral content and areal bone mineral density (aBMD) for femur and lumbar spine (L3-L5), and volumetric BMD measurements by CT for the midshaft and distal femur, femur neck, and fifth lumbar vertebra (L5). A total of 70,000 polymorphic single-nucleotide polymorphisms (SNPs) distributed throughout the genome were selected from genotypes obtained from the Affymetrix rat custom SNPs array for the HS rat population. These SNPs spanned the HS rat genome with a mean linkage disequilibrium coefficient between neighboring SNPs of 0.95. Haplotypes were estimated across the entire genome for each rat using a multipoint haplotype reconstruction method, which calculates the probability of descent for each genotyped locus from each of the eight founder HS strains. The haplotypes were tested for association with each bone density phenotype via a mixed model with covariate adjustment. We identified quantitative trait loci (QTLs) for BMD phenotypes on chromosomes 2, 9, 10, and 13 meeting a conservative genomewide empiric significance threshold (false discovery rate [FDR] = 5%; p < 3 × 10(-6)). Importantly, most QTLs were localized to very small genomic regions (1-3 megabases [Mb]), allowing us to identify a narrow set of potential candidate genes including both novel genes and genes previously shown to have roles in skeletal development and homeostasis.

High-throughput quantification of the mechanical competence of murine femora--a highly automated approach for large-scale genetic studies

Animal models are widely used to gain insight into the role of genetics on bone structure and function. One of the main strategies to map the genes regulating specific traits is called quantitative trait loci (QTL) analysis, which generally requires a very large number of animals (often more than 1000) to reach statistical significance. QTL analysis for mechanical traits has been mainly based on experimental mechanical testing, which, in view of the large number of animals, is time consuming. Hence, the goal of the present work was to introduce an automated method for large-scale high-throughput quantification of the mechanical properties of murine femora. Specifically, our aims were, first, to develop and validate an automated method to quantify murine femoral bone stiffness. Second, to test its high-throughput capabilities on murine femora from a large genetic study, more specifically, femora from two growth hormone (GH) deficient inbred strains of mice (B6-lit/lit and C3.B6-lit/lit) and their first (F1) and second (F2) filial offsprings. Automated routines were developed to convert micro-computed tomography (micro-CT) images of femora into micro-finite element (micro-FE) models. The method was experimentally validated on femora from C57BL/6J and C3H/HeJ mice: for both inbred strains the micro-FE models closely matched the experimentally measured bone stiffness when using a single tissue modulus of 13.06 GPa. The mechanical analysis of the entire dataset (n=1990) took approximately 44 CPU hours on a supercomputer. In conclusion, our approach, in combination with QTL analysis could help to locate genes directly involved in controlling bone mechanical competence.

Identification of candidate gene regions in the rat by co-localization of QTLs for bone density, size, structure and strength

Susceptibility to osteoporotic fracture is influenced by genetic factors that can be dissected by whole-genome linkage analysis in experimental animal crosses. The aim of this study was to characterize quantitative trait loci (QTLs) for biomechanical and two-dimensional dual-energy X-ray absorptiometry (DXA) phenotypes in reciprocal F2 crosses between diabetic GK and normo-glycemic F344 rat strains and to identify possible co-localization with previously reported QTLs for bone size and structure. The biomechanical measurements of rat tibia included ultimate force, stiffness and work to failure while DXA was used to characterize tibial area, bone mineral content (BMC) and areal bone mineral density (aBMD). F2 progeny (108 males, 98 females) were genotyped with 192 genome-wide markers followed by sex- and reciprocal cross-separated whole-genome QTL analyses. Significant QTLs were identified on chromosome 8 (tibial area; logarithm of odds (LOD) = 4.7 and BMC; LOD = 4.1) in males and on chromosome 1 (stiffness; LOD = 5.5) in females. No QTLs showed significant sex-specific interactions. In contrast, significant cross-specific interactions were identified on chromosome 2 (aBMD; LOD = 4.7) and chromosome 6 (BMC; LOD = 4.8) for males carrying F344mtDNA, and on chromosome 15 (ultimate force; LOD = 3.9) for males carrying GKmtDNA, confirming the effect of reciprocal cross on osteoporosis-related phenotypes. By combining identified QTLs for biomechanical-, size- and qualitative phenotypes (pQCT and 3D CT) from the same population, overlapping regions were detected on chromosomes 1, 3, 4, 6, 8 and 10. These are strong candidate regions in the search for genetic risk factors for osteoporosis.

Heterogeneous stock rat: a unique animal model for mapping genes influencing bone fragility

Previously, we demonstrated that skeletal mass, structure and biomechanical properties vary considerably among 11 different inbred rat strains. Subsequently, we performed quantitative trait loci (QTL) analysis in four inbred rat strains (F344, LEW, COP and DA) for different bone phenotypes and identified several candidate genes influencing various bone traits. The standard approach to narrowing QTL intervals down to a few candidate genes typically employs the generation of congenic lines, which is time consuming and often not successful. A potential alternative approach is to use a highly genetically informative animal model resource capable of delivering very high resolution gene mapping such as Heterogeneous stock (HS) rat. HS rat was derived from eight inbred progenitors: ACI/N, BN/SsN, BUF/N, F344/N, M520/N, MR/N, WKY/N and WN/N. The genetic recombination pattern generated across 50 generations in these rats has been shown to deliver ultra-high even gene-level resolution for complex genetic studies. The purpose of this study is to investigate the usefulness of the HS rat model for fine mapping and identification of genes underlying bone fragility phenotypes. We compared bone geometry, density and strength phenotypes at multiple skeletal sites in HS rats with those obtained from five of the eight progenitor inbred strains. In addition, we estimated the heritability for different bone phenotypes in these rats and employed principal component analysis to explore relationships among bone phenotypes in the HS rats. Our study demonstrates that significant variability exists for different skeletal phenotypes in HS rats compared with their inbred progenitors. In addition, we estimated high heritability for several bone phenotypes and biologically interpretable factors explaining significant overall variability, suggesting that the HS rat model could be a unique genetic resource for rapid and efficient discovery of the genetic determinants of bone fragility.

Genetic loci for bone architecture determined by three-dimensional CT in crosses with the diabetic GK rat

The F344 rat carries alleles contributing to bone fragility while the GK rat spontaneously develops type-2 diabetes. These characteristics make F344×GK crosses well suited for the identification of genes related to bone size and allow for future investigation on the association with type-2 diabetes. The aim of this study was to identify quantitative trait loci (QTLs) for bone size phenotypes measured by a new application of three-dimensional computed tomography (3DCT) and to investigate the effects of sex- and reciprocal cross. Tibia from male and female GK and F344 rats, representing the parental, F1 and F2 generations, were examined with 3DCT and analyzed for: total and cortical volumetric BMD, straight and curved length, peri- and endosteal area at mid-shaft. F2 progeny (108 male and 98 female) were genotyped with 192 genome-wide microsatellite markers (average distance 10 cM). Sex- and reciprocal cross-separated QTL analyses were performed for the identification of QTLs linked to 3DCT phenotypes and true interactions were confirmed by likelihood ratio analysis in all F2 animals. Several genome-wide significant QTLs were found in the sex- and reciprocal cross-separated progeny on chromosomes (chr) 1, 3, 4, 9, 10, 14, and 17. Overlapping QTLs for both males and females in the (GK×F344)F2 progeny were located on chr 1 (39-67 cM). This region confirms previously reported pQCT QTLs and overlaps loci for fasting glucose. Sex separated linkage analysis confirmed a male specific QTL on chr 9 (67-82 cM) for endosteal area at the fibula site. Analyses separating the F2 population both by sex and reciprocal cross identified cross specific QTLs on chr 14 (males) and chr 3 and 4 (females). Two loci, chr 4 and 6, are unique to 3DCT and separate from pQCT generated loci. The 3DCT method was highly reproducible and provided high precision measurements of bone size in the rat enabling identification of new sex- and cross-specific loci. The QTLs on chr 1 indicate potential genetic association between bone-related phenotypes and traits affecting type-2 diabetes. The results illustrate the complexity of the genetic architecture of bone size phenotypes and demonstrate the importance of complementary methods for bone analysis.

Genetics of osteoporosis

Osteoporosis is a common disease with a strong genetic component characterized by reduced bone mass, defects in the microarchitecture of bone tissue, and an increased risk of fragility fractures. Twin and family studies have shown high heritability of bone mineral density (BMD) and other determinants of fracture risk such as ultrasound properties of bone, skeletal geometry, and bone turnover. Osteoporotic fractures also have a heritable component, but this reduces with age as environmental factors such as risk of falling come into play. Susceptibility to osteoporosis is governed by many different genetic variants and their interaction with environmental factors such as diet and exercise. Notable successes in identification of genes that regulate BMD have come from the study of rare Mendelian bone diseases characterized by major abnormalities of bone mass where variants of large effect size are operative. Genome-wide association studies have also identified common genetic variants of small effect size that contribute to regulation of BMD and fracture risk in the general population. In many cases, the loci and genes identified by these studies had not previously been suspected to play a role in bone metabolism. Although there has been extensive progress in identifying the genes and loci that contribute to the regulation of BMD and fracture over the past 15 yr, most of the genetic variants that regulate these phenotypes remain to be discovered.

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.

Identification of genes influencing skeletal phenotypes in congenic P/NP rats

We previously showed that alcohol-preferring (P) rats have higher bone density than alcohol-nonpreferring (NP) rats. Genetic mapping in P and NP rats identified a major quantitative trait locus (QTL) between 4q22 and 4q34 for alcohol preference. At the same location, several QTLs linked to bone density and structure were detected in Fischer 344 (F344) and Lewis (LEW) rats, suggesting that bone mass and strength genes might cosegregate with genes that regulate alcohol preference. The aim of this study was to identify the genes segregating for skeletal phenotypes in congenic P and NP rats. Transfer of the NP chromosome 4 QTL into the P background (P.NP) significantly decreased areal bone mineral density (aBMD) and volumetric bone mineral density (vBMD) at several skeletal sites, whereas transfer of the P chromosome 4 QTL into the NP background (NP.P) significantly increased bone mineral content (BMC) and aBMD in the same skeletal sites. Microarray analysis from the femurs using Affymetrix Rat Genome arrays revealed 53 genes that were differentially expressed among the rat strains with a false discovery rate (FDR) of less than 10%. Nine candidate genes were found to be strongly correlated (r(2) > 0.50) with bone mass at multiple skeletal sites. The top three candidate genes, neuropeptide Y (Npy), alpha synuclein (Snca), and sepiapterin reductase (Spr), were confirmed using real-time quantitative PCR (qPCR). Ingenuity pathway analysis revealed relationships among the candidate genes related to bone metabolism involving beta-estradiol, interferon-gamma, and a voltage-gated calcium channel. We identified several candidate genes, including some novel genes on chromosome 4 segregating for skeletal phenotypes in reciprocal congenic P and NP rats.

Genetics of osteoporosis

Osteoporosis is a common disease with a strong genetic component characterized by reduced bone mass and an increased risk of fragility fractures. Twin and family studies have shown that the heritability of bone mineral density and other determinants of fracture risk, such as ultrasound properties of bone, skeletal geometry, and bone turnover, is high, although heritability of fracture is modest. Many different genetic variants contribute to the regulation of these phenotypes. Most are common variants of small effect size, but there is evidence that rare variants of large effect size also contribute in some individuals. Many of the genes that regulate susceptibility to osteoporosis have been identified through studies of rare bone diseases, but genome-wide association studies have also been successful in identifying genes that predispose to osteoporosis. Although there has been extensive progress in this area over the past 10 years, most of the genetic variants that regulate susceptibility to osteoporosis remain to be discovered.

Linkage mapping of femoral material properties in a reciprocal intercross of HcB-8 and HcB-23 recombinant mouse strains

Skeletal fragility is an important health problem with a large genetic component. We performed a 603 animal F2 reciprocal intercross of the recombinant congenic strains HcB-8 and HcB-23 to genetically map quantitative trait loci (QTLs) for tissue-level femoral biomechanical performance. These included elastic and post-yield strain, Young's modulus, elastic and maximum stress, and toughness and were calculated from 3-point bend testing of femora by the application of standard beam equations. We mapped these with R/qtl and QTL Cartographer and established significance levels empirically by permutation testing. Significant QTLs for at least one trait are present on chromosomes 1, 6, and 10 in the full F2 population, with additional QTLs evident in subpopulations defined by sex and cross direction. On chromosome 10, we find a QTL for post-yield strain and toughness, phenotypes that have not been mapped previously. Notably, the HcB-8 allele at this QTL increases post-yield strain and toughness, but decreases bone mineral density (BMD), while the material property QTLs on chromosomes 1, 6, and at a second chromosome 10 QTL are independent of BMD. We find significant sex x QTL and cross direction x QTL interactions. A robust, pleiotropic chromosome 4 QTL that we previously reported at the whole-bone level showed no evidence of linkage at the tissue-level, supporting our interpretation that modeling capacity is its primary phenotype. Our data demonstrate an inverse relationship between femoral perimeter and Young's modulus, with R(2)=0.27, supporting the view that geometric and material bone properties are subject to an integrated set of regulatory mechanisms. Mapping QTLs for tissue-level biomechanical performance advances understanding of the genetic basis of bone quality.

Heritability of lumbar trabecular bone mechanical properties in baboons

Genetic effects on mechanical properties have been demonstrated in rodents, but not confirmed in primates. Our aim was to quantify the proportion of variation in vertebral trabecular bone mechanical properties that is due to the effects of genes. L3 vertebrae were collected from 110 females and 46 male baboons (6-32 years old) from a single extended pedigree. Cranio-caudally oriented trabecular bone cores were scanned with microCT then tested in monotonic compression to determine apparent ultimate stress, modulus, and toughness. Age and sex effects and heritability (h(2)) were assessed using maximum likelihood-based variance components methods. Additive effects of genes on residual trait variance were significant for ultimate stress (h(2)=0.58), toughness (h(2)=0.64), and BV/TV (h(2)=0.55). When BV/TV was accounted for, the residual variance in ultimate stress accounted for by the additive effects of genes was no longer significant. Toughness, however, showed evidence of a non-BV/TV-related genetic effect. Overall, maximum stress and modulus show strong genetic effects that are nearly entirely due to bone volume. Toughness shows strong genetic effects related to bone volume and shows additional genetic effects (accounting for 10% of the total trait variance) that are independent of bone volume. These results support continued use of bone volume as a focal trait to identify genes related to skeletal fragility, but also show that other focal traits related to toughness and variation in the organic component of bone matrix will enhance our ability to find additional genes that are particularly relevant to fatigue-related fractures.

Cross-sectional geometry of the femoral midshaft in baboons is heritable

A great deal of research into the determinants of bone strength has unequivocally demonstrated that variation in bone strength is highly subject to genetic factors. Increasing attention in skeletal genetic studies is being paid to indicators of bone quality that complement studies of BMD, including studies of the genetic control of bone geometry. The aim of this study is to investigate the degree to which normal population-level variation in femoral midshaft geometry in a population of pedigreed baboons (Papio hamadryas spp.) can be attributed to the additive effect of genes. Using 110 baboons (80 females, 30 males), we 1) characterize normal variation in midshaft geometry of the femur with regard to age and sex, and 2) determine the degree to which the residual variation is attributable to additive genetic effects. Cross-sectional area (CSA), minimum (I(MIN)) and maximum (I(MAX)) principal moments of inertia, and polar moment of inertia (J) were calculated from digitized images of transverse midshaft sections. Maximum likelihood-based variance decomposition methods were used to estimate the mean effects of age, sex, and genes. Together age and sex effects account for approximately 56% of the variance in each property. In each case the effect of female sex is negative and that of age is positive, although of a lower magnitude than the effect of female sex. Increased age is associated with decreased mean cross-sectional geometry measures in the oldest females. Residual h(2) values range from 0.36 to 0.50, reflecting genetic effects accounting for 15% to 23% of the total phenotypic variance in individual properties. This study establishes the potential of the baboon model for the identification of genes that regulate bone geometric properties in primates. This model is particularly valuable because it allows for experimental designs, environmental consistency, availability of tissues, and comprehensive assessments of multiple integrated bone phenotypes that are not possible in human populations. The baboon is of particular importance in genetic studies, because it provides results that are likely highly relevant to the human condition due to the phylogenetic proximity of baboons to humans.

Differentially expressed genes strongly correlated with femur strength in rats

The region of chromosome 1q33-q54 harbors quantitative trait loci (QTL) for femur strength in COPxDA and F344xLEW F2 rats. The purpose of this study is to identify the genes within this QTL region that contribute to the variation in femur strength. Microarray analysis was performed using RNA extracted from femurs of COP, DA, F344 and LEW rats. Genes differentially expressed in the 1q33-q54 region among these rat strains were then ranked based on the strength of correlation with femur strength in F2 animals derived from these rats. A total of 214 genes in this QTL region were differentially expressed among all rat strains, and 81 genes were found to be strongly correlated (r(2)>0.50) with femur strength. Of these, 12 candidate genes were prioritized for further validation, and 8 of these genes (Ifit3, Ppp2r5b, Irf7, Mpeg1, Bloc1s2, Pycard, Sec23ip, and Hps6) were confirmed by quantitative PCR (qPCR). Ingenuity Pathway Analysis suggested that these genes were involved in interferon alpha, nuclear factor-kappa B (NFkB), extracellular signal-related kinase (ERK), hepatocyte nuclear factor 4 alpha (HNF4A) and tumor necrosis factor (TNF) pathways.

Serum IGF-1 determines skeletal strength by regulating subperiosteal expansion and trait interactions

Strong correlations between serum IGF-1 levels and fracture risk indicate that IGF-1 plays a critical role in regulating bone strength. However, the mechanism by which serum IGF-1 regulates bone structure and fracture resistance remains obscure and cannot be determined using conventional approaches. Previous analysis of adult liver-specific IGF-1-deficient (LID) mice, which exhibit 75% reductions in serum IGF-1 levels, showed reductions in periosteal circumference, femoral cross-sectional area, cortical thickness, and total volumetric BMD. Understanding the developmental sequences and the resultant anatomical changes that led to this adult phenotype is the key for understanding the complex relationship between serum IGF-1 levels and fracture risk. Here, we identified a unique developmental pattern of morphological and compositional traits that contribute to bone strength. We show that reduced bone strength associated with low levels of IGF-1 in serum (LID mice) result in impaired subperiosteal expansion combined with impaired endosteal apposition and lack of compensatory changes in mineralization throughout growth and aging. We show that serum IGF-1 affects cellular activity differently depending on the cortical surface. Last, we show that chronic reductions in serum IGF-1 indirectly affect bone strength through its effect on the marrow myeloid progenitor cell population. We conclude that serum IGF-1 not only regulates bone size, shape, and composition during ontogeny, but it plays a more fundamental role-that of regulating an individual's ability to adapt its bone structure to mechanical loads during growth and development.

Linkage screen for BMD phenotypes in male and female COP and DA rat strains

Because particular inbred strains of experimental animals are informative for only a subset of the genes underlying variability in BMD, we undertook a genome screen to identify quantitative trait loci (QTLs) in 828 F(2) progeny (405 males and 423 females) derived from the Copenhagen 2331 (COP) and dark agouti (DA) strains of rats. This screen was performed to complement our study in female Fischer 344 (F344) and Lewis (LEW) rats and to further delineate the factors underlying the complex genetic architecture of BMD in the rat model. Microsatellite genotyping was performed using markers at an average density of 20 cM. BMD was measured by pQCT and DXA. These data were analyzed in the R/qtl software to detect QTLs acting in both sexes as well as those having sex-specific effects. A QTL was detected in both sexes on chromosome 18 for midfemur volumetric BMD (vBMD; genome-wide, p < 0.01). On distal chromosome 1, a QTL was found for femur and vertebral aBMD as well as distal femur vBMD, and this QTL appears distinct from the proximal chromosome 1 QTL impacting BMD in our F344/LEW cross. Additional aBMD and vBMD QTLs and several sex-specific QTLs were also detected. These included a male-specific QTL (p < 0.01) on chromosome 8 and a female-specific QTL on chromosomes 7 and 14 (p < 0.01). Few of the QTLs identified showed overlap with the significant QTLs from the F344/LEW cross. These results confirm that the genetic influence on BMD in the rat model is quite complex and would seem to be influenced by a number of different genes, some of which have sex-specific effects.

Genomic expression analysis of rat chromosome 4 for skeletal traits at femoral neck

Hip fracture is the most devastating osteoporotic fracture type with significant morbidity and mortality. Several studies in humans and animal models identified chromosomal regions linked to hip size and bone mass. Previously, we identified that the region of 4q21-q41 on rat chromosome (Chr) 4 harbors multiple femoral neck quantitative trait loci (QTLs) in inbred Fischer 344 (F344) and Lewis (LEW) rats. The purpose of this study is to identify the candidate genes for femoral neck structure and density by correlating gene expression in the proximal femur with the femoral neck phenotypes linked to the QTLs on Chr 4. RNA was extracted from proximal femora of 4-wk-old rats from F344 and LEW strains, and two other strains, Copenhagen 2331 and Dark Agouti, were used as a negative control. Microarray analysis was performed using Affymetrix Rat Genome 230 2.0 arrays. A total of 99 genes in the 4q21-q41 region were differentially expressed (P < 0.05) among all strains of rats with a false discovery rate <10%. These 99 genes were then ranked based on the strength of correlation between femoral neck phenotypes measured in F2 animals, homozygous for a particular strain's allele at the Chr 4 QTL and the expression level of the gene in that strain. A total of 18 candidate genes were strongly correlated (r(2) > 0.50) with femoral neck width and prioritized for further analysis. Quantitative PCR analysis confirmed 14 of 18 of the candidate genes. Ingenuity pathway analysis revealed several direct or indirect relationships among the candidate genes related to angiogenesis (VEGF), bone growth (FGF2), bone formation (IGF2 and IGF2BP3), and resorption (TNF). This study provides a shortened list of genetic determinants of skeletal traits at the hip and may lead to novel approaches for prevention and treatment of hip fracture.

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.

Genetic loci affecting bone structure and strength in inbred COP and DA rats

Previous studies have shown that the Copenhagen 2331 (COP) and Dark Agouti (DA) rats have significant differences in bone structure and strength despite their similar body mass. Thus, these inbred rat strains may provide a unique resource to identify the genetics underlying the phenotypic variation in bone fragility. A sample of 828 (405 males and 423 females) COPxDA F2 progeny had extensive phenotyping for bone structure measures including cortical bone area and polar moment of inertia at the femur midshaft and total, cortical and trabecular bone areas, for the lumbar vertebra 5 (L5). Bone strength phenotypes included ultimate force, stiffness and work to failure of femur and L5. These skeletal phenotypes were measured using peripheral quantitative computed tomography (pQCT) and mechanical testing. A whole-genome screen was conducted in the F2 rats, using microsatellite markers spaced at approximately 20 cM intervals. Genetic marker maps were generated from the F2 data and used for genome-wide linkage analyses to detect linkage to the bone structure and strength phenotypes. Permutation testing was employed to obtain the thresholds for genome-wide significance (p<0.01). Significant QTL for femur structure and strength were identified on chromosome (Chr) 1 with a maximum LOD score of 33.5; evidence of linkage was found in both the male and female rats. In addition, Chrs 6, 7, 10, 13, 15 and 18 were linked to femur midshaft structure. QTL linked to femur strength were identified on Chrs 5 and 10. For L5 vertebrae, Chrs 2, 16, and 18 harbored QTL for cortical structure and trabecular structure for L5 was linked to Chrs 1, 7, 12, and 18. One female-specific QTL for femur ultimate force was identified on Chr 5, and two male-specific QTL for L5 cortical area were found on Chrs 2 and 18. Our study demonstrates strong evidence of linkage for bone structure and strength to multiple rat chromosomes.

Genetics of osteoporosis

Osteoporosis is a common disease with a strong genetic component characterised by reduced bone mass and an increased risk of fragility fractures. Twin and family studies have shown that genetic factors contribute to osteoporosis by influencing bone mineral density (BMD), and other phenotypes that are associated with fracture risk, although the heritability of fracture itself is modest. Linkage studies have identified several quantitative trait loci that regulate BMD but most causal genes remain to be identified. In contrast, linkage studies in monogenic bone diseases have been successful in gene identification, and polymorphisms in many of these genes have been found to contribute to the regulation of bone mass in the normal population. Population-based studies have identified polymorphisms in several candidate genes that have been associated with bone mass or osteoporotic fracture, although individually these polymorphisms only account for a small amount of the genetic contribution to BMD regulation. Environmental factors such as diet and physical activity are also important determinants of BMD, and in some cases specific nutrients have been found to interact with genetic polymorphisms to regulate BMD. From a clinical standpoint, advances in knowledge about the genetic basis of osteoporosis are likely to be important in increasing the understanding of the pathophysiology of the disease; providing new genetic markers with which to assess fracture risk and in identifying genes and pathways that form molecular targets for the design of the next generation of drug treatments.

Genetic effects on bone mechanotransduction in congenic mice harboring bone size and strength quantitative trait loci

The degree to which bone tissue responds to mechanical loading events is partially under genetic control. We assess the contribution of three genetic loci (QTLs linked to bone geometry and strength)--located on mouse Chrs. 1, 8, and 13--to mechanically stimulated bone formation, through in vivo skeletal loading of congenic strains. Bone size was not consistently associated with mechano-responsiveness, indicating that the genetic regulation of mechanotransduction is a complex process that involves a number of genes and is sex-specific.

INTRODUCTION

We showed previously that C57BL/6J (B6) mice are more responsive to mechanical stimulation than C3H/HeJ (C3H) mice and that B6 mice harboring a 40-Mb region of distal C3H Chromosome (Chr.) 4 are more responsive to mechanical stimulation than are fully B6 mice. Here, we assess the contribution of three more genetic loci--located on mouse Chrs. 1, 8, and 1--to mechanically stimulated bone formation.

MATERIALS AND METHODS

Three congenic mouse strains were created in which a region of mouse Chr. 1 (approximately 64 cM; 150 Mb), Chr. 8 (approximately 45 cM; 86 Mb), or Chr. 13 (approximately 24 cM; 42 Mb) was moved from C3H stock to a B6 background through selective breeding over nine generations. The regions moved to the B6 background correspond to three of several quantitative trait loci (QTLs) identified for bone size and strength. The resulting congenic mice were 99% B6, with the remaining genomic DNA comprised of the Chr. 1, 8, or 13 QTLs of interest. Male and female congenic (1T, 8T, and 13B) and B6 control mice were subjected to in vivo loading of the right ulna at one of three different load magnitudes. A separate set of animals from each group had strain gauges applied at the ulnar midshaft to estimate strain at each loading level. Loading was conducted once per day for 3 days (60 cycles/d; 2 Hz). Fluorochrome labels were injected intraperitoneally 4 and 11 days after loading began. Using quantitative histomorphometry, bone formation rates were measured in loaded (right) and control (left) ulnas.

RESULTS

All male congenic mice exhibited significantly reduced mechano-responsiveness compared with male B6 controls, but the same comparison among females yielded no difference from controls, with the exception of the 1T congenics, which showed increased responsiveness to loading. Among the congenic strains, smaller bone size was not consistently associated with reduced mechano-responsiveness.

CONCLUSIONS

Our results indicate that the genetic regulation of mechanotransduction is a complex process that involves a number of genes and is sex-specific. Our data might explain why different individuals can engage in similar exercise protocols yet experience different results in terms of bone mass accrual.

Quantitative trait loci that determine mouse tibial nanoindentation properties in an F2 population derived from C57BL/6J x C3H/HeJ

Use of nanoindentation technology to identify quantitative trait loci (QTL) that regulate bone properties represents a novel approach to improving our understanding of molecular mechanisms that control bone matrix properties. Tibiae for QTL mapping were from an F2 population derived from C57BL/6J and C3H/HeJ. A nanoindenter (Triboindenter; Hysitron, Minneapolis, MN) was used to conduct indentation tests on transverse sections. Genotyping was performed in The Jackson Laboratory. QTL mapping was conducted using software. We found that (1) tibiae from mice at 16 weeks of age were mature and suitable for measurement by a nanoindentor; (2) both stiffness modulus and hardness modulus in the F2 population appeared to have normal distributions, which suggested that multiple genetic factors control the bone properties; and (3) QTL for hardness were identified from five chromosomes (Chr 8, 12, 13, 17, and 19) and for stiffness, from four chromosomes (Chr 3, 8, 12, and 13). Among all detected QTL, one at the same location on Chr 12 was detected for both hardness and stiffness data. It explained the highest percentage of phenotypic variation in bone properties. Using nanoindentation technology to identify QTL that regulate bone properties yielded as many as six different chromosomal regions. Although the actual genes remain to be identified, nanoindentation will contribute to our understanding of molecular mechanisms and normal development processes that control the matrix properties of bone.

Large-scale genome-wide linkage analysis for loci linked to BMD at different skeletal sites in extreme selected sibships

Few genome-wide linkage studies of osteoporosis have been conducted in the Asian population. We performed a genome-wide scan involving 3093 adult siblings with at least one sib-pair extremely concordant or discordant for hip BMD. Our results indicated four genome-wide significant QTLs for BMD. In comparison with 12 previous reported linkage studies, we reveal novel linkage regions that have reaching global significance.

INTRODUCTION

The genetic basis for osteoporosis has been firmly established, but efforts to identify genes associated with this complex trait have been incomplete, especially in Asian populations. The purpose of this study was to identify quantitative trait loci (QTLs) for BMD in a Chinese population.

MATERIALS AND METHODS

We performed a genome-wide scan involving 3093 siblings 25-64 years of age from 941 families, with at least one sib-pair extreme concordant or discordant for total hip BMD from a large community-based cohort (n = 23,327) in Anhui, China. Linkage analysis was performed on BMD residuals adjusted for age, height, weight, occupation, cigarette smoking, physical activity, and alcohol consumption using the revised Haseman-Elston regression-based linkage model.

RESULTS

Our results revealed significant QTLs on chromosome 7p21.2 for femoral neck BMD (LOD = 3.68) and on chromosome 2q24.3 for total hip BMD (LOD = 3.65). Suggestive linkage regions were found to overlap among different skeletal sites on chromosomes 2q, 7p, and 16q. Sex-specific linkage analysis further revealed a significant QTL for lumbar spine BMD on chromosome 13q21.1 (LOD = 3.62) in women only. When performing multivariate linkage analysis by combining BMDs at four skeletal sites (i.e., whole body, total hip, femoral neck, and lumbar spine BMD), an additional significant QTL was found at chromosome 5q21.2 (LOD = 4.56). None of these significant QTLs found in our study overlapped with major QTLs reported by other studies.

CONCLUSIONS

This study reveals four novel QTLs in a Chinese population and suggests that BMD at different skeletal sites may also share common genetic determinants.

Identification of a quantitative trait locus on rat chromosome 4 that is strongly linked to femoral neck structure and strength

Risk factors for osteoporotic hip fracture include reduced bone mineral density and poor structure of the femoral neck, both of which are heritable traits. Previously, we showed that despite similar body size, Fischer 344 (F344) rats have significantly different skeletal traits compared with Lewis (LEW) rats. To identify a gene or genes regulating fracture risk at the femoral neck, we mapped quantitative trait loci (QTL) for femoral neck density and structure phenotypes using a 595 F2 progeny derived from the inbred F344 and LEW strains of rats. Femoral neck phenotypes included volumetric bone mineral density (vBMD), neck width, femoral neck cross-sectional area and polar moment of inertia (Ip). A 20-cM genome-wide scan was performed using 118 microsatellite markers and linkage analysis was conducted to identify chromosomal regions harbor QTL for femoral neck phenotypes. Strong evidence of linkage (P<0.01) to femoral neck vBMD was observed on chromosomes (Chrs) 1, 2, 4, 5, 7, 10 and 15. QTL affecting femoral neck structure and biomechanical properties were detected only on Chr 4 where the F344 alleles were shown to improve femoral neck structure, whereas these alleles had no effect on bone measurements at the lumbar spine and only modest effects at the femoral midshaft. In contrast, QTL on Chrs 1, 2 and 10 affected multiple skeletal sites. Several QTL regions in this study are homologous to human chromosomal regions, where linkage to femoral neck and related phenotypes has been reported previously. These findings represent an important first step in localizing and identifying genes that influence hip fragility.

Genome-wide scan identified QTLs underlying femoral neck cross-sectional geometry that are novel studied risk factors of osteoporosis

A genome-wide screen was conducted using a large white sample to identify QTLs for FNCS geometry. We found significant linkage of FNCS parameters to 20q12 and Xq25, plus significant epistatic interactions and sex-specific QTLs influencing FNCS geometry variation.

INTRODUCTION

Bone geometry, a highly heritable trait, is a critical component of bone strength that significantly determines osteoporotic fracture risk. Specifically, femoral neck cross-sectional (FNCS) geometry is significantly associated with hip fracture risk as well as genetic factors. However, genetic research in this respect is still in its infancy.

MATERIALS AND METHODS

To identify the underlying genomic regions influencing FNCS variables, we performed a remarkably large-scale whole genome linkage scan involving 3998 individuals from 434 pedigrees for four FNCS geometry parameters, namely buckling ratio (BR), cross-sectional area (CSA), cortical thickness (CT), and section modulus (Z). The major statistical approach adopted is the variance component method implemented in SOLAR.

RESULTS

Significant linkage evidence (threshold LOD = 3.72 after correction for tests of multiple phenotypes) was found in the regions of 20q12 and Xq25 for CT (LOD = 4.28 and 3.90, respectively). We also identified eight suggestive linkage signals (threshold LOD = 2.31 after correction for multiple tests) for the respective geometry traits. The above findings were supported by principal component linkage analysis. Of them, 20q12 was of particular interest because it was linked to multiple FNCS geometry traits and significantly interacted with five other genomic loci to influence CSA variation. The effects of 20q12 on FNCS geometry were present in both male and female subgroups. Subgroup analysis also revealed the presence of sex-specific quantitative trait loci (QTLs) for FNCS traits in the regions such as 2p14, 3q26, 7q21 and 15q21.

CONCLUSIONS

Our findings laid a foundation for further replication and fine-mapping studies as well as for positional and functional candidate gene studies, aiming at eventually finding the causal genetic variants and hidden mechanisms concerning FNCS geometry variation and the associated hip fractures.