Abstract
BACKGROUND AND PURPOSE: Histopathologic studies have demonstrated WM damage in primary Sjögren syndrome. The purpose of this study was to evaluate WM microstructural changes by use of DTI-derived parameters in patients with primary Sjögren syndrome.
MATERIALS AND METHODS: DTI was performed in 19 patients with primary Sjögren syndrome (age, 64.73 ± 9.1 years; disease duration, 11.5 ± 7.56 years) and 16 age-matched control subjects. Exclusion criteria were a history of major metabolic, neurologic, or psychiatric disorder and high risk for cardiovascular disease. Data were analyzed by use of tract-based spatial statistics, for which the WM skeleton was created, and a permutation-based inference with 5000 permutations was used with a threshold of P < .01, corrected for multiple comparisons to enable identification of abnormalities in fractional anisotropy, mean diffusivity, radial diffusivity, and axial diffusivity.
RESULTS: Tract-based spatial statistics showed decreased fractional anisotropy in multiple areas in patients with primary Sjögren syndrome compared with control subjects, located mainly in the corticospinal tract, superior longitudinal fasciculus, anterior thalamic radiation, inferior fronto-occipital fasciculus, uncinate fasciculus, and inferior longitudinal fasciculus. Increased mean diffusivity and radial diffusivity and decreased axial diffusivity were observed in most of the fiber tracts of the brain in patients with primary Sjögren syndrome, compared with control subjects.
CONCLUSIONS: Patients with primary Sjögren syndrome show loss of WM microstructural integrity, probably related to both Wallerian degeneration and demyelination.
ABBREVIATIONS:
- AD
- axial diffusivity
- FA
- fractional anisotropy
- MD
- mean diffusivity
- pSS
- primary Sjögren syndrome
- RD
- radial diffusivity
- TBSS
- tract-based spatial statistics
Sjögren syndrome is a chronic systemic autoimmune disease that can be classified as primary Sjögren syndrome (pSS) when presenting in isolation or secondary when related to another connective tissue disease.1 The prevalence of pSS reported in different studies ranges from 0.1–4.8%.2 pSS is characterized by mononuclear infiltration and destruction of the exocrine glands, mainly the lachrymal and salivary glands, but extraglandular manifestations are also reported (eg, arthralgia, pulmonary involvement, renal tubular acidosis, etc).1,3 Involvement of both the peripheral and the CNS has also been reported in pSS.4⇓–6 Although involvement of the peripheral nervous system is a well-documented feature of the disease, the prevalence, the type, and the underlying mechanism of CNS involvement remain unclear.4,5,7 The estimated frequency of CNS involvement ranges from 10–60% in different reports, depending on the parameters studied (eg, patient selection, diagnostic criteria, etc).4⇓–6 Patients with pSS can present with a wide range of focal or diffuse neurologic or psychiatric manifestations, including motor/sensory deficits, transverse myelitis, and cognitive impairment.4⇓–6 The current data from MRI studies support an increased frequency of high signal intensity lesions in the periventricular and/or subcortical WM on FLAIR and T2-weighted imaging, observed mainly in patients with pSS and evidence of CNS disease.8⇓⇓–11 The volumetric analysis of GM and WM by use of the voxel-based morphometry method demonstrated diffuse atrophy in patients with pSS.8 SPECT and PET studies have demonstrated reduced CBF and lowered glucose metabolism in patients with pSS.12,13
DTI is a technique that allows assessment of the preferential direction of the Brownian motion of protons, which in the brain reflects the microscopic architecture of the WM.14,15 Four quantitative diffusion parameters can be derived from DTI data: 1) fractional anisotropy (FA), reflecting the directionality of water diffusion and coherence of WM fiber tracts; 2) mean diffusivity (MD), quantifying the overall magnitude of water diffusion; 3) axial diffusivity (AD), measuring the magnitude of diffusivity along the principal diffusion direction; and 4) radial diffusivity (RD), reflecting the magnitude of diffusivity perpendicular to the principal diffusion direction.14,15 These metrics have been correlated with the microstructural organization of WM and are used to infer structural characteristics of the local tissue.14,15 The DTI-derived parameters have been used to investigate WM microstructure in various disorders by use of an ROI approach, in which structures of interest are manually defined on MR images, but information about changes in brain diffusivity in pSS is scarce.16 Using the ROI method, a decrease in FA and an increase in MD values have been reported in patients with pSS.16 Similar analytic methods that are based on manually selected ROIs have limitations, mainly because they do not examine the whole brain and they are laborious and time-consuming and therefore prone to human error.17 To eliminate the limitations of the ROI-based methods, the so-called voxelwise analysis methods have been developed, which examine the whole brain automatically at a voxel level.17 The technique of tract-based spatial statistics (TBSS) allows voxelwise statistical analysis of DTI-derived data.18 TBSS has been widely used for DTI analysis because of its advanced registration capabilities and its robust nonparametric assessment of local differences in WM integrity between groups.18
The purpose of the present study was to assess the presence and location of WM damage and to elucidate the basis of WM microstructural changes in patients with pSS by analyzing DTI-derived parameters with the automated TBSS method.
Materials and Methods
Study Patients
The study population consisted of 19 consecutively registered, unselected patients with pSS being followed up in the outpatient rheumatology clinic of our hospital, ages 47–78 years (mean ± SD; 64.73 ± 9.1 years), with a disease duration of 5–28 years (mean ± SD; 11.5 ± 7.56 years). The diagnosis of pSS was established according to the American-European Consensus Criteria. Association with other connective tissue diseases was ruled out, and only patients with pSS were included. The demographic and clinical characteristics of the patients are shown in Table 1. The control group consisted of 16 age-matched healthy volunteers, ages 45–76 years (mean ± SD; 62.57 ± 8.3 years). The protocol and the procedure were explained in detail to all patients and control subjects who had the same educational background. The study was performed with the approval of the institutional review board, and all the participants signed a written informed consent agreement.
Demographic and clinical characteristics of 19 patients with pSS
Exclusion criteria were a history or clinical signs of cardiovascular disease, peripheral arterial disease, hepatic dysfunction (serum transaminase levels >1.5 times the upper limit of normal), renal insufficiency (serum creatinine concentration >1.6 mg/dL), proteinuria (> 0.5 g/d), diabetes mellitus (fasting plasma glucose concentration ≥126 mg/dL or use of antidiabetic medication), hypertension (arterial blood pressure >140/90 mm Hg or use of antihypertensive medication), serum level of thyroid-stimulating hormone >5 mU/mL, and treatment with corticosteroids during the previous 6 months. None of the study patients or control subjects had findings suggestive of CNS or psychiatric disorder. A total of 30 patients with pSS were initially evaluated. However, 11 were excluded. More specifically, 3 had hypertension for many years, 2 had atrial fibrillation, 4 had diabetes mellitus and dyslipidemia, and 2 had claustrophobia. Thus, 19 patients were finally included. Routine neurologic examination was performed on all subjects and did not reveal deficits indicative of central or peripheral nervous system involvement. Similarly, patients and caregivers did not report any symptoms indicative of cognitive decline, depression, or fatigue.
Data Acquisition
DTI was performed by use of a 1.5T scanner (Gyroscan ACS NT; Philips Healthcare, Best, The Netherlands). We used a single-shot EPI sequence. Parameters for DTI acquisition were as follows: FOV = 230 × 230 mm, 112 × 128 matrix, section thickness of 3 mm, TE = 131 ms, TR = 9825 ms, number of sections = 42, section gap = 0 mm. We used 16 noncollinear gradient directions, with maximum b = 700 seconds/mm2 and scanning time 4 minutes, 34 seconds. The imaging protocol also included 1) a FLAIR sequence (TR = 6300 ms, TE = 120 ms, FOV = 250 mm, matrix = 256 × 256, section thickness of 6 mm, intersection gap = 0.6, scanning time = 2 minutes, 50 seconds), 2) a T1-weighted, high-resolution (1 × 1 × 1 mm), 3D spoiled gradient-echo sequence (TR = 25 ms, TE = 4.6 ms, acquisition matrix = 256 × 228, FOV = 220 mm, scanning time = 5 minutes, 43 seconds).
The presence of areas of high signal intensity in the WM was assessed for each subject on FLAIR images by 2 neuroradiologists who were blinded to the patient/control status.
Data Preprocessing
Images were processed by use of the FSL (FMRIB Software Library; http://www.fmrib.ox.ac.uk/fsl) software package.19 For each subject, all images including diffusion-weighted and B0 images were corrected for eddy current–induced distortion and subject motion effect by use of the FSL Diffusion Toolbox. Brain mask was created from the first B0 image by use of the FSL Brain Extraction Tool, and Diffusion Toolbox was used to fit the tensor model and to compute the FA, MD, AD, and RD maps.
Tract-Based Spatial Statistics Analysis
Voxelwise analysis was performed by use of TBSS.18 First, every FA image was aligned to every other one; by use of all these comparisons, the software then identified the “most representative” one and used it as the target image. This target image was then affine-aligned into Montreal Neurological Institute 152 standard space. The FA data of all subjects were aligned to this target image by use of the nonlinear registration tool FNIRT, which uses a b-spline representation of the registration warp field. Next, the mean FA image was created and thinned to create a mean FA skeleton, which represents the centers of all tracts common to the group. A threshold of FA >0.2 was applied to the skeleton to include only major fiber bundles. The aligned FA data of each subject were then projected onto this skeleton. By applying the original nonlinear registration of each subject's FA to standard space, the MD, RD, and AD maps were also projected onto the mean FA skeleton. The projected maps were separated into 2 groups (control subjects and patients with pSS) and were used to calculate voxelwise cross-subject diffusion statistics. The localization of all the anatomic information described was based on the Johns Hopkins University WM tractography atlas and the International Consortium for Brain Mapping DTI-81 WM labels (part of the FSL package).
Statistical Analysis
To determine FA, MD, RD, and AD differences between groups, the voxelwise analysis was performed by a permutation-based inference method as implemented in the Randomize FSL tool.20 We used a t test, 5000 permutations, and threshold-free cluster enhancement with a threshold of P < .01, corrected for multiple comparisons by use of family-wise error correction to assess differences in the FA, MD, RD, and AD between the patients with pSS and the healthy control subjects.21 In addition, for the patients with pSS, voxelwise associations between each diffusion metric (FA, MD, AD, and RD) and disease duration were performed by use of a nonparametric, linear regression model. The significance threshold for correlations was set at P < .05, corrected for multiple comparisons by use of family-wise error correction (threshold-free cluster enhancement option in the Randomize permutation-testing tool).
Results
There was no significant difference in age between patients and control subjects. Areas of high signal intensity in WM were observed in 13 of the 19 patients (68.4%) and in 6 of the 16 control subjects (37.5%).
Tract-Based Spatial Statistics
In the voxelwise-based group comparison, multiple WM areas with significant FA decrease (P < .01, family-wise error corrections for multiple comparison) were found bilaterally in patients with pSS compared with control subjects, as shown in Fig 1A, in the corticospinal tract, the superior longitudinal fasciculus, the anterior thalamic radiation, the inferior fronto-occipital fasciculus, the uncinate fasciculus, and the inferior longitudinal fasciculus (Table 2).
Tract-based spatial statistics results demonstrate voxelwise comparisons between 19 patients with pSS and 16 control subjects. Statistical maps (thresholded at threshold-free cluster enhancement, P < .01) are overlaid onto the mean FA skeleton and the Montreal Neurological Institute 152 template. Decreased fractional anisotropy in patients with pSS is shown in red (A), increased mean diffusivity is shown in blue (B), decreased axial diffusivity is shown in orange-yellow (C), and increased radial diffusivity is shown in gray (D). A, Corticospinal tract, superior longitudinal fasciculus, anterior thalamic radiation, inferior fronto-occipital fasciculus, uncinate fasciculus, and inferior longitudinal fasciculus; B, superior longitudinal fasciculus, thalamic radiation, inferior fronto-occipital fasciculus, uncinate fasciculus and inferior longitudinal fasciculus, corticospinal tract, cingulum, genu, and splenium of the corpus callosum; C, superior longitudinal fasciculus, thalamic radiation, inferior fronto-occipital fasciculus, uncinate fasciculus and inferior longitudinal fasciculus, corticospinal tract, cingulum, and genu and splenium of the corpus callosum; D, superior longitudinal fasciculus, thalamic radiation, inferior fronto-occipital fasciculus, uncinate fasciculus, inferior longitudinal fasciculus, corticospinal tract, and cingulum.
Neuroanatomic regions with reduced FA in patients with pSS compared with control subjects
Voxelwise-based group comparison of MD, AD, and RD between patients with pSS and control subjects showed increased MD and RD and decreased AD in the patients in a widespread, diffuse pattern involving most of the major WM tracts throughout the brain (Fig 1B, 1C, and 1D, respectively).
Finally, regression analysis did not show any significant association between the DTI metrics and disease duration.
Discussion
The main findings in the present study were decreased FA values in patients with pSS compared with control subjects in multiple major WM tracts, including the corticospinal tract, the superior longitudinal fasciculus, the anterior thalamic radiation, the inferior fronto-occipital fasciculus, the uncinate fasciculus, and the inferior longitudinal fasciculus. These changes in FA are related to widespread decreased AD values and increased MD and RD values in patients with pSS in comparison with control subjects.
The decreased FA and the increased MD values found in the present study are in accordance with the findings of the single previous study evaluating DTI in a similar number of patients with pSS (n = 19), which used the ROI method to evaluate 2 areas in the frontal lobes.16 The present study adds to these previous findings by use of the TBSS analysis, which permits voxelwise statistical analysis of all DTI data and revealed more extensive changes in patients with pSS.18
FA is a measure of the degree to which water diffusion is constrained in the brain, and its primary determinant is the packing attenuation of axons within a voxel.15,22 Axonal packing attenuation encompasses a variety of microstructural level variables (eg, degree of myelination, axonal diameters, and extracellular space).15,22 The decreased FA values demonstrated in patients with pSS thus indicate loss of WM fiber integrity. MD quantifies the amount of diffusion within a brain voxel, but it lacks directional information, and increased MD values point to an increase overall in directionally nonspecific water diffusivity and suggest tissue breakdown with an increase in brain-water content.15
The MD and FA indices allow for quantitative evaluation of the random translational motion of water molecules and have been shown to reflect a variety of pathologic states in the brain,23,24 but, despite being sensitive, they are lacking in specificity. The changes in diffusion-tensor eigenvalues λ1, λ2, and λ3, which are the source indexes for calculating the MD and FA, may provide further information about the underlying neuropathologic mechanisms.25,26 The λ1 or AD measures the diffusion coefficient along the direction of maximum diffusivity and reflects changes in restrictive barriers along the direction of a tract.25,26 The (λ2+λ3)/2 or RD measures the diffusion coefficient perpendicular to the direction of maximum diffusivity and reflects mainly changes in the axonal membrane and the myelin sheath.25,26 It should be noted that AD and RD, unlike FA, are not rotationally invariant metrics of the diffusion tensor, and therefore their values depend on the orientation of the diffusion tensor ellipsoid at each voxel. Consequently, their usage in multisubject studies requires precise registration that ensures alignment of the different diffusion tensor ellipsoids in every voxel; otherwise, interpretation of changes of the “axial” and “radial” diffusivities on the basis of the underlying tissue structure becomes problematic.27 Experimental studies have demonstrated that axonal damage leads to a marked decrease in AD and modest, often insignificant, decreases in RD, whereas demyelination increases RD without changing AD when these phenomena take place in isolation.28,29 In the present study, the decreased AD and increased RD may be suggestive of decreased organization of tracts or axonal damage, which prevents diffusion along the long (axial) axis, or alternatively, decreased myelination that allows for more room between axons for water molecules to move perpendicular to the tract (ie, radially).
Several pathogenetic mechanisms may account for the findings in the present study: 1) WM involvement could reflect axonal damage through anterograde or Wallerian degeneration. Cerebral small-vessel vasculitis has been reported in patients with pSS, leading to hypoperfusion of the cortex and atrophy.12,13,30,31 With the use of voxel-based morphometry, diffuse GM atrophy has been demonstrated in patients with pSS.8 This cortical atrophy may induce Wallerian degeneration of WM tracts.8 Wallerian degeneration is characterized by a stereotypical course, starting with disintegration of axonal structures within days after injury, followed by fragmentation-degradation of myelin caused by infiltration of macrophages and finally, fibrosis, and atrophy of the affected fiber tracts.32,33 Loss of axonal structure may result in less restricted diffusion perpendicular to the main direction of fibers and consequently give rise to elevated RD.22,34 The membrane disintegration and cellular debris create new diffusion barriers that lead to a decrease in diffusivity parallel to the main fiber direction and therefore reduced AD22,34; 2) WM damage could be caused by a direct insult from antibodies against myelin. Recently, an autoantibody targeting the water channel protein aquaporin-4 (anti–aquaporin-4) has been discovered in neuromyelitis optica, leading to death of oligodendrocytes, demyelination, and axonal loss.35 It is well known that there is a strong association between pSS and neuromyelitis optica, and there is evidence of an increased prevalence of anti–aquaporin-4 antibodies in patients with pSS.36⇓–38 Loss of myelin reduces the barriers that restrict diffusion perpendicular to the WM fibers, leading to increased RD.28 Axonal damage is an integral part of demyelination, and the concurrent axonal transactions and disruptions associated with axonal damage lead to the addition of diffusion barriers parallel to the axon and thus to reduced AD.39 Finally, the underlying process accounting for the WM involvement might be a combination of Wallerian degeneration and demyelination, which is supported by a limited number of autopsy studies that have demonstrated both axonal degeneration and demyelination in the CNS of patients with pSS.40,41 Similarly, decreased FA and AD and increased RD have been demonstrated in patients with neuromyelitis optica.42,43 These changes involve multiple, major WM tracts in the normal-appearing WM, accumulate with increasing disease duration, and are probably related to both demyelination and Wallerian degeneration.42,43
A global decrease in WM integrity throughout the brain was observed in multiple major fiber tracts that control a wide range of brain functions. The CNS manifestations in pSS are heterogeneous, manifested as focal or diffuse involvement.5,6 The wide spectrum of CNS manifestations in pSS includes movement disorders, motor and sensory loss, seizures, cognitive impairment, dementia, psychiatric abnormalities, encephalopathy, optic neuropathies, and others.5,6 This wide range of CNS manifestations in patients with pSS is in accordance with the diffuse WM involvement throughout the brain detected by DTI in the present study.
The present study has some certain limitations that must be noted. First, the sample size was relatively small, which might reduce the power of the statistical significance and the generalization of the findings. Second, the cross-sectional design of the study and the reliance on a single imaging technique do not allow for assessment of the relative time that WM integrity is affected. Third, our DTI acquisition sequence had a limited number of directions that might have limited the precision of the values of the DTI parameters. Furthermore, the TBSS analysis that was used in the present study (as in any other analyses) is not without flaws or limitations. One is the lack of accuracy caused by the low-resolution DTI (partial volume effect). However, this is more an acquisition limitation rather than a postprocessing flaw, and TBSS remains one of the most reliable methods for multisubject DTI analysis available today. Finally, the lack of official cognitive or neuropsychiatric testing and reliable correlations with our imaging findings weaken their clinical interpretation/value. Notwithstanding the problematic issues raised, this study adds to the limited existing literature on CNS involvement in patients with pSS and should be considered preliminary. Therefore, future studies in larger populations, as well as prospective, longitudinal studies with the use of additional advanced MR imaging techniques with histopathologic and clinical correlations, will be useful for further probing the nature and the relationship of WM abnormalities to clinical symptoms in pSS.
Conclusions
Loss of WM microstructural integrity is demonstrated in patients with pSS as reduced FA and AD values and increased MD and RD values. Involvement of the WM might be caused by demyelination or through anterograde Wallerian degeneration. Additional histopathologic and advanced MR imaging studies evaluating larger series of patients will be necessary to better elucidate the role of CNS involvement in patients with pSS and to improve the understanding of the pathologic mechanisms behind the reduction of WM tract integrity in pSS.
References
- Received June 13, 2013.
- Accepted after revision August 10, 2013.
- © 2014 by American Journal of Neuroradiology