Finite element analysis predicts experimental failure patterns in vertebral bodies loaded via intervertebral discs up to large deformation

Vertebral compression fracture is a common medical problem in osteoporotic individuals. The quantitative computed tomography (QCT)-based finite element (FE) method may be used to predict vertebral strength in vivo, but needs to be validated with experimental tests. The aim of this study was to validate a nonlinear anatomy specific QCT-based FE model by using a novel testing setup. Thirty-seven human thoracolumbar vertebral bone slices were prepared by removing cortical endplates and posterior elements. The slices were scanned with QCT and the volumetric bone mineral density (vBMD) was computed with the standard clinical approach. A novel experimental setup was designed to induce a realistic failure in the vertebral slices in vitro. Rotation of the loading plate was allowed by means of a ball joint. To minimize device compliance, the specimen deformation was measured directly on the loading plate with three sensors. A nonlinear FE model was generated from the calibrated QCT images and computed vertebral stiffness and strength were compared to those measured during the experiments. In agreement with clinical observations, most of the vertebrae underwent an anterior wedge-shape fracture. As expected, the FE method predicted both stiffness and strength better than vBMD (R2 improved from 0.27 to 0.49 and from 0.34 to 0.79, respectively). Despite the lack of fitting parameters, the linear regression of the FE prediction for strength was close to the 1:1 relation (slope and intercept close to one (0.86 kN) and to zero (0.72 kN), respectively). In conclusion, a nonlinear FE model was successfully validated through a novel experimental technique for generating wedge-shape fractures in human thoracolumbar vertebrae.

Morphological ripening of fluid inclusions and coupled zone-refining in quartz crystals revealed by cathodoluminescence imaging: Implications for CL-petrography, fluid inclusion analysis and trace-element geothermometry

Cathodoluminescence (CL) studies have previously shown that some secondary fluid inclusions in luminescent quartz are surrounded by dark, non-luminescent patches, resulting from fracture-sealing by late, trace-element-poor quartz. This finding has led to the tacit generalization that all dark CL patches indicate influx of low temperature, late-stage fluids. In this study we have examined natural and synthetic hydrothermal quartz crystals using CL imaging supplemented by in-situ elemental analysis. The results lead us to propose that all natural, liquid-water-bearing inclusions in quartz, whether trapped on former crystal growth surfaces (i.e., of primary origin) or in healed fractures (i.e., of pseudosecondary or secondary origin), are surrounded by three-dimensional, non-luminescent patches. Cross-cutting relations show that the patches form after entrapment of the fluid inclusions and therefore they are not diagnostic of the timing of fluid entrapment. Instead, the dark patches reveal the mechanism by which fluid inclusions spontaneously approach morphological equilibrium and purify their host quartz over geological time. Fluid inclusions that contain solvent water perpetually dissolve and reprecipitate their walls, gradually adopting low-energy euhedral and equant shapes. Defects in the host quartz constitute solubility gradients that drive physical migration of the inclusions over distances of tens of μm (commonly) up to several mm (rarely). Inclusions thus sequester from their walls any trace elements (e.g., Li, Al, Na, Ti) present in excess of equilibrium concentrations, thereby chemically purifying their host crystals in a process analogous to industrial zone refining. Non-luminescent patches of quartz are left in their wake. Fluid inclusions that contain no liquid water as solvent (e.g., inclusions of low-density H2O vapor or other non-aqueous volatiles) do not undergo this process and therefore do not migrate, do not modify their shapes with time, and are not associated with dark-CL zone-refined patches. This new understanding has implications for the interpretation of solids within fluid inclusions (e.g., Ti- and Al-minerals) and for the elemental analysis of hydrothermal and metamorphic quartz and its fluid inclusions by microbeam methods such as LA-ICPMS and SIMS. As Ti is a common trace element in quartz, its sequestration by fluid inclusions and its depletion in zone-refined patches impacts on applications of the Ti-in-quartz geothermometer.

Fluid-related inclusions in Alpine high-pressure peridotite reveal trace element recycling during subduction-zone dehydration of serpentinized mantle (Cima di Gagnone, Swiss Alps).

