Diversity of the proximal femur in humans: morphological variations of the head-neck junction

The proximal femur is a high-diversity region of the human skeleton, especially at the anterior junction between head and neck, where various bony morphologies have been recognized since mid nineteenth century. Classical literature on this topic is chaotic and contradictory, making almost impossible the comparison of data from different researches. Starting from an extensive bibliographic review, the first standardized method to score these traits has been created. This method allows representing both the anatomical diversity of the region already described in literature and a part of variability not considered before, giving few and univocal definitions and allowing to collect comparable data. The method has been applied to three identified and five archaeological European skeletal collections, with the aim of investigating the distribution of these features by sex, age and side, in different places and time periods. It has also been applied to 3D digital reconstructions of femurs from CT scan files of coxo-femoral joints from fresh cadavers. In addition to the osseous traits described in the standardized method, the presence and frequency of some features known as herniation pits have been scored both on bones and on CT scans. The various osseous traits of the proximal femur are present at similar frequencies in skeletal samples from different countries and different historical periods, even if with clear local differentiation. Some of the features examined show significant trends related to their distribution by gender and age. Some hypotheses are proposed about the etiology of these morphologies and their possible implication with the acquisition of bipedalism in Humans. It is therefore highlighted the possible relation of some of these traits with the development of disorders of the hip joint. Moreover, it is not recommended the use of any of these features as a specific activity-related marker.

Mesh-based vs. Image-based Statistical Appearance Model of the Human Femur: a Preliminary Comparison Study for the Creation of Finite Element Meshes

Statistical models have been recently introduced in computational orthopaedics to investigate the bone mechanical properties across several populations. A fundamental aspect for the construction of statistical models concerns the establishment of accurate anatomical correspondences among the objects of the training dataset. Various methods have been proposed to solve this problem such as mesh morphing or image registration algorithms. The objective of this study is to compare a mesh-based and an image-based statistical appearance model approaches for the creation of nite element(FE) meshes. A computer tomography (CT) dataset of 157 human left femurs was used for the comparison. For each approach, 30 finite element meshes were generated with the models. The quality of the obtained FE meshes was evaluated in terms of volume, size and shape of the elements. Results showed that the quality of the meshes obtained with the image-based approach was higher than the quality of the mesh-based approach. Future studies are required to evaluate the impact of this finding on the final mechanical simulations.

Femur Specific Polyaffine Model to Regularize the Log-Domain Demons Registration

Osteoarticular allograft transplantation is a popular treatment method in wide surgical resections with large defects. For this reason hospitals are building bone data banks. Performing the optimal allograft selection on bone banks is crucial to the surgical outcome and patient recovery. However, current approaches are very time consuming hindering an efficient selection. We present an automatic method based on registration of femur bones to overcome this limitation. We introduce a new regularization term for the log-domain demons algorithm. This term replaces the standard Gaussian smoothing with a femur specific polyaffine model. The polyaffine femur model is constructed with two affine (femoral head and condyles) and one rigid (shaft) transformation. Our main contribution in this paper is to show that the demons algorithm can be improved in specific cases with an appropriate model. We are not trying to find the most optimal polyaffine model of the femur, but the simplest model with a minimal number of parameters. There is no need to optimize for different number of regions, boundaries and choice of weights, since this fine tuning will be done automatically by a final demons relaxation step with Gaussian smoothing. The newly developed synthesis approach provides a clear anatomically motivated modeling contribution through the specific three component transformation model, and clearly shows a performance improvement (in terms of anatomical meaningful correspondences) on 146 CT images of femurs compared to a standard multiresolution demons. In addition, this simple model improves the robustness of the demons while preserving its accuracy. The ground truth are manual measurements performed by medical experts.

Minimization of Intra-Operative Shaping of Orthopaedic Fixation Plates: A Population-Based Design

In this paper we present a new population-based method for the design of bone fixation plates. Standard pre-contoured plates are designed based on the mean shape of a certain population. We propose a computational process to design implants while reducing the amount of required intra-operative shaping, thus reducing the mechanical stresses applied to the plate. A bending and torsion model was used to measure and minimize the necessary intra-operative deformation. The method was applied and validated on a population of 200 femurs that was further augmented with a statistical shape model. The obtained results showed substantial reduction in the bending and torsion needed to shape the new design into any bone in the population when compared to the standard mean-based plates.

Personalized X-Ray Reconstruction of the Proximal Femur via Intensity-Based Non-rigid 2D-3D Registration

