Geometry-Aware Multiscale Image Registration via OBBTree-Based Polyaffine Log-Demons

Non-linear image registration is an important tool in many areas of image analysis. For instance, in morphometric studies of a population of brains, free-form deformations between images are analyzed to describe the structural anatomical variability. Such a simple deformation model is justified by the absence of an easy expressible prior about the shape changes. Applying the same algorithms used in brain imaging to orthopedic images might not be optimal due to the difference in the underlying prior on the inter-subject deformations. In particular, using an un-informed deformation prior often leads to local minima far from the expected solution. To improve robustness and promote anatomically meaningful deformations, we propose a locally affine and geometry-aware registration algorithm that automatically adapts to the data. We build upon the log-domain demons algorithm and introduce a new type of OBBTree-based regularization in the registration with a natural multiscale structure. The regularization model is composed of a hierarchy of locally affine transformations via their logarithms. Experiments on mandibles show improved accuracy and robustness when used to initialize the demons, and even similar performance by direct comparison to the demons, with a significantly lower degree of freedom. This closes the gap between polyaffine and non-rigid registration and opens new ways to statistically analyze the registration results.

A MICROFLUIDIC METHOD TO MEASURE DIFFUSION IN HYDROGELS

A novel microfluidic method is proposed for studying diffusion of small molecules in a hydrogel. Microfluidic devices were prepared with semi-permeable microchannels defined by crosslinked poly(ethylene glycol) (PEG). Uptake of dye molecules from aqueous solutions flowing through the microchannels was observedoptically and diffusion of the dye into the hydrogel was quantified. To complement the diffusion measurements from the microfluidic studies, nuclear magnetic resonance(NMR) characterization of the diffusion of dye in the PEG hydrogels was performed. The diffusion of small molecules in a hydrogel is relevant to applications such asdrug delivery and modeling transport for tissue-engineering applications. The diffusion of small molecules in a hydrogel is dependent on the extent of crosslinking within the gel, gel structure, and interactions between the diffusive species and the hydrogel network. These effects were studied in a model environment (semi-infinite slab) at the hydrogelfluid boundary in a microfluidic device. The microfluidic devices containing PEG microchannels were fabricated using photolithography. The unsteady diffusion of small molecules (dyes) within the microfluidic device was monitored and recorded using a digital microscope. The information was analyzed with techniques drawn from digital microscopy and image analysis to obtain concentration profiles with time. Using a diffusion model to fit this concentration vs. position data, a diffusion coefficient was obtained. This diffusion coefficient was compared to those from complementary NMR analysis. A pulsed field gradient (PFG) method was used to investigate and quantify small molecule diffusion in gradient (PFG) method was used to investigate and quantify small molecule diffusion in hydrogels. There is good agreement between the diffusion coefficients obtained from the microfluidic methods and those found from the NMR studies. The microfluidic approachused in this research enables the study of diffusion at length scales that approach those of vasculature, facilitating models for studying drug elution from hydrogels in blood-contacting applications.

Brain areas involved in medial temporal lobe seizures: a principal component analysis of ictal SPECT data

The study describes brain areas involved in medial temporal lobe (mTL) seizures of 12 patients. All patients showed so-called oro-alimentary behavior within the first 20 s of clinical seizure manifestation characteristic of mTL seizures. Single photon emission computed tomography (SPECT) images of regional cerebral blood flow (rCBF) were acquired from the patients in ictal and interictal phases and from normal volunteers. Image analysis employed categorical comparisons with statistical parametric mapping and principal component analysis (PCA) to assess functional connectivity. PCA supplemented the findings of the categorical analysis by decomposing the covariance matrix containing images of patients and healthy subjects into distinct component images of independent variance, including areas not identified by the categorical analysis. Two principal components (PCs) discriminated the subject groups: patients with right or left mTL seizures and normal volunteers, indicating distinct neuronal networks implicated by the seizure. Both PCs were correlated with seizure duration, one positively and the other negatively, confirming their physiological significance. The independence of the two PCs yielded a clear clustering of subject groups. The local pattern within the temporal lobe describes critical relay nodes which are the counterpart of oro-alimentary behavior: (1) right mesial temporal zone and ipsilateral anterior insula in right mTL seizures, and (2) temporal poles on both sides that are densely interconnected by the anterior commissure. Regions remote from the temporal lobe may be related to seizure propagation and include positively and negatively loaded areas. These patterns, the covarying areas of the temporal pole and occipito-basal visual association cortices, for example, are related to known anatomic paths.

An introduction to low-level analysis methods of DNA microarray data

