Summary Statement. Subject specific FE models of human Achilles tendon were developed and optimum material properties were found. Stress concentration occurred at the midsection but dependent on stiffening and thinning of tendon, indicating that they are two major factors for tendon rupture. Introduction. Achilles tendon injuries are common, occurring about 250,000 per year in the US alone, yet the mechanisms of tendinopathy and rupture remain unknown. Most Achilles tendon ruptures occur at 2 to 6 cm above the insertion to the calcaneus bone. Previous angiographic studies have suggested that there is an avascular area in this region. However, it is not understood why that region receives poor blood supply and prone to rupture. The aim of this study is to investigate influence of geometry and material properties on Achilles tendon rupture with mechanical experiment and corresponding subject-specific finite element (FE) analysis. Patients & Methods. Mechanical experiment was performed on 10 fresh human Achilles tendons. High frequency ultrasound images were used to measure cross sectional areas at the midsection of the tendon. Cyclic testing was performed to measure mechanical properties and failure loads. Subject-specific FE models of these tendons were generated with Free Form Deformation (FFD) technique. FE mechanical simulations that mimic the experimental cyclic loading were performed on these subject specific models. Tendon material properties were described as transversely isotropic hyperelastic and the optimum material parameters for the human Achilles tendon were obtained. Linear portion of the cyclic loading data was used as boundary conditions. Measured strains from the experiment were compared with predicted strains from the FE analysis. This process was repeated until optimum parameters were found. The influence of geometry and material properties on the Achilles tendon rupture was then investigated– first with subject-specific geometry with average material properties and then with subject-specific material properties with average geometry. Results. Our results indicate that a significant variation exist in the geometry and material properties in human Achilles tendons. Stress concentrations occurred at the midsection of the tendon, supporting previous studies that reported tendon rupture at the region. In particular the thinning of midsection in geometry is highly correlated with the collagen uncrimpping rate in material properties where thinner midsection leads to faster uncrimpping of collagen fibres. Variations in geometry led to shifts in the location of stress concentration within the midsection while variations in material property led the change in the magnitude of stress concentration. Discussion/Conclusion. Our results indicate that Achilles tendon rupture is highly dependent on subject-specific geometry and material properties. In particular the mid section is the location of stress concentration but depending on the geometrical shape, multiple stress concentrations occur, making the tendon more prone to rupture while the material properties influenced the magnitude of stress concentration. Our results indicate stiffening and thinning of tendon may lead to higher risk for tendon rupture

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

The study aimed to compare trochlear profiles in recent total knee arthroplasty (TKA) models and to determine whether they feature improvements compared to their predecessors. The hypothesis was that recent TKA models have more anatomic trochlear compartments and would display no signs of trochlear dysplasia. The authors analyzed the geometry of the 6 following TKA models using engineering software: PFC and Attune (DePuy), NexGen and Persona (Zimmer), Noetos and KneeTec (Tornier). The mediolateral trochlear profiles were plotted at various flexion angles (0°, 15°, 30° and 45°) to deduce the sulcus angle. Analysis of sulcus angles reveals general convergence of recent designs towards anatomic values. At 0° of flexion, sulcus angles of recent implant models were between 156.0–157.4°, while those of previous generation models between 154.5–165.5°. At 30° of flexion, sulcus angles of recent models also lie within 145.7–148.6°, but those of previous models are between 149.5–152.0°. All three manufacturers deepened their trochlear profile at 30° of flexion in recent models compared to earlier designs. Sulcus angles converge towards anatomic values but still exceed radiologic signs of dysplasia by 2–5°. Recent TKA designs have more anatomic trochlear geometries than earlier TKA models by the same manufacturers, but trochlear compartments still exceed radiologic signs of trochlear dysplasia by 2° to 5°. The hypothesis that recent TKA models display no signs of trochlear dysplasia is therefore refuted. Surgeons should be aware of design limitations to optimize choice of implant and extensor mechanisms alignment. Level of evidence: IV geometric implant analysis

