Quantitative trait loci for periosteal circumference (PC): identification of single loci and epistatic effects in F2 MRL/SJL mice

Widespread differential maternal and paternal genome effects on fetal bone phenotype at mid-gestation

Parent-of-origin-dependent (epi)genetic factors are important determinants of prenatal development that program adult phenotype. However, data on magnitude and specificity of maternal and paternal genome effects on fetal bone are lacking. We used an outbred bovine model to dissect and quantify effects of parental genomes, fetal sex, and nongenetic maternal effects on the fetal skeleton and analyzed phenotypic and molecular relationships between fetal muscle and bone. Analysis of 51 bone morphometric and weight parameters from 72 fetuses recovered at day 153 gestation (54% term) identified six principal components (PC1-6) that explained 80% of the variation in skeletal parameters. Parental genomes accounted for most of the variation in bone wet weight (PC1, 72.1%), limb ossification (PC2, 99.8%), flat bone size (PC4, 99.7%), and axial skeletal growth (PC5, 96.9%). Limb length showed lesser effects of parental genomes (PC3, 40.8%) and a significant nongenetic maternal effect (gestational weight gain, 29%). Fetal sex affected bone wet weight (PC1, p < 0.0001) and limb length (PC3, p < 0.05). Partitioning of variation explained by parental genomes revealed strong maternal genome effects on bone wet weight (74.1%, p < 0.0001) and axial skeletal growth (93.5%, p < 0.001), whereas paternal genome controlled limb ossification (95.1%, p < 0.0001). Histomorphometric data revealed strong maternal genome effects on growth plate height (98.6%, p < 0.0001) and trabecular thickness (85.5%, p < 0.0001) in distal femur. Parental genome effects on fetal bone were mirrored by maternal genome effects on fetal serum 25-hydroxyvitamin D (96.9%, p < 0.001) and paternal genome effects on alkaline phosphatase (90.0%, p < 0.001) and their correlations with maternally controlled bone wet weight and paternally controlled limb ossification, respectively. Bone wet weight and flat bone size correlated positively with muscle weight (r = 0.84 and 0.77, p < 0.0001) and negatively with muscle H19 expression (r = -0.34 and -0.31, p < 0.01). Because imprinted maternally expressed H19 regulates growth factors by miRNA interference, this suggests muscle-bone interaction via epigenetic factors.

Impact of the environment on the skeleton: is it modulated by genetic factors?

The etiology of skeletal disease is driven by genetic and environmental factors. Genome-wide association studies (GWAS) of osteoporotic phenotypes have identified novel candidate genes, but have only uncovered a small proportion of the trait variance explained. This "missing heritability" is caused by several factors, including the failure to consider gene-by-environmental (G*E) interactions. Some G*E interactions have been investigated, but new approaches to integrate environmental data into genomic studies are needed. Advances in genotyping and meta-analysis techniques now allow combining genotype data from multiple studies, but the measurement of key environmental factors in large human cohorts still lags behind, as do the statistical tools needed to incorporate these measures in genome-wide association meta-studies. This review focuses on discussing ways to enhance G*E interaction studies in humans and how the use of rodent models can inform genetic studies. Understanding G*E interactions will provide opportunities to effectively target intervention strategies for individualized therapy.

Skeletal biology over the life span: a view from the surfaces

The biocultural interpretation of skeletal remains is based upon the foundation of skeletal biology. In this review we examine the current state of skeletal biology research outside of the mainstream anthropology literature. The focus is on the structural changes of bone development and growth, and modeling and repair in the four bone surfaces: periosteal, Haversian, endosteal, and trabecular. The pattern of skeletal changes is placed within the framework of the human life span. New perspectives and direction of research on the environmental, biological, and genetic influences on modeling and remodeling processes are discussed chronologically at each bone surface. Implications for biological anthropologists are considered. This approach emphasizes variation in skeletal biology as a dynamic record of development, maturity, and aging.

Genetic analysis of vertebral trabecular bone density and cross-sectional area in older men

We investigated 383 bone candidate genes for associations between single nucleotide polymorphisms and vertebral trabecular volumetric bone mineral density (vBMD) and cross-sectional area (CSA) in 2,018 Caucasian men aged ≥ 65 years. SNPs in TGFBR3, SOST, KL, CALCR, LEP, CSF1R, PTN, GNRH2, FGFR2, and MEPE were associated with vBMD and SNPs in CYP11B1, DVL2, DLX5, WNT4, and PAX7 were associated with CSA in independent study samples (p < 0.005).

INRODUCTION

Vertebral bone mineral density and cross-sectional area are important determinants of vertebral bone strength. Little is known about the specific genetic variants that influence these phenotypes in humans.

METHODS

We investigated the potential genetic variants associated with vertebral trabecular volumetric BMD and CSA measured by quantitative computed tomography. We initially tested for association between these phenotypes and 4608 tagging and potentially functional single nucleotide polymorphisms (SNPs) in 383 candidate genes in 862 community-dwelling Caucasian men aged ≥ 65 years in the Osteoporotic Fractures in Men Study.

RESULTS

SNP associations were then validated by genotyping an additional 1,156 randomly sampled men from the same cohort. We identified 11 SNPs in 10 genes (TGFBR3, SOST, KL, CALCR, LEP, CSF1R, PTN, GNRH2, FGFR2, and MEPE) that were consistently associated with trabecular vBMD and five SNPs in five genes (CYP11B1, DVL2, DLX5, WNT4, and PAX7) that were consistently associated with CSA in both samples (p < 0.005).

CONCLUSION

None of the SNPs associated with trabecular vBMD were associated with CSA. Our findings raise the possibility that at least some of the loci for vertebral trabecular BMD and bone size may be distinct.

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.

A comparison of principal component analysis and factor analysis strategies for uncovering pleiotropic factors

