Limitations of the continuum assumption in cancellous bone

Trabecular bone variation in the gorilla calcaneus

OBJECTIVES

Calcaneal external shape differs among nonhuman primates relative to locomotion. Such relationships between whole-bone calcaneal trabecular structure and locomotion, however, have yet to be studied. Here we analyze calcaneal trabecular architecture in Gorilla gorilla gorilla, Gorilla beringei beringei, and G. b. graueri to investigate general trends and fine-grained differences among gorilla taxa relative to locomotion.

MATERIALS AND METHODS

Calcanei were micro-CT scanned. A three-dimensional geometric morphometric sliding semilandmark analysis was carried out and the final landmark configurations used to position 156 volumes of interest. Trabecular thickness (Tb.Th), trabecular spacing (Tb.Sp), and bone volume fraction (BV/TV) were calculated using the BoneJ plugin for ImageJ and MATLAB. Non-parametric MANOVAs were run to test for significant differences among taxa in parameter raw values and z-scores. Parameter distributions were visualized using color maps and summarized using principal components analysis.

RESULTS

There are no significant differences in raw BV/TV or Tb.Th among gorillas, however G. b. beringei significantly differs in z-scores for both parameters (p = <0.0271). All three taxa exhibit relatively lower BV/TV and Tb.Th in the posterior half of the calcaneus. This gradation is exacerbated in G. b. beringei. G. b. graueri significantly differs from other taxa in Tb.Sp z-scores (p < 0.001) indicating a different spacing distribution.

DISCUSSION

Relatively higher Tb.Th and BV/TV in the anterior calcaneus among gorillas likely reflects higher forces associated with body mass (transmitted through the subtalar joint) relative to forces transferred through the posterior calcaneus. The different Tb.Sp pattern in G. b. graueri may reflect proposed differences in foot positioning during locomotion.

Correlation between tibial and femoral bone and cartilage changes in end-stage knee osteoarthritis

Knee osteoarthritis is a whole joint disease highlighting the coupling of cartilage and bone adaptations. However, the structural properties of the subchondral bone plate (SBP) and underlying subchondral trabecular bone (STB) in the femoral compartment have received less attention compared to the tibial side. Furthermore, how the properties in the femoral compartment relate to those in the corresponding tibial site is unknown. Therefore, this study aimed to quantify the structural bone and cartilage morphology in the femoral compartment and investigate its association with those of the tibial plateau. Specifically, tibial plateaus and femoral condyles were retrieved from 28 patients with end-stage knee-osteoarthritis (OA) and varus deformity. The medial condyle of tibial plateaus and the distal part of the medial femoral condyles were micro-CT scanned (20.1 μm/voxel). Cartilage thickness (Cart.Th), SBP, and STB microarchitecture were quantified. Significant (P < <.001; 0.79 ≤ r ≤ 0.97) correlations with a relative difference within 10% were found between the medial side of the femoral and tibial compartments. The highest correlations were found for SBP porosity (r = 0.97, mean absolute difference of 0.50%, and mean relative difference of 9.41%) and Cart.Th (r = 0.96, mean absolute difference of 0.18 mm, and relative difference of 7.08%). The lowest correlation was found for trabecular thickness (r = 0.79, mean absolute difference of 21.07 μm, and mean relative difference of 5.17%) and trabecular number (r = 0.79, mean absolute difference of 0.18 mm-1, and relative difference of 5.02%). These findings suggest that the distal femur is affected by OA in a similar way as the proximal tibia. Given that bone adaptation is a response to local mechanical forces, our results suggest that varus deformity similarly affects the stress distribution of the medial tibial plateau and the medial distal femur.

Investigation on primary stability of dental implants through considering peri-implant bone damage, caused by small and large deformations: A validated non-linear micro finite element study

Primary stability of a dental implant is defined as its ability to resist the applied load without showing excessive damage in peri-implant bone, which is a prerequisite for secondary stability, and consequently for implantation success. The main goal of this study was to develop a validated micro-finite element (μFE) approach to assess the primary stability of dental implants in terms of stiffness, stiffness reduction, and irreversible displacement of the bone-implant system, subjected to an increasing step-wise quasi-static compressive loading-unloading test. The μFE models were generated based on the μCT images of bone, taken from extracted bovine tibia trabecular bone samples after drilling and implantation. A tissue constitutive model was considered for trabecular bone by describing elasto-plasticity with a modified von Mises yield criterion and element deletion technique to account for trabecular bone damage behavior. Then, the obtained force-displacement curves from the simulation were compared with the in-vitro mechanical test curves to evaluate the validity of the model. The results showed that the proposed μFE model could be properly predict the bone-implant system mechanical response in terms of irreversible displacement (R2 = 0.99), stiffness (R2 = 0.77), and stiffness reduction (R2 = 0.72) of the bone-implant construct for all the applied displacements without a significant difference from the unit slope and zero intercept of the QQ-plot (p-value<0.05). Moreover, a qualitative agreement was seen between the peri-implant bone damage predicted by the μFE model and the observed from μCT images. The adopted methodology used in this study can predict the mechanical failure response of the bone-implant system, which can be employed as a representative tool to study the effects of various dental implant design parameters on the primary stability with the ultimate goal of optimizing dental implants design.

Approach to standardized material characterization of the human lumbopelvic system-Specification, preparation and storage

The complexity of the osseo-ligamentous lumbopelvic system has made it difficult to perform both, the overall preparation as well as specimen harvesting and material testing with a reasonable amount of time and personnel. The logistics of such studies present a hurdle for reproducibility. A structured procedure was developed and proved, which allows all necessary steps to be carried out reproducibly and in a reasonable time. This enables the extraction of 26 soft tissue, 33 trabecular and 32 cortical bone specimens from this anatomical region per cadaver. The integrity of the specimens remains maintained while keeping requirements within manageable limits. The practicability of the intended five-day specimen harvesting and testing procedure could be demonstrated on five test and two pre-test sequences. The intended minimization of physical, biological, and chemical external influences on specimens could be achieved. All protocols, instructions and models of preparation and storage devices are included in the supporting information. The high grade of applicability and reproducibility will lead to better comparability between different biomechanical investigations. This procedure proven on the human pelvis is transferable to other anatomical regions.

Mechanical parameter assessment of fresh human cancellous bone of the femoral head in atraumatic femoral head necrosis and primary coxarthrosis

BACKGROUND