This paper presents a new approach for reconstructing a patient-specific shape model and internal relative intensity distribution of the proximal femur from a limited number (e.g., 2) of calibrated C-arm images or X-ray radiographs. Our approach uses independent shape and appearance models that are learned from a set of training data to encode the a priori information about the proximal femur. An intensity-based non-rigid 2D-3D registration algorithm is then proposed to deformably fit the learned models to the input images. The fitting is conducted iteratively by minimizing the dissimilarity between the input images and the associated digitally reconstructed radiographs of the learned models together with regularization terms encoding the strain energy of the forward deformation and the smoothness of the inverse deformation. Comprehensive experiments conducted on images of cadaveric femurs and on clinical datasets demonstrate the efficacy of the present approach.

Reconstruction of patient-specific 3D bone surface from 2D calibrated fluoroscopic images and point distribution model

Reconstruction of patient-specific 3D bone surface from 2D calibrated fluoroscopic images and a point distribution model is discussed. We present a 2D/3D reconstruction scheme combining statistical extrapolation and regularized shape deformation with an iterative image-to-model correspondence establishing algorithm, and show its application to reconstruct the surface of proximal femur. The image-to-model correspondence is established using a non-rigid 2D point matching process, which iteratively uses a symmetric injective nearest-neighbor mapping operator and 2D thin-plate splines based deformation to find a fraction of best matched 2D point pairs between features detected from the fluoroscopic images and those extracted from the 3D model. The obtained 2D point pairs are then used to set up a set of 3D point pairs such that we turn a 2D/3D reconstruction problem to a 3D/3D one. We designed and conducted experiments on 11 cadaveric femurs to validate the present reconstruction scheme. An average mean reconstruction error of 1.2 mm was found when two fluoroscopic images were used for each bone. It decreased to 1.0 mm when three fluoroscopic images were used.

Comparison of six radiographic projections to assess femoral head/neck asphericity

Early radiographic detection of femoroacetabular impingement might prevent initiation and progression of osteoarthritis. The structural abnormality in femoral-induced femoroacetabular impingement (cam type) is frequently asphericity at the anterosuperior head/neck contour. To determine which of six radiographic projections (anteroposterior, Dunn, Dunn/45 degrees flexion, cross-table/15 degrees internal rotation, cross-table/neutral rotation, and cross-table/15 degrees external rotation) best identifies femoral head/neck asphericity, we studied 21 desiccated femurs; 11 with an aspherical femoral head/neck contour and 10 with a spherical femoral head/neck contour. To radiographically quantify femoral head asphericity, we measured the angle where the femoral head/neck leaves sphericity (angle alpha). The aspherical femoral head/neck contours had a greater maximum angle alpha (70 degrees ) compared with the spherical head/neck contours (50 degrees ). The angle alpha varied depending on the radiographic projection: it was greatest in the Dunn view with 45 degrees hip flexion (71 degrees +/- 10 degrees ) and least in the cross-table view in 15 degrees external rotation (51 degrees +/- 7 degrees ). Diagnosis of a pathologic femoral head/neck contour depends on the radiologic projection. The Dunn view in 45 degrees or 90 degrees flexion or a cross-table projection in internal rotation best show femoral head/neck asphericity, whereas anteroposterior or externally rotated cross-table views are likely to miss asphericity. Level of Evidence: Prognostic study, level II-1 (retrospective study).

Accurate and robust reconstruction of proximal femur from sparse intraoperative data and dense point distribution model for surgical navigation

Constructing a 3D surface model from sparse-point data is a nontrivial task. Here, we report an accurate and robust approach for reconstructing a surface model of the proximal femur from sparse-point data and a dense-point distribution model (DPDM). The problem is formulated as a three-stage optimal estimation process. The first stage, affine registration, is to iteratively estimate a scale and a rigid transformation between the mean surface model of the DPDM and the sparse input points. The estimation results of the first stage are used to establish point correspondences for the second stage, statistical instantiation, which stably instantiates a surface model from the DPDM using a statistical approach. This surface model is then fed to the third stage, kernel-based deformation, which further refines the surface model. Handling outliers is achieved by consistently employing the least trimmed squares (LTS) approach with a roughly estimated outlier rate in all three stages. If an optimal value of the outlier rate is preferred, we propose a hypothesis testing procedure to automatically estimate it. We present here our validations using four experiments, which include 1 leave-one-out experiment, 2 experiment on evaluating the present approach for handling pathology, 3 experiment on evaluating the present approach for handling outliers, and 4 experiment on reconstructing surface models of seven dry cadaver femurs using clinically relevant data without noise and with noise added. Our validation results demonstrate the robust performance of the present approach in handling outliers, pathology, and noise. An average 95-percentile error of 1.7-2.3 mm was found when the present approach was used to reconstruct surface models of the cadaver femurs from sparse-point data with noise added.