Principal component analysis (PCA) and factor analysis (FA) are often used to uncover genetic factors that contribute to complex disease phenotypes. The purpose of such an analysis is to distill a genetic signal from a large number of correlated phenotype measurements. That signal can then be used in genetic analyses (e.g. linkage analysis), presumably leading to greater success at finding genes than one would achieve with any one raw trait. Although both PCA and FA have been used this way, there has been no comparison of their performance in the literature. We compared the ability of these two procedures to extract unobserved underlying genetic components from complex simulated data on nuclear families. We first simulated seven underlying genetic and environmentally determined traits. Then we derived two sets of 50 complex (observed) traits using algebraic combinations of the underlying components. We next performed PCA and FA on the complex traits. We assessed two aspects of the performance of the methods: (1) ability to detect the underlying genetic components; (2) whether the methods worked better when applied to raw traits or to residuals (after regressing out significant environmental covariates). Our results indicate that both the methods behave similarly in most cases, although FA generally produced factors that had stronger correlations with the underlying traits. We also found that using residuals in PCA or FA analyses greatly increased the probability that the PCs or factors detected common genetic components instead of common environmental factors, except if there was statistical interaction between genetic and environmental factors.

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.

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.

Bone morphometry strongly predicts cortical bone stiffness and strength, but not toughness, in inbred mouse models of high and low bone mass

Inbred strains of mice make useful models to study bone properties. Our aim was to compare bone competence and cortical morphometric parameters of two inbred strains to better determine the role of bone structure and geometry in the process of bone failure. Morphometric analysis was performed on 20 murine femora with a low bone mass (C57BL/6J; B6) and 20 murine femora with a high bone mass (C3H/HeJ; C3H) using desktop microCT. The bones were tested under three-point bending to measure their mechanical properties. Results showed that the C3H strain is a more reproducible model regarding bone morphometric and mechanical phenotypes than the B6 strain. Bone strength, stiffness, yield force, yield displacement, and toughness, as well as morphometric traits, were all significantly different between the two strains, whereas postyield displacement was not. It was found that bone volume, cortical thickness, and cross-sectional area predicted almost 80% (p < 0.05) of bone stiffness, strength, and yield force. Nevertheless, cortical bone postyield properties such as bone toughness could not be explained by morphometry, but postyield whitening was observed in that phase. In conclusion, we found that morphometric parameters are strong predictors of preyield but not postyield properties. The lack of morphometric influence on bone competence in the postyield phase in combination with the observed postyield whitening confirmed the important contribution of ultrastructure and microdamage in the process of overall bone failure behavior, especially in the postyield phase.

Evidence for major pleiotropic effects on bone size variation from a principal component analysis of 451 Caucasian families

AIM

To identify pleiotropic quantitative trait loci (QTL) influencing bone size (BS) at different skeletal sites in Caucasians.

METHODS

In a sample containing 3899 Caucasians from 451 pedigrees, 410 microsatellite markers spaced approximately 8.9 cM apart across the human genome were genotyped. Phenotypical and genetic correlations of BS at lumbar spine, hip (femoral neck, trochanter, and intertrochanter regions), and wrist (ultradistal, mid-distal, and one-third distal sites) were determined using bivariate quantitative genetic analysis. A principal component analysis (PCA) was performed to obtain principal component (PC) factors that were then subjected to variance components linkage analysis to identify regions linked to the PC.

RESULTS

Genetic correlations of BS at different skeletal sites ranged from 0.40 to 0.79 (P<0.001). The PCA yielded a PC named PCtotal, which explained up to 76% of the total (co)variation of all the BS at the 7 skeletal sites for the whole sample. We identified a QTL influencing the BS of multiple skeletal sites on chromosome 7 at 140 cM [logarithm of odds (LOD)=2.85] in the overall sample. Sex-specific evidence for linkage was observed on chromosome 11 at 53 cM (LOD =2.82) in the male-only data subset.

CONCLUSION

Our study identified several genomic regions that may have pleiotropic effects on different skeletal sites. These regions may contain genes that play a critical role in overall bone development and osteoporosis at multiple skeletal sites, hence are biologically and clinically important.

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

INTRODUCTION

Study of mutations with large phenotypic effects has allowed the identification of key players in skeletal development. However, the molecular nature of variation in large, phenotypically normal populations tends to be characterized by smaller phenotypic effects that remain undefined.

MATERIALS AND METHODS

We use interval mapping and quantitative trait locus (QTL) mapping techniques in the combined F2-F3 populations (n = 2111) of an LG/J x SM/J mouse intercross to detect QTLs associated with the lengths of the humerus, ulna, femur, and tibia.

RESULTS

Seventy individual trait QTLs affecting long bone lengths were detected, with several chromosomes harboring multiple QTLs. The genetic architecture suggests mainly small, additive effects on long bone length, with roughly one third of the QTLs displaying dominance. Sex interactions were common, and four sex-specific QTLs were observed. Pleiotropy could not be rejected for most of the QTLs identified. Thirty-one epistatic interactions were detected, almost all affecting regions including or immediately adjacent to QTLs.

CONCLUSIONS

A complex regulatory network with many gene interactions modulates bone growth, possibly with integrated skeletal modules that allow fine-tuning of developmental processes present. Candidate genes in the QTL CIs include many genes known to affect endochondral bone growth and genes that have not yet been associated with bone growth or body size but have a strong potential to influence these traits.

Genetics of the musculoskeletal system: a pleiotropic approach

The risk of osteoporotic fracture can be viewed as a function of loading conditions and the ability of the bone to withstand the load. Skeletal loads are dominated by muscle action. Recently, it has become clear that bone and muscle share genetic determinants. Involution of the musculoskeletal system manifests as bone loss (osteoporosis) and muscle wasting (sarcopenia). Therefore, the consideration of pleiotropy is an important aspect in the study of the genetics of osteoporosis and sarcopenia. This Perspective will provide the evidence for a shared genetic influence on bone and muscle. We will start with an overview of accumulating evidence that physical exercise produces effects on the adult skeleton, seeking to unravel some of the contradictory findings published thus far. We will provide indications that there are pleiotropic relationships between bone structure/mass and muscle mass/function. Finally, we will offer some insights and practical recommendations as to the value of studying shared genetic factors and will explore possible directions for future research. We consider several related questions that together comprise the general paradigm of bone responses to mechanical loading and the relationship between muscle strength and bone parameters, including the genetic factors that modulate these responses. We believe that further progress in understanding the common genetic etiology of osteoporosis and sarcopenia will provide valuable insight into important biological underpinnings for both conditions and may translate into new approaches to reduce the burdens of both conditions through improved diagnosis, prevention, and early targeted treatment.

