Segmentation of brain tumor images based on atlas-registration combined with a Markov-Random-Field lesion growth model

We present an automatic method to segment brain tissues from volumetric MRI brain tumor images. The method is based on non-rigid registration of an average atlas in combination with a biomechanically justified tumor growth model to simulate soft-tissue deformations caused by the tumor mass-effect. The tumor growth model, which is formulated as a mesh-free Markov Random Field energy minimization problem, ensures correspondence between the atlas and the patient image, prior to the registration step. The method is non-parametric, simple and fast compared to other approaches while maintaining similar accuracy. It has been evaluated qualitatively and quantitatively with promising results on eight datasets comprising simulated images and real patient data.

Feature-invariant image registration method for quantification of surgical outcomes in patients with craniosynostosis: a preliminary study

Craniosynostosis consists of a premature fusion of the sutures in an infant skull that restricts skull and brain growth. During the last decades, there has been a rapid increase of fundamentally diverse surgical treatment methods. At present, the surgical outcome has been assessed using global variables such as cephalic index, head circumference, and intracranial volume. However, these variables have failed in describing the local deformations and morphological changes that may have a role in the neurologic disorders observed in the patients. This report describes a rigid image registration-based method to evaluate outcomes of craniosynostosis surgical treatments, local quantification of head growth, and indirect intracranial volume change measurements. The developed semiautomatic analysis method was applied to computed tomography data sets of a 5-month-old boy with sagittal craniosynostosis who underwent expansion of the posterior skull with cranioplasty. Quantification of the local changes between pre- and postoperative images was quantified by mapping the minimum distance of individual points from the preoperative to the postoperative surface meshes, and indirect intracranial volume changes were estimated. The proposed methodology can provide the surgeon a tool for the quantitative evaluation of surgical procedures and detection of abnormalities of the infant skull and its development.

Identification of candidate genes and physiological pathways involved in gonad deformation in whitefish (Coregonus spp.) from Lake Thun, Switzerland

In 2000, fishermen reported the appearance of deformed reproductive organs in whitefish (Coregonus spp.) from Lake Thun, Switzerland. Despite intensive investigations, the causes of these abnormalities remain unknown. Using gene expression profiling, we sought to identify candidate genes and physiological processes possibly associated with the observed gonadal deformations, in order to gain insights into potential causes. Using in situ-synthesized oligonucleotide arrays, we compared the expression levels at 21,492 unique transcript probes in liver and head kidney tissue of male whitefish with deformed and normally developed gonads, respectively. The fish had been collected on spawning sites of two genetically distinct whitefish forms of Lake Thun. We contrasted the gene expression profiles of 56 individuals, i.e., 14 individuals of each phenotype and of each population. Gene-by-gene analysis revealed weak expression differences between normal and deformed fish, and only one gene, ictacalcin, was found to be up-regulated in head kidney tissue of deformed fish from both whitefish forms, However, this difference could not be confirmed with quantitative real-time qPCR. Enrichment analysis on the level of physiological processes revealed (i) the involvement of immune response genes in both tissues, particularly those linked to complement activation in the liver, (ii) proteolysis in the liver and (iii) GTPase activity and Ras protein signal transduction in the head kidney. In comparison with current literature, this gene expression pattern signals a chronic autoimmune disease in the testes. Based on the recent observations that gonad deformations are induced through feeding of zooplankton from Lake Thun we hypothesize that a xenobiotic accumulated in whitefish via the plankton triggering autoimmunity as the likely cause of gonad deformations. We propose several experimental strategies to verify or reject this hypothesis.

Simultaneous Multiscale Polyaffine Registration by Incorporating Deformation Statistics

