Optimized quantitative magnetic resonance spectroscopy for clinical routine

Several practical obstacles in data handling and evaluation complicate the use of quantitative localized magnetic resonance spectroscopy (qMRS) in clinical routine MR examinations. To overcome these obstacles, a clinically feasible MR pulse sequence protocol based on standard available MR pulse sequences for qMRS has been implemented along with newly added functionalities to the free software package jMRUI-v5.0 to make qMRS attractive for clinical routine. This enables (a) easy and fast DICOM data transfer from the MR console and the qMRS-computer, (b) visualization of combined MR spectroscopy and imaging, (c) creation and network transfer of spectroscopy reports in DICOM format, (d) integration of advanced water reference models for absolute quantification, and (e) setup of databases containing normal metabolite concentrations of healthy subjects. To demonstrate the work-flow of qMRS using these implementations, databases for normal metabolite concentration in different regions of brain tissue were created using spectroscopic data acquired in 55 normal subjects (age range 6-61 years) using 1.5T and 3T MR systems, and illustrated in one clinical case of typical brain tumor (primitive neuroectodermal tumor). The MR pulse sequence protocol and newly implemented software functionalities facilitate the incorporation of qMRS and reference to normal value metabolite concentration data in daily clinical routine. Magn Reson Med, 2013. © 2012 Wiley Periodicals, Inc.

Development of a morphometric magnetic resonance image parameter suitable for distinguishing between normal dogs and dogs with cerebellar atrophy

Neurodegenerative diseases affect the cerebellum of numerous dog breeds. Although subjective, magnetic resonance (MR) imaging has been used to detect cerebellar atrophy in these diseases, but there are few data available on the normal size range of the cerebellum relative to other brain regions. The purpose of this study was to determine whether the size of the cerebellum maintains a consistent ratio with other brain regions in different ages and breeds of normal dogs and to define a measurement that can be used to identify cerebellar atrophy on MR images. Images from 52 normal and 13 dogs with cerebellar degenerative diseases were obtained. Volume and mid-sagittal cross-sectional area of the forebrain, brainstem, and cerebellum were calculated for each normal dog and compared between different breeds and ages as absolute and relative values. The ratio of the cerebellum to total brain and of the brainstem to cerebellum mid-sagittal cross-sectional area was compared between normal and affected dogs and the sensitivity and specificity of these ratios at distinguishing normal from affected dogs was calculated. The percentage of the brain occupied by the cerebellum in diverse dog breeds between 1 and 5 years of age was not significantly different, and cerebellar size did not change with increasing age. Using a cut off of 89%, the ratio between the brainstem and cerebellum mid-sagittal cross-sectional area could be used successfully to differentiate affected from unaffected dogs with a sensitivity and specificity of 100%, making this ratio an effective tool for identifying cerebellar atrophy on MR images.

COVARIATE-ADJUSTED NONPARAMETRIC ANALYSIS OF MAGNETIC RESONANCE IMAGES USING MARKOV CHAIN MONTE CARLO

Permutation tests are useful for drawing inferences from imaging data because of their flexibility and ability to capture features of the brain that are difficult to capture parametrically. However, most implementations of permutation tests ignore important confounding covariates. To employ covariate control in a nonparametric setting we have developed a Markov chain Monte Carlo (MCMC) algorithm for conditional permutation testing using propensity scores. We present the first use of this methodology for imaging data. Our MCMC algorithm is an extension of algorithms developed to approximate exact conditional probabilities in contingency tables, logit, and log-linear models. An application of our non-parametric method to remove potential bias due to the observed covariates is presented.

High-resolution three-dimensional 3 T magnetic resonance angiography for the evaluation of experimental aneurysm in the rabbit

