Mapping quantitative trait loci for vertebral trabecular bone volume fraction and microarchitecture in mice

Semi-Automated MicroCT Analysis of Bone Anatomy and Mineralization in Mouse Models

The skeletal system mirrors several processes in the vertebrate body that impact developmental malfunctions, hormonal disbalance, malfunction of calcium metabolism and turn over, and inflammation processes such as arthrosis. X-ray micro computed tomography is a useful tool for 3D in situ evaluation of the skeletal system in a time-related manner, but results depend highly on resolution. Here, we provide the methodological background for a graduated evaluation from whole-body analysis of skeletal morphology and mineralization to high-resolution analysis of femoral and vertebral microstructure. We combine an expert-based evaluation with a machine-learning-based computational approach, including pre-setup analytical task lists. © 2024 Wiley Periodicals LLC. Basic Protocol 1: In vivo microCT scanning and skeletal analysis in mice Basic Protocol 2: Ex vivo high-resolution microCT scanning and microstructural analysis of the femur and L4 vertebra.

Assessing patterns of variation in BV/TV in the calcaneus and C2 vertebra of Gorilla gorilla, Pan troglodytes, and populations of Homo sapiens from the Pleistocene and Holocene that differ in physical activity levels

OBJECTIVES

Because trabecular bone volume fraction (BV/TV) is influenced by variations in physical activity recent declines in BV/TV in humans are often attributed to modern sedentary lifestyles. This study tests the hypothesis that presumed variations in mechanical loading between groups can predict the observed BV/TV patterns in humans, chimpanzees and gorillas in two bones: the calcaneus which experiences high and well characterized impact forces, and the C2 vertebrae which experiences reduced locomotor forces.

MATERIALS AND METHODS

BV/TV and other structural variables were quantified from high-resolution microCT scans in gorillas, chimpanzees, and four Homo sapiens populations: Pleistocene, semi-sedentary Natufians; Holocene hunter-gatherers from Point Hope, Alaska; Holocene nomadic pastoralists from medieval Europe; and modern, sedentary Americans.

RESULTS

In the calcaneal tuberosity, Natufian BV/TV was 36, 46, and 46% greater than Alaskans (p = .02), Europeans (p = .005) and modern Americans (p = .002), respectively, but not significantly different from apes. BV/TV was not significantly different between modern Americans and Alaskans or Europeans. In the C2, Natufian BV/TV was 53 and 25% greater than in the Alaskan (p = .0001) and European (p = .048) populations.

DISCUSSION

These results suggest that phenomena other than or in addition to variations in physical activity are needed to explain BV/TV patterns observed in H. sapiens, and point to a systemic decline in H. sapiens BV/TV after the Pleistocene.

Genetic Dissection of Femoral and Tibial Microarchitecture

Our understanding of the genetic control of bone strength has relied mainly on estimates of bone mineral density. Here we have mapped genetic factors that influence femoral and tibial microarchitecture using high‐resolution x‐ray computed tomography (8‐μm isotropic voxels) across a family of 61 BXD strains of mice, roughly 10 isogenic cases per strain and balanced by sex. We computed heritabilities for 25 cortical and trabecular traits. Males and females have well‐matched heritabilities, ranging from 0.25 to 0.75. We mapped 16 genetic loci most of which were detected only in females. There is also a bias in favor of loci that control cortical rather than trabecular bone. To evaluate candidate genes, we combined well‐established gene ontologies with bone transcriptome data to compute bone‐enrichment scores for all protein‐coding genes. We aligned candidates with those of human genome‐wide association studies. A subset of 50 strong candidates fell into three categories: (1) experimentally validated genes already known to modulate bone function (Adamts4, Ddr2, Darc, Adam12, Fkbp10, E2f6, Adam17, Grem2, Ifi204); (2) candidates without any experimentally validated function in bone (eg, Greb1, Ifi202b), but linked to skeletal phenotypes in human cohorts; and (3) candidates that have high bone‐enrichment scores, but for which there is not yet any functional link to bone biology or skeletal system disease (including Ifi202b, Ly9, Ifi205, Mgmt, F2rl1, Iqgap2). Our results highlight contrasting genetic architecture between sexes and among major bone compartments. The alignment of murine and human data facilitates function analysis and should prove of value for preclinical testing of molecular control of bone structure. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

Dietary calcium intake and genetics have site-specific effects on peak trabecular bone mass and microarchitecture in male mice

Trabecular bone (Tb) is used for rapid exchange of calcium (Ca) in times of physiologic need and the site-specific characteristics of Tb may explain why certain sites are more vulnerable to osteoporosis. We hypothesized that peak trabecular bone mass (PTBM) and Tb microarchitecture are differentially regulated by dietary Ca intake, genetics, or Gene-by-Diet (GxD) interactions at the distal femur and the fifth lumbar (L5) vertebra. Male mice from 62 genetically distinct lines were fed basal (0.5%) or low (0.25%) Ca diets from 4 to 12 wks of age. Afterwards, the right femur and L5 vertebra were removed and trabecular bone was analyzed by μCT. In mice fed the basal diet, bone volume fraction (BV/TV), trabecular number (Tb.N), and connectivity density (Conn.D) were significantly higher in the L5 vertebra than femur. Femur Tb had a weaker, more rod-like structure than the L5 vertebrae while mice fed the low Ca diet developed rod-like structures at both sites. Dietary Ca restriction also caused a greater relative reduction of Tb.N and Conn.D in the femur than L5 vertebra, i.e. it was more harmful to the integrity of Tb microarchitecture in femur. Genetics was a major determinant of Tb at both sites, e.g. heritability of BV/TV on the basal diet = 0.65 (femur) and 0.68 (L5 vertebra). However, while GxD interactions altered the impact of dietary Ca restriction on Tb parameters at both sites, the effect was not uniform, e.g. some lines had site-specific responses to Ca restriction. The significance of our work is that there are site-specific effects of dietary Ca restriction and genetics that work independently and interactively to influence the attainment of PTBM and Tb microarchitecture.