Contribution of gender-specific genetic factors to osteoporosis risk

Common diseases result from the complex relationship between genetic and environmental factors. The aim of this review is to provide perspective for a conceptual framework aimed at studying the interplay of gender-specific genetic and environmental factors in the etiology of complex disease, using osteoporosis as an example. In recent years, gender differences in the heritability of the osteoporosis-related phenotypes have been reported and sex-specific quantitative-trait loci were discovered by linkage studies in humans and mice. Results of numerous allelic association studies also differed by gender. In most cases, it was not clear whether or not this phenomenon should be attributed to the effect of sex-chromosomes, sex hormones, or other intrinsic or extrinsic differences between the genders, such as the level of bioavailable estrogen and of physical activity. We conclude that there is need to consider gender-specific genetic and environmental factors in the planning of future association studies on the etiology of osteoporosis and other complex diseases prevalent in the general population.

Mapping of the chromosome 17 BMD QTL in the F(2) male mice of MRL/MpJ x SJL/J

Developing treatment strategies for osteoporosis would be facilitated by identifying genes regulating bone mineral density (BMD). One way to do so is through quantitative trait locus (QTL) mapping. However, there are sex differences in terms of the presence/absence and locations of BMD QTLs. In a previous study, our group identified a BMD QTL on chromosome 17 in the F(2) female mice of the MRL/MpJ x SJL/J cross. Here, we determined whether it was also present in the male mice of the same cross. Furthermore, we also intended to reduce the QTL region by increasing marker density. Interval mapping showed that the same QTL based on chromosomal positions was present in the male mice, with logarithmic odds (LOD) scores of 4.0 for femur BMD and 5.2 for total body BMD. Although there was a body weight QTL at the same location, the BMD QTL was not affected by the adjustment for body weight. Mapping with increased marker density indicated a most likely region of 35-55 Mb for this QTL. There were also co-localized QTLs for femur length, femur periosteal circumference (PC) and total body bone area, suggesting possibility of pleiotropy. Runx2 and VEGFA are strong candidate genes located within this QTL region.

Pleiotropy and heterogeneity in the expression of bone strength-related phenotypes in extended pedigrees

Genetic analysis in 3,535 relative pairs from extended multigenerational families of African heritage showed that volumetric BMD is a highly heritable polygenic trait that is under compartment-specific genetic regulation. The majority of the phenotypic variation in bone size and volumetric BMD also seems to be strongly influenced by distinct genes for each trait.

INTRODUCTION

BMD and bone size contribute to bone strength and the risk of fracture. Little is known about the genetic architecture of QCT measures of volumetric BMD and bone size. We studied the contribution of genes, shared genes (pleiotropy), and shared environment to cortical and trabecular volumetric BMD and bone size using variance components analysis.

MATERIALS AND METHODS

A total of 471 individuals >or=18 yr of age (mean, 43 yr) from eight multigenerational Afro-Caribbean families (mean family size > 50; 3535 relative pairs) underwent a peripheral QCT scan of the radius and tibia and anthropometry.

RESULTS

Strong positive genetic correlations were observed for trabecular or cortical BMD measured at the tibia and radius (rho(G) > 0.82, p < 0.01), but not between trabecular and cortical BMD measured within the same anatomical site. Genetic correlations between volumetric BMD and bone length and circumference were also not statistically significant.

CONCLUSIONS

BMD is a highly heritable polygenic trait that is under compartment-specific genetic regulation. The majority of the phenotypic variation in skeletal size and density seems to be strongly influenced by distinct sets of genes for each trait.

Bivariate linkage study of proximal hip geometry and body size indices: the Framingham study

Femoral geometry and body size are both characterized by substantial heritability. The purpose of this study was to discern whether hip geometry and body size (height and body mass index, BMI) share quantitative trait loci (QTL). Dual-energy X-ray absorptiometric scans of the proximal femur from 1,473 members in 323 pedigrees (ages 31-96 years) from the Framingham Osteoporosis Study were studied. We measured femoral neck length, neck-shaft angle, subperiosteal width (outer diameter), cross-sectional bone area, and section modulus, at the narrowest section of the femoral neck (NN), intertrochanteric (IT), and femoral shaft (S) regions. In variance component analyses, genetic correlations (rho ( G )) between hip geometry traits and height ranged 0.30-0.59 and between hip geometry and BMI ranged 0.11-0.47. In a genomewide linkage scan with 636 markers, we obtained nominally suggestive linkages (bivariate LOD scores > or =1.9) for geometric traits and either height or BMI at several chromosomes (4, 6, 9, 15, and 21). Two loci, on chr. 2 (80 cM, BMI/shaft section modulus) and chr. X (height/shaft outer diameter), yielded bivariate LOD scores > or =3.0; although these loci were linked in univariate analyses with a geometric trait, neither was linked with either height or BMI. In conclusion, substantial genetic correlations were found between the femoral geometric traits, height and BMI. Linkage signals from bivariate linkage analyses of bone geometric indices and body size were similar to those obtained in univariate linkage analyses of femoral geometric traits, suggesting that most of the detected QTL primarily influence geometry of the hip.

Novel loci regulating bone anabolic response to loading: expression QTL analysis in C57BL/6JXC3H/HeJ mice cross