OBJECTIVE: The standard technique of two-dimensional intra-arterial digital subtraction angiography (2D-DSA) for the imaging of experimental rabbit aneurysms is invasive and has considerable surgical risks. Therefore, minimally invasive techniques ideally providing three-dimensional imaging for intervention planning and follow-up are needed. This study evaluates the feasibility and quality of three-dimensional 3-T magnetic resonance angiography (3D-3T-MRA) and compares 3D-3T-MRA with 2D-DSA in experimental aneurysms in the rabbit. METHOD: Three microsurgically created aneurysms in three rabbits were evaluated using 2D-DSA and 3D-3T-MRA. Imaging of the aneurysms was performed 2 weeks after creation using 2D-DSA and contrast-enhanced (CE) MRA. Measurements included aneurysm dome (length and width) and aneurysm neck. Aneurysm volumes were determined using CE-MRA. RESULTS: The measurements of the aneurysms' dimensions and the evaluation of vicinity vessels with both techniques showed a good correlation. The mean aneurysm length, aneurysm width and neck width measured with DSA (6.9, 4.1 and 2.8 mm, respectively) correlated with the measurements performed in 3D-3T-MRA (6.9, 4 and 2.5 mm, respectively). The mean aneurysm volumes measured with CE-MRA was 46.7 mm(3). CONCLUSION: 3D-3T CE-MRA is feasible and less invasive and is a safer imaging alternative to DSA for experimental aneurysm. Additionally, aneurysm technique this precise offers the possibility of repetitive 3D aneurysm volumetry for long-term follow-up studies after endovascular aneurysm occlusion.

Magnetic resonance spectroscopy investigations of functionally defined language areas in schizophrenia patients with and without auditory hallucinations

Background: Cerebral dysfunction occurring in mental disorders can show metabolic disturbances which are limited to circumscribed brain areas. Auditory hallucinations have been shown to be related to defined cortical areas linked to specific language functions. Here, we investigated if the study of metabolic changes in auditory hallucinations requires a functional rather than an anatomical definition of their location and size to allow a reliable investigation by magnetic resonance spectroscopy (MRS). Methods: Schizophrenia patients with (AH; n = 12) and without hallucinations (NH; n = 8) and healthy controls (HC; n = 11) underwent a verbal fluency task in functional MRI (fMRI) to functionally define Broca's and Wernicke's areas. Left and right Heschl's gyri were defined anatomically. Results: The mean distances in native space between the fMRI-defined regions and a corresponding anatomically defined area were 12.4 ± 6.1 mm (range: 2.7–36.1 mm) for Broca's area and 16.8 ± 6.2 mm (range: 4.5–26.4 mm) for Wernicke's area, respectively. Hence, the spatial variance was of similar extent as the size of the investigated regions. Splitting the investigations into a single voxel examination in the frontal brain and a spectroscopic imaging part for the more homogeneous field areas led to good spectral quality for almost all spectra. In Broca's area, there was a significant group effect (p = 0.03) with lower levels of N-acetyl-aspartate (NAA) in NH compared to HC (p = 0.02). There were positive associations of NAA levels in the left Heschl's gyrus with total (p = 0.03) and negative (p = 0.006) PANSS scores. In Broca's area, there was a negative association of myo-inositol levels with total PANSS scores (p = 0.008). Conclusion: This study supports the neurodegenerative hypothesis of schizophrenia only in a frontal region whereas the results obtained from temporal regions are in contrast to the majority of previous studies. Future research should test the hypothesis raised by this study that a functional definition of language regions is needed if neurochemical imbalances are expected to be restricted to functional foci.

Improved cerebellar tissue classification on magnetic resonance images of brain.

PURPOSE: To develop and implement a method for improved cerebellar tissue classification on the MRI of brain by automatically isolating the cerebellum prior to segmentation. MATERIALS AND METHODS: Dual fast spin echo (FSE) and fluid attenuation inversion recovery (FLAIR) images were acquired on 18 normal volunteers on a 3 T Philips scanner. The cerebellum was isolated from the rest of the brain using a symmetric inverse consistent nonlinear registration of individual brain with the parcellated template. The cerebellum was then separated by masking the anatomical image with individual FLAIR images. Tissues in both the cerebellum and rest of the brain were separately classified using hidden Markov random field (HMRF), a parametric method, and then combined to obtain tissue classification of the whole brain. The proposed method for tissue classification on real MR brain images was evaluated subjectively by two experts. The segmentation results on Brainweb images with varying noise and intensity nonuniformity levels were quantitatively compared with the ground truth by computing the Dice similarity indices. RESULTS: The proposed method significantly improved the cerebellar tissue classification on all normal volunteers included in this study without compromising the classification in remaining part of the brain. The average similarity indices for gray matter (GM) and white matter (WM) in the cerebellum are 89.81 (+/-2.34) and 93.04 (+/-2.41), demonstrating excellent performance of the proposed methodology. CONCLUSION: The proposed method significantly improved tissue classification in the cerebellum. The GM was overestimated when segmentation was performed on the whole brain as a single object.

