BMP-2 liberated from biomimetic implant coatings induces and sustains direct ossification in an ectopic rat model

INTRODUCTION: Using a rat model, we evaluated the kinetics and histomorphometry of ectopic bone formation in association with biomimetic implant coatings containing BMP-2. MATERIALS AND METHODS: One experimental and three control groups were set up: titanium-alloy discs coated with a biomimetically co-precipitated layer of calcium phosphate and BMP-2 [1.7 microg per disc (incorporated-BMP group)]; uncoated discs (control); discs biomimetically coated with a layer of calcium phosphate alone (control); and discs biomimetically coated with a layer of calcium phosphate bearing superficially adsorbed BMP-2 [0.98 microg per disc (control)]. Discs (n = 6 per group) were implanted subcutaneously in rats and retrieved at 7-day intervals over a period of 5 weeks for kinetic, histomorphometrical, morphological and histochemical analyses. RESULTS: In the incorporated-BMP-2 group, osteogenic activity was first observed 2 weeks after implantation and thereafter continued unabated until the end of the monitoring period. The net weekly rates of bone formation per disc were 5.8 mm3 at 2 weeks and 3.64 mm3 at 5 weeks. The total volumes of bone formed per disc at these junctures were 5.8 mm3 and 10.3 mm3, respectively. Bone tissue, which was formed by a direct ossification mechanism, was deposited at distances of up to 340 microm from the implant surfaces. The biomimetic coatings were degraded gradually, initially by foreign body giant cells alone and then also by osteoclasts. Forty percent of the coating material (and thus presumably of the incorporated BMP-2) remained at the end of the monitoring period. Hence, 60% of the incorporated BMP-2 had been released. At this 5-week juncture, no bone tissue was associated with any of the control implants. CONCLUSION: BMP-2 incorporated into biomimetic calcium phosphate coatings is capable not only of inducing bone formation at an ectopic site in vivo but also of doing so with a very high potency at a low pharmacological level, and of sustaining this activity for a considerable period of time. The sustainment of osteogenic activity is of great clinical importance for the osseointegration of dental and orthopedic implants.

Land use change and nitrogen feedbacks constrain the trajectory of the land carbon sink

Our understanding of Earth's carbon climate system depends critically upon interactions between rising atmospheric CO2, changing land use, and nitrogen limitation on vegetation growth. Using a global land model, we show how these factors interact locally to generate the global land carbon sink over the past 200 years. Nitrogen constraints were alleviated by N2 fixation in the tropics and by atmospheric nitrogen deposition in extratropical regions. Nonlinear interactions between land use change and land carbon and nitrogen cycling originated from three major mechanisms: (i) a sink foregone that would have occurred without land use conversion; (ii) an accelerated response of secondary vegetation to CO2 and nitrogen, and (iii) a compounded clearance loss from deforestation. Over time, these nonlinear effects have become increasingly important and reduce the present-day net carbon sink by ~40% or 0.4 PgC yr−1.

Long-term injury induced plasticity in {\it Aplysia\/} neurons: Behavioral, electrophysiological, and morphological studies

Nerve injury is known to produce a variety of electrophysiological and morphological neuronal alterations (reviewed by Titmus and Faber, 1990; Bulloch and Ridgeway, 1989; Walters, 1994). Determining if these alterations are adaptive and how they are activated and maintained could provide important insight into basic cellular mechanisms of injury-induced plasticity. Furthermore, characterization of injury-induced plasticity provides a useful assay system for the identification of possible induction signals underlying these neuronal changes. Understanding fundamental mechanisms and underlying induction signals of injury-induced neuronal plasticity could facilitate development of treatment strategies for neural injury and neuropathic pain in humans.^ This dissertation characterizes long-lasting, injury-induced neuronal alterations using the nervous system of Aplysia californica as a model. These changes are examined at the behavioral, electrophysiological, and morphological levels. Injury-induced changes in the electrophysiological properties of neurons were found that increased the signaling effectiveness of the injured neurons. This increase in signalling effectiveness could act to compensate for partial destruction of the injured neuron's peripheral processes. Recovery of a defensive behavioral response which serves to protect the animal from further injury was found within 2 weeks of injury. For the behavioral recovery to occur, new neural pathways must have been formed between the denervated area and the CNS. This was found to be mediated at least in part by new axonal growth which extended from the injured cell back along the original pathway (i.e. into the injured nerve). In addition, injury produced central axonal sprouting into different nerves that do not usually contain the injured neuron's axons. This could be important for (i) finding alternative pathways to the periphery when the original pathways are impassable and (ii) the formation of additional synaptic connections with post-synaptic targets which would further enhance the signalling effectiveness of the injured cell. ^

