Understanding the relationship between interfragmentary movement and callus growth in fracture healing: A novel computational approach

Introduction & Aims Optimising fracture treatments requires a sound understanding of relationships between stability, callus development and healing outcomes. This has been the goal of computational modelling, but discrepancies remain between simulations and experimental results. We compared healing patterns vs fixation stiffness between a novel computational callus growth model and corresponding experimental data. Hypothesis We hypothesised that callus growth is stimulated by diffusible signals, whose production is in turn regulated by mechanical conditions at the fracture site. We proposed that introducing this scheme into computational models would better replicate the observed tissue patterns and the inverse relationship between callus size and fixation stiffness. Method Finite element models of bone healing under stiff and flexible fixation were constructed, based on the parameters of a parallel rat femoral osteotomy study. An iterative procedure was implemented, to simulate the development of callus and its mechanical regulation. Tissue changes were regulated according to published mechano-biological criteria. Predictions of healing patterns were compared between standard models, with a pre-defined domain for callus development, and a novel approach, in which periosteal callus growth is driven by a diffusible signal. Production of this signal was driven by local mechanical conditions. Finally, each model’s predictions were compared to the corresponding histological data. Results Models in which healing progressed within a prescribed callus domain predicted that greater interfragmentary movements would displace early periosteal bone formation further from the fracture. This results from artificially large distortional strains predicted near the fracture edge. While experiments showed increased hard callus size under flexible fixation, this was not reflected in the standard models. Allowing the callus to grow from a thin soft tissue layer, in response to a mechanically stimulated diffusible signal, results in a callus shape and tissue distribution closer to those observed histologically. Importantly, the callus volume increased with increasing interfragmentary movement. Conclusions A novel method to incorporate callus growth into computational models of fracture healing allowed us to successfully capture the relationship between callus size and fixation stability observed in our rat experiments. This approach expands our toolkit for understanding the influence of different fixation strategies on healing outcomes.

Increased intra-cortical porosity reduces bone stiffness and strength in pediatric patients with osteogenesis imperfecta

Osteogenesis imperfecta (OI) is a heritable disease occurring in one out of every 20,000 births. Although it is known that Type I collagen mutation in OI leads to increased bone fragility, the mechanism of this increased susceptibility to fracture is not clear. The aim of this study was to assess the microstructure of cortical bone fragments from patients with osteogenesis imperfecta (OI) using polarized light microscopy, and to correlate microstructural observations with the results of previously performed mechanical compression tests on bone from the same source. Specimens of cortical bone were harvested from the lower limbs of three (3) OI patients at the time of surgery, and were divided into two groups. Group 1 had been subjected to previous micro-mechanical compression testing, while Group 2 had not been subjected to any prior testing. Polarized light microscopy revealed disorganized bone collagen architecture as has been previously observed, as well as a large increase in the areal porosity of the bone compared to typical values for healthy cortical bone, with large (several hundred micron sized), asymmetrical pores. Importantly, the areal porosity of the OI bone samples in Group 1 appears to correlate strongly with their previously measured apparent Young's modulus and compressive strength. Taken together with prior nanoindentation studies on OI bone tissue, the results of this study suggest that increased intra-cortical porosity is responsible for the reduction in macroscopic mechanical properties of OI cortical bone, and therefore that in vivo imaging modalities with resolutions of ~ 100 μm or less could potentially be used to non-invasively assess bone strength in OI patients. Although the number of subjects in this study is small, these results highlight the importance of further studies in OI bone by groups with access to human OI tissue in order to clarify the relationship between increased porosity and reduced macroscopic mechanical integrity.

Toughness and toughening mechanisms of porous thin films

Toughness is the ability of a material to deform plastically and to absorb energy before fracture. The first of its kind, this book covers the most recent developments in the toughening of hard coatings and the methodologies for measuring the toughness of thin films and coatings. The book looks at the present status of toughness for coatings and discusses high-temperature nanocomposite coatings, porous thin films, laser treated surface layers, cracking resistance, indentation techniques, sliding contact fracture, IPN hybrid composites for protection, and adhesion strength.

Genome-wide association study using extreme truncate selection identifies novel genes affecting bone mineral density and fracture risk

