Abstract
BACKGROUND AND PURPOSE: Multi-parametric MRI, provides a variety of biomarkers sensitive to white matter integrity, However, spinal cord MRI data in pediatrics is rare compared to adults. The purpose of this work was 3-fold: 1) to develop a processing pipeline for atlas-based generation of the typically developing pediatric spinal cord WM tracts, 2) to derive atlas-based normative values of the DTI indices for various WM pathways, and 3) to investigate age-related changes in the obtained normative DTI indices along the extracted tracts.
MATERIALS AND METHODS: DTI scans of 30 typically developing subjects (age range, 6–16 years) were acquired on a 3T MR imaging scanner. The data were registered to the PAM50 template in the Spinal Cord Toolbox. Next, the DTI indices for various WM regions were extracted at a single section centered at the C3 vertebral body in all the 30 subjects. Finally, an ANOVA test was performed to examine the effects of the following: 1) laterality, 2) functionality, and 3) age, with DTI-derived indices in 34 extracted WM regions.
RESULTS: A postprocessing pipeline was developed and validated to delineate pediatric spinal cord WM tracts. The results of ANOVA on fractional anisotropy values showed no effect for laterality (P = .72) but an effect for functionality (P < .001) when comparing the 30 primary WM labels. There was a significant (P < .05) effect of age and maturity of the left spinothalamic tract on mean diffusivity, radial diffusivity, and axial diffusivity values.
CONCLUSIONS: The proposed automated pipeline in this study incorporates unique postprocessing steps followed by template registration and quantification of DTI metrics using atlas-based regions. This method eliminates the need for manual ROI analysis of WM tracts and, therefore, increases the accuracy and speed of the measurements.
ABBREVIATIONS:
- AD
- axial diffusivity
- CV
- coefficient of variation
- DF
- dorsal funiculi
- FA
- fractional anisotropy
- LF
- lateral funiculi
- MD
- mean diffusivity
- RD
- radial diffusivity
- VF
- ventral funiculi
In recent years, DTI has shown promise as a noninvasive imaging biomarker for evaluating the tissue microstructure,1 which is sensitive to the integrity of the spinal cord WM. Recently, numerous studies have demonstrated the utility of DTI as a noninvasive tool to assess adult spinal cord WM microstructure.2,3 However, studies related to the pediatric spinal cord have been limited so far.4,5 This issue highlights the importance of studying the pediatric spinal cord at a microstructural level and establishing normative values for various diffusion-derived indices obtained from the WM structures. A limited number of studies in pediatric subjects have demonstrated DTI indices differences between age and cord regions.4,5 However, all of these studies, so far, have used manual delineation of the WM structures in the pediatric spinal cord.5 Recently, methods have been developed that use automated extraction of the WM structures in the adult spinal cord.6 This approach alleviates any subjectivity in identification of WM structures. In this work, we have extended these methods to the pediatric spinal cord.
In summary, the aim of this work was 3-fold: 1) to develop a processing pipeline for atlas-based generation of pediatric spinal cord WM tracts. This provides a faster and more accurate method of tract delineation and eliminates the need for manual WM tract extraction. 2) to derive atlas-based normative values of the DTI parameters for various WM pathways at the C3 level of a typically developing pediatric spinal cord, and 3) to examine age-related changes in the derived normative DTI parameters along the extracted tracts.
MATERIALS AND METHODS
Study Design and Population
All subjects provided written informed consent obtained under a protocol approved by the institutional review board at Thomas Jefferson University. The study used 30 typically developing pediatric subjects (age range, 6–16 years; mean age, 12.38 [SD, 2.81] years). Following the recruitment, all subjects went through a brief assessment, and those with any abnormality of the nervous system were excluded.
MR Imaging Protocol
The scans were obtained using a 3T Verio MR imaging scanner (Siemens) with a 4-channel neck and an 8-channel spine matrix coil. All the subjects underwent the scan protocol, consisting of an initial T2-weighted gradient-echo scout sagittal acquisition of the cervical and thoracic spinal cord. The sagittal scans were used to prescribe axial sections of both the cervical and thoracic spine. Imaging was then followed by obtaining an axial T2-weighted gradient recalled-echo sequence, a sagittal T2-weighted 3D sampling perfection with application-optimized contrasts by using different flip angle evolution (3D-SPACE; Siemens) sequence, a sagittal TSE T1-weighted sequence, a sagittal TSE T2-weighted sequence, and axial diffusion-weighted scans. We used the 3D-SPACE and diffusion-weighted data in this study. Diffusion-weighted images were acquired axially using 2 overlapping slabs to cover the entire cervical and thoracic spinal cord using an inner FOV spin-echo-based EPI pulse sequence with a tilted excitation plane. The inner FOV diffusion-weighted sequence consists of a spatially selective 2D radiofrequency excitation profile.7 This sequence has the advantage of reducing distortion artifacts in diffusion-weighted scans by enabling faster k-space encoding via a reduction of phase-encoding steps. Manual shim and fat-saturation volume adjustments were also performed.
The imaging parameters of the diffusion-weighted acquisition for each slab are as follows: FOV = 164 mm, phase FOV = 28.4% (47 mm), 3 averages of 20 diffusion directions, 6 b = 0 acquisitions, b = 800 s/mm,2 voxel size = 0.8 × 0.8 × 6 mm3, number of slices = 40, TR = 7900 ms, TE = 110 ms, acquisition time = 8 minutes and 49 seconds per slab. The T2-weighted SPACE imaging parameters used were: voxel size = 1 × 1 × 1 mm3, TR = 1500 ms, TE = 122 ms, flip angle = 140°, number of averages = 2, and acquisition time = 3 minutes and 17 seconds (for each overlapping slab).
Image Preprocessing
The diffusion-weighted data from all the individual subjects were motion-corrected and later registered to the template (PAM50) in the Spinal Cord Toolbox (SCT; https://spinalcordtoolbox.com/en/latest/) for performing atlas-based analysis of WM tracts.6 Registration is best done using the high-resolution anatomic image. In this study, all the pre- and postprocessing of the data was performed entirely within the framework of the SCT.8 These processes are described in detail below and illustrated in Fig 1.
Procedure for Motion Correction
All diffusion-weighted data were motion-corrected using a slice-wise motion-correction method.8 Specifically, using the sct_dmri_moco module from SCT, slice-by-slice planar translations within the diffusion-weighted data were estimated, while the regularization was constrained along the z-direction (SliceReg algorithm).
Procedure for Data Registration to the Template
After applying motion correction on the diffusion-weighted data, we registered the data with the PAM50 6 template available in the SCT. The registration procedure was done following a 2-step multimodal registration process within the SCT. The first step involved registration between the PAM50 template and an anatomical image of all the individual subjects. This was followed by a second registration between the template (in anatomical space) and the diffusion-weighted data of the same subjects. At the end of this process, the PAM50 atlas, which contains the various WM tracts, was registered to our own subject-specific diffusion-weighted data. The following steps describe the registration process in detail:
Segmentation of spinal cord: The structural image and diffusion-weighted data of the cervical spine of all subjects were initially segmented using the sct_propseg9module. The goal of this step was to extract a binary mask of the segmented spinal cord. Manual corrections were made on the images if the segmentation at certain locations was not satisfactory.
Labeling the anatomical image: In this step, the segmented structural images obtained from the previous step were labeled using the sct_label_vertebrae10 module. This step generated segmented cord labeled with vertebral levels.
Identifying vertebral levels along the cervical cord: Two vertebral levels (C3 and T1) were identified on the anatomical image using the sct_label_utils module. The main purpose of this step was to enable more accurate registration between the pediatric anatomical image and the PAM50 template by providing 2 reference voxels, one at the C3 level and a second one at the T1 level (orange arrows pointing to the corresponding intervertebral discs shown in Fig 1).
Next, each subject’s individual T2-weighted image was registered to the PAM50 template using the sct_register_to_template module. This process includes straightening of the spinal cord and application of local deformations using a multistep nonlinear registration constrained in the axial plane and regularized in the z-direction. A T2 template in the anatomic space (template2anat) and the corresponding warping field (warp_template2anat) were generated as the output of this step.
A second registration was then performed to register (template2anat); the image from the previous step to each subject’s diffusion data using the sct_register_multimodal module. A template image in each individual’s diffusion native space (template2diff) as well as the corresponding warping field (warp_anat2diff) were generated as the output of this step.
Concatenation of deformation fields: The warping fields from the last 2 steps of registrations were concatenated to create a single global transformation between the template and diffusion data (warp_template2diff).
Obtaining the WM atlas: Finally, the global warping field was used to warp the PAM50 template objects (ie, WM atlas) to the diffusion data.
Quantifying DTI Metrics Using the Atlas
After successful registration was verified visually for each individual subject, DTI metrics were extracted (sct_extract_metric) for each subject at the C3 level using a weighted-average method.11 This step resulted in quantifying DTI metrics of fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) for 34 extracted WM labels from the PAM50 template. Of the 34 labels, 30 represent individual sensory and motor WM tracts, 3 represent combined labels of the dorsal funiculi (DF), lateral funiculi (LF), and ventral funiculi (VF), and 1 label was assigned to represent the entire WM itself. The normative values were obtained by averaging all these measurements across all 30 subjects. The full names of labels are provided in the Online Supplemental Data.
Statistical Analysis
GraphPad Prism software 8.0 (GraphPad Software), a DOS-based scientific graphing application, was used for statistical analysis of the data. For ease of visualization, DTI normative measurements of all the 34 labels were plotted using violin plots. One-way ANOVA was used to study the laterality (left/right), functionality (motor/sensory), and age-related changes among all the tracts. Linear regression analysis was used to compare age-related changes of DF, LF, and VF labels using DTI metrics. For the age-related analysis, the subject population was divided into 2 age groups: one group of younger children from 7–11 years of age and the older group from 12–16 years of age. The age groups were defined on the basis of a previous pediatric study.5 Throughout the study, a P value of .05 was used to determine the statistical significance.
RESULTS
Delineated WM Tracts
Figure 2A depicts a WM atlas of the human spinal cord. Figure 2B illustrates a subset of tracts overlaid on a b=0 image at the C3 level in a representative subject from this study. The labeled WM tracts in this figure are explained in the figure legend. Figure 2C shows a sagittal T2-weighted cross-section of the spinal cord marked at the C3 level.
Normative Values
The developed automated WM pipeline was successful in extracting DTI values from all the various tracts across all 30 subjects (Fig 3). The violin plots were subdivided into the dorsal, lateral, and ventral funiculi as well as combined tracts (WM, DF, LF, and VF). The FA in the dorsal column tracts ranged from 0.61 to 0.62 with the mean 0.62 (SD, 0.09). The corresponding values for lateral and ventral columns were 0.47 (SD, 0.12) and 0.28 (SD, 0.12), respectively.
FA, MD, and RD showed a higher level of dispersion in all the tracts, including dorsal, lateral, and ventral tracts compared with AD. The coefficient of variation (CV) of the measurements in FA ranged from 14.62% in the dorsal column to 45.74% in ventral tracts. Similar measurements in MD were observed (23.75% to 34.06% from the dorsal-to-ventral tracts and in RD from 33.43% to 40.90%, respectively). AD showed less variability across all the tracts, with a CV of 22% in the dorsal columns to 26% in the ventral tracts. As observed, ventral funiculi regions present larger variations compared with other WM regions. WM measurements from “combined tracts” had the lowest variation in FA, with a CV of 18.6%, followed by the second lowest in AD (CV = 20.11%). The VF had higher variation in all the metrics, with the average CV of 33% compared with its counterparts DF and LF (CV(DF) = 22%, CV(LF) = 25%).
The Table compares the diffusion indices obtained for 6 selected tracts shown in Fig 2B. From all the tracts in the Table, the left fasciculus gracilis shows the highest FA of 0.66 averaged among all 30 subjects with the lowest variation (CV = 14.02%). FA was the lowest in the right lateral corticospinal tract with a higher CV of 18.26%.
Effects of Laterality, Functionality, and Age
Results of ANOVA using FA, MD, RD, and AD showed no significant difference (P > .05) between a total of 15 left and 15 right tracts. However, when tested for functionality (i.e., whether a tract is sensory or motor), FA and MD showed a significant difference (PFA = 3.7E−5, PMD = .043), but AD and RD showed no effect.
Examining the effect of age with each single individual tract as well as combined tracts revealed significant differences among the maturation of some tracts and age. Among all 34 selected tracts studied, only the left spinal lemniscus (spinothalamic and spinoreticular tracts) showed an effect with age using MD, RD, and AD (P < .05). Additionally, linear regression analysis showed a positive correlation between FA and age as shown in Fig 4A, while MD, RD, and AD (Fig 4B–D) showed a negative correlation in all 3 combined tracts (dorsal, ventral, and lateral funiculi). MD showed significant differences with age in the VF and DF (P < .05), whereas AD showed an effect in the LF and VF but not the DF. Similar analysis showed no effect for FA but a significant effect for RD in all 3 combined tracts of DF, LF, and VF (P < .05).
DISCUSSION
In this study, we developed a processing pipeline for automatic generation of pediatric spinal cord WM tracts. The acquisition was followed by processing the data with a unique pipeline and registration to an external template. The quantification of DTI indices was facilitated with the atlas incorporated in PAM50,6 followed by establishing DTI normative values for 34 WM pathways. Former studies have addressed the generation of the spinal cord atlas in the adult population.8,11 However, no atlas of WM pathways in the pediatric spinal cord is currently available. This study is the first to propose a pipeline for automated WM tract extraction in pediatrics enabling atlas-based analysis of the spinal cord.
WM Tract-generation Pipeline
The WM atlas of the spinal cord based on DTI enables automated atlas-based analysis of WM pathways. The extracted tracts in this study are aligned with the PAM50 T2-weighted template, which minimizes the interpolation errors during diffusion metrics extraction.
In the SCT, a digitized version of the Gray’s Anatomy spinal cord WM atlas is merged with PAM50. Currently, this atlas is the most detailed atlas of human anatomy.11,12 However, one of the limitations of this atlas is that it does not take into account the between-subject variability of the spinal cord shape, which has been addressed in previous studies.13 We incorporated a multistep registration method to register the pediatric diffusion-weighted data to the adult PAM50 template. To make this process more feasible, we used the sct_label_utils tool from the SCT to create 2 landmarks at the C3 and T1 vertebral levels in the anatomical image of each subject. This method considers the height difference between the adult and pediatric population and adjusts for the discrepancies. We then used a BSplineSyN14 algorithm in the final step of the multimodal registration to preserve the shape of the cord. We visually examined the outcome of the registration in each individual subject, and the results were satisfactory. The suggested method provides a fully automated and accurate WM segmentation compared with whole-cord manual ROI drawings used in similar spinal cord DTI studies such as the work of Reynolds et al.4 The manual ROI drawing performed in the studies of Reynolds et al4 and Saksena et al5 is also prone to partial volume contamination, which can be alleviated with our proposed automated pipeline.15
DTI Metrics Quantification of the Spinal Cord
In this work, we obtained normative values of diffusion data at the C3 vertebral level for ventral, dorsal, and lateral spinal cord funiculi as well as a specific subset of WM tracts in the healthy pediatric population. We believe that these values could be useful when comparing healthy and patient populations.16
Considering the various known spinal cord complications, study of spinal cord WM and characterizing its microstructure noninvasively will help in accurate understanding and diagnosis of disease and pathologies.
Injury to the spinal cord due to traumatic or nontraumatic events and subsequently to white matter tracts in general leads to specific dysfunction and various clinical disabilities. Therefore, identifying the various specialized types of the WM pathways in the cord that may have been damaged is crucial in helping the prognosis of functional recovery and disease management.
Note that our analysis focused only on a single spinal cord level (C3 vertebral level) due to the better performance of proposed segmentation and registration algorithms with lower interrater variability reported by Prados et al.17 Also, as shown by Losseff et al,18 C2/C3 level results are the most reproducible values when studying spinal cord atrophy. In addition, the work of Stevenson et al19 suggests that the measurement variability that occurs using other spinal levels could be eliminated by focusing the analysis on the C2/C3 intervertebral disc. This possibility is due to the excellent contrast of the cord and CSF at the C3 level (i.e., a capacious CSF pool).
When we investigated the function of age and DTI indices, our observed age-related trend of DTI indices matches that of other spinal cord studies.5 However, the lack of a statistical significance effect with a decrease in AD observed in the dorsal columns suggests that the proportion of larger-sized axons increases with the child’s growth.4
In comparing the DTI metrics of the funiculi with age, we found no significant difference between FA values of the LF, DF, and VF. However, MD, RD, and AD values of the VF showed a significant effect with age-related changes. These differences could reflect the variability of the axon diameter among these funicles20 containing different kinds of tracts. Despite the lack of a significant effect observed between the FA and age, the correlation of age with FA specifically and other DTI-derived indices generally was consistent with that in former studies.5,21
While examining the effect of age on each individual tract in this study, we observed that there is a relationship of diffusivity metrics and maturation of some tracts (i.e., the spinothalamic tract). This could be because of the myelination process in different tracts. Myelination constitutes the last stage in the development of WM, which begins in utero in the second trimester of pregnancy and continues up to the end of adolescence.22 There are several factors that affect the rate of myelination, including the type of fiber, where faster rates are observed in sensory compared with motor tracts.23 We believe that this could justify the significant effect of spinothalamic tract (a somatosensory tract) maturation in DTI indices in our study.
We also showed that there is no significant difference in the laterality (left/right) of the tracts. These findings suggest structural similarity between the left- and right-sided tracts. In a study by Duval et al,24 the authors have shown similar average axon diameters in identical tracts from the left and right sides of the spinal cord. Although there was no significant difference in the laterality of the tracts, ANOVA showed a significant difference between sensory (i.e., fasciculus gracilis, fasciculus cuneatus) and motor (i.e., corticospinal tract) pathways. We believe that this finding is related to the microstructural organization of the sensory and motor WM tracts. The anatomic variation of axon diameter is noticeable; whereas the dorsal column fasciculus gracilis contains mainly small-diameter axons, the corticospinal tract in the lateral and ventral columns contains a combination of large- and small-diameter axons.25 Another anatomic variation to consider is the fiber density variation along the entire length of the spinal cord, especially at the exit points of the brachial and lumbosacral plexuses.12 In a study by Ohnishi et al,26 the authors reported a decrease in the number of myelinated fibers of the fasciculus gracilis in humans between C3 and T5. Terao et al27 also reported an association of axonal diameter variation at the vertebral level in the human lateral corticospinal tract. These observations are in accordance with previous histology work28 in humans, emphasizing the importance of broader studies of histologic validation in the ex vivo human spinal cord.
Limitations
Due to the diffusion data, SNR limitations, lack of anatomical landmarks, and the smaller size of the cord in children, there could be higher intrasubject variability resulting from spatial misalignment during registration of diffusion data to structural scans. This limitation emphasizes the need for tools that provide robust preprocessing, normalization, and motion correction for spinal cord diffusion data. Another limitation stems from the intrinsic variability of spinal cord diffusion imaging, which induces inaccuracies in metrics quantification. This issue, along with the selected registration algorithms method of metric extraction and partial volume effects, contributes to intersubject variability, particularly with a smaller cord size. An additional limitation is the small number of subjects as well as the narrow age range of our population as opposed to other pediatric diffusion studies.4 All of the aforementioned limitations could be contributing factors to observed discrepancies in our results compared with other similar studies.
CONCLUSIONS
We demonstrated the feasibility of a fully automated WM tract extraction in the pediatric spinal cord. We also established normative values for DTI-derived indices, looking at specific WM regions of the cord. These results show that the presented values are microstructure-sensitive and evolve in the pediatric spinal cord during development. This finding warrants future work with spinal cord patients (i.e., spinal cord injury) as well as an examination of DTI-derived metrics in the entire cord compared with a single C3 section.
Footnotes
This work was supported by the National Institute of Neurological Disorders of the National Institutes of Health under award Nos. R01NS079635 and R01NS111113.
Disclosures: Laura Krisa—RELATED: Grant: National Institutes of Health.* *Money paid to the institution.
Indicates open access to non-subscribers at www.ajnr.org
References
- Received February 9, 2021.
- Accepted after revision April 19, 2021.
- © 2021 by American Journal of Neuroradiology