Quantitative MRI of spinal cord injury in a rat model: Correlative studies

Serial quantitative and correlative studies of experimental spinal cord injury (SCI) in rats were conducted using three-dimensional magnetic resonance imaging (MRI). Correlative measures included morphological histopathology, neurobehavioral measures of functional deficit, and biochemical assays for N-acetyl-aspartate (NAA), lactate, pyruvate, and ATP. A spinal cord injury device was characterized and provided a reproducible injury severity. Injuries were moderate and consistent to within $\pm$20% (standard deviation). For MRI, a three-dimensional implementation of the single spin-echo FATE (Fast optimum angle, short TE) pulse sequence was used for rapid acquisition, with a 128 x 128 x 32 (x,y,z) matrix size and a 0.21 x 0.21 x 1.5 mm resolution. These serial studies revealed a bimodal characteristic in the evolution in MRI pathology with time. Early and late phases of SCI pathology were clearly visualized in $T\sb2$-weighted MRI, and these corresponded to specific histopathological changes in the spinal cord. Centralized hypointense MRI regions correlated with evidence of hemorrhagic and necrotic tissue, while surrounding hyperintense regions represented edema or myelomalacia. Unexpectedly, $T\sb2$-weighted MRI pathology contrast at 24 hours after injury appeared to subside before peaking at 72 hours after injury. This change is likely attributable to ongoing secondary injury processes, which may alter local $T\sb2$ values or reduce the natural anisotropy of the spinal cord. MRI, functional, and histological measures all indicated that 72 hours after injury was the temporal maximum for quantitative measures of spinal cord pathology. Thereafter, significant improvement was seen only in neurobehavioral scores. Significant correlations were found between quantitated MRI pathology and histopathology. Also, NAA and lactate levels correlated with behavioral measures of the level of function deficit. Asymmetric (rostral/caudal) changes in NAA and lactate due to injury indicate that rostral and caudal segments from the injury site are affected differently by the injury. These studies indicate that volumetric quantitation of MRI pathology from $T\sb2$-weighted images may play an important role in early prediction of neurologic deficit and spinal cord pathology. The loss of $T\sb2$ contrast at 24 hours suggests MR may be able to detect certain delayed mechanisms of secondary injury which are not resolved by histopathology or other radiological modalities. Furthermore, in vivo proton magnetic resonance spectroscopy (MRS) studies of SCI may provide a valuable addition source of information about changes in regional spinal cord lactate and NAA levels, which are indicative of local metabolic and pathological changes. ^

Bone Volume Fraction and Fabric Anisotropy Are Better Determinants of Trabecular Bone Stiffness than Other Morphological Variables

As our population ages, more individuals suffer from osteoporosis. This disease leads to impaired trabecular architecture and increased fracture risk. It is essential to understand how morphological and mechanical properties of the cancellous bone are related. Morphologyelasticity relationships based on bone volume fraction (BV/TV) and fabric anisotropy explain up to 98% of the variation in elastic properties. Yet, other morphological variables such as individual trabeculae segmentation (ITS) and trabecular bone score (TBS) could improve the stiffness predictions. A total of 743 micro-computed tomography reconstructions of cubic trabecular bone samples extracted from femur, radius, vertebrae and iliac crest were analysed. Their morphology was assessed via 25 variables and their stiffness tensor (inline image) was computed from six independent load cases using micro finite element analyses. Variance inflation factors were calculated to evaluate collinearity between morphological variables and decide upon their inclusion in morphology-elasticity relationships. The statistically admissible morphological variables were included in a multi-linear regression modelling the dependent variable inline image. The contribution of each independent variable was evaluated (ANOVA). Our results show that BV/TV is the best determinant of inline image (inline image=0.889), especially in combination with fabric (inline image=0.968). Including the other independent predictors hardly affected the amount of variance explained by the model (inline image=0.975). Across all anatomical sites, BV/TV explained 87% of the variance of the bone elastic properties. Fabric further described 10% of the bone stiffness, but the improvement in variance explanation by adding other independent factors was marginal (<1%). These findings confirm that BV/TV and fabric are the best determinants of trabecular bone stiffness and show, against common belief, that other morphological variables do not bring any further contribution. These overall conclusions remain to be confirmed for specific bone diseases and post-elastic properties.