Variations in the expression levels of bone marker genes among the inbred strains of mice in response to mechanical loading (ML) are largely determined by genetic factors. To explore this, we performed four-point bending on tibiae of 10-week female F2 mice of B6XC3H cross using 9N at 2 Hz, 36 cycles, once per day for 12 days. We collected tibiae from these mice for RNA extraction. We then measured the expression changes of bone marker genes, bone sialoprotein (BSP), alkaline phosphatase (ALP) and housekeeping genes, beta-actin and peptidylprolyl isomerase A (PPIA), by using real-time PCR in both the loaded and the non-loaded tibiae of F2 mice (n=241). A genome-wide scan was performed using 111 micro satellite markers in DNA sample collected from these mice. Mean increase in gene expression, expressed as fold change, ranges from 2.8 to 3.0 for BSP and 2.7 to 2.8 for ALP. Both showed a skewed distribution with a heritability response of 87 to 91%. Absence of significant correlation between the increased gene expression vs. body weight (BW) and bone size (BS) suggests that bone response to loading is independent of BS or BW. Non-parametric mapping (MapQTL program 5) revealed that BSP and ALP expression in response to bending was regulated by several significant and suggestive QTL: Loci regulating both BSP and ALP were located on Chr 8 (60.1 cM), 16 (45.9 cM), 17 (14.2 cM), 18 (38.0 cM) and Chr 19 (3.3 cM); Loci specific to BSP were found on Chrs 1 (LOD score 10.4 at 91.8 cM), 5 (5.2 at 73.2 cM) and 9 (7.0 at 13.1 cM); Loci regulating only ALP were found on Chrs 1 (7.6 at 46 and 75.4 cM), 3 (8.3 at 47 cM) and 4 (5.6 at 54.6 cM). QTLs on Chrs 1, 3, 8, 9, 17 and 18 correspond to QTLs we previously reported by pQCT measurements, thus validating these findings. In addition, we found that the QTL associated with non-loaded tibiae for BSP and ALP on Chrs 4, 16 and 18 was identical to the QTLs associated with ML. This finding suggests that regions on these chromosomes are responsible for natural variation in expression of BSP and ALP as well as for ML. This is the first expression study to provide evidence for the presence of multiple genetic loci regulating bone anabolic response to loading in the B6XC3H intercross and will lead to a better understanding of how exercise improves the skeletal mass.

Quantitative trait locus that determines the cross-sectional shape of the femur in SAMP6 and SAMP2 mice

We segregated a QTL on chromosome 11 that affects femoral cross-sectional shape during growth by generating a congenic strain and an additional 16 subcongenic strains of the senescence-accelerated mouse strain, SAMP6. The QTL region was narrowed down to a 10.0-Mbp region.

INTRODUCTION

Genetic background is known to affect bone characteristics. However, little is known about how polymorphic genes modulate bone shape. In a previous study using SAMP2 and SAMP6 mice, we reported a quantitative trait locus (QTL) on chromosome (Chr) 11 that had significant linkage to peak relative bone mass in terms of cortical thickness index (CTI) in male mice. We named it Pbd1. Here we aimed to clarify the effects of Pbd1 on skeletal phenotype in male mice and to narrow down the QTL region.

MATERIALS AND METHODS

We generated a congenic strain named P6.P2-Pbd1(b), carrying a 39-cM SAMP2-derived Chr11 interval on a SAMP6 genetic background. Sixteen subcongenic strains with smaller overlapping intervals on the SAMP6 background were generated from P6.P2-Pbd1(b) to narrow the region of interest. The effects of Pbd1 on bone properties were determined. Gene expression analysis of all candidate genes in Pbd1 was performed using real-time RT-PCR.

RESULTS

The CTI of strain P6.P2-Pbd1(b) at 16 wk was higher than that of SAMP6. This was not caused by differences in cortical thickness but by cross-sectional shape. Morphological analysis by microCT revealed that the femoral cross-sectional shape of P6.P2-Pbd1(b) (and the other subcongenic strains with higher CTI or bone area fraction [BA/TA]) was more compressed anteroposteriorly than that of SAMP6, which was associated with superior mechanical properties. This feature was formed during bone modeling up to 16 wk of age. Subcongenic strains with a higher CTI showed significant increases in endocortical mineral apposition rate and significant reductions in periosteal mineral apposition rate at 8 wk compared with those of the SAMP6. The Pbd1 locus was successfully narrowed down to a 10.0-Mbp region, and the expression analysis suggested a candidate gene, Cacng4.

CONCLUSIONS

The Pbd1 affects femoral cross-sectional shape by regulating the rate of endocortical and periosteal bone formation of the femur during postnatal growth.

Quantitative trait loci for BMD and bone strength in an intercross between domestic and wildtype chickens

With chicken used as a model species, we used QTL analysis to examine the genetic contribution to bone traits. We report the identification of four QTLs for femoral traits: one for bone strength, one for endosteal circumference, and two affecting mineral density of noncortical bone.

INTRODUCTION

BMD is a highly heritable phenotype, governed by elements at numerous loci. In studies examining the genetic contribution to bone traits, many loci have been identified in humans and in other species. The goal of this study was to identify quantitative trait loci (QTLs) controlling BMD and bone strength in an intercross between wildtype and domestic chickens.

MATERIALS AND METHODS

A set of 164 markers, covering 30 chromosomes (chr.), were used to genotype 337 F2-individuals from an intercross of domesticated white Leghorn and wildtype red junglefowl chicken. DXA and pQCT were used to measure BMD and bone structure. Three-point bending tests and torsional strength tests were performed to determine the biomechanical strength of the bone. QTLs were mapped using forward selection for loci with significant marginal effects.

RESULTS

Four QTLs for femoral bone traits were identified in QTL analysis with body weight included as a covariate. A QTL on chr. 1 affected female noncortical BMD (LOD 4.6) and is syntenic to human 12q21-12q23. Also located on chr. 1, a locus with synteny to human 12q13-14 affected endosteal circumference (LOD 4.6). On chr. 2, a QTL corresponding to human 5p13-p15, 7p12, 18q12, 18q21, and 9q22-9q31 affected BMD in females; noncortical (LOD 4.0) and metaphyseal (LOD 7.0) BMD by pQCT and BMD by DXA (LOD 5.9). A QTL located on chr. 20 (LOD 5.2) affected bone biomechanical strength and had sex-dependent effects. In addition to the significant QTLs, 10 further loci with suggestive linkage to bone traits were identified.

CONCLUSIONS

Four QTLs were identified: two for noncortical BMD, one for endosteal circumference, and one affecting bone biomechanical strength. The future identification of genes responsible for these QTLs will increase the understanding of vertebrate skeletal biology.

Detecting novel bone density and bone size quantitative trait loci using a cross of MRL/MpJ and CAST/EiJ inbred mice