Age-related normal structural and functional ventricular values in cardiac function assessed by magnetic resonance

BACKGROUND The heart is subject to structural and functional changes with advancing age. However, the magnitude of cardiac age-dependent transformation has not been conclusively elucidated. METHODS This retrospective cardiac magnetic resonance (CMR) study included 183 subjects with normal structural and functional ventricular values. End systolic volume (ESV), end diastolic volume (EDV), and ejection fraction (EF) were obtained from the left and the right ventricle in breath-hold cine CMR. Patients were classified into four age groups (20-29, 30-49, 50-69, and ≥70 years) and cardiac measurements were compared using Pearson's rank correlation over the four different groups. RESULTS With advanced age a slight but significant decrease in ESV (r=-0.41 for both ventricles, P<0.001) and EDV (r=-0.39 for left ventricle, r=-0.35 for right ventricle, P<0.001) were observed associated with a significant increase in left (r=0.28, P<0.001) and right (r=0.27, P<0.01) ventricular EF reaching a maximal increase in EF of +8.4% (P<0.001) for the left and +6.1% (P<0.01) for the right ventricle in the oldest compared to the youngest patient group. Left ventricular myocardial mass significantly decreased over the four different age groups (P<0.05). CONCLUSIONS The aging process is associated with significant changes in left and right ventricular EF, ESV and EDV in subjects with no cardiac functional and structural abnormalities. These findings underline the importance of using age adapted values as standard of reference when evaluating CMR studies.

Relationship of vasodilator-induced changes in myocardial oxygenation with the severity of coronary artery stenosis: a study using oxygenation-sensitive cardiovascular magnetic resonance

Aims To explore the impact of the functional severity of coronary artery stenosis on changes in myocardial oxygenation during pharmacological vasodilation, using oxygenation-sensitive cardiovascular magnetic resonance (OS-CMR) imaging and invasive fractional flow reserve (FFR). An FFR is considered a standard of reference for assessing haemodynamic relevance of coronary artery stenosis; yet, the relationship of FFR to changes in myocardial oxygenation during vasodilator stress and thus to an objective marker for ischaemia on the tissue level is not well understood. Methods and results We prospectively recruited 64 patients with suspected/known coronary artery disease undergoing invasive angiography. The FFR was performed in intermediate coronary artery stenosis. OS-CMR images were acquired using a T2*-sensitive sequence before and after adenosine-induced vasodilation, with myocardial segments matched to angiography. Very strict image quality criteria were defined to ensure the validity of results. The FFR was performed in 37 patients. Because of the strict image quality criteria, 41% of segments had to be excluded, leaving 29/64 patients for the blinded OS-CMR analysis. Coronary territories with an associated FFR of <0.80 showed a lack of increase in myocardial oxygenation [mean signal intensity (ΔSI) −0.49%; 95% confidence interval (CI) −3.78 to 2.78 vs. +7.30%; 95% CI 4.08 to 10.64; P < 0.001]. An FFR of <0.54 best predicted a complete lack of a vasodilator-induced oxygenation increase (sensitivity 71% and specificity 75%). An OS-CMR ΔSI <4.78% identified an FFR of <0.8 with a sensitivity of 86% and specificity of 92%. Conclusion An FFR of <0.80 is associated with a lack of an adenosine-inducible increase in oxygenation of the dependent coronary territory, while a complete lack of such an increase was best predicted by an FFR of <0.54. Further studies are warranted to identify clinically meaningful cut-off values for FFR measurements and to assess the utility of OS-CMR as an alternative clinical tool for assessing the functional relevance of coronary artery stenosis.

A dual resonance, image -guided, volume localization technique for magnetic resonance spectroscopy