Abnormal Thyroid Hormone Status Differentially Affects Bone Mass Accrual and Bone Strength in C3H/HeJ Mice: A Mouse Model of Type I Deiodinase Deficiency

C3H/HeJ (C3H) mice are deficient of type I deiodinase (D1), an enzyme that activates thyroid hormone (TH), converting thyroxine (T4) to triiodothyronine (T3). Nevertheless, C3H mice present normal serum T3 and a gross euthyroid phenotype. To investigate if a global D1 deficiency interferes in the TH effects on bone, we compared bone growth, bone mass accrual and bone strength of C3H and C57BL/6J (B6) mice under abnormal TH status. Four-week-old female mice of both strains were grouped as Euthyroid, Hypothyroid (pharmacologically-induced), 1xT4 and 10xT4 (hypothyroid animals receiving 1- or 10-fold the physiological dose of T4 /day/16 weeks). Hypothyroidism and TH excess similarly impaired body weight (BW) gain and body growth in both mice strains. In contrast, whereas hypothyroidism only slightly impaired bone mineral density (BMD) accrual in B6 mice, it severely impaired BMD accrual in C3H mice. No differences were observed in serum and bone concentrations of T3 between hypothyroid animals of both strains. Interestingly, treatment with 10xT4 was less deleterious to BMD accrual in C3H than in B6 mice and resulted in less elevated T3 serum levels in B6 than in C3H mice, which is probably explained by the lower D1 activity in C3H mice. In addition, hypothyroidism decreased bone strength only in C3H but not in B6 mice, while TH excess decreased this parameter in both strains. These findings indicate that D1 deficiency contributes to the TH excess-induced differences in bone mass accrual in C3H vs. B6 mice and suggest that deiodinase-unrelated genetic factors might account for the different skeleton responses to hypothyroidism between strains.

Mouse Models and Online Resources for Functional Analysis of Osteoporosis Genome-Wide Association Studies

Osteoporosis is a complex genetic disease in which the number of loci associated with the bone mineral density, a clinical risk factor for fracture, has increased at an exponential rate in the last decade. The identification of the causative variants and candidate genes underlying these loci has not been able to keep pace with the rate of locus discovery. A large number of tools and data resources have been built around the use of the mouse as model of human genetic disease. Herein, we describe resources available for functional validation of human Genome Wide Association Study (GWAS) loci using mouse models. We specifically focus on large-scale phenotyping efforts focused on bone relevant phenotypes and repositories of genotype-phenotype data that exist for transgenic and mutant mice, which can be readily mined as a first step toward more targeted efforts designed to deeply characterize the role of a gene in bone biology.

Screening Gene Knockout Mice for Variation in Bone Mass: Analysis by μCT and Histomorphometry

PURPOSE OF REVIEW

The international mouse phenotyping consortium (IMPC) is producing defined gene knockout mouse lines. Here, a phenotyping program is presented that is based on micro-computed tomography (μCT) assessment of distal femur and vertebra. Lines with significant variation undergo a computer-based bone histomorphometric analysis.

RECENT FINDINGS

Of the 220 lines examined to date, approximately 15% have a significant variation (high or low) by μCT, most of which are not identified by the IMPC screen. Significant dimorphism between the sexes and bone compartments adds to the complexity of the skeletal findings. The μCT information that is posted at www.bonebase.org can group KOMP lines with similar morphological features. The histological data is presented in a graphic form that associates the cellular features with a specific anatomic group. The web portal presents a bone-centric view appropriate for the skeletal biologist/clinician to organize and understand the large number of genes that can influence skeletal health. Cataloging the relative severity of each variant is the first step towards compiling the dataset necessary to appreciate the full polygenic basis of degenerative bone disease.

Genetic Dissection of Trabecular Bone Structure with Mouse Intersubspecific Consomic Strains

Trabecular bone structure has an important influence on bone strength, but little is known about its genetic regulation. To elucidate the genetic factor(s) regulating trabecular bone structure, we compared the trabecular bone structures of two genetically remote mouse strains, C57BL/6J and Japanese wild mouse-derived MSM/Ms. Phenotyping by X-ray micro-CT revealed that MSM/Ms has structurally more fragile trabecular bone than C57BL/6J. Toward identification of genetic determinants for the difference in fragility of trabecular bone between the two mouse strains, we employed phenotype screening of consomic mouse strains in which each C57BL/6J chromosome is substituted by its counterpart from MSM/Ms. The results showed that many chromosomes affect trabecular bone structure, and that the consomic strain B6-Chr15^MSM, carrying MSM/Ms-derived chromosome 15 (Chr15), has the lowest values for the parameters BV/TV, Tb.N, and Conn.D, and the highest values for the parameters Tb.Sp and SMI. Subsequent phenotyping of subconsomic strains for Chr15 mapped four novel trabecular bone structure-related QTL (Tbsq1-4) on mouse Chr15. These results collectively indicate that genetic regulation of trabecular bone structure is highly complex, and that even in the single Chr15, the combined action of the four Tbsqs controls the fragility of trabecular bone. Given that Tbsq4 is syntenic to human Chr 12q12-13.3, where several bone-related SNPs are assigned, further study of Tbsq4 should facilitate our understanding of the genetic regulation of bone formation in humans.

CK2.1, a novel peptide, induces articular cartilage formation in vivo

Bone morphogenetic protein 2 regulates chondrogenesis and cartilage formation. However, it also induces chondrocyte hypertrophy and cartilage matrix degradation. We recently designed three peptides CK2.1, CK2.2, and CK2.3 that activate the BMP signaling pathways by releasing casein kinase II (CK2) from distinct sites at the bone morphogenetic protein receptor type Ia (BMPRIa). Since BMP2 is a major regulator of chondrogenesis and the peptides activated BMP signaling in a similar way, we evaluated the effect of these peptides on chondrogenesis and cartilage formation. C3H10T1/2 cells were stimulated with CK2.1, CK2.2, and CK2.3 and evaluated for the chondrogenic and osteogenic potential. For chondrogenesis, Alcian blue staining was performed. Additionally, collagen types II and X expression was measured. For osteogenesis, osteocalcin and von Kossa staining were performed. From the three peptides, CK2.1 was the most promising peptide to induce chondrogenesis but not osteogenesis. To investigate the effect of CK2.1 on articular cartilage formation in vivo, we injected CK2.1 into the tail vein of mice. Injection of CK2.1 into the tail vein of mice led to increased articular cartilage formation but not BMD. In sharp contrast, injection of BMP2 led to increased BMD and expression of collagen type X, a marker of chondrocyte hypertrophy. MMP13 expression was unchanged. Our study demonstrates that CK2.1 drives chondrogenesis and cartilage formation without induction of chondrocyte hypertrophy. Peptide CK2.1 may, therefore, be a valuable therapeutic for cartilage degenerative diseases. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:876-885, 2017.

Genetics of aging bone

With aging, the skeleton experiences a number of changes, which include reductions in mass and changes in matrix composition, leading to fragility and ultimately an increase of fracture risk. A number of aspects of bone physiology are controlled by genetic factors, including peak bone mass, bone shape, and composition; however, forward genetic studies in humans have largely concentrated on clinically available measures such as bone mineral density (BMD). Forward genetic studies in rodents have also heavily focused on BMD; however, investigations of direct measures of bone strength, size, and shape have also been conducted. Overwhelmingly, these studies of the genetics of bone strength have identified loci that modulate strength via influencing bone size, and may not impact the matrix material properties of bone. Many of the rodent forward genetic studies lacked sufficient mapping resolution for candidate gene identification; however, newer studies using genetic mapping populations such as Advanced Intercrosses and the Collaborative Cross appear to have overcome this issue and show promise for future studies. The majority of the genetic mapping studies conducted to date have focused on younger animals and thus an understanding of the genetic control of age-related bone loss represents a key gap in knowledge.

Novel Lesions of Bones and Joints Associated with Chikungunya Virus Infection in Two Mouse Models of Disease: New Insights into Disease Pathogenesis

Chikungunya virus is an arbovirus spread predominantly by Aedes aegypti and Ae. albopictus mosquitoes, and causes debilitating arthralgia and arthritis. While these are common manifestations during acute infection and it has been suggested they can recur in patients chronically, gaps in knowledge regarding the pathogenesis still exist. Two established mouse models were utilized (adult IRF 3/7 -/- -/- and wild-type C57BL/6J mice) to evaluate disease manifestations in bones and joints at various timepoints. Novel lesions in C57BL/6J mice consisted of periostitis (91%) and foci of cartilage of necrosis (50% of mice at 21 DPI). Additionally, at 21 DPI, 50% and 75% of mice exhibited periosteal bone proliferation affecting the metatarsal bones, apparent via histology and μCT, respectively. μCT analysis did not reveal any alterations in trabecular bone volume measurements in C57BL/6J mice. Novel lesions demonstrated in IRF 3/7 -/- -/- mice at 5 DPI included focal regions of cartilage necrosis (20%), periosteal necrosis (66%), and multifocal ischemic bone marrow necrosis (100%). Contralateral feet in 100% of mice of both strains had similar, though milder lesions. Additionally, comparison of control IRF 3/7 -/- -/- and wild-type C57BL/6J mice demonstrated differences in cortical bone. These experiments demonstrate novel manifestations of disease similar to those occurring in humans, adding insight into disease pathogenesis, and representing new potential targets for therapeutic interventions. Additionally, results demonstrate the utility of μCT in studies of bone and joint pathology and illustrate differences in bone dynamics between mouse strains.

High-Resolution Peripheral Quantitative Computed Tomography (HR-PQCT) and Dual Energy X-Ray Absorptiometry (DXA) Measurements of Proximal Tibia in Patients Undergoing Total Knee Arthroplasty

The aim of the study was to assess bone morphometric indices of the proximal tibia and compare it with bone mineral density (BMD) at hip and lumbar spine and compare with the pQCT. Fifty consecutive patients who underwent total knee arthroplasty (TKA) had a dual energy X-ray absorptiometry (DXA) scan of the upper femur and lumbar spine. Upper tibial cuts were harvested from the tibial condyles and a DXA of tibial cuts was done during TKA. Bone morphometry studies were carried out using HR-pQCT. The bone mineral density at the hip was 0.54 ± 0.08 g/cm2 and spine was 0.73 ± 0.1 with a T score at the hip −2.23 ± 0.44 and spine −2.61 ± 0.45. The bone mineral density of the tibial cut was 0.356 ± 0.03 g/cm2 (P < 0.001) and T score was −6.58 ± 2.87 (P < 0.001). The average bone volume (BV) was 115.27 ± 40.45 mm3, trabecular number (Tb.N) was 1.45 ± 0.32 (1.009–2.37) and trabecular thickness (Tb.Th) was 0.181 ± 0.03 mm (0.111–0.268). The mean bone mineral density measured was 206.24 ± 50.58 mg HA/ccm. This study shows that there is highly significant difference between BMD measured by DXA and pQCT, and that bone morphometric analysis indicates that there is marked decrease in the mechanical properties of the bone in the proximal tibia due to knee Osteoarthritis and Osteoporosis indicating poor bone architecture and quality.

Intramembranous bone regeneration differs among common inbred mouse strains following marrow ablation