Most previous studies to identify loci involved in bone mineral density (BMD) regulation have used inbred strains with high and low BMD in generating F(2) mice. However, differences in BMD may not be a requirement in selecting parental strains for BMD quantitative trait loci (QTL) studies. In this study, we intended to identify novel QTL using a cross of two strains, MRL/MpJ (MRL) and CAST/EiJ (CAST), both of which exhibit relatively high BMD when compared to previously used strains. In addition, CAST was genetically distinct. We generated 328 MRL x CAST F(2) mice of both sexes and measured femur BMD and periosteal circumference (PC) using peripheral quantitative computed tomography. Whole-genome genotyping was performed with 86 microsatellite markers. A new BMD QTL on chromosome 10 and another suggestive one on chromosome 15 were identified. A significant femur PC QTL identified on chromosome 9 and a suggestive one on chromosome 2 were similar to those detected in MRL x SJL. QTL were also identified for other femur and forearm bone density and bone size phenotypes, some of which were colocalized within the same chromosomal positions as those for femur BMD and femur PC. This study demonstrates the utility of crosses involving inbred strains of mice which exhibit a similar phenotype in QTL identification.

Molecular genetic studies of gene identification for osteoporosis: a 2004 update

This review summarizes comprehensively the most important and representative molecular genetics studies of gene identification for osteoporosis published up to the end of December 2004. It is intended to constitute a sequential update of our previously published review covering the available data up to the end of 2002. Evidence from candidate gene association studies and genome-wide linkage studies in humans, as well as quantitative trait locus mapping animal models are reviewed separately. Studies of transgenic and knockout mice models relevant to osteoporosis are summarized. An important extension of this update is incorporation of functional genomic studies (including DNA microarrays and proteomics) on osteogenesis and osteoporosis, in light of the rapid advances and the promising prospects of the field. Comments are made on the most notable findings and representative studies for their potential influence and implications on our present understanding of genetics of osteoporosis. The format adopted by this review should be ideal for accommodating future new advances and studies.

Identification of genetic loci that regulate bone adaptive response to mechanical loading in C57BL/6J and C3H/HeJ mice intercross

Strain-dependent differences in bone adaptive responses to loading among inbred mouse strains suggest that genetic background contributes significantly to adaptation to exercise. To explore the genetic regulation of response to loading, we performed a genome-wide search for linkage in a cross between two strains, a good responder, C57BL6/J (B6), and a poor responder, C3H/HeJ (C3H). Using a four-point bending model, the right tibia was loaded by applying 9 N force for 36 cycles for 12 days in 10-week-old female B6xC3H F2 mice. Changes in bone density (BMD) and bone size were evaluated in vivo by pQCT. Measurements from non-loaded left tibia were used as an internal control to calculate loading-induced percent increase in BMD and bone size, thus excluding the possibility of identifying background QTL(s) due to natural allelic variation in mapping strains. A genome-wide scan was performed using 111 microsatellite markers in DNA samples collected from 329 F2 mice. Heritability of bone adaptive response to loading was between 70 and 80%. The mean increase, expressed as percent of unloaded tibia, was 5% for BMD, 9% for periosteal circumference (PC), and 14% for cortical thickness in F2 mice (n = 329). All these phenotypes showed normal distributions. Absence of significant correlation between BMD response to four-point bending and body weight or bone size suggested that the bone adaptive response was independent of bone size. Interval mapping revealed that BMD response to four-point bending was influenced by three significant loci on Chrs 1 (log-of-odds ratio score (LOD) 3.4, 91.8 cM), 3 (LOD 3.6, 50.3 cM), and 8 (LOD 4.2, 60.1 cM) and one suggestive QTL on Chr 9 (LOD 2.5, 33.9 cM). Loading-induced increases in PC and Cth were influenced by four significant loci on Chrs 8 (LOD 3.0, 68.9 cM), 9 (LOD 3.0, 13.1 cM), 17 (LOD 3.0, 39.3 cM), and 18 (LOD 3.0, 0 cM) and two suggestive loci on Chr 9 (LOD 2.2, 24 cM) and 11 (LOD 2.1, 69.9 cM). Pairwise analysis showed the presence of several significant and suggestive interactions between loci on Chrs 1, 3, 8, and 13 for BMD trait. This is the first study that provides evidence for the presence of multiple genetic loci regulating bone anabolic responses to loading in the B6xC3H intercross. Knowledge of the genes underlying these loci could provide novel approaches to improve skeletal mass.

QTL analysis of trabecular bone in BXD F2 and RI mice

A sample of 693 mice was used to identify regions of the mouse genome associated with trabecular bone architecture as measured using microCT. QTLs for bone in the proximal tibial metaphysis were identified on several chromosomes indicating regions containing genes that regulate properties of trabecular bone.

INTRODUCTION

Age-related osteoporosis is a condition of major concern because of the morbidity and mortality associated with osteoporotic fractures in humans. Osteoporosis is characterized by reduced bone density, strength, and altered trabecular architecture, all of which are quantitative traits resulting from the actions of many genes working in concert with each other and the environment over the lifespan. microCT gives accurate measures of trabecular bone architecture providing phenotypic data related to bone volume and trabecular morphology. The primary objective of this research was to identify chromosomal regions called quantitative trait loci (QTLs) that contain genes influencing trabecular architecture as measured by microCT.

MATERIALS AND METHODS

The study used crosses between C57BL/6J (B6) and DBA/2J (D2) as progenitor strains of a second filial (F2) generation (n = 141 males and 148 females) and 23 BXD recombinant inbred (RI) strains (n approximately 9 of each sex per strain). The proximal tibial metaphyses of the 200-day-old mice were analyzed by microCT to assess phenotypic traits characterizing trabecular bone, including bone volume fraction, trabecular connectivity, and quantitative measures of trabecular orientation and anisotropy. Heritabilities were calculated and QTLs were identified using composite interval mapping.

RESULTS