We have studied the three-dimensional geometry of the proximal humerus on human cadaver specimens using a digitised measuring device linked to a computer. Our findings demonstrated the variable shape of the proximal humerus as well as its variable dimensions. The articular surface, which is part of a sphere varies individually in its orientation as regards inclination and retroversion, and it has variable medial and posterior offsets. These variations cannot be accommodated by the designs of most contemporary humeral components. Although good clinical results can be achieved with current modular and non-modular components their relatively fixed geometry prevents truly anatomical restoration in many cases. To try to restore the original three-dimensional geometry of the proximal humerus, we have developed a new type of humeral component which is modular and adaptable to the individual anatomy. Such adaptability allows correct positioning of the prosthetic head in relation to an individual anatomical neck, after removal of the marginal osteophytes. The design of this third-generation prosthesis respects the four geometrical variations which have been demonstrated in the present study. These are inclination, retroversion, medial offset and posterior offset

In a study on ten fresh human cadavers we examined the change in the height of the intervertebral disc space, the angle of lordosis and the geometry of the facet joints after insertion of intervertebral total disc replacements. SB III Charité prostheses were inserted at L3-4, L4-5, and L5-S1. The changes studied were measured using computer navigation sofware applied to CT scans before and after instrumentation. After disc replacement the mean lumbar disc height was doubled (p < 0.001). The mean angle of lordosis and the facet joint space increased by a statistically significant extent (p < 0.005 and p = 0.006, respectively). By contrast, the mean facet joint overlap was significantly reduced (p < 0.001). Our study indicates that the increase in the intervertebral disc height after disc replacement changes the geometry at the facet joints. This may have clinical relevance

Trochlear geometry of modern femoral implants is designed for mechanical alignment (MA) technique for TKA. The biomechanical goal is to create a proximalised and more valgus trochlea to better capture the patella and optimize tracking. In contrast, Kinematic alignment (KA) technique for TKA respects the integrity of the soft tissue envelope and therefore aims to restore native articular surfaces, either femoro-tibial or femoro-patellar. Consequently, it is possible that current implant designs are not suitable for restoring patient specific trochlea anatomy when they are implanted using the kinematic technique, this could cause patellar complications, either anterior knee pain, instability or accelerated wear or loosening. The aim of our study is therefore to explore the extent to which native trochlear geometry is restored when the Persona. ®. implant (Zimmer, Warsaw, USA) is kinematically aligned. A retrospective study of a cohort of 15 patients with KA-TKA was performed with the Persona. ®. prosthesis (Zimmer, Warsaw, USA). Preoperative knee MRIs and postoperative knee CTs were segmented to create 3D femoral models. MRI and CT segmentation used Materialise Mimics and Acrobot Modeller software, respectively. Persona. ®. implants were laser scanned to generate 3D implant models. Those implant models have been overlaid on the 3D femoral implant model (generated via segmentation of postoperative CTs) to replicate, in silico, the alignment of the implant on the post-operative bone and to reproduce in the computer models the features of the implant lost due to CT metal artefacts. 3D models generated from post-operative CT and pre-operative MRI were registered to the same coordinate geometry. A custom written planner was used to align the implant, as located on the CT, onto the pre-operative MRI based model. In house software enabled a comparison of trochlea parameters between the native trochlea and the performed prosthetic trochlea. Parameters assessed included 3D trochlear axis and anteroposterior offset from medial facet, central groove, and lateral facet. Sulcus angle at 30% and 40% flexion was also measured. Inter and intra observer measurement variabilities have been assessed. Varus-valgus rotation between the native and prosthetic trochleae was significantly different (p<0.001), with the prosthetic trochlear groove being on average 7.9 degrees more valgus. Medial and lateral facets and trochlear groove were significantly understuffed (3 to 6mm) postoperatively in the proximal two thirds of the trochlear, with greatest understuffing for the lateral facet (p<0.05). The mean medio-lateral translation and internal-external rotation of the groove and the sulcus angle showed no statistical differences, pre and postoperatively. Kinematic alignment of Persona. ®. implants poorly restores native trochlear geometry. Its clinical impact remains to be defined