Various intact and post-injury bone phenotypes are heritable traits. In this study, we sought to determine if intramembranous bone regeneration following marrow ablation differed among common inbred mouse strains and to identify how early the differences appear. We found a ∼four-fold difference in the regenerated bone volume 21 days after marrow ablation in females from four inbred mouse strains: FVB/N (15.7 ± 8.1%, mean and standard deviation), C3H/He (15.5 ± 4.2%), C57BL/6 (12.2 ± 5.2%), and BALB/c (4.0 ± 4.4%); with BALB/c different from FVB/N (p = 0.007) and C3H/He (p = 0.002). A second experiment showed that FVB/N compared to BALB/c mice had more regenerated bone 7 and 14 days after ablation (p < 0.001), while at 21 days FVB/N mice had a greater fraction of mineralizing surface (p = 0.008) without a difference in mineral apposition rate. Thus, differences among strains are evident early during intramembranous bone regeneration following marrow ablation and appear to be associated with differences in osteogenic cell recruitment, but not osteoblast activity. The amount of regenerating bone was not correlated with other heritable traits such as the intact bone phenotype or soft tissue wound healing, suggesting that there may be independent genetic pathways for these traits.

High-throughput phenotyping and genetic linkage of cortical bone microstructure in the mouse

Background

Understanding cellular structure and organization, which plays an important role in biological systems ranging from mechanosensation to neural organization, is a complicated multifactorial problem depending on genetics, environmental factors, and stochastic processes. Isolating these factors necessitates the measurement and sensitive quantification of many samples in a reliable, high-throughput, unbiased manner. In this manuscript we present a pipelined approach using a fully automated framework based on Synchrotron-based X-ray Tomographic Microscopy (SRXTM) for performing a full 3D characterization of millions of substructures.

Results

We demonstrate the framework on a genetic study on the femur bones of in-bred mice. We measured 1300 femurs from a F2 cross experiment in mice without the growth hormone (which can confound many of the smaller structural differences between strains) and characterized more than 50 million osteocyte lacunae (cell-sized hollows in the bone). The results were then correlated with genetic markers in a process called quantitative trait localization (QTL). Our findings provide a mapping between regions of the genome (all 19 autosomes) and observable phenotypes which could explain between 8–40 % of the variance using between 2–10 loci for each trait. This map shows 4 areas of overlap with previous studies looking at bone strength and 3 areas not previously associated with bone.

Conclusions

The mapping of microstructural phenotypes provides a starting point for both structure-function and genetic studies on murine bone structure and the specific loci can be investigated in more detail to identify single gene candidates which can then be translated to human investigations. The flexible infrastructure offers a full spectrum of shape, distribution, and connectivity metrics for cellular networks and can be adapted to a wide variety of materials ranging from plant roots to lung tissue in studies requiring high sample counts and sensitive metrics such as the drug-gene interactions and high-throughput screening.

Electronic supplementary material

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

Systemic injection of CK2.3, a novel peptide acting downstream of bone morphogenetic protein receptor BMPRIa, leads to increased trabecular bone mass

Bone Morphogenetic Protein 2 (BMP2) regulates bone integrity by driving both osteogenesis and osteoclastogenesis. However, BMP2 as a therapeutic has significant drawbacks. We have designed a novel peptide CK2.3 that blocks the interaction of Casein Kinase 2 (CK2) with Bone Morphogenetic Protein Receptor type Ia (BMPRIa), thereby activating BMP signaling pathways in the absence of ligand. Here, we show that CK2.3 induced mineralization in primary osteoblast cultures isolated from calvaria and bone marrow stromal cells (BMSCs) of 8 week old mice. Further, systemic tail vein injections of CK2.3 in 8 week old mice resulted in increased bone mineral density (BMD) and mineral apposition rate (MAR). In situ immunohistochemistry of the femur found that CK2.3 injection induced phosphorylation of extracellular signal-related kinase (ERK), but not Smad in osteocytes and osteoblasts, suggesting that CK2.3 signaling occurred through Smad independent pathway. Finally mice injected with CK2.3 exhibited decreased osteoclast differentiation and osteoclast activity. These data indicate that the novel mimetic peptide CK2.3 activated BMPRIa downstream signaling to enhance bone formation without the increase in osteoclast activity that accompanies BMP 2 stimulation.

Quantitative trait loci that modulate trabecular bone's risk of failure during unloading and reloading

Genetic makeup of an individual is a strong determinant of the morphologic and mechanical properties of bone. Here, in an effort to identify quantitative trait loci (QTLs) for changes in the simulated mechanical parameters of trabecular bone during altered mechanical demand, we subjected 352 second generation female adult (16 weeks old) BALBxC3H mice to 3 weeks of hindlimb unloading followed by 3 weeks of reambulation. Longitudinal in vivo microcomputed tomography (μCT) scans tracked trabecular changes in the distal femur. Tomographies were directly translated into finite element (FE) models and subjected to a uniaxial compression test. Apparent trabecular stiffness and components of the Von Mises (VM) stress distributions were computed for the distal metaphysis and associated with QTLs. At baseline, five QTLs explained 20% of the variation in trabecular peak stresses across the mouse population. During unloading, three QTLs accounted for 14% of the variability in peak stresses. During reambulation, one QTL accounted for 5% of the variability in peak stresses. QTLs were also identified for mechanically induced changes in stiffness, median stress values and skewness of stress distributions. There was little overlap between QTLs identified for baseline and QTLs for longitudinal changes in mechanical properties, suggesting that distinct genes may be responsible for the mechanical response of trabecular bone. Unloading related QTLs were also different from reambulation related QTLs. Further, QTLs identified here for mechanical properties differed from previously identified QTLs for trabecular morphology, perhaps revealing novel gene targets for reducing fracture risk in individuals exposed to unloading and for maximizing the recovery of trabecular bone's mechanical properties during reambulation.

Sclerostin deficient mice rapidly heal bone defects by activating β-catenin and increasing intramembranous ossification

