The potential of 3T high-resolution magnetic resonance imaging for diagnosis, staging, and follow-up of retinoblastoma.

We demonstrate the value of high-resolution magnetic resonance imaging (MRI) in diagnosing, staging, and follow-up of retinoblastoma during eye-saving treatment. We have included informative retinoblastoma cases scanned on a 3T MRI system from a retrospective retinoblastoma cohort from 2009 through 2013. We show that high-resolution MRI has the potential to detect small intraocular seeds, hemorrhage, and metastatic risk factors not visible with fundoscopy (e.g., optic nerve invasion and choroidal invasion), and treatment response. Unfortunately, however, the diagnostic accuracy of high-resolution MRI is not perfect, especially for subtle intraocular seeds or minimal postlaminar optic nerve invasion. The most important application of MRI is the detection of metastatic risk factors, as these cannot be found by fundoscopy and ultrasound.

Cardiovascular magnetic resonance imaging techniques for motion-suppressed multi-dimensional coronary artery imaging with high spatial and temporal resolution

L'imagerie par résonance magnétique (IRM) peut fournir aux cardiologues des informations diagnostiques importantes sur l'état de la maladie de l'artère coronarienne dans les patients. Le défi majeur pour l'IRM cardiaque est de gérer toutes les sources de mouvement qui peuvent affecter la qualité des images en réduisant l'information diagnostique. Cette thèse a donc comme but de développer des nouvelles techniques d'acquisitions des images IRM, en changeant les techniques de compensation du mouvement, pour en augmenter l'efficacité, la flexibilité, la robustesse et pour obtenir plus d'information sur le tissu et plus d'information temporelle. Les techniques proposées favorisent donc l'avancement de l'imagerie des coronaires dans une direction plus maniable et multi-usage qui peut facilement être transférée dans l'environnement clinique. La première partie de la thèse s'est concentrée sur l'étude du mouvement des artères coronariennes sur des patients en utilisant la techniques d'imagerie standard (rayons x), pour mesurer la précision avec laquelle les artères coronariennes retournent dans la même position battement après battement (repositionnement des coronaires). Nous avons découvert qu'il y a des intervalles dans le cycle cardiaque, tôt dans la systole et à moitié de la diastole, où le repositionnement des coronaires est au minimum. En réponse nous avons développé une nouvelle séquence d'acquisition (T2-post) capable d'acquérir les données aussi tôt dans la systole. Cette séquence a été testée sur des volontaires sains et on a pu constater que la qualité de visualisation des artère coronariennes est égale à celle obtenue avec les techniques standard. De plus, le rapport signal sur bruit fourni par la séquence d'acquisition proposée est supérieur à celui obtenu avec les techniques d'imagerie standard. La deuxième partie de la thèse a exploré un paradigme d'acquisition des images cardiaques complètement nouveau pour l'imagerie du coeur entier. La technique proposée dans ce travail acquiert les données sans arrêt (free-running) au lieu d'être synchronisée avec le mouvement cardiaque. De cette façon, l'efficacité de la séquence d'acquisition est augmentée de manière significative et les images produites représentent le coeur entier dans toutes les phases cardiaques (quatre dimensions, 4D). Par ailleurs, l'auto-navigation de la respiration permet d'effectuer cette acquisition en respiration libre. Cette technologie rend possible de visualiser et évaluer l'anatomie du coeur et de ses vaisseaux ainsi que la fonction cardiaque en quatre dimensions et avec une très haute résolution spatiale et temporelle, sans la nécessité d'injecter un moyen de contraste. Le pas essentiel qui a permis le développement de cette technique est l'utilisation d'une trajectoire d'acquisition radiale 3D basée sur l'angle d'or. Avec cette trajectoire, il est possible d'acquérir continûment les données d'espace k, puis de réordonner les données et choisir les paramètres temporel des images 4D a posteriori. L'acquisition 4D a été aussi couplée avec un algorithme de reconstructions itératif (compressed sensing) qui permet d'augmenter la résolution temporelle tout en augmentant la qualité des images. Grâce aux images 4D, il est possible maintenant de visualiser les artères coronariennes entières dans chaque phase du cycle cardiaque et, avec les mêmes données, de visualiser et mesurer la fonction cardiaque. La qualité des artères coronariennes dans les images 4D est la même que dans les images obtenues avec une acquisition 3D standard, acquise en diastole Par ailleurs, les valeurs de fonction cardiaque mesurées au moyen des images 4D concorde avec les valeurs obtenues avec les images 2D standard. Finalement, dans la dernière partie de la thèse une technique d'acquisition a temps d'écho ultra-court (UTE) a été développée pour la visualisation in vivo des calcifications des artères coronariennes. Des études récentes ont démontré que les acquisitions UTE permettent de visualiser les calcifications dans des plaques athérosclérotiques ex vivo. Cepandent le mouvement du coeur a entravé jusqu'à maintenant l'utilisation des techniques UTE in vivo. Pour résoudre ce problème nous avons développé une séquence d'acquisition UTE avec trajectoire radiale 3D et l'avons testée sur des volontaires. La technique proposée utilise une auto-navigation 3D pour corriger le mouvement respiratoire et est synchronisée avec l'ECG. Trois échos sont acquis pour extraire le signal de la calcification avec des composants au T2 très court tout en permettant de séparer le signal de la graisse depuis le signal de l'eau. Les résultats sont encore préliminaires mais on peut affirmer que la technique développé peut potentiellement montrer les calcifications des artères coronariennes in vivo. En conclusion, ce travail de thèse présente trois nouvelles techniques pour l'IRM du coeur entier capables d'améliorer la visualisation et la caractérisation de la maladie athérosclérotique des coronaires. Ces techniques fournissent des informations anatomiques et fonctionnelles en quatre dimensions et des informations sur la composition du tissu auparavant indisponibles. CORONARY artery magnetic resonance imaging (MRI) has the potential to provide the cardiologist with relevant diagnostic information relative to coronary artery disease of patients. The major challenge of cardiac MRI, though, is dealing with all sources of motions that can corrupt the images affecting the diagnostic information provided. The current thesis, thus, focused on the development of new MRI techniques that change the standard approach to cardiac motion compensation in order to increase the efficiency of cardioavscular MRI, to provide more flexibility and robustness, new temporal information and new tissue information. The proposed approaches help in advancing coronary magnetic resonance angiography (MRA) in the direction of an easy-to-use and multipurpose tool that can be translated to the clinical environment. The first part of the thesis focused on the study of coronary artery motion through gold standard imaging techniques (x-ray angiography) in patients, in order to measure the precision with which the coronary arteries assume the same position beat after beat (coronary artery repositioning). We learned that intervals with minimal coronary artery repositioning occur in peak systole and in mid diastole and we responded with a new pulse sequence (T2~post) that is able to provide peak-systolic imaging. Such a sequence was tested in healthy volunteers and, from the image quality comparison, we learned that the proposed approach provides coronary artery visualization and contrast-to-noise ratio (CNR) comparable with the standard acquisition approach, but with increased signal-to-noise ratio (SNR). The second part of the thesis explored a completely new paradigm for whole- heart cardiovascular MRI. The proposed techniques acquires the data continuously (free-running), instead of being triggered, thus increasing the efficiency of the acquisition and providing four dimensional images of the whole heart, while respiratory self navigation allows for the scan to be performed in free breathing. This enabling technology allows for anatomical and functional evaluation in four dimensions, with high spatial and temporal resolution and without the need for contrast agent injection. The enabling step is the use of a golden-angle based 3D radial trajectory, which allows for a continuous sampling of the k-space and a retrospective selection of the timing parameters of the reconstructed dataset. The free-running 4D acquisition was then combined with a compressed sensing reconstruction algorithm that further increases the temporal resolution of the 4D dataset, while at the same time increasing the overall image quality by removing undersampling artifacts. The obtained 4D images provide visualization of the whole coronary artery tree in each phases of the cardiac cycle and, at the same time, allow for the assessment of the cardiac function with a single free- breathing scan. The quality of the coronary arteries provided by the frames of the free-running 4D acquisition is in line with the one obtained with the standard ECG-triggered one, and the cardiac function evaluation matched the one measured with gold-standard stack of 2D cine approaches. Finally, the last part of the thesis focused on the development of ultrashort echo time (UTE) acquisition scheme for in vivo detection of calcification in the coronary arteries. Recent studies showed that UTE imaging allows for the coronary artery plaque calcification ex vivo, since it is able to detect the short T2 components of the calcification. The heart motion, though, prevented this technique from being applied in vivo. An ECG-triggered self-navigated 3D radial triple- echo UTE acquisition has then been developed and tested in healthy volunteers. The proposed sequence combines a 3D self-navigation approach with a 3D radial UTE acquisition enabling data collection during free breathing. Three echoes are simultaneously acquired to extract the short T2 components of the calcification while a water and fat separation technique allows for proper visualization of the coronary arteries. Even though the results are still preliminary, the proposed sequence showed great potential for the in vivo visualization of coronary artery calcification. In conclusion, the thesis presents three novel MRI approaches aimed at improved characterization and assessment of atherosclerotic coronary artery disease. These approaches provide new anatomical and functional information in four dimensions, and support tissue characterization for coronary artery plaques.