Serpentinites release at sub-arc depths volatiles and several fluid-mobile trace elements found in arc magmas. Constraining element uptake in these rocks and defining the trace element composition of fluids released upon serpentinite dehydration can improve our understanding of mass transfer across subduction zones and to volcanic arcs. The eclogite-facies garnet metaperidotite and chlorite harzburgite bodies embedded in paragneiss of the subduction melange from Cima di Gagnone derive from serpentinized peridotite protoliths and are unique examples of ultramafic rocks that experienced subduction metasomatism and devolatilization. In these rocks, metamorphic olivine and garnet trap polyphase inclusions representing the fluid released during high-pressure breakdown of antigorite and chlorite. Combining major element mapping and laser-ablation ICP-MS bulk inclusion analysis, we characterize the mineral content of polyphase inclusions and quantify the fluid composition. Silicates, Cl-bearing phases, sulphides, carbonates, and oxides document post-entrapment mineral growth in the inclusions starting immediately after fluid entrapment. Compositional data reveal the presence of two different fluid types. The first (type A) records a fluid prominently enriched in fluid-mobile elements, with Cl, Cs, Pb, As, Sb concentrations up to 10(3) PM (primitive mantle), similar to 10(2) PM Tit Ba, while Rb, B, Sr, Li, U concentrations are of the order of 10(1) PM, and alkalis are similar to 2 PM. The second fluid (type B) has considerably lower fluid-mobile element enrichments, but its enrichment patterns are comparable to type A fluid. Our data reveal multistage fluid uptake in these peridotite bodies, including selective element enrichment during seafloor alteration, followed by fluid-rock interaction along with subduction metamorphism in the plate interface melange. Here, infiltration of sediment-equilibrated fluid produced significant enrichment of the serpentinites in As, Sb, B, Pb, an enriched trace element pattern that was then transferred to the fluid released at greater depth upon serpentine dehydration (type A fluid). The type B fluid hosted by garnet may record the composition of the chlorite breakdown fluid released at even greater depth. The Gagnone study-case demonstrates that serpentinized peridotites acquire water and fluid-mobile elements during ocean floor hydration and through exchange with sediment-equilibrated fluids in the early subduction stages. Subsequent antigorite devolatilization at subarc depths delivers aqueous fluids to the mantle wedge that can be prominently enriched in sediment-derived components, potentially triggering arc magmatism without the need of concomitant dehydration/melting of metasediments or altered oceanic crust.

Comparison of proximal femur and vertebral body strength improvements in the FREEDOM trial using an alternative finite element methodology

Denosumab reduced the incidence of new fractures in postmenopausal women with osteoporosis by 68% at the spine and 40% at the hip over 36 months compared with placebo in the FREEDOM study. This efficacy was supported by improvements from baseline in vertebral (18.2%) strength in axial compression and femoral (8.6%) strength in sideways fall configuration at 36 months, estimated in Newtons by an established voxel-based finite element (FE) methodology. Since FE analyses rely on the choice of meshes, material properties, and boundary conditions, the aim of this study was to independently confirm and compare the effects of denosumab on vertebral and femoral strength during the FREEDOM trial using an alternative smooth FE methodology. Unlike the previous FE study, effects on femoral strength in physiological stance configuration were also examined. QCT data for the proximal femur and two lumbar vertebrae were analyzed by smooth FE methodology at baseline, 12, 24, and 36 months for 51 treated (denosumab) and 47 control (placebo) subjects. QCT images were segmented and converted into smooth FE models to compute bone strength. L1 and L2 vertebral bodies were virtually loaded in axial compression and the proximal femora in both fall and stance configurations. Denosumab increased vertebral body strength by 10.8%, 14.0%, and 17.4% from baseline at 12, 24, and 36 months, respectively (p < 0.0001). Denosumab also increased femoral strength in the fall configuration by 4.3%, 5.1%, and 7.2% from baseline at 12, 24, and 36 months, respectively (p < 0.0001). Similar improvements were observed in the stance configuration with increases of 4.2%, 5.2%, and 5.2% from baseline (p ≤ 0.0007). Differences between the increasing strengths with denosumab and the decreasing strengths with placebo were significant starting at 12 months (vertebral and femoral fall) or 24 months (femoral stance). Using an alternative smooth FE methodology, we confirmed the significant improvements in vertebral body and proximal femur strength previously observed with denosumab. Estimated increases in strength with denosumab and decreases with placebo were highly consistent between both FE techniques.

Clinical Use of Quantitative Computed Tomography–Based Finite Element Analysis of the Hip and Spine in the Management of Osteoporosis in Adults: the 2015 ISCD Official Positions—Part II

The International Society for Clinical Densitometry (ISCD) has developed new official positions for the clinical use of quantitative computed tomography (QCT)-based finite element analysis of the spine and hip. The ISCD task force for QCT reviewed the evidence for clinical applications and presented a report with recommendations at the 2015 ISCD Position Development Conference. Here we discuss the agreed upon ISCD official positions with supporting medical evidence, rationale, controversy, and suggestions for further study. Parts I and III address the clinical use of QCT of the hip, and the clinical feasibility of existing techniques for opportunistic screening of osteoporosis using CT scans obtained for other diagnosis such as colonography was addressed.

Nonlinear quasi-static finite element simulations predict in vitro strength of human proximal femora assessed in a dynamic sideways fall setup