A number of phenotypes were found to be highly heritable. Heritability values for measured phenotypes using RI strains ranged from 0.15 for degree of anisotropy in females to 0.51 for connectivity density in females and total volume in males. Significant and confirmed QTLs, with LOD scores 4.3 in the F2 cohort and 1.5 in the corresponding RI cohort were found on chromosomes 1 (43 cM), 5 (44 cM), 6 (20 cM), and 8 (49 cM). Other QTLs with LOD scores ranging from 2.8 to 6.9 in the F2 analyses were found on chromosomes 1, 5, 6, 8, 9, and 12. QTLs were identified using data sets comprised of both male and female quantitative traits, suggesting similar genetic action in both sexes, whereas others seemed to be associated exclusively with one sex or the other, suggesting the possibility of sex-dependent effects.

CONCLUSIONS

Identification of the genes underlying these QTLs may lead to improvements in recognizing individuals most at risk for developing osteoporosis and in the design of new therapeutic interventions.

Femur mechanical properties in the F2 progeny of an NZB/B1NJ x RF/J cross are regulated predominantly by genetic loci that regulate bone geometry

Genetic analysis of an NZB/B1NJ x RF/J cross has identified QTLs for femur mechanical, geometric, and densitometric phenotypes. Most mechanical QTLs were associated with geometric QTLs, strongly suggesting common genetic regulation.

INTRODUCTION

Previous studies have shown that bone architecture and BMD are important factors affecting bone strength, and both are genetically regulated. We conducted genetic analyses for loci regulating femur mechanical properties, geometric properties, and BMD in a cohort of F2 mice derived from intercross matings of (NZB/B1NJ x RF/J)F1 parents.

MATERIALS AND METHODS

Femurs were isolated from 662 10-week-old females. Mechanical properties were determined for a femur from each animal by three-point bending. Geometric properties and volumetric BMD (vBMD) were determined by pQCT. Genotype data were obtained by PCR assays for polymorphic markers carried in the genomic DNA of each mouse. Genome-wide scans were carried out for co-segregation of genetic marker data with values from 23 different phenotypes. Quantitative trait loci (QTLs) were identified for mechanical, geometric, and mineral density phenotypes.

RESULTS

QTLs for many phenotypes were significantly refined by covariate analyses using body weight and femur length. Major QTLs for mechanical and geometric phenotypes were found on chromosomes 5, 7, 9, 11, and 12. Nine chromosomal locations were identified with mechanical QTLs and 17 locations with one or more geometric QTLs. The significance of five mechanical and nine geometric QTLs was affected by the inclusion of covariates. These changes included both decreases and increases in significance. The QTLs on chromosomes 5 and 12 were decreased by inclusion of the covariates in the analysis, but QTLs on 7 and 11 were unaffected. Mechanical QTLs were almost always associated with geometric QTLs and less commonly (two of six) with vBMD QTLs.

CONCLUSIONS

Genetic regulation of mechanical properties in the F(2) mice of this NZB/B1NJ x RF/J cross seems to be caused by genes regulating femur geometry.

Genetics of canid skeletal variation: size and shape of the pelvis

The mammalian skeleton presents an ideal system in which to study the genetic architecture of a set of related polygenic traits and the skeleton of the domestic dog (Canis familiaris) is arguably the best system in which to address the relationship between genes and anatomy. We have analyzed the genetic basis for skeletal variation in a population of >450 Portuguese Water Dogs. At this stage of this ongoing project, we have identified >40 putative quantitative trait loci (QTLs) for heritable skeletal phenotypes located on 22 different chromosomes, including the "X." A striking aspect of these is the regulation of suites of traits representing bones located in different parts of the skeleton but related by function. Here we illustrate this by describing genetic variation in postcranial morphology. Two suites of traits are involved. One regulates the size of the pelvis relative to dimensions of the limb bones. The other regulates the shape of the pelvis. Both are examples of trade-offs that may be prototypical of different breeds. For the size of the pelvis relative to limb bones, we describe four QTLs located on autosome CFA 12, 30, 31, and X. For pelvic shape we describe QTLs on autosome CFA 2, 3, 22, and 36. The relation of these polygenic systems to musculoskeletal function is discussed.

Identification of quantitative trait loci that regulate obesity and serum lipid levels in MRL/MpJ x SJL/J inbred mice

The total body fat mass and serum concentration of total cholesterol, HDL cholesterol, and triglyceride (TG) differ between standard diet-fed female inbred mouse strains MRL/MpJ (MRL) and SJL/J (SJL) by 38-120% (P < 0.01). To investigate genetic regulation of obesity and serum lipid levels, we performed a genome-wide linkage analysis in 621 MRLx SJL F2 female mice. Fat mass was affected by two significant loci, D11Mit36 [43.7 cM, logarithm of the odds ratio (LOD) 11.2] and D16Mit51 (50.3 cM, LOD 3.9), and one suggestive locus at D7Mit44 (50 cM, LOD 2.4). TG levels were affected by two novel loci at D1Mit43 (76 cM, LOD 3.8) and D12Mit201 (26 cM, LOD 4.1), and two suggestive loci on chromosomes 5 and 17. HDL and cholesterol concentrations were influenced by significant loci on chromosomes 1, 3, 5, 7, and 17 that were in the regions identified earlier for other strains of mice, except for a suggestive locus on chromosome 14 that was specific to the MRL x SJL cross. In summary, linkage analysis in MRL x SJL F2 mice disclosed novel loci affecting TG, HDL, and fat mass, a measure of obesity. Knowledge of the genes in these quantitative trait loci will enhance our understanding of obesity and lipid metabolism.

Genetic epidemiology of osteoporosis: past, present, and future

Research during the past several decades has unequivocally established a role of heredity in the etiology of osteoporosis. Major efforts are currently underway to identify the genes and allelic variants that confer genetic susceptibility to this common and disabling condition. Genome-wide linkage mapping in families, candidate gene association studies in unrelated individuals, and quantitative trait locus mapping in animal models are the primary strategies being used to search for the genetic contributors to osteoporosis. Genome-wide mapping efforts have identified the low-density lipoprotein receptor-related protein 5, bone morphogenetic protein 2, and 15-lipoxygenase as potential susceptibility genes for osteoporosis in the past few years, providing a rich new base for understanding bone biology. Candidate gene association analyses have also provided evidence for a modest role of allelic variants in several additional genes including collagen type Ialpha1, vitamin D receptor, and estrogen receptor-alpha. With the development of a high-density genome-wide polymorphism and haplotype map and continued improvements in high-throughput and cost-effective genotyping technologies, many more genetic contributors to osteoporosis will probably be identified in the near future. The results of this research should facilitate the development of new methods for diagnosing, preventing, and treating the growing clinical and public health problem of osteoporosis.