Locally affine (polyaffine) image registration methods capture intersubject non-linear deformations with a low number of parameters, while providing an intuitive interpretation for clinicians. Considering the mandible bone, anatomical shape differences can be found at different scales, e.g. left or right side, teeth, etc. Classically, sequential coarse to fine registration are used to handle multiscale deformations, instead we propose a simultaneous optimization of all scales. To avoid local minima we incorporate a prior on the polyaffine transformations. This kind of groupwise registration approach is natural in a polyaffine context, if we assume one configuration of regions that describes an entire group of images, with varying transformations for each region. In this paper, we reformulate polyaffine deformations in a generative statistical model, which enables us to incorporate deformation statistics as a prior in a Bayesian setting. We find optimal transformations by optimizing the maximum a posteriori probability. We assume that the polyaffine transformations follow a normal distribution with mean and concentration matrix. Parameters of the prior are estimated from an initial coarse to fine registration. Knowing the region structure, we develop a blockwise pseudoinverse to obtain the concentration matrix. To our knowledge, we are the first to introduce simultaneous multiscale optimization through groupwise polyaffine registration. We show results on 42 mandible CT images.

Prevalence of keel bone deformities in Swiss laying hens

The goal of this study was to evaluate the prevalence of keel bone deformities of laying hens in Switzerland. The keel bones of 100 end-of-lay hens from each of 39 flocks (3900 in total) were palpated. On average, 25.4% of the hens had moderately or severely deformed keel bones and the overall prevalence including slight deformities was 55%. 3. Variation between flocks was considerable. Thus, the prevalence of moderately or severely deformed keel bones ranged from 6 to 48%, and the overall prevalence including slight deformities ranged from 20 to 83%. Aviary housing was associated with a higher prevalence of total, and severe or moderate deformations, compared with floor pens. There were no significant differences in the number of deformities between the different plumage colours, hybrids or perch materials.

Generation of custom modes in a Nd:YAG laser with a semipassive bimorph adaptive mirror

Custom modes at a wavelength of 1064 nm were generated with a deformable mirror. The required surface deformations of the adaptive mirror were calculated with the Collins integral written in a matrix formalism. The appropriate size and shape of the actuators as well as the needed stroke were determined to ensure that the surface of the controllable mirror matches the phase front of the custom modes. A semipassive bimorph adaptive mirror with five concentric ring-shaped actuators and one defocus actuator was manufactured and characterised. The surface deformation was modelled with the response functions of the adaptive mirror in terms of an expansion with Zernike polynomials. In the experiments the Nd:YAG laser crystal was quasi-CW pumped to avoid thermally induced distortions of the phase front. The adaptive mirror allows to switch between a super-Gaussian mode, a doughnut mode, a Hermite-Gaussian fundamental beam, multi-mode operation or no oscillation in real time during laser operation.

Chondrocyte mechanotransduction: effects of compression on deformation of intracellular organelles and relevance to cellular biosynthesis

OBJECTIVE: The effects of mechanical deformation of intact cartilage tissue on chondrocyte biosynthesis in situ have been well documented, but the mechanotransduction pathways that regulate such phenomena have not been elucidated completely. The goal of this study was to examine the effects of tissue deformation on the morphology of a range of intracellular organelles which play a major role in cell biosynthesis and metabolism. DESIGN: Using chemical fixation, high pressure freezing, and electron microscopy, we imaged chondrocytes within mechanically compressed cartilage explants at high magnification and quantitatively and qualitatively assessed changes in organelle volume and shape caused by graded levels of loading. RESULTS: Compression of the tissue caused a concomitant reduction in the volume of the extracellular matrix (ECM), chondrocyte, nucleus, rough endoplasmic reticulum, and mitochondria. Interestingly, however, the Golgi apparatus was able to resist loss of intraorganelle water and retain a portion of its volume relative to the remainder of the cell. These combined results suggest that a balance between intracellular mechanical and osmotic gradients govern the changes in shape and volume of the organelles as the tissue is compressed. CONCLUSIONS: Our results lead to the interpretive hypothesis that organelle volume changes appear to be driven mainly by osmotic interactions while shape changes are mediated by structural factors, such as cytoskeletal interactions that may be linked to extracellular matrix deformations. The observed volume and shape changes of the chondrocyte organelles and the differential behavior between organelles during tissue compression provide evidence for an important mechanotransduction pathway linking translational and post-translational events (e.g., elongation and sulfation of glycosaminoglycans (GAGs) in the Golgi) to cell deformation.