Osteoporotic proximal femur fractures are caused by low energy trauma, typically when falling on the hip from standing height. Finite element simulations, widely used to predict the fracture load of femora in fall, usually include neither mass-related inertial effects, nor the viscous part of bone's material behavior. The aim of this study was to elucidate if quasi-static non-linear homogenized finite element analyses can predict in vitro mechanical properties of proximal femora assessed in dynamic drop tower experiments. The case-specific numerical models of thirteen femora predicted the strength (R2=0.84, SEE=540 N, 16.2%), stiffness (R2=0.82, SEE=233 N/mm, 18.0%) and fracture energy (R2=0.72, SEE=3.85 J, 39.6%); and provided fair qualitative matches with the fracture patterns. The influence of material anisotropy was negligible for all predictions. These results suggest that quasi-static homogenized finite element analysis may be used to predict mechanical properties of proximal femora in the dynamic sideways fall situation.

Comparison of Mixed and Kinematic Uniform Boundary Conditions in Homogenized Elasticity of Femoral Trabecular Bone Using Microfinite Element Analyses

Mechanical properties of human trabecular bone play an important role in age-related bone fragility and implant stability. Micro-finite element (microFE) analysis allows computing the apparent elastic properties of trabecular bone biopsies, but the results depend on the type of applied boundary conditions (BCs). In this study, 167 femoral trabecular cubic biopsies with a side length of 5.3 mm were analyzed using microFE analysis to compare their stiffness systematically with kinematic uniform boundary conditions (KUBCs) and periodicity-compatible mixed uniform boundary conditions (PMUBCs). The obtained elastic constants were then used in the volume fraction and fabric-based orthotropic Zysset-Curnier model to identify their respective model parameters. As expected, PMUBCs lead to more compliant apparent elastic properties than KUBCs, especially in shear. The differences in stiffness decreased with bone volume fraction and mean intercept length. Unlike KUBCs, PMUBCs were sensitive to heterogeneity of the biopsies. The Zysset-Curnier model predicted apparent elastic constants successfully in both cases with adjusted coefficients of determination of 0.986 for KUBCs and 0.975 for PMUBCs. The role of these boundary conditions in finite element analyses of whole bones and bone-implant systems will need to be investigated in future work.

Biomechanical validation of computer assisted planning of periacetabular osteotomy: A preliminary study based on finite element analysis

Periacetabular Osteotomy (PAO) is a joint preserving surgical intervention intended to increase femoral head coverage and thereby to improve stability in young patients with hip dysplasia. Previously, we developed a CT-based, computer-assisted program for PAO diagnosis and planning, which allows for quantifying the 3D acetabular morphology with parameters such as acetabular version, inclination, lateral center edge (LCE) angle and femoral head coverage ratio (CO). In order to verify the hypothesis that our morphology-based planning strategy can improve biomechanical characteristics of dysplastic hips, we developed a 3D finite element model based on patient-specific geometry to predict cartilage contact stress change before and after morphology-based planning. Our experimental results demonstrated that the morphology-based planning strategy could reduce cartilage contact pressures and at the same time increase contact areas. In conclusion, our computer-assisted system is an efficient tool for PAO planning.

The supergene thorium and rare-earth element deposit at Morro do Ferro, Poços de Caldas, Minas Gerais, Brazil

The thorium and rare-earth element (Th-REE) deposit at Morro do Ferro formed under supergene lateritic weathering conditions. The ore body consists of shallow NW-SE elongated argillaceous lenses that extend from the top of the hill downwards along its south-eastern slope. The deposit is capped by a network of magnetite layers which protected the underlying highly weathered, argillaceous host rock from excessive erosion. The surrounding country rocks comprise a sequence of subvolcanic phonolite intrusions that have been strongly altered by hydrothermal and supergene processes. From petrological, mineralogical and geochemical studies, and mass balance calculations, it is inferred that the highly weathered host rock was originally carbonatitic in composition, initially enriched in Th and REEs compared to the surrounding silicate rocks. The intrusion of the carbonatite caused fenitic alteration in the surrounding phonolites, consisting of early potassic alteration followed by a vein-type Th-REE mineralization with associated fluorite, carbonate, pyrite and zircon. Subsequent weathering has completely decomposed the carbonatite forming a residual supergene enrichment of Th and REEs. Initial weathering of the carbonatite has created a chemical environment that might have been conductive to carbonate and phosphate complexing of the REEs in groundwaters. This may have appreciably restricted the dissolution of primary REE phases. Strongly oxidic weathering has resulted in a fractionation between Ce and the other light rare earth elements (LREEs). Ce3+ is oxidized to Ce4+ and retained together with Th by secondary mineral formation (cerianite, thorianite), and by adsorption on poorly crystalline iron- and aluminium-hydroxides. In contrast, the trivalent LREEs are retained to a lesser degree and are thus more available for secondary mineral formation (Nd-lanthanite) and adsorption at greater depths down the weathering column. Seasonally controlled fluctuations of recharge waters into the weathering column may help to explain the observed repetition of Th-Ce enriched zones underlain by trivalent LREE enriched zones.

