Non-linear methods in statistical shape analysis have become increasingly important in orthopedic research as they allow for more accurate and robust analysis of complex shape data such as articulated joints, bony defects and cartilage loss. These methods involve the use of non-linear transformations to describe shapes, rather than the traditional linear approaches, and have been shown to improve the precision and sensitivity of shape analysis in a variety of applications. In orthopedic research, non-linear methods have been used to study a range of topics, including the analysis of bone shape and structure in relation to osteoarthritis, the assessment of joint deformities and their impact on joint function, and the prediction of patient outcomes following surgical interventions. Overall, the use of non-linear methods in statistical shape analysis has the potential to advance our understanding of the relationship between shape and function in the musculoskeletal system and improve the diagnosis and treatment of orthopedic conditions

Following anterior cruciate ligament reconstruction (ACLR) using a semitendinosus (ST) autograft measures such as length, cross-sectional area, and volume may not fully describe the effects of tendon harvest on muscle morphology as these discrete measures cannot characterize three-dimensional muscle shape. This study aimed to determine between-limb ST shape similarity and regional morphology in individuals with a unilateral history of ACLR using a ST graft, and healthy controls. A secondary analysis of magnetic resonance imaging was undertaken from 18 individuals with unilateral history of ST ACLR and 18 healthy controls. ST muscles were manually segmented, and shape similarity were assessed between limbs and groups using Jaccard index (0-1) and Hausdorff distance (mm). ST length (cm), peak cross-sectional area (CSA) (cm. 2. ), and volume (cm. 3. ) was compared between surgically reconstructed and uninjured contralateral limbs, and between the left and right limbs of control participants with no history of injury. Cohen's d was reported as a measure of effect size. Compared to healthy controls, the ACLR group had significantly (p<0.001, d= −2.33) lower bilateral ST shape similarity. Furthermore, the deviation in muscle shape was significantly (p<0.001, d= 2.12) greater in the ACLR group. Within the ACLR group, maximum Hausdorff distance indicated ST from the ACLR limb deviated (23.1±8.68 mm) from the shape of the healthy contralateral ST, this was observed particularly within the distal region of the muscle. Compared to the uninjured contralateral limb and healthy controls, deficits in peak cross-sectional area and volume in ACLR group were largest in proximal (p<0.001, d= −2.52 to −1.28) and middle (p<0.001, d= −1.81 to −1.04) regions. Findings highlight morphological features in distal ST not identified by traditional discrete morphology measures. ST shape was most different in the distal region of the muscle, despite deficits in CSA and volume being most pronounced in proximal and middle regions. ST shape following ACLR may affect force transmission and distribution within the hamstrings and contribute to persistent deficits in knee flexor and internal rotator strength

Introduction. Kienböck's disease is generally defined as the collapse of the lunate bone, and this may lead to early wrist osteoarthritis. Replacing the collapsed lunate with an implant has regained renewed interest with the advancing technology of additive manufacturing, enabling the design of patient-specific implants. The aims of this project are (1) to determine how accurate it is to use the contralateral lunate shape as a template for patient-specific lunate implants, and (2) to study the effects of shape variations wrist kinematics using 4D-computed tomography (CT) scanning. Methods. A 3D statistical shape model (SSM) of the lunate was built based on bilateral CT scans of 54 individuals. Using SMM, shape variations of the lunate were identified and the intra- and inter-subject shape variations were compared by performing an intraclass correlation analysis. A radiolucent motor-controlled wrist-holder was designed to guide flexion/extension and radial/ulnar deviation of ex vivo wrist specimens under 4D-CT scanning. In this pilot, three shape mode variations were tested per specimen in two specimens were. After post-processing each CT, the scapholunate angle (SLA) and capitolunate angle (CLA) were measured. Results. The shape of the lunate was not symmetrical, defined as exceeding the intra-subject variation in five different shape modes. The FE tests show a generalized increase in scapholunate and capitolunate angle when using lunate implants, and comparing variation of shape modes showed that shape mode 3 has a significant effect on the measured angles (p<0.05). Discussion. The design of patient-specific lunate implants may prove to be challenging using a ‘mirror’-design as it will lead to a degree of shape asymmetry. The pilot study, to determine the effects of those shape variations on wrist kinematics suggest that the degree of shape variation observed indeed may alter the wrist kinematics, although this needs to be further investigated in study using more specimens

