Spatially selective 2D RF inner field of view (iFOV) diffusion kurtosis imaging (DKI) of the pediatric spinal cord

Correlations of diffusion tensor imaging and clinical measures with spinal cord cross-sectional area measurements in pediatric spinal cord injury patients

PURPOSE

The purpose of this work was to employ a semi-automatic method for measuring spinal cord cross-sectional area (SCCSA) and investigate the correlations between diffusion tensor imaging (DTI) metrics and SCCSA for the cervical and thoracic spinal cord for typically developing pediatric subjects and pediatric subject with spinal cord injury.

METHODS

Ten typically developing (TD) pediatric subjects and ten pediatric subjects with spinal cord injury (SCI) were imaged using a Siemens Verio 3 T MR scanner to acquire DTI and high-resolution anatomic scans covering the cervical and thoracic spinal cord (C1-T12). SCCSA was measured using a semi-automated edge detection algorithm for the entire spinal cord. DTI metrics were obtained from whole cord axial ROIs at each vertebral level. SCCSA measures were compared to DTI metrics by vertebral level throughout the entire cord, and above and below the injury site. Correlation analysis was performed to compare SCCSA, DTI and clinical measures as determined by the International Standards for Neurological Classification of Spinal Cord Injury (ISNCSCI) examination.

RESULTS

In subjects with SCI, FA and SCCSA had a positive correlation (r = 0.81, P < 0.01), while RD and SCCSA had a negative correlation (r = -0.68, P = 0.02) for the full spinal cord. FA and SCCSA were correlated above (r = 0.56, P < 0.01) and below (r = 0.54, P < 0.01) the injury site. TD subjects showed negative correlations between AD and SCCSA (r = -0.73, P = 0.01) and RD and SCCSA (r = -0.79, P < 0.01).

CONCLUSION

The ability to quickly and effectively measure SCCSA in subjects with SCI has the potential to allow for a better understanding of the progression of atrophy following a SCI. Correlations between cord cross section and DTI metrics by vertebral level suggest that imaging inferior and superior to lesion may yield useful information for diagnosis and prognosis.

Diffusion Kurtosis Imaging of Neonatal Spinal Cord in Clinical Routine

Diffusion kurtosis imaging (DKI) has undisputed advantages over the more classical diffusion magnetic resonance imaging (dMRI) as witnessed by the fast-increasing number of clinical applications and software packages widely adopted in brain imaging. However, in the neonatal setting, DKI is still largely underutilized, in particular in spinal cord (SC) imaging, because of its inherently demanding technological requirements. Due to its extreme sensitivity to non-Gaussian diffusion, DKI proves particularly suitable for detecting complex, subtle, fast microstructural changes occurring in this area at this early and critical stage of development, which are not identifiable with only DTI. Given the multiplicity of congenital anomalies of the spinal canal, their crucial effect on later developmental outcome, and the close interconnection between the SC region and the brain above, managing to apply such a method to the neonatal cohort becomes of utmost importance. This study will (i) mention current methodological challenges associated with the application of advanced dMRI methods, like DKI, in early infancy, (ii) illustrate the first semi-automated pipeline built on Spinal Cord Toolbox for handling the DKI data of neonatal SC, from acquisition setting to estimation of diffusion measures, through accurate adjustment of processing algorithms customized for adult SC, and (iii) present results of its application in a pilot clinical case study. With the proposed pipeline, we preliminarily show that DKI is more sensitive than DTI-related measures to alterations caused by brain white matter injuries in the underlying cervical SC.

DTI of chronic spinal cord injury in children without MRI abnormalities (SCIWOMR) and with pathology on MRI and comparison to severity of motor impairment

STUDY DESIGN

This investigation was a cohort study that included: 36 typically developing (TD) children and 19 children with spinal cord lesions who underwent spinal cord MRI.

OBJECTIVES

To investigate diffusion tensor imaging (DTI) cervical and thoracic spinal cord changes in pediatric patients that have clinically traumatic and non-traumatic spinal cord injury (SCI) without MR (SCIWOMR) abnormalities.

SETTING

Thomas Jefferson University, Temple University, Shriners Hospitals for Children all in Philadelphia, USA.

METHODS

36 TD children and 19 children with spinal cord lesions that represent either a chronic traumatic acquired SCI or chronic non-traumatic SCI (≥6 months post injury), age range, 6-16 years who underwent cervical and thoracic spinal cord MRI in 2014-2017. Additionally DTI was correlated to clinical American Spinal Injury Association Impairment Scale (AIS).

RESULTS

