Computational and experimental methodology for site-matched investigations of the influence of mineral mass fraction and collagen orientation on the axial indentation modulus of lamellar bone

Relationships between mineralization, collagen orientation and indentation modulus were investigated in bone structural units from the mid-shaft of human femora using a site-matched design. Mineral mass fraction, collagen fibril angle and indentation moduli were measured in registered anatomical sites using backscattered electron imaging, polarized light microscopy and nano-indentation, respectively. Theoretical indentation moduli were calculated with a homogenization model from the quantified mineral densities and mean collagen fibril orientations. The average indentation moduli predicted based on local mineralization and collagen fibers arrangement were not significantly different from the average measured experimentally with nanoindentation (p=0.9). Surprisingly, no substantial correlation of the measured indentation moduli with tissue mineralization and/or collagen fiber arrangement was found. Nano-porosity, micro-damage, collagen cross-links, non-collagenous proteins or other parameters affect the indentation measurements. Additional testing/simulation methods need to be considered to properly understand the variability of indentation moduli, beyond the mineralization and collagen arrangement in bone structural units.

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.

Early changes in biochemical markers of bone formation during teriparatide therapy correlate with improvements in vertebral strength in men with glucocorticoid-induced osteoporosis

Summary Changes of the bone formation marker PINP correlated positively with improvements in vertebral strength in men with glucocorticoid-induced osteoporosis (GIO) who received 18-month treatment with teriparatide, but not with risedronate. These results support the use of PINP as a surrogate marker of bone strength in GIO patients treated with teriparatide. Introduction To investigate the correlations between biochemical markers of bone turnover and vertebral strength estimated by finite element analysis (FEA) in men with GIO. Methods A total of 92 men with GIO were included in an 18-month, randomized, open-label trial of teriparatide (20 μg/day, n = 45) and risedronate (35 mg/week, n = 47). High-resolution quantitative computed tomography images of the 12th thoracic vertebra obtained at baseline, 6 and 18 months were converted into digital nonlinear FE models and subjected to anterior bending, axial compression and torsion. Stiffness and strength were computed for each model and loading mode. Serum biochemical markers of bone formation (amino-terminal-propeptide of type I collagen [PINP]) and bone resorption (type I collagen cross-linked C-telopeptide degradation fragments [CTx]) were measured at baseline, 3 months, 6 months and 18 months. A mixed-model of repeated measures analysed changes from baseline and between-group differences. Spearman correlations assessed the relationship between changes from baseline of bone markers with FEA variables. Results PINP and CTx levels increased in the teriparatide group and decreased in the risedronate group. FEA-derived parameters increased in both groups, but were significantly higher at 18 months in the teriparatide group. Significant positive correlations were found between changes from baseline of PINP at 3, 6 and 18 months with changes in FE strength in the teriparatide-treated group, but not in the risedronate group. Conclusions Positive correlations between changes in a biochemical marker of bone formation and improvement of biomechanical properties support the use of PINP as a surrogate marker of bone strength in teriparatide-treated GIO patients.

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.

High resolution quantitative computed tomography-based assessment of trabecular microstructure and strength estimates by finite-element analysis of the spine, but not DXA, reflects vertebral fracture status in men with glucocorticoid-induced osteoporosis

High-resolution quantitative computed tomography (HRQCT)-based analysis of spinal bone density and microstructure, finite element analysis (FEA), and DXA were used to investigate the vertebral bone status of men with glucocorticoid-induced osteoporosis (GIO). DXA of L1–L3 and total hip, QCT of L1–L3, and HRQCT of T12 were available for 73 men (54.6±14.0years) with GIO. Prevalent vertebral fracture status was evaluated on radiographs using a semi-quantitative (SQ) score (normal=0 to severe fracture=3), and the spinal deformity index (SDI) score (sum of SQ scores of T4 to L4 vertebrae). Thirty-one (42.4%) subjects had prevalent vertebral fractures. Cortical BMD (Ct.BMD) and thickness (Ct.Th), trabecular BMD (Tb.BMD), apparent trabecular bone volume fraction (app.BV/TV), and apparent trabecular separation (app.Tb.Sp) were analyzed by HRQCT. Stiffness and strength of T12 were computed by HRQCT-based nonlinear FEA for axial compression, anterior bending and axial torsion. In logistic regressions adjusted for age, glucocorticoid dose and osteoporosis treatment, Tb.BMD was most closely associated with vertebral fracture status (standardized odds ratio [sOR]: Tb.BMD T12: 4.05 [95% CI: 1.8–9.0], Tb.BMD L1–L3: 3.95 [1.8–8.9]). Strength divided by cross-sectional area for axial compression showed the most significant association with spine fracture status among FEA variables (2.56 [1.29–5.07]). SDI was best predicted by a microstructural model using Ct.Th and app.Tb.Sp (r2=0.57, p<0.001). Spinal or hip DXA measurements did not show significant associations with fracture status or severity. In this cross-sectional study of males with GIO, QCT, HRQCT-based measurements and FEA variables were superior to DXA in discriminating between patients of differing prevalent vertebral fracture status. A microstructural model combining aspects of cortical and trabecular bone reflected fracture severity most accurately.

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.