In vivo imaging of the central and peripheral effects of sleep deprivation and suprachiasmatic nuclei lesion on PERIOD-2 protein in mice.

STUDY OBJECTIVES: That sleep deprivation increases the brain expression of various clock genes has been well documented. Based on these and other findings we hypothesized that clock genes not only underlie circadian rhythm generation but are also implicated in sleep homeostasis. However, long time lags have been reported between the changes in the clock gene messenger RNA levels and their encoded proteins. It is therefore crucial to establish whether also protein levels increase within the time frame known to activate a homeostatic sleep response. We report on the central and peripheral effects of sleep deprivation on PERIOD-2 (PER2) protein both in intact and suprachiasmatic nuclei-lesioned mice. DESIGN: In vivo and in situ PER2 imaging during baseline, sleep deprivation, and recovery. SETTINGS: Mouse sleep-recording facility. PARTICIPANTS: Per2::Luciferase knock-in mice. INTERVENTIONS: N/A. MEASUREMENTS AND RESULTS: Six-hour sleep deprivation increased PER2 not only in the brain but also in liver and kidney. Remarkably, the effects in the liver outlasted those observed in the brain. Within the brain the increase in PER2 concerned the cerebral cortex mainly, while leaving suprachiasmatic nuclei (SCN) levels unaffected. Against expectation, sleep deprivation did not increase PER2 in the brain of arrhythmic SCN-lesioned mice because of higher PER2 levels in baseline. In contrast, liver PER2 levels did increase in these mice similar to the sham and partially lesioned controls. CONCLUSIONS: Our results stress the importance of considering both sleep-wake dependent and circadian processes when quantifying clock-gene levels. Because sleep deprivation alters PERIOD-2 in the brain as well as in the periphery, it is tempting to speculate that clock genes constitute a common pathway mediating the shared and well-known adverse effects of both chronic sleep loss and disrupted circadian rhythmicity on metabolic health.