Both SCIWOMR and MRI positive (+) groups showed abnormal FA and RD DTI values in the adjacent MRI-normal appearing segments of cephalad and caudal spinal cord compared to TD. The FA values demonstrated perilesional abnormal DTI findings in the middle and proximal segments of the cephalad and caudal cord in the SCIWOMR AIS A/B group compared to SCIWOMR AIS C/D group.

CONCLUSIONS

We found DTI changes in children with SCIWOMR with different causes of spinal lesions. We also investigated the relationship between DTI and clinical AIS scores. This study further examined the potential diagnostic value of DTI and should be translatable to adults with spinal cord lesions.

Advanced Diffusion of the Pediatric Brain and Spine

This article discusses new diffusion-weighted imaging (DWI) sequences, diffusion tensor imaging (DTI), and fiber tractography (FT), as well as more advanced diffusion imaging in pediatric brain and spine. Underlying disorder and pathophysiology causing diffusion abnormalities are discussed. Multishot echo planar imaging (EPI) DWI and non-EPI DWI provide higher spatial resolution with less susceptibility artifact and distortion, which are replacing conventional single-shot EPI DWI. DTI and FT have established clinical significance in pediatric brain and spine. This article discusses advanced diffusion imaging, including diffusion kurtosis imaging, neurite orientation dispersion and density imaging, diffusion spectrum imaging, intravoxel incoherent motion, and oscillating-gradient spin-echo.

Reduced Field-of-View Diffusion-Weighted Magnetic Resonance Imaging of the Pancreas With Tilted Excitation Plane: A Preliminary Study

BACKGROUND

Reduced field-of-view diffusion-weighted imaging (rDWI) with tilted two-dimensional radiofrequency (RF) excitation planes has not yet been applied to the imaging of the pancreas although the utility of this technique which allows the acquisition of high-quality images without aliasing artifacts in the phase-encoding direction has been evaluated for brain and spinal cord imaging.

PURPOSE

To evaluate the visual image quality of the pancreas by tilting the excitation plane (tilted rDWI) in comparison to conventional DWI (cDWI) and rDWI without using the tilted excitation plane.

STUDY TYPE

Retrospective.

POPULATION

Thirty-two patients evaluated for suspected pancreatobiliary diseases.

FIELD STRENGTH/SEQUENCE

Echo-planar imaging DWI (cDWI, rDWI, and tilted rDWI) acquired at 3 T.

ASSESSMENT

Images from each DWI sequence were analyzed by five radiologists to compare image quality (conspicuity of pancreatic edges, interslice signal homogeneity, overall image quality, and conspicuity of focal pancreatic lesions) and artifacts (presence of blurring or ghosting artifacts, susceptibility artifacts, and aliasing artifact).

STATISTICAL TESTS

Shapiro-Wilk test was performed to assess whether data were normally distributed. Friedman test followed by Bonferroni-adjusted Wilcoxon signed-rank test for post hoc analysis was performed to compare image quality and artifact scores.

RESULTS

The mean scores for conspicuity of pancreatic edges (3.36 vs. 2.37), interslice signal homogeneity (3.14 vs. 2.81), presence of ghosting artifacts (3.32 vs. 2.66), susceptibility artifacts (3.06 vs. 2.30), and aliasing artifacts (3.90 vs. 2.34), and overall image quality (3.49 vs. 2.36) were significantly higher in the tilted rDWI than in the rDWI (P < 0.017 for all parameters). The conspicuity score for focal pancreatic lesions tended to be higher in tilted rDWI than in rDWI (2.44 vs. 2.00, P = 0.07).

DATA CONCLUSION

Tilted rDWI had better image quality and reduced artifacts relative to cDWI and rDWI techniques in the pancreas.

LEVEL OF EVIDENCE

4 TECHNICAL EFFICACY: Stage 2.

Multi-task convolutional neural network-based design of radio frequency pulse and the accompanying gradients for magnetic resonance imaging