The anatomy of the femur shows a high inter-patient variability, making it challenging to design standard prosthetic devices that perfectly adapt to the geometry of each individual. Over the past decade, Statistical Shape Models (SSMs) have been largely used as a tool to represent an average shape of many three-dimensional objects, as well as their variation in shape. However, no studies of the morphology of the residual femoral canal in patients who have undergone an amputation have been performed. The aim of this study was therefore to evaluate the main modes of variation in the shape of the canal, therefore simulating and analysing different levels of osteotomy. To assess the variability of the femoral canal, 72 CT-scans of the lower limb were selected. A segmentation was performed to isolate the region of interest (ROI), ranging from the lesser tip of the trochanter to the 75% of the length of the femur. The canals were then sized to scale, aligned, and 16 osteotomy levels were simulated, starting from a section corresponding to 25% of the ROI and up to the distal section. For each level, the main modes of variations of the femoral canal were identified through Principal Component Analysis (PCA), thus generating the mean geometry and the extreme shapes (±2 stdev) of the principal modes of variation. The shape of the canals obtained from these geometries was reconstructed every 10 mm, best- fitted with an ellipse and the following parameters were evaluated: i) ellipticity, by looking at the difference between axismax and axismin; ii) curvature of the canal, calculating the arc of circumference passing through the shapes’ centroids; iii) conicity, by looking at the maximum/minimum diameter; iv) mean diameter. To understand the association between the main modes and the shape variance, these parameters were compared, for each level of osteotomy, between the two extreme geometries of the main modes of variation. Results from PCA pointed out that the first three PCs explained more than the 87% of the total variance, for each level of simulated osteotomy. By analysing the extreme geometries for a distal osteotomy (e.g. 80% of the length of the canal), the first PC was associated to a combination of ROC (var%=41%), conicity (var%=28%) and ellipticity (var%=7%). PC2 was still associated with the ROC (var%=16%), while PC3 turned out to be associated with the diameter (var%=38%). Through the SSM presented in this study, a quantitatively evaluation of the deformation of the intramedullary canal has been made possible. By analysing the extreme geometries obtained from the first three modes of variance, it is clear that the first three PCs accounted for the variations in terms of curvature, conicity, ellipticity and diameter of the femoral canal with a different weight, depending on the level of osteotomy. Through this work, it was also possible to parametrize these variations according to the level of excision. The results given for the segment corresponding to the 80% of the length of the canal showed that, at that specified level, the ROC, conicity and ellipticity were the anatomical parameters with the highest range of variability, followed by the variation in terms of diameter. Therefore, the analysis carried out can provide information about the relevance of these parameters depending on the level of osteotomy suffered by the amputee. In this way, optimal strategies for the design and/or customization of osteo-integrated stems can be offered depending on the patient's residual limb

For many years, marker-based systems have been used for motion analysis. However, the emergence of new technologies, such as 4D scanners provide exciting new opportunities for motion analysis. In 4D scanners, the subjects are measured as a dense mesh, which enables the use of shape analysis techniques. In this talk, we will explore how the combination of the rising new motion analysis methods and shape modelling may change the way we think about movement and its analysis

During shoulder arthroplasty the native functionality of the diseased shoulder joint is restored, this functionality is strongly dependent upon the native anatomy of the pre-diseased shoulder joint. Therefore, surgeons often use the healthy contralateral scapula to plan the surgery, however in bilateral diseases such as osteoarthritis this is not always feasible. Virtual reconstructions are then used to reconstruct the pre-diseased anatomy and plan surgery or subject-specific implants. In this project, we develop and validate a statistical shape modeling method to reconstruct the pre-diseased anatomy of eroded scapulae with the aim to investigate the existence of predisposing anatomy for certain shoulder conditions. The training dataset for the statistical shape model consisted of 110 CT images from patients without observable scapulae pathologies as judged by an experienced shoulder surgeon. 3D scapulae models were constructed from the segmented images. An open-source non-rigid B-spline-based registration algorithm was used to obtain point-to-point correspondences between the models. The statistical shape model was then constructed from the dataset using principle component analysis. The cross-validation was performed similarly to the procedure described by Plessers et al. Virtual defects were created on each of the training set models, which closely resemble the morphology of glenoid defects according to the Wallace classification method. The statistical shape model was reconstructed using the leave-one-out method, so the corresponding training set model is no longer incorporated in the shape model. Scapula reconstruction was performed using a Monte Carlo Markov chain algorithm, random walk proposals included both shape and pose parameters, the closest fitting proposal was selected for the virtual reconstruction. Automatic 3D measurements were performed on both the training and reconstructed 3D models, including glenoid version, critical shoulder angle, glenoid offset and glenoid center position. The root-mean-square error between the measurements of the training data and reconstructed models was calculated for the different severities of glenoid defects. For the least severe defect, the mean error on the inclination, version and critical shoulder angle (°) was 2.22 (± 1.60 SD), 2.59 (± 1.86 SD) and 1.92 (± 1.44 SD) respectively. The reconstructed models predicted the native glenoid offset and centre position (mm) an accuracy of 0.87 (± 0.96 SD) and 0.88 (± 0.57 SD) respectively. The overall reconstruction error was 0.71 mm for the reconstructed part. For larger defects each error measurement increased significantly. A virtual reconstruction methodology was developed which can predict glenoid parameters with high accuracy. This tool will be used in the planning of shoulder surgeries and investigation of predisposing scapular morphologies

