Mophological analysis of the human vasculature from medical imaging. Application to cerebral and coranary arteries

In this project, we have investigated new ways of modelling and analysis of human vasculature from Medical images. The research was divided in two main areas: cerebral vasculature analysis and coronary arteries modeling. Regarding cerebral vasculature analysis, we have studed cerebral aneurysms, internal carotid and the Circle of Willis (CoW). Aneurysms are abnormal vessel enlargements that can rupture causing important cerebral damages or death. The understanding of this pathology, together with its virtual treatment, and image diagnosis and prognosis, includes identification and detailed measurement of the aneurysms. In this context, we have proposed two automatic aneurysm isolation method, to separate the abnormal part of the vessel from the healthy part, to homogenize and speed-up the processing pipeline usually employed to study this pathology, [Cardenes2011TMI, arrabide2011MedPhys]. The results obtained from both methods have been also compared and validatied in [Cardenes2012MBEC]. A second important task here the analysis of the internal carotid [Bogunovic2011Media] and the automatic labelling of the CoW, Bogunovic2011MICCAI, Bogunovic2012TMI]. The second area of research covers the study of coronary arteries, specially coronary bifurcations because there is where the formation of atherosclerotic plaque is more common, and where the intervention is more challenging. Therefore, we proposed a novel modelling method from Computed Tomography Angiography (CTA) images, combined with Conventional Coronary Angiography (CCA), to obtain realistic vascular models of coronary bifurcations, presented in [Cardenes2011MICCAI], and fully validated including phantom experiments in [Cardene2013MedPhys]. The realistic models obtained from this method are being used to simulate stenting procedures, and to investigate the hemodynamic variables in coronary bifurcations in the works submitted in [Morlachi2012, Chiastra2012]. Additionally, another preliminary work has been done to reconstruct the coronary tree from rotational angiography, and published in [Cardenes2012ISBI].

Temporal diffeomorphic free-form deformation: application to motion and strain estimation from 3D echocardiography

This paper presents a new registration algorithm, called Temporal Di eomorphic Free Form Deformation (TDFFD), and its application to motion and strain quanti cation from a sequence of 3D ultrasound (US) images. The originality of our approach resides in enforcing time consistency by representing the 4D velocity eld as the sum of continuous spatiotemporal B-Spline kernels. The spatiotemporal displacement eld is then recovered through forward Eulerian integration of the non-stationary velocity eld. The strain tensor iscomputed locally using the spatial derivatives of the reconstructed displacement eld. The energy functional considered in this paper weighs two terms: the image similarity and a regularization term. The image similarity metric is the sum of squared di erences between the intensities of each frame and a reference one. Any frame in the sequence can be chosen as reference. The regularization term is based on theincompressibility of myocardial tissue. TDFFD was compared to pairwise 3D FFD and 3D+t FFD, bothon displacement and velocity elds, on a set of synthetic 3D US images with di erent noise levels. TDFFDshowed increased robustness to noise compared to these two state-of-the-art algorithms. TDFFD also proved to be more resistant to a reduced temporal resolution when decimating this synthetic sequence. Finally, this synthetic dataset was used to determine optimal settings of the TDFFD algorithm. Subsequently, TDFFDwas applied to a database of cardiac 3D US images of the left ventricle acquired from 9 healthy volunteers and 13 patients treated by Cardiac Resynchronization Therapy (CRT). On healthy cases, uniform strain patterns were observed over all myocardial segments, as physiologically expected. On all CRT patients, theimprovement in synchrony of regional longitudinal strain correlated with CRT clinical outcome as quanti ed by the reduction of end-systolic left ventricular volume at follow-up (6 and 12 months), showing the potential of the proposed algorithm for the assessment of CRT.

Is spinal stenosis assessment dependent on slice orientation? A magnetic resonance imaging study.