3D volumetric intensity reconsturction from 2D x-ray images using partial least squares regression

Reconstruction of shape and intensity from 2D x-ray images has drawn more and more attentions. Previously introduced work suffers from the long computing time due to its iterative optimization characteristics and the requirement of generating digitally reconstructed radiographs within each iteration. In this paper, we propose a novel method which uses a patient-specific 3D surface model reconstructed from 2D x-ray images as a surrogate to get a patient-specific volumetric intensity reconstruction via partial least squares regression. No DRR generation is needed. The method was validated on 20 cadaveric proximal femurs by performing a leave-one-out study. Qualitative and quantitative results demonstrated the efficacy of the present method. Compared to the existing work, the present method has the advantage of much shorter computing time and can be applied to both DXA images as well as conventional x-ray images, which may hold the potentials to be applied to clinical routine task such as total hip arthroplasty (THA).

Stromal derived growth factor alpha (SDF-1α) induces the differentiation of bone marrow progenitor cells (BMCs) into pericytes by regulating platelet derived growth factor B (PDGF-B) transcription

We previously demonstrated that bone marrow cells (BMCs) migrate to TC71 and A4573 Ewing’s sarcoma tumors where they can differentiate into endothelial cells (ECs) and pericytes and, participate in the tumor vascular development. This process of neo-vascularization, known as vasculogenesis, is essential for Ewing’s sarcoma growth with the soluble vascular endothelial growth factor, VEGF165, being the chemotactic factor for BMC migration to the tumor site. Inhibiting VEGF165 in TC71 tumors (TC/siVEGF7-1) inhibited BMC infiltration to the tumor site and tumor growth. Introducing the stromal-derived growth factor (SDF-1α) into the TC/siVEGF7-1 tumors partially restored vasculogenesis with infiltration of BMCs to a perivascular area where they differentiated into pericytes and rescued tumor growth. RNA collected from the SDF-1α-treated TC/siVEGF7-1 tumors also revealed an increase in platelet-derived growth factor B (PDGF-B) mRNA levels. PDGF-B expression is elevated in several cancer types and the role of PDGF-B and its receptor, PDGFR-β, has been extensively described in the process of pericyte maturation. However, the mechanisms by which PDGF-B expression is up-regulated during vascular remodeling and the process by which BMCs differentiate into pericytes during tumor vasculogenesis remain areas of investigation. In this study, we are the first to demonstrate that SDF-1α regulates the expression of PDGF-B via a transcriptional mechanism which involves binding of the ELK-1 transcription factor to the pdgf-b promoter. We are also first to validate the critical role of the SDF-1α/PDGF-B pathway in the differentiation of BMCs into pericytes both in vitro and in vivo. SDF-1α up-regulated PDGF-B expression in both TC/siVEGF7-1 and HEK293 cells. In contrast, down-regulating SDF-1α, down-regulated PDGF-B. We cloned the 2 kb pdgf-b promoter fragment into the pGL3 reporter vector and showed that SDF-1α induced pdgf-b promoter activity. We used chromatin immunoprecipitation (ChIP) and demonstrated that the ELK-1 transcription factor bound to the pdgf-b promoter in response to SDF-1α stimulation in both TC/siVEGF7-1 and HEK293 cells. We collected BMCs from the hind femurs of mice and cultured the cells in medium containing SDF-1α and PDGF-B and found that PDGFR-β+ BMCs differentiated into NG2 and desmin positive pericytes in vitro. In contrast, inhibiting SDF-1α and PDGF-B abolished this differentiation process. In vivo, we injected TC71 or A4573 tumor-bearing mice with the SDF-1α antagonist, AMD3100 and found that inhibiting SDF-1α signaling in the tumor microenvironment decreased the tumor microvessel density, decreased the tumor blood vessel perfusion and, increased tumor cell apoptosis. We then analyzed the effect of AMD3100 on vasculogenesis of Ewing’s sarcoma and found that BMCs migrated to the tumor site where they differentiated into ECs but, they did not form thick perivascular layers of NG2 and desmin positive pericytes. Finally, we stained the AMD3100-treated tumors for PDGF-B and showed that inhibiting SDF-1α signaling also inhibited PDGF-B expression. All together, these findings demonstrated that the SDF-1α/PDGF-B pathway plays a critical role in the formation of BM-derived pericytes during vasculogenesis of Ewing’s sarcoma tumors.

