Oxidative metabolism is impaired by phosphate deficiency during fracture healing and is mechanistically related to BMP induced chondrocyte differentiation.

Prior studies of acute phosphate restriction during the endochondral phase of fracture healing showed delayed chondrocyte differentiation was mechanistically linked to decreased bone morphogenetic protein signaling. In the present study, transcriptomic analysis of fracture callus gene expression in three strains of mice was used to identify differentially expressed (FDR = q ≤ 0.05) genes in response to phosphate (Pi) restriction. Ontology and pathway analysis of these genes showed that independent of genetic background, a Pi-deficient diet downregulated (p = 3.16 × 10-23) genes associated with mitochondrial oxidative phosphorylation pathways as well as multiple other pathways of intermediate metabolism. Temporal clustering was used to identify co-regulation of these specific pathways. This analysis showed that specific Ox/Phos, tricarboxylic acid cycle, pyruvate dehydrogenase. Arginine, proline metabolism genes, and prolyl 4-hydroxylase were all co-regulated in response to dietary Pi restriction. The murine C3H10T½ mesenchymal stem cell line was used to assess the functional relationships between BMP2-induced chondrogenic differentiation, oxidative metabolism and extracellular matrix formation. BMP2-induced chondrogenic differentiation of C3H10T½ was carried out in culture media in the absence or presence of ascorbic acid, the necessary co-factor for proly hydroxylation, and in media with normal and 25 % phosphate levels. BMP2 treatment led to decreased proliferation, increased protein accumulation and increased collagen and aggrecan gene expression. Across all conditions, BMP2 increased total oxidative activity and ATP synthesis. Under all conditions, the presence of ascorbate further increased total protein accumulation, proly-hydroxylation and aggrecan gene expression, oxidative capacity and ATP production. Lower phosphate levels only diminished aggrecan gene expression with no other effects of metabolic activity being observed. These data suggest that dietary phosphate restriction controls endochondral growth in vivo indirectly through BMP signaling, which upregulates oxidative activity that is linked to overall protein production and collagen hydroxylation.

Dosimetry of [212Pb]VMT01, a MC1R-Targeted Alpha Therapeutic Compound, and Effect of Free 208Tl on Tissue Absorbed Doses.

[212Pb]VMT01 is a melanocortin 1 receptor (MC1R) targeted theranostic ligand in clinical development for alpha particle therapy for melanoma. 212Pb has an elementally matched gamma-emitting isotope 203Pb; thus, [203Pb]VMT01 can be used as an imaging surrogate for [212Pb]VMT01. [212Pb]VMT01 human serum stability studies have demonstrated retention of the 212Bi daughter within the chelator following beta emission of parent 212Pb. However, the subsequent alpha emission from the decay of 212Bi into 208Tl results in the generation of free 208Tl. Due to the 10.64-hour half-life of 212Pb, accumulation of free 208Tl in the injectate will occur. The goal of this work is to estimate the human dosimetry for [212Pb]VMT01 and the impact of free 208Tl in the injectate on human tissue absorbed doses. Human [212Pb]VMT01 tissue absorbed doses were estimated from murine [203Pb]VMT01 biodistribution data, and human biodistribution values for 201Tl chloride (a cardiac imaging agent) from published data were used to estimate the dosimetry of free 208Tl. Results indicate that the dose-limiting tissues for [212Pb]VMT01 are the red marrow and the kidneys, with estimated absorbed doses of 1.06 and 8.27 mGyRBE = 5/MBq. The estimated percent increase in absorbed doses from free 208Tl in the injectate is 0.03% and 0.09% to the red marrow and the kidneys, respectively. Absorbed doses from free 208Tl result in a percent increase of no more than 1.2% over [212Pb]VMT01 in any organ or tissue. This latter finding indicates that free 208Tl in the [212Pb]VMT01 injectate will not substantially impact estimated tissue absorbed doses in humans.

Temporal and Quantitative Transcriptomic Differences Define Sexual Dimorphism in Murine Postnatal Bone Aging.

