Embedding of human vertebral bodies leads to higher ultimate load and altered damage localisation under axial compression

Computer tomography (CT)-based finite element (FE) models of vertebral bodies assess fracture load in vitro better than dual energy X-ray absorptiometry, but boundary conditions affect stress distribution under the endplates that may influence ultimate load and damage localisation under post-yield strains. Therefore, HRpQCT-based homogenised FE models of 12 vertebral bodies were subjected to axial compression with two distinct boundary conditions: embedding in polymethylmethalcrylate (PMMA) and bonding to a healthy intervertebral disc (IVD) with distinct hyperelastic properties for nucleus and annulus. Bone volume fraction and fabric assessed from HRpQCT data were used to determine the elastic, plastic and damage behaviour of bone. Ultimate forces obtained with PMMA were 22% higher than with IVD but correlated highly (R² = 0.99). At ultimate force, distinct fractions of damage were computed in the endplates (PMMA: 6%, IVD: 70%), cortex and trabecular sub-regions, which confirms previous observations that in contrast to PMMA embedding, failure initiated underneath the nuclei in healthy IVDs. In conclusion, axial loading of vertebral bodies via PMMA embedding versus healthy IVD overestimates ultimate load and leads to distinct damage localisation and failure pattern.     

Calibration of hyperelastic material properties of the human lumbar intervertebral disc under fast dynamic compressive loads

Under fast dynamic loading conditions (e.g. high-energy impact), the load rate dependency of the intervertebral disc (IVD) material properties may play a crucial role in the biomechanics of spinal trauma. However, most finite element models (FEM) of dynamic spinal trauma uses material properties derived from quasi-static experiments, thus neglecting this load rate dependency. The aim of this study was to identify hyperelastic material properties that ensure a more biofidelic simulation of the IVD under a fast dynamic compressive load. A hyperelastic material law based on a first-order Mooney-Rivlin formulation was implemented in a detailed FEM of a L2-L3 functional spinal unit (FSU) to represent the mechanical behavior of the IVD. Bony structures were modeled using an elasto-plastic Johnson-Cook material law that simulates bone fracture while ligaments were governed by a viscoelastic material law. To mimic experimental studies performed in fast dynamic compression, a compressive loading velocity of 1 m/s was applied to the superior half of L2, while the inferior half of L3 was fixed. An exploratory technique was used to simulate dynamic compression of the FSU using 34 sets of hyperelastic material constants randomly selected using an optimal Latin hypercube algorithm and a set of material constants derived from quasi-static experiments. Selection or rejection of the sets of material constants was based on compressive stiffness and failure parameters criteria measured experimentally. The two simulations performed with calibrated hyperelastic constants resulted in nonlinear load-displacement curves with compressive stiffness (7335 and 7079 N/mm), load (12,488 and 12,473 N), displacement (1.95 and 2.09 mm) and energy at failure (13.5 and 14.7 J) in agreement with experimental results (6551 ± 2017 N/mm, 12,411 ± 829 N, 2.1 ± 0.2 mm and 13.0 ± 1.5 J respectively). The fracture pattern and location also agreed with experimental results. The simulation performed with constants derived from quasi-static experiments showed a failure energy (13.2 J) and a fracture pattern and location in agreement with experimental results, but a compressive stiffness (1580 N/mm), a failure load (5976 N) and a displacement to failure (4.8 mm) outside the experimental corridors. The proposed method offers an innovative way to calibrate the hyperelastic material properties of the IVD and to offer a more realistic simulation of the FSU in fast dynamic compression.     

Finite element analyses of human vertebral bodies embedded in polymethylmethalcrylate or loaded via the hyperelastic intervertebral disc models provide equivalent predictions of experimental strength