(MA)LDI MS Imaging at High Specificity and Sensitivity

(Matrix-assisted) laser desorption/ionization ((MA)LDI) mass spectrometry imaging (MSI) has been driven by remarkable technological developments in the last couple of years. Although molecular information of a wide range of molecules including peptides, lipids, metabolites, and xenobiotics can be mapped, (MA)LDI MSI only leads to the detection of the most abundant soluble molecules in the cells and, consequently, does not provide access to the least expressed species, which can be very informative in the scope of disease research. Within a short period of time, numerous protocols and concepts have been developed and introduced in order to increase MSI sensitivity, including in situ tissue chemistry and solvent-free matrix depositions. In this chapter, we will discuss some of the latest developments in the field of high-sensitivity MSI using solvent-free matrix depositions and will detail protocols of two methods with their capability of enriching molecular MSI signal as demonstrated within our laboratory.

Microtubule-associated protein-2 immunoreactivity: a useful tool in the differential diagnosis of low-grade neuroepithelial tumors.

Complex and variable morphological phenotypes pose a major challenge to the histopathological classification of neuroepithelial tumors. This applies in particular for low-grade gliomas and glio-neuronal tumors. Recently, we and others have identified microtubule-associated protein-2 (MAP2) as an immunohistochemical marker expressed in the majority of glial tumors. Characteristic cell morphologies can be recognized by MAP2 immunoreactivity in different glioma entities, i.e., process sparse oligodendroglial versus densely ramified astrocytic elements. Here, we describe MAP2-immunoreactivity patterns in a large series of various neuroepithelial tumors and related neoplasms (n = 960). Immunohistochemical analysis led to the following conclusions: (1) specific pattern of MAP2-positive tumor cells can be identified in 95% of glial neoplasms; (2) ependymal tumors do not express MAP2 in their rosette-forming cell component; (3) tumors of the pineal gland as well as malignant embryonic tumors are also characterized by abundant MAP2 immunoreactivity; (4) virtually no MAP2 expression can be observed in the neoplastic glial component of glio-neuronal tumors, i.e. gangliogliomas; (5) malignant glial tumor variants (WHO grade III or IV) exhibit different and less specific MAP2 staining patterns compared to their benign counterparts (WHO grade I or II); (6) with the exception of melanomas and small cell lung cancers, MAP2 expression is very rare in metastatic and non-neuroepithelial tumors; (7) glial MAP2 expression was not detected in 56 non-neoplastic lesions. These data point towards MAP2 as valuable diagnostic tool for pattern recognition and differential diagnosis of low-grade neuroepithelial tumors.

