Coupling biomechanics to a cellular level model: An approach to patient-specific image driven multi-scale and multi-physics tumor simulation

Modeling of tumor growth has been performed according to various approaches addressing different biocomplexity levels and spatiotemporal scales. Mathematical treatments range from partial differential equation based diffusion models to rule-based cellular level simulators, aiming at both improving our quantitative understanding of the underlying biological processes and, in the mid- and long term, constructing reliable multi-scale predictive platforms to support patient-individualized treatment planning and optimization. The aim of this paper is to establish a multi-scale and multi-physics approach to tumor modeling taking into account both the cellular and the macroscopic mechanical level. Therefore, an already developed biomodel of clinical tumor growth and response to treatment is self-consistently coupled with a biomechanical model. Results are presented for the free growth case of the imageable component of an initially point-like glioblastoma multiforme tumor. The composite model leads to significant tumor shape corrections that are achieved through the utilization of environmental pressure information and the application of biomechanical principles. Using the ratio of smallest to largest moment of inertia of the tumor material to quantify the effect of our coupled approach, we have found a tumor shape correction of 20\% by coupling biomechanics to the cellular simulator as compared to a cellular simulation without preferred growth directions. We conclude that the integration of the two models provides additional morphological insight into realistic tumor growth behavior. Therefore, it might be used for the development of an advanced oncosimulator focusing on tumor types for which morphology plays an important role in surgical and/or radio-therapeutic treatment planning.

Pathology hinting as the combination of automatic segmentation with a statistical shape model

With improvements in acquisition speed and quality, the amount of medical image data to be screened by clinicians is starting to become challenging in the daily clinical practice. To quickly visualize and find abnormalities in medical images, we propose a new method combining segmentation algorithms with statistical shape models. A statistical shape model built from a healthy population will have a close fit in healthy regions. The model will however not fit to morphological abnormalities often present in the areas of pathologies. Using the residual fitting error of the statistical shape model, pathologies can be visualized very quickly. This idea is applied to finding drusen in the retinal pigment epithelium (RPE) of optical coherence tomography (OCT) volumes. A segmentation technique able to accurately segment drusen in patients with age-related macular degeneration (AMD) is applied. The segmentation is then analyzed with a statistical shape model to visualize potentially pathological areas. An extensive evaluation is performed to validate the segmentation algorithm, as well as the quality and sensitivity of the hinting system. Most of the drusen with a height of 85.5 microm were detected, and all drusen at least 93.6 microm high were detected.

Exogenous surfactant application in a rat lung ischemia reperfusion injury model: effects on edema formation and alveolar type II cells

BACKGROUND: Prophylactic exogenous surfactant therapy is a promising way to attenuate the ischemia and reperfusion (I/R) injury associated with lung transplantation and thereby to decrease the clinical occurrence of acute lung injury and acute respiratory distress syndrome. However, there is little information on the mode by which exogenous surfactant attenuates I/R injury of the lung. We hypothesized that exogenous surfactant may act by limiting pulmonary edema formation and by enhancing alveolar type II cell and lamellar body preservation. Therefore, we investigated the effect of exogenous surfactant therapy on the formation of pulmonary edema in different lung compartments and on the ultrastructure of the surfactant producing alveolar epithelial type II cells. METHODS: Rats were randomly assigned to a control, Celsior (CE) or Celsior + surfactant (CE+S) group (n = 5 each). In both Celsior groups, the lungs were flush-perfused with Celsior and subsequently exposed to 4 h of extracorporeal ischemia at 4 degrees C and 50 min of reperfusion at 37 degrees C. The CE+S group received an intratracheal bolus of a modified natural bovine surfactant at a dosage of 50 mg/kg body weight before flush perfusion. After reperfusion (Celsior groups) or immediately after sacrifice (Control), the lungs were fixed by vascular perfusion and processed for light and electron microscopy. Stereology was used to quantify edematous changes as well as alterations of the alveolar epithelial type II cells. RESULTS: Surfactant treatment decreased the intraalveolar edema formation (mean (coefficient of variation): CE: 160 mm3 (0.61) vs. CE+S: 4 mm3 (0.75); p < 0.05) and the development of atelectases (CE: 342 mm3 (0.90) vs. CE+S: 0 mm3; p < 0.05) but led to a higher degree of peribronchovascular edema (CE: 89 mm3 (0.39) vs. CE+S: 268 mm3 (0.43); p < 0.05). Alveolar type II cells were similarly swollen in CE (423 microm3(0.10)) and CE+S (481 microm3(0.10)) compared with controls (323 microm3(0.07); p < 0.05 vs. CE and CE+S). The number of lamellar bodies was increased and the mean lamellar body volume was decreased in both CE groups compared with the control group (p < 0.05). CONCLUSION: Intratracheal surfactant application before I/R significantly reduces the intraalveolar edema formation and development of atelectases but leads to an increased development of peribronchovascular edema. Morphological changes of alveolar type II cells due to I/R are not affected by surfactant treatment. The beneficial effects of exogenous surfactant therapy are related to the intraalveolar activity of the exogenous surfactant.

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.

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.

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.

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.

Type A aortic dissection model to improve endovascular research and technologies.

OBJECTIVE Type A aortic dissection is a life-threatening disease requiring immediate surgical treatment. With emerging catheter-based technologies, endovascular stent-graft implantation to treat aneurysms and dissections has become a standardized procedure. However, endovascular treatment of the ascending aorta remains challenging. Thus we designed an ascending aortic dissection model to allow simulation of endovascular treatment. METHODS Five formalin-fixed human aortas were prepared. The ascending aorta was opened semicircularly in the middle portion and the medial layer was separated from the intima. The intimal tube was readapted using running monofilament sutures. The preparations were assessed by 128-slice computed tomography. A bare-metal stent was implanted for thoracic endovascular aortic repair in 4 of the aortic dissection models. RESULTS Separation of the intimal and medial layer of the aorta was considered to be sufficient because computed tomography showed a clear image of the dissection membrane in each aorta. The dissection was located 3.9 ± 1.4 cm proximally from the aortic annulus, with a length of 4.6 ± 0.9 cm. Before stent implantation, the mean distance from the intimal flap to the aortic wall was measured as 0.63 ± 0.163 cm in the ascending aorta. After stent implantation, this distance decreased to 0.26 ± 0.12 cm. CONCLUSION This model of aortic dissection of the ascending human aorta was reproducible with a comparable pathological and morphological appearance. The technique and model can be used to evaluate new stent-graft technologies to treat type A dissection and facilitate training for surgeons.