We investigated the influence of the osteocyte protein, sclerostin, on fracture healing by examining the dynamics and mechanisms of repair of single-cortex, stabilized femoral defects in sclerostin knockout (Sost(-/-); KO) and sclerostin wild-type (Sost(+/+); WT) mice. Fourteen days following generation of bone defects, Sost KO mice had significantly more bone in the healing defect than WT mice. The increase in regenerating bone was due to an increase in the thickness of trabecularized spicules, osteoblast numbers and surfaces within the defect. Enhanced healing of bone defects in Sost KO mice was associated with significantly more activated β-catenin expression than observed in WT mice. The findings were similar to those observed in Axin2(-/-) mice, in which β-catenin signaling is known to be enhanced to facilitate bone regeneration. Taken together, these data indicate that enhanced β-catenin signaling is present in Sost(-/-) mice that demonstrate accelerated healing of bone defects, suggesting that modulation of β-catenin signaling in bone could be used to promote fracture repair.

Comparative anatomy of mouse and human nail units

Recent studies of mice with hair defects have resulted in major contributions to the understanding of hair disorders. To use mouse models as a tool to study nail diseases, a basic understanding of the similarities and differences between the human and mouse nail unit is required. In this study we compare the human and mouse nail unit at the macroscopic and microscopic level and use immunohistochemistry to determine the keratin expression patterns in the mouse nail unit. Both species have a proximal nail fold, cuticle, nail matrix, nail bed, nail plate, and hyponychium. Distinguishing features are the shape of the nail and the presence of an extended hyponychium in the mouse. Expression patterns of most keratins are similar. These findings indicate that the mouse nail unit shares major characteristics with the human nail unit and overall represents a very similar structure, useful for the investigation of nail diseases and nail biology.

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.

Genotype imputation accuracy in a F2 pig population using high density and low density SNP panels

Background

F₂ resource populations have been used extensively to map QTL segregating between pig breeds. A limitation associated with the use of these resource populations for fine mapping of QTL is the reduced number of founding individuals and recombinations of founding haplotypes occurring in the population. These limitations, however, become advantageous when attempting to impute unobserved genotypes using within family segregation information. A trade-off would be to re-type F₂ populations using high density SNP panels for founding individuals and low density panels (tagSNP) in F₂ individuals followed by imputation. Subsequently a combined meta-analysis of several populations would provide adequate power and resolution for QTL mapping, and could be achieved at relatively low cost. Such a strategy allows the wealth of phenotypic information that has previously been obtained on experimental resource populations to be further mined for QTL identification. In this study we used experimental and simulated high density genotypes (HD-60K) from an F₂ cross to estimate imputation accuracy under several genotyping scenarios.

Results

Selection of tagSNP using physical distance or linkage disequilibrium information produced similar imputation accuracies. In particular, tagSNP sets averaging 1 SNP every 2.1 Mb (1,200 SNP genome-wide) yielded imputation accuracies (IA) close to 0.97. If instead of using custom panels, the commercially available 9K chip is used in the F₂, IA reaches 0.99. In order to attain such high imputation accuracy the F₀ and F₁ generations should be genotyped at high density. Alternatively, when only the F₀ is genotyped at HD, while F₁ and F₂ are genotyped with a 9K panel, IA drops to 0.90.

Conclusions

Combining 60K and 9K panels with imputation in F₂ populations is an appealing strategy to re-genotype existing populations at a fraction of the cost.

Sclerostin alters serum vitamin D metabolite and fibroblast growth factor 23 concentrations and the urinary excretion of calcium

Inactivating mutations of the SOST (sclerostin) gene are associated with overgrowth and sclerosis of the skeleton. To determine mechanisms by which increased amounts of calcium and phosphorus are accreted to enable enhanced bone mineralization in the absence of sclerostin, we measured concentrations of calciotropic and phosphaturic hormones, and urine and serum calcium and inorganic phosphorus in mice in which the sclerostin (sost) gene was replaced by the β-D-galactosidase (lacZ) gene in the germ line. Knockout (KO) (sost(-/-)) mice had increased bone mineral density and content, increased cortical and trabecular bone thickness, and greater net bone formation as a result of increased osteoblast and decreased osteoclast surfaces compared with wild-type (WT) mice. β-Galactosidase activity was detected in osteocytes of sost KO mice but was undetectable in WT mice. Eight-week-old, male sost KO mice had increased serum 1α,25-dihydroxyvitamin D, decreased 24,25-dihydroxyvitamin D, decreased intact fibroblast growth factor 23, and elevated inorganic phosphorus concentrations compared with age-matched WT mice. 25-Hydroxyvitamin D 1α-hydroxylase cytochrome P450 (cyp27B1) mRNA was increased in kidneys of sost KO mice compared with WT mice. Treatment of cultured proximal tubule cells with mouse recombinant sclerostin decreased cyp27B1 mRNA transcripts. Urinary calcium and renal fractional excretion of calcium were decreased in sost KO mice compared with WT mice. Sost KO and WT mice had similar serum calcium and parathyroid hormone concentrations. The data show that sclerostin not only alters bone mineralization, but also influences mineral metabolism by altering concentrations of hormones that regulate mineral accretion.

Weak genetic relationship between trabecular bone morphology and obesity in mice

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

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

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

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.

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.

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.

Fluoride's effects on the formation of teeth and bones, and the influence of genetics

Fluorides are present in the environment. Excessive systemic exposure to fluorides can lead to disturbances of bone homeostasis (skeletal fluorosis) and enamel development (dental/enamel fluorosis). The severity of dental fluorosis is also dependent upon fluoride dose and the timing and duration of fluoride exposure. Fluoride's actions on bone cells predominate as anabolic effects both in vitro and in vivo. More recently, fluoride has been shown to induce osteoclastogenesis in mice. Fluorides appear to mediate their actions through the MAPK signaling pathway and can lead to changes in gene expression, cell stress, and cell death. Different strains of inbred mice demonstrate differential physiological responses to ingested fluoride. Genetic studies in mice are capable of identifying and characterizing fluoride-responsive genetic variations. Ultimately, this can lead to the identification of at-risk human populations who are susceptible to the unwanted or potentially adverse effects of fluoride action and to the elucidation of fundamental mechanisms by which fluoride affects biomineralization.