Modern MRI systems usually load the predesigned RFs and the accompanying gradients during clinical scans, with minimal adaption to the specific requirements of each scan. Here, we describe a neural network-based method for real-time design of excitation RF pulses and the accompanying gradients' waveforms to achieve spatially two-dimensional selectivity. Nine thousand sets of radio frequency (RF) and gradient waveforms with two-dimensional spatial selectivity were generated as the training dataset using the Shinnar-Le Roux (SLR) method. Neural networks were created and trained with five strategies (TS-1 to TS-5). The neural network-designed RF and gradients were compared with their SLR-designed counterparts and underwent Bloch simulation and phantom imaging to investigate their performances in spin manipulations. We demonstrate a convolutional neural network (TS-5) with multi-task learning to yield both the RF pulses and the accompanying two channels of gradient waveforms that comply with the SLR design, and these design results also provide excitation spatial profiles comparable with SLR pulses in both simulation (normalized root mean square error [NRMSE] of 0.0075 ± 0.0038 over the 400 sets of testing data between TS-5 and SLR) and phantom imaging. The output RF and gradient waveforms between the neural network and SLR methods were also compared, and the joint NRMSE, with both RF and the two channels of gradient waveforms considered, was 0.0098 ± 0.0024 between TS-5 and SLR. The RF and gradients were generated on a commercially available workstation, which took ~130 ms for TS-5. In conclusion, we present a convolutional neural network with multi-task learning, trained with SLR transformation pairs, that is capable of simultaneously generating RF and two channels of gradient waveforms, given the desired spatially two-dimensional excitation profiles.

Application of Color Transformation Techniques in Pediatric Spinal Cord MR Images: Typically Developing and Spinal Cord Injury Population

The purpose of this study was to evaluate an improved and reliable visualization method for pediatric spinal cord MR images in healthy subjects and patients with spinal cord injury (SCI). A total of 15 pediatric volunteers (10 healthy subjects and 5 subjects with cervical SCI) with a mean age of 11.41 years (range 8-16 years) were recruited and scanned using a 3.0T Siemens Verio MR scanner. T2-weighted axial images were acquired covering entire cervical spinal cord level C1 to C7. These gray-scale images were then converted to color images by using five different techniques including hue-saturation-value (HSV), rainbow, red-green-blue (RGB), and two enhanced RGB techniques using automated contrast stretching and intensity inhomogeneity correction. Performance of these techniques was scored visually by two neuroradiologists within three selected cervical spinal cord intervertebral disk levels (C2-C3, C4-C5, and C6-C7) and quantified using signal to noise ratio (SNR) and contrast to noise ratio (CNR). Qualitative and quantitative evaluation of the color images shows consistent improvement across all the healthy and SCI subjects over conventional gray-scale T2-weighted gradient echo (GRE) images. Inter-observer reliability test showed moderate to strong intra-class correlation (ICC) coefficients in the proposed techniques (ICC > 0.73). The results suggest that the color images could be used for quantification and enhanced visualization of the spinal cord structures in addition to the conventional gray-scale images. This would immensely help towards improved delineation of the gray/white and CSF structures and further aid towards accurate manual or automatic drawings of region of interests (ROIs).

Diffusion MRI in acute nervous system injury

Diffusion weighted magnetic resonance imaging (DWI) and related techniques such as diffusion tensor imaging (DTI) are uniquely sensitive to the microstructure of the brain and spinal cord. In the acute aftermath of nervous system injury, for example, DWI reveals changes caused by injury that remains invisible on other MRI contrasts such as T₂-weighted imaging. This ability has led to a demonstrated clinical utility in cerebral ischemia. However, despite strong promise in preclinical models and research settings, DWI has not been as readily adopted for other acute injuries such as traumatic spinal cord, brain, or peripheral nerve injury. Furthermore, the precise biophysical mechanisms that underlie DWI and DTI changes are not fully understood. In this report, we review the DWI and DTI changes that occur in acute neurological injury of cerebral ischemia, spinal cord injury, traumatic brain injury, and peripheral nerve injury. Their associations with the underlying biology are examined with an emphasis on the role of acute axon and dendrite beading. Lastly, emerging DWI techniques to overcome the limitations of DTI are discussed as these may offer the needed improvements to translate to clinical settings.

Age related diffusion and tractography changes in typically developing pediatric cervical and thoracic spinal cord

Background and objective

Diffusion tensor imaging (DTI) and diffusion tensor tractography (DTT) are two techniques that can measure white matter integrity of the spinal cord. Recently, DTI indices have been shown to change with age. The purpose of this study is (a) to evaluate the maturational states of the entire pediatric spinal cord using DTI and DTT indices including fractional anisotropy (FA), mean diffusivity (MD), mean length of white matter fiber tracts and tract density and (b) to analyze the DTI and DTT parameters along the entire spinal cord as a function of spinal cord levels and age.

Method