Atraumatic femoral head necrosis is a rare pathological change of the femoral head. It is characterized by local necrosis of the cancellous bone as a result of reduced blood supply to the bone. Even today it remains unclear how to assess the hardness of the necrosis, whether it is soft tissue that is easily removed, or hard tissue that is difficult to resect.

METHODS

Femoral heads with primary coxarthrosis were selected as a comparison group. For this purpose, 49 femoral heads obtained during total hip arthroplasty surgery with either condition (23 femoral head necrosis, 26 coxarthrosis) were transferred to the testing laboratory in fresh condition. Cylindrical specimens were obtained using a tenon cutter along the main trabecular load direction in the subchondral region of the femoral head. Additionally, thin bone slices were extracted proximal and distal to the specimens for density measurements. Brass plates were glued to the circular surfaces of the specimens. After curing of the adhesive, the specimens were mounted in the testing machine and destructive uniaxial compression tests were conducted.

FINDINGS

The recorded mean compressive strengths and elastic moduli were almost identical for both groups, but the necrosis group showed significantly higher data scattering and range regarding the elastic modulus. The mean density of the coxarthrosis specimens was significantly higher than that of the necrotic specimens.

INTERPRETATION

The mechanical properties of cancellous bone vary considerably in the presence of femoral head necrosis. The existence of hard necrosis implies a potential challenge regarding the clinical resection of these tissues.

Effects of loading rate, age, and morphology on the material properties of human rib trabecular bone

The material and morphometric properties of trabecular bone have been studied extensively in bones bearing significant weight, such as the appendicular long bones and spine. Less attention has been devoted to the ribs, where quantification of material properties is vital to understanding thoracic injury. The objective of this study was to quantify the compressive material properties of human rib trabecular bone and assess the effects of loading rate, age, and morphology on the material properties. Material properties were quantified via uniaxial compression tests performed on trabecular bone samples at two loading rates: 0.005 s^-1 and 0.5 s^-1. Morphometric parameters of each sample were quantified before testing using micro-computed tomography. Rib trabecular bone material properties were lower on average compared to trabecular bone from other anatomical locations. Morphometric parameters indicated an anisotropic structure with low connectivity and a sparser density of trabeculae in the rib compared to other locations. No significant differences in material properties were observed between the tested loading rates. Material properties were only significantly correlated with age at the 0.005 s^-1 loading rate, and no morphometric parameter was significantly correlated with age. Trabecular separation and thickness were most strongly correlated with the material properties, indicating the sparser trabecular matrix likely contributed to the lower material property values compared to other sites. The novel trabecular bone material properties reported in this study can be used to improve the thoracic response and injury prediction of computational models.

Comparison of bone formation mediated by bone morphogenetic protein delivered by nanoclay gels with clinical techniques (autograft and InductOs®) in an ovine bone model

Development of a growth factor delivery vehicle providing appropriate temporal-spatial release together with an appropriate preclinical large animal model to evaluate bone formation is critical in the development of delivery strategies for bone tissue regeneration. Smectite nanoclays such as LAPONITE™ possess unique thixotropic and protein retention properties offering promise for use in growth factor delivery in bone repair and regeneration. This study has examined bone formation mediated by a clinically approved growth factor delivery system (InductOs®) in combination with Laponite gel in an aged female ovine femoral condyle defect preclinical model (10 weeks). Two different designs, one containing a low volume of Laponite gel (LLG) in combination with the InductOs® absorbable collagen sponge (ACS), the other in which Laponite gel formed the implant (HLG), were compared against InductOs® alone and an autograft positive control. Thus, five groups: (i) empty defect, (ii) autograft, (iii) BMP2 + ACS, (iv) BMP2 + ACS + LLG and (v) BMP2 + HLG + ACS were examined in 9 mm × 12 mm defects performed bilaterally in the medial femoral condyles of 24 aged (>5 years) sheep. Bone formation within the defect was assessed using micro-computed tomography (micro-CT), digital volume correlation (DVC) for biomechanical characterisation as well as histology. The autograft and InductOs® mediated enhanced bone formation (p < 0001) compared to blank controls, while no significant differences were observed between the Laponite/Collagen/BMP delivery vehicles. However, the current study illustrated the excellent biocompatibility of Laponite and its ability to deliver localised active BMP-2, with the opportunity for improved efficacy with further optimisation. Interestingly, DVC-computed strain distributions indicated that the regenerated bone structure is mechanically adapted to bear external loads from the early remodelling stages of the bone reparation cascade. The current studies of selected nanoclay delivery platforms for BMP, assessed in a clinically relevant large animal model auger well for the development of bone fracture therapeutics for an ageing population.

A shrewd inspection of vertebral regionalization in large shrews (Soricidae: Crocidurinae)

The regionalization of the mammalian spinal column is an important evolutionary, developmental, and functional hallmark of the clade. Vertebral column regions are usually defined using transitions in external bone morphology, such as the presence of transverse foraminae or rib facets, or measurements of vertebral shape. Yet the internal structure of vertebrae, specifically the trabecular (spongy) bone, plays an important role in vertebral function, and is subject to the same variety of selective, functional, and developmental influences as external bone morphology. Here we investigated regionalization of external and trabecular bone morphology in the vertebral column of a group of shrews (family Soricidae). The primary goals of this study were to: 1) determine if vertebral trabecular bone morphology is regionalized in large shrews, and if so, in what configuration relative to external morphology; 2) assess correlations between trabecular bone regionalization and functional or developmental influences; and 3) determine if external and trabecular bone regionalization patterns provide clues about the function of the highly modified spinal column of the hero shrew Scutisorex.

Trabecular bone is regionalized along the soricid vertebral column, but the configuration of trabecular bone regions does not match that of the external vertebral morphology, and is less consistent across individuals and species. The cervical region has the most distinct and consistent trabecular bone morphology, with dense trabeculae indicative of the ability to withstand forces in a variety of directions. Scutisorex exhibits an additional external morphology region compared to unmodified shrews, but this region does not correspond to a change in trabecular architecture.

Although trabecular bone architecture is regionalized along the soricid vertebral column, and this regionalization is potentially related to bone functional adaptation, there are likely aspects of vertebral functional regionalization that are not detectable using trabecular bone morphology. For example, the external morphology of the Scutisorex lumbar spine shows signs of an extra functional region that is not apparent in trabecular bone analyses. It is possible that body size and locomotor mode affect the degree to which function is manifest in trabecular bone, and broader study across mammalian size and ecology is warranted to understand the relationship between trabecular bone morphology and other measures of vertebral function such as intervertebral range of motion.

