Phantomless calibration of CT scans for hip fracture risk prediction in silico: Comparison with phantom-based calibration.

Finite element models built from quantitative computed tomography images rely on element-wise mapping of material properties starting from Hounsfield Units (HU), which can be converted into mineral densities upon calibration. While calibration is preferably carried out by scanning a phantom with known-density components, conducting phantom-based calibration may not always be possible. In such cases, a phantomless procedure, where the scanned subject's tissues are used as a phantom, is an interesting alternative. The aim of this study was to compare a phantom-based and a phantomless calibration method on 41 postmenopausal women. The proposed phantomless calibration utilized air, adipose, and muscle tissues, with reference equivalent mineral density values of -797, -95, and 38 mg/cm3, extracted from a previously performed phantom-based calibration. A 9-slice volume of interest (VOI) centred between the femoral head and knee rotation centres was chosen. Reference HU values for air, adipose, and muscle tissues were extracted by identifying HU distribution peaks within the VOI, and patient-specific calibration was performed using linear regression. Comparison of FE models calibrated with the two methods showed average relative differences of 1.99% for Young's modulus1.30% for tensile and 1.34% for compressive principal strains. Excellent correlations (R2 > 0.99) were identified for superficial maximum tensile and minimum compressive strains. Maximum normalised root mean square relative error (RMSRE) values settled at 4.02% for Young's modulus, 2.99% for tensile, and 3.22% for compressive principal strains, respectively. The good agreement found between the two methods supports the adoption of the proposed methodology when phantomless calibration is needed.

A statistical shape analysis for the assessment of the main geometrical features of the distal femoral medullary canal.

Statistical Shape Models (SSMs) are widely used in orthopedics to extract the main shape features from bone regions (e.g., femur). This study aims to develop an SSM of the femoral medullary canal, investigate its anatomical variability, and assess variations depending on canal length. The canals were isolated from 72 CT femur scans, through a threshold-based segmentation. A region of interest (ROI) was selected; sixteen segments were extracted from the ROI, ranging from 25% of the full length down to the most distal segment. An SSM was developed to identify the main modes of variation for each segment. The number of Principal Components (PCs) needed to explain at least 90% of the shape variance were three/four based on the length of the canal segment. The study examined the relationship between the identified PCs and geometric parameters like length, radius of curvature, ellipticity, mean diameter, and conicity, reporting range and percentage variation of these parameters for each segment. The SSMs provide insights into the anatomical variability of the femoral canal, emphasizing the importance of considering different segments to capture shape variations at various canal length. These findings can contribute for the design of personalized orthopedic implants involving the distal femur.

Development and validation of a semi-automated and unsupervised method for femur segmentation from CT.

Quantitative computed tomography (QCT)-based in silico models have demonstrated improved accuracy in predicting hip fractures with respect to the current gold standard, the areal bone mineral density. These models require that the femur bone is segmented as a first step. This task can be challenging, and in fact, it is often almost fully manual, which is time-consuming, operator-dependent, and hard to reproduce. This work proposes a semi-automated procedure for femur bone segmentation from CT images. The proposed procedure is based on the bone and joint enhancement filter and graph-cut algorithms. The semi-automated procedure performances were assessed on 10 subjects through comparison with the standard manual segmentation. Metrics based on the femur geometries and the risk of fracture assessed in silico resulting from the two segmentation procedures were considered. The average Hausdorff distance (0.03 ± 0.01 mm) and the difference union ratio (0.06 ± 0.02) metrics computed between the manual and semi-automated segmentations were significantly higher than those computed within the manual segmentations (0.01 ± 0.01 mm and 0.03 ± 0.02). Besides, a blind qualitative evaluation revealed that the semi-automated procedure was significantly superior (p < 0.001) to the manual one in terms of fidelity to the CT. As for the hip fracture risk assessed in silico starting from both segmentations, no significant difference emerged between the two (R2 = 0.99). The proposed semi-automated segmentation procedure overcomes the manual one, shortening the segmentation time and providing a better segmentation. The method could be employed within CT-based in silico methodologies and to segment large volumes of images to train and test fully automated and supervised segmentation methods.