Microfibril-associated glycoprotein-1, an extracellular matrix regulator of bone remodeling

MAGP1 is an extracellular matrix protein that, in vertebrates, is a ubiquitous component of fibrillin-rich microfibrils. We previously reported that aged MAGP1-deficient mice (MAGP1Delta) develop lesions that are the consequence of spontaneous bone fracture. We now present a more defined bone phenotype found in MAGP1Delta mice. A longitudinal DEXA study demonstrated age-associated osteopenia in MAGP1Delta animals and muCT confirmed reduced bone mineral density in the trabecular and cortical bone. Further, MAGP1Delta mice have significantly less trabecular bone, the trabecular microarchitecture is more fragmented, and the diaphyseal cross-sectional area is significantly reduced. The remodeling defect seen in MAGP1Delta mice is likely not due to an osteoblast defect, because MAGP1Delta bone marrow stromal cells undergo osteoblastogenesis and form mineralized nodules. In vivo, MAGP1Delta mice exhibit normal osteoblast number, mineralized bone surface, and bone formation rate. Instead, our findings suggest increased bone resorption is responsible for the osteopenia. The number of osteoclasts derived from MAGP1Delta bone marrow macrophage cells is increased relative to the wild type, and osteoclast differentiation markers are expressed at earlier time points in MAGP1Delta cells. In vivo, MAGP1Delta mice have more osteoclasts lining the bone surface. RANKL (receptor activator of NF-kappaB ligand) expression is significantly higher in MAGP1Delta bone, and likely contributes to enhanced osteoclastogenesis. However, bone marrow macrophage cells from MAGP1Delta mice show a higher propensity than do wild-type cells to differentiate to osteoclasts in response to RANKL, suggesting that they are also primed to respond to osteoclast-promoting signals. Together, our findings suggest that MAGP1 is a regulator of bone remodeling, and its absence results in osteopenia associated with an increase in osteoclast number.

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

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

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

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

Candidate gene analysis of femoral neck trabecular and cortical volumetric bone mineral density in older men

In contrast to conventional dual-energy X-ray absorptiometry, quantitative computed tomography separately measures trabecular and cortical volumetric bone mineral density (vBMD). Little is known about the genetic variants associated with trabecular and cortical vBMD in humans, although both may be important for determining bone strength and osteoporotic risk. In the current analysis, we tested the hypothesis that there are genetic variants associated with trabecular and cortical vBMD at the femoral neck by genotyping 4608 tagging and potentially functional single-nucleotide polymorphisms (SNPs) in 383 bone metabolism candidate genes in 822 Caucasian men aged 65 years or older from the Osteoporotic Fractures in Men Study (MrOS). Promising SNP associations then were tested for replication in an additional 1155 men from the same study. We identified SNPs in five genes (IFNAR2, NFATC1, SMAD1, HOXA, and KLF10) that were robustly associated with cortical vBMD and SNPs in nine genes (APC, ATF2, BMP3, BMP7, FGF18, FLT1, TGFB3, THRB, and RUNX1) that were robustly associated with trabecular vBMD. There was no overlap between genes associated with cortical vBMD and trabecular vBMD. These findings identify novel genetic variants for cortical and trabecular vBMD and raise the possibility that some genetic loci may be unique for each bone compartment.

Nephroblastoma overexpressed (Nov) inactivation sensitizes osteoblasts to bone morphogenetic protein-2, but nov is dispensable for skeletal homeostasis

Overexpression of nephroblastoma overexpressed (Nov), a member of the Cyr 61, connective tissue growth factor, Nov family of proteins, inhibits osteoblastogenesis and causes osteopenia. The consequences of Nov inactivation on osteoblastogenesis and the postnatal skeleton are not known. To study the function of Nov, we inactivated Nov by homologous recombination. Nov null mice were maintained in a C57BL/6 genetic background after the removal of the neomycin selection cassette and compared with wild-type controls of identical genetic composition. Nov null mice were identified by genotyping and absent Nov mRNA in calvarial extracts and osteoblast cultures. Nov null mice did not exhibit developmental skeletal abnormalities or postnatal changes in weight, femoral length, body fat, or bone mineral density and appeared normal. Bone volume and trabecular number were decreased only in 1-month-old female mice. In older mice, after 7 months of age, osteoblast surface and bone formation were increased in females, and osteoclast and eroded surfaces were increased in male Nov null mice. Calvarial osteoblasts from Nov null mice displayed enhanced alkaline phosphatase activity, alkaline phosphatase mRNA, and transactivation of a bone morphogenetic protein (BMP)/phosphorylated mothers against decapentaplegic reporter construct in response to BMP-2. Similar results were obtained after the down-regulation of Nov by RNA interference in ST-2 stromal and MC3T3 cells. Osteoclast number was increased in marrow stromal cell cultures from Nov null mice. Surface plasmon resonance demonstrated direct interactions between Nov and BMP-2. In conclusion, Nov sensitizes osteoblasts to BMP-2, but Nov is dispensable for the maintenance of bone mass.

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.

Hierarchical microimaging of bone structure and function