An interleaved, dual resonance, volume localization technique for $\sp1$H/$\sp{31}$P magnetic resonance spectroscopy has been designed, implemented on a 2 T imager/spectrometer, and verified with phantom studies.^ Localization techniques, including several single voxel techniques and spectroscopic imaging, were implemented, and studies were performed to compare the efficiency of each sequence of $\sp1$H/$\sp{31}$P spectral acquisitions. The sequence chosen was a hybrid of the stimulated echo single voxel technique and the spectroscopic imaging technique.^ Water suppression during the $\sp1$H spectral acquisitions was accomplished by the use of three narrow bandwidth RF saturation pulses in combination with three spoiler gradients. The spoiler gradient amplitudes were selected on the basis of a numerical solution of the Bloch equations. A post-acquisition water suppression algorithm was used to minimize any residual water signal.^ For interleaved $\sp1$H/$\sp{31}$P acquisitions, a dual resonance RF coil was constructed and interfaced to the existing RF detection system via a custom-designed dual resonance transcoupler and switching system. Programmable attenuators were incorporated to allow for changes in receiver and transmitter attenuation "on the fly".^ To provide the rapidly switched gradient fields required for the $\sp1$H/$\sp{31}$P acquisitions, an actively screened gradient coil system was designed and implemented. With this system, gradient field rise times on the order of 100 $\mu$s were obtained. These rapid switching times were necessary for minimizing intrasequence delays and for improving localization quality and water suppression efficiency.^ The interleaved $\sp1$H/$\sp{31}$P volume localization technique was tested using a two-compartment phantom. Analysis of the data showed that the spectral contamination was less than three percent. One-to-one spatial correspondence of the $\sp1$H and $\sp{31}$P spectra was verified and allowed for direct correlation of the spectral data with a standard magnetic resonance image. ^

Brain tissue segmentation on Diffusion Weighted Magnetic Resonance data

Introduction Diffusion weighted Imaging (DWI) techniques are able to measure, in vivo and non-invasively, the diffusivity of water molecules inside the human brain. DWI has been applied on cerebral ischemia, brain maturation, epilepsy, multiple sclerosis, etc. [1]. Nowadays, there is a very high availability of these images. DWI allows the identification of brain tissues, so its accurate segmentation is a common initial step for the referred applications. Materials and Methods We present a validation study on automated segmentation of DWI based on the Gaussian mixture and hidden Markov random field models. This methodology is widely solved with iterative conditional modes algorithm, but some studies suggest [2] that graph-cuts (GC) algorithms improve the results when initialization is not close to the final solution. We implemented a segmentation tool integrating ITK with a GC algorithm [3], and a validation software using fuzzy overlap measures [4]. Results Segmentation accuracy of each tool is tested against a gold-standard segmentation obtained from a T1 MPRAGE magnetic resonance image of the same subject, registered to the DWI space. The proposed software shows meaningful improvements by using the GC energy minimization approach on DTI and DSI (Diffusion Spectrum Imaging) data. Conclusions The brain tissues segmentation on DWI is a fundamental step on many applications. Accuracy and robustness improvements are achieved with the proposed software, with high impact on the application’s final result.

Noninvasive monitoring of serial changes in pulmonary vascular resistance and acute vasodilator testing using cardiac magnetic resonance

Objectives The study sought to evaluate the ability of cardiac magnetic resonance (CMR) to monitor acute and long-term changes in pulmonary vascular resistance (PVR) noninvasively. Background PVR monitoring during the follow-up of patients with pulmonary hypertension (PH) and the response to vasodilator testing require invasive right heart catheterization. Methods An experimental study in pigs was designed to evaluate the ability of CMR to monitor: 1) an acute increase in PVR generated by acute pulmonary embolization (n = 10); 2) serial changes in PVR in chronic PH (n = 22); and 3) changes in PVR during vasodilator testing in chronic PH (n = 10). CMR studies were performed with simultaneous hemodynamic assessment using a CMR-compatible Swan-Ganz catheter. Average flow velocity in the main pulmonary artery (PA) was quantified with phase contrast imaging. Pearson correlation and mixed model analysis were used to correlate changes in PVR with changes in CMR-quantified PA velocity. Additionally, PVR was estimated from CMR data (PA velocity and right ventricular ejection fraction) using a formula previously validated. Results Changes in PA velocity strongly and inversely correlated with acute increases in PVR induced by pulmonary embolization (r = –0.92), serial PVR fluctuations in chronic PH (r = –0.89), and acute reductions during vasodilator testing (r = –0.89, p ≤ 0.01 for all). CMR-estimated PVR showed adequate agreement with invasive PVR (mean bias –1.1 Wood units,; 95% confidence interval: –5.9 to 3.7) and changes in both indices correlated strongly (r = 0.86, p < 0.01). Conclusions CMR allows for noninvasive monitoring of acute and chronic changes in PVR in PH. This capability may be valuable in the evaluation and follow-up of patients with PH.