Assessment of the Bone Mineral Density and Microstructure of the Human Femoral Head according to Different Tip-apex Distances Can Guide the Treatment of Intertrochanteric Hip Fractures

Purpose

We analyzed the microstructure and bone mineral density (BMD) of the trabecular bone in the femoral head of patients with osteoporosis.

Materials and Methods

Sixteen femoral heads with osteoporotic femoral neck fractures underwent micro-computed tomography scanning. In each tip-apex distance (TAD) of 15, 20, and 25 mm, five regions of interest (ROIs) were extracted from the central, anterior, posterior, superior, and inferior sections. A total of 15 ROIs were extracted from TADs of 15, 20, and 25 mm. The measurement parameters included BMD, percent bone volume: bone volume/total volume (BV/TV), trabecular thickness (Tb.Th), trabecular number (Tb.N), structural model index (SMI), and degree of anisotropy (DOA).

Results

The lowest BMD and BV/TV values were observed in the inferior region and differed significantly from those in other regions (P<0.05). Lower Tb.Th and Tb.N values were observed in the inferior region compared with those in the central region (P<0.05). The highest SMI value was observed in the inferior region (P<0.05). With TAD of 15 and 20 mm, the DOA values in the inferior region were lower than those in the anterior region (P<0.05). Lower BMD and BV/TV values were observed in the anterior, central, and inferior regions of TAD of 15 mm compared with those in the corresponding regions of TAD of 25 mm (P<0.05).

Conclusion

Positioning the lag screw between TAD of 20 to 25 mm and in the inferior region is recommended, and TAD of less than 15 mm is not recommended.

The effect of Staphylococcus aureus exposure on white-tailed deer trabecular bone stiffness and yield

With a growing number of osteomyelitis diagnoses, many of which are linked to Staphylococcus aureus (S. aureus), it is imperative to understand the pathology of S. aureus in relation to bone to provide better diagnostics and patient care. While the cellular mechanisms of S. aureus and osteomyelitis have been studied, little information exists on the biomechanical effects of such infections. The aim of this study was to determine the effect of S. aureus exposure on the stiffness and yield of trabecular bone tissue. S. aureus-ATCC-12600, a confirmed biofilm producer, along with one hundred and three trabecular cubes (5 × 5 × 5 mm) from the proximal tibiae of Odocoileus virginianus (white-tailed deer) were used in this experiment. Bone cubes were disinfected and then swabbed to confirm no residual living microbes or endospore contamination before inoculation with S. aureus (test group) or sterile nutrient broth (control group) for 72 h. All cubes were then tested in compression until yield using an Instron 5942 Single-Column machine. Structural stiffness (N/mm) and yield (MPa) were calculated and compared between the two groups. Our results revealed that acute exposure to S. aureus, within the context of our deer tibia model, does not significantly decrease trabecular bone stiffness or yield. The results of this study may be of value clinically when assessing fracture risks for osteomyelitis or other patients whose cultures test positive for S. aureus.

Characterization of Ultralow Density Cellular Solids: Lessons from 30 years of Bone Biomechanics Research

Advances in additive manufacturing techniques have enabled the development of micro-architectured materials displaying a combination of low-density and lightweight structures with high specific strength and toughness. The mechanical performance of micro-architectured materials can be assessed using standard techniques; however, when studying low- and ultralow density micro-architectured materials, standard characterization techniques can be subject to experimental artifacts. Additionally, quantitative assessment and comparisons of microarchitectures with distinct lattice patterns is not always straightforward. Cancellous bone is a natural, ultralow density (porosity often exceeding 90%), irregular, cellular solid that has been thoroughly characterized in terms of micro-architecture and mechanical performance over the past 30 years. However, most the literature on cancellous bone mechanical properties and micro-structure-function relationships is in the medical literature and is not immediately accessible to materials designers. Here we provide a brief review of state-of-the-art approaches for characterizing the micro-architecture and mechanical performance of ultralow density cancellous bone, including methods of addressing experimental artifacts during mechanical characterization of ultralow density cellular solids, methods of quantifying microarchitecture, and currently understood structure-function relationships.

Influence of osteoporosis on the compressive properties of femoral cancellous bone and its dependence on various density parameters

Data collection of mechanical parameters from compressive tests play a fundamental role in FE modelling of bone tissues or the developing and designing of bone implants, especially referring to osteoporosis or other forms of bone loss. A total of 43 cylindrical samples (Ø8 × 16 mm) were taken from 43 freshly frozen proximal femora using a tenon cutter. All femora underwent BMD measurement and additionally apparent- and relative- and bulk density (ρ_app, ρ_r, ρ_b) were determined using samples bordering the compressive specimen on the proximal and distal regions. All samples were classified as "normal", "osteopenia" and "osteoporosis" based on the DEXA measurements. Distal apparent density was most suitable for predicting bone strength and BMD. One novel aspect is the examination of the plateau stress as it describes the stress at which the failure of spongious bone progresses. No significant differences in mechanical properties (compressive modulus E; compressive stress σ_max and plateau stress σ_p) were found between osteopenic and osteoporotic bone. The results suggest that already in the case of a known osteopenia, actions should be taken as they are applied in the case of osteoporosis A review of the literature regarding extraction and testing methods illustrates the urgent need for standardized biomechanical compressive material testing.

Fluid-structure interaction (FSI) modeling of bone marrow through trabecular bone structure under compression

The present study has sought to investigate the fluid characteristic and mechanical properties of trabecular bone using fluid-structure interaction (FSI) approach under different trabecular bone orientations. This method imposed on trabecular bone structure at both longitudinal and transverse orientations to identify effects on shear stress, permeability, stiffness and stress regarded to the trabeculae. Sixteen FSI models were performed on different range trabecular cubes of 27 mm³ with eight models developed for each longitudinal and transverse direction. Results show that there was a moderate correlation between permeability and porosity, and surface area in the longitudinal and transverse orientations. For the longitudinal orientation, the permeability values varied between 3.66 × 10^-8 and 1.9 × 10^-7 and the sheer stress values varied between 0.05 and 1.8 Pa, whilst for the transverse orientation, the permeability values varied between 5.95 × 10^-10 and 1.78 × 10^-8 and the shear stress values varied between 0.04 and 3.1 Pa. Here, transverse orientation limits the fluid flow from passing through the trabeculae due to high shear stress disturbance generated within the trabecular bone region. Compared to physiological loading direction (longitudinal orientation), permeability is higher within the range known to trigger a response in bone cells. Additionally, shear stresses also increase with bone surface area. This study suggests the shear stress within bone marrow in real trabecular architecture could provide the mechanical signal to marrow cells that leads to bone anabolism and can depend on trabecular orientation.

