Automated multiclass segmentation, quantification, and visualization of the diseased aorta on hybrid PET/CT-SEQUOIA

BACKGROUND

Cardiovascular disease is the most common cause of death worldwide, including infection and inflammation related conditions. Multiple studies have demonstrated potential advantages of hybrid positron emission tomography combined with computed tomography (PET/CT) as an adjunct to current clinical inflammatory and infectious biochemical markers. To quantitatively analyze vascular diseases at PET/CT, robust segmentation of the aorta is necessary. However, manual segmentation is extremely time-consuming and labor-intensive.

PURPOSE

To investigate the feasibility and accuracy of an automated tool to segment and quantify multiple parts of the diseased aorta on unenhanced low-dose computed tomography (LDCT) as an anatomical reference for PET-assessed vascular disease.

METHODS

A software pipeline was developed including automated segmentation using a 3D U-Net, calcium scoring, PET uptake quantification, background measurement, radiomics feature extraction, and 2D surface visualization of vessel wall calcium and tracer uptake distribution. To train the 3D U-Net, 352 non-contrast LDCTs from (2-[18 F]FDG and Na[18 F]F) PET/CTs performed in patients with various vascular pathologies with manual segmentation of the ascending aorta, aortic arch, descending aorta, and abdominal aorta were used. The last 22 consecutive scans were used as a hold-out internal test set. The remaining dataset was randomly split into training (n = 264; 80%) and validation (n = 66; 20%) sets. Further evaluation was performed on an external test set of 49 PET/CTs. The dice similarity coefficient (DSC) and Hausdorff distance (HD) were used to assess segmentation performance. Automatically obtained calcium scores and uptake values were compared with manual scoring obtained using clinical softwares (syngo.via and Affinity Viewer) in six patient images. intraclass correlation coefficients (ICC) were calculated to validate calcium and uptake values.

RESULTS

Fully automated segmentation of the aorta using a 3D U-Net was feasible in LDCT obtained from PET/CT scans. The external test set yielded a DSC of 0.867 ± 0.030 and HD of 1.0 [0.6-1.4] mm, similar to an open-source model with a DSC of 0.864 ± 0.023 and HD of 1.4 [1.0-1.8] mm. Quantification of calcium and uptake values were in excellent agreement with clinical software (ICC: 1.00 [1.00-1.00] and 0.99 [0.93-1.00] for calcium and uptake values, respectively).

CONCLUSIONS

We present an automated pipeline to segment the ascending aorta, aortic arch, descending aorta, and abdominal aorta on LDCT from PET/CT and to accurately provide uptake values, calcium scores, background measurement, radiomics features, and a 2D visualization. We call this algorithm SEQUOIA (SEgmentation, QUantification, and visualizatiOn of the dIseased Aorta) and is available at https://github.com/UMCG-CVI/SEQUOIA. This model could augment the utility of aortic evaluation at PET/CT studies tremendously, irrespective of the tracer, and potentially provide fast and reliable quantification of cardiovascular diseases in clinical practice, both for primary diagnosis and disease monitoring.

Evaluation of the impact of strain correction on the orientation of cardiac diffusion tensors with in vivo and ex vivo porcine hearts

PURPOSE

To evaluate the importance of strain-correcting stimulated echo acquisition mode echo-planar imaging cardiac diffusion tensor imaging.

METHODS

Healthy pigs (n = 11) were successfully scanned with a 3D cine displacement-encoded imaging with stimulated echoes and a monopolar-stimulated echo-planar imaging diffusion tensor imaging sequence at 3 T during diastasis, peak systole, and strain sweet spots in a midventricular short-axis slice. The same diffusion tensor imaging sequence was repeated ex vivo after arresting the hearts in either a relaxed (KCl-induced) or contracted (BaCl₂ -induced) state. The displacement-encoded imaging with stimulated echoes data were used to strain-correct the in vivo cardiac diffusion tensor imaging in diastole and systole. The orientation of the primary (helix angles) and secondary (E2A) diffusion eigenvectors was compared with and without strain correction and to the strain-free ex vivo data.