The etiology of cleft palate formation in BMP7-deficient mice

Palatogenesis is a complex process implying growth, elevation and fusion of the two lateral palatal shelves during embryogenesis. This process is tightly controlled by genetic and mechanistic cues that also coordinate the growth of other orofacial structures. Failure at any of these steps can result in cleft palate, which is a frequent craniofacial malformation in humans. To understand the etiology of cleft palate linked to the BMP signaling pathway, we studied palatogenesis in Bmp7-deficient mouse embryos. Bmp7 expression was found in several orofacial structures including the edges of the palatal shelves prior and during their fusion. Bmp7 deletion resulted in a general alteration of oral cavity morphology, unpaired palatal shelf elevation, delayed shelf approximation, and subsequent lack of fusion. Cell proliferation and expression of specific genes involved in palatogenesis were not altered in Bmp7-deficient embryos. Conditional ablation of Bmp7 with Keratin14-Cre or Wnt1-Cre revealed that neither epithelial nor neural crest-specific loss of Bmp7 alone could recapitulate the cleft palate phenotype. Palatal shelves from mutant embryos were able to fuse when cultured in vitro as isolated shelves in proximity, but not when cultured as whole upper jaw explants. Thus, deformations in the oral cavity of Bmp7-deficient embryos such as the shorter and wider mandible were not solely responsible for cleft palate formation. These findings indicate a requirement for Bmp7 for the coordination of both developmental and mechanistic aspects of palatogenesis.

Topological lattice actions for the 2d XY model

We consider the 2d XY Model with topological lattice actions, which are invariant against small deformations of the field configuration. These actions constrain the angle between neighbouring spins by an upper bound, or they explicitly suppress vortices (and anti-vortices). Although topological actions do not have a classical limit, they still lead to the universal behaviour of the Berezinskii-Kosterlitz-Thouless (BKT) phase transition — at least up to moderate vortex suppression. In the massive phase, the analytically known Step Scaling Function (SSF) is reproduced in numerical simulations. However, deviations from the expected universal behaviour of the lattice artifacts are observed. In the massless phase, the BKT value of the critical exponent ηc is confirmed. Hence, even though for some topological actions vortices cost zero energy, they still drive the standard BKT transition. In addition we identify a vortex-free transition point, which deviates from the BKT behaviour.

Heterotic T-folds with a small number of neutral moduli

We discuss non-geometric supersymmetric heterotic string models in D=4, in the framework of the free fermionic construction. We perform a systematic scan of models with four a priori left-right asymmetric Z2 projections and shifts. We analyze some 220 models, identifying 18 inequivalent classes and addressing variants generated by discrete torsions. They do not contain geometrical or trivial neutral moduli, apart from the dilaton. However, we show the existence of flat directions in the form of exactly marginal deformations and identify patterns of symmetry breaking where product gauge groups, realized at level one, are broken to their diagonal at higher level. We also describe an “inverse Gepner map” from Heterotic to Type II models that could be used, in certain non geometric settings, to define “effective” topological invariants.

Quantification of Popliteal Artery Deformation During Leg Flexion in Subjects With Peripheral Artery Disease: A Pilot Study