DXA-based statistical models of shape and intensity outperform aBMD hip fracture prediction: A retrospective study.

Areal bone mineral density (aBMD) currently represents the clinical gold standard for hip fracture risk assessment. Nevertheless, it is characterised by a limited prediction accuracy, as about half of the people experiencing a fracture are not classified as at being at risk by aBMD. In the context of a progressively ageing population, the identification of accurate predictive tools would be pivotal to implement preventive actions. In this study, DXA-based statistical models of the proximal femur shape, intensity (i.e., density) and their combination were developed and employed to predict hip fracture on a retrospective cohort of post-menopausal women. Proximal femur shape and pixel-by-pixel aBMD values were extracted from DXA images and partial least square (PLS) algorithm adopted to extract corresponding modes and components. Subsequently, logistic regression models were built employing the first three shape, intensity and shape-intensity PLS components, and their ability to predict hip fracture tested according to a 10-fold cross-validation procedure. The area under the ROC curves (AUC) for the shape, intensity, and shape-intensity-based predictive models were 0.59 (95%CI 0.47-0.69), 0.80 (95%CI 0.70-0.90) and 0.83 (95%CI 0.73-0.90), with the first being significantly lower than the latter two. aBMD yielded an AUC of 0.72 (95%CI 0.59-0.82), found to be significantly lower than the shape-intensity-based predictive model. In conclusion, a methodology to assess hip fracture risk uniquely based on the clinically available imaging technique, DXA, is proposed. Our study results show that hip fracture risk prediction could be enhanced by taking advantage of the full set of information DXA contains.

Credibility assessment of computational models according to ASME V&V40: Application to the Bologna Biomechanical Computed Tomography solution.

BACKGROUND AND OBJECTIVE: When a computational model aims to be adopted beyond research purposes, e.g. to inform a clinical or regulatory decision, trust must be placed in its predictive accuracy. This practically translates into the need to demonstrate its credibility. In fact, prior to its adoption for regulatory purposes, an in silico methodology should be proven credible enough for the scope. This has become especially relevant as, although evidence of the safety and efficacy of new medical products or interventions has been traditionally provided to the regulator experimentally, i.e., in vivo or ex vivo, recently the idea to inform a regulatory decision in silico has made its way in the regulatory scenario. While a harmonised technical standard is currently missing in the EU regulatory system, in 2018 the ASME issued V&V40-2018, where a risk-based framework to assess the credibility of a computational model through the performance of predefined credibility activities is provided. The credibility framework is here applied to Bologna Biomechanical Computed Tomography (BBCT) solution, which predicts the absolute risk of fracture at the femur for a subject. BBCT has recently been the object of a qualification advice request to the European Medicine Agency. METHODS: The full implementation of ASME V&V40-2018 framework on BBCT is shown. Starting from BBCT proposed context of use the whole credibility plan is presented and the credibility activities (Verification, Validation, Applicability) described together with the achieved credibility levels. RESULTS: BBCT risk is judged medium, and the credibility levels achieved considered acceptable. The uncertainties intrinsically present in the material properties assignment affected BBCT predictions to the highest extent. CONCLUSIONS: This work provides the practical application of the ASME V&V40-2018 risk-based credibility assessment framework, which could be applied to demonstrate model credibility in any field and support future regulatory submissions and foster the adoption of In Silico Trials.

Experimental validation of a subject-specific finite element model of lumbar spine segment using digital image correlation.