Abstract. Objectives. The objective of this study is to investigate the effect of solvents and rheological properties of PCL/Hydroxyapatite ink on the shape fidelity of the 3D printed scaffolds for bone tissue engineering. Methods. A series of inks were made consisting of 50% (w/v) of polycaprolactone (PCL) filled with 0%, 3.5% and 12.5% (w/V) of hydroxyapatite (HA) in dichloromethane (DCM) and chloroform (CHF). Steady and oscillatory shear rheological tests were performed on a rheometer (Discovery HR-3). Solvent-cast direct ink writing was performed with a custom-made 3D printer for the fabrication of PCL/HA scaffold structures with 2–8 layers. Optical microscope and scanning electron microscopy (SEM) were used to assess the shape fidelity. Results. Shape fidelity of the inks was quantitively assessed on the 3D printed scaffold structures allowing subjective comparisons. The addition of HA particles increased zero-shear viscosity by up to 900%. For oscillatory tests, plateau of storage modulus was observed in the low-frequency region which is attributed to good dispersion of the HA particles inside the matrix that leads to the formation of filler networks, resulting in pseudo-solid behavior and shape fidelity improvement. As the HA concentration increases, the plateau becomes more pronounced and the shape fidelity increases. With the same concentration, all DCM inks also show higher viscosity (from 10% to 200%) and better shape fidelity than CHF inks. As DCM has a lower boiling point (39.6 °C) than CHF (61.2°C), DCM evaporates quicker reducing the fusion and diffusion of deposited ink filaments before solidification which is observed in SEM images. Conclusions. This study reveals insights into using rheological characterizations as a tool for evaluation of shape fidelity of solvent-based DIW inks and also provides fundamental information on the influence of different solvents on the fidelity of 3D printed scaffolds. Declaration of Interest. (b) declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported:I declare that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research project

Pre-operative 3D glenoid planning improves component placement in terms of version, inclination, offset and orientation. Version and inclination measurements require the position of the inferior angle. As a consequence, current planning tools require a 3D model of the full scapula to accurately determine the glenoid parameters. Statistical shape models (SSMs) can be used to reconstruct the missing anatomy of bones. Therefore, the objective of this study is to develop and validate an SSM for the reconstruction of the inferior scapula, hereby reducing the irradiation exposure for patients. The training dataset for the statistical shape consisted of 110 CT images from patients without observable scapulae pathologies as judged by an experienced shoulder surgeon. 3D scapulae models were constructed from the segmented images. An open-source non-rigid B-spline-based registration algorithm was used to obtain point-to-point correspondences between the models. A statistical shape model was then constructed from the dataset using principal component analysis. Leave-one-out cross-validation was performed to evaluate the accuracy of the predicted glenoid parameters from virtual partial scans. Five types of virtual partial scans were created on each of the training set models, where an increasing amount of scapular body was removed to mimic a partial CT scan. The statistical shape model was reconstructed using the leave-one-out method, so the corresponding training set model is no longer incorporated in the shape model. Reconstruction was performed using a Monte Carlo Markov chain algorithm, random walk proposals included both shape and pose parameters, the closest fitting proposal was selected for the virtual reconstruction. Automatic 3D measurements were performed on both the training and reconstructed 3D models, including glenoid version, inclination, glenoid centre point position and glenoid offset. In terms of inclination and version we found a mean absolute difference between the complete model and the different virtual partial scan models of 0.5° (SD 0.4°). The maximum difference between models was 3° for inclination and 2° for version. For offset and centre point position the mean absolute difference was 0 mm with an absolute maximum of 1 mm. The magnitude of the mean and maximum differences for all anatomic measurements between the partial scan and complete models is smaller than the current surgical accuracy. Considering these findings, we believe a SSM based reconstruction technique can be used to accurately reconstruct the glenoid parameters from partial CT scans