Osteoporotic fracture is a major cause of morbidity and mortality worldwide. Low bone mineral density (BMD) is a major predisposing factor to fracture and is known to be highly heritable. Site-, gender-, and age-specific genetic effects on BMD are thought to be significant, but have largely not been considered in the design of genome-wide association studies (GWAS) of BMD to date. We report here a GWAS using a novel study design focusing on women of a specific age (postmenopausal women, age 55-85 years), with either extreme high or low hip BMD (age- and gender-adjusted BMD z-scores of +1.5 to +4.0, n = 1055, or -4.0 to -1.5, n = 900), with replication in cohorts of women drawn from the general population (n = 20,898). The study replicates 21 of 26 known BMD-associated genes. Additionally, we report suggestive association of a further six new genetic associations in or around the genes CLCN7, GALNT3, IBSP, LTBP3, RSPO3, and SOX4, with replication in two independent datasets. A novel mouse model with a loss-of-function mutation in GALNT3 is also reported, which has high bone mass, supporting the involvement of this gene in BMD determination. In addition to identifying further genes associated with BMD, this study confirms the efficiency of extreme-truncate selection designs for quantitative trait association studies. © 2011 Duncan et al.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ~2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2×10-9). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3×10-12, and -0.16 SD per G allele, P = 1.2×10-15, respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3×10-9), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9×10-6 and rs2707466: OR = 1.22, P = 7.2×10-6). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16-/- mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5×10-13<P<5.9×10-4) at both femur and tibia, compared with their wild-type littermates. Natural variation in humans and targeted disruption in mice demonstrate that WNT16 is an important determinant of CBT, BMD, bone strength, and risk of fracture. © 2012 Zheng et al.

Fracture healing via periosteal callus formation requires macrophages for both initiation and progression of early endochondral ossification

The distribution, phenotype, and requirement of macrophages for fracture-associated inflammation and/or early anabolic progression during endochondral callus formation were investigated. A murine femoral fracture model [internally fixed using a flexible plate (MouseFix)] was used to facilitate reproducible fracture reduction. IHC demonstrated that inflammatory macrophages (F4/80+Mac-2+) were localized with initiating chondrification centers and persisted within granulation tissue at the expanding soft callus front. They were also associated with key events during soft-to-hard callus transition. Resident macrophages (F4/80+Mac-2neg), including osteal macrophages, predominated in the maturing hard callus. Macrophage Fas-induced apoptosis transgenic mice were used to induce macrophage depletion in vivo in the femoral fracture model. Callus formation was completely abolished when macrophage depletion was initiated at the time of surgery and was significantly reduced when depletion was delayed to coincide with initiation of early anabolic phase. Treatment initiating 5 days after fracture with the pro-macrophage cytokine colony stimulating factor-1 significantly enhanced soft callus formation. The data support that inflammatory macrophages were required for initiation of fracture repair, whereas both inflammatory and resident macrophages promoted anabolic mechanisms during endochondral callus formation. Overall, macrophages make substantive and prolonged contributions to fracture healing and can be targeted as a therapeutic approach for enhancing repair mechanisms. Thus, macrophages represent a viable target for the development of pro-anabolic fracture treatments with a potentially broad therapeutic window...

Risedronate in adults with osteogenesis imperfecta type I: Increased bone mineral density and decreased bone turnover, but high fracture rate persists

Summary Bisphosphonates can increase bone mineral density (BMD) in children with osteogenesis imperfecta (OI). In this study of adults with OI type I, risedronate increased BMD at lumbar spine (but not total hip) and decreased bone turnover. However, the fracture rate in these patients remained high. Introduction Intravenous bisphosphonates given to children with OI can increase BMD and reduce fracture incidence. Oral and/or intravenous bisphosphonates may have similar effects in adults with OI. We completed an observational study of the effect of risedronate in adults with OI type I. Methods Thirty-two adults (mean age, 39 years) with OI type I were treated with risedronate (total dose, 35 mg weekly) for 24 months. Primary outcome measures were BMD changes at lumbar spine (LS) and total hip (TH). Secondary outcome measures were fracture incidence, bone pain, and change in bone turnover markers (serum procollagen type I aminopropeptide (P1NP) and bone ALP). A meta-analysis of published studies of oral bisphosphonates in adults and children with OI was performed. Results Twenty-seven participants (ten males and seventeen females) completed the study. BMD increased at LS by 3.9% (0.815 vs. 0.846 g/cm 2, p=0.007; mean Z-score, -1.93 vs. -1.58, p=0.002), with no significant change at TH. P1NP fell by 37% (p=0.00041), with no significant change in bone ALP (p=0.15). Bone pain did not change significantly (p=0.6). Fracture incidence remained high, with 25 clinical fractures and 10 major fractures in fourteen participants (0.18 major fractures per person per year), with historical data of 0.12 fractures per person per year. The meta-analysis did not demonstrate a significant difference in fracture incidence in patients with OI treated with oral bisphosphonates. Conclusions Risedronate in adults with OI type I results in modest but significant increases in BMD at LS, and decreased bone turnover. However, this may be insufficient to make a clinically significant difference to fracture incidence.