INTRODUCTION: Lumbar spinal stenosis (LSS) treatment is based primarily on the clinical criteria providing that imaging confirms radiological stenosis. The radiological measurement more commonly used is the dural sac cross-sectional area (DSCA). It has been recently shown that grading stenosis based on the morphology of the dural sac as seen on axial T2 MRI images, better reflects severity of stenosis than DSCA and is of prognostic value. This radiological prospective study investigates the variability of surface measurements and morphological grading of stenosis for varying degrees of angulation of the T2 axial images relative to the disc space as observed in clinical practice. MATERIALS AND METHODS: Lumbar spine TSE T2 three-dimensional (3D) MRI sequences were obtained from 32 consecutive patients presenting with either suspected spinal stenosis or low back pain. Axial reconstructions using the OsiriX software at 0°, 10°, 20° and 30° relative to the disc space orientation were obtained for a total of 97 levels. For each level, DSCA was digitally measured and stenosis was graded according to the 4-point (A-D) morphological grading by two observers. RESULTS: A good interobserver agreement was found in grade evaluation of stenosis (k = 0.71). DSCA varied significantly as the slice orientation increased from 0° to +10°, +20° and +30° at each level examined (P < 0.0001) (-15 to +32% at 10°, -24 to +143% at 20° and -29 to +231% at 30° of slice orientation). Stenosis definition based on the surface measurements changed in 39 out of the 97 levels studied, whereas the morphology grade was modified only in two levels (P < 0.01). DISCUSSION: The need to obtain continuous slices using the classical 2D MRI acquisition technique entails often at least a 10° slice inclination relative to one of the studied discs. Even at this low angulation, we found a significantly statistical difference between surface changes and morphological grading change. In clinical practice, given the above findings, it might therefore not be necessary to align the axial cuts to each individual disc level which could be more time-consuming than obtaining a single series of axial cuts perpendicular to the middle of the lumbar spine or to the most stenotic level. In conclusion, morphological grading seems to offer an alternative means of assessing severity of spinal stenosis that is little affected by image acquisition technique.

Paediatric cardiac CT examinations: impact of the iterative reconstruction method ASIR on image quality--preliminary findings.

BACKGROUND: Radiation dose exposure is of particular concern in children due to the possible harmful effects of ionizing radiation. The adaptive statistical iterative reconstruction (ASIR) method is a promising new technique that reduces image noise and produces better overall image quality compared with routine-dose contrast-enhanced methods. OBJECTIVE: To assess the benefits of ASIR on the diagnostic image quality in paediatric cardiac CT examinations. MATERIALS AND METHODS: Four paediatric radiologists based at two major hospitals evaluated ten low-dose paediatric cardiac examinations (80 kVp, CTDI(vol) 4.8-7.9 mGy, DLP 37.1-178.9 mGy·cm). The average age of the cohort studied was 2.6 years (range 1 day to 7 years). Acquisitions were performed on a 64-MDCT scanner. All images were reconstructed at various ASIR percentages (0-100%). For each examination, radiologists scored 19 anatomical structures using the relative visual grading analysis method. To estimate the potential for dose reduction, acquisitions were also performed on a Catphan phantom and a paediatric phantom. RESULTS: The best image quality for all clinical images was obtained with 20% and 40% ASIR (p < 0.001) whereas with ASIR above 50%, image quality significantly decreased (p < 0.001). With 100% ASIR, a strong noise-free appearance of the structures reduced image conspicuity. A potential for dose reduction of about 36% is predicted for a 2- to 3-year-old child when using 40% ASIR rather than the standard filtered back-projection method. CONCLUSION: Reconstruction including 20% to 40% ASIR slightly improved the conspicuity of various paediatric cardiac structures in newborns and children with respect to conventional reconstruction (filtered back-projection) alone.

Combined transesophageal and surface MRI provides optimal imaging in aortic atherosclerosis.

OBJECTIVE: Surface magnetic resonance imaging (MRI) for aortic plaque assessment is limited by the trade-off between penetration depth and signal-to-noise ratio (SNR). For imaging the deep seated aorta, a combined surface and transesophageal MRI (TEMRI) technique was developed 1) to determine the individual contribution of TEMRI and surface coils to the combined signal, 2) to measure the signal improvement of a combined surface and TEMRI over surface MRI, and 3) to assess for reproducibility of plaque dimension analysis. METHODS AND RESULTS: In 24 patients six black blood proton-density/T2-weighted fast-spin echo images were obtained using three surface and one TEMRI coil for SNR measurements. Reproducibility of plaque dimensions (combined surface and TEMRI) was measured in 10 patients. TEMRI contributed 68% of the signal in the aortic arch and descending aorta, whereas the overall signal gain using the combined technique was up to 225%. Plaque volume measurements had an intraclass correlation coefficient of as high as 0.97. CONCLUSION: Plaque volume measurements for the quantification of aortic plaque size are highly reproducible for combined surface and TEMRI. The TEMRI coil contributes considerably to the aortic MR signal. The combined surface and TEMRI approach improves aortic signal significantly as compared to surface coils alone. CONDENSED ABSTRACT: Conventional MRI aortic plaque visualization is limited by the penetration depth of MRI surface coils and may lead to suboptimal image quality with insufficient reproducibility. By combining a transesophageal MRI (TEMRI) with surface MRI coils we enhanced local and overall image SNR for improved image quality and reproducibility.

GLUT8 is dispensable for embryonic development but influences hippocampal neurogenesis and heart function.