A transposable element insertion in APOB causes cholesterol deficiency in Holstein cattle

Cholesterol deficiency, a new autosomal recessive inherited genetic defect in Holstein cattle, has been recently reported to have an influence on the rearing success of calves. The affected animals show unresponsive diarrhea accompanied by hypocholesterolemia and usually die within the first weeks or months of life. Here, we show that whole genome sequencing combined with the knowledge about the pedigree and inbreeding status of a livestock population facilitates the identification of the causative mutation. We resequenced the entire genomes of an affected calf and a healthy partially inbred male carrying one copy of the critical 2.24-Mb chromosome 11 segment in its ancestral state and one copy of the same segment with the cholesterol deficiency mutation. We detected a single structural variant, homozygous in the affected case and heterozygous in the non-affected carrier male. The genetic makeup of this key animal provides extremely strong support for the causality of this mutation. The mutation represents a 1.3kb insertion of a transposable LTR element (ERV2-1) in the coding sequence of the APOB gene, which leads to truncated transcripts and aberrant splicing. This finding was further supported by RNA sequencing of the liver transcriptome of an affected calf. The encoded apolipoprotein B is an essential apolipoprotein on chylomicrons and low-density lipoproteins, and therefore, the mutation represents a loss of function mutation similar to autosomal recessive inherited familial hypobetalipoproteinemia-1 (FHBL1) in humans. Our findings provide a direct gene test to improve selection against this deleterious mutation in Holstein cattle.

Nitinol Stent Oversizing in Patients Undergoing Popliteal Artery Revascularization: A Finite Element Study

Nitinol stent oversizing is frequently performed in peripheral arteries to ensure a desirable lumen gain. However, the clinical effect of mis-sizing remains controversial. The goal of this study was to provide a better understanding of the structural and hemodynamic effects of Nitinol stent oversizing. Five patient-specific numerical models of non-calcified popliteal arteries were developed to simulate the deployment of Nitinol stents with oversizing ratios ranging from 1.1 to 1.8. In addition to arterial biomechanics, computational fluid dynamics methods were adopted to simulate the physiological blood flow inside the stented arteries. Results showed that stent oversizing led to a limited increase in the acute lumen gain, albeit at the cost of a significant increase in arterial wall stresses. Furthermore, localized areas affected by low Wall Shear Stress increased with higher oversizing ratios. Stents were also negatively impacted by the procedure as their fatigue safety factors gradually decreased with oversizing. These adverse effects to both the artery walls and stents may create circumstances for restenosis. Although the ideal oversizing ratio is stent-specific, this study showed that Nitinol stent oversizing has a very small impact on the immediate lumen gain, which contradicts the clinical motivations of the procedure.

Variable effects of the conserved RNA hairpin element upon 3' end processing of histone pre-mRNA in vitro.

We have studied the requirements for efficient histone-specific RNA 3' processing in nuclear extract from mammalian tissue culture cells. Processing is strongly impaired by mutations in the pre-mRNA spacer element that reduce the base-pairing potential with U7 RNA. Moreover, by exchanging the hairpin and spacer elements of two differently processed H4 genes, we find that this difference is exclusively due to the spacer element. Finally, processing is inhibited by the addition of competitor RNAs, if these contain a wild-type spacer sequence, but not if their spacer element is mutated. Conversely, the importance of the hairpin for histone RNA 3' processing is highly variable: A hairpin mutant of the H4-12 gene is processed with almost wild-type efficiency in extract from K21 mouse mastocytoma cells but is strongly affected in HeLa cell extract, whereas an identical hairpin mutant of the H4-1 gene is affected in both extracts. The hairpin defect of H4-12-specific RNA in HeLa cells can be overcome by a compensatory mutation that increases the base complementarity to U7 snRNA. Very similar results were also obtained in RNA competition experiments: processing of H4-12-specific RNA can be competed by RNA carrying a wild-type hairpin element in extract from HeLa, but not K21 cells, whereas processing of H4-1-specific RNA can be competed in both extracts. With two additional histone genes we obtained results that were in one case intermediate and in the other similar to those obtained with H4-1. These results suggest that hairpin binding factor(s) can cooperatively support the ability of U7 snRNPs to form an active processing complex, but is(are) not directly involved in the processing mechanism.