Inter-subject variability is inherently present in patient anatomy and is apparent in differences in shape, size and relative alignment of the bony structures. Understanding the variability in patient anatomy is useful for distinguishing between pathologies and to assist in surgical planning. With the aim of supporting the development of stratified orthopaedic interventions, this work introduces an Articulated Statistical Shape Model (ASSM) of the lower limb. The model captures inter-subject variability and allows reconstructing ‘virtual’ knee joints of the lower limb shape while considering pose. A training dataset consisting of 173 lower limbs from CT scans of 110 subjects (77 male, 33 female) was used to construct the ASSM of the lower limb. Each bone of the lower limb was segmented using ScanIP (Simpleware Ltd., UK), reconstructed into 3D surface meshes, and a SSM of each bone was created. A series of sizing and positioning procedures were carried out to ensure all the lower limbs were in full extension, had the same femoral length and that the femora were aligned with a coincident centre. All articulated lower limbs were represented as: (femur scale factor) × (full extension articulated lower limb + relative transformation of tibia, fibula and patella to femur). Articulated lower limbs were in full extension were used to construct a statistical shape model, representing the variance of lower limb morphology. Relative transformations of the tibia, fibula and patella versus the femur were used to form a statistical pose model. Principal component analysis (PCA) was used to extract the modes of changes in the model. The first 30 modes of the shape model covered 90% of the variance in shape and the first 10 modes of the pose model covered 90% of the pose variance. The first mode captures changes of the femoral CCD angle and the varus/valgus alignment of the knee. The second mode represents the changes in the ratio of femur to tibia length. The third mode reflects change of femoral shaft diameter and patella size. The first mode characterising pose captures the medial/lateral translation between femur and tibia. The second mode represents variation in knee flexion. The third mode reflects variation in tibio-femoral joint space. An articulated statistical modelling approach was developed to characterize inter-subject variability in lower limb morphology for a set of training specimens. This model can generate large sets of lower limbs to systematically study the effect of anatomical variability on joint replacement performance. Moreover, if a series of images of the lower limb during a dynamic activity are used as training data, this method can be applied to analyse variance of lower limb motion across a population

As stem cells and primary cells hold potential for improving disease outcomes and patient lives, methods for steering cell fate are of considerable importance. In this context, an emerging method is directing cell function through controlling cellular shape. The talk will discuss how cell functions are based on mechano-transduction events related to the balance of intra- and extracellular forces. The talk will explore the multiple biophysical cues that affect cell shape and present methods for directly generating cell shape, e. g. micro-contact printing used for directing the differentiation lineage of stem cells. Based on our own work, the talk will introduce the novel concept that specific biomaterial types and stiffnesses can be chosen for generating specific cellular “baseline shapes” and associated function. As our cells are exposed to continuously changing biomechanical forces, the talk will also report how specific forces can be used for engineering shape. The talk will explore how biomaterial stiffness and biomechanical forces act together on cellular shape, and whether one of the two stimuli is able to override the other. The novel insights reported here are fundamental for designing cell shape-instructive 3D biomaterials in the context of steering cell function in situ for regenerative medicine

Summary. Our statistical shape analysis showed that size is the primary geometrical variation factor in the medial meniscus. Shape variations are primarily focused in the posterior horn, suggesting that these variations could influence cartilage contact pressures. Introduction. Variations in meniscal geometry are known to influence stresses and strains inside the meniscus and the articulating cartilage surfaces. This geometry-dependent functioning emphasizes that understanding the natural variation in meniscus geometry is essential for a correct selection of allograft menisci and even more crucial for the definition of different sizes for synthetic meniscal implants. Moreover, the design of such implants requires a description of 3D meniscus geometry. Therefore, the aim of this study was to quantify 3D meniscus geometry and to determine whether variation in medial meniscus geometry is size or shape driven. Patients & Methods. Sagittal knee MR images (n=35; 15 males, 20 females, aged 33±12) were acquired at 3 Tesla using a 3D SPACE sequence with isotropic resolution of 0.5×0.5×0.5mm. 3D models were generated by manual segmentation of the medial menisci from the MR scans. The surface of a reference meniscus was then described by 250 landmarks. Using an affine iterative closest point transformation, these landmarks were registered onto the full set of 3D models. Based on the set of corresponding landmarks, a point distribution model was created using the Shapeworks software (NITRC, University of Utah), an open source algorithm for constructing correspondence-based statistical models of sets of similar shapes. Several modules from Shapeworks and the Arthron software (UCD, Dublin) were used to perform principal component analysis (PCA) upon the set of landmarks. The results of the PCA enabled quantification and visualisation of the primary modes of variation in meniscal geometry. Results. The majority (77%) of variation in medial meniscus geometry was found to be due to sizing (principal component (PC) 1). Including the shape-related PC's 2 to 4, increased the cumulative percentage of represented geometry variation to over 90%. The independent shape variations described by PCs 2–4 all display larger variations in geometry of the posterior meniscal horn than the anterior section. Discussion. From this study, we can conclude that geometry variation of the medial meniscus is mainly determined by differences in size. However, since the posterior aspect of the medial meniscus experiences higher loads during daily activities than the anterior part, the shape variations described by PCs 2–4 may have a significant influence on cartilage contact pressures. Therefore, PCA alone does not provide sufficient information to define the number of implant sizes to cover a majority of the population. Analysis of the sensitivity of cartilage contact pressures to the shape variations identified in this analysis could provide the additional information needed

