Model-based Segmentation and Fusion of 3D Computed Tomography and 3D Ultrasound of the Eye for Radiotherapy Planning

Computed Tomography (CT) represents the standard imaging modality for tumor volume delineation for radiotherapy treatment planning of retinoblastoma despite some inherent limitations. CT scan is very useful in providing information on physical density for dose calculation and morphological volumetric information but presents a low sensitivity in assessing the tumor viability. On the other hand, 3D ultrasound (US) allows a highly accurate definition of the tumor volume thanks to its high spatial resolution but it is not currently integrated in the treatment planning but used only for diagnosis and follow-up. Our ultimate goal is an automatic segmentation of gross tumor volume (GTV) in the 3D US, the segmentation of the organs at risk (OAR) in the CT and the registration of both modalities. In this paper, we present some preliminary results in this direction. We present 3D active contour-based segmentation of the eye ball and the lens in CT images; the presented approach incorporates the prior knowledge of the anatomy by using a 3D geometrical eye model. The automated segmentation results are validated by comparing with manual segmentations. Then, we present two approaches for the fusion of 3D CT and US images: (i) landmark-based transformation, and (ii) object-based transformation that makes use of eye ball contour information on CT and US images.

Visuo-motor pathways in humans revealed by event-related fMRI.

Whether different brain networks are involved in generating unimanual responses to a simple visual stimulus presented in the ipsilateral versus contralateral hemifield remains a controversial issue. Visuo-motor routing was investigated with event-related functional magnetic resonance imaging (fMRI) using the Poffenberger reaction time task. A 2 hemifield x 2 response hand design generated the "crossed" and "uncrossed" conditions, describing the spatial relation between these factors. Both conditions, with responses executed by the left or right hand, showed a similar spatial pattern of activated areas, including striate and extrastriate areas bilaterally, SMA, and M1 contralateral to the responding hand. These results demonstrated that visual information is processed bilaterally in striate and extrastriate visual areas, even in the "uncrossed" condition. Additional analyses based on sorting data according to subjects' reaction times revealed differential crossed versus uncrossed activity only for the slowest trials, with response strength in infero-temporal cortices significantly correlating with crossed-uncrossed differences (CUD) in reaction times. Collectively, the data favor a parallel, distributed model of brain activation. The presence of interhemispheric interactions and its consequent bilateral activity is not determined by the crossed anatomic projections of the primary visual and motor pathways. Distinct visuo-motor networks need not be engaged to mediate behavioral responses for the crossed visual field/response hand condition. While anatomical connectivity heavily influences the spatial pattern of activated visuo-motor pathways, behavioral and functional parameters appear to also affect the strength and dynamics of responses within these pathways.

Positive contrast MR-lymphography using inversion recovery with ON-resonant water suppression (IRON).

PURPOSE: To investigate the utility of inversion recovery with ON-resonant water suppression (IRON) to create positive signal in normal lymph nodes after injection of superparamagnetic nanoparticles. MATERIALS AND METHODS: Experiments were conducted on six rabbits, which received a single bolus injection of 80 mumol Fe/kg monocrystalline iron oxide nanoparticle (MION-47). Magnetic resonance imaging (MRI) was performed at baseline, 1 day, and 3 days after MION-47 injection using conventional T(1)- and T(2)*-weighted sequences and IRON. Contrast-to-noise ratios (CNR) were measured in blood and in paraaortic lymph nodes. RESULTS: On T(2)*-weighted images, as expected, signal attenuation was observed in areas of paraaortic lymph nodes after MION-47 injection. However, using IRON the paraaortic lymph nodes exhibited very high contrast enhancement, which remained 3 days after injection. CNR with IRON was 2.2 +/- 0.8 at baseline, increased markedly 1 day after injection (23.5 +/- 5.4, P < 0.01 vs. baseline), and remained high after 3 days (21.8 +/- 5.7, *P < 0.01 vs. baseline). CNR was also high in blood 1 day after injection (42.7 +/- 7.2 vs. 1.8 +/- 0.7 at baseline, P < 0.01) but approached baseline after 3 days (1.9 +/- 1.4, P = NS vs. baseline). CONCLUSION: IRON in conjunction with superparamagnetic nanoparticles can be used to perform 'positive contrast' MR-lymphography, particularly 3 days after injection of the contrast agent, when signal is no longer visible within blood vessels. The proposed method may have potential as an adjunct for nodal staging in cancer screening.

