3D scaffolds co-seeded with human endothelial progenitor and mesenchymal stem cells: evidence of prevascularisation within 7 days

Blood supply is a critical issue in most tissue engineering approaches for large defect healing. As vessel ingrowth from surrounding tissues is proven to be insufficient, current strategies are focusing on the neo-vascularisation process. In the present study, we developed an in vitro pre-vascularised construct using 3D polyurethane (PU) scaffolds, based on the association of human Endothelial Progenitor Cells (EPC, CD34+ and CD133+) with human Mesenchymal Stem Cells (MSC). We showed the formation of luminal tubular structures in the co-seeded scaffolds as early as day 7 in culture. These tubular structures were proven positive for endothelial markers von Willebrand Factor and PECAM-1. Of special significance in our constructs is the presence of CD146-positive cells, as a part of the neovasculature scaffolding. These cells, coming from the mesenchymal stem cells population (MSC or EPC-depleted MSC), also expressed other markers of pericyte cells (NG2 and αSMA) that are known to play a pivotal function in the stabilisation of newly formed pre-vascular networks. In parallel, in co-cultures, osteogenic differentiation of MSCs occurred earlier when compared to MSCs monocultures, suggesting the close cooperation between the two cell populations. The presence of angiogenic factors (from autologous platelet lysates) in association with osteogenic factors seems to be crucial for both cell populations' cooperation. These results are promising for future clinical applications, as all components (cells, growth factors) can be prepared in an autologous way.

Human intervertebral disc stiffness correlates better with the Otsu threshold computed from axial T2 map of its posterior annulus fibrosus than with clinical classifications

Degeneration of the intervertebral disc, sometimes associated with low back pain and abnormal spinal motions, represents a major health issue with high costs. A non-invasive degeneration assessment via qualitative or quantitative MRI (magnetic resonance imaging) is possible, yet, no relation between mechanical properties and T2 maps of the intervertebral disc (IVD) has been considered, albeit T2 relaxation time values quantify the degree of degeneration. Therefore, MRI scans and mechanical tests were performed on 14 human lumbar intervertebral segments freed from posterior elements and all soft tissues excluding the IVD. Degeneration was evaluated in each specimen using morphological criteria, qualitative T2 weighted images and quantitative axial T2 map data and stiffness was calculated from the load-deflection curves of in vitro compression, torsion, lateral bending and flexion/extension tests. In addition to mean T2, the OTSU threshold of T2 (TOTSU), a robust and automatic histogram-based method that computes the optimal threshold maximizing the distinction of two classes of values, was calculated for anterior, posterior, left and right regions of each annulus fibrosus (AF). While mean T2 and degeneration schemes were not related to the IVDs' mechanical properties, TOTSU computed in the posterior AF correlated significantly with those classifications as well as with all stiffness values. TOTSU should therefore be included in future degeneration grading schemes.

Tissue properties of the human vertebral body sub-structures evaluated by means of microindentation

Purpose The better understanding of vertebral mechanical properties can help to improve the diagnosis of vertebral fractures. As the bone mechanical competence depends not only from bone mineral density (BMD) but also from bone quality, the goal of the present study was to investigate the anisotropic indentation moduli of the different sub-structures of the healthy human vertebral body and spondylophytes by means of microindentation. Methods Six human vertebral bodies and five osteophytes (spondylophytes) were collected and prepared for microindentation test. In particular, indentations were performed on bone structural units of the cortical shell (along axial, circumferential and radial directions), of the endplates (along the anterio-posterior and lateral directions), of the trabecular bone (along the axial and transverse directions) and of the spondylophytes (along the axial direction). A total of 3164 indentations down to a maximum depth of 2.5 µm were performed and the indentation modulus was computed for each measurement. Results The cortical shell showed an orthotropic behavior (indentation modulus, Ei, higher if measured along the axial direction, 14.6±2.8 GPa, compared to the circumferential one, 12.3±3.5 GPa, and radial one, 8.3±3.1 GPa). Moreover, the cortical endplates (similar Ei along the antero-posterior, 13.0±2.9 GPa, and along the lateral, 12.0±3.0 GPa, directions) and the trabecular bone (Ei= 13.7±3.4 GPa along the axial direction versus Ei=10.9±3.7 GPa along the transverse one) showed transversal isotropy behavior. Furthermore, the spondylophytes showed the lower mechanical properties measured along the axial direction (Ei=10.5±3.3 GPa). Conclusions The original results presented in this study improve our understanding of vertebral biomechanics and can be helpful to define the material properties of the vertebral substructures in computational models such as FE analysis.

