The I-PREDICT 50th Percentile Male Warfighter Finite Element Model: Development and Validation of the Thoracolumbar Spine

Military personnel are commonly at risk of lower back pain and thoracolumbar spine injury. Human volunteers and postmortem human subjects have been used to understand the scenarios where injury can occur and the tolerance of the warfighter to these loading regimes. Finite element human body models (HBMs) can accurately simulate the mechanics of the human body and are a useful tool for understanding injury. In this study, a HBM thoracolumbar spine was developed and hierarchically validated as part of the Incapacitation Prediction for Readiness in Expeditionary Domains: an Integrated Computational Tool (I-PREDICT) program. Constitutive material models were sourced from literature and the vertebrae and intervertebral discs were hexahedrally meshed from a 50th percentile male CAD dataset. Ligaments were modeled through attaching beam elements at the appropriate anatomical insertion sites. 94 simulations were replicated from experimental PMHS tests at the vertebral body, functional spinal unit (FSU), and regional lumbar spine levels. The BioRank (BRS) biofidelity ranking system was used to assess the response of the I-PREDICT model. At the vertebral body level, the I-PREDICT model showed good agreement with experimental results. The I-PREDICT FSUs showed good agreement in tension and compression and had comparable stiffness values in flexion, extension, and axial rotation. The regional lumbar spine exhibited "good" biofidelity when tested in tension, compression, extension, flexion, posterior shear, and anterior shear (BRS regional average = 1.05). The validated thoracolumbar spine of the I-PREDICT model can be used to better understand and mitigate injury risk to the warfighter.

Super-resolution of clinical CT: Revealing microarchitecture in whole bone clinical CT image data

Osteoporotic fractures, prevalent in the elderly, pose a significant health and economic burden. Current methods for predicting fracture risk, primarily relying on bone mineral density, provide only modest accuracy. If better spatial resolution of trabecular bone in a clinical scan were available, a more complete assessment of fracture risk would be obtained using microarchitectural measures of bone (i.e. trabecular thickness, trabecular spacing, bone volume fraction, etc.). However, increased resolution comes at the cost of increased radiation or can only be applied at small volumes of distal skeletal locations. This study explores super-resolution (SR) technology to enhance clinical CT scans of proximal femurs and better reveal the trabecular microarchitecture of bone. Using a deep-learning-based (i.e. subset of artificial intelligence) SR approach, low-resolution clinical CT images were upscaled to higher resolution and compared to corresponding MicroCT-derived images. SR-derived 2-dimensional microarchitectural measurements, such as degree of anisotropy, bone volume fraction, trabecular spacing, and trabecular thickness were within 16 % error compared to MicroCT data, whereas connectivity density exhibited larger error (as high as 1094 %). SR-derived 3-dimensional microarchitectural metrics exhibited errors <18 %. This work showcases the potential of SR technology to enhance clinical bone imaging and holds promise for improving fracture risk assessments and osteoporosis detection. Further research, including larger datasets and refined techniques, can advance SR's clinical utility, enabling comprehensive microstructural assessment across whole bones, thereby improving fracture risk predictions and patient-specific treatment strategies.

Mobility Solutions After a Lower Extremity Fracture and Applicability to Battlefield and Wilderness Medicine

The potential for delayed evacuation of injured Service members from austere environments highlights the need to develop solutions that can stabilize a wound and enable mobility during these prolonged casualty care (PCC) scenarios. Lower extremity fractures have traditionally been treated by immobilization (splinting) followed by air evacuation - a paradigm not practical in PCC scenarios. In the civilian sector, treatment of extremity injuries sustained during remote recreational activities have similar challenges, particularly when adverse weather or terrain precludes early ground or air rescue. This review examines currently available fracture treatment solutions to include splinting, orthotic devices, and biological interventions and evaluates their feasibility: 1) for prolonged use in austere environments and 2) to enable patient mobilization. This review returned three common types of splints to include: a simple box splint, pneumatic splints, and traction splints. None of these splinting techniques allowed for ambulation. However, fixed facility-based orthotic interventions that include weight-bearing features may be combined with common splinting techniques to improve mobility. Biologically-focused technologies to stabilize a long bone fracture are still in their infancy. Integrating design features across these technologies could generate advanced treatments which would enable mobility, thus maximizing survivability until patient evacuation is feasible.

Biomechanical characteristics of proficient free-throw shooters-markerless motion capture analysis

The winning game outcome in basketball is partially contingent on the team's ability to secure and make more free-throw shooting attempts, especially close to the end of the game. Thus, the purpose of the present study was to perform a comprehensive biomechanical analysis of the free-throw shooting motion to examine differences between (a) proficient (≥70%) and non-proficient shooters (<70%) and (b) made and missed free-throw shoots within the proficient group of shooters. Thirty-four recreationally active males with previous basketball playing experience attempted ten consecutive free-throw shots (4.57 m), with a 10-15 s rest interval between each shot. An innovative three-dimensional markerless motion capture system (SwRI Enable, San Antonio, TX, USA) composed of nine high-definition cameras recording at 120 Hz was used to capture and analyze the biomechanical parameters of interest. Independent t-tests and Mann-Whitney U tests were used to examine a presence of statistically significant differences. The findings of the present study reveal that proficient free-throw shooters performed the shooting motion in a more controlled manner by having significantly lower knee and center of mass peak and mean angular velocities. Also, proficient shooters attained a significantly greater release height and had less forward trunk lean when compared to non-proficient shooters at the time point of the ball release. Moreover, despite being beneficial for improvements in shooting accuracy, our findings suggest that overemphasizing the release height may be in certain instances counterproductive, as it may lead to more missed than made free-throw shots within the proficient group of shooters.

Response of Thoraco-Abdominal Tissue in High-Rate Compression

Body armor protects the human from penetrating injuries. However, in the process of defeating a projectile, the back face of the armor can deform into the wearer at extremely high rates. This deformation can cause a variety of soft and hard tissue injuries. Finite element modeling represents one of the best tools to predict injuries from this high-rate compression mechanism. However, the validity of a model is reliant on accurate material properties for biological tissues. In this study, we measured the stress-strain response of thoraco-abdominal tissue during high-rate compression (1000 and 1900 s-1) using a split Hopkinson pressure bar. Using this method, high-rate material properties of porcine adipose, heart, spleen, and stomach tissue were characterized. At a strain rate of 1000 s-1, adipose (E=4.7MPa) was the most compliant, followed by spleen (E=9.6MPa), and then heart (E=13.6MPa) tissue. At a strain rate of 1900 s-1, adipose (E=7.3MPa) was most compliant, followed by spleen (E=10.7MPa), heart (E=14.1MPa), and stomach (E=32.6MPa) tissue. Only adipose tissue demonstrated a consistent rate dependence across these rates, with a stiffer response at 1900 s-1. However, comparison of these tissues to previously published quasi-static and intermediate dynamic experiments revealed a strong rate dependence with increasing stress response from quasi-static to dynamic to high strain rates. Together, these findings can be used to develop finite element models for high-rate compression injuries.

