1. Introduction
Osteoporosis (OP) is a degenerative bone condition that is thought to be responsible for 8.9 million fractures per year [
1]. It is estimated that one in two women and one in five men over the age of 50 will suffer a fragility fracture, which is defined as a fracture caused by a fall from standing height or less. These fractures are typically associated with or attributed to osteoporosis or osteopenia. In the UK alone there were approximately 527,000 new fragility fractures in 2019, estimated to increase by 26.2% to 665,000 in 2034, the cost similarly rising constantly by GBP~100 million/year from its present level of GBP~5.0 billion/year [
2]. At present, OP is defined as having a bone mass 2.5 standard deviations below the young adult reference mean [
3]. Another prevalent condition that affects bone tissue is osteoarthritis. Osteoarthritis is normally considered only for its impact on the articular cartilage of the synovial joints; the knock-on effects of the compromised joints causes structural changes to occur in the subchondral bone [
4]. Osteoarthritis affects 8.75 million people in the UK, and it is estimated that 33% of the population over the age of 45 have sought treatment for osteoarthritis. The joints most affected by the condition are the knee and hip affecting 4.7 and 2.46 million people, respectively.
Current protocol in determining a patient’s fracture risk and whether they are osteoporotic is based on dual energy X-ray absorptiometry (DEXA). This assessment using DEXA gives an indication of the patient’s bone mineral density (BMD) which is the product of both the porosity and density of the mineralized bone tissue; this is usually taken at the hip [
5]. The DEXA results are assessed using the fracture risk assessment tool as recommended by the World Health Organization. While this provides valuable data on an individual’s fracture risk, advancements in medical imaging technology allow for development of more robust and accurate risk assessment tools [
6].
The primary role of bone in the body is as a structural material and the cancellous regions can be considered as a cellular solid [
7,
8,
9,
10,
11]. As such, the mechanical properties of cancellous bone are impacted by the base material properties of the structure and the micro-architecture of the structure. All variations in the micro-structural properties of the tissue, from the quantity of bone tissue to the orientation of individual trabecular architecture, will impact the resultant mechanical properties of the tissue. The current DEXA protocol, however, fails to consider the architecture of the individual trabeculae. The most common mechanical property that is investigated is the compressive strength of the bone tissue, which fails to consider the ability of the tissue to resist fracture, an extremely important consideration when assessing the ability of bone to carry out its daily tasks, specifically its ability not to fracture under load. This has been considered by a previous study [
11], in which the fracture toughness of discs and beams of cancellous bone were measured, conforming to ASTM standards.
Multiple regression represents an advancement beyond traditional linear regression, because it allows the utilization of multiple predictors to estimate the value of a variable based on the values of two or more predictors. It also assesses the collective impact of multiple predictors on determining the outcome, providing a comprehensive understanding of the overall fit. It is a tool rarely used to predict the mechanical response of bone based on its architecture and has never before been used to predict the fracture toughness of bone in these terms. The authors recognize that in the application of multiple linear regression, the resultant models are not prescriptive of the underlying mechanisms but rather a descriptive method to ascertain the relationships within the sample set.
In a series of recent studies, we have demonstrated the importance and impact of changes in the micro-architecture and material properties in cancellous bone mechanics [
5,
12,
13,
14,
15,
16,
17], and these will provide a basis for the work that is presented here. In this study, we have the two following primary objectives: (a) investigate the use of predictive models to help in the prediction of fracture toughness, and (b) investigate if there are any significant differences between the models produced from samples loaded in the across (A
C) and along (A
L) loading configurations.
3. Results
Table 2 shows a descriptive and statistical comparison of the A
L and A
C as separate groups, as well as the average parameters collected for the entire cohort. Values measured between the groups were not statistically significantly different, with the exception of DA, which may be an artefact of the cutting and selection process. Even with this consideration in mind, it shows that the differences between the subsequent correlations and regression analysis is due to the contribution of the parameters to the loading in the specified direction. A comparison of the morphological data between males and females and with other studies was reported in a previous study [
5].
As well as considering the different loading conditions, differences between the OP and OA groups and the relationship they have with fracture toughness are also considered. Therefore, the relationships between the architectural properties of OP and OA bone and their corresponding fracture toughness are presented (
Table 3). Within the OP sample set, Kc had a higher correlation with trabecular spacing than observed in the OA group, whereas BV/TV and TbN correlations in the OA groups were much higher than in the OP. There was a consistency in parameters that correlated significantly between the OP and OA groups, except for connectivity density (Conn. D), which was found to be significant (
p < 0.05) in the OA group but not in the OP group.
Table 3 also includes a comparison of how the morphological data collected here compares to previous studies.
3.1. Micro-Architecture
Table 4 shows the correlations between the micro-architectural parameters and fracture toughness, with R
2 and
p-values given, in the A
L and A
C groups as well as in the combined groups. The parameter with the highest R
2 value in the A
L and combined groups was BV/TV whilst in the A
C group it was the TbN and BS/TV. Most of the parameters measured were found to impact upon fracture toughness, except for DA across the entire cohort and Conn. D in the A
C loading group. In
Table 5, the correlations within the OP and OA separate groups are shown. When considering the entire cohort, the BV/TV had the highest R
2 value in both groups. The DA and Conn. D did not correlate significantly in either group. Additionally, the Dmat was seen not to be significant with the OP/OA separation.
3.2. Regression Analysis
For multiple regressions, BV/TV was taken to be the base predictor as it was consistently the parameter that correlated highest with fracture toughness. Additionally, BV/TV is very closely linked to the metric currently used in the assessment of OP, as results from DEXA are mostly influenced by the quantity of bone rather than the density of the material itself [
5]. The performance of predictions produced by multiple regressions are shown in
Figure 2,
Figure 3 and
Figure 4. For A
C samples (
Figure 2) the inclusion of additional mArch parameters added to the predictive power of BV/TV alone. However, as shown in
Table 6, stepwise regression was unable to identify any additional mArch parameters that could significantly improve the R
2 value for the A
L group beyond BV/TV alone. The R
2 values for the A
L group for BV/TV and stepwise best fits were the same. The final best fit for the A
C group is given in
Table 7 and for the entire cohort in
Table 8.
The best models utilizing as many as possible variants, at a
p-value < 0.05, were step-3 (
Table 7) for the A
C group, and step-4 (
Table 8) for the entire cohort. There are two technical aspects of applying the stepwise regressions that are worth noting: (i) to add or subtract a parameter to the model, the
p-value was set to 0.15; and (ii) applying multiple regressions to the entire cohort adhered to the ‘rule of ten’ which suggests a minimum of ten samples for every predictor in the model. However, in the A
L groups model this was not maintained. Whilst it has been suggested that this is not necessary, maintaining a high number of predictors to samples is advantageous and helps reduce the effects of over fitting [
35].
4. Discussion
The research presented in this article outlines the fundamental relationships between the fracture toughness of cancellous bone and the material quality factors measured by μ-CT, and implements the use of statistical models to predict the mechanical properties of the samples. The collection of samples, which have been used in previous studies, are unique in that they are the only instance of measuring cancellous fracture toughness considering the start of growth of a major crack [
4,
11], as opposed to the total work under the load/deformation curve [
36]. Previous μ-CT research on this cohort has investigated the micro-architecture and material quality whilst looking at differences between the male and female samples in the cohort, and treated the samples loaded in different configurations (A
L and A
C) indiscriminately [
5]. Here, we have taken the opposite approach and treated the male/female samples indiscriminately and separated the A
L and A
C loading configurations. This research also has the inclusion of OA samples which represent perhaps the opposite of OP, in that the effects of OA tend to lead to a thickening of the subchondral bone. The use of μ-CT imaging presents an opportunity to assess the skeleton not currently found across the array of medical machines available. Current OP diagnosis by DEXA assesses the BMD which is a representation of the density of cancellous architecture and the material density of the bone itself. Medical CT scanners can also be used to assess the structure of the skeleton; however, the voxel size and image resolution currently obtainable from these systems is nowhere near as great as that which can be achieved in μ-CT. Therefore, all data and its associated methodology presented here represent what could potentially be assessed in the future, and are precursors to future non-invasive assessment of bone fracture toughness in diagnostic clinics, if we could only develop the ability to assess these same characteristics in vivo.
As previously mentioned, there is a real danger of over fitting data in a multiple regression analysis, which would produce models that claim to predict better than they are capable of. Here, we have taken every care to include the fewest number of predictors and to ensure that the predictors are independent of each other. In bone, however, this is very difficult due to the dependence of parameters on other physical characteristics, including both the obvious links between BV/TV and apparent density, and the less apparent links between the material density and the BV/TV [
34]. Multiple regression analysis was not carried out in the OP and OA subgroups due to a very small samples size of the OA group. The SMI values reported in this study were negative; this is due to the samples containing a significant number of concave surfaces. In the SMI calculation, it is assumed that the number of concave surfaces is negligible [
37]. Therefore, SMI was excluded from the multiple regression models due to lack of suitability but was included to demonstrate that the number of concave surfaces in cancellous bone are significant.
The comparison of OP and OA subgroups has supported the notion that OP leads to the loss of bone, shown by the significant differences between the BV/TV of the two groups (
p < 0.01). The average BV/TV of the OA group is still within the range previously reported in the literature (
Table 3), suggesting that in the OA condition, there is no extreme deposition of new bone tissue within the cancellous regions. The measured BS/BV and calculated BS/TV are within the literature ranges for both OP and OA groups. The differences between BS/BV for the OP and OA groups suggest that there are more surfaces available within the OP groups, which is consistent with the notion that remodeling is a surface effect [
29,
38]. Therefore, greater rates of remodeling could lead to a loss of bone which is typically associated with OP [
39,
40]. The trabecular number and thickness are higher in the OA group, which is typically consistent with the increased BV/TV and consistent with increases in mechanical strength [
9,
11,
38,
41]. There were no significant differences in the morphology measured between the A
L and A
C groups suggesting that any differences between correlations with the architectural parameters and any differences in the multiple regression models produced are due to the contributions of the individual parameters in the different loading directions.
When looking at the entire sample set, BV/TV was seen to have the highest correlation with fracture toughness, enforcing the assertion that the quantity of bone is the biggest contributor to bone strength. However, in the division of the A
L and A
C subgroups this only held true for the A
L group, whilst in the A
C group the TbN was seen to have the highest correlation. This suggests that a denser trabecular packing may have a bigger impact on cancellous bones’ resistance to fracture in the A
C loading configuration than in the A
L. The significant Dmat correlation across all the groups suggests that the material composition of the bone tissue plays an important role in the ability of the tissue to resist fracture, which supports previously found differences between the physio-chemistry of normal and OP bone tissue [
16]. However, the effect of Dmat is clearly not as important as the structural properties of the tissue as evidenced by the much lower R
2 values. Between the OP and OA subgroups, the parameters that impacted on fracture toughness followed the same trends, except for connectivity density (Conn. D), which was a significant contributor in the OA group but not the OP. This is perhaps due to the connectivity density being significantly higher in the OA group.
Multiple Regressions
Using multiple linear regressions, we were able to demonstrate that multiple morphological parameters impact upon the fracture toughness of bone when loaded in the AC direction or when loading direction is not considered. By accounting for these parameters within the model, it is possible to better predict the fracture toughness of bone than by consideration of multiple parameters. However, in the AL group, the use of multiple regression was unable to identify any parameter that would significantly improve the model. This has very profound implications on the understanding of bone fracture toughness and suggests that in the AL loading direction, the only parameter that resists fracture is the quantity of bone available, and that other parameters such as the average thickness of trabeculae do not develop in such a way to resist fracture. In the AC direction, however, other parameters had a significant effect on the ability of material to resist fracture. This is consistent with basic underpinning mechanisms of bone remodeling suggested by Wolff, whereby bone is responsive and adapts to the loads applied to it. The samples in the AC groups are orientated across the primary direction of loading in the hip, so the bone will have adapted to resist fracture in this direction and as such, this adaptation has led to reorientation of the trabeculae to achieve this. The AL group were orientated perpendicular to the primary loading of bone and as such, the trabecular structure has not adapted in micro-orientation to resist fracture.
The two primary aims of this study have been addressed as follows: (a) we have shown that across the entire cohort, consideration of multiple morphological parameters can help produce models that can inform on bone quality and can perhaps be used to predict fracture toughness with further development; and (b) separation of the models produced between the AL and AC groups was found to be revealing, in that for the AL group, no additional parameter was seen to improve the predictive ability over and above the use of BV/TV. This is incredibly surprising and has implications on our comprehension of how bone at the hip remodels to help resist fractures. To conclude, the use of multiple regressions represents a real opportunity to develop models to predict the likelihood of a patient’s fracture using bone micro-architecture, and there is a clear case for investigation into the remodeling of bone at the hip.
5. Conclusions
This study has considered the impact of the micro-architecture of cancellous bone on the fracture toughness of the tissue. We have been able to use a relatively large cohort of samples collected from patients undergoing hip replacement surgery and determined to be either osteoporotic or osteoarthritic. The findings support the currently used DEXA model, whereby a significantly reduced bone mass leads to a reduction in the mechanical competency of the tissue. It has additionally supported previous reports that multiple structural parameters such as TbTh, TbSp, TbN, and BS/BV also contribute significantly to the fracture toughness. We also employed the use of a statistical tool, multiple regression analysis, to demonstrate that the combination of multiple structural parameters can lead to an improved model of fracture toughness that may provide a basis to predict the fracture risk of a patient. The use of multiple regressions also highlighted that in the AL loading condition, the quantity of bone is the biggest contributor to fracture toughness and that the inclusion of additional parameters did not significantly improve the predictive power. The same cannot be said for the AC group, which showed a marked improvement with the addition of multiple parameters. This further proves Wolff’s law, or at least the principle, that bone truly remodels to its loading, and in this case, to resist fracture at the hip. The use of multiple regression is not without its limitations; in this study, from a statistical perspective, the sample size is relatively small; however, from a study on human bone samples perspective, it can be considered relatively large. Further work is required to investigate if these architectural parameters can be included alongside the currently collected BMD to improve the prediction of patients’ fracture risk.