Relationship between motion of coronary arteries and diaphragm during free breathing: lessons from real-time MR imaging.

OBJECTIVE: Diaphragmatic navigators are frequently used in free-breathing coronary MR angiography, either to gate or prospectively correct slice position or both. For such approaches, a constant relationship between coronary and diaphragmatic displacement throughout the respiratory cycle is assumed. The purpose of this study was to evaluate the relationship between diaphragmatic and coronary artery motion during free breathing. SUBJECTS AND METHODS: A real-time echoplanar MR imaging sequence was used in 12 healthy volunteers to obtain 30 successive images each (one per cardiac cycle) that included the left main coronary artery and the domes of both hemidiaphragms. The coronary artery and diaphragm positions (relative to isocenter) were determined and analyzed for effective diaphragmatic gating windows of 3, 5, and 7 mm (diaphragmatic excursions of 0-3, 0-5, and 0-7 mm from the end-expiratory position, respectively). RESULTS: Although the mean slope correlating the displacement of the right diaphragm and the left main coronary artery was approximately 0.6 for all diaphragmatic gating windows, we also found great variability among individual volunteers. Linear regression slopes varied from 0.17 to 0.93, and r2 values varied from .04 to .87. CONCLUSION: Wide individual variability exists in the relationship between coronary and diaphragmatic respiratory motion during free breathing. Accordingly, coronary MR angiographic approaches that use diaphragmatic navigator position for prospective slice correction may benefit from patient-specific correction factors. Alternatively, coronary MR angiography may benefit from a more direct assessment of the respiratory displacement of the heart and coronary arteries, using left ventricular navigators.

Three-dimensional ordered-subset expectation maximization iterative protocol for evaluation of left ventricular volumes and function by quantitative gated SPECT: a dynamic phantom study.

The purposes of this study were to characterize the performance of a 3-dimensional (3D) ordered-subset expectation maximization (OSEM) algorithm in the quantification of left ventricular (LV) function with (99m)Tc-labeled agent gated SPECT (G-SPECT), the QGS program, and a beating-heart phantom and to optimize the reconstruction parameters for clinical applications. METHODS: A G-SPECT image of a dynamic heart phantom simulating the beating left ventricle was acquired. The exact volumes of the phantom were known and were as follows: end-diastolic volume (EDV) of 112 mL, end-systolic volume (ESV) of 37 mL, and stroke volume (SV) of 75 mL; these volumes produced an LV ejection fraction (LVEF) of 67%. Tomographic reconstructions were obtained after 10-20 iterations (I) with 4, 8, and 16 subsets (S) at full width at half maximum (FWHM) gaussian postprocessing filter cutoff values of 8-15 mm. The QGS program was used for quantitative measurements. RESULTS: Measured values ranged from 72 to 92 mL for EDV, from 18 to 32 mL for ESV, and from 54 to 63 mL for SV, and the calculated LVEF ranged from 65% to 76%. Overall, the combination of 10 I, 8 S, and a cutoff filter value of 10 mm produced the most accurate results. The plot of the measures with respect to the expectation maximization-equivalent iterations (I x S product) revealed a bell-shaped curve for the LV volumes and a reverse distribution for the LVEF, with the best results in the intermediate range. In particular, FWHM cutoff values exceeding 10 mm affected the estimation of the LV volumes. CONCLUSION: The QGS program is able to correctly calculate the LVEF when used in association with an optimized 3D OSEM algorithm (8 S, 10 I, and FWHM of 10 mm) but underestimates the LV volumes. However, various combinations of technical parameters, including a limited range of I and S (80-160 expectation maximization-equivalent iterations) and low cutoff values (< or =10 mm) for the gaussian postprocessing filter, produced results with similar accuracies and without clinically relevant differences in the LV volumes and the estimated LVEF.