Semiautomatic Cochleostomy Target and Insertion Trajectory Planning for Minimally Invasive Cochlear Implantation

A major component of minimally invasive cochlear implantation is atraumatic scala tympani (ST) placement of the electrode array. This work reports on a semiautomatic planning paradigm that uses anatomical landmarks and cochlear surface models for cochleostomy target and insertion trajectory computation. The method was validated in a human whole head cadaver model (n = 10 ears). Cochleostomy targets were generated from an automated script and used for consecutive planning of a direct cochlear access (DCA) drill trajectory from the mastoid surface to the inner ear. An image-guided robotic system was used to perform both, DCA and cochleostomy drilling. Nine of 10 implanted specimens showed complete ST placement. One case of scala vestibuli insertion occurred due to a registration/drilling error of 0.79 mm. The presented approach indicates that a safe cochleostomy target and insertion trajectory can be planned using conventional clinical imaging modalities, which lack sufficient resolution to identify the basilar membrane.

A novel animal model for external anal sphincter insufficiency

Abstract PURPOSE: Reliable animal models are essential to evaluate future therapeutic options like cell-based therapies for external anal sphincter insufficiency. The goal of our study was to describe the most reliable model for external sphincter muscle insufficiency by comparing three different methods to create sphincter muscle damage. METHODS: In an experimental animal study, female Lewis rats (200-250 g) were randomly assigned to three treatment groups (n = 5, each group). The external sphincter muscle was weakened in the left dorsal quadrant by microsurgical excision, cryosurgery, or electrocoagulation by diathermy. Functional evaluation included in vivo measurements of resting pressure, spontaneous muscle contraction, and contraction in response to electrical stimulation of the afferent nerve at baseline and at 2, 4, and 6 weeks after sphincter injury. Masson's trichrome staining and immunofluorescence for skeletal muscle markers was performed for morphological analysis. RESULTS: Peak contraction after electrical stimulation was significantly decreased after sphincter injury in all groups. Contraction forces recovered partially after cryosurgery and electrocoagulation but not after microsurgical excision. Morphological analysis revealed an incomplete destruction of the external sphincter muscle in the cryosurgery and electrocoagulation groups compared to the microsurgery group. CONCLUSIONS: For the first time, three different models of external sphincter muscle insufficiency were directly compared. The animal model using microsurgical sphincter destruction offers the highest level of consistency regarding tissue damage and sphincter insufficiency, and therefore represents the most reliable model to evaluate future therapeutic options. In addition, this study represents a novel model to specifically test the external sphincter muscle function.

Comparison of morphological and rheological conditions between conventional and eversion carotid endarterectomy using computational fluid dynamics - a pilot study

PURPOSE To compare postoperative morphological and rheological conditions after eversion carotid endarterectomy versus conventional carotid endarterectomy using computational fluid dynamics. BASIC METHODS Hemodynamic metrics (velocity, wall shear stress, time-averaged wall shear stress and temporal gradient wall shear stress) in the carotid arteries were simulated in one patient after conventional carotid endarterectomy and one patient after eversion carotid endarterectomy by computational fluid dynamics analysis based on patient specific data. PRINCIPAL FINDINGS Systolic peak of the eversion carotid endarterectomy model showed a gradually decreased pressure along the stream path, the conventional carotid endarterectomy model revealed high pressure (about 180 Pa) at the carotid bulb. Regions of low wall shear stress in the conventional carotid endarterectomy model were much larger than that in the eversion carotid endarterectomy model and with lower time-averaged wall shear stress values (conventional carotid endarterectomy: 0.03-5.46 Pa vs. eversion carotid endarterectomy: 0.12-5.22 Pa). CONCLUSIONS Computational fluid dynamics after conventional carotid endarterectomy and eversion carotid endarterectomy disclosed differences in hemodynamic patterns. Larger studies are necessary to assess whether these differences are consistent and might explain different rates of restenosis in both techniques.