Biglycan and chondroitin sulfate play pivotal roles in bone toughness via retaining bound water in bone mineral matrix

Recent in vitro evidence shows that glycosaminoglycans (GAGs) and proteoglycans (PGs) in bone matrix may functionally be involved in the tissue-level toughness of bone. In this study, we showed the effect of biglycan (Bgn), a small leucine-rich proteoglycan enriched in extracellular matrix of bone and the associated GAG subtype, chondroitin sulfate (CS), on the toughness of bone in vivo, using wild-type (WT) and Bgn deficient mice. The amount of total GAGs and CS in the mineralized compartment of Bgn KO mouse bone matrix decreased significantly, associated with the reduction of the toughness of bone, in comparison with those of WT mice. However, such differences between WT and Bgn KO mice diminished once the bound water was removed from bone matrix. In addition, CS was identified as the major subtype in bone matrix. We then supplemented CS to both WT and Bgn KO mice to test whether supplemental GAGs could improve the tissue-level toughness of bone. After intradermal administration of CS, the toughness of WT bone was greatly improved, with the GAGs and bound water amount in the bone matrix increased, while such improvement was not observed in Bgn KO mice or with supplementation of dermatan sulfate (DS). Moreover, CS supplemented WT mice exhibited higher bone mineral density and reduced osteoclastogenesis. Interestingly, Bgn KO bone did not show such differences irrespective of the intradermal administration of CS. In summary, the results of this study suggest that Bgn and CS in bone matrix play a pivotal role in imparting the toughness to bone most likely via retaining bound water in bone matrix. Moreover, supplementation of CS improves the toughness of bone in mouse models.

Validity and reliability of an accelerometer-based player tracking device

This study aimed to determine the intra- and inter-device accuracy and reliability of wearable athletic tracking devices, under controlled laboratory conditions. A total of nineteen portable accelerometers (Catapult OptimEye S5) were mounted to an aluminum bracket, bolted directly to an Unholtz Dickie 20K electrodynamic shaker table, and subjected to a series of oscillations in each of three orthogonal directions (front-back, side to side, and up-down), at four levels of peak acceleration (0.1g, 0.5g, 1.0g, and 3.0g), each repeated five times resulting in a total of 60 tests per unit, for a total of 1140 records. Data from each accelerometer was recorded at a sampling frequency of 100Hz. Peak accelerations recorded by the devices, Catapult PlayerLoad™, and calculated player load (using Catapult’s Cartesian formula) were used for the analysis. The devices demonstrated excellent intradevice reliability and mixed interdevice reliability. Differences were found between devices for mean peak accelerations and PlayerLoad™ for each direction and level of acceleration. Interdevice effect sizes ranged from a mean of 0.54 (95% CI: 0.34–0.74) (small) to 1.20 (95% CI: 1.08–1.30) (large) and ICCs ranged from 0.77 (95% CI: 0.62–0.89) (very large) to 1.0 (95% CI: 0.99–1.0) (nearly perfect) depending upon the magnitude and direction of the applied motion. When compared to the player load determined using the Cartesian formula, the Catapult reported PlayerLoad™ was consistently lower by approximately 15%. These results emphasize the need for industry wide standards in reporting validity, reliability and the magnitude of measurement errors. It is recommended that device reliability and accuracy are periodically quantified.

Determination of sex differences of human cadaveric mandibular condyles using statistical shape and trait modeling

The objective of this study was to elucidate sex differences in the anatomy of human temporomandibular joint mandibular condyles using a statistical shape and trait model (SSTM). Mandibles were obtained from 16 human cadavers (79±13years). The condyles were scanned using micro-computed tomography with 27μm resolution. An image processing algorithm was used to segment the bone, determine the border of the entire mandibular condyle and trabecular bone compartments, and create triangulated meshes of the compartments. One subject was chosen as the template and was registered to the other individuals using a coherence point drift algorithm. This process positioned all vertices at corresponding anatomic locations. For the trabecular bone region, around each vertex position, the average bone image intensity, which is proportional to bone density, and microstructural traits, including trabecular bone volume fraction, thickness, separation, connectivity, and connectivity density were calculated. For the entire mandibular condyle mesh, the surface vertices were extracted to represent the overall anatomy of the condyle. Using a SSTM, the shape and trait information was reduced to a small set of independent and uncorrelated variables for each individual. Wilcoxon rank sum tests were used to test for differences in the variables between sexes. A lasso approach was used to determine a set of variables that differentiate between sexes. Male condyles were on average larger than female condyles, with complex differences in the microstructural traits. Two out of 15 principal components were statistically different between males and females (p<0.1). The lasso approach determined a set of 7 principal components that fully described the complex shape and trait differences between males and females. A SSTM was able to determine sex-dependent differences in the shape of the mandibular condyle. These differences may alter the biomechanics of the joint and contribute to the development of temporomandibular joint disease.

Sex dependent mechanical properties of the human mandibular condyle

