In vivo morphometric and mechanical characterization of trabecular bone from high resolution magnetic resonance imaging

La osteoporosis es una enfermedad ósea que se manifiesta con una menor densidad ósea y el deterioro de la arquitectura del hueso esponjoso. Ambos factores aumentan la fragilidad ósea y el riesgo de sufrir fracturas óseas, especialmente en mujeres, donde existe una alta prevalencia. El diagnóstico actual de la osteoporosis se basa en la cuantificación de la densidad mineral ósea (DMO) mediante la técnica de absorciometría dual de rayos X (DXA). Sin embargo, la DMO no puede considerarse de manera aislada para la evaluación del riesgo de fractura o los efectos terapéuticos. Existen otros factores, tales como la disposición microestructural de las trabéculas y sus características que es necesario tener en cuenta para determinar la calidad del hueso y evaluar de manera más directa el riesgo de fractura. Los avances técnicos de las modalidades de imagen médica, como la tomografía computarizada multidetector (MDCT), la tomografía computarizada periférica cuantitativa (HR-pQCT) y la resonancia magnética (RM) han permitido la adquisición in vivo con resoluciones espaciales elevadas. La estructura del hueso trabecular puede observarse con un buen detalle empleando estas técnicas. En particular, el uso de los equipos de RM de 3 Teslas (T) ha permitido la adquisición con resoluciones espaciales muy altas. Además, el buen contraste entre hueso y médula que proporcionan las imágenes de RM, así como la utilización de radiaciones no ionizantes sitúan a la RM como una técnica muy adecuada para la caracterización in vivo de hueso trabecular en la enfermedad de la osteoporosis. En la presente tesis se proponen nuevos desarrollos metodológicos para la caracterización morfométrica y mecánica del hueso trabecular en tres dimensiones (3D) y se aplican a adquisiciones de RM de 3T con alta resolución espacial. El análisis morfométrico está compuesto por diferentes algoritmos diseñados para cuantificar la morfología, la complejidad, la topología y los parámetros de anisotropía del tejido trabecular. En cuanto a la caracterización mecánica, se desarrollaron nuevos métodos que permiten la simulación automatizada de la estructura del hueso trabecular en condiciones de compresión y el cálculo del módulo de elasticidad. La metodología desarrollada se ha aplicado a una población de sujetos sanos con el fin de obtener los valores de normalidad del hueso esponjoso. Los algoritmos se han aplicado también a una población de pacientes con osteoporosis con el fin de cuantificar las variaciones de los parámetros en la enfermedad y evaluar las diferencias con los resultados obtenidos en un grupo de sujetos sanos con edad similar.Los desarrollos metodológicos propuestos y las aplicaciones clínicas proporcionan resultados satisfactorios, presentando los parámetros una alta sensibilidad a variaciones de la estructura trabecular principalmente influenciadas por el sexo y el estado de enfermedad. Por otra parte, los métodos presentan elevada reproducibilidad y precisión en la cuantificación de los valores morfométricos y mecánicos. Estos resultados refuerzan el uso de los parámetros presentados como posibles biomarcadores de imagen en la enfermedad de la osteoporosis.Alberich Bayarri, Á. (2010). In vivo morphometric and mechanical characterization of trabecular bone from high resolution magnetic resonance imaging [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/8981Palanci

Rod-like trabeculae extraction from cancellous bone microstructure using topological analysis

The study of the microstructure in cancellous bone is important in osteoporosis research. In this paper, we focus in the extraction of rod-like trabeculae in cancellous hone, whose topological and geometrical properties are known to be strongly related to bone strength. Three dimensional (3D) cancellous bone volume data obtained from micro Computed Tomography (micro-CT) scans is used in our analysis. Because of the complexity of the cancellous bone microstructure, the mere use of existing 3D digital topology method's is unable to yield satisfactory results for trabeculae extraction. We devise several enhancements which are specific to the identification of rod-like trabeculae. Experimental results show that the accuracy of trabeculae extraction is greatly improved. This application study may help clinicians to solve problems in future. © 2005 IEEE.published_or_final_versio

Bone mechanical properties in healthy and diseased states

The mechanical properties of bone are fundamental to the ability of our skeletons to support movement and to provide protection to our vital organs. As such, deterioration in mechanical behavior with aging and/or diseases such as osteoporosis and diabetes can have profound consequences for individuals’ quality of life. This article reviews current knowledge of the basic mechanical behavior of bone at length scales ranging from hundreds of nanometers to tens of centimeters. We present the basic tenets of bone mechanics and connect them to some of the arcs of research that have brought the field to recent advances. We also discuss cortical bone, trabecular bone, and whole bones, as well as multiple aspects of material behavior, including elasticity, yield, fracture, fatigue, and damage. We describe the roles of bone quantity (e.g., density, porosity) and bone quality (e.g., cross-linking, protein composition), along with several avenues of future research.Author manuscrip

Synchrotron Radiation Micro-CT Imaging of Bone Tissue

International audienc

Long-term prediction of three-dimensional bone architecture in simulations of pre-, peri- and post-menopausal microstructural bone remodeling

The mechanical behavior of trabecular bone depends on the internal bone structure. It is generally accepted now that the trabecular bone structure is a result of a load adaptive bone remodeling. The mathematical laws that relate bone remodeling to the local state of stress and strain, however, are still under investigation. The aim of this project was to investigate if changes in the trabecular architecture as observed with age-related bone loss and osteoporosis can be predicted from a computer model that simulates bone resorption after hormone depletion based on realistic models of trabecular microstructure using micro-computed tomography (μCT). A compact desktop μCT providing a nominal isotropic resolution of 14μm was used to measure two groups of seven trabecular bone specimens from pre-menopausal and post-menopausal women respectively. A novel algorithm was developed to simulate age-related bone loss for the specimens in the first group. The algorithm, also referred to as simulated bone atrophy (SIBA), describes a truly three-dimensional approach and is based directly on cellular bone remodeling with an underlying realistic time frame. Bone resorption is controlled by osteoclastic penetration depth and bone formation is governed by the efficiency level of the osteoblasts. The simulation itself describes an iterative process with a cellular remodeling cycle of 197 days. Activation frequency is controllable and can be adjusted for the different phases of pre-, peri- and post-menopause. For our simulations, osteoblastic and osteoclastic activities were in balance until the onset of menopause, set to be at the age of 50 years. In that period, the structure remained almost constant. After the onset of menopause an imbalance in the cell activities was modeled resulting in a net bone loss. The doubling of the activation frequency in the peri-menopausal phase caused a pronounced loss. Using advanced animation tools and quantitative bone morphometry, the changes in bone architecture associated with the bone loss were monitored over an average observation time of 43 years until the age of 80 years. In that time, bone volume density decreased monotonously with the progression of the simulation for all specimens. Right after the onset of menopause, bone was lost fast, where with the progression of age losses slowed down. The structures at the end-point of the simulations were then compared qualitatively and quantitatively to the structures of the post-menopausal group with all morphometric indices being within a narrow margin of error. These results suggest the feasibility of transforming "normal” to "osteopenic” bone on a microstructural level yielding in realistic bone models similar in appearance as well as in structural behavior if compared to a post-menopausal group of wome

Mechanical and morphometric characterization of cancellous bone

[EN] Bone fracture is a social health problem of increasing magnitude because of its prevalence in aged population due to osteoporosis. Bone quality is often characterized by bone mineral density (BMD) measured at cancellous bone regions using dual-energy X-ray absorptiometry (DXA). However, BMD alone cannot predict several cases because not only density is important, but also microstructure plays an important role in cancellous bone strength. The mechanical properties can be used as indicators of bone integrity as a function of age, disease or treatment. Therefore, cancellous bone fracture characterization and its relationship to microstructure has not been completely solved in the literature and is relevant to improve fracture prediction. In this thesis, we aim at characterizing cancellous bone morphometry and mechanical behavior. Morphometry is estimated through the analysis of micro-computed tomography (micro-CT) images of vertebral cancellous bone specimens. With regards to the mechanical behavior, we calculate elastic, yield and failure properties at the apparent and tissue levels. To determine them, we followed different approaches: compression tests, finite element models and micro-CT phantoms. We have developed finite element models that reproduce the elastic and failure response of cancellous bone under compression conditions. We modeled failure as a combination of continuum damage mechanics and the element deletion technique. The numerical models permitted to estimate elastic and failure properties. Failure properties were consistent with results reported in the literature. Specifically, our results revealed that yield strain is relatively constant (0.7 %) over a range of apparent densities, while failure strain presents a wider range of variation. A single strain parameter (equivalent strain) was found as an accurate descriptor of cancellous bone compression failure. Image-based numerical models usually need for the action of a technician to segment the images. Therefore, we studied the sensitivity to variations of the segmentation threshold on the morphometry and the elastic properties of vertebral cancellous bone specimens of different bone volume fractions. The apparent modulus is highly sensitive to the segmentation threshold. We report variations between 45 and 120 % for a ± 15 % threshold variation. Other parameters, such as BS/BV, BS/TV, Tb.Sp, Tb.N, Conn.D and fractal dimension were influenced significantly. Digital image correlation (DIC) was applied to images taken during compression testing to analyze displacement fields at failure and characterize them. Some variables were explored to describe failure and a study is done about how DIC parameters influence the strain field obtained. Facet and step sizes have a relevant effect on the failure strain estimation, and an increment of both parameters reduces the strain estimation up to 40 %. Besides, several parameters combination led to correct failure pattern detection, so values reported in the literature should be referred to the parameters used. Furthermore, we explored if cancellous bone microstructure acts (non-speckle/texture approach) as a proper pattern to calculate displacements using DIC technique. As regards relationships between microstructure and mechanics, single and multiple parameter analysis were performed to assess the morphometric variables that control the explanation of mechanical properties variation. Bone volume fraction (BV/TV), bone surface to volume ratio (BS/BV), mean trabecular thickness (Tb.Th) and fractal dimension (D) presented the best linear correlations to the elastic properties, while both the yield and failure strains did not show correlation to any morphometric parameter. The regressions obtained permit to estimate those mechanical properties that describe the state of a specimen.[ES] Las fracturas óseas constituyen un problema social de salud con magnitud creciente por su prevalencia en la población de edad avanzada debido a la osteoporosis. La calidad del hueso suele caracterizarse mediante la estimación de la densidad mineral ósea (DMO) en regiones de hueso trabecular, utilizando absorciometría de rayos X de energía dual (DXA). No obstante, la DMO por si sola no es capaz de predecir numerosos casos de fractura porque no solo importa la pérdida de densidad, sino que la microestructura también tiene un papel principal en la resistencia del hueso. Las propiedades mecánicas del hueso pueden usarse como indicadores de su integridad en función de la edad, enfermedad o tratamiento. Por lo tanto, la caracterización de la fractura de hueso trabecular y su relación con la microestructura no se ha resuelto de forma completa en la literatura y es relevante para mejorar las predicciones de fractura. En esta tesis, nuestro principal objetivo es caracterizar la morfometría y el comportamiento mecánico del hueso trabecular. Estimamos la morfometría a través del análisis de imágenes obtenidas por micro tomografía computerizada (micro-CT) de muestras de hueso trabecular vertebral de cerdo. Respecto al comportamiento mecánico, calculamos propiedades elásticas, de plasticidad y fractura a escala aparente y de tejido. Para determinar esas propiedades, hemos seguido diferentes procedimientos: ensayos a compresión, modelos de elementos finitos y fantomas de calibración micro-CT. Los modelos de elementos finitos desarrollados reproducen la respuesta elástica y de fallo bajo condiciones de compresión en hueso trabecular, modelando el fallo como combinación de mecánica del daño contínuo y la técnica de eliminación de elementos. Los modelos numéricos desarrollados han permitido estimar propiedades elásticas y de fallo. En concreto, las deformaciones de inicio de fallo estimadas son relativamente constantes para las muestras analizadas (0.7 %), mientras que las deformaciones últimas de fallo presentan un rango de variación mayor. Por otro lado, encontramos que la deformación equivalente es el descriptor más preciso del fallo a compresión del hueso trabecular. Normalmente, los modelos numéricos basados en imágenes suelen necesitar la acción de un técnico para segmentar las imágenes. En este sentido, estudiamos la sensibilidad de la morfometría y la estimación de propiedades elásticas ante variaciones en el umbral de segmentación en muestras con distinta fracción en volumen. Hemos obtenido que la rigidez aparente es muy sensible a cambios en el umbral de segmentación, con variaciones entre 45 y 120 % para una variación de ± 15 % del umbral de segmentación. Otros parámetros, como BS/BV, BS/TV, Tb.Sp, Tb.N, Conn.D y la dimensión fractal se ven afectados significativamente. Por otro lado, hemos aplicado la técnica correlación digital por imagen (DIC) para caracterizar campos de desplazamientos en el fallo a compresión del hueso trabecular, a partir del análisis de imágenes tomadas durante el ensayo de las muestras. Además, estudiamos la influencia de algunos parámetros de la técnica DIC en el campo de deformaciones obtenido. También, hemos explorado la aplicación DIC sin el uso de moteado, utilizando como patrón de reconocimiento la propia microestructura trabecular. En relación al estudio de la influencia de la microestructura en la respuesta mecánica, hemos calculado correlaciones de uno y varios parámetros para analizar qué variables morfométricas explican la variación de las propiedades mecánicas. La fracción en volumen de hueso (BV/TV), la relación entre el área y el volumen de hueso (BS/BV), el espesor trabecular medio (Tb.Th) y la dimensión fractal (D) presentan las mejores correlaciones lineales respecto a las propiedades elásticas, mientras que las deformaciones de inicio de plasticidad y fractura no mostraron correlación con ningún parámetro morfométrico.[CA] Les fractures òssies constitueixen un problema social de salut amb magnitud creixent per la seua prevalença en la població d'edat avançada a causa de l'osteoporosi. La qualitat de l'os sol caracteritzar-se mitjançant l'estimació de la densitat mineral òssia (DMO) en regions d'os trabecular, utilitzant absorciometria de raigs X d'energia dual (DXA). No obstant això, la DMO per si sola no és capaç de predir nombrosos casos de fractura perquè no sols importa la pèrdua de densitat, sinó que la microestructura també té un paper principal en la resistència de l'os. Les propietats mecàniques de l'os poden usar-se com a indicadors de la seua integritat en funció de l'edat, malaltia o tractament. Per tant, la caracterització de la fractura d'os trabecular i la seua relació amb la microestructura no s'ha resolt de manera completa en la literatura i és rellevant per a millorar les prediccions de fractura. En aquesta tesi, el nostre principal objectiu és caracteritzar la morfometria i el comportament mecànic de l'os trabecular. Estimem la morfometria a través de l'anàlisi d'imatges obtingudes per micro tomografia automatitzada (micro-CT) de mostres d'os trabecular vertebral de porc. Respecte al comportament mecànic, calculem propietats elàstiques, de plasticitat i fractura a escala aparent i de teixit. Per a determinar aqueixes propietats, hem seguit diferents procediments: assajos a compressió, models d'elements finits i fantomas de calibratge micro-CT. Hem desenvolupat models d'elements finits que reprodueixen la resposta elàstica i de fallada sota condicions de compressió en os trabecular, modelant la fallada com a combinació de mecànica del dany continu i la tècnica d'eliminació d'elements. Els models numèrics desenvolupats han permés estimar propietats elàstiques i de fallada. Les nostres estimacions respecte a propietats de fallada són consistents amb valors reportats en la literatura. En concret, les deformacions d'inici de fallada estimades són relativament constants per a les mostres analitzades (0.7 %), mentre que les deformacions últimes de fallada presenten un rang de variació major. D'altra banda, trobem que la deformació equivalent és el descriptor més precís de la fallada a compressió de l'os trabecular. Els models numèrics basats en imatges solen necessitar l'acció d'un tècnic per a segmentar les imatges. En aquest sentit, estudiem la sensibilitat de la morfometria i l'estimació de propietats elàstiques davant variacions en el llindar de segmentació en mostres amb diferent fracció en volum. Hem obtingut que la rigidesa aparent és molt sensible a canvis en el llindar de segmentació, amb variacions entre 45 i 120 % per a una variació de ± 15 % del llindar de segmentació. Altres paràmetres, com BS/BV, BS/TV, Tb.Sp, Tb.N, Conn.D i la dimensió fractal es veuen afectats significativament. D'altra banda, hem aplicat la tècnica correlació digital per imatge (DIC) per a caracteritzar camps de desplaçaments en la fallada a compressió de l'os trabecular, a partir de l'anàlisi d'imatges preses durant l'assaig de les mostres. A més, estudiem la influència d'alguns paràmetres de la tècnica DIC en el camp de deformacions obtingut. També, hem explorat l'aplicació DIC sense l'ús de clapejat, utilitzant com a patró de reconeixement la pròpia microestructura trabecular. En relació a l'estudi de la influència de la microestructura en la resposta mecànica, hem calculat correlacions d'un i diversos paràmetres per a analitzar quines variables morfomètriques expliquen la variació de les propietats mecàniques. La fracció en volum d'os (BV/TV), la relació entre l'àrea i el volum d'os (BS/BV), la espessor trabecular mitjà (Tb.th) i la dimensió fractal (D) presenten les millors correlacions lineals respecte a les propietats elàstiques, mentre que les deformacions d'inici de plasticitat i fractura no van mostrar correlació amb cap paràmetre morfomètric.Belda González, R. (2020). Mechanical and morphometric characterization of cancellous bone [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/149376TESI

Diffusion tensor imaging and magnetic resonance spectroscopy assessment of cancellous bone quality in femoral neck of healthy, osteopenic and osteoporotic subjects at 3T: Preliminary experience.

We assessed the potential of diffusion tensor imaging (DTI) in combination with proton magnetic resonance spectroscopy ((1)H-MRS), in cancellous bone quality evaluation of the femoral neck in postmenopausal women. INTRODUCTION: DTI allows for non-invasive microarchitectural characterization of heterogeneous tissue. In this work we hypothesized that DTI parameters mean diffusivity (MD) and fractional anisotropy (FA) of bone marrow water, can provide information about microstructural changes that occur with the development of osteoporosis disease. Because osteoporosis is associated with increased bone marrow fat content, which in principal can alter DTI parameters, the goal of this study was to examine the potential of MD and FA, in combination with bone marrow fat fraction (FF), to discriminate between healthy, osteopenic and osteoporotic subjects, classified according to DXA criteria. MATERIALS AND METHODS: Forty postmenopausal women (mean age, 68.7years; range 52-81years), underwent a Dual-energy X-ray absorptiometry (DXA) examination in femoral neck, to be classified as healthy (n=12), osteopenic (n=14) and osteoporotic (n=14) subjects. (1)H-MRS and DTI (with b value=2500s/mm(2)) of femoral neck were obtained in each subject at 3T. The study protocol was approved by local Ethics Committee. MD, FA, FF and MD/FF, FA/FF were obtained and compared among the three bone-density groups. One-way ANOVA with multiple comparisons Bonferroni test and Pearson correlation analysis were applied. Receiver operating characteristic (ROC) curve analysis was also performed. RESULTS: Reproducibility of DTI measures was satisfactory. CV was approximately 2%-3% for MD and 4%-5% for FA measurements. Moreover, no significant difference was found in both MD and FA measurements between two separate sessions (median 34days apart) comprised of six healthy volunteers. FF was able to discriminate between healthy and osteoporotic subjects only. Conversely MD and FA were able to discriminate healthy from osteopenic and healthy from osteoporotic subjects, but they were not able to discriminate between osteopenic and osteoporotic patients. A significant correlation between MD and FF was observed in healthy group only. A moderate correlation was found between MD and T-score when all groups together are considered. No significant correlation was found between MD and T-score within groups. A significant positive correlation between FA and FF was found in both osteopenic and osteoporotic groups. Vice-versa no correlation between FA and FF was observed in healthy group. A high significant positive correlation was found between FA and T-score in all groups together, in healthy and in osteoporotic groups. MD/FF and FA/FF are characterized by a higher sensitivity and specificity compared to MD and FA in the discrimination between healthy, and osteoporotic subjects. MD/FF vs FA/FF graph extracted from femoral neck, identify all healthy individuals according to DXA results. CONCLUSION: DTI-(1)H-MRS protocol performed in femoral neck seems to be highly sensitive and specific in identifying healthy subjects. A MR exam is more expensive when compared to a DXA investigation. However, even though DXA BMD evaluation has been the accepted standard for osteoporosis diagnosis, DXA result has a low predictive value on patients' risk for future fractures. Thus, new approaches for examining patients at risk for developing osteoporosis would be desirable. Preliminary results showed here suggest that future studies on a larger population based on DTI assessment in the femoral neck, in combination with (1)H-MRS investigations, might allow screening of high-risk populations and the establishment of cut-off values of normality, with potential application of the method to single subjects

Design and testing of additively manufactured lattice structures for musculoskeletal applications

Additive manufacturing (AM) methods present a new frontier in engineering, allowing the fabrication of porous lattice structures with tailored mechanical properties. AM structures can be made using bio-inert metals, creating controlled stiffness biomaterials. As bone formation is strain dependent, these AM biomaterials can be used in implants to optimise the strain in surrounding trabecular bone for peak bone formation. However, the behaviour of AM lattices varies and is subject to manufacturing constraints. The aim of this PhD was to investigate the mechanical behaviour of AM lattices, and maximise the clinical benefits of AM for musculoskeletal applications. Lattice architecture was shown to affect the anisotropy of an AM lattice biomaterial, increasing the stiffness in directions not often tested in the literature. The mechanical and morphological properties of individual struts within powder bed fusion (PBF) lattices were also shown to vary depending on the orientation of the struts to the build direction. The ultimate tensile strength of titanium alloy (Ti6Al4V) struts more than doubled when built at a low angle versus perpendicular to the build platform, and other properties were substantially lower than for the bulk material. Geometric imperfections were found for struts built at low angles. As such, a low stiffness modified stochastic lattice was designed and tested which avoided the problems found with struts built at low angles. The resulting lattice had improved stiffness isotropy and could be used for musculoskeletal applications, tuned to match the mechanical properties in local trabecular bone and enhancing bone formation.Open Acces