Summary Statement. Wear of total knee replacement (TKR) is a clinical concern. This study demonstrated low-conformity moderately cross-linked-polyethylene fixed bearing TKRs showed lower volumetric wear than conventional-polyethylene curved fixed bearing TKRs highlighting potential improvement in TKR performance through design and material selection. Introduction. Wear of total knee replacement (TKR) continues to be a significant factor in the clinical performance of the implants. Historically, failure due to delamination and fatigue directed implant design towards more conforming implants to reduce contact stress. However, the new generations of more oxidatively-stable polyethylene have improved the long-term mechanical properties of the material, and therefore allowed more flexibility in the bearing design. The purpose of this study was to investigate the effect of insert conformity and material on the wear performance of a fixed bearing total knee replacement through experimental simulation. Methods. The wear of TKR bearings were investigated using a physiological six station Prosim knee wear simulator (Simulator Solutions, UK). Six samples of each test configuration (Sigma CR fixed bearing knees (DePuy Synthes, UK) were studied, and compared with previously reported data, tested under identical conditions (1, 2). The central axis of the implant was offset from the aligned axes of applied load and tibial rotation to replicate a right knee. High kinematics, under anterior-posterior displacement control was used for this study (3). The lubricant was 25% (v/v) calf serum supplemented with 0.03% (v/v) sodium azide solution in deionised water, as an antibacterial agent, and was changed approximately every 0.33Mc. Wear was assessed gravimetrically and moisture uptake accounted for using unloaded soak controls. Results. The wear rates for the moderately cross-linked inserts (XLK) were significantly lower than the conventional polyethylene (GVF) for all geometries (ANOVA, p<0.05). There was a significant reduction in wear rate as the insert geometry became less conforming for both materials (ANOVA, p<0.05). The wear scars areas were comparable in size and shape between materials, within a geometry group. The size of the wear scar changed with conformity, with the curved inserts showing the largest scars in both anterior-posterior and medial-lateral dimensions, and the flat inserts showing the smallest wear scars. Discussion/Conclusion. The introduction of a moderately cross-linked polyethylene insert was shown to significantly reduce the wear of a fixed bearing total knee replacement compared with a conventional material. There was a trend for reducing wear rate with reducing conformity for both materials, suggesting that reduced conformity results in higher contact pressures and reduced contact area, leading to a reduced surface for wear to occur. Both material and conformity were shown to have a significant impact on the wear of a fixed bearing TKR, and therefore provide opportunity for enhancing wear performance through material and design selection

This study investigates the relationships between Intervertebral Disc (IVD) morphology and biomechanics using patient-specific (PS) finite element (FE) models and poromechanical simulations.

169 3D lumbar IVD shapes from the European project MySpine (FP7-269909), spanning healthy to Pfirrmann grade 4 degeneration, were obtained from MRIs. A Bayesian Coherent Point Drift algorithm aligned meshes to a previously validated structural FE mesh of the IVD. After mesh quality analyses and Hausdorff distance measurements, mechanical simulations were performed: 8 and 16 hours of sleep and daytime, respectively, applying 0.11 and 0.54 MPa of pressure on the upper cartilage endplate (CEP). Simulation results were extracted from the anterior (ATZ) and posterior regions (PTZ) and the center of the nucleus pulposus (CNP). Data mining was performed using Linear Regression, Support Vector Machine, and eXtreme Gradient Boosting techniques. Mechanical variables of interest in DD, such as pore fluid velocity (FLVEL), water content, and swelling pressure, were examined. The morphological variables of the simulated discs were used as input features.

Local morphological variables significantly impacted the local mechanical response. The local disc heights, respectively in the mid (mh), anterior (ah), and posterior (ph) regions, were key factors in general. Additionally, fluid transport, reflected by FLVEL, was greatly influenced (r2 0.69) by the shape of the upper and lower cartilage endplates (CEPs).

This study suggests that disc morphology affects Mechanical variables of interest in DD. Attention should be paid to the antero-posterior distribution and local effects of disc heights. Surprisingly, the CEP morphology remotely affected the fluid transport in NP volumes around mid-height, and mechanobiological implications shall be explored. In conclusion, patient-specific IVD modeling has strong potential to unravel important correlations between IVD phenotypes and local tissue regulation.

Acknowledgments: European Commission: Disc4All-MSCA-2020-ITN-ETN GA: 955735; O-Health-ERC-CoG-2021-101044828

Summary Statement. Repetitive concavities threaded on the surface of bone implants have been already demonstrated to be effective on ectopic bone formation in vivo. The aim of this study was to investigate the effect of concavity on the mineralization process in vitro. Introduction. The role of implant surface geometry in bone formation has been extensively investigated. Ripamonti and co. investigated the possibility to induce bone formation by threading concavities on the surface of calcium phosphate implants, without the need for exogenous osteogenic soluble factors. The underlying hypothesis was that this geometry, by resembling the hemi-osteon trench observable during osteoclastogenesis, was able to activate the ripple-like cascade of bone tissue induction and morphogenesis. Despite several studies indicating a positive effect of concavities on bone induction, so far no attempts have rationalised this phenomenon by means of in vitro tests. Consequently, this study aimed to evaluate the effect of surface concavities on the mineralization of hydroxyapatite (HA) and beta-tricalciumphosphate (b-TCP) ceramics in vitro. Our hypothesis was that concavities could effectively guide the mineralization process in vitro. Materials and Methods. Different-size concavities were threaded into the surface of HA and b-TCP semi-sintered disks: 1.8 (large concavity, LC), 0.8 (medium concavity, MC) and 0.4 mm (small concavity, SC). Disks were fully sintered at 1200 or 1100 °C and soaked in simulated body fluid (SBF) up to 28 days. Distinct experiments were carried out in order to assess the role of chemical composition, sintering temperature and concavity size on the extent of mineralization in vitro. The mineralization process was followed by SEM, EDS, XRD and Ca assay tests. Results and Discussion. Massive mineralization occurred exclusively at the surface of the HA disks treated at 1200 °C. Firstly, aggregates of spherical-like amorphous calcium phosphate nucleated specifically within concavities and not on the planar surface. At a later stage, a flake-like apatitic phase replaced the spherical-like apatite and spread out the concavities until covering the entire disk surface. Instead, a lower extent of mineralization occurred on HA disk treated at 1100 °C, with no formation of flake-like phase. In contrast, no significative mineralization was observed for b-TCP disks, irrespective of sintering temperature. Finally, the extent of mineralization on disks exhibiting SC on the surface was found to be about 120- and 10-fold higher than mineralization disks threaded with LC and MC, respectively. Conclusions. The main results of this study are: i) the in vitro mineralization process of CaP ceramics with surface concavities starts preferentially within the concavities and not on the planar surface; ii) concavity size is an extremely effective parameter for controlling the extent of mineralization in vitro. Finally, the results reported in this study suggest correlation between the positive effect of concavities on mineralization in vitro and on ectopic bone formation in vivo