RESULTS

Strain correction reduces systolic E2A significantly when compared without strain correction and ex vivo (median absolute E2A = 34.3° versus E2A = 57.1° (P = 0.01), E2A = 60.5° (P = 0.006), respectively). The systolic distribution of E2A without strain correction is closer to the contracted ex vivo distribution than with strain correction, root mean square deviation of 0.027 versus 0.038.

CONCLUSIONS

The current strain-correction model amplifies the contribution of microscopic strain to diffusion resulting in an overcorrection of E2A. Results show that a new model that considers cellular rearrangement is required. Magn Reson Med 79:2205-2215, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Effect of inversion time on the precision of myocardial late gadolinium enhancement quantification evaluated with synthetic inversion recovery MR imaging

OBJECTIVES

To evaluate the influence of inversion time (TI) on the precision of myocardial late gadolinium enhancement (LGE) quantification using synthetic inversion recovery (IR) imaging in patients with myocardial infarction (MI).

METHODS

Fifty-three patients with suspected prior MI underwent 1.5-T cardiac MRI with conventional magnitude (MagIR) and phase-sensitive IR (PSIR) LGE imaging and T1 mapping at 15 min post-contrast. T1-based synthetic MagIR and PSIR images were calculated with a TI ranging from -100 to +150 ms at 5-ms intervals relative to the optimal TI (TI₀). LGE was quantified using a five standard deviation (5SD) and full width at half-maximum (FWHM) thresholds. Measurements were compared using one-way analysis of variance.

RESULTS

The MagIR_sy technique provided precise assessment of LGE area at TIs ≥ TI₀, while precision was decreased below TI₀. The LGE area showed significant differences at ≤ -25 ms compared to TI₀ using 5SD (P < 0.001) and at ≤ -65 ms using the FWHM approach (P < 0.001). LGE measurements did not show significant difference over the analysed TI range in the PSIR_sy images using either of the quantification methods.

CONCLUSIONS

T1 map-based PSIR_sy images provide precise quantification of MI independent of TI at the investigated time point post-contrast. MagIR_sy-based MI quantification is precise at TI₀ and at longer TIs while showing decreased precision at TI values below TI₀.

KEY POINTS

• Synthetic IR imaging retrospectively generates LGE images at any theoretical TI • Synthetic IR imaging can simulate the effect of TI on LGE quantification • Fifteen minutes post-contrast MagIR _sy accurately quantifies infarcts from TI ₀ to TI ₀   + 150 ms • Fifteen minutes post-contrast PSIR _sy provides precise infarct size independent of TI • Synthetic IR imaging has further advantages in reducing operator dependence.

Validation of T2* in-line analysis for tissue iron quantification at 1.5 T

BACKGROUND

There is a need for improved worldwide access to tissue iron quantification using T2* cardiovascular magnetic resonance (CMR). One route to facilitate this would be simple in-line T2* analysis widely available on MR scanners. We therefore compared our clinically validated and established T2* method at Royal Brompton Hospital (RBH T2*) against a novel work-in-progress (WIP) sequence with in-line T2* measurement from Siemens (WIP T2*).

METHODS

Healthy volunteers (n = 22) and patients with iron overload (n = 78) were recruited (53 males, median age 34 years). A 1.5 T study (Magnetom Avanto, Siemens) was performed on all subjects. The same mid-ventricular short axis cardiac slice and transaxial slice through the liver were used to acquire both RBH T2* images and WIP T2* maps for each participant. Cardiac white blood (WB) and black blood (BB) sequences were acquired. Intraobserver, interobserver and interstudy reproducibility were measured on the same data from a subset of 20 participants.

RESULTS

Liver T2* values ranged from 0.8 to 35.7 ms (median 5.1 ms) and cardiac T2* values from 6.0 to 52.3 ms (median 31 ms). The coefficient of variance (CoV) values for direct comparison of T2* values by RBH and WIP were 6.1-7.8 % across techniques. Accurate delineation of the septum was difficult on some WIP T2* maps due to artefacts. The inability to manually correct for noise by truncation of erroneous later echo times led to some overestimation of T2* using WIP T2* compared with the RBH T2*. Reproducibility CoV results for RBH T2* ranged from 1.5 to 5.7 % which were better than the reproducibility of WIP T2* values of 4.1-16.6 %.