A dosimetric study of radionuclide therapy for bone marrow ablation

In a phase I clinical trial, six multiple myeloma patients, who were non-responsive to conventional therapy and were scheduled for bone marrow transplantation, received Holmium-166 ($\sp{166}$Ho) labeled to a bone seeking agent, DOTMP (1,4,7,10-tetraazacyclododecane-1,4,7,10-tetramethylene-phosphonic acid), for the purpose of bone marrow ablation. The specific aims of my research within this protocol were to evaluate the toxicity and efficacy of $\sp{166}$Ho DOTMP by quantifying the in vivo pharmacokinetics and radiation dosimetry, and by correlating these results to the biologic response observed. The reproducibility of pharmacokinetics from multiple injections of $\sp{166}$Ho DOTMP administered to these myeloma patients was demonstrated from both blood and whole body retention. The skeletal concentration of $\sp{166}$Ho DOTMP was heterogenous in all six patients: high in the ribs, pelvis, and lumbar vertebrae regions, and relatively low in the femurs, arms, and head.^ A novel technique was developed to calculate the radiation dose to the bone marrow in each skeletal ROI, and was applied to all six $\sp{166}$Ho DOTMP patients. Radiation dose estimates for the bone marrow calculated using the standard MIRD "S" factors were compared with the average values derived from the heterogenous distribution of activity in the skeleton (i.e., the regional technique). The results from the two techniques were significantly different; the average of the dose estimates from the regional technique were typically 30% greater. Furthermore, the regional technique provided a range of radiation doses for the entire marrow volume, while the MIRD "S" factors only provided a single value. Dose volume histogram analysis of data from the regional technique indicated a range of dose estimates that varied by a factor of 10 between the high dose and low dose regions. Finally, the observed clinical response of cells and abnormal proteins measured in bone marrow aspirates and peripheral blood samples were compared with radiation dose estimates for the bone marrow calculated from the standard and regional technique. The results showed the regional technique values correlated more closely to several clinical response parameters. (Abstract shortened by UMI.) ^

Characterization of osteopontin charge forms produced by osteoblasts

Osteopontin (OPN) is a highly-phosphorylated extracellular matrix protein localized in bone, kidney, placenta, T-lymphocytes, macrophages, smooth muscle of the vascular system, milk, urine, and plasma. In ROS 17/2.8 osteoblast-like osteosarcoma cells, 1,25-dihydroxyvitamin D3 [1,25(OH)2D 3] regulates OPN at the transcriptional level resulting in increased steady state mRNA levels and increased production of OPN protein, maximal at 48 hours. Using ROS 17/2.8 cells as an osteoblast model, OPN was purified from culture medium after three hour treatments of either vehicle (ethanol) or 1,25(OH)2D3 via barium citrate precipitation followed by immunoaffinity chromatography. ^ Here, further evidence of regulation of OPN by 1,25(OH)2D 3 at the posttranslational level is presented. Prior to the up-regulation of OPN at the transcriptional level, 1,25(OH)2D3 induces a shift in OPN isoelectric point (pI) detected on two-dimensional gels from pI 4.6 to pI 5.1. Loading equal amounts of [32P]-labeled OPN recovered from ROS 17/2.8 cells exposed to 1,25(OH)2D3 or vehicle alone for three hours reveals that the shift from pI 4.6 to 5.1 is the result of reduced phosphorylation. Using structural analogs to 1,25(OH) 2D3, analog AT [25-(OH)-16-ene-23-yne-D3], which triggers Ca2+ influx through voltage sensitive Ca2+ channels but does not bind to the vitamin D receptor, mimicked the OPN pI shift while analog BT [1,25(OH)2-22-ene-24-cyclopropyl-D 3], which binds to the vitamin D receptor but does not allow Ca 2+ influx, did not. Inclusion of the Ca2+ channel blocker nifedipine also blocks the charge shift conversion of OPN. Further analysis of the signaling pathway initiated by 1,25(OH)2D3 reveals that inhibition of the cyclic 3′,5′ -adenosine monophosphate-dependent kinase, protein kinase A, or inhibition of the cyclic 3′,5′-guanine monophosphate-dependent kinase, protein kinase G, also prevents the charge shift conversion. ^ Isolation of OPN from rat femurs and tibiae provides evidence for the existence of these two OPN charge forms in vivo, evidenced by differential migration on isoelectric focusing gels and sodium dodecyl sulfate-polyacrylamide gels. Peptide sequencing of rat long bone fractions revealed the presence of a presumed dentin specific protein, dentin matrix protein-1 (DMP-1). Western blot analysis confirmed the existence of DMP-1 in these fractions. ^ Using the OPN charge forms in functional assays, it was determined that the charge forms have differential roles in both cell surface and mineralization functions. In cell attachment assays and Ca2+ influx assays using PC-3 prostate cancer cells, the pI 5.1 charge form of OPN was found to permit binding and increase intracellular Ca2+ concentrations of PC-3 cells. The increase in intracellular Ca2+ concentration was found to be integrin αvβ3-dependent. In mineralization assays, the pI 4.6 charge form of OPN promoted hydroxyapatite formation, while the pI 5.1 charge form had improved Ca2+ binding ability. ^ In conclusion, these findings suggest that 1,25(OH) 2D3 regulates OPN not only at the transcriptional level, but also plays a role in determination of the OPN phosphorylation state. The latter involves a short term (less than three hours) treatment and is associated with membrane-initiated Ca2+ influx. Functional assays utilizing the two OPN charge forms reveal the dependence of OPN post-translational state on its function. ^