Summary. Anatomical variations in hip joint anatomy are associated with both the presence and location of tibiofemoral osteoarthritis (OA). Introduction. Variations in hip joint anatomy can alter the moment-generating capacity of the hip abductor muscles, possibly leading to changes in the magnitude and direction of ground reaction force and altered loading at the knee. Through analysis of full-limb anteroposterior radiographs, this study explored the hypothesis that knees with lateral and medial knee OA demonstrate hip geometry that differs from that of control knees without OA. Patients and Methods. This cross-sectional study is an ancillary to the Multicenter Osteoarthritis Study (MOST), an observational cohort study of incident and progressive knee OA in community-dwelling men and women, ages 50–79 years. We report on 160 knees with lateral OA (LOA), 168 knees with medial OA (MOA), and 336 controls. All participants with LOA at the baseline MOST visit were included. An equal number of knees with MOA, and twice the number of control knees were then randomly selected. In participants with bilateral eligibility, a single knee was randomly selected so that all participants contributed only one case or one control knee to the analysis. Case knees were identified as having Kellgren/Lawrence (K/L) ≥ 2 with joint space narrowing (JSN) ≥ 1 in the specified compartment with no JSN in the adjoining compartment. Controls had no radiographic OA (K/L=0 or 1 and JSN=0) in either compartment. Hip joint anatomy parameters were assessed from full-limb standing radiographs using custom OsiriX software by an author (AB) blinded to knee OA status, and unreadable radiographs (N = 8) were discarded prior to unblinding. We measured parameters that influence the abductor moment arm of the hip, including: abductor lever arm, femoral offset, femoral neck length, femoral neck-shaft angle, height of hip centre, body weight lever arm, acetabular version, and abductor angle. All hip measurements were taken from the ipsilateral side of the knee in interest. Each variable was then compared independently among the three groups via analysis of covariance (ANCOVA), controlling for age, sex, and body mass index (BMI), and followed up with a post-hoc Bonferroni analysis to distinguish pairwise group differences. Results. The ANCOVA analysis showed a significant difference in height of hip centre (p=0.001), femoral neck-shaft angle (p=0.009), and abductor angle (p=0.001). Compared to controls, knees with LOA had an increased height of hip centre (p=0.001) and knees with MOA had a decreased abductor angle (p=0.046). Compared to knees with MOA, those with LOA had a more valgus neck-shaft angle (p=0.007) and increased abductor angle (p=0.001). Conclusion. Our study demonstrates that variations in hip geometry that affect the moment-generating capacity of the hip abductors are associated with the presence and compartmental distribution of tibiofemoral OA in the ipsilateral knee. Anatomical arrangements that reduce the hip abductor moment arm are more strongly associated with LOA than with MOA