Deviant trajectories of cortical maturation in 22q11.2 deletion syndrome (22q11DS): a cross-sectional and longitudinal study.

22q11.2 deletion syndrome (22q11DS) is associated with an increased susceptibility to develop schizophrenia. Despite a large body of literature documenting abnormal brain structure in 22q11DS, cerebral changes associated with brain maturation in 22q11DS remained largely unexplored. To map cortical maturation from childhood to adulthood in 22q11.2 deletion syndrome, we used cerebral MRI from 59 patients with 22q11DS, aged 6 to 40, and 80 typically developing controls; three year follow-up assessments were also available for 32 patients and 31 matched controls. Cross-sectional cortical thickness trajectories during childhood and adolescence were approximated in age bins. Repeated-measures were also conducted with the longitudinal data. Within the group of patients with 22q11DS, exploratory measures of cortical thickness differences related to COMT polymorphism, IQ, and schizophrenia were also conducted. We observed deviant trajectories of cortical thickness changes with age in patients with 22q11DS. In affected preadolescents, larger prefrontal thickness was observed compared to age-matched controls. Afterward, we observed greater cortical loss in 22q11DS with a convergence of cortical thickness values by the end of adolescence. No compelling evidence for an effect of COMT polymorphism on cortical maturation was observed. Within 22q11DS, significant differences in cortical thickness were related to cognitive level in children and adolescents, and to schizophrenia in adults. Deviant trajectories of cortical thickness from childhood to adulthood provide strong in vivo cues for a defect in the programmed synaptic elimination, which in turn may explain the susceptibility of patients with 22q11DS to develop psychosis.

Correction for heart rate variability during 3D whole heart MR coronary angiography.

PURPOSE: To evaluate the effect of a real-time adaptive trigger delay on image quality to correct for heart rate variability in 3D whole-heart coronary MR angiography (MRA). MATERIALS AND METHODS: Twelve healthy adults underwent 3D whole-heart coronary MRA with and without the use of an adaptive trigger delay. The moment of minimal coronary artery motion was visually determined on a high temporal resolution MRI. Throughout the scan performed without adaptive trigger delay, trigger delay was kept constant, whereas during the scan performed with adaptive trigger delay, trigger delay was continuously updated after each RR-interval using physiological modeling. Signal-to-noise, contrast-to-noise, vessel length, vessel sharpness, and subjective image quality were compared in a blinded manner. RESULTS: Vessel sharpness improved significantly for the middle segment of the right coronary artery (RCA) with the use of the adaptive trigger delay (52.3 +/- 7.1% versus 48.9 +/- 7.9%, P = 0.026). Subjective image quality was significantly better in the middle segments of the RCA and left anterior descending artery (LAD) when the scan was performed with adaptive trigger delay compared to constant trigger delay. CONCLUSION: Our results demonstrate that the use of an adaptive trigger delay to correct for heart rate variability improves image quality mainly in the middle segments of the RCA and LAD.

Coronary artery imaging by magnetic resonance.

Non-invasive visualization of the coronary arteries represents a major challenge in modern cardiology, but this goal may be achieved in the near future by MR angiography. Possible applications are non-invasive diagnosis of coronary artery disease, and follow-up examinations for therapy control after PTCA, in order to detect restenosis at an early stage. A multiple slice technique (2 mm slice thickness, with a spatial resolution of 1 x 1 mm, Philips Gyroscan ACS-II, 1.5 Tesla) was used. Ten volunteers were imaged and 10 patients with coronary artery disease were examined before and after PTCA. MR measurements were validated by quantitative coronary angiography. The diameters of the proximal coronary arteries as measured by both methods were compared, and a good correlation was found (r = 0.76). Thus, it is concluded that non-invasive visualization of the coronary arteries is possible before and after PTCA and allows to determine potential restenoses. However, patient cooperation is essential for good image quality. Moreover, limited spatial image resolution and breathing artifacts restrict MR coronary angiography today to be used as a routine diagnostic tool for the diagnosis of coronary artery disease.