Quantitative computer tomography (QCT)-based finite element (FE) models of vertebral body provide better prediction of vertebral strength than dual energy X-ray absorptiometry. However, most models were validated against compression of vertebral bodies with endplates embedded in polymethylmethalcrylate (PMMA). Yet, loading being as important as bone density, the absence of intervertebral disc (IVD) affects the strength. Accordingly, the aim was to assess the strength predictions of the classic FE models (vertebral body embedded) against the in vitro and in silico strengths of vertebral bodies loaded via IVDs. High resolution peripheral QCT (HR-pQCT) were performed on 13 segments (T11/T12/L1). T11 and L1 were augmented with PMMA and the samples were tested under a 4° wedge compression until failure of T12. Specimen-specific model was generated for each T12 from the HR-pQCT data. Two FE sets were created: FE-PMMA refers to the classical vertebral body embedded model under axial compression; FE-IVD to their loading via hyperelastic IVD model under the wedge compression as conducted experimentally. Results showed that FE-PMMA models overestimated the experimental strength and their strength prediction was satisfactory considering the different experimental set-up. On the other hand, the FE-IVD models did not prove significantly better (Exp/FE-PMMA: R²=0.68; Exp/FE-IVD: R²=0.71, p=0.84). In conclusion, FE-PMMA correlates well with in vitro strength of human vertebral bodies loaded via real IVDs and FE-IVD with hyperelastic IVDs do not significantly improve this correlation. Therefore, it seems not worth adding the IVDs to vertebral body models until fully validated patient-specific IVD models become available.     

Noninvasive prediction of vertebral body compressive strength using nonlinear finite element method and an image based technique

Noninvasive prediction of vertebral body strength under compressive loading condition is a valuable tool for the assessment of clinical fractures. This paper presents an effective specimen-specific approach for noninvasive prediction of human vertebral strength using a nonlinear finite element (FE) model and an image based parameter based on the quantitative computed tomography (QCT). Nine thoracolumbar vertebrae excised from three cadavers with an average age of 42 years old were used as the samples. The samples were scanned using the QCT. Then, a segmentation technique was performed on each QCT sectional image. The segmented images were then converted into three-dimensional FE models for linear and nonlinear analyses. A new material model was implemented in our nonlinear model being more compatible with real mechanical behavior of trabecular bone. A new image based MOS (Mechanic of Solids) parameter named minimum sectional strength ((σ_uA)_min) was used for the ultimate compressive strength prediction. Subsequently, the samples were destructively tested under uniaxial compression and their experimental ultimate compressive strengths were obtained. Results indicated that our new implemented FE model can predict ultimate compressive strength of human vertebra with a correlation coefficient (R² = 0.94) better than usual linear and nonlinear FE models (R² = 0.83 and 0.85 respectively). The image based parameter introduced in this study ((σ_uA)_min) was also correlated well with the experimental results (R² = 0.86). Although nonlinear FE method with new implemented material model predicts compressive strength better than the (σ_uA)_min, this parameter is clinically more feasible due to its simplicity and lower computational costs. This can make future applications of the (σ_uA)_min more justified for human vertebral body compressive strength prediction.     

A combined numerical and experimental technique for estimation of the forces and moments in the lumbar intervertebral disc

Evaluation of the loads on lumbar intervertebral discs (IVD) is critically important since it is closely related to spine biomechanics, pathology and prosthesis design. Non-invasive estimation of the loads in the discs remains a challenge. In this study, we proposed a new technique to estimate in vivo loads in the IVD using a subject-specific finite element (FE) model of the disc and the kinematics of the disc endplates as input boundary conditions. The technique was validated by comparing the forces and moments in the discs calculated from the FE analyses to the in vitro experiment measurements of three corresponding lumbar discs. The results showed that the forces and moments could be estimated within an average error of 20%. Therefore, this technique can be a promising tool for non-invasive estimation of the loads in the discs and may be extended to be used on living subjects.     

Harmonizing finite element modelling for non-invasive strength estimation by high-resolution peripheral quantitative computed tomography