Pathologies such as cancer metastasis and osteoporosis strongly affect the mechanical properties of the vertebral bone and increase the risk of fragility fractures. The prediction of the fracture risk with a patient-specific model, directly generated from the diagnostic images of the patient, could help the clinician in the choice of the correct therapy to follow. But before such models can be used to support any clinical decision, their credibility must be demonstrated through verification, validation, and uncertainty quantification. In this study we describe a procedure for the generation of such patient-specific finite element models and present a first validation of the kinematics of the spine segment. Quantitative computed tomography images of a cadaveric lumbar spine segment presenting vertebral metastatic lesions were used to generate the model. The applied boundary conditions replicated a specific experimental test where the spine segment was loaded in compression-flexion. Model predictions in terms of vertebral surface displacements were compared against the full-field experimental displacements measured with Digital Image Correlation. A good agreement was obtained from the local comparison between experimental data and simulation results (R2 > 0.9 and RMSE% <8%). In conclusion, this work demonstrates the possibility to apply the developed modelling pipeline to predict the displacement field of human spine segment under physiological loading conditions, which is a first fundamental step in the credibility assessment of these clinical decision-support technology.

Toward A Regulatory Pathway for the Use of in Silico Trials in the CE Marking of Medical Devices.

In Silico Trials methodologies will play a growing and fundamental role in the development and de-risking of new medical devices in the future. While the regulatory pathway for Digital Patient and Personal Health Forecasting solutions is clear, it is more complex for In Silico Trials solutions, and therefore deserves a deeper analysis. In this position paper, we investigate the current state of the art towards the regulatory system for in silico trials applied to medical devices while exploring the European regulatory system toward this topic. We suggest that the European regulatory system should start a process of innovation: in principle to limit distorted quality by different internal processes within notified bodies, hence avoiding that the more innovative and competitive companies focus their attention on the needs of other large markets, like the USA, where the use of such radical innovations is already rapidly developing.

Personalised 3D Assessment of Trochanteric Soft Tissues Improves HIP Fracture Classification Accuracy.

Passive soft tissues surrounding the trochanteric region attenuate fall impact forces and thereby control hip fracture risk. The degree of attenuation is related to Soft Tissue Thickness (STT). STT at the neutral hip impact orientation, estimated using a regression relation in body mass index (BMI), was previously shown to influence the current absolute risk of hip fracture (ARF0) and its fracture classification accuracy. The present study investigates whether fracture classification using ARF0 improves when STT is determined from the subject's Computed-Tomography (CT) scans (i.e. personalised) in an orientation-specific (i.e. 3D) manner. STT is calculated as the shortest distance along any impact orientation between a semi-automatically segmented femur surface and an automatically segmented soft tissue/air boundary. For any subject, STT along any of the 33 impact orientations analysed always exceeds the value estimated using BMI. Accuracy of fracture classification using ARF0 improves when using personalised 3D STT estimates (AUC = 0.87) instead of the BMI-based STT estimate (AUC = 0.85). The improvement is smaller (AUC = 0.86) when orientation-specificity of CT-based STT is suppressed and is nil when personalisation is suppressed instead. Thus, fracture classification using ARF0 improves when CT is used to personalise STT estimates and improves further when, in addition, the estimates are orientation specific.

Improving the Hip Fracture Risk Prediction with a Statistical Shape-and-Intensity Model of the Proximal Femur.

Severe predictions have been made regarding osteoporotic fracture incidence for the next years, with major economic and social impacts in a worldwide greying society. However, the performance of the currently adopted gold standard for fracture risk prediction, the areal Bone Mineral Density (aBMD), remains moderate. To overcome current limitations, the construction of statistical models of the proximal femur, based on three-dimensional shape and intensity (a hallmark of bone density), is here proposed for predicting hip fracture in a Caucasian postmenopausal cohort. Partial Least Square (PLS)-based statistical models of the shape, intensity and their combination were developed, and the corresponding modes and components were identified. Logistic regression models using the first two shape, intensity and shape-intensity PLS components were implemented and tested within a 10-fold cross-validation procedure as predictors of hip fracture. It emerged that (1) intensity components were superior to shape components in stratifying patients according to their fracture status, and that (2) a combination of intensity and shape improved patients risk stratification. The area under the ROC curve was 0.64, 0.85 and 0.92 for the models based on shape, intensity and shape-intensity combination respectively, against a 0.72 value for the aBMD standard approach. Based on these findings, the presented methodology turns out to be promising in tackling the need for an enhanced fracture risk assessment.