Submillimeter three-dimensional coronary MR angiography with real-time navigator correction: comparison of navigator locations.

Three-dimensional free-breathing coronary magnetic resonance angiography was performed in eight healthy volunteers with use of real-time navigator technology. Images acquired with the navigator localized at the right hemidiaphragm and at the left ventricle were objectively compared. The diaphragmatic navigator was found to be superior for vessel delineation of middle to distal portions of the coronary arteries.

Quantification of the local heartwall motion by magnetic resonance myocardial tagging.

Sophisticated magnetic resonance tagging techniques provide powerful tools for the non-invasive assessment of the local heartwall motion towards a deeper fundamental understanding of local heart function. For the extraction of motion data from the time series of magnetic resonance tagged images and for the visualization of the local heartwall motion a new image analysis procedure has been developed. New parameters have been derived which allows quantification of the motion patterns and are highly sensitive to any changes in these patterns. The new procedure has been applied for heart motion analysis in healthy volunteers and in patient collectives with different heart diseases. The achieved results are summarized and discussed.

Molecular genetics of ocular diseases

The eye is a complex organ, which provides one of our most important senses, sight. The retina is the neuronal component of the eye and represents the connection with the central nervous system for the transmission of the information that leads to image processing. Retinitis pigmentosa (RP) is one of the most common forms of inherited retinal degeneration, in which the primary death of rods, resulting in night blindness, is always followed by the loss of cones, which leads to legal blindness. Clinical and genetic heterogeneity in retinitis pigmentosa is not only due to different mutations in different genes, but also to different effects of the same mutation in different individuals, sometimes even within the same family. My thesis work has been mainly focused on an autosomal dominant form of RP linked to mutations in the PRPF31 gene, which often shows reduced penetrance. Our study has led to the identification of the major regulator of the penetrance of PRPF31 mutations, the CNOT3 protein, and to the characterization of its mechanism of action. Following the same rationale of investigating molecular mechanisms that are responsible for clinical and genetic heterogeneity of retinitis pigmentosa, we studied a recessive form of the disease associated with mutations in the recently-identified gene FAMI61 A, where mutations in the same gene give rise to variable clinical manifestations. Our data have increased the knowledge of the relationship between genotype and phenotype in this form of the disease. Whole genome sequencing technique was also tested as a strategy for disease gene identification in unrelated patients with recessive retinitis pigmentosa and proved to be effective in identifying disease-causing variants that might have otherwise failed to be detected with other screening methods. Finally, for the first time we reported a choroidal tumor among the clinical manifestations of PTEN hamartoma tumor syndrome, a genetic disorder caused by germline mutations of the tumor suppressor gene PTEN. Our study has highlighted the heterogeneity of this choroidal tumor, showing that genetic and/or epigenetic alterations in different genes may contribute to the tumor development and growth. - L'oeil est un organe complexe, à l'origine d'un de nos sens les plus importants, la vue. La rétine est la composante neuronale de l'oeil qui constitue la connexion avec le système nerveux central pour la transmission de l'information et qui conduit à la formation des images. La rétinite pigmentaire (RP) est une des formes les plus courantes de dégénérescence rétinienne héréditaire, dans laquelle la mort primaire de bâtonnets, entraînant la cécité nocturne, est toujours suivie par la perte de cônes qui conduit à la cécité complète. L'hétérogénéité clinique et génétique dans la rétinite pigmentaire n'est pas seulement due aux différentes mutations dans des gènes différents, mais aussi à des effets différents de la même mutation chez des individus différents, parfois même dans la même famille. Mon travail de thèse s'est principalement axé sur une forme autosomique dominante de RP liée à des mutations dans le gène PRPF31, associées souvent à une pénétrance réduite, me conduisant à l'identification et à la caractérisation du mécanisme d'action du régulateur principal de la pénétrance des mutations: la protéine CNOT3. Dans la même logique d'étude des mécanismes moléculaires responsables de l'hétérogénéité clinique et génétique de la RP, nous avons étudié une forme récessive de la maladie associée à des mutations dans le gène récemment identifié FAMI61 A, dont les mutations dans le même gène donnent lieu à des manifestations cliniques différentes. Nos données ont ainsi accru la connaissance de la relation entre le génotype et le phénotype dans cette forme de maladie. La technique de séquençage du génome entier a été ensuite testée en tant que stratégie pour l'identification du gène de la maladie chez les patients atteints de RP récessive. Cette approche a montré son efficacité dans l'identification de variantes pathologiques qui n'auraient pu être détectées avec d'autres méthodes de dépistage. Enfin, pour la première fois, nous avons identifié une tumeur choroïdienne parmi les manifestations cliniques du PTEN hamartoma tumor syndrome, une maladie génétique causée par des mutations germinales du gène suppresseur de tumeur PTEN. Notre étude a mis en évidence l'hétérogénéité de cette tumeur choroïdienne, montrant que les altérations génétiques et/ou épigénétiques dans les différents gènes peuvent contribuer au développement et à la croissance tumorale.