Time is a central element of the sexual dimorphic patterns of development, pathology, and aging of the skeleton. Because the transcriptome is a representation of the phenome, we hypothesized that both sex and sex-specific temporal, transcriptomic differences in bone tissues over an 18-month period would be informative to the underlying molecular processes that lead to postnatal sexual dimorphism. Regardless of age, sex-associated changes of the whole bone transcriptomes were primarily associated not only with bone but also vascular and connective tissue ontologies. A pattern-based approach used to screen the entire Gene Expression Omnibus (GEO) database against those that were sex-specific in bone identified two coordinately regulated gene sets: one related to high phosphate-induced aortic calcification and one induced by mechanical stimulation in bone. Temporal clustering of the transcriptome identified two skeletal tissue-associated, sex-specific patterns of gene expression. One set of genes, associated with skeletal patterning and morphology, showed peak expression earlier in females. The second set of genes, associated with coupled remodeling, had quantitatively higher expression in females and exhibited a broad peak between 3 to 12 months, concurrent with the animals' reproductive period. Results of phenome-level structural assessments of the tibia and vertebrae, and in vivo and in vitro analysis of cells having osteogenic potential, were consistent with the existence of functionally unique, skeletogenic cell populations that are separately responsible for appositional growth and intramedullary functions. These data suggest that skeletal sexual dimorphism arises through sex-specific, temporally different processes controlling morphometric growth and later coupled remodeling of the skeleton during the reproductive period of the animal. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Tributyltin protects against ovariectomy-induced trabecular bone loss in C57BL/6J mice with an attenuated effect in high fat fed mice.

Risk factors for poor bone quality include estrogen loss at menopause, a high fat diet and exposures to drugs/chemicals that activate peroxisome proliferator activated receptor gamma (PPARγ). We previously reported that the PPARγ and retinoid X receptor dual ligand, tributyltin (TBT), repressed periosteal bone formation but enhanced trabecular bone formation in vivo. Here, we examined the interaction of diet, ovariectomy (OVX) and TBT exposure on bone structure. C57BL/6J mice underwent either sham surgery or OVX at 10 weeks of age. At 12 weeks of age, they were placed on a low (10% kcal) or high (45% kcal) fat, sucrose-matched diet and treated with vehicle or TBT (1 or 5 mg/kg) for 14 weeks. OVX increased body weight gain in mice on either diet. TBT enhanced body weight gain in intact mice fed a high fat diet, but decreased weight gain in OVX mice. Elemental tin concentrations increased dose-dependently in bone. TBT had marginal effects on cortical and trabecular bone in intact mice fed either diet. OVX caused a reduction in cortical and trabecular bone, regardless of diet. In high fat fed OVX mice, TBT further reduced cortical thickness, bone area and total area. Interestingly, TBT protected against OVX-induced trabecular bone loss in low fat fed mice. The protective effect of TBT was nullified by the high fat. These results show that TBT protects against trabecular bone loss, even in the presence of a strongly resorptive environment, at an even lower level of exposure than we showed repressed homeostatic resorption.

Bayesian modeling of factorial time-course data with applications to a bone aging gene expression study.

Many scientific studies, especially in the biomedical sciences, generate data measured simultaneously over a multitude of units, over a period of time, and under different conditions or combinations of factors. Often, an important question of interest asked relates to which units behave similarly under different conditions, but measuring the variation over time complicates the analysis significantly. In this article we address such a problem arising from a gene expression study relating to bone aging, and develop a Bayesian statistical method that can simultaneously detect and uncover signals on three levels within such data: factorial, longitudinal, and transcriptional. Our model framework considers both cluster and time-point-specific parameters and these parameters uniquely determine the shapes of the temporal gene expression profiles, allowing the discovery and characterization of latent gene clusters based on similar underlying biological mechanisms. Our methodology was successfully applied to discover transcriptional networks in a microarray data set comparing the transcriptomic changes that occurred during bone aging in male and female mice expressing one or both copies of the bromodomain (Brd2) gene, a transcriptional regulator which exhibits an age-dependent sex-linked bone loss phenotype.

Are Narrative Letters of Recommendation for Medical Students Interpreted as Intended by Orthopaedic Surgery Residency Programs?