Genetic determinants of bone density and fracture risk: State of the art and future directions

Context: Osteoporosis is a common, highly heritable condition that causes substantial morbidity and mortality, the etiopathogenesis of which is poorly understood. Genetic studies are making increasingly rapid progress in identifying the genes involved. Evidence Acquisition and Synthesis: In this review, we will summarize the current understanding of the genetics of osteoporosis based on publications from PubMed from the year 1987 onward. Conclusions: Most genes involved in osteoporosis identified to date encode components of known pathways involved in bone synthesis or resorption, but as the field progresses, new pathways are being identified. Only a small proportion of the total genetic variation involved in osteoporosis has been identified, and new approaches will be required to identify most of the remaining genes.

Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture

Bone mineral density (BMD) is the most widely used predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and east Asian ancestry. We tested the top BMD-associated markers for replication in 50,933 independent subjects and for association with risk of low-trauma fracture in 31,016 individuals with a history of fracture (cases) and 102,444 controls. We identified 56 loci (32 new) associated with BMD at genome-wide significance (P < 5 × 10−8). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal stem cell differentiation, endochondral ossification and Wnt signaling pathways. However, we also discovered loci that were localized to genes not known to have a role in bone biology. Fourteen BMD-associated loci were also associated with fracture risk (P < 5 × 10−4, Bonferroni corrected), of which six reached P < 5 × 10−8, including at 18p11.21 (FAM210A), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.

The critical role of water in spider silk and its consequence for protein mechanics

Due to its remarkable mechanical and biological properties, there is considerable interest in understanding, and replicating, spider silk's stress-processing mechanisms and structure-function relationships. Here, we investigate the role of water in the nanoscale mechanics of the different regions in the spider silk fibre, and their relative contributions to stress processing. We propose that the inner core region, rich in spidroin II, retains water due to its inherent disorder, thereby providing a mechanism to dissipate energy as it breaks a sacrificial amide-water bond and gains order under strain, forming a stronger amide-amide bond. The spidroin I-rich outer core is more ordered under ambient conditions and is inherently stiffer and stronger, yet does not on its own provide high toughness. The markedly different interactions of the two proteins with water, and their distribution across the fibre, produce a stiffness differential and provide a balance between stiffness, strength and toughness under ambient conditions. Under wet conditions, this balance is destroyed as the stiff outer core material reverts to the behaviour of the inner core.

A three-dimensional hybrid smoothed finite element method (H-SFEM) for nonlinear solid mechanics problems

This paper presents a novel three-dimensional hybrid smoothed finite element method (H-SFEM) for solid mechanics problems. In 3D H-SFEM, the strain field is assumed to be the weighted average between compatible strains from the finite element method (FEM) and smoothed strains from the node-based smoothed FEM with a parameter α equipped into H-SFEM. By adjusting α, the upper and lower bound solutions in the strain energy norm and eigenfrequencies can always be obtained. The optimized α value in 3D H-SFEM using a tetrahedron mesh possesses a close-to-exact stiffness of the continuous system, and produces ultra-accurate solutions in terms of displacement, strain energy and eigenfrequencies in the linear and nonlinear problems. The novel domain-based selective scheme is proposed leading to a combined selective H-SFEM model that is immune from volumetric locking and hence works well for nearly incompressible materials. The proposed 3D H-SFEM is an innovative and unique numerical method with its distinct features, which has great potential in the successful application for solid mechanics problems.

Prediction of masonry compressive behaviour using a damage mechanics inspired modelling method

Masonry under compression is affected by the properties of its constituents and their interfaces. In spite of extensive investigations of the behaviour of masonry under compression, the information in the literature cannot be regarded as comprehensive due to ongoing inventions of new generation products – for example, polymer modified thin layer mortared masonry and drystack masonry. As comprehensive experimental studies are very expensive, an analytical model inspired by damage mechanics is developed and applied to the prediction of the compressive behaviour of masonry in this paper. The model incorporates a parabolic progressively softening stress-strain curve for the units and a progressively stiffening stress-strain curve until a threshold strain for the combined mortar and the unit-mortar interfaces is reached. The model simulates the mutual constraints imposed by each of these constituents through their respective tensile and compressive behaviour and volumetric changes. The advantage of the model is that it requires only the properties of the constituents and considers masonry as a continuum and computes the average properties of the composite masonry prisms/wallettes; it does not require discretisation of prism or wallette similar to the finite element methods. The capability of the model in capturing the phenomenological behaviour of masonry with appropriate elastic response, stiffness degradation and post peak softening is presented through numerical examples. The fitting of the experimental data to the model parameters is demonstrated through calibration of some selected test data on units and mortar from the literature; the calibrated model is shown to predict the responses of the experimentally determined masonry built using the corresponding units and mortar quite well. Through a series of sensitivity studies, the model is also shown to predict the masonry strength appropriately for changes to the properties of the units and mortar, the mortar joint thickness and the ratio of the height of unit to mortar joint thickness. The unit strength is shown to affect the masonry strength significantly. Although the mortar strength has only a marginal effect, reduction in mortar joint thickness is shown to have a profound effect on the masonry strength. The results obtained from the model are compared with the various provisions in the Australian Masonry Structures Standard AS3700 (2011) and Eurocode 6.

Impact of malnutrition on 12-month mortality following acute hip fracture

Background Studies investigating the relationship between malnutrition and post-discharge mortality following acute hip fracture yield conflicting results. This study aimed to determine whether malnutrition independently predicted 12-month post-fracture mortality after adjusting for clinically relevant covariates. Methods An ethics approved, prospective, consecutive audit was undertaken for all surgically treated hip fracture inpatients admitted to a dedicated orthogeriatric unit (November 2010–October 2011). The 12-month mortality data were obtained by a dual search of the mortality registry and Queensland Health database. Malnutrition was evaluated using the Subjective Global Assessment. Demographic (age, gender, admission residence) and clinical covariates included fracture type, time to surgery, anaesthesia type, type of surgery, post-surgery time to mobilize and post-operative complications (delirium, pulmonary and deep vein thrombosis, cardiac complications, infections). The Charlson Comorbidity Index was retrospectively applied. All diagnoses were confirmed by the treating orthogeriatrician. Results A total of 322 of 346 patients were available for audit. Increased age (P = 0.004), admission from residential care (P < 0.001), Charlson Comorbidity Index (P = 0.007), malnutrition (P < 0.001), time to mobilize >48 h (P < 0.001), delirium (P = 0.003), pulmonary embolism (P = 0.029) and cardiovascular complication (P = 0.04) were associated with 12-month mortality. Logistic regression analysis demonstrated that malnutrition (odds ratio (OR) 2.4 (95% confidence interval (CI) 1.3–4.7, P = 0.007)), in addition to admission from residential care (OR 2.6 (95% CI 1.3–5.3, P = 0.005)) and pulmonary embolism (OR 11.0 (95% CI 1.5–78.7, P = 0.017)), independently predicted 12-month mortality. Conclusions Findings substantiate malnutrition as an independent predictor of 12-month mortality in a representative sample of hip fracture inpatients. Effective strategies to identify and treat malnutrition in hip fracture should be prioritized.

Development of a transient gradient enhanced non local continuum damage mechanics model for masonry

Due to the advent of varied types of masonry systems a comprehensive failure mechanism of masonry essential for the understanding of its behaviour is impossible to be determined from experimental testing. As masonry is predominantly used in wall structures a biaxial stress state dominates its failure mechanism. Biaxial testing will therefore be necessary for each type of masonry, which is expensive and time consuming. A computational method would be advantageous; however masonry is complex to model which requires advanced computational modelling methods. This thesis has formulated a damage mechanics inspired modelling method and has shown that the method effectively determines the failure mechanisms and deformation characteristics of masonry under biaxial states of loading.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 x 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 x 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 x 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.