CONCLUSIONS

Iron estimation using the T2* CMR sequence in combination with Siemens' in-line data processing is generally satisfactory and may help facilitate global access to tissue iron assessment. The current automated T2* map technique is less good for tissue iron assessment with noisy data at low T2* values.

Saturation pulse design for quantitative myocardial T1 mapping

BACKGROUND

Quantitative saturation-recovery based T1 mapping sequences are less sensitive to systematic errors than the Modified Look-Locker Inversion recovery (MOLLI) technique but require high performance saturation pulses. We propose to optimize adiabatic and pulse train saturation pulses for quantitative T1 mapping to have <1 % absolute residual longitudinal magnetization (|MZ/M0|) over ranges of B0 and [Formula: see text] (B1 scale factor) inhomogeneity found at 1.5 T and 3 T.

METHODS

Design parameters for an adiabatic BIR4-90 pulse were optimized for improved performance within 1.5 T B0 (±120 Hz) and [Formula: see text] (0.7-1.0) ranges. Flip angles in hard pulse trains of 3-6 pulses were optimized for 1.5 T and 3 T, with consideration of T1 values, field inhomogeneities (B0 = ±240 Hz and [Formula: see text]=0.4-1.2 at 3 T), and maximum achievable B1 field strength. Residual MZ/M0 was simulated and measured experimentally for current standard and optimized saturation pulses in phantoms and in-vivo human studies. T1 maps were acquired at 3 T in human subjects and a swine using a SAturation recovery single-SHot Acquisition (SASHA) technique with a standard 90°-90°-90° and an optimized 6-pulse train.

RESULTS

Measured residual MZ/M0 in phantoms had excellent agreement with simulations over a wide range of B0 and [Formula: see text]. The optimized BIR4-90 reduced the maximum residual |MZ/M0| to <1 %, a 5.8× reduction compared to a reference BIR4-90. An optimized 3-pulse train achieved a maximum residual |MZ/M0| <1 % for the 1.5 T optimization range compared to 11.3 % for a standard 90°-90°-90° pulse train, while a 6-pulse train met this target for the wider 3 T ranges of B0 and [Formula: see text]. The 6-pulse train demonstrated more uniform saturation across both the myocardium and entire field of view than other saturation pulses in human studies. T1 maps were more spatially homogeneous with 6-pulse train SASHA than the reference 90°-90°-90° SASHA in both human and animal studies.

CONCLUSIONS

Adiabatic and pulse train saturation pulses optimized for different constraints found at 1.5 T and 3 T achieved <1 % residual |MZ/M0| in phantom experiments, enabling greater accuracy in quantitative saturation recovery T1 imaging.

Myocardial Late Gadolinium Enhancement: Accuracy of T1 Mapping-based Synthetic Inversion-Recovery Imaging

PURPOSE

To compare the accuracy of detection and quantification of myocardial late gadolinium enhancement (LGE) with a synthetic inversion-recovery (IR) approach with that of conventional IR techniques.

MATERIALS AND METHODS

This prospective study was approved by the institutional review board and compliant with HIPAA. All patients gave written informed consent. Between June and November 2014, 43 patients (25 men; mean age, 54 years ± 16) suspected of having previous myocardial infarction underwent magnetic resonance (MR) imaging, including contrast material-enhanced LGE imaging and T1 mapping. Synthetic magnitude and phase-sensitive IR images were generated on the basis of T1 maps. Images were assessed by two readers. Differences in the per-patient and per-segment LGE detection rates between the synthetic and conventional techniques were analyzed with the McNemar test, and the accuracy of LGE quantification was calculated with the paired t test and Bland-Altman statistics. Interreader agreement for the detection and quantification of LGE was analyzed with κ and Bland-Altman statistics, respectively.