BACKGROUND: Narrative letters of recommendation are an important component of the residency application process. However, because narrative letters of recommendation are almost always positive, it is unclear whether those reviewing the letters understand the writer's intended strength of support for a given applicant. QUESTIONS/PURPOSES: (1) Is the perception of letter readers for narrative letters of recommendation consistent with the intention of the letter's author? (2) Is there inter-reviewer consistency in selection committee members' perceptions of the narrative letters of recommendation? METHODS: Letter writers who wrote two or more narrative letters of recommendation for applicants to one university-based orthopaedic residency program for the 2014 to 2015 application cycle were sent a survey linked to a specific letter of recommendation they authored to assess the intended meaning regarding the strength of an applicant. A total of 247 unstructured letters of recommendation and accompanying surveys were sent to their authors, and 157 surveys were returned and form the basis of this study (response percentage 64%). The seven core members of the admissions committee (of 22 total reviewers) at a university-based residency program were sent a similar survey regarding their perception of the letter. To answer our research question about whether letter readers' perceptions about a candidate were consistent with the letter writer's intention, we used kappa values to determine agreement for survey questions involving discrete variables and Spearman correlation coefficients (SCCs) to determine agreement for survey questions involving continuous variables. To answer our research question regarding inter-reviewer consistency among the seven faculty members, we compared the letter readers' responses to each survey question using intraclass correlation coefficients (ICCs). RESULTS: There was a negligible to moderate correlation between the intended and perceived strength of the letters (SCC 0.26 to 0.57), with only one of seven letter readers scoring in the moderate correlation category. When stratifying the applicants into thirds, there was only slight agreement (kappa 0.07 to 0.19) between the writers and reviewers. There were similarly low kappa values for agreement about how the writers and readers felt regarding the candidate matching into their program (kappa 0.14 to 0.30). The ICC for each question among the seven faculty reviewers ranged from poor to moderate (ICC 0.42 to 0.52). CONCLUSION: Our results demonstrate that the reader's perception of narrative letters of recommendation did not correlate well with the letter writer's intended meaning and was not consistent between letter readers at a single university-based urban orthopaedic surgery residency program. CLINICAL RELEVANCE: Given the low correlation between the intended strength of the letter writers and the perceived strength of those letters, we believe that other options such as a slider bar or agreed-upon wording as is used in many dean's letters may be helpful.

Association of vertebral endplate microstructure with bone strength in men and women.

Epidemiological and biomechanical evidence indicates that the risk of vertebral fracture differs between men and women, and that vertebral fracture frequently involves failure of the endplate region. The goal of this study was to compare the bone microstructure of the endplate region-defined as the (bony) vertebral endplate and underlying subchondral trabecular bone-between sexes and to determine whether any such sex differences are associated with vertebral strength. The bone density (volume fraction, apparent density and tissue mineral density) of the superior-most 2 mm of the vertebra, and the bone density and trabecular architecture of the next 5 mm were quantified using micro-computed tomography in human T8 (12 female, 16 male) and L1 (13 female, 12 male) vertebrae. Average density of the vertebra (integral bone mineral density (BMD)) was determined by quantitative computed tomography and compressive strength by mechanical testing. Few differences were found between male and female vertebrae in the density of the endplate region; none were found in trabecular architecture. However, whereas endplate volume fraction was positively correlated with integral BMD in male vertebrae (r = 0.654, p < .001), no correlation was found in the female vertebrae (r = 0.157, p = .455). Accounting for the density of the endplate region improved predictions of vertebral strength (p < .034) and eliminated sex-specificity in the strength prediction that was based on integral BMD alone. These results suggest that the density of the endplate region influences vertebral fracture and that non-invasive assessment of this region's density can contribute to predictions of vertebral strength in men and women.

Hypophosphatemia Regulates Molecular Mechanisms of Circadian Rhythm.