A total of 23 typically developing (TD) pediatric subjects ranging in age from 6 to 16 years old (11.94 ± 3.26 (mean ± standard deviation), 13 females and 10 males) were recruited, and scanned using 3.0 T MR scanner. Reduced FOV diffusion tensor images were acquired axially in the same anatomical location prescribed for the T2-weighted images to cover the entire spinal cord (C1-mid L1 levels). To mitigate motion induced artifacts, diffusion directional images were aligned with the reference image (b0) using a rigid body registration algorithm performed by in-house software developed in Matlab (MathWorks, Natick, Massachusetts). Diffusion tensor maps (FA and MD) and streamline deterministic tractography were then generated from the motion corrected DTI dataset. DTI and DTT parameters were calculated by using ROIs drawn to encapsulate the whole cord along the entire spinal cord by an independent board certified neuroradiologist. These indices then were compared between two age groups (age group A = 6-11 years (n = 11) and age group B = 12-16 years (n = 12)) based on similar standards and age definitions used for reporting spinal cord injury in the pediatric population. Standard least squared linear regression based on a restricted maximum likelihood (REML) method was used to evaluate the relationship between age and DTI and DTT parameters.

Results

An increase in FA (group A = 0.42 ± 0.097, group B = 0.49 ± 0.116), white matter tract density (group A = 368.01 ± 236.88, group B = 440.13 ± 245.24) and mean length of fiber tracts (group A = 48.16 ± 20.48 mm, group B = 60.28 ± 23.87 mm) and a decrease in MD (group A = 1.06 ± 0.23 × 10-3 mm2/s, group B = 0.82 ± 0.24 × 10-3 mm2/s) were observed with age along the entire spinal cord. Statistically significant increases have been shown in FA (p = 0.004, R2 = 0.57), tract density (p = 0.0004, R2 = 0.58), mean length of fiber tracts (p < 0.001, R2 = 0.5) and a significant decrease has been shown in MD (p = 0.002, R2 = 0.59) between group A and group B. Also, it has been shown DTI and DTT parameters vary along the spinal cord as a function of intervertebral disk and mid-vertebral body level.

Conclusion

This study provides an initial understanding of age related changes of DTI values as well as DTT metrics of the spinal cord. The results show significant differences in DTI and DTT parameters which may result from decreasing water content, myelination of fiber tracts, and the thickening diameter of fiber tracts during the maturation process. Consequently, when quantitative DTI and DTT of the spinal cord is undertaken in the pediatric population an age and level matched normative dataset should be used to accurately interpret the quantitative results.

Reduced Field of View Diffusion Tensor Imaging and Fiber Tractography of the Pediatric Cervical and Thoracic Spinal Cord Injury

The aim of this study is to assess the utility and effectiveness of diffusion tensor imaging (DTI) and diffusion tensor tractography (DTT) of the entire pediatric cervical and thoracic spinal cord toward discrimination of typically developing (TD) controls and subjects with spinal cord injury (SCI). A total of 43 pediatric subjects, including 23 TD subjects ranging in age from 6 to 16 years old and 20 subjects with SCI ranging in age from 7 to 16 years, were recruited and scanned using a 3.0 Tesla magnetic resonance scanner. Reduced field of view diffusion tensor images were acquired axially to cover the entire spinal cord across two slabs. For DTI analysis, motion correction was performed by coregistration of the diffusion-weighted images to the reference image (b0). Streamline deterministic tractography results were generated from the preprocessed data. DTI and DTT parameters of the whole cord, including fractional anisotropy (FA), mean diffusivity (MD), tract length, and tract density, were calculated, averaged across the whole spinal cord, and compared between the TD and SCI groups. Statistically significant decreases have been shown in FA (TD = 0.46 ± 0.11; SCI = 0.37 ± 0.09; p < 0.0001) and tract density (TD = 405.93 ± 243.84; SCI = 268.90 ± 270.34; p < 0.0001). However, the mean length of tracts and MD did not show significant differences. When investigating differences in DTI and DTT parameters above and below the injury site, it was shown that the FA and tract density in patients with cervical SCI decreased significantly in the thoracic region. An identical trend was observed in the cervical region for patients with thoracic SCI as well. When comparing TD and SCI subjects, FA and tract density were the most sensitive parameters in detecting functional changes of the spinal cord in chronic pediatric SCI. The results show that both DTI and DTT have the potential to be imaging biomarkers in the diagnosis of SCI.

Identification of ghost artifact using texture analysis in pediatric spinal cord diffusion tensor images

PURPOSE

Ghost artifacts are a major contributor to degradation of spinal cord diffusion tensor images. A multi-stage post-processing pipeline was designed, implemented and validated to automatically remove ghost artifacts arising from reduced field of view diffusion tensor imaging (DTI) of the pediatric spinal cord.