Parametric investigation of the effects of load level on fatigue crack growth in trabecular bone based on artificial neural network computation

This study reports the development of an artificial neural network computation model to predict the accumulation of crack density and crack length in cancellous bone under a cyclic load. The model was then applied to conduct a parametric investigation into the effects of load level on fatigue crack accumulation in cancellous bone. The method was built in three steps: (1) conducting finite element simulations to predict fatigue growth of different three-dimensional micro-computed tomography cancellous bone specimens considering input combinations based on a factorial experimental design; (2) performing a training stage of an artificial neural network based on the results of step 1; and (3) applying the trained artificial neural network to rapidly predict the crack density and the crack length growth for cancellous bone under a cyclic loading for a given applied apparent strain, cycle frequency, bone volume fraction, bone density and apparent elastic modulus.

Deciphering an extreme morphology: bone microarchitecture of the hero shrew backbone (Soricidae: Scutisorex)

Biological structures with extreme morphologies are puzzling because they often lack obvious functions and stymie comparisons to homologous or analogous features with more typical shapes. An example of such an extreme morphotype is the uniquely modified vertebral column of the hero shrew Scutisorex, which features numerous accessory intervertebral articulations and massively expanded transverse processes. The function of these vertebral structures is unknown, and it is difficult to meaningfully compare them to vertebrae from animals with known behavioural patterns and spinal adaptations. Here, we use trabecular bone architecture of vertebral centra and quantitative external vertebral morphology to elucidate the forces that may act on the spine of Scutisorex and that of another large shrew with unmodified vertebrae (Crocidura goliath). X-ray micro-computed tomography (µCT) scans of thoracolumbar columns show that Scutisorex thori is structurally intermediate between C. goliath and S. somereni internally and externally, and both Scutisorex species exhibit trabecular bone characteristics indicative of higher in vivo axial compressive loads than C. goliath. Under compressive load, Scutisorex vertebral morphology is adapted to largely restrict bending to the sagittal plane (flexion). Although these findings do not solve the mystery of how Scutisorex uses its byzantine spine in vivo, our work suggests potentially fruitful new avenues of investigation for learning more about the function of this perplexing structure.

Baby steps towards linking calcaneal trabecular bone ontogeny and the development of bipedal human gait

Trabecular bone structure in adulthood is a product of a process of modelling during ontogeny and remodelling throughout life. Insight into ontogeny is essential to understand the functional significance of trabecular bone structural variation observed in adults. The complex shape and loading of the human calcaneus provides a natural experiment to test the relationship between trabecular morphology and locomotor development. We investigated the relationship between calcaneal trabecular bone structure and predicted changes in loading related to development of gait and body size in growing children. We sampled three main trabecular regions of the calcanei using micro-computed tomography scans of 35 individuals aged between neonate to adult from the Norris Farms #36 site (1300 AD, USA) and from Cambridge (1200-1500 AD, UK). Trabecular properties were calculated in volumes of interest placed beneath the calcaneocuboid joint, plantar ligaments, and posterior talar facet. At birth, thin trabecular struts are arranged in a dense and relatively isotropic structure. Bone volume fraction strongly decreases in the first year of life, whereas anisotropy and mean trabecular thickness increase. Dorsal compressive trabecular bands appear around the onset of bipedal walking, although plantar tensile bands develop prior to predicted propulsive toe-off. Bone volume fraction and anisotropy increase until the age of 8, when gait has largely matured. Connectivity density gradually reduces, whereas trabeculae gradually thicken from birth until adulthood. This study demonstrates that three different regions of the calcaneus develop into distinct adult morphologies through varying developmental trajectories. These results are similar to previous reports of ontogeny in human long bones and are suggestive of a relationship between the mechanical environment and trabecular bone architecture in the human calcaneus during growth. However, controlled experiments combined with more detailed biomechanical models of gait maturation are necessary to establish skeletal markers linking growth to loading. This has the potential to be a novel source of information for understanding loading levels, activity patterns, and perhaps life history in the fossil record.

Global and site-specific analysis of bone in a rat model of spinal cord injury-induced osteoporosis

Micro-Computed Tomography bone analysis is the gold standard method for assessing trabecular and cortical bone microarchitecture in small animal bones. This technique reports morphometric parameters as averages over selected volumes of interest (VOIs). This study proposes the introduction of an additional global 2D morphometric step into the analysis process, that provides a survey of the underlying morphometric variation present throughout both trabecular and cortical bone. The visualisation of these morphometric distributions provides a systematic approach to VOI selection that provides rationale and adds confidence to subsequent 3D morphometric analysis. To test the applicability and value of this methodological addition it was applied to the distal femur of a rat model of spinal cord injury (SCI)-induced osteoporosis. The 2D morphometric variation of both trabecular and cortical bone was quantified as a function of bone length. SCI-induced osteoporosis was localised in i) trabecular bone, where metaphyseal bone was more severely affected than epiphyseal bone, and there was a significant reduction in Distal Femoral Trabecular Extent, a new parameter defined here that quantifies how far trabecular bone penetrates in to the marrow cavity, ii) cortical bone, where diaphyseal bone underwent significant lowering of both cortical area and thickness, while distal-metaphyseal bone did not. Theses site-specific changes were validated, further elucidated and compared with follow-up conventional 3D analysis. The techniques applied here are equally applicable to other long bones (tibia, humerus, radius, ulna), other types of imaging modality and other types of experimental design including the effects of rehabilitation, aging, loading, gene knockout and pharmacological intervention.

•

2D morphological surveying identifies regions warranting further 3D investigation.

•

Trabecular microarchitecture site-specifically varies in the distal femur.

•

SCI-induced osteoporosis changes metaphyseal more than epiphyseal trabecular bone.

•

SCI-induced osteoporosis reduced the extent of metaphyseal trabecular bone.

Analysis of Trabecular Bone Mechanics Using Machine Learning

"Bone remodeling" is a dynamic process, and mutliphase analysis incorporated with the forecasting algorithm can help the biologists and orthopedics to interpret the laboratory generated results and to apply them in improving applications in the fields of "drug design, treatment, and therapy" of diseased bones. The metastasized bone microenvironment has always remained a challenging puzzle for the researchers. A multiphase computational model is interfaced with the artificial intelligence algorithm in a hybrid manner during this research. Trabecular surface remodeling is presented in this article, with the aid of video graphic footage, and the associated parametric thresholds are derived from artificial intelligence and clinical data.

Cancellous bone and theropod dinosaur locomotion. Part I-an examination of cancellous bone architecture in the hindlimb bones of theropods

This paper is the first of a three-part series that investigates the architecture of cancellous ('spongy') bone in the main hindlimb bones of theropod dinosaurs, and uses cancellous bone architectural patterns to infer locomotor biomechanics in extinct non-avian species. Cancellous bone is widely known to be highly sensitive to its mechanical environment, and has previously been used to infer locomotor biomechanics in extinct tetrapod vertebrates, especially primates. Despite great promise, cancellous bone architecture has remained little utilized for investigating locomotion in many other extinct vertebrate groups, such as dinosaurs. Documentation and quantification of architectural patterns across a whole bone, and across multiple bones, can provide much information on cancellous bone architectural patterns and variation across species. Additionally, this also lends itself to analysis of the musculoskeletal biomechanical factors involved in a direct, mechanistic fashion. On this premise, computed tomographic and image analysis techniques were used to describe and analyse the three-dimensional architecture of cancellous bone in the main hindlimb bones of theropod dinosaurs for the first time. A comprehensive survey across many extant and extinct species is produced, identifying several patterns of similarity and contrast between groups. For instance, more stemward non-avian theropods (e.g. ceratosaurs and tyrannosaurids) exhibit cancellous bone architectures more comparable to that present in humans, whereas species more closely related to birds (e.g. paravians) exhibit architectural patterns bearing greater similarity to those of extant birds. Many of the observed patterns may be linked to particular aspects of locomotor biomechanics, such as the degree of hip or knee flexion during stance and gait. A further important observation is the abundance of markedly oblique trabeculae in the diaphyses of the femur and tibia of birds, which in large species produces spiralling patterns along the endosteal surface. Not only do these observations provide new insight into theropod anatomy and behaviour, they also provide the foundation for mechanistic testing of locomotor hypotheses via musculoskeletal biomechanical modelling.

Trabecular Bone Structural Variation in the Proximal Sacrum Among Primates

The sacrum occupies a functionally important anatomical position as part of the pelvic girdle and vertebral column. Sacral orientation and external morphology in modern humans are distinct from those in other primates and compatible with the demands of habitual bipedal locomotion. Among nonhuman primates, however, how sacral anatomy relates to positional behaviors is less clear. As an alternative to evaluation of the sacrum's external morphology, this study assesses if the sacrum's internal morphology (i.e., trabecular bone) differs among extant primates. The primary hypothesis tested is that trabecular bone parameters with established functional relevance will differ in the first sacral vertebra (S1) among extant primates that vary in positional behaviors. Results for analyses of individual variables demonstrate that bone volume fraction, degree of anisotropy, trabecular number, and size-corrected trabecular thickness differ among primates grouped by positional behaviors to some extent, but not always in ways consistent with functional expectations. When examined as a suite, these trabecular parameters distinguish obligate bipeds from other positional behavior groups; and, the latter three trabecular bone variables further distinguish knuckle-walking terrestrial quadrupeds from manual suspensor-brachiators, vertical clingers and leapers, and arboreal quadrupeds, as well as between arboreal and terrestrial quadrupeds. As in other regions of the skeleton in modern humans, trabecular bone in S1 exhibits distinctively low bone volume fraction. Results from this study of extant primate S1 trabecular bone structural variation provide a functional context for interpretations concerning the positional behaviors of extinct primates based on internal sacral morphology. Anat Rec, 302:1354-1371, 2019. © 2018 Wiley Periodicals, Inc.

Comparison of patient-specific computational models vs. clinical follow-up, for adjacent segment disc degeneration and bone remodelling after spinal fusion

Spinal fusion is a standard surgical treatment for patients suffering from low back pain attributed to disc degeneration. However, results are somewhat variable and unpredictable. With fusion the kinematic behaviour of the spine is altered. Fusion and/or stabilizing implants carrying considerable load and prevent rotation of the fused segments. Associated with these changes, a risk for accelerated disc degeneration at the adjacent levels to fusion has been demonstrated. However, there is yet no method to predict the effect of fusion surgery on the adjacent tissue levels, i.e. bone and disc. The aim of this study was to develop a coupled and patient-specific mechanoregulated model to predict disc generation and changes in bone density after spinal fusion and to validate the results relative to patient follow-up data. To do so, a multiscale disc mechanoregulation adaptation framework was developed and coupled with a previously developed bone remodelling algorithm. This made it possible to determine extra cellular matrix changes in the intervertebral disc and bone density changes simultaneously based on changes in loading due to fusion surgery. It was shown that for 10 cases the predicted change in bone density and degeneration grade conforms reasonable well to clinical follow-up data. This approach helps us to understand the effect of surgical intervention on the adjacent tissue remodelling. Thereby, providing the first insight for a spine surgeon as to which patient could potentially be treated successfully by spinal fusion and in which patient has a high risk for adjacent tissue changes.

Mapping anisotropy improves QCT-based finite element estimation of hip strength in pooled stance and side-fall load configurations

Hip fractures are one of the most severe consequences of osteoporosis. Compared to the clinical standard of DXA-based aBMD at the femoral neck, QCT-based FEA delivers a better surrogate of femoral strength and gains acceptance for the calculation of hip fracture risk when a CT reconstruction is available. Isotropic, homogenised voxel-based, finite element (hvFE) models are widely used to estimate femoral strength in cross-sectional and longitudinal clinical studies. However, fabric anisotropy is a classical feature of the architecture of the proximal femur and the second determinant of the homogenised mechanical properties of trabecular bone. Due to the limited resolution, fabric anisotropy cannot be derived from clinical CT reconstructions. Alternatively, fabric anisotropy can be extracted from HR-pQCT images of cadaveric femora. In this study, fabric anisotropy from HR-pQCT images was mapped onto QCT-based hvFE models of 71 human proximal femora for which both HR-pQCT and QCT images were available. Stiffness and ultimate load computed from anisotropic hvFE models were compared with previous biomechanical tests in both stance and side-fall configurations. The influence of using the femur-specific versus a mean fabric distribution on the hvFE predictions was assessed. Femur-specific and mean fabric enhance the prediction of experimental ultimate force for the pooled, i.e. stance and side-fall, (isotropic: r²=0.81, femur-specific fabric: r²=0.88, mean fabric: r²=0.86,p<0.001) but not for the individual configurations. Fabric anisotropy significantly improves bone strength prediction for the pooled configurations, and mapped fabric provides a comparable prediction to true fabric. The mapping of fabric anisotropy is therefore expected to help generate more accurate QCT-based hvFE models of the proximal femur for personalised or multiple load configurations.