With recent advances in molecular medicine and disease treatment in osteoporosis, quantitative image processing of three-dimensional bone structures is critical in the context of bone quality assessment. Biomedical imaging technology such as MRI or CT is readily available, but few attempts have been made to expand the capabilities of these systems by integrating quantitative analysis tools and by exploring structure-function relationships in a hierarchical fashion. Nevertheless, such quantitative end points are an important factor for success in basic research and in the development of novel therapeutic strategies. CT is key to these developments, as it images and quantifies bone in three dimensions and provides multiscale biological imaging capabilities with isotropic resolutions of a few millimeters (clinical CT), a few tens of micrometers (microCT) and even as high as 100 nanometers (nanoCT). The technology enables the assessment of the relationship between microstructural and ultrastructural measures of bone quality and certain diseases or therapies. This Review focuses on presenting strategies for three-dimensional approaches to hierarchical biomechanical imaging in the study of microstructural and ultrastructural bone failure. From this Review, it can be concluded that biomechanical imaging is extremely valuable for the study of bone failure mechanisms at different hierarchical levels.

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.

Molecular imaging of vessels in mouse models of disease

Vascular imaging of angiogenesis in mouse models of disease requires multi modal imaging hardware capable of targeting both structure and function at different physical scales. The three dimensional (3D) structure and function vascular information allows for accurate differentiation between biological processes. For example, image analysis of vessel development in angiogenesis vs. arteriogenesis enables more accurate detection of biological variation between subjects and more robust and reliable diagnosis of disease. In the recent years a number of micro imaging modalities have emerged in the field as preferred means for this purpose. They provide 3D volumetric data suitable for analysis, quantification, validation, and visualization of results in animal models. This review highlights the capabilities of microCT, ultrasound and microPET for multimodal imaging of angiogenesis and molecular vascular targets in a mouse model of tumor angiogenesis. The basic principles of the imaging modalities are described and experimental results are presented.

A three-dimensional microcomputed tomographic study of site-specific variation in trabecular microarchitecture in the human second metacarpal

Variation in trabecular microarchitecture is widely accepted as being regulated by both functional (mechanical loading) and genetic parameters, although the relative influence of each is unclear. Studies reporting inter-site differences in trabecular morphology (volume, number and structure) reveal a complex interaction at the gene-environment interface. We report inter- and intra-site variation in trabecular anatomy using a novel model of contralateral (left vs right) and ipsilateral (head vs base) comparisons for the human second metacarpal in a sample of n = 29 historically known 19th century EuroCanadians. Measures of bone volume fraction, structure model index, connectivity, trabecular number, spacing and thickness as well as degree of anisotropy were obtained from 5-mm volumes of interest using three-dimensional microcomputed tomography. We hypothesized that: (i) the more diverse loading environment of metacarpal heads should produce a more robust trabecular architecture than corresponding bases within sides and (ii) the ipsilateral differences between epiphyses will be larger on the right side than on the left side, as a function of handedness. Analysis of covariance (Side x Epiphysis) with Age as covariate revealed a clear dichotomy between labile and constrained architectures within and among anatomical sites. The predicted variation in loading was accommodated by changes in trabecular volume, whereas trabecular structure did not vary significantly by side or by epiphysis within sides. Age was a significant covariate only for females. We conclude that environmental and genetic regulation of bone adaptation may act through distinct pathways and local anatomies to ensure an integrated lattice of sufficient mass to meet normal functional demands.

Phenotypic integration among trabecular and cortical bone traits establishes mechanical functionality of inbred mouse vertebrae

Conventional approaches to identifying quantitative trait loci (QTLs) regulating bone mass and fragility are limited because they examine cortical and trabecular traits independently. Prior work examining long bones from young adult mice and humans indicated that skeletal traits are functionally related and that compensatory interactions among morphological and compositional traits are critical for establishing mechanical function. However, it is not known whether trait covariation (i.e., phenotypic integration) also is important for establishing mechanical function in more complex, corticocancellous structures. Covariation among trabecular, cortical, and compositional bone traits was examined in the context of mechanical functionality for L(4) vertebral bodies across a panel of 16-wk-old female AXB/BXA recombinant inbred (RI) mouse strains. The unique pattern of randomization of the A/J and C57BL/6J (B6) genome among the RI panel provides a powerful tool that can be used to measure the tendency for different traits to covary and to study the biology of complex traits. We tested the hypothesis that genetic variants affecting vertebral size and mass are buffered by changes in the relative amounts of cortical and trabecular bone and overall mineralization. Despite inheriting random sets of A/J and B6 genomes, the RI strains inherited nonrandom sets of cortical and trabecular bone traits. Path analysis, which is a multivariate analysis that shows how multiple traits covary simultaneously when confounding variables like body size are taken into consideration, showed that RI strains that tended to have smaller vertebrae relative to body size achieved mechanical functionality by increasing mineralization and the relative amounts of cortical and trabecular bone. The interdependence among corticocancellous traits in the vertebral body indicated that variation in trabecular bone traits among inbred mouse strains, which is often thought to arise from genetic factors, is also determined in part by the adaptive response to variation in traits describing the cortical shell. The covariation among corticocancellous traits has important implications for genetic analyses and for interpreting the response of bone to genetic and environmental perturbations.

Bone health and back pain: what do we know and where should we go?

Bone health is generally not considered in patients who present with chronic back pain. Nonetheless, bone health and back pain share common genetic and environmental correlates suggesting a co-dependence. Evidence exists for a relationship between back pain and impaired bone health. Here we present the evidence, theoretic framework and clinical relevance. Bone health and back pain are important determinants of musculoskeletal health. Back pain experienced in youth is a risk factor for future back pain, while suboptimal bone health during development increases the risk of skeletal fragility in later life. Generally, bone health is not considered in patients with chronic back pain who do not demonstrate other well-recognised bone health risk factors or associated conditions. Nonetheless, evidence suggests that back pain and impaired bone health share common environmental and genetic correlates, indicating that bone health ought to be considered in the context of back pain in otherwise healthy individuals. This review describes the likely mechanisms explaining the relationship between back pain and impaired bone health, evidence concerning the relationship and suggestions for future research. A narrative literature search was conducted using CINAHL, Medline, PubMed and Web of Science electronic databases. A history of back pain is associated with decreased bone mineral density in adults, yet this tends to be site-specific. No studies were identified examining this association in youth, yet the negative effects of childhood skeletal trauma and obesity on bone and spinal health provide indirect evidence for an association. Further research is required to clarify the impact of back pain on bone health at different lifespan stages using prospective cohort designs.