DXA predictions of human femoral mechanical properties depend on the load configuration

The aim of this study was to evaluate the ability of dual energy X-rays absorptiometry (DXA) areal bone mineral density (aBMD) measured in different regions of the proximal part of the human femur for predicting the mechanical properties of matched proximal femora tested in two different loading configurations. 36 pairs of fresh frozen femora were DXA scanned and tested until failure in two loading configurations: a fall on the side or a one-legged standing. The ability of the DXA output from four different regions of the proximal femur in predicting the femoral mechanical properties was measured and compared for the two loading scenarios. The femoral neck DXA BMD was best correlated to the femoral ultimate force for both configurations and predicted significantly better femoral failure load (R2=0.80 vs. R2=0.66, P<0.05) when simulating a side than when simulating a standing configuration. Conversely, the work to failure was predicted similarly for both loading configurations (R2=0.54 vs. R2=0.53, P>0.05). Therefore, neck BMD should be considered as one of the key factors for discriminating femoral fracture risk in vivo. Moreover, the better predictive ability of neck BMD for femoral strength if tested in a fall compared to a one-legged stance configuration suggests that DXA's clinical relevance may not be as high for spontaneous femoral fractures than for fractures associated to a fall.

Finite element based estimation of contact areas and pressures of the human scaphoid in various functional positions of the hand

The scaphoid is the most frequently fractured carpal bone. When investigating fixation stability, which may influence healing, knowledge of forces and moments acting on the scaphoid is essential. The aim of this study was to evaluate cartilage contact forces acting on the intact scaphoid in various functional wrist positions using finite element modeling. A novel methodology was utilized as an attempt to overcome some limitations of earlier studies, namely, relatively coarse imaging resolution to assess geometry, assumption of idealized cartilage thicknesses and neglected cartilage pre-stresses in the unloaded joint. Carpal bone positions and articular cartilage geometry were obtained independently by means of high resolution CT imaging and incorporated into finite element (FE) models of the human wrist in eight functional positions. Displacement driven FE analyses were used to resolve inter-penetration of cartilage layers, and provided contact areas, forces and pressure distribution for the scaphoid bone. The results were in the range reported by previous studies. Novel findings of this study were: (i) cartilage thickness was found to be heterogeneous for each bone and vary considerably between carpal bones; (ii) this heterogeneity largely influenced the FE results and (iii) the forces acting on the scaphoid in the unloaded wrist were found to be significant. As major limitations, accuracy of the method was found to be relatively low, and the results could not be compared to independent experiments. The obtained results will be used in a following study to evaluate existing and recently developed screws used to fix scaphoid fractures.

A nonlinear QCT-based finite element model validation study for the human femur tested in two configurations in vitro

Purpose Femoral fracture is a common medical problem in osteoporotic individuals. Bone mineral density (BMD) is the gold standard measure to evaluate fracture risk in vivo. Quantitative computed tomography (QCT)-based homogenized voxel finite element (hvFE) models have been proved to be more accurate predictors of femoral strength than BMD by adding geometrical and material properties. The aim of this study was to evaluate the ability of hvFE models in predicting femoral stiffness, strength and failure location for a large number of pairs of human femora tested in two different loading scenarios. Methods Thirty-six pairs of femora were scanned with QCT and total proximal BMD and BMC were evaluated. For each pair, one femur was positioned in one-legged stance configuration (STANCE) and the other in a sideways configuration (SIDE). Nonlinear hvFE models were generated from QCT images by reproducing the same loading configurations imposed in the experiments. For experiments and models, the structural properties (stiffness and ultimate load), the failure location and the motion of the femoral head were computed and compared. Results In both configurations, hvFE models predicted both stiffness (R2=0.82 for STANCE and R2=0.74 for SIDE) and femoral ultimate load (R2=0.80 for STANCE and R2=0.85 for SIDE) better than BMD and BMC. Moreover, the models predicted qualitatively well the failure location (66% of cases) and the motion of the femoral head. Conclusions The subject specific QCT-based nonlinear hvFE model cannot only predict femoral apparent mechanical properties better than densitometric measures, but can additionally provide useful qualitative information about failure location.