The finite element (FE) method based on high-resolution peripheral quantitative computed tomography (HR-pQCT) use a variety of tissue constitutive properties and boundary conditions at different laboratories making comparison of mechanical properties difficult. Furthermore, the advent of a second-generation HR-pQCT poses challenges due to improved resolution and a larger region of interest (ROI). This study addresses the need to harmonize results across FE models. The aims are to establish the relationship between FE results as a function of boundary conditions and a range of tissue properties for the first-generation HR-pQCT system, and to determine appropriate model parameters for the second-generation HR-pQCT system. We implemented common boundary conditions and tissue properties on a large cohort (N = 1371), and showed the relationships were highly linear (R² > 0.99) for yield strength and reaction force between FE models. Cadaver radii measured on both generation HR-pQCT with matched ROIs were used to back-calculate a tissue modulus that accounts for the increased resolution (61 µm versus 82 µm), resulting in a modulus of 8748 MPa for second-generation HR-pQCT to produce bone yield strength and reaction force equivalent to using 6829 MPa for first-generation HR-pQCT. Finally, in vivo scans (N = 61) conducted on both generations demonstrated that the larger ROI in the second-generation system results in stronger bone outcome measures, suggesting it is not advisable to convert FE results across HR-pQCT generations without matching ROIs. Together, these findings harmonize FE results by providing a means to compare findings with different boundary conditions and tissue properties, and across scanner generations.     

In vitro extracellular matrix accumulation of nasal and articular chondrocytes for intervertebral disc repair

Chondrocyte based regenerative therapies for intervertebral disc repair such as Autologous Disc Cell Transplantation (ADCT, CODON) and allogeneic juvenile chondrocyte implantation (NuQu^®, ISTO Technologies) have demonstrated good outcomes in clinical trials. However concerns remain with the supply demand reconciliation and issues surrounding immunoreactivity which exist for allogeneic-type technologies. The use of stem cells is challenging due to high growth factor requirements, regulatory barriers and differentiation towards a stable phenotype. Therefore, there is a need to identify alternative non-disc cell sources for the development and clinical translation of next generation therapies for IVD regeneration. In this study, we compared Nasal Chondrocytes (NC) as a non-disc alternative chondrocyte source with Articular Chondrocytes (AC) in terms of cell yield, morphology, proliferation kinetics and ability to produce key extracellular matrix components under 5% and 20% oxygen conditions, with and without exogenous TGF-β supplementation.Results indicated that NC maintained proliferative capacity with high amounts of sGAG and lower collagen accumulation in the absence of TGF-β supplementation under 5% oxygen conditions. Importantly, osteogenesis and calcification was inhibited for NC when cultured in IVD-like microenvironmental conditions. The present study provides a rationale for the exploration of nasal chondrocytes as a promising, potent and clinically feasible autologous cell source for putative IVD repair strategies.     

Loss of tenomodulin expression is a risk factor for age‐related intervertebral disc degeneration

The intervertebral disc (IVD) degeneration is thought to be closely related to ingrowth of new blood vessels. However, the impact of anti‐angiogenic factors in the maintenance of IVD avascularity remains unknown. Tenomodulin (Tnmd) is a tendon/ligament‐specific marker and anti‐angiogenic factor with abundant expression in the IVD. It is still unclear whether Tnmd contributes to the maintenance of IVD homeostasis, acting to inhibit vascular ingrowth into this normally avascular tissue. Herein, we investigated whether IVD degeneration could be induced spontaneously by the absence of Tnmd. Our results showed that Tnmd was expressed in an age‐dependent manner primarily in the outer annulus fibrous (OAF) and it was downregulated at 6 months of age corresponding to the early IVD degeneration stage in mice. Tnmd knockout (Tnmd^?/?) mice exhibited more rapid progression of age‐related IVD degeneration. These signs include smaller collagen fibril diameter, markedly lower compressive stiffness, reduced multiple IVD‐ and tendon/ligament‐related gene expression, induced angiogenesis, and macrophage infiltration in OAF, as well as more hypertrophic‐like chondrocytes in the nucleus pulposus. In addition, Tnmd and chondromodulin I (Chm1, the only homologous gene to Tnmd) double knockout (Tnmd^?/?Chm1^?/?) mice displayed not only accelerated IVD degeneration, but also ectopic bone formation of IVD. Lastly, the absence of Tnmd in OAF‐derived cells promoted p65 and matrix metalloproteinases upregulation, and increased migratory capacity of human umbilical vein endothelial cells. In sum, our data provide clear evidences that Tnmd acts as an angiogenic inhibitor in the IVD homeostasis and protects against age‐related IVD degeneration. Targeting Tnmd may represent a novel therapeutic strategy for attenuating age‐related IVD degeneration.     