Purpose: To quantify the in vivo deformations of the popliteal artery during leg flexion in subjects with clinically relevant peripheral artery disease (PAD). Methods: Five patients (4 men; mean age 69 years, range 56–79) with varying calcification levels of the popliteal artery undergoing endovascular revascularization underwent 3-dimensional (3D) rotational angiography. Image acquisition was performed with the leg straight and with a flexion of 70°/20° in the knee/hip joints. The arterial centerline and the corresponding branches in both positions were segmented to create 3D reconstructions of the arterial trees. Axial deformation, twisting, and curvatures were quantified. Furthermore, the relationships between the calcification levels and the deformations were investigated. Results: An average shortening of 5.9%±2.5% and twist rate of 3.8±2.2°/cm in the popliteal artery were observed. Maximal curvatures in the straight and flexed positions were 0.12±0.04 cm−1 and 0.24±0.09 cm−1, respectively. As the severity of calcification increased, the maximal curvature in the straight position increased from 0.08 to 0.17 cm−1, while an increase from 0.17 to 0.39 cm−1 was observed for the flexed position. Axial elongations and arterial twisting were not affected by the calcification levels. Conclusion: The popliteal artery of patients with symptomatic PAD is exposed to significant deformations during flexion of the knee joint. The severity of calcification directly affects curvature, but not arterial length or twisting angles. This pilot study also showed the ability of rotational angiography to quantify the 3D deformations of the popliteal artery in patients with various levels of calcification.

Integrated Segmentation of Brain Tumor Images for Radiotherapy and Neurosurgery

Image-based modeling of tumor growth combines methods from cancer simulation and medical imaging. In this context, we present a novel approach to adapt a healthy brain atlas to MR images of tumor patients. In order to establish correspondence between a healthy atlas and a pathologic patient image, tumor growth modeling in combination with registration algorithms is employed. In a first step, the tumor is grown in the atlas based on a new multi-scale, multi-physics model including growth simulation from the cellular level up to the biomechanical level, accounting for cell proliferation and tissue deformations. Large-scale deformations are handled with an Eulerian approach for finite element computations, which can operate directly on the image voxel mesh. Subsequently, dense correspondence between the modified atlas and patient image is established using nonrigid registration. The method offers opportunities in atlasbased segmentation of tumor-bearing brain images as well as for improved patient-specific simulation and prognosis of tumor progression.

O(N) Models with Topological Lattice Actions

A variety of lattice discretisations of continuum actions has been considered, usually requiring the correct classical continuum limit. Here we discuss “weird” lattice formulations without that property, namely lattice actions that are invariant under most continuous deformations of the field configuration, in one version even without any coupling constants. It turns out that universality is powerful enough to still provide the correct quantum continuum limit, despite the absence of a classical limit, or a perturbative expansion. We demonstrate this for a set of O(N) models (or non-linear σ-models). Amazingly, such “weird” lattice actions are not only in the right universality class, but some of them even have practical benefits, in particular an excellent scaling behaviour.

Optical Tolerance Analysis of the Multi-Beam Limb Viewing Instrument STEAMR

In this paper, we report on an optical tolerance analysis of the submillimeter atmospheric multi-beam limb sounder, STEAMR. Physical optics and ray-tracing methods were used to quantify and separate errors in beam pointing and distortion due to reflector misalignment and primary reflector surface deformations. Simulations were performed concurrently with the manufacturing of a multi-beam demonstrator of the relay optical system which shapes and images the beams to their corresponding receiver feed horns. Results from Monte Carlo simulations show that the inserts used for reflector mounting should be positioned with an overall accuracy better than 100 μm (~ 1/10 wavelength). Analyses of primary reflector surface deformations show that a deviation of magnitude 100 μm can be tolerable before deployment, whereas the corresponding variations should be less than 30 μm during operation. The most sensitive optical elements in terms of misalignments are found near the focal plane. This localized sensitivity is attributed to the off-axis nature of the beams at this location. Post-assembly mechanical measurements of the reflectors in the demonstrator show that alignment better than 50 μm could be obtained.

FEM/FD Immersed Boundary FSI Simulations

Immersed boundary simulations have been under development for physiological flows, allowing for elegant handling of fluid-structure interaction modelling with large deformations due to retained domain-specific meshing. We couple a structural system in Lagrangian representation that is formulated in a weak form with a Navier-Stokes system discretized through a finite differences scheme. We build upon a proven highly scalable imcompressible flow solver that we extend to handle FSI. We aim at applying our method to investigating the hemodynamics of Aortic Valves. The code is going to be extended to conform to the new hybrid-node supercomputers.