This article gives an overview over the methods used in the low--level analysis of gene expression data generated using DNA microarrays. This type of experiment allows to determine relative levels of nucleic acid abundance in a set of tissues or cell populations for thousands of transcripts or loci simultaneously. Careful statistical design and analysis are essential to improve the efficiency and reliability of microarray experiments throughout the data acquisition and analysis process. This includes the design of probes, the experimental design, the image analysis of microarray scanned images, the normalization of fluorescence intensities, the assessment of the quality of microarray data and incorporation of quality information in subsequent analyses, the combination of information across arrays and across sets of experiments, the discovery and recognition of patterns in expression at the single gene and multiple gene levels, and the assessment of significance of these findings, considering the fact that there is a lot of noise and thus random features in the data. For all of these components, access to a flexible and efficient statistical computing environment is an essential aspect.

Quantitative analysis of OCT characteristics in patients with achromatopsia and blue-cone monochromatism

PURPOSE: To quantify optical coherence tomography (OCT) images of the central retina in patients with blue-cone monochromatism (BCM) and achromatopsia (ACH) compared with healthy control individuals. METHODS: The study included 15 patients with ACH, 6 with BCM, and 20 control subjects. Diagnosis of BCM and ACH was established by visual acuity testing, morphologic examination, color vision testing, and Ganzfeld ERG recording. OCT images were acquired with the Stratus OCT 3 (Carl Zeiss Meditec AG, Oberkochen, Germany). Foveal OCT images were analyzed by calculating longitudinal reflectivity profiles (LRPs) from scan lines. Profiles were analyzed quantitatively to determine foveal thickness and distances between reflectivity layers. RESULTS: Patients with ACH and BCM had a mean visual acuity of 20/200 and 20/60, respectively. Color vision testing results were characteristic of the diseases. The LRPs of control subjects yielded four peaks (P1-P4), presumably representing the RPE (P1), the ovoid region of the photoreceptors (P2), the external limiting membrane (ELM) (P3), and the internal limiting membrane (P4). In patients with ACH, P2 was absent, but foveal thickness (P1-P4) did not differ significantly from that in the control subjects (187 +/- 20 vs. 192 +/- 14 microm, respectively). The distance from P1 to P3 did not differ significantly (78 +/- 10 vs. 82 +/- 5 microm) between ACH and controls subjects. In patients with BCM, P3 was lacking, and P2 advanced toward P1 compared with the control subjects (32 +/- 6 vs. 48 +/- 4 microm). Foveal thickness (153 +/- 16 microm) was significantly reduced compared with that in control subjects and patients with ACH. CONCLUSIONS: Quantitative OCT image analysis reveals distinct patterns for controls subjects and patients with ACH and BCM, respectively. Quantitative analysis of OCT imaging can be useful in differentiating retinal diseases affecting photoreceptors. Foveal thickness is similar in both normal subjects and patients with ACH but is decreased in patients with BCM.

BEYOND ROOTS ALONE: NOVEL METHODOLOGIES FOR ANALYZING COMPLEX SOIL AND MINIRHIZOTRON IMAGERY USING IMAGE PROCESSING AND GIS TOOLS

Quantifying belowground dynamics is critical to our understanding of plant and ecosystem function and belowground carbon cycling, yet currently available tools for complex belowground image analyses are insufficient. We introduce novel techniques combining digital image processing tools and geographic information systems (GIS) analysis to permit semi-automated analysis of complex root and soil dynamics. We illustrate methodologies with imagery from microcosms, minirhizotrons, and a rhizotron, in upland and peatland soils. We provide guidelines for correct image capture, a method that automatically stitches together numerous minirhizotron images into one seamless image, and image analysis using image segmentation and classification in SPRING or change analysis in ArcMap. These methods facilitate spatial and temporal root and soil interaction studies, providing a framework to expand a more comprehensive understanding of belowground dynamics.

Computer-aided diagnosis of carotid atherosclerosis based on ultrasound image statistics, laws' texture and neural networks

Quantitative characterisation of carotid atherosclerosis and classification into symptomatic or asymptomatic is crucial in planning optimal treatment of atheromatous plaque. The computer-aided diagnosis (CAD) system described in this paper can analyse ultrasound (US) images of carotid artery and classify them into symptomatic or asymptomatic based on their echogenicity characteristics. The CAD system consists of three modules: a) the feature extraction module, where first-order statistical (FOS) features and Laws' texture energy can be estimated, b) the dimensionality reduction module, where the number of features can be reduced using analysis of variance (ANOVA), and c) the classifier module consisting of a neural network (NN) trained by a novel hybrid method based on genetic algorithms (GAs) along with the back propagation algorithm. The hybrid method is able to select the most robust features, to adjust automatically the NN architecture and to optimise the classification performance. The performance is measured by the accuracy, sensitivity, specificity and the area under the receiver-operating characteristic (ROC) curve. The CAD design and development is based on images from 54 symptomatic and 54 asymptomatic plaques. This study demonstrates the ability of a CAD system based on US image analysis and a hybrid trained NN to identify atheromatous plaques at high risk of stroke.

Do mature pulmonary lobes grow after transplantation into an immature recipient?

BACKGROUND The use of reduced-size adult lung transplants could help solve the profound pediatric donor lung shortage. However, adequate long-term function of the mature grafts requires growth in proportion to the recipient's development. METHODS Mature left lower lobes from adult mini-pigs (age: 7 months; mean body weight: 30 kg) were transplanted into 14-week-old piglets (mean body weight: 15 kg). By the end of the 14-week holding period, lungs of the recipients (n = 4) were harvested. After volumetric measurements, the lung morphology was studied using light microscopy, scanning, and transmission electron microscopy. Changes of alveolar airspace volume were determined using a computer aided image analysis system. Comparisons were made to age- and weight-matched controls. RESULTS Volumetric studies showed no significant differences (p = 0.49) between the specific volume (mL/kg body weight) of lobar grafts and left lower lobes of adult controls. Morphologic studies showed marked structural differences between the grafts and the right native lungs of the recipients, with increased average alveolar diameter of the grafts. On light microscopy and scanning electron microscopy, alveoli appeared dilated and rounded compared to the normal polygonal shape in the controls. The computer generated semi-quantitative data of relative alveolar airspace volume tended to be higher in transplanted lobes. CONCLUSIONS The mature pulmonary lobar grafts have filled the growing left hemithorax of the developing recipient. Emphysema-like alterations of the grafts were observed without evidence of alveolar growth in the mature lobar transplants. Thus, it can be questioned whether mature pulmonary grafts can guarantee sufficient long-term gas exchange in growing recipients.

Intact flagellar motor of Borrelia burgdorferi revealed by cryo-electron tomography: evidence for stator ring curvature and rotor/C-ring assembly flexion.

The bacterial flagellar motor is a remarkable nanomachine that provides motility through flagellar rotation. Prior structural studies have revealed the stunning complexity of the purified rotor and C-ring assemblies from flagellar motors. In this study, we used high-throughput cryo-electron tomography and image analysis of intact Borrelia burgdorferi to produce a three-dimensional (3-D) model of the in situ flagellar motor without imposing rotational symmetry. Structural details of B. burgdorferi, including a layer of outer surface proteins, were clearly visible in the resulting 3-D reconstructions. By averaging the 3-D images of approximately 1,280 flagellar motors, a approximately 3.5-nm-resolution model of the stator and rotor structures was obtained. flgI transposon mutants lacked a torus-shaped structure attached to the flagellar rod, establishing the structural location of the spirochetal P ring. Treatment of intact organisms with the nonionic detergent NP-40 resulted in dissolution of the outermost portion of the motor structure and the C ring, providing insight into the in situ arrangement of the stator and rotor structures. Structural elements associated with the stator followed the curvature of the cytoplasmic membrane. The rotor and the C ring also exhibited angular flexion, resulting in a slight narrowing of both structures in the direction perpendicular to the cell axis. These results indicate an inherent flexibility in the rotor-stator interaction. The FliG switching and energizing component likely provides much of the flexibility needed to maintain the interaction between the curved stator and the relatively symmetrical rotor/C-ring assembly during flagellar rotation.

Do binucleate cardiomyocytes have a role in myocardial repair? Insights using isolated rodent myocytes and cell culture.

Neonatal and adult cardiomyocytes were isolated from rat hearts. Some of the adult myocytes were cultured to allow for cell dedifferentiation, a phenomenon thought to mimic cell changes that occur in stressed myocardium, with myocytes regressing to a fetal pattern of metabolism and stellate neonatal shape.Using fluorescence deconvolution microscopy, cells were probed with fluorescent markers and scanned for a number of proteins associated with ion control, calcium movements and cardiac function. Image analysis of deconvoluted image stacks and sequential real-time image recordings of calcium transients of cells were made.All three myocyte groups were predominantly comprised of binucleate cells. Clustering of proteins to a single nucleus was a common observation, suggesting that one nucleus is active in protein synthesis pathways, while the other nucleus assumes a 'dormant' or different role and that cardiomyocytes might be mitotically active even in late development, or specific protein syntheses could be targeted and regulated for reintroduction into the cell cycle.Such possibilities would extend cardiac disease associated stem cell research and therapy options, while producing valuable insights into developmental and death pathways of binucleate cardiomyocytes (word count 183).