Detection of neuritic plaques in Alzheimer’s disease by magnetic resonance microscopy

Magnetic resonance microscopy (MRM) theoretically provides the spatial resolution and signal-to-noise ratio needed to resolve neuritic plaques, the neuropathological hallmark of Alzheimer’s disease (AD). Two previously unexplored MR contrast parameters, T2* and diffusion, are tested for plaque-specific contrast to noise. Autopsy specimens from nondemented controls (n = 3) and patients with AD (n = 5) were used. Three-dimensional T2* and diffusion MR images with voxel sizes ranging from 3 × 10−3 mm3 to 5.9 × 10−5 mm3 were acquired. After imaging, specimens were cut and stained with a microwave king silver stain to demonstrate neuritic plaques. From controls, the alveus, fimbria, pyramidal cell layer, hippocampal sulcus, and granule cell layer were detected by either T2* or diffusion contrast. These structures were used as landmarks when correlating MRMs with histological sections. At a voxel resolution of 5.9 × 10−5 mm3, neuritic plaques could be detected by T2*. The neuritic plaques emerged as black, spherical elements on T2* MRMs and could be distinguished from vessels only in cross-section when presented in three dimension. Here we provide MR images of neuritic plaques in vitro. The MRM results reported provide a new direction for applying this technology in vivo. Clearly, the ability to detect and follow the early progression of amyloid-positive brain lesions will greatly aid and simplify the many possibilities to intervene pharmacologically in AD.

Development of a magnetic resonance compatible wrist device for the analysis of movement encoding in the brain

Tese de mestrado integrado em Engenharia Biomédica e Biofísica, apresentada à Universidade de Lisboa, através da Faculdade de Ciências, 2016

3-D magnetic resonance angiography versus conventional angiography in peripheral arterial disease: Pilot study

Angiography is usually performed as the preoperative road map for those requiring revascularization for lower extremity peripheral arterial disease (PAD). The alternative investigations are ultrasound, 3-D magnetic resonance angiography (3-D MRA) and computed tomography angiography. This pilot study aimed to assess whether 3-D MRA could replace the gold standard angiography in preoperative planning. Eight patients considered for aortoiliac or infrainguinal arterial bypass surgery were recruited. All underwent both imaging modalities within 7 days. A vascular surgeon and a radiologist each reported on the images from both the 3-D MRA and the angiography, with blinding to patient details and each others reports. Comparisons were made between the reports for the angiographic and the 3-D MRA images, and between the reports of the vascular surgeon and the radiologist. Compared to the gold standard angiogram, 3-D MRA had a sensitivity of 77% and specificity of 94% in detecting occlusion, and a sensitivity of 72% and specificity of 90% in differentiating high grade (> 50%) versus low grade (< 50%) stenoses. There was an overall concordance of 78% between the two investigations with a range of 62% in the peroneal artery to 94% in the aorta. 3-D MRA showed flow in 23% of cases where conventional angiography showed no flow. In the present pilot study, 3-D MRA had reasonable concordance with the gold standard angiography, depending on the level of the lesion. At times it showed vessel flow where occlusion was shown on conventional angiogram. 3-D MRA in peripheral vascular disease is challenging the gold standard, but is inconsistent at present.

Tensor Completion for Multidimensional Inverse Problems with Applications to Magnetic Resonance Relaxometry

This thesis deals with tensor completion for the solution of multidimensional inverse problems. We study the problem of reconstructing an approximately low rank tensor from a small number of noisy linear measurements. New recovery guarantees, numerical algorithms, non-uniform sampling strategies, and parameter selection algorithms are developed. We derive a fixed point continuation algorithm for tensor completion and prove its convergence. A restricted isometry property (RIP) based tensor recovery guarantee is proved. Probabilistic recovery guarantees are obtained for sub-Gaussian measurement operators and for measurements obtained by non-uniform sampling from a Parseval tight frame. We show how tensor completion can be used to solve multidimensional inverse problems arising in NMR relaxometry. Algorithms are developed for regularization parameter selection, including accelerated k-fold cross-validation and generalized cross-validation. These methods are validated on experimental and simulated data. We also derive condition number estimates for nonnegative least squares problems. Tensor recovery promises to significantly accelerate N-dimensional NMR relaxometry and related experiments, enabling previously impractical experiments. Our methods could also be applied to other inverse problems arising in machine learning, image processing, signal processing, computer vision, and other fields.