Fourier transform infrared imaging microspectroscopy and tissue-level mechanical testing reveal intraspecies variation in mouse bone mineral and matrix composition

Fracture susceptibility is heritable and dependent upon bone morphology and quality. However, studies of bone quality are typically overshadowed by emphasis on bone geometry and bone mineral density. Given that differences in mineral and matrix composition exist in a variety of species, we hypothesized that genetic variation in bone quality and tissue-level mechanical properties would also exist within species. Sixteen-week-old female A/J, C57BL/6J (B6), and C3H/HeJ (C3H) inbred mouse femora were analyzed using Fourier transform infrared imaging and tissue-level mechanical testing for variation in mineral composition, mineral maturity, collagen cross-link ratio, and tissue-level mechanical properties. A/J femora had an increased mineral-to-matrix ratio compared to B6. The C3H mineral-to-matrix ratio was intermediate of A/J and B6. C3H femora had reduced acid phosphate and carbonate levels and an increased collagen cross-link ratio compared to A/J and B6. Modulus values paralleled mineral-to-matrix values, with A/J femora being the most stiff, B6 being the least stiff, and C3H having intermediate stiffness. In addition, work-to-failure varied among the strains, with the highly mineralized and brittle A/J femora performing the least amount of work-to-failure. Inbred mice are therefore able to differentially modulate the composition of their bone mineral and the maturity of their bone matrix in conjunction with tissue-level mechanical properties. These results suggest that specific combinations of bone quality and morphological traits are genetically regulated such that mechanically functional bones can be constructed in different ways.

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.

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.

Trabecular bone deterioration in col9a1+/- mice associated with enlarged osteoclasts adhered to collagen IX-deficient bone

INTRODUCTION

Short collagen IX, the exclusive isoform expressed by osteoblasts, is synthesized through alternative transcription of the col9a1 gene. The function of short collagen IX in bone was characterized in col9a1-null mutant mice.

MATERIALS AND METHODS

Trabecular bone morphometry of lumbar bones and tibias was evaluated by muCT and nondecalcified histology. Osteoblastic and osteoclastic activities were evaluated by PCR- and microarray-based gene expression assays and TRACP-5b and C-terminal telopeptide (CTX) assays, as well as in vitro using bone marrow stromal cells and splenocytes. The effect of col9a1(+/-) mutation on osteoclast morphology was evaluated using RAW264.7-derived osteoclastic cells cultured on the mutant or wildtype calvarial bone substrates.

RESULTS

Col9a1 knockout mutation caused little effects on the skeletal development; however, young adult female col9a1(-/-) and col9a1(+/-) mice exhibited significant loss of trabecular bone. The trabecular bone architecture was progressively deteriorated in both male and female heterozygous col9a1(+/-) mice while aging. The aged mutant mice also exhibited signs of thoracic kyphosis and weight loss, resembling the clinical signs of osteoporosis. The col9a1(+/-) osteoblasts synthesized short col9a1 transcripts at decreased rates. Whereas bone formation activities in vitro and in vivo were not affected, the mutant osteoblast expressed the elevated ratio of RANKL/osteoprotegerin. Increased serum TRACP-5b and CTX levels were found in col9a1(+/-) mice, whose bone surface was associated with osteoclastic cells that were abnormally flattened and enlarged. The mutant and wildtype splenocytes underwent similar osteoclastogenesis in vitro; however, RAW264.7-derived osteoclastic cells, when cultured on the col9a1(+/-) calvaria, widely spread over the bone surface and formed large resorption pits. The surface of col9a1(+/-) calvaria was found to lack the typical nanotopography.

CONCLUSIONS

The mineralized bone matrix deficient of short collagen IX may become susceptible to osteoclastic bone resorption, possibly through a novel non-cell-autonomous mechanism. The data suggest the involvement of bone collagen IX in the pathogenesis of osteoporosis.

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.

Resolution dependence of the non-metric trabecular structure indices

Non-metric indices of topological features of trabecular bone structure, such as structure model index (SMI), connectivity density (Conn.D), and degree of anisotropy (DA), provide unique information relevant to bone quality. With recent technological advancement, in vivo assessment of these indices may be possible from images acquired using high-resolution imaging techniques such as high-resolution peripheral quantitative computed tomography (HR-pQCT). However, more detailed investigation of the dependence of non-metric indices on spatial resolution is needed to determine their applicability. The purpose of this study was to determine whether these three non-metric indices are affected by the spatial resolution of CT images. First, the SMI, Conn.D, and DA were calculated for trabecular bone specimens with varying plate-like and rod-like structures from resampled muCT images across a range of spatial resolutions and compared to the reference values. To account for differences in size across different species and anatomical sites, the results are reported in normalized resolution units. Next, the impact of resolution on the non-metric indices for cores of human distal tibia trabecular bone from clinical HR-pQCT images was evaluated to determine the applicability of the non-metric indices to in vivo imaging. We found that the non-metric indices of trabecular bone structure were affected by spatial resolution of CT images. Particularly, the SMI deviated from the high-resolution muCT reference value depending on the structure type, whether plate-like or rod-like. Both Conn.D and DA were underestimated in the images obtained at an in vivo resolution. It is not trivial to determine absolute threshold for validity of these non-metric indices without considering a specific study design (e.g. relative resolution, the size of the treatment effect to detect, and specimen type). The results of this study provide an upper bound for the accuracy of the non-metric indices under limited resolution scenarios.