METHOD

A total of 12 pediatric subjects including 7 healthy subjects (mean age=11.34years) with no evidence of spinal cord injury or pathology and 5 patients (mean age=10.96years) with cervical spinal cord injury were studied. Ghost/true cords, labeled as region of interests (ROIs), in non-diffusion weighted b0 images were segmented automatically using mathematical morphological processing. Initially, 21 texture features were extracted from each segmented ROI including 5 first-order features based on the histogram of the image (mean, variance, skewness, kurtosis and entropy) and 16s-order feature vector elements, incorporating four statistical measures (contrast, correlation, homogeneity and energy) calculated from co-occurrence matrices in directions of 0°, 45°, 90° and 135°. Next, ten features with a high value of mutual information (MI) relative to the pre-defined target class and within the features were selected as final features which were input to a trained classifier (adaptive neuro-fuzzy interface system) to separate the true cord from the ghost cord.

RESULTS

The implemented pipeline was successfully able to separate the ghost artifacts from true cord structures. The results obtained from the classifier showed a sensitivity of 91%, specificity of 79%, and accuracy of 84% in separating the true cord from ghost artifacts.

CONCLUSION

The results show that the proposed method is promising for the automatic detection of ghost cords present in DTI images of the spinal cord. This step is crucial towards development of accurate, automatic DTI spinal cord post processing pipelines.

Spinal cord microstructure integrating phase-sensitive inversion recovery and diffusional kurtosis imaging

PURPOSE

The aim of this prospective study was to determine the feasibility in terms of repeatability and reproducibility of diffusional kurtosis imaging (DKI) for microstructural assessment of the normal cervical spinal cord (cSC) using a phase-sensitive inversion recovery (PSIR) sequence as the anatomical reference for accurately defining white-matter (WM) and gray-matter (GM) regions of interests (ROIs).

METHODS

Thirteen young healthy subjects were enrolled to undergo DKI and PSIR sequences in the cSC. The repeatability and reproducibility of kurtosis metrics and fractional anisotropy (FA) were calculated in GM, WM, and cerebral-spinal-fluid (CSF) ROIs drawn by two independent readers on PSIR images of three different levels (C1-C4). The presence of statistically significant differences in DKI metrics for levels, ROIs (GM, WM, and CSF) repeatability, reproducibility, and inter-reader agreement was evaluated.

RESULTS

Intra-class correlation coefficients between the two readers ranged from good to excellent (0.75 to 0.90). The inferior level consistently had the highest concordance. The lower values of scan-rescan variability for all DKI parameters were found for the inferior level. Statistically significant differences in kurtosis values were not found in the lateral white-matter bundles of the spinal cord.

CONCLUSION

The integration of DKI and PSIR sequences in a clinical MR acquisition to explore the regional microstructure of the cSC in healthy subjects is feasible, and the results obtainable are reproducible. Further investigation will be required to verify the possibility to translate this method to a clinical setting to study patients with SC involvement especially in the absence of MRI abnormalities on standard sequences.

Fast imaging of mean, axial and radial diffusion kurtosis

Diffusion kurtosis imaging (DKI) is being increasingly reported to provide sensitive biomarkers of subtle changes in tissue microstructure. However, DKI also imposes larger data requirements than diffusion tensor imaging (DTI), hence, the widespread adaptation and exploration of DKI would benefit from more efficient acquisition and computational methods. To meet this demand, we recently developed a method capable of estimating mean kurtosis with only 13 diffusion weighted images. This approach was later shown to provide very accurate mean kurtosis estimates and to be more efficient in terms of contrast to noise per unit time. However, insofar, the computation of two other critical DKI parameters, radial and axial kurtosis, has required the estimation of all 22 variables parameterizing the full DKI signal expression. Here, we present two strategies for estimating all of DKI's principal parameters - mean kurtosis, radial kurtosis, and axial kurtosis - using only 19 diffusion weighted images, compared to the current state-of-the-art acquisitions typically requiring about 60 images. The first approach is based on axially symmetric diffusion and kurtosis tensors, presented here for the first time, and referred to as axially symmetric DKI. The second approach is applicable in tissues with a priori known principal diffusion direction, and does not require fitting of any kind. The approaches are evaluated in human brain in vivo as well as in fixed rat spinal cord, and are demonstrated to provide metrics in good agreement with their full DKI counterparts estimated with nonlinear least squares. For small data sets and in white matter, axially symmetric DKI provides more accurate and robust estimates than unconstrained DKI. The significant acceleration achieved further paves the way to routine application of the technique.