End caps prevent nail migration in elastic stable intramedullary nailing in paediatric femoral fractures: a biomechanical study using synthetic and cadaveric bones.

End caps are intended to prevent nail migration (push-out) in elastic stable intramedullary nailing. The aim of this study was to investigate the force at failure with and without end caps, and whether different insertion angles of nails and end caps would alter that force at failure. Simulated oblique fractures of the diaphysis were created in 15 artificial paediatric femurs. Titanium Elastic Nails with end caps were inserted at angles of 45°, 55° and 65° in five specimens for each angle to create three study groups. Biomechanical testing was performed with axial compression until failure. An identical fracture was created in four small adult cadaveric femurs harvested from two donors (both female, aged 81 and 85 years, height 149 cm and 156 cm, respectively). All femurs were tested without and subsequently with end caps inserted at 45°. In the artificial femurs, maximum force was not significantly different between the three groups (p = 0.613). Push-out force was significantly higher in the cadaveric specimens with the use of end caps by an up to sixfold load increase (830 N, standard deviation (SD) 280 vs 150 N, SD 120, respectively; p = 0.007). These results indicate that the nail and end cap insertion angle can be varied within 20° without altering construct stability and that the risk of elastic stable intramedullary nailing push-out can be effectively reduced by the use of end caps.

Osteotomy of the distal femur. Surgical technique using the LCP Pediatric Condylar Plate 90°.

OBJECTIVE Correction of all kind of deformities at the distal part of the femur (supracondylar). INDICATIONS Flexion, extension osteotomies, and varus or valgus, and external or internal rotation osteotomies, and shortening osteotomies of the distal femur or combined surgical procedures (e.g., extension and de-rotation osteotomy). CONTRAINDICATIONS Osteotomy through unknown bony process. SURGICAL TECHNIQUE LCP system provides angular stable fixation. POSTOPERATIVE MANAGEMENT Without concomitant surgical procedures of soft tissue (e.g., patellar tendon shortening), early functional rehabilitation is possible with immediate weight bearing (35 kg for small fragment plates and 70 kg for large fragment plates). RESULTS The surgical procedure is safe and is associated with few complications. Overall complication rate in this series of patients was 3%.

Fully Automatic Segmentation of Hip CT Images

Automatic segmentation of the hip joint with pelvis and proximal femur surfaces from CT images is essential for orthopedic diagnosis and surgery. It remains challenging due to the narrowness of hip joint space, where the adjacent surfaces of acetabulum and femoral head are hardly distinguished from each other. This chapter presents a fully automatic method to segment pelvic and proximal femoral surfaces from hip CT images. A coarse-to-fine strategy was proposed to combine multi-atlas segmentation with graph-based surface detection. The multi-atlas segmentation step seeks to coarsely extract the entire hip joint region. It uses automatically detected anatomical landmarks to initialize and select the atlas and accelerate the segmentation. The graph based surface detection is to refine the coarsely segmented hip joint region. It aims at completely and efficiently separate the adjacent surfaces of the acetabulum and the femoral head while preserving the hip joint structure. The proposed strategy was evaluated on 30 hip CT images and provided an average accuracy of 0.55, 0.54, and 0.50 mm for segmenting the pelvis, the left and right proximal femurs, respectively.