Damage to articular cartilage is a common injury, for which there is no effective treatment. Our aims were to investigate the temporal sequence of the repair of articular cartilage and to define a critical-size defect.

Full-thickness defects were made in adult male New Zealand white rabbits. The diameter (1 to 4 mm) of the defects was varied in order to determine the effect that the size and depth of the defect had on its healing. The defects were made in the femoral groove of the knee with one defect per knee and eight knees per group. The tissues were fixed in formalin at days 3, 7, 14, 21, 28, 42, 84 and 126 after operation and the sections stained with Toluidine Blue. These were then examined and evaluated for several parameters including the degree of metachromasia and the amount of subchondral bone which had reformed in the defect.

The defects had a characteristic pattern of healing which differed at different days and for different sizes of defect. Specifically, the defects of 1 mm first peaked in terms of metachromasia at day 21, those of 2 mm at day 28, followed by defects of 3 mm and 4 mm. The healing of the subchondral bone was slowest in defects of 1 mm.

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

Abstract. Objectives. The fidelity of a 3D model created using image segmentation must be precisely quantified and evaluated for the model to be trusted for use in subsequent biomechanical studies such as finite element analysis. The bones within the ankle joint vary significantly in size and shape. The purpose of this study was to test the hypothesis that the accuracy and reliability of a segmented bone geometry is independent of the particular bone being measured. Methods. Computed tomography (CT) scan data (slice thickness 1 mm, pixel size 808±7 µm) from three anonymous patients was used for the development of the ankle geometries (consisting of the tibia, fibula, talus, calcaneus, and navicular bones) using Simpleware Scan IP software (Synopsys, Exeter, UK). Each CT scan was segmented 4 times by an inexperienced undergraduate, resulting in a total of 12 geometry assemblies. An experienced researcher segmented each scan once, and this was used as the ‘gold standard’ to quantify the accuracy. The solid bone geometries were imported into CAD software (Inventor 2023, Autodesk, CA, USA) for measurement of the surface area and volume of each bone, and the distances between bones (tibia to talus, talus to navicular, talus to calcaneus, and tibia to fibula) were carried out. The intra-class coefficient (ICC) was used to assess intra-observer reliability. Bland Altman plots were employed as a statistical measure for criteria validity (accuracy) [1]. Results. The average ICC score was 0.93, which is regarded as a high reliability score for an inexperienced user. The talus to navicular and talus to tibia separations, which had the smallest distances, showed a slight decrease in reliability and this was observed for all separations shorter than 2 mm. According to the Bland-Altman plots, more than 95% of the data points were inside the borders of agreement, which is an excellent indication of accuracy. The bias percentage (average error percentage) varied between 1% and 4% and was constant across all parameters, with the proportion rising for short distance separations. Conclusions. The current study demonstrates that an inexperienced undergraduate, with access to software manuals, can segment an ankle CT scan with excellent reliability. The present study also concluded that all five bones were segmented with high levels of accuracy, and this was not influenced by bone volume or type. The only factor found to influence the reliability was the magnitude of distance between bones, where if this was smaller than 2 mm it reduced the reliability, indicating the influence of CT scan resolution on the segmentation reliability. 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