Transcriptomic analysis showed that the central circadian pathway genes had significantly altered expression in fracture calluses from mice fed a low phosphate diet. This led us to hypothesize that phosphate deficiency altered the circadian cycle in peripheral tissues. Analysis of the expression of the central clock genes over a 24-36 hour period in multiple peripheral tissues including fracture callus, proximal tibia growth plate and cardiac tissues after 12 days on a low phosphate diet showed higher levels of gene expression in the hypophosphatemia groups (p < 0.001) and a 3 to 6 hour elongation of the circadian cycle. A comparative analysis of the callus tissue transcriptome genes that were differentially regulated by hypophosphatemia with published data for the genes in bone that are diurnally regulated identified 1879 genes with overlapping differential regulation, which were shown by ontology assessment to be associated with oxidative metabolism and apoptosis. Network analysis of the central circadian pathway genes linked their expression to the up regulated expression of the histone methyltransferase gene EZH2, a gene that when mutated in both humans and mice controls overall skeletal growth. These data suggest that phosphate is an essential metabolite that controls circadian function in both skeletal and non skeletal peripheral tissues and associates its levels with the overall oxidative metabolism and skeletal growth of animals.

Bone Mechanical Properties in Healthy and Diseased States.

The mechanical properties of bone are fundamental to the ability of our skeletons to support movement and to provide protection to our vital organs. As such, deterioration in mechanical behavior with aging and/or diseases such as osteoporosis and diabetes can have profound consequences for individuals' quality of life. This article reviews current knowledge of the basic mechanical behavior of bone at length scales ranging from hundreds of nanometers to tens of centimeters. We present the basic tenets of bone mechanics and connect them to some of the arcs of research that have brought the field to recent advances. We also discuss cortical bone, trabecular bone, and whole bones, as well as multiple aspects of material behavior, including elasticity, yield, fracture, fatigue, and damage. We describe the roles of bone quantity (e.g., density, porosity) and bone quality (e.g., cross-linking, protein composition), along with several avenues of future research.

Carbonic anhydrase III protects osteocytes from oxidative stress.

Osteocytes are master orchestrators of bone remodeling; they control osteoblast and osteoclast activities both directly via cell-to-cell communication and indirectly via secreted factors, and they are the main postnatal source of sclerostin and RANKL (receptor activator of NF-kB ligand), two regulators of osteoblast and osteoclast function. Despite progress in understanding osteocyte biology and function, much remains to be elucidated. Recently developed osteocytic cell lines-together with new genome editing tools-has allowed a closer look at the biology and molecular makeup of these cells. By using single-cell cloning, we identified genes that are associated with high Sost/sclerostin expression and analyzed their regulation and function. Unbiased transcriptome analysis of high- vs. low-Sost/sclerostin-expressing cells identified known and novel genes. Dmp1 (dentin matrix protein 1), Dkk1 (Dickkopf WNT signaling pathway inhibitor 1), and Phex were among the most up-regulated known genes, whereas Srpx2, Cd200, and carbonic anhydrase III (CAIII) were identified as novel markers of differentiated osteocytes. Aspn, Enpp2, Robo2, Nov, and Serpina3g were among the transcripts that were most significantly suppressed in high-Sost cells. Considering that CAII was recently identified as being regulated by Sost/sclerostin and capable of controlling mineral homeostasis, we focused our attention on CAIII. Here, we report that CAIII is highly expressed in osteocytes, is regulated by parathyroid hormone both in vitro and in vivo, and protects osteocytes from oxidative stress.-Shi, C., Uda, Y., Dedic, C., Azab, E., Sun, N., Hussein, A. I., Petty, C. A., Fulzele, K., Mitterberger-Vogt, M. C., Zwerschke, W., Pereira, R., Wang, K., Divieti Pajevic, P. Carbonic anhydrase III protects osteocytes from oxidative stress.

Correlation between RUST assessments of fracture healing to structural and biomechanical properties.

Radiographic Union Score for Tibia (RUST) and modified RUST (mRUST) are radiographic tools for quantitatively evaluating fracture healing using a cortical scoring system. This tool has high intra-class correlation coefficients (ICCs); however, little evidence has evaluated the scores against the physical properties of bone healing. Closed, stabilized fractures were made in the femora of C3H/HeJ male mice (8-12 week-old) of two dietary groups: A control and a phosphate restricted diet group. Micro-computed tomography (µCT) and torsion testing were carried out at post-operative days (POD) 14, 21, 35, and 42 (n = 10-16) per group time-point. Anteroposterior and lateral radiographic views were constructed from the µCT scans and scored by five raters. The raters also indicated if the fracture were healed. ICCs were 0.71 (mRUST) and 0.63 (RUST). Both RUST scores were positively correlated with callus bone mineral density (BMD) (r = 0.85 and 0.80, p < 0.001) and bone volume fraction (BV/TV) (r = 0.86 and 0.80, p < 0.001). Both RUST scores positively correlated with callus strength (r = 0.35 and 0.26, p < 0.012) and rigidity (r = 0.50 and 0.39, p < 0.001). Radiographically healed calluses had a mRUST ≥13 and a RUST ≥10 and had excellent relationship to structural and biomechanical metrics. Effect of delayed healing due to phosphate dietary restrictions was found at later time points with all mechanical properties (p < 0.011), however no differences found in the RUST scores (p > 0.318). Clinical relevance of this study is both RUST scores showed high correlation to physical properties of healing and generally distinguished healed vs. non-healed fractures. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:945-953, 2018.