Combining shape and intensity dxa-based statistical approaches for osteoporotic HIP fracture risk assessment.

Aiming to improve osteoporotic hip fracture risk detection, factors other than the largely adopted Bone Mineral Density (BMD) have been investigated as potential risk predictors. In particular Hip Structural Analysis (HSA)-derived parameters accounting for femur geometry, extracted from Dual-energy X-ray Absorptiometry (DXA) images, have been largely considered as geometric risk factors. However, HSA-derived parameters represent discrete and cross-correlated quantities, unable to describe proximal femur geometry as a whole and tightly related to BMD. Focusing on a post-menopausal cohort (N = 28), in this study statistical models of bone shape and BMD distribution have been developed to investigate their possible role in fracture risk. Due to unavailable retrospective patient-specific fracture risk information, here a surrogate fracture risk based on 3D computer simulations has been employed for the statistical framework construction. When considered separately, BMD distribution performed better than shape in explaining the surrogate fracture risk variability for the analysed cohort. However, the combination of BMD and femur shape quantities in a unique statistical model yielded better results. In detail, the first shape-intensity combined mode identified using a Partial Least Square (PLS) algorithm was able to explain 70% of the surrogate fracture risk variability, thus suggesting that a more effective patients stratification can be obtained applying a shape-intensity combination approach, compared to T-score. The findings of this study strongly advocate future research on the role of a combined shape-BMD statistical framework in fracture risk determination.

Finite Element Modeling Application in Forensic Practice: A Periprosthetic Femoral Fracture Case Study.

The incidence of periprosthetic fractures has rapidly increased in the last two decades and has been the cause of a large number of revision surgeries and permanent physical disability for many patients, as well as a significant socioeconomic burden for many nations. This research deals with a periprosthetic femur fracture real event, occurred following a total hip arthroplasty and treated with one of the most widespread internal fixation methods: the implant of a periprosthetic femur plate system. A Finite Element analysis was performed to investigate the implanted femur plate break after a short follow-up and to understand the plate break causes. Such events are currently object of forensic debate as more and more often hospitals, surgeons, and medical device manufacturers are denounced by patients to whom similar events occur. In this work, different load situations acting on the femur during daily and incidental activities were simulated, in order to validate the correct behavior of the plate, according to the intended use recommended by the manufacturer. The analysis demonstrates that the plate failure can occur in situations of unconventional loading such as that caused by stumbling and in presence of incomplete bone healing.

Decellularized Human Dermal Matrix as a Biological Scaffold for Cardiac Repair and Regeneration.

The complex and highly organized environment in which cells reside consists primarily of the extracellular matrix (ECM) that delivers biological signals and physical stimuli to resident cells. In the native myocardium, the ECM contributes to both heart compliance and cardiomyocyte maturation and function. Thus, myocardium regeneration cannot be accomplished if cardiac ECM is not restored. We hypothesize that decellularized human skin might make an easily accessible and viable alternate biological scaffold for cardiac tissue engineering (CTE). To test our hypothesis, we decellularized specimens of both human skin and human myocardium and analyzed and compared their composition by histological methods and quantitative assays. Decellularized dermal matrix was then cut into 600-µm-thick sections and either tested by uniaxial tensile stretching to characterize its mechanical behavior or used as three-dimensional scaffold to assess its capability to support regeneration by resident cardiac progenitor cells (hCPCs) in vitro. Histological and quantitative analyses of the dermal matrix provided evidence of both effective decellularization with preserved tissue architecture and retention of ECM proteins and growth factors typical of cardiac matrix. Further, the elastic modulus of the dermal matrix resulted comparable with that reported in literature for the human myocardium and, when tested in vitro, dermal matrix resulted a comfortable and protective substrate promoting and supporting hCPC engraftment, survival and cardiomyogenic potential. Our study provides compelling evidence that dermal matrix holds promise as a fully autologous and cost-effective biological scaffold for CTE.