Space-qualified laser system for the BepiColombo Laser Altimeter

The space-qualified design of a miniaturized laser for pulsed operation at a wavelength of 1064 nm and at repetition rates up to 10 Hz is presented. This laser consists of a pair of diode-laser pumped, actively q-switched Nd:YAG rod oscillators hermetically sealed and encapsulated in an environment of dry synthetic air. The system delivers at least 300 million laser pulses with 50 mJ energy and 5 ns pulse width (FWHM). It will be launched in 2017 aboard European Space Agency’s Mercury Planetary Orbiter as part of the BepiColombo Laser Altimeter, which, after a 6-years cruise, will start recording topographic data from orbital altitudes between 400 and 1500 km above Mercury’s surface.

Analysis of polygonal cracking patterns in chloride-bearing terrains on Mars: Indicators of ancient playa settings

The ancient southern highlands on Mars (~3.5 Gyr old) contain > 600 regions that display spectral evidence in the infrared for the presence of chloride-bearing materials. Many of these locations were previously reported to display polygonal cracking patterns. We studied more than 80 of the chloride-bearing terrains using high-resolution (0.25-0.5 m/pixel) images, as well as near-infrared spectral data, to characterize the surface textures and the associated cracking patterns and mineralogies. Our study indicates that ~75% of the studied locations display polygonal cracks that resemble desiccation cracks, while some resemble salt expansion/thrust polygons. Furthermore, we detect, spectrally, the presence of smectites in association with ~30% of the studied fractured terrains. We note that smectites are a special class of swelling clay minerals that can induce formation of large desiccation cracks. As such, we suggest that the cracking patterns are indicative of the presence of smectite phyllosilicates even in the absence of spectral confirmation. Our results suggest that many chloride-bearing terrains have a lacustrine origin and a geologic setting similar to playas on Earth. Such locations would have contained ephemeral lakes that may have undergone repeated cycles of desiccation and recharging by a near-surface fluctuating water table in order to account for the salt-phyllosilicates associations. These results have notable implications for the ancient hydrology of Mars. We propose that the morphologies and sizes of the polygonal cracks can be used as paleoenvironmental, as well as lithological, indicators that could be helpful in planning future missions.

Morphology and evolution of the ejecta of Hale crater in Argyre basin, Mars: Results from high resolution mapping

We use various data sets, including images from the High Resolution Imaging Science Experiment camera (HiRISE), to examine the ejecta of the generally fresh-looking Hale crater that occurs in the rugged mountain terrain of Nereidum Montes in the northern rim materials of the Argyre impact structure on Mars. Our investigation reveals that the distal parts of the Hale crater ejecta and other basin deposits behave like viscous flows, which we attribute to the secondary flow of ejecta mixed with water–ice-rich basin materials. Consistent with water-enrichment of the basin materials, our mapping further reveals occasionally deformed surfaces, including highly conspicuous features such as mounds and fractured plateaus that we interpret to be a result of periglacial modification, subsequent (including possibly present-day) to the transient localized melting and fluvial erosion caused by Hale-impact-generated heating. In particular, our morphometric analysis of a well-defined valley system west of Hale crater suggests that it may have been formed through hydrologic/glacial activity prior to the Hale impact, with additional modification resulting from the impact and subsequent geologic and hydrologic phenomena including glacial and periglacial activity.

Empirical chemosensitivity testing in a spheroid model of ovarian cancer using a microfluidics-based multiplex platform.

The use of biomarkers to infer drug response in patients is being actively pursued, yet significant challenges with this approach, including the complicated interconnection of pathways, have limited its application. Direct empirical testing of tumor sensitivity would arguably provide a more reliable predictive value, although it has garnered little attention largely due to the technical difficulties associated with this approach. We hypothesize that the application of recently developed microtechnologies, coupled to more complex 3-dimensional cell cultures, could provide a model to address some of these issues. As a proof of concept, we developed a microfluidic device where spheroids of the serous epithelial ovarian cancer cell line TOV112D are entrapped and assayed for their chemoresponse to carboplatin and paclitaxel, two therapeutic agents routinely used for the treatment of ovarian cancer. In order to index the chemoresponse, we analyzed the spatiotemporal evolution of the mortality fraction, as judged by vital dyes and confocal microscopy, within spheroids subjected to different drug concentrations and treatment durations inside the microfluidic device. To reflect microenvironment effects, we tested the effect of exogenous extracellular matrix and serum supplementation during spheroid formation on their chemotherapeutic response. Spheroids displayed augmented chemoresistance in comparison to monolayer culturing. This resistance was further increased by the simultaneous presence of both extracellular matrix and high serum concentration during spheroid formation. Following exposure to chemotherapeutics, cell death profiles were not uniform throughout the spheroid. The highest cell death fraction was found at the center of the spheroid and the lowest at the periphery. Collectively, the results demonstrate the validity of the approach, and provide the basis for further investigation of chemotherapeutic responses in ovarian cancer using microfluidics technology. In the future, such microdevices could provide the framework to assay drug sensitivity in a timeframe suitable for clinical decision making.