The mandibular condyle consists of articular cartilage and subchondral bone that play an important role in bearing loads at the temporomandibular joint (TMJ) during static occlusion and dynamic mastication. The objective of the current study was to examine effects of sex and cartilage on 1) static and dynamic mechanical analysis (DMA) based dynamic energy storage and dissipation for the cartilage-subchondral bone construct of the human mandibular condyle, and 2) their correlations with the tissue mineral density and trabecular morphological parameters of subchondral bone. Cartilage-subchondral bone constructs were obtained from 16 individual human cadavers (9 males, 7 females, 79.00±13.10 years). After scanning with micro-computed tomography, the specimens were subjected to a non-destructive compressive static loading up to 7N and DMA using a cyclic loading profile (-5±2N at 2Hz). After removing the cartilage from the same specimen, the series of loading experiments were repeated. Static stiffness (K) and energy dissipation (W), and dynamic storage (K'), loss (K'') stiffness, and energy dissipation (tan δ) were assessed. Gray values, which are proportional to degree of bone mineralization, and trabecular morphological parameters of the subchondral bone were also measured. After removal of the cartilage, static energy dissipation significantly decreased (p<0.009) but dynamic energy dissipation was not influenced (p>0.064). Many subchondral bone properties were significantly correlated with the overall mechanical behavior of the cartilage-subchondral bone constructs for males (p<0.047) but not females (p>0.054). However, after removal of cartilage from the constructs, all of the significant correlations were no longer found (p>0.057). The current findings indicate that the subchondral bone is responsible for bearing static and dynamic loading in males but not in females. This result indicates that the female condyle may have a mechanically disadvantageous TMJ loading environment.

Finite Element Study of a Lumbar Intervertebral Disc Nucleus Replacement Device

Nucleus replacement technologies are a minimally invasive alternative to spinal fusion and total disc replacement that have the potential to reduce pain and restore motion for patients with degenerative disc disease. Finite element modeling can be used to determine the biomechanics associated with nucleus replacement technologies. The current study focuses on a new nucleus replacement device designed as a conforming silicone implant with an internal void. A validated finite element model of the human lumbar L3-L4 motion segment was developed and used to investigate the influence of the nucleus replacement device on spine biomechanics. In addition, the effect of device design changes on biomechanics was determined. A 3D, L3-L4 finite element model was constructed from medical imaging data. Models were created with the normal intact nucleus, the nucleus replacement device, and a solid silicone implant. Probabilistic analysis was performed on the normal model to provide quantitative validation metrics. Sensitivity analysis was performed on the silicone Shore A durometer of the device. Models were loaded under axial compression followed by flexion/extension, lateral bending, or axial rotation. Compressive displacement, endplate stresses, reaction moment, and annulus stresses were determined and compared between the different models. The novel nucleus replacement device resulted in similar compressive displacement, endplate stress, and annulus stress and slightly higher reaction moment compared with the normal nucleus. The solid implant resulted in decreased displacement, increased endplate stress, decreased annulus stress, and decreased reaction moment compared with the novel device. With increasing silicone durometer, compressive displacement decreased, endplate stress increased, reaction moment increased, and annulus stress decreased. Finite element analysis was used to show that the novel nucleus replacement device results in similar biomechanics compared with the normal intact nucleus.

Development and validation of a statistical shape modeling-based finite element model of the cervical spine under low-level multiple direction loading conditions

Cervical spinal injuries are a significant concern in all trauma injuries. Recent military conflicts have demonstrated the substantial risk of spinal injury for the modern warfighter. Finite element models used to investigate injury mechanisms often fail to examine the effects of variation in geometry or material properties on mechanical behavior. The goals of this study were to model geometric variation for a set of cervical spines, to extend this model to a parametric finite element model, and, as a first step, to validate the parametric model against experimental data for low-loading conditions. Individual finite element models were created using cervical spine (C3–T1) computed tomography data for five male cadavers. Statistical shape modeling (SSM) was used to generate a parametric finite element model incorporating variability of spine geometry, and soft-tissue material property variation was also included. The probabilistic loading response of the parametric model was determined under flexion-extension, axial rotation, and lateral bending and validated by comparison to experimental data. Based on qualitative and quantitative comparison of the experimental loading response and model simulations, we suggest that the model performs adequately under relatively low-level loading conditions in multiple loading directions. In conclusion, SSM methods coupled with finite element analyses within a probabilistic framework, along with the ability to statistically validate the overall model performance, provide innovative and important steps toward describing the differences in vertebral morphology, spinal curvature, and variation in material properties. We suggest that these methods, with additional investigation and validation under injurious loading conditions, will lead to understanding and mitigating the risks of injury in the spine and other musculoskeletal structures.

Fracture risk predictions based on statistical shape and density modeling of the proximal femur

Increased risk of skeletal fractures due to bone mass loss is a major public health problem resulting in significant morbidity and mortality, particularly in the case of hip fractures. Current clinical methods based on two-dimensional measures of bone mineral density (areal BMD or aBMD) are often unable to identify individuals at risk of fracture. We investigated predictions of fracture risk based on statistical shape and density modeling (SSDM) methods using a case-cohort sample of individuals from the Osteoporotic Fractures in Men (MrOS) study. Baseline quantitative computed tomography (QCT) data of the right femur were obtained for 513 individuals, including 45 who fractured a hip during follow-up (mean 6.9 year observation, validated by physician review). QCT data were processed for 450 individuals (including 40 fracture cases) to develop individual models describing three-dimensional bone geometry and density distribution. Comparison of mean fracture and non-case models indicated complex structural differences that appear to be responsible for resistance to hip fracture. Logistic regressions were used to model the relation of baseline hip BMD and SSDM weighting factors to the occurrence of hip fracture. Area under the receiver operating characteristic (ROC) curve (AUC) for a prediction model based on weighting factors and adjusted by age was significantly greater than AUC for a prediction model based on aBMD and age (0.94 versus 0.83, respectively). The SSDM-based prediction model adjusted by age correctly identified 55% of the fracture cases (and 94.7% of the non-cases), whereas the clinical standard aBMD correctly identified 10% of the fracture cases (and 91.3% of the non-cases). SSDM identifies subtle changes in combinations of structural bone traits (eg, geometric and BMD distribution traits) that appear to indicate fracture risk. Investigation of important structural differences in the proximal femur between fracture and no-fracture cases may lead to improved prediction of those at risk for future hip fracture.

Baboons as a model to study genetics and epigenetics of human disease

A major challenge for understanding susceptibility to common human diseases is determining genetic and environmental factors that influence mechanisms underlying variation in disease-related traits. The most common diseases afflicting the US population are complex diseases that develop as a result of defects in multiple genetically controlled systems in response to environmental challenges. Unraveling the etiology of these diseases is exceedingly difficult because of the many genetic and environmental factors involved. Studies of complex disease genetics in humans are challenging because it is not possible to control pedigree structure and often not practical to control environmental conditions over an extended period of time. Furthermore, access to tissues relevant to many diseases from healthy individuals is quite limited. The baboon is a well-established research model for the study of a wide array of common complex diseases, including dyslipidemia, hypertension, obesity, and osteoporosis. It is possible to acquire tissues from healthy, genetically characterized baboons that have been exposed to defined environmental stimuli. In this review, we describe the genetic and physiologic similarity of baboons with humans, the ability and usefulness of controlling environment and breeding, and current genetic and genomic resources. We discuss studies on genetics of heart disease, obesity, diabetes, metabolic syndrome, hypertension, osteoporosis, osteoarthritis, and intrauterine growth restriction using the baboon as a model for human disease. We also summarize new studies and resources under development, providing examples of potential translational studies for targeted interventions and therapies for human disease.

Characterization of complex, co-adapted skeletal biomechanics phenotypes: a needed paradigm shift in the genetics of bone structure and function

The genetic architecture of skeletal biomechanical performance has tremendous potential to advance our knowledge of the biological mechanisms that drive variation in skeletal fragility and osteoporosis risk. Research using traditional approaches that focus on specific gene pathways is increasing our understanding of how and to what degree those pathways may affect population-level variation in fracture susceptibility, and shows that known pathways may affect bone fragility through unsuspected mechanisms. Non-traditional approaches that incorporate a new appreciation for the degree to which bone traits co-adapt to functional loading environments, using a wide variety of redundant compensatory mechanisms to meet both physiological and mechanical demands, represent a radical departure from the dominant reductionist paradigm and have the potential to rapidly advance our understanding of bone fragility and identification of new targets for therapeutic intervention.

Platelet-derived growth-factor-releasing aligned collagen-nanoparticle fibers promote the proliferation and tenogenic differentiation of adipose-derived stem cells

In order to enhance the healing potential of an injured tendon, we have prepared a novel biomimetic aligned collagen-nanoparticle (NP) composite fiber using an electrochemical process. The aligned collagen-NP composite fiber is designed to affect the cellular activity of adipose-derived stem cells (ADSCs) through two different ways: (i) topographic cues from the alignment of collagen fibril and (ii) controlled release of platelet-derived growth factors (PDGFs) from the NPs. PDGF released from collagen-NP fibers significantly enhanced the proliferation of ADSCs when tested for up to 7 days. Moreover, compared to random collagen fibers with PDGFs, aligned collagen-NP fibers significantly promoted the desirable tenogenic differentiation of ADSCs, as evidenced by an increased level of tendon markers such as tenomodulin and scleraxis. On the other hand, no undesirable osteogenic differentiation, as measured by the unchanged level of alkaline phosphatase and osteocalcin, was observed. Together, these results indicate that the aligned collagen-NP composite fiber induced the tenogenic differentiation of ADSCs through both a topographic cue (aligned collagen fibril) and a chemical cue (PDGF released from NPs). Thus, our novel aligned collagen-NP composite fiber has a significant potential to be used for tendon tissue engineering and regeneration.

Intracortical bone remodeling variation shows strong genetic effects

Intracortical microstructure influences crack propagation and arrest within bone cortex. Genetic variation in intracortical remodeling may contribute to mechanical integrity and, therefore, fracture risk. Our aim was to determine the degree to which normal population-level variation in intracortical microstructure is due to genetic variation. We examined right femurs from 101 baboons (74 females, 27 males; aged 7-33 years) from a single, extended pedigree to determine osteon number, osteon area (On.Ar), haversian canal area, osteon population density, percent osteonal bone (%On.B), wall thickness (W.Th), and cortical porosity (Ct.Po). Through evaluation of the covariance in intracortical properties between pairs of relatives, we quantified the contribution of additive genetic effects (heritability [h (2)]) to variation in these traits using a variance decomposition approach. Significant age and sex effects account for 9 % (Ct.Po) to 21 % (W.Th) of intracortical microstructural variation. After accounting for age and sex, significant genetic effects are evident for On.Ar (h (2) = 0.79, p = 0.002), %On.B (h (2) = 0.82, p = 0.003), and W.Th (h (2) = 0.61, p = 0.013), indicating that 61-82 % of the residual variation (after accounting for age and sex effects) is due to additive genetic effects. This corresponds to 48-75 % of the total phenotypic variance. Our results demonstrate that normal, population-level variation in cortical microstructure is significantly influenced by genes. As a critical mediator of crack behavior in bone cortex, intracortical microstructural variation provides another mechanism through which genetic variation may affect fracture risk.

Computational anatomy in the study of bone structure

Osteoporosis is a major public health threat for millions of Americans with billions of dollars per year of national direct costs for osteoporotic fractures. Osteoporosis results in a decrease in overall bone mass and subsequent increase in the risk of bone fracture. Bone strength arises from the combination of bone size and shape, the distribution of bone mass throughout the structure, and the quality of the bone material. Advances in medical imaging have enabled a comprehensive assessment of bone structure through the analysis of high-resolution scans of relevant anatomical sites, eg, the proximal femur. However, conventional imaging analysis techniques use predefined regions of interest that do not take full advantage of such scans. Recently, computational anatomy, a set of imaging-based analysis algorithms, has emerged as a promising technique in studies of osteoporosis. Computational anatomy enables analyses that are not biased to one particular region and provide a more complete assessment of the whole structure. In this article, we review studies that have used computational anatomy to investigate the structure of the proximal femur in relation to age, fracture, osteoporotic treatment, and spaceflight effects.

Measurement and estimation of osteocyte mechanical strain

Osteocytes are the most abundant cell type in bone and are responsible for sensing mechanical strain and signaling bone (re)modeling, making them the primary mechanosensors within bone. Under aging and osteoporotic conditions, bone is known to be less responsive to loading (exercise), but it is unclear why. Perhaps, the levels of mechanical strain required to initiate these biological events are not perceived by the osteocytes embedded within the bone tissue. In this review we examine the methods used to measure and estimate the strains experienced by osteocytes in vivo as well as the results of related published experiments. Although the physiological levels of strain experienced by osteocytes in vivo are still under investigation, through computational modeling and laboratory experiments, it has been shown that there is significant amplification of average bone strain at the level of the osteocyte lacunae. It has also been proposed that the material properties of the perilacunar region surrounding the osteocyte can have significant effects of the strain perceived by the embedded osteocyte. These facts have profound implications for studies involving osteoporotic bone where the material properties are known to become stiffer.

Addressing the gaps: sex differences in osteoarthritis of the knee

An introduction to the accompanying three papers.

Hormonal modulation of connective tissue homeostasis and sex differences in risk for osteoarthritis of the knee

Young female athletes experience a higher incidence of ligament injuries than their male counterparts, females experience a higher incidence of joint hypermobility syndrome (a risk factor for osteoarthritis development), and post-menopausal females experience a higher prevalence of osteoarthritis than age-matched males. These observations indicate that fluctuating sex hormone levels in young females and loss of ovarian sex hormone production due to menopause likely contribute to observed sex differences in knee joint function and risk for loss of function. In studies of osteoarthritis, however, there is a general lack of appreciation for the heterogeneity of hormonal control in both women and men. Progress in this field is limited by the relatively few preclinical osteoarthritis models, and that most of the work with established models uses only male animals. To elucidate sex differences in osteoarthritis, it is important to examine sex hormone mechanisms in cells from knee tissues and the sexual dimorphism in the role of inflammation at the cell, tissue, and organ levels. There is a need to determine if the risk for loss of knee function and integrity in females is restricted to only the knee or if sex-specific changes in other tissues play a role. This paper discusses these gaps in knowledge and suggests remedies.

Mechanical contributors to sex differences in idiopathic knee osteoarthritis

The occurrence of knee osteoarthritis (OA) increases with age and is more common in women compared with men, especially after the age of 50 years. Recent work suggests that contact stress in the knee cartilage is a significant predictor of the risk for developing knee OA. Significant gaps in knowledge remain, however, as to how changes in musculoskeletal traits disturb the normal mechanical environment of the knee and contribute to sex differences in the initiation and progression of idiopathic knee OA. To illustrate this knowledge deficit, we summarize what is known about the influence of limb alignment, muscle function, and obesity on sex differences in knee OA. Observational data suggest that limb alignment can predict the development of radiographic signs of knee OA, potentially due to increased stresses and strains within the joint. However, these data do not indicate how limb alignment could contribute to sex differences in either the development or worsening of knee OA. Similarly, the strength of the knee extensor muscles is compromised in women who develop radiographic and symptomatic signs of knee OA, but the extent to which the decline in muscle function precedes the development of the disease is uncertain. Even less is known about how changes in muscle function might contribute to the worsening of knee OA. Conversely, obesity is a stronger predictor of developing knee OA symptoms in women than in men. The influence of obesity on developing knee OA symptoms is not associated with deviation in limb alignment, but BMI predicts the worsening of the symptoms only in individuals with neutral and valgus (knock-kneed) knees. It is more likely, however, that obesity modulates OA through a combination of systemic effects, particularly an increase in inflammatory cytokines, and mechanical factors within the joint. The absence of strong associations of these surrogate measures of the mechanical environment in the knee joint with sex differences in the development and progression of knee OA suggests that a more multifactorial and integrative approach in the study of this disease is needed. We identify gaps in knowledge related to mechanical influences on the sex differences in knee OA.

Neural and psychosocial contributions to sex differences in knee osteoarthritic pain

People with osteoarthritis (OA) can have significant pain that interferes with function and quality of life. Women with knee OA have greater pain and greater reductions in function and quality of life than men. In many cases, OA pain is directly related to sensitization and activation of nociceptors in the injured joint and correlates with the degree of joint effusion and synovial thickening. In some patients, however, the pain does not match the degree of injury and continues after removal of the nociceptors with a total joint replacement. Growth of new nociceptors, activation of nociceptors in the subchondral bone exposed after cartilage degradation, and nociceptors innervating synovium sensitized by inflammatory mediators could all augment the peripheral input to the central nervous system and result in pain. Enhanced central excitability and reduced central inhibition could lead to prolonged and enhanced pain that does not directly match the degree of injury. Psychosocial variables can influence pain and contribute to pain variability. This review explores the neural and psychosocial factors that contribute to knee OA pain with an emphasis on differences between the sexes and gaps in knowledge.

Micromechanical modeling of R-curve behaviors in human cortical bone

The risk of bone fracture increases with age because of a variety of factors that include, among others, decreasing bone quantity and quality due to increasing porosity and crack density with age. Experimental evidence has indicated that changes in bone microstructure and trace mineralization with age can result in different crack-tip strain field and fracture response, leading to different fracture mechanisms and R-curve behaviors. In this paper, a micromechanical modeling approach is developed to predict the R-curve response of bone tissue by delineating fracture mechanisms that lead to microdamage and ligament bridging by incorporating the influence of increasing porosity and crack density with age. The effects of age on fracture of human femur cortical bone due to porosity (bone quantity) and bone quality (crack density) with age are then examined via the micromechanical model.

Load/strain distribution between ulna and radius in the mouse forearm compression loading model

Finite element analysis (FEA) of the mouse forearm compression loading model is used to relate strain distributions with downstream changes in bone formation and responses of bone cells. The objective of this study was to develop two FEA models - the first one with the traditional ulna only and the second one in which both the ulna and radius are included, in order to examine the effect of the inclusion of the radius on the strain distributions in the ulna. The entire mouse forearm was scanned using microCT and images were converted into FEA tetrahedral meshes using a suite of software programs. The performance of both linear and quadratic tetrahedral elements and coarse and fine meshes were studied. A load of 2N was applied to the ulna/radius model and a 1.3N load (based on previous investigations of load sharing between the ulna and radius in rats) was applied to the ulna only model for subsequent simulations. The results showed differences in the cross sectional strain distributions and magnitude within the ulna for the combined ulna/radius model versus the ulna only model. The maximal strain in the combined model occurred about 4mm toward the distal end from the ulna mid-shaft in both models. Results from the FEA model simulations were also compared to experimentally determined strain values. We conclude that inclusion of the radius in FE models to predict strains during in vivo forearm loading increases the magnitude of the estimated ulna strains compared to those predicted from a model of the ulna alone but the distribution was similar. This has important ramifications for future studies to understand strain thresholds needed to activate bone cell responses to mechanical loading.

Dark horse in osteocyte biology: Glycocalyx around the dendrites is critical for osteocyte mechanosensing

Osteocytes are considered as the major mechanosensory cells of the bone tissue that control the bone remodeling process. Since osteocytes are buried inside mineralized matrix, they maintain a strong communication network with other cells. Long dendritic processes of the osteocytes act as communication cables, conveying mechanical signals to the neighboring osteocytes and the cells on the bone surface; like osteoblasts and osteoclasts. Gap junctions and hemichannels formed by Connexin (Cx) 43 are observed to be involved in responding to the mechanical stimulus and in communicating the mechano-responsive biochemical signals. The contrast in the arrangement of the osteocyte cell body and the dendrites raises an important question of how these parts of the osteocyte respond to mechanical stimulation. We addressed this issue in our recent report through the stimulation of either osteocyte cell body or dendrites and our findings suggest that the osteocyte dendritic processes are sensitive to mechanical stimulation in comparison to the cell body. Most importantly, we observed that the dendritic processes are capable of conveying the mechanical signals to the cell body. Our findings also suggested that the glycocalyx surrounding the dendrites is required for sensing and conveying the mechanical signals. Degradation of the glycocalyx also leads to poor integrin attachment, thereby, affecting dendritic stiffness. These results suggest that the osteocyte dendritic processes are highly responsive towards mechanical stimulation and the glycocalyx surrounding the dendrites is critical in transducing these mechanical signals.

Development of a parametric finite element model of the proximal femur using statistical shape and density modelling

Skeletal fractures associated with bone mass loss are a major clinical problem and economic burden, and lead to significant morbidity and mortality in the ageing population. Clinical image-based measures of bone mass show only moderate correlative strength with bone strength. However, engineering models derived from clinical image data predict bone strength with significantly greater accuracy. Currently, image-based finite element (FE) models are time consuming to construct and are non-parametric. The goal of this study was to develop a parametric proximal femur FE model based on a statistical shape and density model (SSDM) derived from clinical image data. A small number of independent SSDM parameters described the shape and bone density distribution of a set of cadaver femurs and captured the variability affecting proximal femur FE strength predictions. Finally, a three-dimensional FE model of an 'unknown' femur was reconstructed from the SSDM with an average spatial error of 0.016 mm and an average bone density error of 0.037 g/cm(3).

Non-destructive characterization of microdamage in cortical bone using low field pulsed NMR

The microcracking and damage accumulation process in human cortical bone was characterized by performing cyclic loading under four-point bending at ambient temperature. A non-destructive nuclear magnetic resonance (NMR) spin-spin (T(2)) relaxation technique was applied to quantify the apparent changes in bone porosity as a function of cyclic loading and prior damage accumulation, first to unloaded cortical bone to quantify the initial porosity and then to fatigued cortical bone that was subjected to cyclic loading to various levels of modulus degradation and microdamage in the form of microcracks. The NMR T(2) relaxation time and amplitude data of the fatigued bone were compared against the undamaged state. The difference in the T(2) relaxation time data was taken as a measure of the increase in pore size, bone porosity or microcrack density due to microdamage induced by cyclic loading. A procedure was developed to deduce the number and size distributions of microcracks formed in cortical bone. Serial sectioning of the fatigued bone showed the formation of microcracks along the cement lines or within the interstitial tissue. The results on the evolution of microdamage derived from NMR measurements were verified by independent experimental measurements of microcrack density using histological characterization techniques. The size distribution and population of the microcracks were then utilized in conjunction with an analytical model to predict the degradation of the elastic modulus of cortical bone as a function of damage accumulation.

Higher doses of bisphosphonates further improve bone mass, architecture, and strength but not the tissue material properties in aged rats

We report the results of a series of experiments designed to determine the effects of ibandronate (Ibn) and risedronate (Ris) on a number of bone quality parameters in aged osteopenic rats to explain how bone material and bone mass may be affected by the dose of bisphosphonates (BP) and contribute to their anti-fracture efficacy. Eighteen-month old female rats underwent either ovariectomy or sham surgery. The ovariectomized (OVX) groups were left untreated for 2 months to develop osteopenia. Treatments started at 20 months of age as follows: sham and OVX control (treated with saline), OVX + risedronate 30 and 90 (30 or 90 microg/kg/dose), and OVX + ibandronate 30 and 90 (30 or 90 microg/kg/dose). The treatments were given monthly for 4 months by subcutaneous injection. At sacrifice at 24 months of age the 4th lumbar vertebra was used for microCT scans (bone mass, architecture, and degree of mineralization of bone, DMB) and histomorphometry, and the 6th lumbar vertebra, tibia, and femur were collected for biomechanical testing to determine bone structural and material strength, cortical fracture toughness, and tissue elastic modulus. The compression testing of the vertebral bodies (LVB6) was simulated using finite-element analysis (FEA) to also estimate the bone structural stiffness. Both Ibn and Ris dose-dependently increased bone mass and improved vertebral bone microarchitecture and mechanical properties compared to OVX control. Estimates of vertebral maximum stress from FEA were correlated with vertebral maximum load (r=0.5, p<0.001) and maximum stress (r=0.4, p<0.005) measured experimentally. Tibial bone bending modulus and cortical strength increased compared to OVX with both BP but no dose-dependent effect was observed. DMB and elastic modulus of trabecular bone were improved with Ibn 30 compared to OVX but were not affected in other BP-treated groups. DMB of tibial cortical bone showed no change with BP treatments. The fracture toughness examined in midshaft femurs did not change with BP even with the higher doses. In summary, the anti-fracture efficacy of BP is largely due to their preservation of bone mass and while the higher doses further improve the bone structural properties do not improve the localized bone material characteristics such as tissue strength, elastic modulus, and cortical toughness.

Statistical shape modeling describes variation in tibia and femur surface geometry between Control and Incidence groups from the osteoarthritis initiative database

We hypothesize that variability in knee subchondral bone surface geometry will differentiate between patients at risk and those not at risk for developing osteoarthritis (OA) and suggest that statistical shape modeling (SSM) methods form the basis for developing a diagnostic tool for predicting the onset of OA. Using a subset of clinical knee MRI data from the osteoarthritis initiative (OAI), the objectives of this study were to (1) utilize SSM to compactly and efficiently describe variability in knee subchondral bone surface geometry and (2) determine the efficacy of SSM and rigid body transformations to distinguish between patients who are not expected to develop osteoarthritis (i.e. Control group) and those with clinical risk factors for OA (i.e. Incidence group). Quantitative differences in femur and tibia surface geometry were demonstrated between groups, although differences in knee joint alignment measures were not statistically significant, suggesting that variability in individual bone geometry may play a greater role in determining joint space geometry and mechanics. SSM provides a means of explicitly describing complete articular surface geometry and allows the complex spatial variation in joint surface geometry and joint congruence between healthy subjects and those with clinical risk of developing or existing signs of OA to be statistically demonstrated.

Heritability of lumbar trabecular bone mechanical properties in baboons

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

Correlation of cell strain in single osteocytes with intracellular calcium, but not intracellular nitric oxide, in response to fluid flow

Osteocytes compose 90-95% of all bone cells and are the mechanosensors of bone. In this study, the strain experienced by individual osteocytes resulting from an applied fluid flow shear stress was quantified and correlated to two biological responses measured in real-time within the same individual osteocytes: (1) the upregulation of intracellular calcium and (2) changes in intracellular nitric oxide. Osteocyte-like MLO-Y4 cells were loaded with Fluo-4 AM and DAR-4M and exposed to uniform laminar fluid flow shear stresses of 2, 8, or 16 dyn/cm(2). Intracellular calcium and nitric oxide changes were determined by measuring the difference in fluorescence intensity from the cell's basal level prior to fluid flow and the level immediately following exposure. Individual cell strains were calculated using digital image correlation. MLO-Y4 cells showed a linear increase in cell strain, intracellular calcium concentration, and nitric oxide concentration with an increase in applied fluid flow rate. The increase in intracellular calcium was well correlated to the strain that each cell experienced. This study shows that osteocytes exposed to the same fluid flow experienced a range of individual strains and changes in intracellular calcium and nitric oxide concentrations, and the changes in intracellular calcium were correlated with cell strain. These results are among the first to establish a relationship between the strain experienced by osteocytes in response to fluid flow shear and a biological response at the single cell level. Mechanosensing and chemical signaling in osteocytes has been hypothesized to occur at the single cell level, making it imperative to understand the biological response of the individual cell.

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.

Relating crack-tip deformation to mineralization and fracture resistance in human femur cortical bone

The risk of bone fracture increases with age because of a variety of factors that include, among others, decreasing bone quantity and quality. Despite recent advances, the roles of bone microstructure and trace mineralization in the fracture process are not well understood. In this study, we utilize a combination of in-situ fracture toughness testing, digital strain mapping, and X-ray photoelectron spectroscopy techniques to characterize the near-tip strain field, fracture toughness, and chemical elements on the fracture surface of bone specimens from donors of two ages (48-year-old and 78-year-old females). We show that age-related embrittlement of bone fracture is associated with higher near-tip strains by lamellar shear and crack deflection at lamellar interfaces in the young bone and their absence in the old bone. The different near-tip deformation behaviors may be associated with the presence of Si and Zn in the young bone but more Ca and P and the lack of Si and Zn in the old bone.

Measurements of mobile and bound water by nuclear magnetic resonance correlate with mechanical properties of bone

Since clinical measures of bone mineral density do not necessarily predict whether a person will fracture a bone without an intervention, there is a need to find supplementary tools for assessing bone quality. Presently, we hypothesized that measures of mobile and bound water by a Nuclear Magnetic Resonance (NMR) technique are correlated with bone strength and toughness, respectively. To test this, bending specimens from the mid-diaphysis of 18 human femurs were collected from 18 male donors and divided into middle aged and elderly groups. After NMR measurements of each hydrated specimen, an inversion technique was used to convert the free induction decay data into a distribution of spin-spin (T2) relaxation rates. Then, the distribution resolved into three distinct components that likely represent solid hydrogen, water bound to bone tissue, and mobile water that occupy microscopic pores within the bone specimen. The integrated signal intensities of the bound and mobile components were normalized by the wet mass of the specimen. Following NMR measurements, three point bending tests were conducted to determine the modulus of elasticity, flexure strength, and work to fracture of each specimen. Next, the porosity, mineral-to-collagen ratio, and pentosidine concentration were measured. In this sample of human cortical bone, there was no age-related difference in the amount of mobile water, but the decrease in the amount of bound water with increasing age was statistically significant. Moreover, bound water was associated with both strength and work to fracture of bone, while mobile water was correlated with modulus of elasticity and appeared to quantify the level of microscopic pores within bone. On the other hand, bound water was correlated with the concentration of non-enzymatic collagen cross-links. The results of this study indicate that quantifying mobile and bound water with magnetic resonance techniques could potentially serve as indicators of bone quality.

Tissue strain amplification at the osteocyte lacuna: a microstructural finite element analysis

A parametric finite element model of an osteocyte lacuna was developed to predict the microstructural response of the lacuna to imposed macroscopic strains. The model is composed of an osteocyte lacuna, a region of perilacunar tissue, canaliculi, and the surrounding bone tissue. A total of 45 different simulations were modeled with varying canalicular diameters, perilacunar tissue material moduli, and perilacunar tissue thicknesses. Maximum strain increased with a decrease in perilacunar tissue modulus and decreased with an increase in perilacunar tissue modulus, regardless of the thickness of the perilacunar region. An increase in the predicted maximum strain was observed with an increase in canalicular diameter from 0.362 to 0.421 microm. In response to the macroscopic application of strain, canalicular diameters increased 0.8% to over 1.0% depending on the perilacunar tissue modulus. Strain magnification factors of over 3 were predicted. However, varying the size of the perilacunar tissue region had no effect on the predicted perilacunar tissue strain. These results indicate that the application of average macroscopic strains similar to strain levels measured in vivo can result in significantly greater perilacunar tissue strains and canaliculi deformations. A decrease in the perilacunar tissue modulus amplifies the perilacunar tissue strain and canaliculi deformation while an increase in the local perilacunar tissue modulus attenuates this effect.

The effect of three-dimensional shape optimization on the probabilistic response of a cemented femoral hip prosthesis

Probabilistic analyses allow the effect of uncertainty in system parameters on predicted model performance measures to be determined. Furthermore, using performance functions to describe a failure event, the probability of failure can be quantified. The effect of three-dimensional prosthesis shape optimization on the probabilistic response and failure probability of a cemented hip prosthesis system is investigated. Random variables include joint and muscle loading, cortical and cancellous bone and PMMA bone cement elastic properties, and strength parameters describing failure of the bone cement and the prosthesis-bone cement interface. Several performance functions describing the bone cement and prosthesis-cement interface are used to compute the probability of failure. When evaluated deterministically, most performance functions indicated a safe design, with the exception of interface tensile failure. However, when evaluated probabilistically, finite probabilities of failure were computed, some significant. The most likely mode of failure before shape optimization was prosthesis-bone cement interface tensile failure with a predicted probability of failure of 97.9%. Deterministic prosthesis shape optimization reduced the probability of failure for all performance functions and reduced prosthesis-bone cement interface tensile failure by 31.7%. Probability sensitivity factors indicate that the uncertainty in the joint loading, cement strength, and implant-cement interface strength have the greatest effect on the computed probability of failure. Implant shape optimization results in a more robust implant design that is less sensitive to uncertainties in joint loading, which cannot be easily controlled, and more sensitive to cement and interface properties, which are easier to modify.

Microstrain fields for cortical bone in uniaxial tension: Optical analysis method

This study employed an optical strain measurement method, called microdisplacements by machine vision photogrammetry (DISMAP), to measure both the global and local strain fields in microtensile specimens of cortical bone subjected to controlled uniaxial tension. The variation of local maximum principal strains was measured within the gauge region of samples as a function of applied tensile stress during testing. High gradients of local strain appeared around microstructural features in stressed bone even while the global strain for the entire gauge region showed a strong linear correlation with increasing tensile stress (r2 = 0.98, p < 0.0001). The highest local strain around micro-structural features in bone was 11.5-79.5 times higher than the global strain.

Measurement of microstructural strain in cortical bone

It is well known that mechanical factors affect bone remodeling such that increased mechanical demand results in net bone formation, whereas decreased demand results in net bone resorption. Current theories suggest that bone modeling and remodeling is controlled at the cellular level through signals mediated by osteocytes. The objective of this study was to investigate how macroscopically applied bone strains similar in magnitude to those that occur in vivo are manifest at the microscopic level in the bone matrix. Using a digital image correlation strain measurement technique, experimentally determined bone matrix strains around osteocyte lacuna resulting from macroscopic strains of approximately 2,000 microstrain (0.2%) reach levels of over 30,000 microstrain (3%) over fifteen times greater than the applied macroscopic strain. Strain patterns were highly heterogeneous and in some locations similar to observed microdamage around osteocyte lacuna indicating the resulting strains may represent the precursors to microdamage. This information may lead to a better understanding of how bone cells are affected by whole bone functional loading.

Osteocyte lacunae tissue strain in cortical bone

Current theories suggest that bone modeling and remodeling are controlled at the cellular level through signals mediated by osteocytes. However, the specific signals to which bone cells respond are still unknown. Two primary theories are: (1) osteocytes are stimulated via the mechanical deformation of the perilacunar bone matrix and (2) osteocytes are stimulated via fluid flow generated shear stresses acting on osteocyte cell processes within canaliculi. Recently, much focus has been placed on fluid flow theories since in vitro experiments have shown that bone cells are more responsive to analytically estimated levels of fluid shear stress than to direct mechanical stretching using macroscopic strain levels measured on bone in vivo. However, due to the complex microstructural organization of bone, local perilacunar bone tissue strains potentially acting on osteocytes cannot be reliably estimated from macroscopic bone strain measurements. Thus, the objective of this study was to quantify local perilacunar bone matrix strains due to macroscopically applied bone strains similar in magnitude to those that occur in vivo. Using a digital image correlation strain measurement technique, experimentally measured bone matrix strains around osteocyte lacunae resulting from macroscopic strains of approximately 2000 microstrain are significantly greater than macroscopic strain on average and can reach peak levels of over 30,000 microstrain locally. Average strain concentration factors ranged from 1.1 to 3.8, which is consistent with analytical and numerical estimates. This information should lead to a better understanding of how bone cells are affected by whole bone functional loading.

Microstructural strain near osteocyte lacuna in cortical bone in vitro

Mechanical factors affect bone remodeling such that increased mechanical demand results in net bone formation, whereas decreased demand results in net bone resorption. Two proposed mechanical signals are stress-generated fluid flow forces acting on cells and bone matrix deformation itself. A prominent current theory is that bone cells are more responsive to fluid flow than to mechanical strain. Recent experiments support this conclusion: bone cells increase their production of osteopontin (OPN) mRNA, prostaglandin (PGE(2)), and nitric oxide (NO) in response to fluid flow in contrast to cells stimulated by mechanical strain levels similar to those measured in vivo. However, when cells are subjected to substrate strains levels many times greater than those measured in vivo, increased biological activity again results. We assert that it is neither fluid flow nor matrix deformation per se, but rather the resulting cell deformation that causes cell biological response. Machined specimens of undamaged bovine cortical bone were subjected to increasing levels of macroscopic strain while observed under an optical microscope at 220X. Continuum level strain was measured using a standard foil strain gauge attached to the back of the specimen and ranged from 500 to 6,000 microstrain. Images of the specimen surface at each strain level were captured. To determine the level of osteocyte deformation that results from fluid flow in vitro, MLO-Y4 cells were cultured on collagen coated 190 cm2 plastic sheets and subjected to steady fluid flow at 16 dynes/cm(2). Images representing the initial undisturbed cell configuration and the configuration of the cells after ten minutes of fluid flow were acquired from a videotape of the flow experiment. The captured unloaded vs. loaded image pairs were analyzed to determine the local deformation and strain fields using a digital stereoimaging system. When subjected to a nominal continuum strain level approximately equal to that measured in humans in vivo during rigorous activity (2,000 microstrain), the local, osteocyte level strains can be as high as 12,000 to 15,000 microstrain (1.2% to 1.5%). Average osteocyte strains due to fluid flow in vitro increase from 7,972 microstrains after 16 seconds of flow to 22,856 microstrains after 64 seconds of flow. In contrast, maximum strains measured in vivo are approximately 1,800 microstrain in humans and up to 3,000 microstrain in other species. These data may help to explain why bone cells are more sensitive to fluid flow than substrate strain; fluid forces result in cell deformations much higher than those considered to be "physiological".

Primary adult human bone cells do not respond to tissue (continuum) level strains

Bone adapts to its mechanical environment, and, since the late 1800s, investigators have presumed that this adaptation relates to strain magnitude. Indeed, overwhelming evidence supports the view that either strain or some strain-related quantity stimulates bone adaptation or remodeling. Virtually all investigators, implicitly or explicitly, assume that the level of strain magnitude responsible for bone adaptation is that measured by strain gauges in vivo (i.e., 100-2500 microstrain) and that bone cells are directly deformed by strained matrix. We present evidence that bone cell deformation in this range does not cause bone adaptation. First, bone cells in vitro typically do not respond to average (continuum) levels of strain magnitude. Second, bone cells in vitro do respond to fluid flow-induced shear stresses in these ostensible physiological ranges. Third, in vivo strain magnitudes presumed to stimulate remodeling reflect only averages, and not local peaks, which are 2-15 times higher. Thus, we hypothesize that sensing cells do not respond to levels of strain presumed to be physiological.

Machine vision photogrammetry: a technique for measurement of microstructural strain in cortical bone

Understanding local microstructural deformations and strains in cortical bone may lead to a better understanding of cortical bone damage development, fracture, and remodeling. Traditional experimental techniques for measuring deformation and strain do not allow characterization of these quantities at the microstructural level in cortical bone. This study describes a technique based on digital stereoimaging used to measure the microstructural strain fields in cortical bone. The technique allows the measurement of material surface displacements and strains by comparing images acquired from a specimen at two distinct stress states. The accuracy of the system is investigated by analyzing an undeformed image set; the test image is identical to the reference image but translated by a known pixel amount. An increase in the correlation sub-image train parameter results in an increase in displacement measurement accuracy from 0.049 to 0.012 pixels. Errors in strain calculated from the measured displacement field were between 39 and 564 microstrain depending upon the sub-image train size and applied image displacement. The presence of a microcrack in cortical bone results in local strain at the crack tip reaching 0.030 (30,000 microstrain) and 0.010 (10,000 microstrain) near osteocyte lacunae. It is expected that the use of this technique will allow a greater understanding of bone strength and fracture as well as bone mechanotransduction.