Harmonization of Longitudinal Diffusion Tensor Imaging Data of the Pediatric Cervical and Thoracic Spinal Cord at 3T Using Longitudinal ComBat

Diffusion tensor imaging (DTI) of the spinal cord has been extensively used to identify biomarkers for spinal cord pathology. Previously, the longitudinal ComBat (longComBat) technique was examined to reduce scanner effects in multi-site, multi-scanner spinal cord DTI data. This study aimed to assess its effectiveness on longitudinal scans using a single-scanner pediatric dataset, including healthy and spinal cord injury (SCI) subjects. Two identical datasets were collected from 42 healthy and 27 SCI subjects with a 2-hour interval between scans on a 3T Siemens MRI scanner. Axial DTI images of the entire cervical and thoracic spinal cord were obtained, and various average diffusion tensor metrics (FA, MD, RD, & AD) were measured at each vertebral level. Pearson correlation and intraclass correlation coefficients were used to evaluate inter- and intra-subject agreement pre- and post-harmonization. Minimal improvement in agreement was observed with the mean square residual (MSR) model, while the restricted maximum likelihood estimator (REML) model showed reduced intra-subject agreement in all the tensor metrics. The significant variability between longitudinal DTI scans within a single scanner was likely due to physiological motion rather than scanner effects. Post-harmonization using the longComBat MSR model showed limited improvement in agreement.