Renal arteries: navigator-gated balanced fast field-echo projection MR angiography with aortic spin labeling: initial experience.

A cardiac-triggered free-breathing three-dimensional balanced fast field-echo projection magnetic resonance (MR) angiographic sequence with a two-dimensional pencil-beam aortic labeling pulse was developed for the renal arteries. For data acquisition during free breathing in eight healthy adults and seven consecutive patients with renal artery disease, real-time navigator technology was implemented. This technique allows high-spatial-resolution and high-contrast renal MR angiography and visualization of renal artery stenosis without exogenous contrast agent or breath hold. Initial promising results warrant larger clinical studies.

A fast 3D approach for coronary MRA.

Two-dimensional (2D)-breath-hold coronary magnetic resonance angiography (MRA) has been shown to be a fast and reliable method to depict the proximal coronary arteries. Recent developments, however, allow for free-breathing navigator gated and navigator corrected three-dimensional (3D) coronary MRA. These 3D approaches have potential for improved signal-to-noise ratio (SNR) and allow for the acquisition of adjacent thin slices without the misregistration problems known from 2D approaches. Still, a major impediment of a 3D acquisition is the increased scan time. The purpose of this study was the implementation of a free-breathing navigator gated and corrected ultra-fast 3D coronary MRA technique, which allows for scan times of less than 5 minutes. Twelve healthy adult subjects were examined in the supine position using a navigator gated and corrected ECG triggered ultra-fast 3D interleaved gradient echo planar imaging sequence (TFE-EPI). A 3D slab, consisting of 20 slices with a reconstructed slice thickness of 1.5 mm, was acquired with free-breathing. The diastolic TFE-EPI acquisition block was preceded by a T2prep pre-pulse, a diaphragmatic navigator pulse, and a fat suppression pre-pulse. With a TR of 19 ms and an effective TE of 5.4 ms, the duration of the data acquisition window duration was 38 ms. The in-plane spatial resolution was 1.0-1.3 mm*1.5-1.9 mm. In all cases, the entire left main (LM) and extensive portions of the left anterior descending (LAD) and right coronary artery (RCA) could be visualized with an average scan time for the entire 3D-volume data set of 2:57 +/- 0:51 minutes. Average contiguous vessel length visualized was 53 +/- 11 mm (range: 42 to 75 mm) for the LAD and 84 +/- 14 mm (range: 62 to 112 mm) for the RCA. Contrast-to-noise between coronary blood and myocardium was 5.0 +/- 2.3 for the LM/LAD and 8.0 +/- 2.9 for the RCA, resulting in an excellent suppression of myocardium. We present a new approach for free-breathing 3D coronary MRA, which allows for scan times superior to corresponding 2D coronary MRA approaches, and which takes advantage of the enhanced SNR of 3D acquisitions and the post-processing benefits of thin adjacent slices. The robust image quality and the short average scanning time suggest that this approach may be useful for screening the major coronary arteries or identification of anomalous coronary arteries. J. Magn. Reson. Imaging 1999;10:821-825.