An active contour-based atlas registration model applied to automatic subthalamic nucleus targeting on MRI: method and validation.

This paper presents a new non parametric atlas registration framework, derived from the optical flow model and the active contour theory, applied to automatic subthalamic nucleus (STN) targeting in deep brain stimulation (DBS) surgery. In a previous work, we demonstrated that the STN position can be predicted based on the position of surrounding visible structures, namely the lateral and third ventricles. A STN targeting process can thus be obtained by registering these structures of interest between a brain atlas and the patient image. Here we aim to improve the results of the state of the art targeting methods and at the same time to reduce the computational time. Our simultaneous segmentation and registration model shows mean STN localization errors statistically similar to the most performing registration algorithms tested so far and to the targeting expert's variability. Moreover, the computational time of our registration method is much lower, which is a worthwhile improvement from a clinical point of view.

Mental imagery for full and upper human bodies: common right hemisphere activations and distinct extrastriate activations.

The processing of human bodies is important in social life and for the recognition of another person's actions, moods, and intentions. Recent neuroimaging studies on mental imagery of human body parts suggest that the left hemisphere is dominant in body processing. However, studies on mental imagery of full human bodies reported stronger right hemisphere or bilateral activations. Here, we measured functional magnetic resonance imaging during mental imagery of bilateral partial (upper) and full bodies. Results show that, independently of whether a full or upper body is processed, the right hemisphere (temporo-parietal cortex, anterior parietal cortex, premotor cortex, bilateral superior parietal cortex) is mainly involved in mental imagery of full or partial human bodies. However, distinct activations were found in extrastriate cortex for partial bodies (right fusiform face area) and full bodies (left extrastriate body area). We propose that a common brain network, mainly on the right side, is involved in the mental imagery of human bodies, while two distinct brain areas in extrastriate cortex code for mental imagery of full and upper bodies.