Relationships between in vivo dynamic knee joint loading, static alignment and tibial subchondral bone microarchitecture in end-stage knee osteoarthritis

OBJECTIVE

To study, in end-stage knee osteoarthritis (OA) patients, relationships between indices of in vivo dynamic knee joint loads obtained pre-operatively using gait analysis, static knee alignment, and the subchondral trabecular bone (STB) microarchitecture of their excised tibial plateau quantified with 3D micro-CT.

DESIGN

Twenty-five knee OA patients scheduled for total knee arthroplasty underwent pre-operative gait analysis. Mechanical axis deviation (MAD) was determined radiographically. Following surgery, excised tibial plateaus were micro-CT-scanned and STB microarchitecture analysed in four subregions (anteromedial, posteromedial, anterolateral, posterolateral). Regional differences in STB microarchitecture and relationships between joint loading and microarchitecture were examined.

RESULTS

STB microarchitecture differed among subregions (P < 0.001), anteromedially exhibiting highest bone volume fraction (BV/TV) and lowest structure model index (SMI). Anteromedial BV/TV and SMI correlated strongest with the peak external rotation moment (ERM; r = -0.74, r = 0.67, P < 0.01), despite ERM being the lowest (by factor of 10) of the moments considered, with majority of ERM measures below accuracy thresholds; medial-to-lateral BV/TV ratios correlated with ERM, MAD, knee adduction moment (KAM) and internal rotation moment (|r|-range: 0.54-0.74). When controlling for walking speed, KAM and MAD, the ERM explained additional 11-30% of the variations in anteromedial BV/TV and medial-to-lateral BV/TV ratio (R² = 0.59, R² = 0.69, P < 0.01).

CONCLUSIONS

This preliminary study suggests significant associations between tibial plateau STB microarchitecture and knee joint loading indices in end-stage knee OA patients. Particularly, anteromedial BV/TV correlates strongest with ERM, whereas medial-to-lateral BV/TV ratio correlates strongest with indicators of medial-to-lateral joint loading (MAD, KAM) and rotational moments. However, associations with ERM should be interpreted with caution.

Trabecular architecture in the sciuromorph femoral head: allometry and functional adaptation

BACKGROUND

Sciuromorpha (squirrels and close relatives) are diverse in terms of body size and locomotor behavior. Individual species are specialized to perform climbing, gliding or digging behavior, the latter being the result of multiple independent evolutionary acquisitions. Each lifestyle involves characteristic loading patterns acting on the bones of sciuromorphs. Trabecular bone, as part of the bone inner structure, adapts to such loading patterns. This network of thin bony struts is subject to bone modeling, and therefore reflects habitual loading throughout lifetime. The present study investigates the effect of body size and lifestyle on trabecular structure in Sciuromorpha.

METHODS

Based upon high-resolution computed tomography scans, the femoral head 3D inner microstructure of 69 sciuromorph species was analyzed. Species were assigned to one of the following lifestyle categories: arboreal, aerial, fossorial and semifossorial. A cubic volume of interest was selected in the center of each femoral head and analyzed by extraction of various parameters that characterize trabecular architecture (degree of anisotropy, bone volume fraction, connectivity density, trabecular thickness, trabecular separation, bone surface density and main trabecular orientation). Our analysis included evaluation of the allometric signals and lifestyle-related adaptation in the trabecular parameters.

RESULTS

We show that bone surface density, bone volume fraction, and connectivity density are subject to positive allometry, and degree of anisotropy, trabecular thickness, and trabecular separation to negative allometry. The parameters connectivity density, bone surface density, trabecular thickness, and trabecular separation show functional signals which are related to locomotor behavior. Aerial species are distinguished from fossorial ones by a higher trabecular thickness, lower connectivity density and lower bone surface density. Arboreal species are distinguished from semifossorial ones by a higher trabecular separation.

CONCLUSION

This study on sciuromorph trabeculae supplements the few non-primate studies on lifestyle-related functional adaptation of trabecular bone. We show that the architecture of the femoral head trabeculae in Sciuromorpha correlates with body mass and locomotor habits. Our findings provide a new basis for experimental research focused on functional significance of bone inner microstructure.

A novel in silico method to quantify primary stability of screws in trabecular bone

Insufficient primary stability of screws in bone leads to screw loosening and failure. Unlike conventional continuum finite-element models, micro-CT based finite-element analysis (micro-FE) is capable of capturing the patient-specific bone micro-architecture, providing accurate estimates of bone stiffness. However, such in silico models for screws in bone highly overestimate the apparent stiffness. We hypothesized that a more accurate prediction of primary implant stability of screws in bone is possible by considering insertion-related bone damage. We assessed two different screw types and loading scenarios in 20 trabecular bone specimens extracted from 12 cadaveric human femoral heads (N = 5 for each case). In the micro-FE model, we predicted specimen-specific Young's moduli of the peri-implant bone damage region based on morphometric parameters such that the apparent stiffness of each in silico model matched the experimentally measured stiffness of the corresponding in vitro specimen as closely as possible. The standard micro-FE models assuming perfectly intact peri-implant bone overestimated the stiffness by over 330%. The consideration of insertion related damaged peri-implant bone corrected the mean absolute percentage error down to 11.4% for both loading scenarios and screw types. Cross-validation revealed a mean absolute percentage error of 14.2%. We present the validation of a novel micro-FE modeling technique to quantify the apparent stiffness of screws in trabecular bone. While the standard micro-FE model overestimated the bone-implant stiffness, the consideration of insertion-related bone damage was crucial for an accurate stiffness prediction. This approach provides an important step toward more accurate specimen-specific micro-FE models. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:2415-2424, 2017.

Calculation of cancellous bone elastic properties with the polarization-based FFT iterative scheme

The Fast Fourier Transform-based method, originally introduced by Moulinec and Suquet in 1994 has gained popularity for computing homogenized properties of composites. In this work, the method is used for the computational homogenization of the elastic properties of cancellous bone. To the authors' knowledge, this is the first study where the Fast Fourier Transform scheme is applied to bone mechanics. The performance of the method is analyzed for artificial and natural bone samples of 2 species: bovine femoral heads and implanted femurs of Hokkaido rats. Model geometries are constructed using data from X-ray tomographies, and the bone tissue elastic properties are measured using microindentation and nanoindentation tests. Computed results are in excellent agreement with those available in the literature. The study shows the suitability of the method to accurately estimate the fully anisotropic elastic response of cancellous bone. Guidelines are provided for the construction of the models and the setting of the algorithm.