Differences in Trabecular Microarchitecture and Simplified Boundary Conditions Limit the Accuracy of Quantitative Computed Tomography-Based Finite Element Models of Vertebral Failure.

Vertebral fractures are common in the elderly, but efforts to reduce their incidence have been hampered by incomplete understanding of the failure processes that are involved. This study's goal was to elucidate failure processes in the lumbar vertebra and to assess the accuracy of quantitative computed tomography (QCT)-based finite element (FE) simulations of these processes. Following QCT scanning, spine segments (n = 27) consisting of L1 with adjacent intervertebral disks and neighboring endplates of T12 and L2 were compressed axially in a stepwise manner. A microcomputed tomography scan was performed at each loading step. The resulting time-lapse series of images was analyzed using digital volume correlation (DVC) to quantify deformations throughout the vertebral body. While some diversity among vertebrae was observed on how these deformations progressed, common features were large strains that developed progressively in the superior third and, concomitantly, in the midtransverse plane, in a manner that was associated with spatial variations in microstructural parameters such as connectivity density. Results of FE simulations corresponded qualitatively to the measured failure patterns when boundary conditions were derived from DVC displacements at the endplate. However, quantitative correspondence was often poor, particularly when boundary conditions were simplified to uniform compressive loading. These findings suggest that variations in trabecular microstructure are one cause of the differences in failure patterns among vertebrae and that both the lack of incorporation of these variations into QCT-based FE models and the oversimplification of boundary conditions limit the accuracy of these models in simulating vertebral failure.

Serum proteomic assessment of the progression of fracture healing.

A targeted proteomic analysis of murine serum over a 35-day course of fracture healing was carried out to determine if serum proteomic changes could be used to monitor the biological progression of fracture healing. Transverse, closed femoral fractures where generated and stabilized with intramedullary fixation. A single stranded DNA aptamer-based multiplexed proteomic approach was used to assay 1,310 proteins. The transcriptomic profiles for genes matching the 1,310 proteins were obtained by microarray analysis of callus mRNA. Of the 1,310 proteins analyzed, 850 proteins showed significant differences among the time points (p-value <0.05). Ontology assessment associated these proteins with osteoblasts, monocyte/macrophage lineages, mesenchymal stem cell lines, hepatic tissues, and lymphocytes. Temporal clustering of these data identified proteins associated with inflammation, cartilage formation and bone remodeling stages of healing. VEGF, Wnt, and TGF-ßsignaling pathways were restricted to the period of cartilage formation. Comparison of the proteomic and transcriptomic profiles showed that 87.5% of proteins in serum had concordant expression to their mRNA expression in the callus, while 12.5% of the protein and mRNA expression patterns were discordant. The discordant proteins that were elevated in the serum but down regulated in callus mRNA expression were related to clotting functions, allograft rejection, and complement function. While proteins down regulated in the serum and elevated in callus mRNA were associated with osteoblast function, NF-ĸb, and activin signaling. These data show the serum proteome may be used to monitor the different biological stages of fracture healing and have translational potential in assessing human fracture healing. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:1153-1163, 2018.

Acute Phosphate Restriction Impairs Bone Formation and Increases Marrow Adipose Tissue in Growing Mice.