Cells typically respond to a variety of geometrical cues in their environment, ranging from nanoscale surface topography to mesoscale surface curvature. The ability to control cellular organisation and fate by engineering the shape of the extracellular milieu offers exciting opportunities within tissue engineering. Despite great progress, however, many questions regarding geometry-driven tissue growth remain unanswered. Here, we combine mathematical surface design, high-resolution microfabrication, in vitro cell culture, and image-based characterization to study spatiotemporal cell patterning and bone tissue formation in geometrically complex environments. Using concepts from differential geometry, we rationally designed a library of complex mesostructured substrates (10. 1. -10. 3. µm). These substrates were accurately fabricated using a combination of two-photon polymerisation and replica moulding, followed by surface functionalisation. Subsequently, different cell types (preosteoblasts, fibroblasts, mesenchymal stromal cells) were cultured on the substrates for varying times and under varying osteogenic conditions. Using imaging-based methods, such as fluorescent confocal microscopy and second harmonic generation imaging, as well as quantitative image processing, we were able to study early-stage spatiotemporal cell patterning and late-stage extracellular matrix organisation. Our results demonstrate clear geometry-dependent cell patterning, with cells generally avoiding convex regions in favour of concave domains. Moreover, the formation of multicellular bridges and collective curvature-dependent cell orientation could be observed. At longer time points, we found clear and robust geometry-driven orientation of the collagenous extracellular matrix, which became apparent with second harmonic generation imaging after ∼2 weeks of culture. Our results highlight a key role for geometry as a cue to guide spatiotemporal cell and tissue organisation, which is relevant for scaffold design in tissue engineering applications. Our ongoing work aims at understanding the underlying principles of geometry-driven tissue growth, with a focus on the interactions between substrate geometry and mechanical forces

Intervertebral discs (IVD) provide flexibility to the back and ensure functional distributions of the spinal loads. They are avascular, and internal diffusion-dependent metabolic transport is vital to supply nutrients to disc cells1, but interactions with personalized IVD shapes and mechanics remain poorly explored. Poromechanical finite element models of seven personalized lumbar IVD geometries, with mean heights ranging from 8 to 16 mm were coupled with a reactive oxygen, glucose and lactate transport model linked with tissue deformations and osmosis . In previous studies, reduced formulations of the divergence of the solute flux (∇ .J = ∇ . (D∇ C) = ∇ D. ∇ C +D∇ 2C) ignored the dependence of the diffusion on the deformation gradients, ∇ D. ∇C. We simulated this phenomenon to explore its significance in mechano-metabolic -transport couplings, in the different geometries, over 24h of simulated rest (8h) and physical activity (16h). ∇ D. ∇ C affected the daily variations of glucose concentrations in IVD thinner than 12 mm but with neglectable variation ranges, while not considering ∇ D. ∇ C in taller discs only slightly overestimated the glucose concentration. Most importantly, tall IVD had nearly 60% less glucose than thin IVD, with local drops below the concentration of 0.5 mM, considered to be critical for disc cells3, in the anterior nucleus pulposus. On the one hand, previous reduced formulations for mechanometabolic-transport models of the IVD seem acceptable, even for patient-specific modelling. On the other hand, tall IVD might suffer from unfortunate combinations of deformation-dependent solute diffusion and large diffusion distances, which may favor early. Acknowledgements: Catalan Government and European Commission (2020 BP 00282; ERC-2021-CoG-O-Health-101044828)