Pattern dystrophy with high intrafamilial variability associated with Y141C mutation in the peripherin/RDS gene and successful treatment of subfoveal CNV related to multifocal pattern type with anti-VEGF (ranibizumab) intravitreal injections.

OBJECTIVE: To identify disease causing mutation in three generations of a Swiss family with pattern dystrophy and high intrafamilial variability of phenotype. To assess the effect of intravitreal ranibizumab injections in the treatment of subfoveal choroidal neovascularization associated with pattern dystrophy in one patient. METHODS: Affected family members were ascertained for phenotypic and genotypic characterization. Ophthalmic evaluations included fundus photography, autofluorescence imaging, optical coherence tomography, and International Society for Clinical Electrophysiology of Vision standard full-field electroretinography. When possible family members had genetic testing. The proband presented with choroidal neovascularization and had intravitreal injections as needed according to visual acuity and optical coherence tomography. RESULTS: Proband had a multifocal type pattern dystrophy, and his choroidal neovascularization regressed after four intravitreal injections. The vision improved from 0.8 to 1.0, and optical coherence tomography showed complete anatomical restoration. A butterfly-shaped pattern was observed in her cousin, whereas a fundus pulverulentus pattern was seen in a second cousin. Aunt had a multifocal atrophic appearance, simulating geographic atrophy in age-related macular degeneration. The Y141C mutation was identified in the peripherin/RDS gene and segregated with disease in the family. CONCLUSION: This is the first report of marked intrafamilial variation of pattern dystrophy because of peripherin/RDS Y141C mutation. Intravitreal ranibizumab injections might be a valuable treatment for associated subfoveal choroidal neovascularization.

Comparison and validation of tissue modelization and statistical classification methods in T1-weighted MR brain images.

This paper presents a validation study on statistical nonsupervised brain tissue classification techniques in magnetic resonance (MR) images. Several image models assuming different hypotheses regarding the intensity distribution model, the spatial model and the number of classes are assessed. The methods are tested on simulated data for which the classification ground truth is known. Different noise and intensity nonuniformities are added to simulate real imaging conditions. No enhancement of the image quality is considered either before or during the classification process. This way, the accuracy of the methods and their robustness against image artifacts are tested. Classification is also performed on real data where a quantitative validation compares the methods' results with an estimated ground truth from manual segmentations by experts. Validity of the various classification methods in the labeling of the image as well as in the tissue volume is estimated with different local and global measures. Results demonstrate that methods relying on both intensity and spatial information are more robust to noise and field inhomogeneities. We also demonstrate that partial volume is not perfectly modeled, even though methods that account for mixture classes outperform methods that only consider pure Gaussian classes. Finally, we show that simulated data results can also be extended to real data.

Representing diffusion MRI in 5D for segmentation of white matter tracts with a level set method.

We present a method for segmenting white matter tracts from high angular resolution diffusion MR. images by representing the data in a 5 dimensional space of position and orientation. Whereas crossing fiber tracts cannot be separated in 3D position space, they clearly disentangle in 5D position-orientation space. The segmentation is done using a 5D level set method applied to hyper-surfaces evolving in 5D position-orientation space. In this paper we present a methodology for constructing the position-orientation space. We then show how to implement the standard level set method in such a non-Euclidean high dimensional space. The level set theory is basically defined for N-dimensions but there are several practical implementation details to consider, such as mean curvature. Finally, we will show results from a synthetic model and a few preliminary results on real data of a human brain acquired by high angular resolution diffusion MRI.

Early processing in the human lateral occipital complex is highly responsive to illusory contours but not to salient regions.