Morphology–elasticity relationships using decreasing fabric information of human trabecular bone from three major anatomical locations

With improving clinical CT scanning technology, the accuracy of CT-based finite element (FE) models of the human skeleton may be ameliorated by an enhanced description of apparent level bone mechanical properties. Micro-finite element (μFE) modeling can be used to study the apparent elastic behavior of human cancellous bone. In this study, samples from the femur, radius and vertebral body were investigated to evaluate the predictive power of morphology–elasticity relationships and to compare them across different anatomical regions. μFE models of 701 trabecular bone cubes with a side length of 5.3 mm were analyzed using kinematic boundary conditions. Based on the FE results, four morphology–elasticity models using bone volume fraction as well as full, limited or no fabric information were calibrated for each anatomical region. The 5 parameter Zysset–Curnier model using full fabric information showed excellent predictive power with coefficients of determination ( r2adj ) of 0.98, 0.95 and 0.94 of the femur, radius and vertebra data, respectively, with mean total norm errors between 14 and 20%. A constant orthotropy model and a constant transverse isotropy model, where the elastic anisotropy is defined by the model parameters, yielded coefficients of determination between 0.90 and 0.98 with total norm errors between 16 and 25%. Neglecting fabric information and using an isotropic model led to r2adj between 0.73 and 0.92 with total norm errors between 38 and 49%. A comparison of the model regressions revealed minor but significant (p<0.01) differences for the fabric–elasticity model parameters calibrated for the different anatomical regions. The proposed models and identified parameters can be used in future studies to compute the apparent elastic properties of human cancellous bone for homogenized FE models.

Impact of dams, dam removal and dam-related river engineering structures on sediment connectivity and channel morphology of the Fugnitz and the Kaja Rivers

In terms of changing flow and sediment regimes of rivers, dams are often regarded as the most dominant form of human impact on fluvial systems. Dams can decrease the flux of water and sediments leading to channel changes such as upstream aggradation and downstream degradation. The opposite effects occur when dams are removed. Channel degradation often requires further intervention in terms of river bed and bank protection works. The situation evolves more complex in river systems that are impacted by a series of dams due to feedback processes between the different system compartments. A number of studies have recently investigated geomorphic systems using connectivity approaches to improve the understanding of geomorphic system response to change. This paper presents a case study investigating the impact of dam construction, dam removal and dam-related river bed and bank protection measures on the sediment connectivity and channel morphology of the Fugnitz and the Kaja Rivers using a combination of DEM analyses, field surveys and landscape evolution modelling. For both river systems the results revealed low sediment connectivity accompanied by a fine river bed sediment facies in river sections upstream of active dams and of removed dams with protection measures. Contrarily, high sediment connectivity which was accompanied by a coarse river bed sediment facies was observed in river sections either located downstream of active dams or of removed dams with upstream protection. In terms of channel changes, significant channel degradation was examined at locations downstream of active dams and of removed dams. Channel bed and bank protection measures prevent erosion and channel slope recovery after dam removal. Landscape evolution modeling revealed a complex geomorphic response to dam construction and dam removal as sediment output rates and therefore geomorphic processes have been shown to act in a non-linear manner. These insights are deemed to have major implications for river management and conservation, as quality and state of riverine habitats are determined by channel morphology and river bed sediment composition.

Cell and matrix morphology in articular cartilage from adult human knee and ankle joints suggests depth-associated adaptations to biomechanical and anatomical roles