Introduction. The human wrist is a highly complex joint, offering extensive motion across various planes. This study investigates scapholunate ligament (SLL) injuries’ impact on wrist stability and arthritis risks using cadaveric experiments and the finite element (FE) method. It aims to validate experimental findings with FE analysis results. Method. The study utilized eight wrist specimens on a custom rig to investigate Scapho-Lunate dissociation. Contact pressure and flexion were measured using sensors. A CT-based 3D geometry reconstruction approach was used to create the geometries needed for the FE analysis. The study used the Friedman test with pairwise comparisons to assess if differences between testing conditions were statistically significant. Result. The study found significant variations in scaphoid and lunate bone movement based on ligament condition. Full tears increased scapholunate distance in the distal-proximal direction and decreased in the medial-lateral direction. Lunate angles shifted from flexion to extension with fully torn ligaments. Conversely, the scaphoid shifted significantly from extension to flexion with full tears. A proximal movement was observed in the distal-proximal direction in all groups, with significant differences in the partial tear group. Lateral deviation of the scaphoid and lunate occurred with ligament damage, being more pronounced in the partial tear group. All groups exhibited statistically significant movement in the volar direction, with the full tear group showing the least movement. Also, radiocarpal joint and finger contact pressure and contact area were studied. Whereas the differences in contact area were not significant, scapholunate ligament tears resulted in significantly decreased finger contact pressures. FEA confirmed these findings, showing notable peak radiocarpal contact pressure differences between intact and fully torn ligaments. Conclusion. Our study found that SLL damage alters wrist stability, potentially leading to early arthritis. The FEA model confirmed these findings, indicating the potential for the clinical use of computer models from CT scans for treatment planning

There is currently no commercially available and clinically successful treatment for scapholunate interosseous ligament rupture, the latter leading to the development of hand-wrist osteoarthritis. We have created a novel biodegradable implant which fixed the dissociated scaphoid and lunate bones and encourages regeneration of the ruptured native ligament. To determine if scaphoid and lunate kinematics in cadaveric specimens were maintained during robotic manipulation, when comparing the native wrist with intact ligament and when the implant was installed. Ten cadaveric experiments were performed with identical conditions, except for implant geometry that was personalised to the anatomy of each cadaveric specimen. Each cadaveric arm was mounted upright in a six degrees of freedom robot using k-wires drilled through the radius, ulna, and metacarpals. Infrared markers were attached to scaphoid, lunate, radius, and 3rd metacarpal. Cadaveric specimens were robotically manipulated through flexion-extension and ulnar-radial deviation by ±40° and ±30°, respectively. The cadaveric scaphoid and lunate kinematics were examined with 1) intact native ligament, 2) severed ligament, 3) and installed implant. Digital wrist models were generated from computed tomography scans and included implant geometry, orientation, and location. Motion data were filtered and aligned relative to neutral wrist in the digital models of each specimen using anatomical landmarks. Implant insertion points in the scaphoid and lunate over time were then calculated using digital models, marker data, and inverse kinematics. Root mean squared distance was compared between severed and implant configurations, relative to intact. Preliminary data from five cadaveric specimens indicate that the implant reduced distance between scaphoid and lunate compared to severed configuration for all but three trials. Preliminary results indicate our novel implant reduced scapho-lunate gap caused by ligament transection. Future analysis will reveal if the implant can achieve wrist kinematics similar to the native intact wrist