A nonlocal constitutive model for trabecular bone softening in compression

Using the three-dimensional morphological data provided by computed tomography, finite element (FE) models can be generated and used to compute the stiffness and strength of whole bones. Three-dimensional constitutive laws capturing the main features of bone mechanical behavior can be developed and implemented into FE software to enable simulations on complex bone structures. For this purpose, a constitutive law is proposed, which captures the compressive behavior of trabecular bone as a porous material with accumulation of irreversible strain and loss of stiffness beyond its yield point and softening beyond its ultimate point. To account for these features, a constitutive law based on damage coupled with hardening anisotropic elastoplasticity is formulated using density and fabric-based tensors. To prevent mesh dependence of the solution, a nonlocal averaging technique is adopted. The law has been implemented into a FE software and some simple simulations are first presented to illustrate its behavior. Finally, examples dealing with compression of vertebral bodies clearly show the impact of softening on the localization of the inelastic process.     

Design and numerical implementation of a 3-D non-linear viscoelastic constitutive model for brain tissue during impact

Finite Element (FE) head models are often used to understand mechanical response of the head and its contents during impact loading in the head. Current FE models do not account for non-linear viscoelastic material behavior of brain tissue. We developed a new non-linear viscoelastic material model for brain tissue and implemented it in an explicit FE code. To obtain sufficient numerical accuracy for modeling the nearly incompressible brain tissue, deviatoric and volumetric stress contributions are separated. Deviatoric stress is modeled in a non-linear viscoelastic differential form. Volumetric behavior is assumed linearly elastic. Linear viscoelastic material parameters were derived from published data on oscillatory experiments, and from ultrasonic experiments. Additionally, non-linear parameters were derived from stress relaxation (SR) experiments at shear strains up to 20%. The model was tested by simulating the transient phase in the SR experiments not used in parameter determination (strains up to 20%, strain rates up to 8s(-1)). Both time- and strain-dependent behavior were predicted accurately (R2>0.96) for strain and strain rates applied. However, the stress was overestimated systematically by approximately 31% independent of strain(rate) applied. This is probably caused by limitations of the experimental data at hand.     

Specimen-specific predictions of contact stress under physiological loading in the human hip: validation and sensitivity studies

Hip osteoarthritis may be initiated and advanced by abnormal cartilage contact mechanics, and finite element (FE) modeling provides an approach with the potential to allow the study of this process. Previous FE models of the human hip have been limited by single specimen validation and the use of quasi-linear or linear elastic constitutive models of articular cartilage. The effects of the latter assumptions on model predictions are unknown, partially because data for the instantaneous behavior of healthy human hip cartilage are unavailable. The aims of this study were to develop and validate a series of specimen-specific FE models, to characterize the regional instantaneous response of healthy human hip cartilage in compression, and to assess the effects of material nonlinearity, inhomogeneity and specimen-specific material coefficients on FE predictions of cartilage contact stress and contact area. Five cadaveric specimens underwent experimental loading, cartilage material characterization and specimen-specific FE modeling. Cartilage in the FE models was represented by average neo-Hookean, average Veronda Westmann and specimen- and region-specific Veronda Westmann hyperelastic constitutive models. Experimental measurements and FE predictions compared well for all three cartilage representations, which was reflected in average RMS errors in contact stress of less than 25 %. The instantaneous material behavior of healthy human hip cartilage varied spatially, with stiffer acetabular cartilage than femoral cartilage and stiffer cartilage in lateral regions than in medial regions. The Veronda Westmann constitutive model with average material coefficients accurately predicted peak contact stress, average contact stress, contact area and contact patterns. The use of subject- and region-specific material coefficients did not increase the accuracy of FE model predictions. The neo-Hookean constitutive model underpredicted peak contact stress in areas of high stress. The results of this study support the use of average cartilage material coefficients in predictions of cartilage contact stress and contact area in the normal hip. The regional characterization of cartilage material behavior provides the necessary inputs for future computational studies, to investigate other mechanical parameters that may be correlated with OA and cartilage damage in the human hip. In the future, the results of this study can be applied to subject-specific models to better understand how abnormal hip contact stress and contact area contribute to OA.     