RESULTS

Seventeen of the 43 patients (39%) had LGE patterns consistent with myocardial infarction. The sensitivity and specificity of synthetic magnitude and phase-sensitive IR techniques in the detection of LGE were 90% and 95%, respectively, with patient-based analysis and 94% and 99%, respectively, with segment-based analysis. The area of LGE measured with synthetic IR techniques showed excellent agreement with that of conventional techniques (4.35 cm(2) ± 1.88 and 4.14 cm(2)± 1.62 for synthetic magnitude and phase-sensitive IR, respectively, compared with 4.25 cm(2) ± 1.92 and 4.22 cm(2) ± 1.86 for conventional magnitude and phase-sensitive IR, respectively; P > .05). Interreader agreement was excellent for the detection (κ > 0.81) and quantification (bias range, -0.34 to 0.40; P > .05) of LGE.

CONCLUSION

The accuracy of the T1 map-based synthetic IR approach in the detection and quantification of myocardial LGE in patients with previous myocardial infarction was similar to that of conventional IR techniques. The use of T1 mapping to derive synthetic LGE images may reduce imaging times and operator dependence in future T1 mapping protocols with full left ventricular coverage.

White matter integrity of the cerebellar peduncles as a mediator of effects of prenatal alcohol exposure on eyeblink conditioning

Fetal alcohol spectrum disorders (FASD) are characterized by a range of neurodevelopmental deficits that result from prenatal exposure to alcohol. These can include cognitive, behavioural, and neurological impairment, as well as structural and functional brain damage. Eyeblink conditioning (EBC) is among the most sensitive endpoints affected in FASD. The cerebellar peduncles, large bundles of myelinated nerve fibers that connect the cerebellum to the brainstem, constitute the principal white matter element of the EBC circuit. Diffusion tensor imaging (DTI) is used to assess white matter integrity in fibre pathways linking brain regions. DTI scans of 54 children with FASD and 23 healthy controls, mean age 10.1 ± 1.0 years, from the Cape Town Longitudinal Cohort were processed using voxelwise group comparisons. Prenatal alcohol exposure was related to lower fractional anisotropy (FA) bilaterally in the superior cerebellar peduncles and higher mean diffusivity (MD) in the left middle peduncle, effects that remained significant after controlling for potential confounding variables. Lower FA and higher MD in these regions were associated with poorer EBC performance. Moreover, effects of alcohol exposure on EBC decreased significantly after inclusion of these DTI measures in regression models, suggesting that these white matter deficits partially mediate the relation of prenatal alcohol exposure to EBC. The associations of greater alcohol consumption with these DTI measures are largely attributable to greater radial diffusivity, possibly indicating poorer myelination. Thus, these data suggest that fetal alcohol-related deficits in EBC are attributable, in part, to poorer myelination in key regions of the cerebellar peduncles.

Free-breathing T2* mapping using respiratory motion corrected averaging

BACKGROUND

Pixel-wise T2* maps based on breath-held segmented image acquisition are prone to ghost artifacts in instances of poor breath-holding or cardiac arrhythmia. Single shot imaging is inherently immune to ghost type artifacts. We propose a free-breathing method based on respiratory motion corrected single shot imaging with averaging to improve the signal to noise ratio.

METHODS

Images were acquired using a multi-echo gradient recalled echo sequence and T2* maps were calculated at each pixel by exponential fitting. For 40 subjects (2 cohorts), two acquisition protocols were compared: (1) a breath-held, segmented acquisition, and (2) a free-breathing, single-shot multiple repetition respiratory motion corrected average. T2* measurements in the interventricular septum and liver were compared for the 2-methods in all studies with diagnostic image quality.

RESULTS