In uncemented total hip arthroplasty (THA), the optimal femoral component should allow both maximum cortical contact with proximal load transfer and accurate restoration of individual joint biomechanics. This is often compromised due to a high variability in proximal femoral anatomy. The aim of this on-going study is to assess the variation in proximal femoral canal shape and its association with geometric and anthropometric parameters in primary hip OA. In a retrospective cohort study, AP-pelvis radiographs of 98 consecutive patients (42 males, 56 females, mean age 61 (range:45-74) years, BMI 27.4 (range:20.3-44.6) kg/m2) who underwent THA for primary hip OA were reviewed. All radiographs were calibrated and femoral offset (FO) and neck-shaft-angle (NSA) were measured using a validated custom programme. Point-based active shape modelling (ASM) was performed to assess the shape of the inner cortex of the proximal femoral meta- and diaphysis. Independent shape modes were identified using principal component analysis (PCA). Hierarchical cluster analysis of the shape modes was performed to identify natural groupings of patients. Differences in geometric measures of the proximal femur (FO, NSA) and demographic parameters (age, height, weight, BMI) between the clusters were evaluated using Kruskal-Wallis one-way-ANOVA or Chi-square tests, as appropriate. In the entire cohort, mean FO was 39.0 mm, mean NSA was 131 degrees. PCA identified 10 independent shape modes accounting for over 90% of variation in proximal femoral canal shape within the dataset. Cluster Analysis revealed 6 shape clusters for which all 10 shape modes demonstrated a significantly different distribution (p-range:0.000-0.015). We observed significant differences in age (p=0.032), FO (p<0.001) and NSA (p<0.001) between the clusters. No significant differences with regard to gender or BMI were seen. Our preliminary analysis has identified 6 different patterns of proximal femoral canal shape which are associated with significant differences in femoral offset, neck-shaft-angle and age at time of surgery. We are currently evaluating the entire dataset of 345 patients which will allow a comprehensive classification of variation in proximal femoral shape and joint geometry. The present data may optimise preoperative planning and improve future implant design in THA

Summary Statement. A large proportion of knee arthroplasty patients are dissatisfied with their replacement. Significant differences exist between preoperative, postoperative and normal kinematics. A better understanding of the inter-relationships between kinematics, shape and prosthesis placement could lead to improved quality of life. Introduction. Knee kinematics are altered by total knee arthroplasty (TKA) both intentionally and unintentionally. Knowledge of how and why kinematics change may improve patient outcome and satisfaction through improved implant design, implant placement or rehabilitation. Comparing preoperative to postoperative kinematics and shape of the natural and replaced joint will allow an investigation of the inter-relationships between knee shape, prosthesis placement, knee kinematics and quality of life. Patients & Methods. Using a sequential-biplanar radiographic protocol that allowed imaging the preoperative and postoperative patellofemoral (PF) and tibiofemoral (TF) joints under weightbearing throughout the range of motion, we imaged and compared the 6 degree-of-freedom PF and TF kinematics of 9 pre-TKA subjects to those of 15 post-TKA subjects (Zimmer NexGen Legacy Posterior Stabilised Gender Solutions (GS) components). Using a novel computed tomography (CT) protocol, we obtained the femoral, tibial and patellar knee shapes, plus component placement after TKA. The same 9 pre-TKA subjects have now been re-imaged a minimum of one year postoperatively (DePuy Sigma Mobile Bearing cruciate-sacrificing components) to determine their changes in knee geometry and kinematics; full analysis is in progress. Results. Clear, statistically significant differences were seen between the kinematics of the pre-TKA and post-TKA groups. For the TF joint, the tibia was more posterior and inferior in the post-TKA group compared to the pre-TKA group (max 20 mm and 15 mm, respectively) (p<0.001). Subjects had neutral alignment in the post-TKA group compared to varus alignment (max 9°) in the pre-TKA group (p<0.001). For the PF joint, the patella was shifted more posteriorly and less laterally postoperatively and was tilted neutrally compared to laterally (p<0.001). Our preliminary analysis of the matched preop-postop subjects likewise shows a more posterior and inferior tibia and neutral versus valgus alignment. Greater tibial rotations were seen postoperatively due to the mobile bearing. The patella was more posterior and less lateral postoperatively, as seen with the two groups. Discussion/Conclusion. The kinematic differences seen are likely due to a combination of surgical, implant and patient factors. Both groups showed differences from normal kinematics, based on previous studies in the literature. In the future, by comparing the preoperative and postoperative kinematics, shape and quality of life for the same subjects (i.e. the 9 pre-TKA subjects in this study), and analyzing the interrelationships amongst these, we aim to determine if a different implant shape or different component positioning could create more normal kinematics, resulting in a better clinical outcome