GLUT8 is a glucose transporter isoform expressed at high levels in testis; at intermediate levels in the brain, including the hippocampus; and at lower levels in the heart and several other tissues. GLUT8 is located in an intracellular compartment and does not appear to translocate to the cell surface, except in blastocysts, where insulin has been reported to induce its surface expression. Here, we generated mice with inactivation of the glut8 gene. We showed that expression of GLUT8 was not required for normal embryonic development and that glut8-/- mice had normal postnatal development, glucose homeostasis, and response to mild stress. Adult glut8-/- mice showed increased proliferation of hippocampal cells but no defect in memory acquisition and retention. Absence of GLUT8 from the heart did not alter heart size and morphology but led to an increase in P-wave duration, which was not associated with abnormal Nav1.5 Na+ channel or connexin expression. Thus, absence of GLUT8 expression in the mouse caused complex but mild physiological alterations.

Label-free cytotoxicity screening assay by digital holographic microscopy

Abstract We introduce a label-free technology based on digital holographic microscopy (DHM) with applicability for screening by imaging, and we demonstrate its capability for cytotoxicity assessment using mammalian living cells. For this first high content screening compatible application, we automatized a digital holographic microscope for image acquisition of cells using commercially available 96-well plates. Data generated through both label-free DHM imaging and fluorescence-based methods were in good agreement for cell viability identification and a Z'-factor close to 0.9 was determined, validating the robustness of DHM assay for phenotypic screening. Further, an excellent correlation was obtained between experimental cytotoxicity dose-response curves and known IC values for different toxic compounds. For comparable results, DHM has the major advantages of being label free and close to an order of magnitude faster than automated standard fluorescence microscopy.

Trinocular System for 3D motion and Dense Stucture Estimation

The relief of the seafloor is an important source of data for many scientists. In this paper we present an optical system to deal with underwater 3D reconstruction. This system is formed by three cameras that take images synchronously in a constant frame rate scheme. We use the images taken by these cameras to compute dense 3D reconstructions. We use Bundle Adjustment to estimate the motion ofthe trinocular rig. Given the path followed by the system, we get a dense map of the observed scene by registering the different dense local reconstructions in a unique and bigger one

Small scale underwater change detection

In this paper, we present a method to deal with the constraints of the underwater medium for finding changes between sequences of underwater images. One of the main problems of underwater medium for automatically detecting changes is the low altitude of the camera when taking pictures. This emphasise the parallax effect between the images as they are not taken exactly at the same position. In order to solve this problem, we are geometrically registering the images together taking into account the relief of the scene

Identification of geometrically consistent interest points for 3D scene reconstruction

We propose an algorithm that extracts image features that are consistent with the 3D structure of the scene. The features can be robustly tracked over multiple views and serve as vertices of planar patches that suitably represent scene surfaces, while reducing the redundancy in the description of 3D shapes. In other words, the extracted features will off er good tracking properties while providing the basis for 3D reconstruction with minimum model complexity

Copper acquisition by the SenC protein regulates aerobic respiration in Pseudomonas aeruginosa PAO1.

Aerobic respiration of Pseudomonas aeruginosa involves four terminal oxidases belonging to the heme-copper family (that is, three cytochrome c oxidases and one quinol oxidase) plus one copper-independent, cyanide-insensitive quinol oxidase (CIO). The PA0114 gene encoding an SCO1/SenC-type protein, which is known to be important for copper delivery to cytochrome c in yeast, Rhodobacter spp. and Agrobacterium tumefaciens, was found to be important for copper acquisition and aerobic respiration in P. aeruginosa. A PA0114 (senC) mutant grew poorly in low-copper media and had low cytochrome cbb(3)-type oxidase activity, but expressed CIO at increased levels, by comparison with the wild-type PAO1. Addition of copper reversed these phenotypes, suggesting that periplasmic copper capture by the SenC protein helps P. aeruginosa to adapt to copper deprivation.

Cognition and brain function in schizotypy : a selective review

Schizotypy refers to a set of personality traits thought to reflect the subclinical expression of the signs and symptoms of schizophrenia. Here, we review the cognitive and brain functional profile associated with high questionnaire scores in schizotypy. We discuss empirical evidence from the domains of perception, attention, memory, imagery and representation, language, and motor control. Perceptual deficits occur early and across various modalities. Whilst the neural mechanisms underlying visual impairments may be linked to magnocellular dysfunction, further effects may be seen downstream in higher cognitive functions. Cognitive deficits are observed in inhibitory control, selective and sustained attention, incidental learning and memory. In concordance with the cognitive nature of many of the aberrations of schizotypy, higher levels of schizotypy are associated with enhanced vividness and better performance on tasks of mental rotation. Language deficits seem most pronounced in higher-level processes. Finally, higher levels of schizotypy are associated with reduced performance on oculomotor tasks, resembling the impairments seen in schizophrenia. Some of these deficits are accompanied by reduced brain activation, akin to the pattern of hypoactivations in schizophrenia spectrum individuals. We conclude that schizotypy is a construct with apparent phenomenological overlap with schizophrenia and stable inter-individual differences that covary with performance on a wide range of perceptual, cognitive and motor tasks known to be impaired in schizophrenia. The importance of these findings lies not only in providing a fine-grained neurocognitive characterisation of a personality constellation known to be associated with real-life impairments, but also in generating hypotheses concerning the aetiology of schizophrenia.