Phosphate plays a critical role in chondrocyte maturation and skeletal mineralization. Studies examining the consequences of dietary phosphate restriction in growing mice demonstrated not only the development of rickets, but also a dramatic decrease in bone accompanied by increased marrow adipose tissue (MAT). Thus studies were undertaken to determine the effects of dietary phosphate restriction on bone formation and bone marrow stromal cell (BMSC) differentiation. Acute phosphate restriction of 28-day-old mice profoundly inhibited bone formation within 48 hours. It also resulted in increased mRNA expression of the early osteolineage markers Sox9 and Runt-related transcription factor 2 (Runx2), accompanied by decreased expression of the late osteolineage markers Osterix and Osteocalcin in BMSCs and osteoblasts, suggesting that phosphate restriction arrests osteoblast differentiation between Runx2 and Osterix. Increased expression of PPARγ and CEBPα, key regulators of adipogenic differentiation, was observed within 1 week of dietary phosphate restriction and was followed by a 13-fold increase in MAT at 3 weeks of phosphate restriction. In vitro phosphate restriction did not alter BMSC osteogenic or adipogenic colony formation, implicating aberrant paracrine or endocrine signaling in the in vivo phenotype. Because BMP signaling regulates the transition between Runx2 and Osterix, this pathway was interrogated. A dramatic decrease in pSmad1/5/9 immunoreactivity was observed in the osteoblasts of phosphate-restricted mice on day 31 (d31) and d35. This was accompanied by attenuated expression of the BMP target genes Id1, KLF10, and Foxc2, the latter of which promotes osteogenic and angiogenic differentiation while impairing adipogenesis. A decrease in expression of the Notch target gene Hey1, a BMP-regulated gene that governs angiogenesis, was also observed in phosphate-restricted mice, in association with decreased metaphyseal marrow vasculature. Whereas circulating phosphate levels are known to control growth plate maturation and skeletal mineralization, these studies reveal novel consequences of phosphate restriction in the regulation of bone formation and osteoblast differentiation. © 2016 American Society for Bone and Mineral Research.

Quantitative, 3D Visualization of the Initiation and Progression of Vertebral Fractures Under Compression and Anterior Flexion.

The biomechanical mechanisms leading to vertebral fractures are not well understood. Clinical and laboratory evidence suggests that the vertebral endplate plays a key role in failure of the vertebra as a whole, but how this role differs for different types of vertebral loading is not known. Mechanical testing of human thoracic spine segments, in conjunction with time-lapsed micro-computed tomography, enabled quantitative assessment of deformations occurring throughout the entire vertebral body under axial compression combined with anterior flexion ("combined loading") and under axial compression only ("compression loading"). The resulting deformation maps indicated that endplate deflection was a principal feature of vertebral failure for both loading modes. Specifically, the onset of endplate deflection was temporally coincident with a pronounced drop in the vertebra's ability to support loads. The location of endplate deflection, and also vertebral strength, were associated with the porosity of the endplate and the microstructure of the underlying trabecular bone. However, the location of endplate deflection and the involvement of the cortex differed between the two types of loading. Under the combined loading, deflection initiated, and remained the largest, at the anterior central endplate or the anterior ring apophysis, depending in part on health of the adjacent intervertebral disc. This deflection was accompanied by outward bulging of the anterior cortex. In contrast, the location of endplate deflection was more varied in compression loading. For both loading types, the earliest progression to a mild fracture according to a quantitative morphometric criterion occurred only after much of the failure process had occurred. The outcomes of this work indicate that for two physiological loading modes, the vertebral endplate and underlying trabecular bone are critically involved in vertebral fracture. These outcomes provide a strong biomechanical rationale for clinical methods, such as algorithm-based qualitative (ABQ) assessment, that diagnose vertebral fracture on the basis of endplate depression. © 2015 American Society for Bone and Mineral Research.

Elastic Anisotropy of Trabecular Bone in the Elderly Human Vertebra.