Identification of novel genetic loci for bone size and mechanosensitivity in an ENU mutant exhibiting decreased bone size

Using a dominant ENU mutagenesis screen in C57BL/6J (B6) mice to reveal gene function, we identified a mutant, 917M, with a reduced bone size phenotype, which is expressed only in males. We show that mutation results in osteoblasts with reduced proliferation, increased apoptosis, and an impaired response to in vitro mechanical load. The mutation is mapped to a novel locus (LOD score of 7.9 at 10.5 cM) on chromosome 4.

INTRODUCTION

Using a dominant ENU mutagenesis screen in C57BL/6J (B6) mice to reveal gene function, we identified a mutant, 917M, with a reduced bone size phenotype, which is expressed only in males. In this report, we show the chromosomal location of this mutation using linkage analysis and cellular characterization of the mutant phenotype.

MATERIALS AND METHODS

The mutant mouse was bred to wildtype B6 to produce progeny for characterization of the bone size phenotype. Periosteal osteoblasts isolated from the tibia and femur of mutant and wildtype mice were studied for proliferation, differentiation, and apoptosis potential. To determine the chromosomal location of the mutation, a low-resolution linkage map was established by completing a genome-wide scan in B6C3H F2 male mice generated from intercross breeding of mutant mice.

RESULTS AND CONCLUSIONS

Mutant progeny (16 weeks old) displayed a total body bone area that was 10-13% lower and a periosteal circumference that was 5-8% lower at the femur and tibia midshaft compared with wildtype B6 mice. Periosteal osteoblasts from mutant mice showed 17-27% reduced cell proliferation and 23% increased apoptosis compared with wildtype controls. In addition, osteoblasts from mutant mice showed an impaired response to shear stress-induced proliferation rate, an in vitro model for mechanical loading. Interval mapping in B6C3H F2 males (n = 69) indicated two major loci affecting bone size on chromosome 1 at 45 cM (LOD 4.9) and chromosome 4 at 10.5 cM (LOD 7.9, genome-wide p < 0.01). Interval mapping using body weight as covariate revealed only one significant interval at chromosome 4 (LOD 6.8). Alleles of the chromosome 4 interval inherited from the B6 mutant strain contributed to a significantly lower bone size than those inherited from C3H. A pairwise interaction analysis showed evidence for a significant interaction between loci on chromosome 1 with the chromosome 4 quantitative trait loci. The 917M locus on chromosome 4 seems to be novel because it does not correspond with those loci previously associated with bone size on chromosome 4 in B6 and C3H/HeJ mice or other crosses.

A gene(s) for all-trans-retinoic acid-induced forelimb defects mapped and confirmed to murine chromosome 11

All-trans-retinoic acid (RA) induces various anatomical limb dysmorphologies in mice dependent on the time of exposure. During early limb development, RA induces forelimb ectrodactyly (digital absence) with varying susceptibilities for different inbred mouse strains; C57BL/6N are highly susceptible while SWV are resistant. To isolate the genetic basis of this defect, a full-genome scan was performed in 406 backcross fetuses of F(1) males to C57BL/6N females. Fetuses were exposed via a maternal injection of 75 mg of RA per kilogram of body weight on gestational day 9.25. The genome-wide analysis revealed significant linkage to a chromosome 11 locus near D11Mit39 with a maximum LOD score of 9.0 and to a chromosome 4 locus near D4Mit170. An epistatic interaction was detected between loci on chromosome 11 (D11Mit39) and chromosome 18 (D18Mit64). Linkage to the chromosome 11 locus (D11Mit39) was confirmed in RA-treated backcross fetuses of F(1) females to C57BL/6N males. Loci associated with bone density/mass in both human and mouse were previously detected in the same region, suggesting a mechanistic linkage with bone homeostasis. The human syntenic region of this locus has been previously linked to Meckel syndrome; the phenotype includes postaxial polydactyly, an ectopic digital defect hypothesized to be induced by a common molecular pathway with ectrodactyly.

Quantitative trait loci analysis of structural and material skeletal phenotypes in C57BL/6J and DBA/2 second-generation and recombinant inbred mice

QTL analyses identified several chromosomal regions influencing skeletal phenotypes of the femur and tibia in BXD F2 and BXD RI populations of mice. QTLs for skeletal traits co-located with each other and with correlated traits such as body weight and length, adipose mass, and serum alkaline phosphatase.

INTRODUCTION

Past research has shown substantial genetic influence on bone quality, and the impact of reduced bone mass on our aging population has heightened the interest in skeletal genetic research.

MATERIALS AND METHODS

Quantitative trait loci (QTL) analyses were performed on morphologic measures and structural and material properties of the femur and tibia in 200-day-old C57BL/6J x DBA/2 (BXD) F2 (second filial generation; n = 400) and BXD recombinant inbred (RI; n = 23 strains) populations of mice. Body weight, body length, adipose mass, and serum alkaline phosphatase were correlated phenotypes included in the analyses.

RESULTS

Skeletal QTLs for morphologic bone measures such as length, width, cortical thickness, and cross-sectional area mapped to nearly every chromosome. QTLs for both structural properties (ultimate load, yield load, or stiffness) and material properties (stress and straincharacteristics and elastic modulus) mapped to chromosomes 4, 6, 9, 12, 13, 15, and 18. QTLs that were specific to structural properties were identified on chromosomes 1, 2, 3, 7, 8, and 17, and QTLs that were specific to skeletal material properties were identified on chromosomes 5, 11, 16, and 19. QTLs for body size (body weight, body length, and adipose mass) often mapped to the same chromosomal regions as those identified for skeletal traits, suggesting that several QTLs identified as influencing bone could be mediated through body size.

CONCLUSION