Improving the Hip Fracture Risk Prediction Through 2D Finite Element Models From DXA Images: Validation Against 3D Models.

Osteoporotic fracture incidence represents a major social and economic concern in the modern society, where the progressive graying of the population involves an highly increased fracture occurrence. Although the gold standard to diagnose osteoporosis is represented by the T-score measurement, estimated from the Bone Mineral Density (BMD) using Dual-energy X-ray Absorptiometry (DXA), the identification of the subjects at high risk of fracture still remains an issue. From this perspective, the purpose of this work is to investigate the role that DXA-based two-dimensional patient-specific finite element (FE) models of the proximal femur, in combination with T-score, could play in enhancing the risk of fracture estimation. With this aim, 2D FE models were built from DXA images of the 28 post-menopausal female subjects involved. A sideways fall condition was reproduced and a Risk of Fracture ( RF ^ ) was computed on the basis of principal strains criteria. The identified RF ^ was then compared to that derived from the CT-based models developed in a previous study. The 2D and 3D RF ^ turned out to be significantly correlated (Spearman's ρ = 0.66, p < 0.001), highlighting the same patients as those at higher risk. Moreover, the 2D RF ^ resulted significantly correlated with the T-score (Spearman's ρ = -0.69, p < 0.001), and managed to better differentiate osteopenic patients, drawing the attention to some of them. The Hip Structural Analysis (HSA) variables explaining the majority of the variance of the 2D and 3D fracture risk were the same as well, i.e., neck-shaft angle and narrow neck buckling ratio. In conclusion, DXA-based FE models, developable from currently available clinical data, appear promising in supporting and integrating the present diagnostic procedure.

Osteoporotic Hip Fracture Prediction: Is T-Score-Based Criterion Enough? A Hip Structural Analysis-Based Model.

At present, the current gold-standard for osteoporosis diagnosis is based on bone mineral density (BMD) measurement, which, however, has been demonstrated to poorly estimate fracture risk. Further parameters in the hands of the clinicians are represented by the hip structural analysis (HSA) variables, which include geometric information of the proximal femur cross section. The purpose of this study was to investigate the suitability of HSA parameters as additional hip fracture risk predictors. With this aim, twenty-eight three-dimensional patient-specific models of the proximal femur were built from computed tomography (CT) images and a sideways fall condition was reproduced by finite element (FE) analyses. A tensile or compressive predominance based on minimum and maximum principal strains was determined at each volume element and a risk factor (RF) was calculated. The power of HSA variables combinations to predict the maximum superficial RF values was assessed by multivariate linear regression analysis. The optimal regression model, identified through the Akaike information criterion (AIC), only comprises two variables: the buckling ratio (BR) and the neck-shaft angle (NSA). In order to validate the study, the model was tested on two additional patients who suffered a hip fracture after a fall. The results classified the patients in the high risk level, confirming the prediction power of the adopted model.

Implementation and validation of constitutive relations for human dermis mechanical response.

Finite element models in conjunction with adequate constitutive relations are pivotal in several physiological and medical applications related to both native and engineered tissues, allowing to predict the tissue response under various loading states. In order to get reliable results, however, the validation of the constitutive models is crucial. Therefore, the main purpose of this work is to provide an experimental-computational approach to the biomechanical investigation of soft tissues such as the dermis. This is accomplished by implementing and validating three widely adopted hyperelastic constitutive models (the Ogden, the Holzapfel, and the Gasser-Ogden-Holzapfel laws) supposed to be adequate to reproduce human reticular dermis mechanical behavior. Biaxial experimental data have represented the basis for the determination of the respective material parameters identified thanks to the definition of a cost function accounting for the discrepancy between experimental and predicted data. Afterwards, the experimental tests have been reproduced through finite element simulations. Hence, the constitutive laws have been validated comparing experimental and numerical outcomes in terms of displacements of four reference points and stress-strain relations. Hence, an experimental-numerical framework is proposed for the investigation of collagenous tissues, which could become more accurate with larger and independent experimental datasets. Graphical abstract ᅟ.