Stereoselective uptake and degradation of (±)‐o,p‐DDD pesticide stereomers in water‐sediment system

The environmental stereoselective uptake and degradation of (±)‐o,p‐DDD pesticide stereomers in water‐sediment system are described. The results were analyzed by artificial neural network model. The optimized experimental parameters were concentration of o,p‐DDD streamers (7.0 μg L^?1), experimental time (60 min), pH (6), dose (5.0 g L^?1), and temperature (25°C). The maximum uptake and degradation were 87% and 85% and 33.0% and 30.5% for (?)‐ and (+)‐stereomers of o,p‐DDD in 15‐day time. Both uptake and degraded phenomenon showed first‐order rate reaction. Thermodynamic variables indicated exothermic nature of uptake and degradation processes. The uptake and degradation were slightly higher for (?)‐stereomer than (+)‐stereomer of o,p‐DDD. It was assumed that both uptake and degradation processes are accountable for the removal of the streomers of o,p‐DDD from earth's ecosystem, but the uptake is responsible for major contribution. The magnitudes of relative errors obtained by artificial neural network model were in the range of ±0.2 to 3.5, indicating good applicability of the experimental data. The results are very useful to control the environmental contamination due to the chiral o,p‐DDD pesticide as its two enantiomers have different ecological toxicities.     

Animal models of regenerative medicine for biological treatment approaches of degenerative disc diseases

Degenerative disc disease (DDD) is a painful, chronic and progressive disease, which is characterized by inflammation, structural and biological deterioration of the intervertebral disc (IVD) tissues. DDD is specified as cell-, age-, and genetic-dependent degenerative process that can be accelerated by environmental factors. It is one of the major causes of chronic back pain and disability affecting millions of people globally. Current treatment options, such as physical rehabilitation, pain management, and surgical intervention, can provide only temporary pain relief. Different animal models have been used to study the process of IVD degeneration and develop therapeutic options that may restore the structure and function of degenerative discs. Several research works have depicted considerable progress in understanding the biological basis of disc degeneration and the therapeutic potentials of cell transplantation, gene therapy, applications of supporting biomaterials and bioactive factors, or a combination thereof. Since animal models play increasingly significant roles in treatment approaches of DDD, we conducted an electronic database search on Medline through June 2020 to identify, compare, and discuss publications regarding biological therapeutic approaches of DDD that based on intradiscal treatment strategies. We provide an up-to-date overview of biological treatment strategies in animal models including mouse, rat, rabbit, porcine, bovine, ovine, caprine, canine, and primate models. Although no animal model could profoundly reproduce the clinical conditions in humans; animal models have played important roles in specifying our knowledge about the pathophysiology of DDD. They are crucial for developing new therapy approaches for clinical applications.     

The contribution of the glenoid labrum to glenohumeral stability under physiological joint loading using finite element analysis

The labrum contributes to passive glenohumeral joint stability. Cadaveric studies have demonstrated that this has position and load dependency, which has not been quantified under physiological loads. This study aims to validate subject-specific finite element (FE) models against in vitro measurements of joint stability and to utilise the FE models to predict joint stability under physiological loads. The predicted stability values were within ± one standard deviation of experimental data and the FE models showed a reduction in stability of 10–15% with high, physiological, loads. The developed regression equations provide the first representation of passive glenohumeral stability and will aid surgical decision-making.     