Joint loading and proximal tibia subchondral trabecular bone microarchitecture differ with walking gait patterns in end-stage knee osteoarthritis

OBJECTIVES

To (1) stratify patient subgroups according to their distinct walking gait patterns in end-stage knee osteoarthritis (OA); (2) compare measures of joint loading and proximal tibia subchondral trabecular bone (STB) microarchitecture among these gait subgroups.

DESIGN

Twenty-five knee OA patients undergoing total knee arthroplasty (TKA) had pre-operative gait analysis. Following surgery, excised tibial plateaus were micro-CT-scanned and STB microarchitecture analysed in four tibial condylar regions of interest. Peak knee moments were input to k-means cluster analysis, to identify subgroups with homogeneous gait patterns. Joint loading and STB microarchitecture parameters were compared among gait subgroups (Kruskal-Wallis, Bonferroni-corrected Mann-Whitney U tests).

RESULTS

Three gait subgroups were revealed: biphasics (n = 7), flexors (n = 9), counter-rotators (n = 9). Peak knee adduction moment (KAM) and KAM impulse were significantly higher (P < 0.05) in biphasics than in flexors and counter-rotators (KAM = -0.65, -0.40 and -0.21 Nm/kg, respectively), suggesting a higher medial-to-lateral tibiofemoral load ratio in biphasics. Interestingly, STB medial-to-lateral bone volume fraction (BV/TV) ratio was also significantly higher (more than double) in biphasics and flexors than in counter-rotators (2.24, 2.00 and 1.00, respectively), whereas in biphasics it was only 10% higher than in flexors and not significantly so.

CONCLUSIONS

Within the confines of the limited sample size, data suggests that different mechanisms between the biphasic and flexor gait subroups may generate comparable loads upon the tibial plateau and corresponding bony responses, despite significantly lower KAM indices in flexors. Hence, in flexor gait OA patients, conservative treatments designed to reduce KAM, may not be appropriate. Understanding joint loading among walking gait patterns and relationships to bone microarchitecture may aid at identifying/improving management of persons at risk for developing knee OA.

Systematic mapping of the subchondral bone 3D microarchitecture in the human tibial plateau: Variations with joint alignment

Tibial subchondral bone plays an important role in knee osteoarthritis (OA). Microarchitectural characterization of subchondral bone plate (SBP), underlying subchondral trabecular bone (STB) and relationships between these compartments, however, is limited. The aim of this study was to characterize the spatial distribution of SBP thickness, SBP porosity and STB microarchitecture, and relationships among them, in OA tibiae of varying joint alignment. Twenty-five tibial plateaus from end-stage knee-OA patients, with varus (n = 17) or non-varus (n = 8) alignment were micro-CT scanned (17 μm/voxel). SBP and STB microarchitecture was quantified via a systematic mapping in 22 volumes of interest per knee (11 medial, 11 lateral). Significant within-condylar and between-condylar (medial vs. lateral) differences (p < 0.05) were found. In varus, STB bone volume fraction (BV/TV) was consistently high throughout the medial condyle, whereas in non-varus, medially, it was more heterogeneously distributed. Regions of high SBP thickness were co-located with regions of high STB BV/TV underneath. In varus, BV/TV was significantly higher medially than laterally, however, not so in non-varus. Moreover, region-specific significant associations between the SBP thickness and SBP porosity and the underlying STB microarchitecture were detected, which in general were not captured when considering the values averaged for each condyle. As subchondral bone changes reflect responses to local mechanical and biochemical factors within the joint, our results suggest that joint alignment influences both the medial-to-lateral and the within-condyle distribution of force across the tibia, generating corresponding local bony responses (adaptation) of both the subchondral bone plate and underlying subchondral trabecular bone microarchitecture. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:1927-1941, 2017.

Influence of meniscus on cartilage and subchondral bone features of knees from older individuals: A cadaver study

Objective

Cartilage and subchondral bone form a functional unit. Here, we aimed to examine the effect of meniscus coverage on the characteristics of this unit in knees of older individuals.

Methods

We assessed the hyaline cartilage, subchondral cortical plate (SCP), and subchondral trabecular bone in areas covered or uncovered by the meniscus from normal cadaver knees (without degeneration). Bone cores harvested from the medial tibial plateau at locations uncovered (central), partially covered (posterior), and completely covered (peripheral) by the meniscus were imaged by micro-CT. The following were measured on images: cartilage volume (Cart.Vol, mm³) and thickness (Cart.Th, mm); SCP thickness (SCP.Th, μm) and porosity (SCP.Por, %); bone volume to total volume fraction (BV/TV, %); trabecular thickness (Tb.Th, μm), spacing (Tb.Sp, μm), and number (Tb.N, 1/mm); structure model index (SMI); trabecular pattern factor (Tb.Pf); and degree of anisotropy (DA).

Results

Among the 28 specimens studied (18 females) from individuals with mean age 82.8±10.2 years, cartilage and SCP were thicker at the central site uncovered by the meniscus than the posterior and peripheral sites, and Cart.Vol was greater. SCP.Por was highest in posterior samples. In the upper 1–5 mm of subchondral bone, central samples were characterized by higher values for BV/TV, Tb.N, Tb.Th, and connectivity (Tb.Pf), a more plate-like trabecular structure and lower anisotropy than with other samples. Deeper down, at 6–10 mm, the differences were slightly higher for Tb.Th centrally, DA peripherally and SMI posteriorly.

Conclusions

The coverage or not by meniscus in the knee of older individuals is significantly associated with Cart.Th, SCP.Th, SCP.Por and trabecular microarchitectural parameters in the most superficial 5 mm and to a lesser extent the deepest area of subchondral trabecular bone. These results suggest an effect of differences in local loading conditions. In subchondral bone uncovered by the meniscus, the trabecular architecture resembles that of highly loaded areas.

Effect of including damage at the tissue level in the nonlinear homogenisation of trabecular bone