OBJECTIVE Marked differences exist between human knee and ankle joints regarding risks and progression of osteoarthritis (OA). Pathomechanisms of degenerative joint disease may therefore differ in these joints, due to differences in tissue structure and function. Focussing on structural issues which are design goals for tissue engineering, we compared cell and matrix morphologies in different anatomical sites of adult human knee and ankle joints. METHODS Osteochondral explants were acquired from knee and ankle joints of deceased persons aged 20 to 40 years and analyzed for cell, matrix and tissue morphology using confocal and electron microscopy and unbiased stereological methods. Variations associated with joint (knee versus ankle) and biomechanical role (convex versus concave articular surfaces) were identified by 2-way analysis of variance and post-hoc analysis. RESULTS Knee cartilage exhibited higher cell densities in the superficial zone than ankle cartilage. In the transitional zone, higher cell densities were observed in association with convex versus concave articular surfaces, without significant differences between knee and ankle cartilage. Highly uniform cell and matrix morphologies were evident throughout the radial zone in the knee and ankle, regardless of tissue biomechanical role. Throughout the knee and ankle cartilage sampled, chondron density was remarkably constant at approximately 4.2×10(6) chondrons/cm(3). CONCLUSION Variation of cartilage cell and matrix morphologies with changing joint and biomechanical environments suggests that tissue structural adaptations are performed primarily by the superficial and transitional zones. Data may aid the development of site-specific cartilage tissue engineering, and help identify conditions where OA is likely to occur.

Finite Element Based Nonlinear Normalization of Human Lumbar Intervertebral Disc Stiffness to Account for Its Morphology

Disc degeneration, usually associated with low back pain and changes of intervertebral stiffness, represents a major health issue. As the intervertebral disc (IVD) morphology influences its stiffness, the link between mechanical properties and degenerative grade is partially lost without an efficient normalization of the stiffness with respect to the morphology. Moreover, although the behavior of soft tissues is highly nonlinear, only linear normalization protocols have been defined so far for the disc stiffness. Thus, the aim of this work is to propose a nonlinear normalization based on finite elements (FE) simulations and evaluate its impact on the stiffness of human anatomical specimens of lumbar IVD. First, a parameter study involving simulations of biomechanical tests (compression, flexion/extension, bilateral torsion and bending) on 20 FE models of IVDs with various dimensions was carried out to evaluate the effect of the disc's geometry on its compliance and establish stiffness/morphology relations necessary to the nonlinear normalization. The computed stiffness was then normalized by height (H), cross-sectional area (CSA), polar moment of inertia (J) or moments of inertia (Ixx, Iyy) to quantify the effect of both linear and nonlinear normalizations. In the second part of the study, T1-weighted MRI images were acquired to determine H, CSA, J, Ixx and Iyy of 14 human lumbar IVDs. Based on the measured morphology and pre-established relation with stiffness, linear and nonlinear normalization routines were then applied to the compliance of the specimens for each quasi-static biomechanical test. The variability of the stiffness prior to and after normalization was assessed via coefficient of variation (CV). The FE study confirmed that larger and thinner IVDs were stiffer while the normalization strongly attenuated the effect of the disc geometry on its stiffness. Yet, notwithstanding the results of the FE study, the experimental stiffness showed consistently higher CV after normalization. Assuming that geometry and material properties affect the mechanical response, they can also compensate for one another. Therefore, the larger CV after normalization can be interpreted as a strong variability of the material properties, previously hidden by the geometry's own influence. In conclusion, a new normalization protocol for the intervertebral disc stiffness in compression, flexion, extension, bilateral torsion and bending was proposed, with the possible use of MRI and FE to acquire the discs' anatomy and determine the nonlinear relations between stiffness and morphology. Such protocol may be useful to relate the disc's mechanical properties to its degree of degeneration.

Experimental Validation of Finite Element Analysis of Human Vertebral Collapse Under Large Compressive Strains