Knowledge of the nature of the elastic symmetry of trabecular bone is fundamental to the study of bone adaptation and failure. Previous studies have classified human vertebral trabecular bone as orthotropic or transversely isotropic but have typically obtained samples from only selected regions of the centrum. In this study, the elastic symmetry of human vertebral trabecular bone was characterized using microfinite element (µFE) analyses performed on 1019 cubic regions of side length equal to 5 mm, obtained via thorough sampling of the centrums of 18 human L1 vertebrae (age = 81.17 ± 7.7 yr; eight males and ten females). An optimization procedure was used to find the closest orthotropic representation of the resulting stiffness tensor for each cube. The orthotropic elastic constants and orientation of the principal elastic axes were then recorded for each cube and were compared to the constants predicted from Cowin's fabric-based constitutive model (Cowin, 1985, "The Relationship Between the Elasticity Tensor and the Fabric Tensor," Mech. Mater., 4(2), pp. 137-147.) and the orientation of the principal axes of the fabric tensor, respectively. Deviations from orthotropy were quantified by the "orthotropic error" (van Rietbergen et al., 1996, "Direct Mechanics Assessment of Elastic Symmetries and Properties of Trabecular Bone Architecture," J. Biomech., 29(12), pp. 1653-1657), and deviations from transverse isotropy were determined by statistical comparison of the secondary and tertiary elastic moduli. The orthotropic error was greater than 50% for nearly half of the cubes, and the secondary and tertiary moduli differed from one another (p < 0.0001). Both the orthotropic error and the difference between secondary and tertiary moduli decreased with increasing bone volume fraction (BV/TV; p ≤ 0.007). Considering only the cubes with an orthotropic error less than 50%, only moderate correlations were observed between the fabric-based and the µFE-computed elastic moduli (R2 ≥ 0.337; p < 0.0001). These results indicate that when using a criterion of 5 mm for a representative volume element (RVE), transverse isotropy or orthotropy cannot be assumed for elderly human vertebral trabecular bone. Particularly at low values of BV/TV, this criterion does not ensure applicability of theories of continuous media. In light of the very sparse and inhomogeneous microstructure found in the specimens analyzed in this study, further work is needed to establish guidelines for selecting a RVE within the aged vertebral centrum.

Endplate deflection is a defining feature of vertebral fracture and is associated with properties of the underlying trabecular bone.

Endplate deflection frequently occurs with vertebral failure, but the relationship between the two remains poorly defined. This study examined associations between endplate deflection under compressive loading and characteristics of the neighboring subchondral bone and intervertebral disc (IVD). Ten L1 vertebrae with adjacent IVDs were dissected, compressed axially in a stepwise manner to failure, and imaged with micro-computed tomography before each loading step. From the images, deflection was measured across the surface of each endplate at each step. Trabecular microstructure and endplate volume fraction were evaluated in 5 mm regions just under the superior endplate. IVDs were assessed using computed tomography and histology. A marked increase in superior endplate deflection coincided with a drop in the load-displacement curve. Endplate deflection was higher in regions with less robust bone microstructure (p < 0.009), though these associations tended to weaken as loading progressed. Immediately following the ultimate point, endplate deflection was higher in regions underlying the nucleus pulposus versus annulus fibrosus (p = 0.035), irrespective of disc grade (p = 0.346). These results indicate that a sudden increase in endplate deflection signals that the mechanical competence of the vertebra has been compromised. The mechanisms of endplate failure likely relate to anatomical features of the endplate, neighboring trabecular bone, and IVD.

Effect of specimen-specific anisotropic material properties in quantitative computed tomography-based finite element analysis of the vertebra.

Intra- and inter-specimen variations in trabecular anisotropy are often ignored in quantitative computed tomography (QCT)-based finite element (FE) models of the vertebra. The material properties are typically estimated solely from local variations in bone mineral density (BMD), and a fixed representation of elastic anisotropy ("generic anisotropy") is assumed. This study evaluated the effect of incorporating specimen-specific, trabecular anisotropy on QCT-based FE predictions of vertebral stiffness and deformation patterns. Orthotropic material properties estimated from microcomputed tomography data ("specimen-specific anisotropy"), were assigned to a large, columnar region of the L1 centrum (n = 12), and generic-anisotropic material properties were assigned to the remainder of the vertebral body. Results were compared to FE analyses in which generic-anisotropic properties were used throughout. FE analyses were also performed on only the columnar regions. For the columnar regions, the axial stiffnesses obtained from the two categories of material properties were uncorrelated with each other (p = 0.604), and the distributions of minimum principal strain were distinctly different (p ≤ 0.022). In contrast, for the whole vertebral bodies in both axial and flexural loading, the stiffnesses obtained using the two categories of material properties were highly correlated (R2 > 0.82, p < 0.001) with, and were no different (p > 0.359) from, each other. Only moderate variations in strain distributions were observed between the two categories of material properties. The contrasting results for the columns versus vertebrae indicate a large contribution of the peripheral regions of the vertebral body to the mechanical behavior of this bone. In companion analyses on the effect of the degree of anisotropy (DA), the axial stiffnesses of the trabecular column (p < 0.001) and vertebra (p = 0.007) increased with increasing DA. These findings demonstrate the need for accurate modeling of the peripheral regions of the vertebral body in analyses of the mechanical behavior of the vertebra.