Being able to predict bone fracture or implant stability needs a proper constitutive model of trabecular bone at the macroscale in multiaxial, non-monotonic loading modes. Its macroscopic damage behaviour has been investigated experimentally in the past, mostly with the restriction of uniaxial cyclic loading experiments for different samples, which does not allow for the investigation of several load cases in the same sample as damage in one direction may affect the behaviour in other directions. Homogenised finite element models of whole bones have the potential to assess complicated scenarios and thus improve clinical predictions. The aim of this study is to use a homogenisation-based multiscale procedure to upscale the damage behaviour of bone from an assumed solid phase constitutive law and investigate its multiaxial behaviour for the first time. Twelve cubic specimens were each submitted to nine proportional strain histories by using a parallel code developed in-house. Evolution of post-elastic properties for trabecular bone was assessed for a small range of macroscopic plastic strains in these nine load cases. Damage evolution was found to be non-isotropic, and both damage and hardening were found to depend on the loading mode (tensile, compression or shear); both were characterised by linear laws with relatively high coefficients of determination. It is expected that the knowledge of the macroscopic behaviour of trabecular bone gained in this study will help in creating more precise continuum FE models of whole bones that improve clinical predictions.

Optimizing finite element predictions of local subchondral bone structural stiffness using neural network-derived density-modulus relationships for proximal tibial subchondral cortical and trabecular bone

BACKGROUND

Quantitative computed tomography based subject-specific finite element modeling has potential to clarify the role of subchondral bone alterations in knee osteoarthritis initiation, progression, and pain. However, it is unclear what density-modulus equation(s) should be applied with subchondral cortical and subchondral trabecular bone when constructing finite element models of the tibia. Using a novel approach applying neural networks, optimization, and back-calculation against in situ experimental testing results, the objective of this study was to identify subchondral-specific equations that optimized finite element predictions of local structural stiffness at the proximal tibial subchondral surface.

METHODS

Thirteen proximal tibial compartments were imaged via quantitative computed tomography. Imaged bone mineral density was converted to elastic moduli using multiple density-modulus equations (93 total variations) then mapped to corresponding finite element models. For each variation, root mean squared error was calculated between finite element prediction and in situ measured stiffness at 47 indentation sites. Resulting errors were used to train an artificial neural network, which provided an unlimited number of model variations, with corresponding error, for predicting stiffness at the subchondral bone surface. Nelder-Mead optimization was used to identify optimum density-modulus equations for predicting stiffness.

FINDINGS

Finite element modeling predicted 81% of experimental stiffness variance (with 10.5% error) using optimized equations for subchondral cortical and trabecular bone differentiated with a 0.5g/cm³ density.

INTERPRETATION

In comparison with published density-modulus relationships, optimized equations offered improved predictions of local subchondral structural stiffness. Further research is needed with anisotropy inclusion, a smaller voxel size and de-blurring algorithms to improve predictions.

The effects of cracks on the quantification of the cancellous bone fabric tensor in fossil and archaeological specimens: a simulation study

Cancellous bone is very sensitive to its prevailing mechanical environment, and study of its architecture has previously aided interpretations of locomotor biomechanics in extinct animals or archaeological populations. However, quantification of architectural features may be compromised by poor preservation in fossil and archaeological specimens, such as post mortem cracking or fracturing. In this study, the effects of post mortem cracks on the quantification of cancellous bone fabric were investigated through the simulation of cracks in otherwise undamaged modern bone samples. The effect on both scalar (degree of fabric anisotropy, fabric elongation index) and vector (principal fabric directions) variables was assessed through comparing the results of architectural analyses of cracked vs. non-cracked samples. Error was found to decrease as the relative size of the crack decreased, and as the orientation of the crack approached the orientation of the primary fabric direction. However, even in the best-case scenario simulated, error remained substantial, with at least 18% of simulations showing a > 10% error when scalar variables were considered, and at least 6.7% of simulations showing a > 10° error when vector variables were considered. As a 10% (scalar) or 10° (vector) difference is probably too large for reliable interpretation of a fossil or archaeological specimen, these results suggest that cracks should be avoided if possible when analysing cancellous bone architecture in such specimens.

Nonlinear homogenisation of trabecular bone: Effect of solid phase constitutive model

Micro-finite element models have been extensively employed to evaluate the elastic properties of trabecular bone and, to a limited extent, its yield behaviour. The macroscopic stiffness tensor and yield surface are of special interest since they are essential in the prediction of bone strength and stability of implants at the whole bone level. While macroscopic elastic properties are now well understood, yield and post-yield properties are not. The aim of this study is to shed some light on what the effect of the solid phase yield criterion is on the macroscopic yield of trabecular bone for samples with different microstructure. Three samples with very different density were subjected to a large set of apparent load cases (which is important since physiological loading is complex and can have multiple components in stress or strain space) with two different solid phase yield criteria: Drucker-Prager and eccentric-ellipsoid. The study found that these two criteria led to small differences in the macroscopic yield strains for most load cases except for those that were compression-dominated; in these load cases, the yield strains for the Drucker-Prager criterion were significantly higher. Higher density samples resulted in higher differences between the two criteria. This work provides a comprehensive assessment of the effect of two different solid phase yield criteria on the macroscopic yield strains of trabecular bone, for a wide range of load cases, and for samples with different morphology.

An evolutionary model of osteoarthritis including articular cartilage damage, and bone remodeling in a computational study

With osteoarthritis, a complex set of progressive chemical, biological, and mechanical changes occur in both cartilage and bone. The aim of this study is to develop a high-fidelity computational model of the complete bone-cartilage unit to study the evolution of osterarthritis-induced articular cartilage (AC) damage and remodeling of subchondral cortical bone (SCB) and subchondral trabecular bone (STB). A finite element model of spherical indentation was developed with a depth-dependent anisotropic model of degenerating articular cartilage, a calcified cartilage (CC) zone, and SCB and STB remodeling regions. Calcified tissue (CC, SCB, and STB) and AC material regions were integrated to form an evolutionary bone-cartilage unit model. Results indicate that with indentation loading, articular cartilage damage occurs at the articular surface. Furthermore, bone remodeling was predicted to occur with a net stiffening of the subchondral bone plate. Changes in indentation force were minimal (<2%) between initial and final peak indentation loading. However, additional degradation and wear of AC and/or alterations in loading may have more pronounced effects on the mechanical response of the bone-cartilage unit. Bone remodeling and articular cartilage damage predictions are consistent with experimental observations that cartilage damage begins at the articular surface and subchondral bone experiences a thickening (i.e., stiffening) response with osteoarthritis. Our results provide insight into the early-term initiation behavior of osteoarthritis; the potential consequences of evolutions in AC, SCB, and STB with disease progression; and may guide future experimental and computational studies to elucidate mechanisms of osteoarthritis progression.