Osteoporosis-related vertebral fractures represent a major health problem in elderly populations. Such fractures can often only be diagnosed after a substantial deformation history of the vertebral body. Therefore, it remains a challenge for clinicians to distinguish between stable and progressive potentially harmful fractures. Accordingly, novel criteria for selection of the appropriate conservative or surgical treatment are urgently needed. Computer tomography-based finite element analysis is an increasingly accepted method to predict the quasi-static vertebral strength and to follow up this small strain property longitudinally in time. A recent development in constitutive modeling allows us to simulate strain localization and densification in trabecular bone under large compressive strains without mesh dependence. The aim of this work was to validate this recently developed constitutive model of trabecular bone for the prediction of strain localization and densification in the human vertebral body subjected to large compressive deformation. A custom-made stepwise loading device mounted in a high resolution peripheral computer tomography system was used to describe the progressive collapse of 13 human vertebrae under axial compression. Continuum finite element analyses of the 13 compression tests were realized and the zones of high volumetric strain were compared with the experiments. A fair qualitative correspondence of the strain localization zone between the experiment and finite element analysis was achieved in 9 out of 13 tests and significant correlations of the volumetric strains were obtained throughout the range of applied axial compression. Interestingly, the stepwise propagating localization zones in trabecular bone converged to the buckling locations in the cortical shell. While the adopted continuum finite element approach still suffers from several limitations, these encouraging preliminary results towardsthe prediction of extended vertebral collapse may help in assessing fracture stability in future work.

Image-based vs. mesh-based statistical appearance models of the human femur: Implications for finite element simulations.

Statistical appearance models have recently been introduced in bone mechanics to investigate bone geometry and mechanical properties in population studies. The establishment of accurate anatomical correspondences is a critical aspect for the construction of reliable models. Depending on the representation of a bone as an image or a mesh, correspondences are detected using image registration or mesh morphing. The objective of this study was to compare image-based and mesh-based statistical appearance models of the femur for finite element (FE) simulations. To this aim, (i) we compared correspondence detection methods on bone surface and in bone volume; (ii) we created an image-based and a mesh-based statistical appearance models from 130 images, which we validated using compactness, representation and generalization, and we analyzed the FE results on 50 recreated bones vs. original bones; (iii) we created 1000 new instances, and we compared the quality of the FE meshes. Results showed that the image-based approach was more accurate in volume correspondence detection and quality of FE meshes, whereas the mesh-based approach was more accurate for surface correspondence detection and model compactness. Based on our results, we recommend the use of image-based statistical appearance models for FE simulations of the femur.

Image-based biomechanical assessment of vertebral body and intervertebral disc in the human lumbar spine

Life expectancy continuously increases but our society faces age-related conditions. Among musculoskeletal diseases, osteoporosis associated with risk of vertebral fracture and degenerative intervertebral disc (IVD) are painful pathologies responsible for tremendous healthcare costs. Hence, reliable diagnostic tools are necessary to plan a treatment or follow up its efficacy. Yet, radiographic and MRI techniques, respectively clinical standards for evaluation of bone strength and IVD degeneration, are unspecific and not objective. Increasingly used in biomedical engineering, CT-based finite element (FE) models constitute the state-of-art for vertebral strength prediction. However, as non-invasive biomechanical evaluation and personalised FE models of the IVD are not available, rigid boundary conditions (BCs) are applied on the FE models to avoid uncertainties of disc degeneration that might bias the predictions. Moreover, considering the impact of low back pain, the biomechanical status of the IVD is needed as a criterion for early disc degeneration. Thus, the first FE study focuses on two rigid BCs applied on the vertebral bodies during compression test of cadaver vertebral bodies, vertebral sections and PMMA embedding. The second FE study highlights the large influence of the intervertebral disc’s compliance on the vertebral strength, damage distribution and its initiation. The third study introduces a new protocol for normalisation of the IVD stiffness in compression, torsion and bending using MRI-based data to account for its morphology. In the last study, a new criterion (Otsu threshold) for disc degeneration based on quantitative MRI data (axial T2 map) is proposed. The results show that vertebral strength and damage distribution computed with rigid BCs are identical. Yet, large discrepancies in strength and damage localisation were observed when the vertebral bodies were loaded via IVDs. The normalisation protocol attenuated the effect of geometry on the IVD stiffnesses without complete suppression. Finally, the Otsu threshold computed in the posterior part of annulus fibrosus was related to the disc biomechanics and meet objectivity and simplicity required for a clinical application. In conclusion, the stiffness normalisation protocol necessary for consistent IVD comparisons and the relation found between degeneration, mechanical response of the IVD and Otsu threshold lead the way for non-invasive evaluation biomechanical status of the IVD. As the FE prediction of vertebral strength is largely influenced by the IVD conditions, this data could also improve the future FE models of osteoporotic vertebra.