Effect of surface restoring on subsurface variability in a climate model during 1949–2005

Initializing the ocean for decadal predictability studies is a challenge, as it requires reconstructing the little observed subsurface trajectory of ocean variability. In this study we explore to what extent surface nudging using well-observed sea surface temperature (SST) can reconstruct the deeper ocean variations for the 1949–2005 period. An ensemble made with a nudged version of the IPSLCM5A model and compared to ocean reanalyses and reconstructed datasets. The SST is restored to observations using a physically-based relaxation coefficient, in contrast to earlier studies, which use a much larger value. The assessment is restricted to the regions where the ocean reanalyses agree, i.e. in the upper 500 m of the ocean, although this can be latitude and basin dependent. Significant reconstruction of the subsurface is achieved in specific regions, namely region of subduction in the subtropical Atlantic, below the thermocline in the equatorial Pacific and, in some cases, in the North Atlantic deep convection regions. Beyond the mean correlations, ocean integrals are used to explore the time evolution of the correlation over 20-year windows. Classical fixed depth heat content diagnostics do not exhibit any significant reconstruction between the different existing bservation-based references and can therefore not be used to assess global average time-varying correlations in the nudged simulations. Using the physically based average temperature above an isotherm (14°C) alleviates this issue in the tropics and subtropics and shows significant reconstruction of these quantities in the nudged simulations for several decades. This skill is attributed to the wind stress reconstruction in the tropics, as already demonstrated in a perfect model study using the same model. Thus, we also show here the robustness of this result in an historical and observational context.

Temperature Prediction Model for Bone Drilling Based on Density Distribution and In Vivo Experiments for Minimally Invasive Robotic Cochlear Implantation

Surgical robots have been proposed ex vivo to drill precise holes in the temporal bone for minimally invasive cochlear implantation. The main risk of the procedure is damage of the facial nerve due to mechanical interaction or due to temperature elevation during the drilling process. To evaluate the thermal risk of the drilling process, a simplified model is proposed which aims to enable an assessment of risk posed to the facial nerve for a given set of constant process parameters for different mastoid bone densities. The model uses the bone density distribution along the drilling trajectory in the mastoid bone to calculate a time dependent heat production function at the tip of the drill bit. Using a time dependent moving point source Green's function, the heat equation can be solved at a certain point in space so that the resulting temperatures can be calculated over time. The model was calibrated and initially verified with in vivo temperature data. The data was collected in minimally invasive robotic drilling of 12 holes in four different sheep. The sheep were anesthetized and the temperature elevations were measured with a thermocouple which was inserted in a previously drilled hole next to the planned drilling trajectory. Bone density distributions were extracted from pre-operative CT data by averaging Hounsfield values over the drill bit diameter. Post-operative [Formula: see text]CT data was used to verify the drilling accuracy of the trajectories. The comparison of measured and calculated temperatures shows a very good match for both heating and cooling phases. The average prediction error of the maximum temperature was less than 0.7 °C and the average root mean square error was approximately 0.5 °C. To analyze potential thermal damage, the model was used to calculate temperature profiles and cumulative equivalent minutes at 43 °C at a minimal distance to the facial nerve. For the selected drilling parameters, temperature elevation profiles and cumulative equivalent minutes suggest that thermal elevation of this minimally invasive cochlear implantation surgery may pose a risk to the facial nerve, especially in sclerotic or high density mastoid bones. Optimized drilling parameters need to be evaluated and the model could be used for future risk evaluation.

Pinus nigra (European black pine) as the dominant species of the last glacial pinewoods in south-western to central Iberia: a morphological study of modern and fossil pollen