New QTLs, not previously reported in the literature, were identified for structural and material properties and morphological measures of the mouse femur and tibia. Body weight and length, adipose mass, and serum alkaline phosphatase were correlated phenotypes that mapped in close proximity of skeletal chromosomal loci. The more specific measures of bone quality included in this investigation enhance our understanding of the functional significance of previously identified QTLs.

Genetic variation in femur extrinsic strength in 29 different inbred strains of mice is dependent on variations in femur cross-sectional geometry and bone density

The femurs from groups of mice from 29 different inbred strains were characterized to study the genetic variations in bone parameters. For these analyses, we used peripheral quantitative computed tomography to assess bone size and density in addition to three-point bend testing to assess bone mechanical properties. Highly significant differences between inbred strains were found for all size, density, and mechanical parameters measured (P < 0.0001). Correcting femoral cross-sectional geometry values or bone mechanical properties values for body weight or femur length reduced but did not eliminate the variations in bone geometry or bone mechanical properties. Mice of similar body size had as much as a 40% difference in the midshaft total area of the femur. Regression analysis suggested that 50.9% of the variation in maximum load among strains was related to variations in section modulus, i.e., cross-sectional geometry, 21.5% was related to variations in material bone density, and 27.7% to variations in quality. These components were further analyzed to show that 3.9-27.8% of the variation in maximum load was related to adaptation to mechanical stress. These findings indicate that there is a significant genetic variation in the femur cross-sectional area, density, and mechanical properties between inbred mouse strains. These studies identify inbred mouse strains suitable for future studies identifying genes regulating bone geometry and mechanical properties.

Periosteum: biology, regulation, and response to osteoporosis therapies

Periosteum contains osteogenic cells that regulate the outer shape of bone and work in coordination with inner cortical endosteum to regulate cortical thickness and the size and position of a bone in space. Induction of periosteal expansion, especially at sites such as the lumbar spine and femoral neck, reduces fracture risk by modifying bone dimensions to increase bone strength. The cell and molecular mechanisms that selectively and specifically activate periosteal expansion, as well as the mechanisms by which osteoporosis drugs regulate periosteum, remain poorly understood. We speculate that an alternate strategy to protect human bones from fracture may be through targeting of the periosteum, either using current or novel agents. In this review, we highlight current concepts of periosteal cell biology, including their apparent differences from endosteal osteogenic cells, discuss the limited data regarding how the periosteal surface is regulated by currently approved osteoporosis drugs, and suggest one potential means through which targeting periosteum may be achieved. Improving our understanding of mechanisms controlling periosteal expansion will likely provide insights necessary to enhance current and develop novel interventions to further reduce the risk of osteoporotic fractures.

Mapping quantitative trait loci that influence blood levels of alkaline phosphatase in MRL/MpJ and SJL/J mice

To examine the hypothesis that serum alkaline phosphatase (ALP) levels have a heritable component, we analyzed blood from two inbred strains of mice, MRL/MpJ and SJL, which exhibit 90% difference in total serum ALP activity (268+/-26 vs. 140+/-15 U/l, respectively, P<0.001). A genome-wide scan was carried out using 137 polymorphic markers in 518 F2 female mice. Serum ALP activity in the F2 progeny showed a normal distribution with an estimated heritability of 56%. Genome-wide scan for cosegregation of genetic marker data with serum ALP activity revealed three major quantitative trait loci (QTL), one each on chromosomes 2 (LOD score 3.8), chromosome 6 (LOD score 12.0), and chromosome 14 (LOD score 3.7). In addition, there was one suggestive QTL on chromosome 2 (LOD score of 3.3). In aggregate, these QTLs explain 22.5% of variance in serum ALP between these two strains. Serum ALP showed a moderate but significant correlation with body weight adjusted total body bone mineral density (r=0.12, P=0.0108) and periosteal circumference at midshaft tibia (r=0.15, P=0.0006) in F2 mice. The chromosome 6 locus harboring the major serum ALP QTL also contains a major BMD and bone size QTL, identified earlier, between these two strains of mice; in addition, this QTL is also close to the locus that regulates IGF-I levels (LOD score 8-9) in C3HB6 F2 mice. These common QTLs indicate that the observed difference in ALP and BMD or bone size may be regulated by same loci (or genes). Accordingly, the osteoblast cells isolated from femur and tibia of MRL mice showed a significantly higher number of ALP +ve cells/colony and two- to threefold higher ALP activity (P<0.001) as compared to the cells isolated from SJL mice, thus suggesting that differences in serum ALP between MRL and SJL reflect difference in ALP expression from osteoblasts from these strains of mice. These data suggest that serum ALP levels are genetically determined and correlate with cellular mechanisms that differentiate BMD accrual in these two strains of mice. The findings that ALP and BMD traits share the same loci on chromosome 6 suggest a role for genetic determinants of bone formation in overall BMD accretion.

HSV-1 amplicon-mediated transfer of 128-kb BMP-2 genomic locus stimulates osteoblast differentiation in vitro

In previous studies, we developed mouse genetic models and discovered genetic components of quantitative trait loci on mouse chromosomes that contribute to phenotypes such as bone size, bone density, and fracture healing. However, these regions contain dozens of genes in several overlapping bacterial artificial chromosomes (BACs) and are difficult to clone by physical cloning strategies. A feasible and efficient approach of identifying candidate genes is to transfer the genomic loci in BAC clones into mammalian cells for functional studies. In this study, we retrofitted a BAC construct into herpes simplex virus-1 amplicon and packaged it into an infectious BAC (iBAC) to test gene function in a cell-based system, using a 128-kb clone containing the complete bone morphogenetic protein-2 (BMP-2) gene. We transduced MC3T3-E1 cells with the iBAC bearing BMP-2 gene and examined transgene expression and function. Our results have demonstrated that an iBAC can efficiently deliver a BMP-2 genomic locus into preosteoblast cells and express functional BMP-2 protein, inducing a phenotype of cell differentiation, as indicated by an increase in alkaline phosphatase activity. Therefore, this experimental system provides a rapid, efficient cell-based model of high-throughput phenotypic screening to identify the BAC clones from physically mapped regions that are important for osteoblast differentiation. It also illustrates the potential of iBAC technology in functional testing of single nucleotide polymorphisms located in the distal promoter or/and intron regions responsible for low bone density.