Next Article in Journal
Exploring the Efficacy of Novel Therapeutic Strategies for Periodontitis: A Literature Review
Previous Article in Journal
Analysing the Influence of Speed and Jumping Performance Metrics on the Percentage Change of Direction Deficit in Adolescent Female Soccer Players
Previous Article in Special Issue
Parametrization of the Calcaneus and Medial Cuneiform to Aid Potential Advancements in Flatfoot Surgery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting the Fracture Toughness of Human Cancellous Bone in Fractured Neck of Femur Patients Using Bone Volume and Micro-Architecture

1
Cranfield Forensic Institute, Cranfield University, Cranfield MK43 0AL, UK
2
nCATS, School of Engineering Science, University of Southampton, Southampton SO17 1BJ, UK
3
Biomedical Engineering Research Group, School of Engineering, University of Hull, Kingston upon Hull HU6 7RX, UK
*
Author to whom correspondence should be addressed.
Life 2024, 14(4), 467; https://doi.org/10.3390/life14040467
Submission received: 18 December 2023 / Revised: 20 March 2024 / Accepted: 29 March 2024 / Published: 3 April 2024
(This article belongs to the Special Issue Hard and Soft Tissue Biomechanics ‘In Translation’)

Abstract

:
The current protocol used to determine if an individual is osteoporotic relies on assessment of the individual’s bone mineral density (BMD), which allows clinicians to judge the condition of a patient with respect to their peers. This, in essence, evaluates a person’s fracture risk, because BMD is a good surrogate measure for strength and stiffness. In recent studies, the authors were the first to produce fracture toughness (FT) data from osteoporotic (OP) and osteoarthritic (OA) patients, by using a testing technique which basically analyzes the prerequisite stress conditions for the onset of growth of a major crack through cancellous bone tissue. FT depends mainly on bone quantity (BV/TV, bone volume/tissue volume), but also on bone micro-architecture (mArch), the inner trabecular design of the bone. The working research hypothesis of the present study is that mArch offers added prediction power to BV/TV in determining FT parameters. Consequently, our aim was to investigate the use of predictive models for fracture toughness and also to investigate if there are any significant differences between the models produced from samples loaded across (AC, transverse to) the main trabecular orientation and along (AL, in parallel) the trabeculae. In multilinear regression analysis, we found that the strength of the relationships varied for a crack growing in these two orthogonal directions. Adding mArch variables in the Ac direction helped to increase the R2 to 0.798. However, in the AL direction, adding the mArch parameters did not add any predictive power to using BV/TV alone; BV/TV on its own could produce R2 = 0.730. The present results also imply that the anisotropic layout of the trabeculae makes it more difficult for a major crack to grow transversely across them. Cancellous bone models and remodels itself in a certain way to resist fracture in a specific direction, and thus, we should be mindful that architectural quality as well as bone quantity are needed to understand the resistance to fracture.

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 (AC) and along (AL) loading configurations.

2. Materials and Methods

2.1. Bone Specimens

A sample set of femoral heads were collected from 37 osteoporotic (OP) and 8 osteoarthritic (OA) patients who had received a total hip replacement surgery due to suffering fragility fractures at the femoral neck (OP) or elective surgery for arthroplasty (OA). During the surgery, specialist surgeons were able to remove the femoral heads intact. The femoral head was used in this study due to the increased availability of tissue at the site compared to the femoral neck, where fracture typically occurs, whilst being physically close to the femoral neck. This cohort allowed for the collection of a number of samples, covering a broad range of porosity values, with OP covering the lower end of the spectrum and OA the normal and higher end of porosities. A good overlap was shown between the two groups. Population characteristics are provided in Table 1. Following removal, all samples were kept at −20 °C prior to sample preparation. Ethical approval for the collection and use of these specimens was provided by Gloucestershire NHS trust REC (acknowledgments).

2.2. Specimen Preparation

Specimen preparation (including sectioning from the femoral head and cleaning) has previously been described in detail [4,11]. Single Edge Notched Disc (SEND) samples were prepared to conform to an adjusted ASTM standard E399-90, in order to assess the necessary stress conditions to instigate crack growth from a man-made notch (Figure 1). Samples were divided into two subsets; with 34 samples orientated along (AL) the primary direction of the trabecular and 27 orientated across (AC) the primary orientation of the trabecular structure. Due to it being a cellular solid comprising struts and beams, directionality is a natural feature of cancellous bone. When the struts and beams are not isotropically laid out, they result in an anisotropic structure, where directions such as along and across the majority of the beams can be discerned. In a few femoral heads, samples were produced in both AL and AC directions. All specimens were stored at −20 °C following a defatting process detailed in [4,11]. The sectioning was performed by using a metallurgical saw (Struers® Accutom-2, Rotherham, UK). The samples were then sanded and polished by using progressively finer grades of carbide paper (400–2500 grit) to the dimensions required for material testing. Specimens were manufactured in the shape of discs, diameter 20 mm and thickness 7.5 mm, for mechanical material testing as SEND. Sample preparation was performed under constant water irrigation, to prevent the production of micro-cracks or other damage to the specimens.