High b-value diffusion-weighted imaging: A sensitive method to reveal white matter differences in schizophrenia.

Over the last 10 years, diffusion-weighted imaging (DWI) has become an important tool to investigate white matter (WM) anomalies in schizophrenia. Despite technological improvement and the exponential use of this technique, discrepancies remain and little is known about optimal parameters to apply for diffusion weighting during image acquisition. Specifically, high b-value diffusion-weighted imaging known to be more sensitive to slow diffusion is not widely used, even though subtle myelin alterations as thought to happen in schizophrenia are likely to affect slow-diffusing protons. Schizophrenia patients and healthy controls were scanned with a high b-value (4000s/mm(2)) protocol. Apparent diffusion coefficient (ADC) measures turned out to be very sensitive in detecting differences between schizophrenia patients and healthy volunteers even in a relatively small sample. We speculate that this is related to the sensitivity of high b-value imaging to the slow-diffusing compartment believed to reflect mainly the intra-axonal and myelin bound water pool. We also compared these results to a low b-value imaging experiment performed on the same population in the same scanning session. Even though the acquisition protocols are not strictly comparable, we noticed important differences in sensitivities in the favor of high b-value imaging, warranting further exploration.

Genetic exchange in an arbuscular mycorrhizal fungus results in increased rice growth and altered mycorrhiza-specific gene transcription.

Arbuscular mycorrhizal fungi (AMF) are obligate symbionts with most terrestrial plants. They improve plant nutrition, particularly phosphate acquisition, and thus are able to improve plant growth. In exchange, the fungi obtain photosynthetically fixed carbon. AMF are coenocytic, meaning that many nuclei coexist in a common cytoplasm. Genetic exchange recently has been demonstrated in the AMF Glomus intraradices, allowing nuclei of different Glomus intraradices strains to mix. Such genetic exchange was shown previously to have negative effects on plant growth and to alter fungal colonization. However, no attempt was made to detect whether genetic exchange in AMF can alter plant gene expression and if this effect was time dependent. Here, we show that genetic exchange in AMF also can be beneficial for rice growth, and that symbiosis-specific gene transcription is altered by genetic exchange. Moreover, our results show that genetic exchange can change the dynamics of the colonization of the fungus in the plant. Our results demonstrate that the simple manipulation of the genetics of AMF can have important consequences for their symbiotic effects on plants such as rice, which is considered the most important crop in the world. Exploiting natural AMF genetic variation by generating novel AMF genotypes through genetic exchange is a potentially useful tool in the development of AMF inocula that are more beneficial for crop growth.

Seismic reflectivity of the sediment-covered seafloor: Effects of velocity gradients and fine-scale layering

Knowledge of the reflectivity of the sediment-covered seabed is of significant importance to marine seismic data acquisition and interpretation as it governs the generation of reverberations in the water layer. In this context pertinent, but largely unresolved, questions concern the importance of the typically very prominent vertical seismic velocity gradients as well as the potential presence and magnitude of anisotropy in soft surficial seabed sediments. To address these issues, we explore the seismic properties of granulometric end-member-type clastic sedimentary seabed models consisting of sand, silt, and clay as well as scale-invariant stochastic layer sequences of these components characterized by realistic vertical gradients of the P- and S-wave velocities. Using effective media theory, we then assess the nature and magnitude of seismic anisotropy associated with these models. Our results indicate that anisotropy is rather benign for P-waves, and that the S-wave velocities in the axial directions differ only slightly. Because of the very high P- to S-wave velocity ratios in the vicinity of the seabed our models nevertheless suggest that S-wave triplications may occur at very small incidence angles. To numerically evaluate the P-wave reflection coefficient of our seabed models, we apply a frequency-slowness technique to the corresponding synthetic seismic wavefields. Comparison with analytical plane-wave reflection coefficients calculated for corresponding isotropic elastic half-space models shows that the differences tend to be most pronounced in the vicinity of the elastic equivalent of the critical angle as well as in the post-critical range. We also find that the presence of intrinsic anisotropy in the clay component of our layered models tends to dramatically reduce the overall magnitude of the P-wave reflection coefficient as well as its variation with incidence angle.