Theoretical models of planetary system formation: mass vs. semi-major axis

Context. Planet formation models have been developed during the past years to try to reproduce what has been observed of both the solar system and the extrasolar planets. Some of these models have partially succeeded, but they focus on massive planets and, for the sake of simplicity, exclude planets belonging to planetary systems. However, more and more planets are now found in planetary systems. This tendency, which is a result of radial velocity, transit, and direct imaging surveys, seems to be even more pronounced for low-mass planets. These new observations require improving planet formation models, including new physics, and considering the formation of systems. Aims: In a recent series of papers, we have presented some improvements in the physics of our models, focussing in particular on the internal structure of forming planets, and on the computation of the excitation state of planetesimals and their resulting accretion rate. In this paper, we focus on the concurrent effect of the formation of more than one planet in the same protoplanetary disc and show the effect, in terms of architecture and composition of this multiplicity. Methods: We used an N-body calculation including collision detection to compute the orbital evolution of a planetary system. Moreover, we describe the effect of competition for accretion of gas and solids, as well as the effect of gravitational interactions between planets. Results: We show that the masses and semi-major axes of planets are modified by both the effect of competition and gravitational interactions. We also present the effect of the assumed number of forming planets in the same system (a free parameter of the model), as well as the effect of the inclination and eccentricity damping. We find that the fraction of ejected planets increases from nearly 0 to 8% as we change the number of embryos we seed the system with from 2 to 20 planetary embryos. Moreover, our calculations show that, when considering planets more massive than ~5 M⊕, simulations with 10 or 20 planetary embryos statistically give the same results in terms of mass function and period distribution.

Highly accurate isotope composition measurements by a miniature laser ablation mass spectrometer designed for in situ investigations on planetary surface

An experimental procedure for precise and accurate measurements of isotope abundances by a miniature laser ablation mass spectrometer for space research is described. The measurements were conducted on different untreated NIST standards and galena samples by applying pulsed UV laser radiation (266 nm, 3 ns and 20 Hz) for ablation, atomisation, and ionisation of the sample material. Mass spectra of released ions are measured by a reflectron-type time-of-flight mass analyser. A computer controlled performance optimiser was used to operate the system at maximum ion transmission and mass resolution. At optimal experimental conditions, the best relative accuracy and precision achieved for Pb isotope compositions are at the per mill level and were obtained in a range of applied laser irradiances and a defined number of accumulated spectra. A similar relative accuracy and precision was achieved in the study of Pb isotope compositions in terrestrial galena samples. The results for the galena samples are similar to those obtained with a thermal ionisation mass spectrometer (TIMS). The studies of the isotope composition of other elements yielded relative accuracy and precision at the per mill level too, with characteristic instrument parameters for each element. The relative accuracy and precision of the measurements is degrading with lower element/isotope concentration in a sample. For the elements with abundances below 100 ppm these values drop to the percent level. Depending on the isotopic abundances of Pb in minerals, 207Pb/206Pb ages with accuracy in the range of tens of millions of years can be achieved.

The Science of Exoplanets and Their Systems

A scientific forum on “The Future Science of Exoplanets and Their Systems,” sponsored by Europlanet* and the International Space Science Institute (ISSI)† and co-organized by the Center for Space and Habitability (CSH)‡ of the University of Bern, was held during December 5 and 6, 2012, in Bern, Switzerland. It gathered 24 well-known specialists in exoplanetary, Solar System, and stellar science to discuss the future of the fast-expanding field of exoplanetary research, which now has nearly 1000 objects to analyze and compare and will develop even more quickly over the coming years. The forum discussions included a review of current observational knowledge, efforts for exoplanetary atmosphere characterization and their formation, water formation, atmospheric evolution, habitability aspects, and our understanding of how exoplanets interact with their stellar and galactic environment throughout their history. Several important and timely research areas of focus for further research efforts in the field were identified by the forum participants. These scientific topics are related to the origin and formation of water and its delivery to planetary bodies and the role of the disk in relation to planet formation, including constraints from observations as well as star-planet interaction processes and their consequences for atmosphere-magnetosphere environments, evolution, and habitability. The relevance of these research areas is outlined in this report, and possible themes for future ISSI workshops are identified that may be proposed by the international research community over the coming 2–3 years.