The shape of the vertebral bodies from L1 to L4 was assessed from lateral dual-energy x-ray absorptiometry (DXA) images using an active shape model. The output from the model was compared to measurements of areal bone mineral density in L1 to L4 (aBMD) using a stepwise linear regression model. A significant relationship was found between aBMD and vertebral shape that suggests that the method may be useful for correcting artefacts such as osteophytes

Summary. This work proposes a novel, automatic method to obtain an anatomical reconstruction for 3D segmented bones with large acetabular defects. The method works through the fitting of a Statistical Shape Model to the non-defect parts of the bone. Introduction. Patient-specific implants can be used to treat patients with large acetabular bone defects (IIa-c, IIIb, Paprosky 1994). These implants require a full 3D preoperative planning that includes segmentation of volumetric images (CT or MRI), extraction of the 3D shape, reconstruction of the bone defect into its anatomic (non-defect) state, design of an implant with a perfect fit and optimal placement of the screws. The anatomic reconstruction of the bone defect will play a key role in diagnosing the amount of bone loss and in the design of the implant. Previous reconstruction methods rely on a healthy contralateral (Gelaude 2007); however this is not always available (e.g. partial scan or implant present). Statistical shape models (SSM) of healthy bones can help to increase the accuracy and usability of the reconstruction and will decrease the manual labor and user dependency. Skadlubowicz (2009) illustrated the use of an SSM to reconstruct pelvic bones with tumor defects; however tumors generally affect a smaller region of the bone so that the reconstruction will be easier than in large acetabular bone defects. Also, the tumor reconstruction method uses 80 manually indicated landmarks, while the proposed method only uses 14. Patients & Methods. CT-scans from subjects with a healthy hemi-pelvis (15 male, 33 female, mean age: 69±20) were used to generate an SSM. The CT-scans were segmented using Mimics (Materialise NV, Belgium) to create a triangulated mesh. Preprocessing of the meshes ensured that the triangulation was smooth and uniform to help solve the corresponding point problem. An algorithm based on Redert (1999) was used to morph the template hemi-pelvis onto each dataset entity, creating a dataset with corresponding points. From this dataset the SSM was calculated using principal component analysis, so that the principal components serve as parameters for the mathematical model of the hemi-pelvis. To fit the SSM to a new defect hemi-pelvis, a matching algorithm was used. The algorithm varies the Principal Components independently optimizing the distance of the non-defect parts of the defect hemi-pelvis to the SSM sample. To validate the reconstruction method, 6 healthy bone meshes were used to generate a synthetic defect in the acetabular region. The original mesh was used as ‘golden standard’ to measure the reconstruction error. To illustrate the clinical use of the reconstruction method, one hemi-pelvis with a substantial defect was reconstructed. Results. The correspondence error for the morphing algorithm was 4.68±0.78 mm. The leave-one-out error for the SSM was 1.30±0.96 mm. The reconstruction error for the non-defect part was 1.44±1.13mm and for the reconstructed part 2.15±1.53mm. Discussion/Conclusion. The proposed method performs comparable to the contralateral method and the tumor reconstruction method, without the need of a healthy contralateral geometry. Consequently, the validation and the clinical illustration show that the proposed method is promising for automatic reconstruction of large acetabular defects