Presence of intervertebral discs alters observed stiffness and failure mechanisms in the vertebra.

Ex vivo mechanical testing is an essential tool for study of vertebral mechanics. However, the common method of testing vertebral bodies in the absence of adjacent intervertebral discs (IVDs) may limit the physiological relevance of the results. The goal of this study was to determine the influence of IVDs on vertebral mechanical properties and failure mechanisms. Rabbit thoracic vertebral bodies were tested with and without IVDs in a stepwise fashion that incorporated a micro-computed tomography scan at each loading step. The image sequences were analyzed using digital volume correlation to quantify deformations throughout the vertebral body. The observed deformation patterns differed substantially between the groups. Specimens tested with IVDs exhibited a slow increase in strain in the inferior and posterior regions, followed by a sudden increase in strain in the anterior cortex right at the yield point. In contrast, the highest strains in the isolated vertebral bodies were in the posterior regions throughout the test. Specimens tested with IVDs had lower stiffness (507.49±184.73N/mm vs. 845.61±296.09N/mm; p=0.044), higher ultimate displacement (2.00±0.68mm vs. 1.17±0.54mm; p=0.043), and higher maximum shear strains (e.g. top 25th percentile: 0.19±0.11 vs. 0.06±0.07mm/mm; p<0.0458), and tended to have lower ultimate force (690.28±160.25N vs. 873.81±131.48N; p=0.056). Similar work to failure (648.15±317.86N-mm vs. 603.49±437.95 N-mm; p=0.844) was observed between the two groups. These results indicate that testing vertebral bodies in the absence of IVDs can elicit artifactual failure mechanisms. These artifacts may be more prominent than the effects on vertebral strength and toughness.

Vascular development during distraction osteogenesis proceeds by sequential intramuscular arteriogenesis followed by intraosteal angiogenesis.

Vascular formation is intimately associated with bone formation during distraction osteogenesis (DO). While prior studies on this association have focused on vascular formation locally within the regenerate, we hypothesized that this vascular formation, as well as the resulting osteogenesis, relies heavily on the response of the vascular network in surrounding muscular compartments. To test this hypothesis, the spatiotemporal sequence of vascular formation was assessed in both muscular and osseous compartments in a murine model of DO and was compared to the progression of osteogenesis. Micro-computed tomography (µCT) scans were performed sequentially, before and after demineralization, on specimens containing contrast-enhanced vascular casts. Image registration and subtraction procedures were developed to examine the co-related, spatiotemporal patterns of vascular and osseous tissue formation. Immunohistochemistry was used to assess the contributory roles of arteriogenesis (formation of large vessels) and angiogenesis (formation of small vessels) to overall vessel formation. Mean vessel thickness showed an increasing trend during the period of active distraction (p=0.068), whereas vessel volume showed maximal increases during the consolidation period (p=0.009). The volume of mineralized tissue in the regenerate increased over time (p<0.039), was correlated with vessel volume (r=0.59; p=0.025), and occurred primarily during consolidation. Immunohistological data suggested that: 1) the period of active distraction was characterized primarily by arteriogenesis in the surrounding muscle; 2) during consolidation, angiogenesis predominated in the intraosteal region; and 3) vessel formation proceeded from the surrounding muscle into the regenerate. These data show that formation of vascular tissue occurs in both muscular and osseous compartments during DO and that periods of intense osteogenesis are concurrent with those of angiogenesis. The results further suggest the presence of morphogenetic factors that coordinate the development of vascular tissues from the intramuscular compartment into the regions of osseous regeneration.