Lumbar spine finite element model for healthy subjects: development and validation

Finite element (FE) method is a proven powerful and efficient tool to study the biomechanics of the human lumbar spine. However, due to the large inter-subject variability of geometries and material properties in human lumbar spines, concerns existed on the accuracy and predictive power of one single deterministic FE model with one set of spinal geometry and material properties. It was confirmed that the combined predictions (median or mean value) of several distinct FE models can be used as an improved prediction of behavior of human lumbar spine under identical loading and boundary conditions. In light of this improved prediction, five FE models (L1-L5 spinal levels) of the human lumbar spine were developed based on five healthy living subjects with identical modeling method. The five models were extensively validated through experimental and computational results in the literature. Mesh convergence and material sensitivity analysis were also conducted. We have shown that the results from the five FE models developed in this paper were consistent with the experimental data and simulation results from the existing literature. The validated modeling method introduced in this study can be used in modeling dysfunctional lumber spines such as disc degeneration and scoliosis in future work.     

Finite element formulation of biphasic poroviscoelastic model for articular cartilage.

The purpose of the present study was to develop a computationally efficient finite element model that could be useful for parametric analysis of the biphasic poroviscoelastic (BPVE) behavior of articular cartilage under various loading conditions. The articular cartilage was modeled as the BPVE mixture of a porous, linear viscoelastic, and incompressible solid and an inviscid and incompressible fluid. A finite element (FE) formulation of the BPVE model was developed using two different algorithms, the continuous and discrete spectrum relaxation functions for the viscoelasticity of the solid matrix. These algorithms were applied to the creep and stress relaxation responses to the confined compression of articular cartilage, and a comparison of their performances was made. It was found that the discrete spectrum algorithm significantly saved CPU time and memory, as compared to the continuous spectrum algorithm. The consistency analysis for the present FE formulation was performed in comparison with the IMSL, a commercially available numerical software package. It was found that the present FE formulation yielded consistent results in predicting model behavior, whereas the IMSL subroutine produced inconsistent results in the velocity field, and thereby in the strain calculation.     

Stiffness of photocrosslinkable gelatin hydrogel influences nucleus pulposus cell propertiesin vitro

A key early sign of degenerative disc disease (DDD) is the loss of nucleus pulposus (NP) cells (NPCs). Accordingly, NPC transplantation is a treatment strategy for intervertebral disc (IVD) degeneration. However, in advanced DDD, due to structural damage of the IVD and scaffold mechanical properties, the transplanted cells are less viable and secrete less extracellular matrix, and thus, are unable to efficiently promote NP regeneration. In this study, we evaluated the encapsulation of NPCs in a photosensitive hydrogel made of collagen hydrolysate gelatin and methacrylate (GelMA) to improve NP regeneration. By adjusting the concentration of GelMA, we prepared hydrogels with different mechanical properties. After examining the mechanical properties, cell compatibility and tissue engineering indices of the GelMA-based hydrogels, we determined the optimal hydrogel concentration of the NPC-encapsulating GelMA hydrogel for NP regeneration as 5%. NPCs effectively combined with GelMA and proliferated. As the concentration of the GelMA hydrogel increased, the survival, proliferation and matrix deposition of the encapsulated NPCs gradually decreased, which is the opposite of NPCs grown on the surface of the hydrogel. The controllability of the GelMA hydrogels suggests that these NPC-encapsulating hydrogels are promising candidates to aid in NP tissue engineering and repairing endogenous NPCs.     

A computational study of invariant I5 in a nearly incompressible transversely isotropic model for white matter