The in vitro mimicking of bone microenvironment for the study of pathologies is a challenging field that requires the design of scaffolds with suitable morphological, structural and cytocompatible properties. During last years, 3D in vitro tumour models have been developed to reproduce mechanical, biochemical and structural bone microenvironment elements, allowing cells to behave as in vivo. In this work, gas foamed polyether urethane foams (PUF) and 3D printed thermoplastic polyether urethane (3DP-PU) designed with different patterns are proposed as scaffolds for in vitro model of bone tissue. Surface coatings for a biomimetic behaviour of the 3D scaffold models were also investigated. Morphological, chemico-physical, mechanical properties, and biological in vitro behaviour were investigated. PUFs for metastases investigation. The suitability of PUF as 3D in vitro model to study the interactions between bone tumour initiating cells and the bone microenvironment was investigated. PUF open porosity (>70%) appeared suitable to mimic trabecular bone structure. Human adipose derived stem cells (ADSC) were cultured and differentiated into osteoblast lineage on the PU foam, as confirmed by Alizarin Red staining and RT-PCR, thus offering a bone biomimetic microenvironment to the further co-culture with bone derived tumour-initiating cells (MCFS). Tumour aggregates were observed after three weeks of co-culture by e-cadherin staining and SEM; modification in CaP distribution was identified by SEM-EDX and associated to the presence of tumour cells. 3DP-PU as tumour bone model. 3D printed scaffolds have pores with a precise and regular geometry (0°-90°, 0°-45°-90°-135°, 0°-60°-120°). PU scaffold porosity evidenced values from 55 to 67%, values that belong to the porosity range of the trabecular bone tissue (30-90%). The compressive modulus varied between 2 and 4 MPa, depending on the printed pattern. Biomimetic nanostructured coating was performed on 0-90° 3DP-PU by Ionized Jet Deposition. Coatings had a submicrometric thickness, variable tuning deposition time, nanostructured surface morphology and biomimetic composition. Coating on 3DP-PU promoted cells colonization of the whole porous scaffolds, compared to the controls, where cells concentrated mostly on the outer layers. In conclusion, based on the obtained results, scaffolds with different geometries have been successfully produced. Morphological and structural properties of the scaffolds here presented are suitable for mimicking the bone tissue, in order to produce a 3D in vitro model useful for bone pathologies research

Geometric deep learning is a relatively new field that combines the principles of deep learning with techniques from geometry and topology to analyze data with complex structures, such as graphs and manifolds. In orthopedic research, geometric deep learning has been applied to a variety of tasks, including the analysis of imaging data to detect and classify abnormalities, the prediction of patient outcomes following surgical interventions, and the identification of risk factors for degenerative joint disease. This review aims to summarize the current state of the field and highlight the key findings and applications of geometric deep learning in orthopedic research. The review also discusses the potential benefits and limitations of these approaches and identifies areas for future research. Overall, the use of geometric deep learning in orthopedic research has the potential to greatly advance our understanding of the musculoskeletal system and improve patient care

Abstract. OBJECTIVES. Dual mobility (DM) total hip replacements (THRs) were introduced to reduce the risk of hip dislocation in at-risk patients. DM THRs have shown good overall survivorship and low rates of dislocation, however, the mechanisms which describe how these bearings function in-vivo are not fully understood. This is partly due to a lack of suitable characterisation methodologies which are appropriate for the novel geometry and function of DM polyethylene liners, whereby both surfaces are subject to articulation. This study aimed to develop a novel semi-quantitative geometric characterisation methodology to assess the wear/deformation of DM liners. METHODS. Three-dimensional coordinate data of the internal and external surfaces of 14 in-vitro tested DM liners was collected using a Legex 322 coordinate measuring machine. Data was input into a custom Matlab script, whereby the unworn reference geometry was determined using a sphere fitting algorithm. The analysis method determined the geometric variance of each point from the reference surface and produced surface deviation heatmaps to visualise areas of wear/deformation. Repeatability of the method was also assessed. RESULTS. Semi-quantitative analysis of the surface deviation heatmaps revealed circumferential damage patterns similar to those reported in the literature. Additionally, the location of the damaged regions corresponded between the internal and external surfaces. Comparing five repeat measurements of the same liner, the maximum geometric variance of each surface varied by 1 µm (standard deviation) suggesting a high repeatability of the method. CONCLUSIONS. This study presents an effective and highly repeatable characterisation methodology to semi-quantitatively assess the wear/deformation of in-vitro tested DM liners. This method is suitable for the analysis of retrieved DM liners whereby no pre-service information is available, which may provide information about the complex in-vivo kinematics and mechanical failure mechanisms of these bearings