Background. The cartilaginous growth plate (GP) is a zonal structure, in which chondrocytes are organized into columns and drive the longitudinal elongation of the endochondral bones. In the proliferative zone (PZ), cells exhibit high mitotic activity, are flattened and oriented along the mediolateral (ML) axis of the GP. Mitotic figures in the elongated chondrocytes lie perpendicular to the proximo-distal (PD) direction of growth, while cell divisions occur parallel to the columns followed by a gliding movement of the daughter cells. The mechanisms responsible for the geometrical anisotropy and columnar arrangement of PZ chondrocytes are poorly understood. Here, we assessed the function of the adhesive receptor β1 integrins on spindle and division geometry in chondrocytes using mouse genetics. Methods. GP slices were prepared from wild type (wt) and β1fl/fl-Col2a1cre mice. Primary rib chondrocytes were cultured on substrate-coated glass slides and fluorescently stained with anti-alpha-tubulin and anti-pericentrin antibodies, with phalloidin and DAPI. Confocal stacks were built from images acquired by confocal microscopy. Cell and spindle orientation relative to the PD axis (in vivo) or to the substrate plane (SP) in vitro were determined by geometric functions. Shape indexes (SI) were calculated as the ratio of the length and height of the cell. Results. During GP morphogenesis, aggregating mesenchymal cells (E11) were polygonal (SI=1.43) and nonoriented. Upon chondrogenic differentiation at E12, wt GP chondrocytes showed moderate flattening (SI=1.93) and tend to align perpendicular to the PD axis. At E13, PZ chondrocytes further flattened (SI=3.41) and largely organized into lines along the ML axis. At E14, the first vertical stacks formed, which were gradually elongated along the PD axis at later stages and composed of extremely flat (SI=4.91), ML-oriented chondrocytes. Early spindles were randomly oriented at all stages, whereas late spindles were aligned along the long axis of the cell. In contrast, β1 integrin null PZ chondrocytes were rounded (SI=1.37), displayed random orientation of both early and late spindles, and failed to organize into columns. On collagen II, both wt and β1-null primary chondrocytes remained rounded and displayed random spindle orientation relative to the SP. On fibronectin (FN), wt chondrocytes were flattened (SI=3.26) and their spindle was aligned parallel to SP. Mutant cells lacking the major FN receptor α5β1 integrin did not spread, were rounded (SI=1.73) and aligned their spindle randomly. On vitronectin, both cell types spread and flattened well (SI=4.11 and SI=4.20) accompanied by parallel spindle orientation to the SP. Conclusion. We propose that β1 integrins provide adhesive cues for chondrocyte geometry and orientation, which in turn orients the mitotic spindle and determines the division axis; and for chondrocyte intercalation to shape the three-dimensional structure of the GP

We report a study of the shapes of the tibial and femoral articular surfaces in sagittal, frontal and coronal planes which was performed on cadaver knees using two techniques, MRI and computer interpolation of sections of the articular surfaces acquired by a three-dimensional digitiser. The findings using MRI, confirmed in a previous study by dissection, were the same as those using the digitiser. Thus both methods appear to be valid anatomical tools. The tibial and femoral articular surfaces can be divided into anterior segments, contacting from 0° to 20 ± 10° of flexion, and posterior segments, contacting from 20 ± 10° to 120° of flexion. The medial and lateral compartments are asymmetrical, particularly anteriorly. Posteromedially, the femur is spherical and is located in a conforming, but partly deficient, tibial socket. Posterolaterally, it is circular only in the sagittal section and the tibia is flat centrally, sloping downwards both anteriorly and posteriorly to receive the meniscal horns. Anteromedially, the femur is convex with a sagittal radius larger than that posteriorly, while the tibia is flat sloping upwards and forwards. Anterolaterally, both the femoral and tibial surfaces are largely deficient. These shapes suggest that medially the femur can rotate on the tibia through three axes intersecting in the middle of the femoral sphere, but that the sphere can only translate anteroposteriorly and even then to a limited extent. Laterally, the femur can freely translate anteroposteriorly, but can only rotate around a transverse axis for that part of the arc, i.e., near extension, during which it comes into contact with the tibia through its flattened distal/medial surface as against its spherical posterior surface