In cohort 1 (N = 28) with age 51.4 ± 17.6 (m ± SD) including 1 subject with severe myocardial iron overload, there were 8 non-diagnostic breath-held studies due to poor image quality resulting from ghost artifacts caused by respiratory motion or arrhythmias. In cohort 2 (N = 12) with age 30.9 ± 7.5 (m ± SD), including 7 subjects with severe myocardial iron overload and 4 subjects with mild iron overload, a single subject was unable to breath-hold. Free-breathing motion corrected T2* maps were of diagnostic quality in all 40 subjects. T2* measurements were in excellent agreement (In cohort #1, T2*FB = 0.95 x T2*BH + 0.41, r2 = 0.93, N = 39 measurements, and in cohort #2, T2*FB = 0.98 x T2*BH + 0.05, r2 > 0.99, N = 22 measurements).

CONCLUSIONS

A free-breathing approach to T2* mapping is demonstrated to produce consistently good quality maps in the presence of respiratory motion and arrhythmias.

Optimized saturation recovery protocols for T1-mapping in the heart: influence of sampling strategies on precision

BACKGROUND

T1-mapping has the potential to detect and quantify diffuse processes such as interstitial fibrosis. Detection of disease at an early stage by measurement of subtle changes requires a high degree of reproducibility. Initial implementation of saturation recovery (SR) T1-mapping employed 3-parameter fitting which was highly accurate but was quite sensitive to noise; 2-parameter fitting greatly reduced the sensitivity to noise at the expense of a small degree of systematic bias. A recently introduced implementation that uses a variable readout flip angle greatly reduces systematic errors in T1-measurement thereby making it feasible to use SR methods with 2-parameter fitting with improved accuracy and precision. SR T1 mapping techniques with multi-heartbeat recovery times have been proposed to better sample the T1 recovery curve, but have not been evaluated for 2-parameter fitting.

METHODS

An analytic formulation for calculating the standard deviation (SD) for SR T1-mapping with 2-parameter fitting is developed and validated using Monte-Carlo simulation. The coefficient of variation is compared for a brute force optimization of sampling and for several previously described sampling schemes for T1 measurement over several uncertainty ranges. Experimental validation is performed in phantoms over a range of T1, and in-vivo both native and post-contrast. Pixel-wise SD maps are calculated for SR T1-mapping.

RESULTS

Sampling schemes that use a non-saturated anchor image and multiple (N) measurements at a single fixed saturation delay are found to be near optimum for the case of known T1 and are close to the brute force optimized solution over wide ranges of native and post-contrast T1 values. The fixed delay sampling scheme is simple to implement and provides an improvement over uniformly distributed schemes.

CONCLUSIONS

Sampling strategies for saturation recovery methods for myocardial T1-mapping have been optimized and validated experimentally. Reduced SD, or improved precision, may be achieved by using fixed saturation delays when considering native myocardium and post-contrast T1 ranges. Pixel-wise estimates of T1 mapping errors have been formulated and validated for SR fitting methods.

Semi-automated left ventricular segmentation based on a guide point model approach for 3D cine DENSE cardiovascular magnetic resonance

BACKGROUND

The most time consuming and limiting step in three dimensional (3D) cine displacement encoding with stimulated echoes (DENSE) MR image analysis is the demarcation of the left ventricle (LV) from its surrounding anatomical structures. The aim of this study is to implement a semi-automated segmentation algorithm for 3D cine DENSE CMR using a guide point model approach.

METHODS

A 3D mathematical model is fitted to guide points which were interactively placed along the LV borders at a single time frame. An algorithm is presented to robustly propagate LV epicardial and endocardial surfaces of the model using the displacement information encoded in the phase images of DENSE data. The accuracy, precision and efficiency of the algorithm are tested.

RESULTS

The model-defined contours show good accuracy when compared to the corresponding manually defined contours as similarity coefficients Dice and Jaccard consist of values above 0.7, while false positive and false negative measures show low percentage values. This is based on a measure of segmentation error on intra- and inter-observer spatial overlap variability. The segmentation algorithm offers a 10-fold reduction in the time required to identify LV epicardial and endocardial borders for a single 3D DENSE data set.

CONCLUSION

A semi-automated segmentation method has been developed for 3D cine DENSE CMR. The algorithm allows for contouring of the first cardiac frame where blood-myocardium contrast is almost nonexistent and reduces the time required to segment a 3D DENSE data set significantly.

Micro-structurally detailed model of a therapeutic hydrogel injectate in a rat biventricular cardiac geometry for computational simulations

Biomaterial injection-based therapies have showed cautious success in restoration of cardiac function and prevention of adverse remodelling into heart failure after myocardial infarction (MI). However, the underlying mechanisms are not well understood. Computational studies utilised simplified representations of the therapeutic myocardial injectates. Wistar rats underwent experimental infarction followed by immediate injection of polyethylene glycol hydrogel in the infarct region. Hearts were explanted, cryo-sectioned and the region with the injectate histologically analysed. Histological micrographs were used to reconstruct the dispersed hydrogel injectate. Cardiac magnetic resonance imaging data from a healthy rat were used to obtain an end-diastolic biventricular geometry which was subsequently adjusted and combined with the injectate model. The computational geometry of the injectate exhibited microscopic structural details found the in situ. The combination of injectate and cardiac geometry provides realistic geometries for multiscale computational studies of intra-myocardial injectate therapies for the rat model that has been widely used for MI research.

Preoperative three-dimensional model creation of magnetic resonance brain images as a tool to assist neurosurgical planning

BACKGROUND

Neurosurgeons regularly plan their surgery using magnetic resonance imaging (MRI) images, which may show a clear distinction between the area to be resected and the surrounding healthy brain tissue depending on the nature of the pathology. However, this distinction is often unclear with the naked eye during the surgical intervention, and it may be difficult to infer depth and an accurate volumetric interpretation from a series of MRI image slices.

OBJECTIVES

In this work, MRI data are used to create affordable patient-specific 3-dimensional (3D) scale models of the brain which clearly indicate the location and extent of a tumour relative to brain surface features and important adjacent structures.

METHODS

This is achieved using custom software and rapid prototyping. In addition, functionally eloquent areas identified using functional MRI are integrated into the 3D models.

RESULTS

Preliminary in vivo results are presented for 2 patients. The accuracy of the technique was estimated both theoretically and by printing a geometrical phantom, with mean dimensional errors of less than 0.5 mm observed.

CONCLUSIONS

This may provide a practical and cost-effective tool which can be used for training, and during neurosurgical planning and intervention.

Relationship between neurocognition and regional brain volumes in traumatized adolescents with and without posttraumatic stress disorder

OBJECTIVES

Studies using convergent neurocognitive and structural imaging paradigms in adolescent posttraumatic stress disorder (PTSD) are limited; in the current study we used both voxel-based morphometry (VBM) to obtain between-group volumetric differences, and Freesurfer to examine the relationship between cognition and regional brain volumes.

METHODS

Participants were 21 traumatized adolescents with PTSD matched with 32 traumatized adolescents without PTSD. Magnetic resonance images were obtained on a 1.5-Tesla MAGNETOM Siemens Symphony scanner. VBM implemented on FSL was then used to compare between-group grey matter volumes, after which Freesurfer was used to obtain global volume and thickness measurements in different brain regions.

RESULTS

Significant between-group neurocognitive differences were found for tests of attention, delayed recall and visual reconstruction. On VBM, reduced grey matter was found in three regions in the PTSD group: left insula, right precuneus and right cingulate gyrus, using uncorrected values (p < 0.001), while no statistically significant between-group differences were found on the initial Freesurfer stream. Further Freesurfer analysis on Qdec revealed significant reductions in the insula for the PTSD group. In addition, volumetric changes in the corpus callosum and insula were significantly associated with deficits in logical memory and visual reproduction on Freesurfer analysis.

CONCLUSIONS

Trauma exposure of itself may be sufficient to cause structural changes in adolescents regardless of PTSD development.

A semi-automated method for measuring thickness and white matter integrity of the corpus callosum

Aim. Diseases affecting cerebral white matter may lead to left-right asymmetries and atrophy of interhemispheric connections, i.e. the corpus callosum (CC). Our aim was to describe and test a semi-automated system that divides the midline CC into a number of segments and determines thickness at each, then performs fibre tracking from these segments. Methods. Six normal female volunteers (average age 25.8 ±6.7 years) and a female patient with diagnosed multiple sclerosis (age 26 years) were scanned on a 3T MRI. We performed diffusion-weighted imaging in 12 directions, and calculated diffusion tensors and fractional anisotropy (FA) maps from this pre-processed data. Fibre tracking from a region-of-interest encompassing the entire CC was done. This fibre data, together with FA maps and the unweighted diffusion tensor imaging (DTI) image (b = 0 s/mm2), were imported into a custom tool written in MATLAB. The midline sagittal position was carefully defined by selecting multiple midline points in coronal and axial views and rotating the image volume and fibre co-ordinates accordingly. Using the customised tool, dorsal and ventral CC contours were manually drawn on the mid-sagittal FA image, initiating automated calculation of a contour midway between these manually drawn lines. The programme was designed to then divide the midline contour into a pre-selected number of segments; from each segment border, perpendicular spokes were projected until they intersected with the dorsal and ventral contours. This technique divided the CC into a pre-set amount of segments, the number of which was limited by the spatial resolution. It was decided to set the number at 40 to ensure that each segment depicted a contiguous strip of voxels across the CC from the dorsal to the ventral contour. The system allows these segments to then be used as seeds for separate fibre tracking in each cerebral hemisphere, and various parameters are automatically plotted as a function of distance along the midline contour. The following parameters are measurable: midline CC thickness; midline FA; fibre volume for each hemisphere (represented as a left/right ratio centred on zero) and mean fibre FA for each hemisphere (also represented as a left/right ratio centred on zero). Results. The tool proved successful in measuring and plotting CC midline thickness and FA, but was not sensitive for peripheral white matter lesions. Conclusions. The technique successfully determined values of CC midline thickness, FA and interhemispheric differences. Future research will determine normal values for age and compare CC thickness with peripheral white matter volume loss in large groups of patients, using the semiautomated technique.

Functional MRI in pre-surgical planning: Case study and cautionary notes

Background. Since its inception almost 20 years ago, functional magnetic resonance imaging (fMRI) has greatly advanced our knowledge of human brain function. Although the clinical applications of fMRI are still limited, there have recently been encouraging advances for its use in pre-operative functional cortical mapping to identify potentially eloquent areas prior to neurosurgery. Objectives. We explore the potential use of this emerging technique by presenting a neurosurgical case study, as performed at the Cape Universities Brain Imaging Centre (CUBIC), Tygerberg, Cape Town. We conclude with a brief summary of the potential pitfalls of this technique, as well as cautionary guidelines based on our experience. Methods and results. A 22-year-old male patient from Tygerberg Hospital underwent the successful resection of an anaplastic astrocytoma after fMRI presurgical planning at our facility. The subject was able to leave the ward unassisted. Conclusion. If consideration is given to the many limitations of this emerging technique, fMRI can be useful in aiding the neurosurgeon in pre-operative planning of his surgical approach.

Mapping right ventricular myocardial mechanics using 3D cine DENSE cardiovascular magnetic resonance

BACKGROUND

The mechanics of the right ventricle (RV) are not well understood as studies of the RV have been limited. This is, in part, due to the RV's thin wall, asymmetric geometry and irregular motion. However, the RV plays an important role in cardiovascular function. This study aims to describe the complex mechanics of the healthy RV using three dimensional (3D) cine displacement encoding with stimulated echoes (DENSE) cardiovascular magnetic resonance (CMR).

METHODS

Whole heart 3D cine DENSE data were acquired from five healthy volunteers. Tailored post-processing algorithms for RV mid-wall tissue tracking and strain estimation are presented. A method for sub-dividing the RV into four regions according to anatomical land marks is proposed, and the temporal evolution of strain was assessed in these regions.

RESULTS

The 3D cine DENSE tissue tracking methods successfully capture the motion and deformation of the RV at a high spatial resolution in all volunteers. The regional Lagrangian peak surface strain and time to peak values correspond with previous studies using myocardial tagging, DENSE and strain encoded CMR. The inflow region consistently displays lower peak strains than the apical and outflow regions, and the time to peak strains suggest RV mechanical activation in the following order: inflow, outflow, mid, then apex.

CONCLUSIONS

Model-free techniques have been developed to study the myocardial mechanics of the RV at a high spatial resolution using 3D cine DENSE CMR. The consistency of the regional RV strain patterns across healthy subjects is encouraging and the techniques may have clinical utility in assessing disrupted RV mechanics in the diseased heart.

Comprehensive cardiovascular magnetic resonance of myocardial mechanics in mice using three-dimensional cine DENSE

BACKGROUND

Quantitative noninvasive imaging of myocardial mechanics in mice enables studies of the roles of individual genes in cardiac function. We sought to develop comprehensive three-dimensional methods for imaging myocardial mechanics in mice.

METHODS

A 3D cine DENSE pulse sequence was implemented on a 7T small-bore scanner. The sequence used three-point phase cycling for artifact suppression and a stack-of-spirals k-space trajectory for efficient data acquisition. A semi-automatic 2D method was adapted for 3D image segmentation, and automated 3D methods to calculate strain, twist, and torsion were employed. A scan protocol that covered the majority of the left ventricle in a scan time of less than 25 minutes was developed, and seven healthy C57Bl/6 mice were studied.

RESULTS

Using these methods, multiphase normal and shear strains were measured, as were myocardial twist and torsion. Peak end-systolic values for the normal strains at the mid-ventricular level were 0.29 ± 0.17, -0.13 ± 0.03, and -0.18 ± 0.14 for E(rr), E(cc), and E(ll), respectively. Peak end-systolic values for the shear strains were 0.00 ± 0.08, 0.04 ± 0.12, and 0.03 ± 0.07 for E(rc), E(rl), and E(cl), respectively. The peak end-systolic normalized torsion was 5.6 ± 0.9°.

CONCLUSIONS

Using a 3D cine DENSE sequence tailored for cardiac imaging in mice at 7 T, a comprehensive assessment of 3D myocardial mechanics can be achieved with a scan time of less than 25 minutes and an image analysis time of approximately 1 hour.

Diffusion tensor imaging of the cerebellum and eyeblink conditioning in fetal alcohol spectrum disorder

BACKGROUND

Prenatal alcohol exposure is related to a wide range of neurocognitive effects. Eyeblink conditioning (EBC), which involves temporal pairing of a conditioned with an unconditioned stimulus, has been shown to be a potential biomarker of fetal alcohol exposure. A growing body of evidence suggests that white matter may be a specific target of alcohol teratogenesis, and the neural circuitry underlying EBC is known to involve the cerebellar peduncles. Diffusion tensor imaging (DTI) is a magnetic resonance imaging (MRI) technique that has proven useful for assessing central nervous system white matter integrity. This study used DTI to examine the degree to which the fetal alcohol-related deficit in EBC may be mediated by structural impairment in the cerebellar peduncles.

METHODS

Thirteen children with fetal alcohol spectrum disorder (FASD) and 12 matched controls were scanned using DTI and structural MRI sequences. The DTI data were processed using a voxelwise technique, and the structural data were used for volumetric analyses. Prenatal alcohol exposure group and EBC performance were examined in relation to brain volumes and outputs from the DTI analysis.

RESULTS

Fractional anisotropy (FA) and perpendicular diffusivity group differences between alcohol-exposed and nonexposed children were identified in the left middle cerebellar peduncle. Alcohol exposure correlated with lower FA and greater perpendicular diffusivity in this region, and these correlations remained significant even after controlling for total brain and cerebellar volumes. Conversely, trace conditioning performance was related to higher FA and lower perpendicular diffusivity in the left middle peduncle. The effect of prenatal alcohol exposure on trace conditioning was partially mediated by lower FA in this region.

CONCLUSIONS

This study extends recent findings that have used DTI to reveal microstructural deficits in white matter in children with FASD. This is the first DTI study to demonstrate mediation of a fetal alcohol-related effect on neuropsychological function by deficits in white matter integrity.