2.3. Micro-Computed Tomography

The samples’ micro-architecture was imaged using cone beam micro-computed tomography (μ-CT, μ-CBCT). Each sample was imaged using a Nikon CT H225 (X-Tek Systems Ltd., Tring, Hertfordshire, UK) cone beam μ-CT (μ-CBCT) scanner. Samples were imaged at 50 kV and 65 μA with a 1000 ms exposure. The resultant voxel size of the scan was ~24 μm. All scans were manually reconstructed using CT Pro 3D. During reconstruction, conditions were optimized to reduce beam hardening and noise, and the noise and beam hardening corrections were standardized across all the samples to ensure that the results were comparable. Image analysis and visualization were carried out using VG Studio Max 2.2 (Volume Graphics GmbH, Heidelberg, Germany). Firstly, the samples’ structural properties were determined, and these parameters included the following: trabecular thickness (TbTh), spacing (TbSp) and number (TbN), surface area (BS), material volume (BV), and total volume (TV). The density of the samples was also determined using a QRM MicroCT-HA (QRM GmbH, Möhrendorf, Germany) calibration phantom. This uses hydroxyapatite of different known concentrations to produce a calibration curve of gray value versus density. Using this calibration, the density of the samples can be determined, which is often referred to as tissue mineral density or material density (Dmat). Following the determination of the Dmat from the average gray value, the apparent density (Dapp) of each sample was determined using Equation (1). The Dapp is often referred to as bone mineral density.
Dapp = Dmat × BV/TV
In this article, Dmat and Dapp are used to indicate that they have been calculated on a volumetric basis using μ-CT data. This is important to distinguish, as bone mineral density measured in DEXA is an areal representation of this bone mArch parameter. BoneJ© [http://bonej.org/; http://rsbweb.nih.gov/ij/] (30 June 2018) was employed at a second stage to calculate additional micro-architectural parameters such as structure model index (SMI), degree of anisotropy (DA), connectivity density (Conn. D) and Euler characteristic (Euler ch.) [18].

2.4. Mechanical Testing

The SEND samples were mechanically characterized, as in previous papers, for fracture toughness using a linear elastic FM approach [4,11]. The Kc (critical stress intensity factor) values were derived for the load at a point where the man-made notch started growing (following extensive yielding and bending of the trabeculae ahead of the notch), caused by snapping of one or more trabeculae in the first instance. The deformation was measured by a miniature extensometer (Model 3442-006M-050ST, DWE Scientific Ltd., Brackley, UK) attached at the mouth of the notch. The dimensions and other restrictions that were followed complied with the usual material testing standards such as ASTM E399-90, as reported in [11]. The mechanical testing was undertaken using a DARTEC® Series HC25 (Zwick Roell, Worcester, UK) materials testing machine driven by a 9610 series controller unit and operated using Workshop 96© software. The load was monitored using a 500 N load cell (RDP Electronics Ltd., Wolverhampton, UK) whilst the gauge length of the crack mouth opening displacement measured by the extensometer was 6 mm. The loading rate during fracture toughness testing was 0.05 mm s−1 (3 mm min−1), with data acquisition at a capture rate of 1000 points per minute. Unlike the more common compression studies, which have tested cancellous bone in cylinders or cubes, these tests were the first ever to attempt a quantification of the necessary loading conditions that would allow a crack to start growing from stability into an unstable fracture mode. In this respect, the mechanical data offer a novel and invaluable way of assessing the structural integrity and loading ability of these samples, in a way that resembles the conditions in FNF situations in a more biofidelic manner [11].

2.5. Statistical Analysis

Statistical analysis was carried out in two steps. Firstly, comparisons were made among the multiple subgroups within the cohort, such as the AL and AC loading configurations, and the OP and OA groups. Pearson’s correlation coefficients were used throughout the study. In the second stage of statistical analysis, multiple linear regressions were employed to explore to what degree BV/TV and the full set of mArch parameters are able to predict Kc. The multiple regression was a stepwise selection in MINITAB© (v.17). The stepwise selection combined forward selection and backward elimination, where added variables were deleted if their contribution to the model is not determined to be significant. It must be considered that the values for trabecular thickness (TbTh) and trabecular number (TbN) are calculated interdependently; therefore, in the development of statistical models, inclusion of both has been avoided.

3. Results

Table 2 shows a descriptive and statistical comparison of the AL and AC 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 R2 and p-values given, in the AL and AC groups as well as in the combined groups. The parameter with the highest R2 value in the AL and combined groups was BV/TV whilst in the AC 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 AC 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 R2 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 AC 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 R2 value for the AL group beyond BV/TV alone. The R2 values for the AL group for BV/TV and stepwise best fits were the same. The final best fit for the AC 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 AC 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 AL 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 (AL and AC) indiscriminately [5]. Here, we have taken the opposite approach and treated the male/female samples indiscriminately and separated the AL and AC 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 AL and AC 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 AL and AC subgroups this only held true for the AL group, whilst in the AC 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 AC loading configuration than in the AL. 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 R2 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.

Author Contributions

Conceptualization: G.J.A. and P.Z.; methodology: G.J.A., R.B.C., M.G. and P.Z.; formal analysis: G.J.A., MG. and P.Z.; writing—original draft preparation: G.J.A.; writing—review and editing: R.B.C. and P.Z.; supervision: P.Z.; project administration: P.Z.; funding acquisition: P.Z. All authors have worked towards the results of the manuscript and agreed to the contents. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge the support provided by the UK Department of Transport under the BOSCOS (Bone Scanning for Occupant Safety) program and the Engineering and Physical Science Research Council (EPSRC under GR/N33225; GR/N33102; GR/M59167; and EP/K020196/1).

Institutional Review Board Statement

The study was approved by Cranfield University Research and Ethics system (CURES/2549/2017, 24 April 2017; Bone mechanobiology research studies group: ageing and diseased bone project); the human specimens were obtained with informed consent in the Gloucester and Cheltenham NHS Trust hospitals under ethical consent (BOSCOS—Mr. Curwen CI REC ref. 01/179G).

Data Availability Statement

Data for this manuscript when the article is in print will be available through the Cranfield University CORD data depository and preservation system at https://cranfield.figshare.com, or through the corresponding author [email protected].

Acknowledgments

The work and reporting were carried out in the Biomechanics and Analytical laboratories of the Cranfield Forensic Institute of Cranfield University and the Biomedical Engineering Research Group of University of Hull, UK. This paper is dedicated to our friend and colleague Mike Gibson, whose untimely death is a great loss to us all.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Johnell, O.; Kanis, J.A. An estimate of the worldwide prevalence and disability associated with osteoporotic fractures. Osteoporos. Int. 2006, 17, 1726–1733. [Google Scholar] [CrossRef] [PubMed]
  2. International Osteoporosis Foundation. Facts and Statistics|International Osteoporosis Foundation. 2015. Available online: https://www.osteoporosis.foundation/facts-statistics (accessed on 30 June 2023).
  3. Summers, G.D. Osteoporosis in men. Radiography 2001, 7, 119–123. [Google Scholar] [CrossRef]
  4. Cook, R.B. Non-Invasively Assessed Skeletal Bone Status and Its Relationship to the Biomechanical Properties and Condition of Cancellous Bone. Ph.D. Thesis, Cranfield University, Cranfield, Bedford, UK, 2005. [Google Scholar]
  5. Greenwood, C.; Clement, J.; Dicken, A.; Evans, J.; Lyburn, I.; Martin, R.; Rogers, K.; Stone, N.; Adams, G.; Zioupos, P. The micro-architecture of human cancellous bone from fracture neck of femur patients in relation to the structural integrity and fracture toughness of the tissue. Bone Rep. 2015, 3, 67–75. [Google Scholar] [CrossRef]
  6. van den Bergh, J.P.W.; van Geel, T.A.C.M.; Lems, W.F.; Geusens, P.P. Assessment of Individual Fracture Risk: FRAX and Beyond. Curr. Osteoporos. Rep. 2010, 8, 131–137. [Google Scholar] [CrossRef]
  7. Huang, J.S.; Gibson, L.J. Fracture toughness of brittle foams. Acta Metall. Mater. 1991, 39, 1627–1636. [Google Scholar] [CrossRef]
  8. Maiti, S.K.; Ashby, M.F.; Gibson, L.J. Fracture toughness of brittle cellular solids. Scr. Metall. 1984, 18, 213–217. [Google Scholar] [CrossRef]
  9. Gibson, L.J. The mechanical behaviour of cancellous bone. J. Biomech. 1985, 18, 317–328. [Google Scholar] [CrossRef] [PubMed]
  10. Brezny, R.; Green, D.J. The effect of cell size on the mechanical behavior of cellular materials. Acta Metall. Mater. 1990, 38, 2517–2526. [Google Scholar] [CrossRef]
  11. Cook, R.B.; Zioupos, P. The fracture toughness of cancellous bone. J. Biomech. 2009, 42, 2054–2060. [Google Scholar] [CrossRef]
  12. Adams, G.J.; Cook, R.B.; Hutchinson, J.R.; Zioupos, P. Microarchitecture and morphology of mammalian bone tissue over a wide range of apparent densities assessed by micro-computed tomography and three different threshold backgrounds. Med. Eng. Phys. 2022, 106, 103828. [Google Scholar]
  13. Adams, G.J.; Cook, R.B.; Hutchinson, J.R.; Zioupos, P. Bone apparent and material densities examined by cone beam computed tomography and the Archimedes technique: Comparison of the two methods and their results. Front. Mech. Eng. 2018, 3, 23. [Google Scholar] [CrossRef]
  14. Adams, G.J. Quality Factors at the Material and Structural Level That Affect the Toughness of Human Cancellous Bone. Ph.D. Thesis, Cranfield University, Cranfield, Bedford, UK, 2017. [Google Scholar]
  15. Greenwood, C.; Clement, J.; Dicken, A.; Evans, P.; Lyburn, I.; Martin, R.M.; Stone, N.; Zioupos, P.; Rogers, K. Age-related changes in femoral head trabecular microarchitecture. Aging Dis. 2018, 9, 976–987. [Google Scholar]
  16. Greenwood, C.; Clement, J.; Dicken, A.; Evans, J.; Lyburn, I.; Martin, R.; Rogers, K.; Stone, N.; Zioupos, P. Towards new material biomarkers for fracture risk. Bone 2016, 93, 55–63. [Google Scholar] [CrossRef] [PubMed]
  17. McGivern, H.; Greenwood, C.; Márquez-Grant, N.; Kranioti, E.F.; Xhemali, B.; Zioupos, P. Age-related trends in the trabecular micro-architecture of the medial clavicle: Is it of use in forensic science? Front. Bioeng. Biotechnol. 2020, 7, 467. [Google Scholar] [CrossRef] [PubMed]
  18. Doube, M.; Kłosowski, M.M.; Wiktorowicz-Conroy, A.M.; Hutchinson, J.R.; Shefelbine, S.J. Trabecular bone scales allometrically in mammals and birds. Proc. Biol. Sci. 2011, 278, 3067–3073. [Google Scholar] [CrossRef] [PubMed]
  19. Link, T.M.; Majumdar, S.; Lin, J.C.; Newitt, D.; Augat, P.; Ouyang, X.; Mathur, A.; Genant, H.K. A Comparative Study of Trabecular Bone Properties in the Spine and Femur Using High Resolution MRI and CT. J. Bone Miner. Res. 1998, 13, 122–132. [Google Scholar] [CrossRef] [PubMed]
  20. Thomsen, J.S.; Jensen, M.V.; Niklassen, A.S.; Ebbesen, E.N.; Brüel, A. Age-related changes in vertebral and iliac crest 3D bone microstructure—Differences and similarities. Osteoporos. Int. 2015, 26, 219–228. [Google Scholar] [CrossRef] [PubMed]
  21. Wu, Z.; Laneve, A.J.; Niebur, G.L. In vivo microdamage is an indicator of susceptibility to initiation and propagation of microdamage in human femoral trabecular bone. Bone 2013, 55, 208–215. [Google Scholar] [CrossRef] [PubMed]
  22. Perilli, E.; Baleani, M.; Ohman, C.; Fognani, R.; Baruffaldi, F.; Viceconti, M. Dependence of mechanical compressive strength on local variations in microarchitecture in cancellous bone of proximal human femur. J. Biomech. 2008, 41, 438–446. [Google Scholar] [CrossRef] [PubMed]
  23. Macdonald, H.M.; Nishiyama, K.K.; Kang, J.; Hanley, D.A.; Boyd, S.K. Age-related patterns of trabecular and cortical bone loss differ between sexes and skeletal sites: A population-based HR-pQCT study. J. Bone Miner. Res. 2011, 26, 50–62. [Google Scholar] [CrossRef]
  24. Green, J.O.; Nagaraja, S.; Diab, T.; Vidakovic, B.; Guldberg, R.E. Age-related changes in human trabecular bone: Relationship between microstructural stress and strain and damage morphology. J. Biomech. 2011, 44, 2279–2285. [Google Scholar] [CrossRef]
  25. Mazurkiewicz, A.; Topoliński, T. Relationships between structure, density and strength of human trabecular bone. Acta Bioeng. Biomech. 2009, 11, 55–61. [Google Scholar] [PubMed]
  26. Zhang, Z.-M.; Li, Z.-C.; Jiang, L.-S.; Jiang, S.-D.; Dai, L.-Y. Micro-CT and mechanical evaluation of subchondral trabecular bone structure between postmenopausal women with osteoarthritis and osteoporosis. Osteoporos. Int. 2010, 21, 1383–1390. [Google Scholar] [CrossRef] [PubMed]
  27. Milovanovic, P.; Djonic, D.; Marshall, R.P.; Hahn, M.; Nikolic, S.; Zivkovic, V.; Amling, M.; Djuric, M. Micro-structural basis for particular vulnerability of the superolateral neck trabecular bone in the postmenopausal women with hip fractures. Bone 2012, 50, 63–68. [Google Scholar] [CrossRef] [PubMed]
  28. Fyhrie, D.P.; Lang, S.M.; Hoshaw, S.J.; Schaffler, M.B.; Kuo, R.F. Human vertebral cancellous bone surface distribution. Bone 1995, 17, 287–291. [Google Scholar] [CrossRef] [PubMed]
  29. Martin, R.B. Porosity and specific surface of bone. Crit. Rev. Biomed. Eng. 1984, 10, 179–222. [Google Scholar] [PubMed]
  30. Ding, M.; Hvid, I. Quantification of age-related changes in the structure model type and trabecular thickness of human tibial cancellous bone. Bone 2000, 26, 291–295. [Google Scholar] [CrossRef] [PubMed]
  31. Yan, Y.-B.; Qi, W.; Wang, J.; Liu, L.-F.; Teo, E.-C.; Tianxia, Q.; Ba, J.-J.; Lei, W. Relationship between architectural parameters and sample volume of human cancellous bone in micro-CT scanning. Med. Eng. Phys. 2011, 33, 764–769. [Google Scholar] [CrossRef] [PubMed]
  32. Saers, J.P.P.; Cazorla-Bak, Y.; Shaw, C.N.; Stock, J.T.; Ryan, T.M. Trabecular bone structural variation throughout the human lower limb. J. Hum. Evol. 2016, 97, 97–108. [Google Scholar] [CrossRef]
  33. Galante, J.; Rostoker, W.; Ray, R.D. Physical properties of trabecular bone. Calcif. Tissue Res. 1970, 5, 236–246. [Google Scholar] [CrossRef] [PubMed]
  34. Zioupos, P.; Cook, R.B.; Hutchinson, J.R. Some basic relationships between density values in cancellous and cortical bone. J. Biomech. 2008, 41, 1961–1968. [Google Scholar] [CrossRef]
  35. Austin, P.C.; Steyerberg, E.W. The number of subjects per variable required in linear regression analyses. J. Clin. Epidemiol. 2015, 68, 627–636. [Google Scholar] [CrossRef] [PubMed]
  36. Garrison, J.G.; Slaboch, C.L.; Niebur, G.L. Density and architecture have greater effects on the toughness of trabecular bone than damage. Bone 2009, 44, 924–929. [Google Scholar] [CrossRef] [PubMed]
  37. Salmon, P.L.; Ohlsson, C.; Shefelbine, S.J.; Doube, M. Structure Model Index Does Not Measure Rods and Plates in Trabecular Bone. Front. Endocrinol. 2015, 6, 162. [Google Scholar] [CrossRef] [PubMed]
  38. Currey, J.D. Bones: Structure and Mechanics; Princeton University Press: Princeton, NJ, USA, 2002. [Google Scholar]
  39. Genant, H.K.; Engelke, K.; Prevrhal, S. Advanced CT bone imaging in osteoporosis. Rheumatology 2008, 47 (Suppl. S4), iv9–iv16. [Google Scholar] [CrossRef] [PubMed]
  40. Buenzli, P.R.; David, C.; Thomas, L.; Clement, J.G.; Pivonka, P. Endocortical bone loss in osteoporosis: The role of bone surface availability. Int. J. Numer. Meth. Biomed. Eng. 2013, 29, 1307–1322. [Google Scholar] [CrossRef]
  41. Kopperdahl, D.L.; Keaveny, T.M. Yield strain behavior of trabecular bone. J. Biomech. 1998, 31, 601–608. [Google Scholar] [CrossRef]
Figure 1. (a) Cancellous bone samples have their trabeculae oriented, in the majority of cases, along the direction of the principal axes of stress and thus they are anisotropic. A propagating crack will find its way either advancing between the trabeculae (along direction AL) or transversely across the trabeculae (across direction AC); (b) design template for preparing the disk-shaped fracture toughness specimens [4]; (c) example of a cancellous bone sample prepared according to the template; (d) sample mounted on the grip, under continuous physiological fluid irrigation and with extensometer attached at the mouth of the crack; (e,f) μ-CT scanned section for the portion of the sample along the path of the growing crack (magnified 5.25), this particular sample is an AL oriented one; where (e) single slice along a sagittal plane through the sample; and (f) a 3D rendering of the sample surface using all the μ-CT data.
Figure 1. (a) Cancellous bone samples have their trabeculae oriented, in the majority of cases, along the direction of the principal axes of stress and thus they are anisotropic. A propagating crack will find its way either advancing between the trabeculae (along direction AL) or transversely across the trabeculae (across direction AC); (b) design template for preparing the disk-shaped fracture toughness specimens [4]; (c) example of a cancellous bone sample prepared according to the template; (d) sample mounted on the grip, under continuous physiological fluid irrigation and with extensometer attached at the mouth of the crack; (e,f) μ-CT scanned section for the portion of the sample along the path of the growing crack (magnified 5.25), this particular sample is an AL oriented one; where (e) single slice along a sagittal plane through the sample; and (f) a 3D rendering of the sample surface using all the μ-CT data.
Life 14 00467 g001
Figure 2. Plots for the AC group (crosses) of (a) BV/TV vs. measured fracture toughness (KC); and (b) best stepwise regression model (step-3, Table 7) for predicted KC vs. measured KC. Regression lines with their 95% confidence (CI) and prediction (PI) intervals.
Figure 2. Plots for the AC group (crosses) of (a) BV/TV vs. measured fracture toughness (KC); and (b) best stepwise regression model (step-3, Table 7) for predicted KC vs. measured KC. Regression lines with their 95% confidence (CI) and prediction (PI) intervals.
Life 14 00467 g002
Figure 3. Plots for the AL group (diamonds) of (a) BV/TV vs. measured fracture toughness (KC); and (b) predicted KC vs. measured KC. The performance is the same because the best predictor for Kc (AL) is in fact BV/TV. Regression lines with their 95% confidence (CI) and prediction (PI) intervals.
Figure 3. Plots for the AL group (diamonds) of (a) BV/TV vs. measured fracture toughness (KC); and (b) predicted KC vs. measured KC. The performance is the same because the best predictor for Kc (AL) is in fact BV/TV. Regression lines with their 95% confidence (CI) and prediction (PI) intervals.
Life 14 00467 g003
Figure 4. Plots for the entire cohort (OP and OA) and for both AC (crosses), AL (diamonds) in the following: (a) BV/TV vs. measured fracture toughness (KC); and (b) predicted KC vs. measured KC by using the best stepwise regression model (step-4, Table 8) for the entire cohort. Regression lines with their 95% confidence (CI) and prediction (PI) intervals.
Figure 4. Plots for the entire cohort (OP and OA) and for both AC (crosses), AL (diamonds) in the following: (a) BV/TV vs. measured fracture toughness (KC); and (b) predicted KC vs. measured KC by using the best stepwise regression model (step-4, Table 8) for the entire cohort. Regression lines with their 95% confidence (CI) and prediction (PI) intervals.
Life 14 00467 g004
Table 1. Anthropometric and demographic data of OP and OA groups.
Table 1. Anthropometric and demographic data of OP and OA groups.
OPOA
Donors378
Male/Female7/305/3
Number of specimens6019
Age range (years)59–9653–76
Age mean (years)82.3 (SD = 6.8)66 (SD = 7.3)
Weight range (kg)41.3–82.668–108
Weight mean (kg)64.2 (SD = 10.5)84.5 (SD = 12.96)
Height range (m)1.55–1.801.65–1.83
Height mean (m)1.67 (SD = 0.08)1.76 (SD = 0.074)
SD: standard deviation.
Table 2. Average micro-architecture properties for samples loaded as AL and AC. Standard deviation (SD) and p-values denoting significant difference (p < 0.05) between the two groups are included. (Abbreviations for the parameters are explained in the micro-computer tomography section in Materials and Methods.)
Table 2. Average micro-architecture properties for samples loaded as AL and AC. Standard deviation (SD) and p-values denoting significant difference (p < 0.05) between the two groups are included. (Abbreviations for the parameters are explained in the micro-computer tomography section in Materials and Methods.)
All SamplesAL (N = 34)AC (N = 27)
ParameterMeanSDMeanSDMeanSDp-Values
BV/TV0.210.0790.230.090.2020.0730.14
BS/BV14.532.6213.992.4315.002.720.09
BS/TV2.940.573.030.622.860.510.18
TbTh0.140.0290.150.030.140.030.14
TbN1.440.341.500.291.380.380.10
Tbsp0.600.240.540.140.660.290.06
DA2.370.682.550.812.210.490.02
Conn. D2.181.142.201.252.161.040.13
SMI−5.874.65−4.923.73−6.705.240.10
Dmat1.790.081.810.071.770.360.88
Dapp0.390.160.410.170.080.140.09
Table 3. Average micro-architecture properties for the OP and OA groups. Standard deviation (SD) and p-values denoting significant difference (p < 0.05) between the two groups are included. Micro-architecture values from other studies are also provided for comparison (modified from [4]) (Abbreviations for the parameters are explained in the micro-computer tomography section in Materials and Methods).
Table 3. Average micro-architecture properties for the OP and OA groups. Standard deviation (SD) and p-values denoting significant difference (p < 0.05) between the two groups are included. Micro-architecture values from other studies are also provided for comparison (modified from [4]) (Abbreviations for the parameters are explained in the micro-computer tomography section in Materials and Methods).
OPOA
ParameterMeanSDMeanSDp-ValuesRange
in the Literature
References
Kc0.290.140.420.280.02--
BV/TV0.200.050.270.094<0.010.07–0.30[19,20,21,22,23,24,25]
BS/BV (mm−1)14.622.3712.792.250.018.70–22.5[25,26,27]
BS/TV (mm−1)2.860.473.330.75<0.010.59–5.00[28,29]
TbTh (mm−1)0.140.020.160.03<0.010.09–0.25[20,21,22,24,25,30]
TbN (mm−1)1.460.251.630.340.040.76–2.52[22,24,25,29,31]
Tbsp (mm−1)0.570.130.500.160.150.30–1.22[21,22,24,25,31]
DA2.410.622.250.810.161.73–2.00[21,24]
Conn. D1.830.883.360.94<0.011.96–5.62[32]
SMI−6.153.48−1.862.89<0.010.50–2.61[21,24,27,30]
Dmat (g cm−3)1.800.091.800.050.741.40–2.00[9,11,33,34]
Dapp (g cm−3)0.360.090.500.18<0.010.12–0.37[23,25]
Table 4. R2 and p-values for correlations between architectural parameters and fracture toughness (bold typeface for p < 0.05). (Abbreviations for the parameters are explained in the micro-computer tomography section in Materials and Methods).
Table 4. R2 and p-values for correlations between architectural parameters and fracture toughness (bold typeface for p < 0.05). (Abbreviations for the parameters are explained in the micro-computer tomography section in Materials and Methods).
AllALAC
ParameterR2p-ValueR2p-valueR2p-Value
BV/TV0.66<0.010.74<0.010.67<0.01
BS/BV0.25<0.010.25<0.010.35<0.01
BS/TV0.57<0.010.51<0.010.72<0.01
TbTh0.34<0.010.370.0250.39<0.01
TbN0.56<0.010.51<0.010.72<0.01
Tbsp0.40<0.010.49<0.010.41<0.01
DA0.000.740.010.600.100.11
Conn. D0.16<0.010.26<0.010.090.13
Dmat0.110.010.130.040.180.025
Dapp0.64<0.010.73<0.010.65<0.01
Table 5. R2 and p-values for correlations between architectural parameters and fracture toughness with OP and OA groups separated (bold typeface for p < 0.05). (Abbreviations for the parameters are explained in the micro-computer tomography section in Materials and Methods).
Table 5. R2 and p-values for correlations between architectural parameters and fracture toughness with OP and OA groups separated (bold typeface for p < 0.05). (Abbreviations for the parameters are explained in the micro-computer tomography section in Materials and Methods).
OPOA
ParameterR2p-ValueR2p-Value
BV/TV0.58<0.010.69<0.01
BS/BV0.17<0.010.370.02
BS/TV0.52<0.010.56<0.01
TbTh0.21<0.010.400.02
TbN0.44<0.010.71<0.01
Tbsp0.48<0.010.290.05
DA0.010.490.010.72
Conn. D0.030.280.280.05
Dmat0.110.220.230.08
Dapp0.54<0.010.68<0.01
Table 6. R2, Adjusted R2 (bold) and Significance-F multiple linear regressions for all samples as well as a separation of the AL and AC loading groups.
Table 6. R2, Adjusted R2 (bold) and Significance-F multiple linear regressions for all samples as well as a separation of the AL and AC loading groups.
AllALAC
R2Adj-R2Signif-FR2Adj-R2Signif-FR2Adj-R2Signif-F
BV/TV0.6700.656<0.0010.7380.730<0.0010.6740.661<0.001
Stepwise selection0.7590.741<0.0010.7380.730<0.0010.7980.771<0.001
Table 7. Stepwise regression steps for the AC group using BV/TV as the base predictor with the addition of degree of anisotropy (DA), trabecular thickness (TbTh), connectivity density (Conn. D), and trabecular spacing (Tbsp) (Alpha to add 0.05), the final step with all predictors being significant (p < 0.05) is given in bold. (Abbreviations for the parameters are explained in the micro-computer tomography section in Materials and Methods).
Table 7. Stepwise regression steps for the AC group using BV/TV as the base predictor with the addition of degree of anisotropy (DA), trabecular thickness (TbTh), connectivity density (Conn. D), and trabecular spacing (Tbsp) (Alpha to add 0.05), the final step with all predictors being significant (p < 0.05) is given in bold. (Abbreviations for the parameters are explained in the micro-computer tomography section in Materials and Methods).
Step
12345
Constant−0.133−0.436−0.1259.565 × 10−6−0.262
BV/TV2.312.263.493.974.74
T-value7.187.776.216.466.19
p-value<0.001<0.001<0.001<0.001<0.001
DA 0.600.620.500.53
T-value 2.572.932.272.47
p-value 0.0170.0080.0330.022
TbTh −4.1−4.7−5.4
T-value −2.48−2.85−3.26
p-value 0.0210.0090.004
Conn. D −0.040−0.044
T-value −1.65−1.85
p-value 0.1140.079
Tbsp 0.32
T-value 1.59
p-value 0.126
S0.1150.1040.0940.0910.088
R20.6730.7440.7900.8200.840
R2 (adjusted)0.6600.7220.7710.7870.801
Table 8. Stepwise regression steps for the entire cohort using BV/TV as the base predictor with the addition TbN, Tbsp, BS/TV, and DA (Alpha to add 0.15), the final step with all predictors being significant (p < 0.05) is given in bold. (Abbreviations for the parameters are explained in the micro-computer tomography section in Materials and Methods).
Table 8. Stepwise regression steps for the entire cohort using BV/TV as the base predictor with the addition TbN, Tbsp, BS/TV, and DA (Alpha to add 0.15), the final step with all predictors being significant (p < 0.05) is given in bold. (Abbreviations for the parameters are explained in the micro-computer tomography section in Materials and Methods).
Step
12345
Constant−0.1712−0.3650−1.0018−1.0062−1.1267
BV/TV2.271.371.601.811.74
T-value10.833.484.094.714.54
p-value<0.001<0.001<0.001<0.001<0.001
TbN 0.2620.4580.9831.006
T-value 2.643.653.964.10
p-value 0.0110.001<0.001<0.001
Tbsp 0.530.540.55
T-value 2.402.532.62
p-value 0.0200.0140.011
BS/TV −0.28−0.27
T-value −2.42−2.36
p-value 0.0190.022
Degree of Anisotropy 0.028
T-value 1.54
p-value 0.129
S0.1090.1040.0990.0960.095
R20.6690.7050.7330.7580.769
R2 (adjusted)0.6630.6950.7180.7410.747
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Adams, G.J.; Cook, R.B.; Gibson, M.; Zioupos, P. Predicting the Fracture Toughness of Human Cancellous Bone in Fractured Neck of Femur Patients Using Bone Volume and Micro-Architecture. Life 2024, 14, 467. https://doi.org/10.3390/life14040467

AMA Style

Adams GJ, Cook RB, Gibson M, Zioupos P. Predicting the Fracture Toughness of Human Cancellous Bone in Fractured Neck of Femur Patients Using Bone Volume and Micro-Architecture. Life. 2024; 14(4):467. https://doi.org/10.3390/life14040467

Chicago/Turabian Style

Adams, George J., Richard B. Cook, Michael Gibson, and Peter Zioupos. 2024. "Predicting the Fracture Toughness of Human Cancellous Bone in Fractured Neck of Femur Patients Using Bone Volume and Micro-Architecture" Life 14, no. 4: 467. https://doi.org/10.3390/life14040467

APA Style

Adams, G. J., Cook, R. B., Gibson, M., & Zioupos, P. (2024). Predicting the Fracture Toughness of Human Cancellous Bone in Fractured Neck of Femur Patients Using Bone Volume and Micro-Architecture. Life, 14(4), 467. https://doi.org/10.3390/life14040467

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop