Bone is a biomaterial undergoing continuous renewal. The renewal process is known as bone remodelling and is operated by bone-resorbing cells (osteoclasts) and bone-forming cells (osteoblasts). Both biochemical and biomechanical regulatory mechanisms have been identified in the interaction between osteoclasts and osteoblasts. Here we focus on an additional and poorly understood potential regulatory mechanism of bone cells, that involves the morphology of the microstructure of bone. Bone cells can only remove and replace bone at a bone surface. However, the microscopic availability of bone surface depends in turn on the ever-changing bone microstructure. The importance of this geometrical dependence is unknown and difficult to quantify experimentally. Therefore, we develop a sophisticated mathematical model of bone cell interactions that takes into account biochemical, biomechanical and geometrical regulations. We then investigate numerically the influence of bone surface availability in bone remodelling within a representative bone tissue sample. The interdependence between the bone cells' activity, which modifies the bone microstructure, and changes in the microscopic bone surface
For bone segmentation, the classical geodesic active contour model is usually limited by its indiscriminate feature extraction, and then struggles to handle the phenomena of edge obstruction, edge leakage and bone fracture. Thus, we propose a fracture interactive geodesic active contour algorithm tailored for bone segmentation, which can better capture bone features and perform robustly to the presence of bone fractures and soft tissues. Inspired by orthopedic knowledge, we construct a novel edge-detector function that combines the intensity and gradient norm, which guides the contour towards bone edges without being obstructed by other soft tissues and therefore reduces mis-segmentation. Furthermore, distance information, where fracture prompts can be embedded, is introduced into the contour evolution as an adaptive step size to stabilize the evolution and help the contour stop at bone edges and fractures. This embedding provides a way to interact with bone fractures and improves the accuracy in the fracture regions. Experiments in pelvic and ankle segmentation demonstrate the effectiveness on addressing the aforementioned problems and show an accurate, stable and consistent perfo
Osteoarthritis (OA) is a multifaceted joint disease which poses significant socioeconomic burdens and remains a significant clinical challenge. Evidence suggests that structural and mechanical changes in subchondral bone influence the pathogenesis and development of OA, leading to diminished bone quality and cartilage degeneration. While changes in microstructure and tissue scale elastic properties are well reported, the tissue yield response of subchondral bone in OA and their correlation with compositional changes have not been investigated. Here, we performed quasistatic micropillar compression and nanoindentation within the subchondral bone plate and trabeculae of hydrated non-diseased (ND) and OA affected specimens retrieved from the distal tibia in vivo. The micropillars, extracted by laser ablation, exhibited a taper angle which mandated the use of an in silico micropillar compression routine to back-calculate elastic modulus and strength of the bone tissue that comprised each micropillar. Elastic modulus remained unchanged between ND and OA subchondral bone, whereas strength increased from 46.0 MPa to 57.3 MPa in OA subchondral trabecular bone but not in the bone plate. Mic
Ultrasound imaging of the cortex of long bones may enable the measurement of the cortical thickness and the ultrasound wave speed in cortical bone tissue. However, with bone loss, the cortical porosity and the size of the pores increase, resulting in strong ultrasound diffuse scattering whose magnitude can exceed that of the specular reflection from the bone-marrow (endosteal) interface. In this study we adapt to bone a specular beamforming technique proposed to better image a needle in soft tissue. Our approach takes into account both wave refraction and specular reflection physics to enhance the contrast of bone surfaces and reduce speckle from intracortical pores. In vivo ultrasound data were acquired at the center of the human tibia in a plane normal to the bone axis of 11 young healthy volunteeers. Ex vivo ultrasound data were acquired from 16 regions of interest from the femoral diaphysis of three elderly donors (donors 66-98 y.o.) using a 2.5 MHz US transducer. A single-element trans mission synthetic aperture imaging sequence was implemented on a research ultrasound system with a 2.5MHz phased array transducer. Image reconstruction was performed: (A) a delay-and-sum (DAS) a
Current approaches in 3D human pose estimation primarily focus on regressing 3D joint locations, often neglecting critical physical constraints such as bone length consistency and body symmetry. This work introduces a recurrent neural network architecture designed to capture holistic information across entire video sequences, enabling accurate prediction of bone lengths. To enhance training effectiveness, we propose a novel augmentation strategy using synthetic bone lengths that adhere to physical constraints. Moreover, we present a bone length adjustment method that preserves bone orientations while substituting bone lengths with predicted values. Our results demonstrate that existing 3D human pose estimation models can be significantly enhanced through this adjustment process. Furthermore, we fine-tune human pose estimation models using inferred bone lengths, observing notable improvements. Our bone length prediction model surpasses the previous best results, and our adjustment and fine-tuning method enhance performance across several metrics on the Human3.6M dataset.
Age-related bone loss and postmenopausal osteoporosis are disorders of bone remodelling, in which less bone is reformed than resorbed. Yet, this dysregulation of bone remodelling does not occur equally in all bone regions. Loss of bone is more pronounced near and at the endocortex, leading to cortical wall thinning and medullary cavity expansion, a process sometimes referred to as "trabecularisation" or "cancellisation". Cortical wall thinning is of primary concern in osteoporosis due to the strong deterioration of bone mechanical properties that it is associated with. In this paper, we examine the possibility that the non-uniformity of microscopic bone surface availability could explain the non-uniformity of bone loss in osteoporosis. We use a computational model of bone remodelling in which microscopic bone surface availability influences bone turnover rate and simulate the evolution of the bone volume fraction profile across the midshaft of a long bone. We find that bone loss is accelerated near the endocortical wall where the specific surface is highest. Over time, this leads to a substantial reduction of cortical wall thickness from the endosteum. The associated expansion of t
This work proposes a gradient-based method to design bone implants using triply-periodic minimal surfaces (TPMS) of spatially varying thickness to maximize bone in-growth. Bone growth into the implant is estimated using a finite element based mechanobiological model considering the magnitude and frequency of in vivo loads, as well as the density distribution of the surrounding bone. The wall thicknesses of the implant unit cells are determined via linear interpolation of the thicknesses over a user defined grid of control points, avoiding mesh dependency and providing control over the sensitivity computation costs. The TPMS structure is modeled as a homogenized material to reduce computational cost. Local properties of the implant are determined at run-time on an element-by-element basis using a pre-constructed surrogate model of the TPMS's physical and geometric properties as a function of the local wall thickness and the density of in-grown bone. Design sensitivities of the bone growth within the implant are computed using the direct sensitivity method. The methodology is demonstrated on a cementless hip, optimizing the implant for bone growth subject to wall thickness constraint
Bone metastasis analysis is a significant challenge in pathology and plays a critical role in determining patient quality of life and treatment strategies. The microenvironment and specific tissue structures are essential for pathologists to predict the primary bone cancer origins and primary bone cancer subtyping. By digitizing bone tissue sections into whole slide images (WSIs) and leveraging deep learning to model slide embeddings, this analysis can be enhanced. However, tumor metastasis involves complex multivariate interactions with diverse bone tissue structures, which traditional WSI analysis methods such as multiple instance learning (MIL) fail to capture. Moreover, graph neural networks (GNNs), limited to modeling pairwise relationships, are hard to represent high-order biological associations. To address these challenges, we propose a dynamic hypergraph neural network (DyHG) that overcomes the edge construction limitations of traditional graph representations by connecting multiple nodes via hyperedges. A low-rank strategy is used to reduce the complexity of parameters in learning hypergraph structures, while a Gumbel-Softmax-based sampling strategy optimizes the patch di
Photoacoustic (PA) technology can provide information on both the physical structure and chemical composition of bone, showing great potential in bone assessment. However, due to the complex composition and porous structure of cancellous bone, the PA signals generated and propagated in cancellous bone are complex and difficult to be directly used in cancellous bone analysis. In this paper, a photoacoustic differential attenuation spectrum (PA-DAS) method is proposed. By eliminating the PA spectrum of the optical absorption sources, the propagation attenuation characteristics of cancellous bone are studied theoretically and experimentally. An analytical solution for the propagation attenuation of broadband ultrasound waves in cancellous bone is given by applying high-frequency and viscous corrections to Biot's theory. An experimental system of PA-DAS with an eccentric excitation differential detection system is established to obtain the PA-DAS of cancellous bone and its acoustic propagation characteristic on the rabbit osteoporosis model. The PA-DAS quantization parameter slope is further extracted to quantify the attenuation of high and low frequency components. The results show th
Bone remodeling involves the coordinated removal of bone by osteoclasts and addition of bone by osteoblasts, a process that is modulated by the prevailing mechanical environment. In this paper a fully coupled model of bone remodeling is developed, based on coupling a bone cell population model with a micromechanical homogenization scheme of bone stiffness. While the former model considers biochemical regulatory mechanisms between bone cells such as the RANK-RANKL-OPG pathway and action of TGF-β, the latter model allows for accurate upscaling of the mechanical properties of bone. Importantly, we consider bone remodeling as being controlled proportionally to the microscopic strain energy density, on the observation scale where the sensing of the mechanical loading takes place, estimated by means of continuum micromechanics-based strain concentration. This approach allows to address two fundamental questions of bone biology: (i) How do biochemical changes influence bone remodeling and so affect the composition and mechanical properties of bone? and (ii) What mechanisms are responsible for mechanoregulation of bone remodeling? Numerical studies highlight the conceptual advantage of thi
The skeletal region is one of the common sites of metastatic spread of cancer in the breast and prostate. CT is routinely used to measure the size of lesions in the bones. However, they can be difficult to spot due to the wide variations in their sizes, shapes, and appearances. Precise localization of such lesions would enable reliable tracking of interval changes (growth, shrinkage, or unchanged status). To that end, an automated technique to detect bone lesions is highly desirable. In this pilot work, we developed a pipeline to detect bone lesions (lytic, blastic, and mixed) in CT volumes via a proxy segmentation task. First, we used the bone lesions that were prospectively marked by radiologists in a few 2D slices of CT volumes and converted them into weak 3D segmentation masks. Then, we trained a 3D full-resolution nnUNet model using these weak 3D annotations to segment the lesions and thereby detected them. Our automated method detected bone lesions in CT with a precision of 96.7% and recall of 47.3% despite the use of incomplete and partial training data. To the best of our knowledge, we are the first to attempt the direct detection of bone lesions in CT via a proxy segmentat
Short and long-term stabilities of cementless implants are strongly determined by the interfacial load transfer between implants and bone tissue. Stress-shielding effects arise from shear stresses due to the difference of material properties between bone and the implant. It remains difficult to measure the stress field in periprosthetic bone tissue. This study proposes to investigate the dependence of the stress field in periprosthetic bone tissue on i) the implant surface roughness, ii) material properties of bone and of the implant, iii) the bone-implant contact ratio. To do so, a microscale 2-D finite element model of an osseointegrated boneimplant interface was developed where the surface roughness was modeled by a sinusoidal surface. The results show that the isostatic pressure is not affected by the presence of the bone-implant interface while shear stresses arise due to the combined effects of a geometrical singularity (for low surface roughness) and of shear stresses at the bone-implant interface (for high surface roughness). Stress-shielding effects are likely to be more important when the bone-implant contact ratio value is low, which corresponds to a case of relatively l
Bone health studies are crucial in medical practice for the early detection and treatment of Osteopenia and Osteoporosis. Clinicians usually make a diagnosis based on densitometry (DEXA scans) and patient history. The applications of AI in this field are ongoing research. Most successful methods rely on deep learning models that use vision alone (DEXA/X-ray imagery) and focus on prediction accuracy, while explainability is often disregarded and left to post hoc assessments of input contributions. We propose ProtoMedX, a multi-modal (multimodal) model that uses both DEXA scans of the lumbar spine and patient records. ProtoMedX's prototype-based architecture is explainable by design, which is crucial for medical applications, especially in the context of the upcoming EU AI Act, as it allows explicit analysis of model decisions, including incorrect ones. ProtoMedX demonstrates state-of-the-art performance in bone health classification while also providing explanations that can be visually understood by clinicians. Using a dataset of 4,160 real NHS patients, the proposed ProtoMedX achieves 87.58% accuracy in vision-only tasks and 89.8% in its multi-modal variant, both surpassing existi
For a long time, bone marrow cell morphology examination has been an essential tool for diagnosing blood diseases. However, it is still mainly dependent on the subjective diagnosis of experienced doctors, and there is no objective quantitative standard. Therefore, it is crucial to study a robust bone marrow cell detection algorithm for a quantitative automatic analysis system. Currently, due to the dense distribution of cells in the bone marrow smear and the diverse cell classes, the detection of bone marrow cells is difficult. The existing bone marrow cell detection algorithms are still insufficient for the automatic analysis system of bone marrow smears. This paper proposes a bone marrow cell detection algorithm based on the YOLOv5 network, trained by minimizing a novel loss function. The classification method of bone marrow cell detection tasks is the basis of the proposed novel loss function. Since bone marrow cells are classified according to series and stages, part of the classes in adjacent stages are similar. The proposed novel loss function considers the similarity between bone marrow cell classes, increases the penalty for prediction errors between dissimilar classes, and
The formation of new bone involves both the deposition of bone matrix, and the formation of a network of cells embedded within the bone matrix, called osteocytes. Osteocytes derive from bone-synthesising cells (osteoblasts) that become buried in bone matrix during bone deposition. The generation of osteocytes is a complex process that remains incompletely understood. Whilst osteoblast burial determines the density of osteocytes, the expanding network of osteocytes regulates in turn osteoblast activity and osteoblast burial. In this paper, a spatiotemporal continuous model is proposed to investigate the osteoblast-to-osteocyte transition. The aims of the model are (i) to link dynamic properties of osteocyte generation with properties of the osteocyte network imprinted in bone, and (ii) to investigate Marotti's hypothesis that osteocytes prompt the burial of osteoblasts when they become covered with sufficient bone matrix. Osteocyte density is assumed in the model to be generated at the moving bone surface by a combination of osteoblast density, matrix secretory rate, rate of entrapment, and curvature of the bone substrate, but is found to be determined solely by the ratio of the ins
Until recently many studies of bone remodeling at the cellular level have focused on the behavior of mature osteoblasts and osteoclasts, and their respective precursor cells, with the role of osteocytes and bone lining cells left largely unexplored. This is particularly true with respect to the mathematical modeling of bone remodeling. However, there is increasing evidence that osteocytes play important roles in the cycle of targeted bone remodeling, in serving as a significant source of RANKL to support osteoclastogenesis, and in secreting the bone formation inhibitor sclerostin. Moreover, there is also increasing interest in sclerostin, an osteocyte-secreted bone formation inhibitor, and its role in regulating local response to changes in the bone microenvironment. Here we develop a cell population model of bone remodeling that includes the role of osteocytes, sclerostin, and allows for the possibility of RANKL expression by osteocyte cell populations. This model extends and complements many of the existing mathematical models for bone remodeling but can be used to explore aspects of the process of bone remodeling that were previously beyond the scope of prior modeling work. Thro
The Ishango Bone is a prehistoric artifact dated to approximately 20,000 years ago, discovered near the Semliki River in what is now the Democratic Republic of Congo, and has been the subject of scholarly debate for decades. The artifact displays sixteen groups of notches organised into three distinct columns, a structure that permits relational analysis, yet its precise function remains debated. This study identifies previously undescribed mathematical patterns across all three columns. Two columns comprise all prime and odd numbers between 9 and 21, with the sole exception of 15, itself the arithmetic mean of this set. Each column sums to 60 and subdivides into two internal groupings, each summing to 30. Exploratory positional adjustments, each uniquely constrained by the data, appear to reveal a consistent grouping rule and arithmetic relationships, spanning all three columns. Five structural properties are evaluated simultaneously through a global permutation test. The fully adjusted configuration achieves a score that is not observed among the 1,000,000 random rearrangements of the same values. The findings in this study support the hypothesis that the Ishango Bone may have fu
Bone drilling is a universal surgical procedure commonly used for internal fracture fixation, implant placement, or reconstructive surgery in orthopedics and dentistry. The increased temperature during such treatment increases the risk of thermal penetration of the bone, which may delay healing or compromise the fixation's integrity. Thus, avoiding penetration during bone drilling is critical to ensuring the implant's stability, which needs surgical drills with an optimized design. Bovine femur and mandible bones are chosen as the work material since human bones are not available and are the closest animal bone to human bone in terms of properties. In the present study, the Taguchi fractional factorial approach was used to determine the best design of surgical drills by comparing the drilling properties (i.e., signal-to-noise ratio and temperature rise). The control factors (spindle speed, drill bit diameter, drill site depth, and their levels) were arranged in an L9 orthogonal array. Drilling experiments were done using nine experimental drills with three repetitions. The findings of this study indicate that the ideal values of the surgical drill's three parameters combination (S1
Segmenting both bone surface and the corresponding acoustic shadow are fundamental tasks in ultrasound (US) guided orthopedic procedures. However, these tasks are challenging due to minimal and blurred bone surface response in US images, cross-machine discrepancy, imaging artifacts, and low signal-to-noise ratio. Notably, bone shadows are caused by a significant acoustic impedance mismatch between the soft tissue and bone surfaces. To leverage this mutual information between these highly related tasks, we propose a single end-to-end network with a shared transformer-based encoder and task independent decoders for simultaneous bone and shadow segmentation. To share complementary features, we propose a cross task feature transfer block which learns to transfer meaningful features from decoder of shadow segmentation to that of bone segmentation and vice-versa. We also introduce a correspondence consistency loss which makes sure that network utilizes the inter-dependency between the bone surface and its corresponding shadow to refine the segmentation. Validation against expert annotations shows that the method outperforms the previous state-of-the-art for both bone surface and shadow s
Continuum bone remodelling is an important tool for predicting the effects of mechanical stimuli on bone density evolution. While the modelling of only cancellous bone is considered in many studies based on continuum bone remodelling, this work presents an approach of modelling also cortical bone and the interaction of both bone types. The distinction between bone types is made by introducing an initial volume fraction. A simple point-wise example is used to study the behaviour of novel model options, as well as a proximal femur example, where the interaction of both bone types is demonstrated using initial density distributions. The results of the proposed model options indicate that the consideration of cortical bone remarkably changes the density evolution of cancellous bone, and should therefore not be neglected.