Introduction. The proteoglycan aggrecan is a major component of the cartilaginous matrices which provides resistance against compressive forces. Spontaneously occurring functional null mutations in the aggrecan gene (Acan) in various species lead to perinatal chondrodysplasia. The aim of the present study was to investigate the cellular and biomechanical properties of the cartilaginous growth plate, and the development of intervertebral disc in a novel, experimentally induced aggrecan mutant mouse strain carrying an insertion in exon 5 of the Acan gene. Methods. The novel aggrecan mutant mice were generated by inserting a loxP site into exon 5 (E5i) by homologous recombination in ES cells. Wild type and homozygous mutant (Acan-E5i/E5i) mice were analyzed by skeletal staining, histology and immunohistochemistry. Proliferation and survival were assessed by phosphorylated histone H3 immunostaining and TUNEL assay, respectively. Shape index (SI) in the proliferative zone (PZ) of the growth plate (GP) was calculated as a ratio of the long and short axes of the cells. Orientation of the PZ chondrocytes was characterized by the angle between the cell long axis and longitudinal direction of the bone growth. Imaging and stiffness measurements were performed by atomic force microscopy (AFM). Results. Acan-E5i/E5i mice are characterized by severe dwarfism, short snout, protruding tongue, cleft palate, and die at birth due to respiratory failure. On sections the cartilage of mutant mice appeared as tightly packed chondrocytes surrounded by a compressed matrix. At E18.5 and E14.5, the mutant PZ consisted of rounded (SI=1.71 at E18.5; SI=1.72 at E14.5) non-oriented chondrocytes, compared to the wild type PZ with flattened (SI=3.92 at E18.5; SI=3.90 at E14.5), columnar cells oriented with right angle to the longitudinal axis of the growth. At E13, the shape and orientation of mutant chondrocytes were similar to control. AFM at E14.5 and E18.5 demonstrated a stiffer matrix with denser collagen network in the mutant compared to wild type. The mutant cartilage had increased apoptosis and reduced proliferation rate at E18.5. The IVDs development appeared normal at E13.5-E14.5, however, the IVD was severely malformed at E18.5. Discussion. We have shown that aggrecan deficiency impairs cartilage biomechanics and results in a stiffer matrix. The altered mechanical properties might be responsible for the disorganization of mutant GP and compression of the IVD at around birth. Interestingly, the altered matrix mechanics is dispensable for early flattening and orientation of GP proliferative chondrocytes. In summary, aggrecan is essential for proper cartilage cytoarchitecture and morphogenesis by ensuring the suitable mechanical environment

Introduction. Osteoporotic vertebral fractures can cause severe vertebral wedging and kyphotic deformity. This study tested the hypothesis that kyphoplasty restores vertebral height, shape and mechanical function to a greater extent than vertebroplasty following severe wedge fractures. Methods. Pairs of thoracolumbar “motion segments” from seventeen cadavers (70–97 yrs) were compressed to failure in moderate flexion and then cyclically loaded to create severe wedge deformity. One of each pair underwent vertebroplasty and the other kyphoplasty. Specimens were then creep loaded at 1.0kN for 1 hour. At each stage of the experiment the following parameters were measured: vertebral height and wedge angle from radiographs, motion segment compressive stiffness, and stress distributions within the intervertebral discs. The latter indicated intra-discal pressure (IDP) and neural arch load-bearing (F. N. ). Results. Fracture and cyclic loading reduced anterior vertebral height by 34%, increased wedge angle from 5.0° to 11.4°, increased F. N. by 58% and reduced IDP and compressive stiffness by 96% and 44% respectively. Kyphoplasty restored anterior height to a greater extent than vertebroplasty (p<0.001), by 96% versus 59% immediately after augmentation, and by 79% versus 47% after subsequent creep loading. Wedge angle was also reduced to a greater extent following kyphoplasty than vertebroplasty (p<0.02) by 7.2° vs 4.2° after augmentation and 6.6° vs 4.0° after creep loading. IDP, F. N. and compressive stiffness were restored to a similar extent by both procedures. Conclusion. Kyphoplasty and vertebroplasty were equally effective in restoring mechanical function following severe wedge fractures, but kyphoplasty was better able to correct deformity by restoring vertebral height and reducing wedging. No conflicts of interest. Sources of funding: Funding was provided by a Royal College of Surgeons of England Research Fellowship and the Gloucestershire Arthritis Trust. Materials were provided by Medtronic and Depuy. This abstract has not been previously published in whole or substantial part nor has been presented previously at a national meeting

In six unloaded cadaver knees we used MRI to determine the shapes of the articular surfaces and their relative movements. These were confirmed by dissection. Medially, the femoral condyle in sagittal section is composed of the arcs of two circles and that of the tibia of two angled flats. The anterior facets articulate in extension. At about 20° the femur ‘rocks’ to articulate through the posterior facets. The medial femoral condyle does not move anteroposteriorly with flexion to 110°. Laterally, the femoral condyle is composed entirely, or almost entirely, of a single circular facet similar in radius and arc to the posterior medial facet. The tibia is roughly flat. The femur tends to roll backwards with flexion. The combination during flexion of no antero-posterior movement medially (i.e., sliding) and backward rolling (combined with sliding) laterally equates to internal rotation of the tibia around a medial axis with flexion. About 5° of this rotation may be obligatory from 0° to 10° flexion; thereafter little rotation occurs to at least 45°. Total rotation at 110° is about 20°, most if not all of which can be suppressed by applying external rotation to the tibia at 90°