Mouse models and quantitative trait loci of susceptibility to bleomycin- and radiation-induced pulmonary fibrosis

Radiotherapy involving the thoracic cavity and chemotherapy with the drug bleomycin are both dose limited by the development of pulmonary fibrosis. From evidence that there is variation in the population in susceptibility to pulmonary fibrosis, and animal data, it was hypothesized that individual variation in susceptibility to bleomycin-induced, or radiation-induced, pulmonary fibrosis is, in part, genetically controlled. In this thesis a three generation mouse genetic model of C57BL/6J (fibrosis prone) and C3Hf/Kam (fibrosis resistant) mouse strains and F1 and F2 (F1 intercross) progeny derived from the parental strains was developed to investigate the genetic basis of susceptibility to fibrosis. In the bleomycin studies the mice received 100 mg/kg (125 for females) of bleomycin, via mini osmotic pump. The animals were sacrificed at eight weeks following treatment or when their breathing rate indicated respiratory distress. In the radiation studies the mice were given a single dose of 14 or 16 Gy (Co$\sp{60})$ to the whole thorax and were sacrificed when moribund. The phenotype was defined as the percent of fibrosis area in the left lung as quantified with image analysis of histological sections. Quantitative trait loci (QTL) mapping was used to identify the chromosomal location of genes which contribute to susceptibility to bleomycin-induced pulmonary fibrosis in C57BL/6J mice compared to C3Hf/Kam mice and to determine if the QTL's which influence susceptibility to bleomycin-induced lung fibrosis in these progenitor strains could be implicated in susceptibility to radiation-induced lung fibrosis. For bleomycin, a genome wide scan revealed QTL's on chromosome 17, at the MHC, (LOD = 11.7 for males and 7.2 for females) accounting for approximately 21% of the phenotypic variance, and on chromosome 11 (LOD = 4.9), in male mice only, adding 8% of phenotypic variance. The bleomycin QTL on chromosome 17 was also implicated for susceptibility to radiation-induced fibrosis (LOD = 5.0) and contributes 7% of the phenotypic variance in the radiation study. In conclusion, susceptibility to both bleomycin-induced and radiation-induced pulmonary fibrosis are heritable traits, and are influenced by a genetic factor which maps to a genomic region containing the MHC. ^

A computer-assisted microscopic analysis of bone tissue developed inside a polyactive plymer implanted into an equine articular surface

One of the most promising applications for the restoration of small or moderately sized focal articular lesions is mosaicplasty (MP). Although recurrent hemarthrosis is a rare complication after MP, recently, various strategies have been designed to find an effective filling material to prevent postoperative bleeding from the donor site. The porous biodegradable polymer Polyactive (PA; a polyethylene glycol terephthalate - polybutylene terephthalate copolymer) represents a promising solution in this respect. A histological evaluation of the longterm PA-filled donor sites obtained from 10 experimental horses was performed. In this study, attention was primarily focused on the bone tissue developed in the plug. A computer-assisted image analysis and quantitative polarized light microscopic measurements of decalcified, longitudinally sectioned, dimethylmethylene blue (DMMB)- and picrosirius red (PS) stained sections revealed that the coverage area of the bone trabecules in the PA-filled donor tunnels was substantially (25%) enlarged compared to the neighboring cancellous bone. For this quantification, identical ROIs (regions of interest) were used and compared. The birefringence retardation values were also measured with a polarized light microscope using monochromatic light. Identical retardation values could be recorded from the bone trabeculae developed in the PA and in the neighboring bone, which indicates that the collagen orientation pattern does not differ significantly among these bone trabecules. Based on our new data, we speculate that PA promotes bone formation, and some of the currently identified degradation products of PA may enhance osteo-conduction and osteoinduction inside the donor canal.

Reconstruction of 3D Vertebral Models from a Single 2D Lateral Fluoroscopic Image

Accurate three-dimensional (3D) models of lumbar vertebrae are required for image-based 3D kinematics analysis. MRI or CT datasets are frequently used to derive 3D models but have the disadvantages that they are expensive, time-consuming or involving ionizing radiation (e.g., CT acquisition). In this chapter, we present an alternative technique that can reconstruct a scaled 3D lumbar vertebral model from a single two-dimensional (2D) lateral fluoroscopic image and a statistical shape model. Cadaveric studies are conducted to verify the reconstruction accuracy by comparing the surface models reconstructed from a single lateral fluoroscopic image to the ground truth data from 3D CT segmentation. A mean reconstruction error between 0.7 and 1.4 mm was found.