Aim Our aim was to discriminate different species of Pinus via pollen analysis in order to assess the responses of particular pine species to orbital and millennial-scale climate changes, particularly during the last glacial period. Location Modern pollen grains were collected from current pine populations along transects from the Pyrenees to southern Iberia and the Balearic Islands. Fossil pine pollen was recovered from the south-western Iberian margin core MD95-2042. Methods We measured a set of morphological traits of modern pollen from the Iberian pine species Pinus nigra, P. sylvestris, P. halepensis, P. pinea and P. pinaster and of fossil pine pollen from selected samples of the last glacial period and the early to mid-Holocene. Classification and regression tree (CART) analysis was used to establish a model from the modern dataset that discriminates pollen from the different pine species and allows identification of fossil pine pollen at the species level. Results The CART model was effective in separating pollen of P. nigra and P. sylvestris from that of the Mediterranean pine group (P. halepensis, P. pinea and P. pinaster). The pollen of Pinus nigra diverged from that of P. sylvestris by having a more flattened corpus. Predictions using this model suggested that fossil pine pollen is mainly from P. nigra in all the samples analysed. Pinus sylvestris was more abundant in samples from Greenland stadials than Heinrich stadials, whereas Mediterranean pines increased in samples from Greenland interstadials and during the early to mid-Holocene. Main conclusions Morphological parameters can be successfully used to increase the taxonomic resolution of fossil pine pollen at the species level for the highland pines (P. nigra and P. sylvestris) and at the group of species level for the Mediterranean pines. Our study indicates that P. nigra was the dominant component of the last glacial south-western/central Iberian pinewoods, although the species composition of these woodlands varied in response to abrupt climate changes.

Terrain model of the eastern slope of the Porcupine Seabight (50 m grid size)

The topography of the eastern margin of the Porcupine Seabight was surveyed in June 2000 utilizing swath bathymetry. The survey was carried out during RV Polarstern cruise ANT XVII/4 as part of the GEOMOUND project. The main objective was to map and investigate the seafloor topography of this region. The investigated area contains a variability of morphological features such as deep sea channels and giant mounds. The survey was planned and realized on the basis of existing data so as to guarantee the complete coverage of the margin. In order to achieve a resolution of the final digital terrain model (DTM) that meets the project demands, data processing was adjusted accordingly. The grid spacing of the DTM was set to 50 m and an accuracy better than 1% of the water depth was achieved for 96% of the soundings.

Terrain model of the Håkon Mosby Mud Volcano (10 m grid size)

The Håkon Mosby Mud Volcano is a natural laboratory to study geological, geochemical, and ecological processes related to deep-water mud volcanism. High resolution bathymetry of the Håkon Mosby Mud Volcano was recorded during RV Polarstern expedition ARK-XIX/3 utilizing the multibeam system Hydrosweep DS-2. Dense spacing of the survey lines and slow ship speed (5 knots) provided necessary point density to generate a regular 10 m grid. Generalization was applied to preserve and represent morphological structures appropriately. Contour lines were derived showing detailed topography at the centre of the Håkon Mosby Mud Volcano and generalized contours in the vicinity. We provide a brief introduction to the Håkon Mosby Mud Volcano area and describe in detail data recording and processing methods, as well as the morphology of the area. Accuracy assessment was made to evaluate the reliability of a 10 m resolution terrain model. Multibeam sidescan data were recorded along with depth measurements and show reflectivity variations from light grey values at the centre of the Håkon Mosby Mud Volcano to dark grey values (less reflective) at the surrounding moat.

Medem Channel profiles based on stand-alone model simulation and corrections with nudging

State-of-the-art process-based models have shown to be applicable to the simulation and prediction of coastal morphodynamics. On annual to decadal temporal scales, these models may show limitations in reproducing complex natural morphological evolution patterns, such as the movement of bars and tidal channels, e.g. the observed decadal migration of the Medem Channel in the Elbe Estuary, German Bight. Here a morphodynamic model is shown to simulate the hydrodynamics and sediment budgets of the domain to some extent, but fails to adequately reproduce the pronounced channel migration, due to the insufficient implementation of bank erosion processes. In order to allow for long-term simulations of the domain, a nudging method has been introduced to update the model-predicted bathymetries with observations. The model-predicted bathymetry is nudged towards true states in annual time steps. Sensitivity analysis of a user-defined correlation length scale, for the definition of the background error covariance matrix during the nudging procedure, suggests that the optimal error correlation length is similar to the grid cell size, here 80-90 m. Additionally, spatially heterogeneous correlation lengths produce more realistic channel depths than do spatially homogeneous correlation lengths. Consecutive application of the nudging method compensates for the (stand-alone) model prediction errors and corrects the channel migration pattern, with a Brier skill score of 0.78. The proposed nudging method in this study serves as an analytical approach to update model predictions towards a predefined 'true' state for the spatiotemporal interpolation of incomplete morphological data in long-term simulations.