The aligned axonal fiber bundles in white matter make it suitable to be modeled as a transversely isotropic material. Recent experimental studies have shown that a minimal form, nearly incompressible transversely isotropic (MITI) material model, is capable of describing mechanical anisotropy of white matter. Here, we used a finite element (FE) computational approach to demonstrate the significance of the fifth invariant (I₅) when modeling the anisotropic behavior of white matter in the large-strain regime. We first implemented and validated the MITI model in an FE simulation framework for large deformations. Next, we applied the model to a plate-hole structural problem to highlight the significance of the invariant I₅ by comparing with the standard fiber reinforcement (SFR) model. We also compared the two models by fitting the experiment data of asymmetric indentation, shear test, and uniaxial stretch of white matter. Our results demonstrated the significance of I₅ in describing shear deformation/anisotropy, and illustrated the potential of the MITI model to characterize transversely isotropic white matter tissues in the large-strain regime.     

Elucidating Turnover Pathways of Bioactive Small Molecules by Isotopomer Analysis: The Persistent Organic Pollutant DDT

The persistent organic pollutant DDT (1,1,1-trichloro-2,2-bis(4-chlorophenyl)ethane) is still indispensable in the fight against malaria, although DDT and related compounds pose toxicological hazards. Technical DDT contains the dichloro congener DDD (1-chloro-4-[2,2-dichloro-1-(4-chlorophenyl)ethyl]benzene) as by-product, but DDD is also formed by reductive degradation of DDT in the environment. To differentiate between DDD formation pathways, we applied deuterium NMR spectroscopy to measure intramolecular deuterium distributions (²H isotopomer abundances) of DDT and DDD. DDD formed in the technical DDT synthesis was strongly deuterium-enriched at one intramolecular position, which we traced back to ²H/¹H fractionation of a chlorination step in the technical synthesis. In contrast, DDD formed by reductive degradation was strongly depleted at the same position, which was due to the incorporation of ²H-depleted hydride equivalents during reductive degradation. Thus, intramolecular isotope distributions give mechanistic information on reaction pathways, and explain a puzzling difference in the whole-molecule ²H/¹H ratio between DDT and DDD. In general, our results highlight that intramolecular isotope distributions are essential to interpret whole-molecule isotope ratios. Intramolecular isotope information allows distinguishing pathways of DDD formation, which is important to identify polluters or to assess DDT turnover in the environment. Because intramolecular isotope data directly reflect isotope fractionation of individual chemical reactions, they are broadly applicable to elucidate transformation pathways of small bioactive molecules in chemistry, physiology and environmental science.     

Nucleus pulposus repair with cultured autologous elastic cartilage derived chondrocytes

Low back pain is one of the most common medical conditions in the Western world. Disc degeneration, an inevitable process of ageing, is one of the major causes of low back pain. Autologous chondrocyte transplantation (ACT) is an increasingly popular method of addressing pathological disorders of cartilage. The purpose of our study was to determine whether autologous chondrocytes from elastic cartilage could survive and synthesise a cartilage specific matrix in the intervertebral disc of rabbits. Sixteen lumbar intervertebral discs (IVD) of New Zealand White rabbits were analysed. In 6 IVD, the nucleus pulposus was evacuated and replaced with tissue engineered autologous chondrocytes from auricular cartilage. In the second group, only the nucleus pulposus was evacuated from 6 IVD, with no chondrocytes implantation. Four non-operated IVD were used as a control. Six months after the operation, the animals were euthanized and the IVD were analysed histologically. Autologous cartilage implants were well tolerated by the host for up to six months in vivo. There was only hyaline-like cartilage in the place of the nucleus pulposus. We could not detect any elastic fibres in the new cartilage matrix. In IVD from which only the nucleus pulposus was evacuated and no chondrocytes were implanted, just fibrous tissue was found instead of nucleus pulposus. The overall histological analysis of new cartilage produced after implantation in our study confirmed the hypothesis that ACT from auricular cartilage can be implanted into the IVD instead of the nucleus pulposus and that a significant percentage of implanted chondrocytes survive and produce hyaline-like cartilage.