Human electrophysiological studies support a model whereby sensitivity to so-called illusory contour stimuli is first seen within the lateral occipital complex. A challenge to this model posits that the lateral occipital complex is a general site for crude region-based segmentation, based on findings of equivalent hemodynamic activations in the lateral occipital complex to illusory contour and so-called salient region stimuli, a stimulus class that lacks the classic bounding contours of illusory contours. Using high-density electrical mapping of visual evoked potentials, we show that early lateral occipital cortex activity is substantially stronger to illusory contour than to salient region stimuli, whereas later lateral occipital complex activity is stronger to salient region than to illusory contour stimuli. Our results suggest that equivalent hemodynamic activity to illusory contour and salient region stimuli probably reflects temporally integrated responses, a result of the poor temporal resolution of hemodynamic imaging. The temporal precision of visual evoked potentials is critical for establishing viable models of completion processes and visual scene analysis. We propose that crude spatial segmentation analyses, which are insensitive to illusory contours, occur first within dorsal visual regions, not the lateral occipital complex, and that initial illusory contour sensitivity is a function of the lateral occipital complex.

New perspectives in the use of ink evidence in forensic science: Part II. Development and testing of mathematical algorithms for the automatic comparison of ink samples analysed by HPTLC

In the first part of this research, three stages were stated for a program to increase the information extracted from ink evidence and maximise its usefulness to the criminal and civil justice system. These stages are (a) develop a standard methodology for analysing ink samples by high-performance thin layer chromatography (HPTLC) in reproducible way, when ink samples are analysed at different time, locations and by different examiners; (b) compare automatically and objectively ink samples; and (c) define and evaluate theoretical framework for the use of ink evidence in forensic context. This report focuses on the second of the three stages. Using the calibration and acquisition process described in the previous report, mathematical algorithms are proposed to automatically and objectively compare ink samples. The performances of these algorithms are systematically studied for various chemical and forensic conditions using standard performance tests commonly used in biometrics studies. The results show that different algorithms are best suited for different tasks. Finally, this report demonstrates how modern analytical and computer technology can be used in the field of ink examination and how tools developed and successfully applied in other fields of forensic science can help maximising its impact within the field of questioned documents.

Atlas-based segmentation of pathological MR brain images using a model of lesion growth.

We propose a method for brain atlas deformation in the presence of large space-occupying tumors, based on an a priori model of lesion growth that assumes radial expansion of the lesion from its starting point. Our approach involves three steps. First, an affine registration brings the atlas and the patient into global correspondence. Then, the seeding of a synthetic tumor into the brain atlas provides a template for the lesion. The last step is the deformation of the seeded atlas, combining a method derived from optical flow principles and a model of lesion growth. Results show that a good registration is performed and that the method can be applied to automatic segmentation of structures and substructures in brains with gross deformation, with important medical applications in neurosurgery, radiosurgery, and radiotherapy.

A surface-based approach to quantify local cortical gyrification.

The high complexity of cortical convolutions in humans is very challenging both for engineers to measure and compare it, and for biologists and physicians to understand it. In this paper, we propose a surface-based method for the quantification of cortical gyrification. Our method uses accurate 3-D cortical reconstruction and computes local measurements of gyrification at thousands of points over the whole cortical surface. The potential of our method to identify and localize precisely gyral abnormalities is illustrated by a clinical study on a group of children affected by 22q11 Deletion Syndrome, compared to control individuals.

ZHARP: three-dimensional motion tracking from a single image plane.

Three-dimensional imaging and quantification of myocardial function are essential steps in the evaluation of cardiac disease. We propose a tagged magnetic resonance imaging methodology called zHARP that encodes and automatically tracks myocardial displacement in three dimensions. Unlike other motion encoding techniques, zHARP encodes both in-plane and through-plane motion in a single image plane without affecting the acquisition speed. Postprocessing unravels this encoding in order to directly track the 3-D displacement of every point within the image plane throughout an entire image sequence. Experimental results include a phantom validation experiment, which compares zHARP to phase contrast imaging, and an in vivo study of a normal human volunteer. Results demonstrate that the simultaneous extraction of in-plane and through-plane displacements from tagged images is feasible.