|Home | About | Journals | Submit | Contact Us | Français|
Chronic fatigue syndrome (CFS) is a complex illness, which is often misdiagnosed as a psychiatric illness. In two previous reports, using 1H MRSI, we found significantly higher levels of ventricular cerebrospinal fluid (CSF) lactate in patients with CFS relative to those with generalized anxiety disorder and healthy volunteers (HV), but not relative to those with major depressive disorder (MDD). In this third independent cross-sectional neuroimaging study, we investigated a pathophysiological model which postulated that elevations of CSF lactate in patients with CFS might be caused by increased oxidative stress, cerebral hypoperfusion and/or secondary mitochondrial dysfunction. Fifteen patients with CFS, 15 with MDD and 13 HVs were studied using the following modalities: (i) 1H MRSI to measure CSF lactate; (ii) single-voxel 1H MRS to measure levels of cortical glutathione (GSH) as a marker of antioxidant capacity; (iii) arterial spin labeling (ASL) MRI to measure regional cerebral blood flow (rCBF); and (iv) 31P MRSI to measure brain high-energy phosphates as objective indices of mitochondrial dysfunction. We found elevated ventricular lactate and decreased GSH in patients with CFS and MDD relative to HVs. GSH did not differ significantly between the two patient groups. In addition, we found lower rCBF in the left anterior cingulate cortex and the right lingual gyrus in patients with CFS relative to HVs, but rCBF did not differ between those with CFS and MDD. We found no differences between the three groups in terms of any high-energy phosphate metabolites. In exploratory correlation analyses, we found that levels of ventricular lactate and cortical GSH were inversely correlated, and significantly associated with several key indices of physical health and disability. Collectively, the results of this third independent study support a pathophysiological model of CFS in which increased oxidative stress may play a key role in CFS etiopathophysiology.
Chronic fatigue syndrome (CFS) is a medically unexplained illness, diagnosed only after alternative medical and neuropsychiatric pathologies have been excluded. The pathognomonic symptom of CFS is debilitating fatigue lasting for at least 6 months, which is accompanied by at least four of eight infectious, rheumatologic or neuropsychiatric symptoms (1,2). The diagnosis remains controversial because of the disorder's heterogeneous phenomenology, extensive overlap with neuropsychiatric symptoms and lack of validated diagnostic tests. As a result, the discovery of CFS-specific biomarkers is currently an area of intense research.
In two previous studies, we used 1H MRSI to compare neurometabolites in CFS with those in generalized anxiety disorder (3) and major depressive disorder (MDD) (4), common neuro-psychiatric disorders with extensive symptom overlap with CFS. In these reports, patients with CFS showed significantly elevated ventricular cerebrospinal fluid (CSF) lactate relative to healthy controls (3,4) and to patients with generalized anxiety disorder (3), whereas no differences were found between patients with CFS and MDD (4). Importantly, our replicated finding of significant elevations of ventricular lactate in CFS (3,4) suggested a potential illness-associated biomarker, whose understanding could shed new light onto the pathophysiology of the illness.
In this study, we investigated a pathophysiological model of CFS, which postulates that sustained oxidative stress (5) and associated oxidant damage (6,7) lead to cerebral hypoperfusion (8–11) and/or to secondary mitochondrial dysfunction (5) which could potentially explain our observed cross-sectional elevations of ventricular lactate (3,4). Specifically, this study had two primary objectives: (i) to use 1H MRSI to replicate, in a new cohort, our finding of cross-sectional elevations of ventricular lactate in CFS (3,4); and (ii) to determine whether the postulated (5) and experimentally documented (6,7) oxidative stress increases in the disorder are associated with antioxidant capacity deficit using 1H MRS to measure in vivo brain levels of glutathione (GSH), the most abundant antioxidant in the central nervous system (12). Our secondary aims were as follows: (i) to use arterial spin labeling (ASL) MRI (13) to replicate previous observations of decreased regional cerebral blood flow (rCBF) in CFS (8–11), which may explain the observed lactate elevations; (ii) to use 31P MRSI to measure regional brain levels of high-energy phosphates (e.g. ATP) as indices of a possible secondary mitochondrial dysfunction in CFS (5), whose presence might also be associated with elevations in lactate; and (iii) to assess whether these neuroimaging markers correlate with clinical characteristics in CFS.
All neuroimaging studies were conducted in a single 90–120-min examination, enabling the acquisition of temporally concordant data for potentially more meaningful comparisons between the outcome measures. As in our previous studies (3,4), age- and sex-matched healthy volunteers (HVs) and patients with MDD served as normal controls and ‘disease controls’, respectively.
Subjects enrolled into the study provided written informed consent before all procedures, and were compensated for their participation. Diagnostic and laboratory assessments for patients with MDD and HVs were conducted at Mount Sinai School of Medicine, and those for patients with CFS were conducted at Beth Israel Medical Center. All neuroimaging scans and additional day-of-scan clinical evaluations were performed at Weill Cornell Medical College. Approval to conduct the studies was granted by the institutional review boards of the participating institutions.
For all participants, general eligibility requirements included age between 18 and 45 years, negative urine toxicology at screening and on the day of scan and, for females of reproductive age, the use of an effective birth control method and negative day-of-scan pregnancy test. In addition, all participants were instructed to refrain from the consumption of alcohol for at least 48 h prior to the neuroimaging scans, as the measured lactate CH3 group has nearly the same MR frequency as the ethanol CH3 group at 1.3 ppm. Exclusion criteria consisted of any unstable medical or neurological illness or any condition precluding MRI exposure (e.g. pacemaker, metallic prosthesis). All participants were assessed for psychiatric disorders using the Structured Clinical Interview for the Diagnostic and Statistical Manual of Mental Disorders, 4th edn. (14).
Patients with CFS were recruited via clinician referrals and media advertisements and were diagnosed by a board-certified neurologist with extensive CFS research experience (BHN), using modified US Centers for Disease Control and Prevention (CDC) guidelines (1). Patients with MDD were recruited through similar methods and were evaluated by a board-certified psychiatrist (SJM or JWM). All participants had been psychotropic medication free for 2 weeks or more (≥ 4 weeks for fluoxetine) prior to the scans, had been free of substance abuse/dependence for 6 months or more, had no lifetime history of psychotic disorder, mania or hypomania, had no pervasive developmental disorder or mental retardation, and had no current eating disorder. HVs, recruited through local media advertisement, did not meet the criteria for CFS or any psychiatric disorder, and were group matched for sex and age to the two patient groups.
Within 1 week of the scan, participants underwent a standardized clinical assessment battery, which included the following: (i) CDC Symptom Inventory for CFS symptoms (1); (ii) fatigue severity using the 20-item Multidimensional Fatigue Inventory (MFI) (15); (iii) physical and mental health using the 36-Item Short Form Health Survey by the RAND Corporation, Santa Monica, CA, USA (RAND SF-36) (16); (iv) functional impairment and disability using the Sheehan Disability Scale (SDS) (17); (v) pain and interference of pain with activities of daily life using the Wisconsin Brief Pain Inventory (18); (vi) sleep quality using the Pittsburgh Sleep Quality Index (19); (vii) depression using the Quick Inventory of Depressive Symptomatology – Self Report (QIDS-SR) (20); and (viii) early life traumatic events using the Childhood Trauma Questionnaire (21).
All neuroimaging studies were conducted on a research-dedicated, multinuclear, General Electric 3.0-T EXCITE MR system at Weill Cornell Medical College.
The in vivo levels of ventricular lactate and the lateral ventricular volume (Fig. 1) which contributed signal to the measured lactate were obtained in all subjects using a commercial, double-quadrature, dual-tuned, 1H/31P head coil (Clinical MR Solutions, LLC, Brookfield, WI, USA), employing a multislice 1H MRSI technique (22) and volumetric analysis of structural MRI data acquired with a standard three-dimensional, T1-weighted, spoiled gradient-recalled echo technique, respectively, as described previously (3,4). The lactate levels presented are the mean values obtained for all voxels in the ventricular space (Fig. 1A, C), expressed in institutional units (i.u.) as ratios relative to the root-mean-square (rms) of the background noise in the spectra, and corrected for the contribution of finite noise at the lactate frequency to ensure a mean lactate value of 0.0 ± SD for the HV group (4). Negative lactate values thus simply re flect random noise fluctuations at the lactate frequency (4).
31P MRSI data were recorded with the same dual-tuned, double-quadrature, 1H/31P volume head coil as used for the 1H MRSI scans from a single 30-mm slice, prescribed to be co-axial with the central slice in the 1H MRSI lactate scan (Fig. 1A), thereby ensuring overlap of the regions of interest in the two scans for enhanced correlation. A pulse sequence (‘FIDCSI’) supplied by the instrument manufacturer and consisting of a single slice-selective radiofrequency (RF) pulse, followed by a slice refocusing gradient and phase-encoding gradients, was implemented with 14 × 14 phase-encoding steps, a field of view of 420 mm, 2048 sample points, a spectral width of 5 kHz, TR = 1000 ms and eight excitations per phase-encoding step, to yield voxels with a nominal size of 3.0 × 3.0 × 3.0 cm3 in a total scan time of 26 min. The field homogeneity in the 31P MRSI slice was automatically optimized by the host computer on the water 1H signal from the same slice detected with the 1H coil prior to initiating the 31P acquisition. The raw data matrix size was zero filled to 32 × 32 spatial points prior to standard three-dimensional Fourier transformation to yield a grid of spectral voxels (Fig. 2) which was selected for further post-processing and statistical analysis. Although our use of a relatively short TR (1000 ms) probably resulted in substantial saturation of all 31P resonances, this did not adversely affect the signal-to-noise ratio (Fig. 2).
Following the 1H and 31P MRSI scans, the dual-tuned 1H/31P head coil was replaced with an eight-channel phased-array head coil to measure levels of GSH in 15 min from a single 3 × 3 × 2-cm3 occipital cortex (OCC) voxel (Fig. 3A), a region in which we had previously found abnormal levels of γ-aminobutyric acid (GABA) in patients with MDD (32), and is thus of interest as depressive symptoms are common among CFS sufferers. It should be noted that, because of its low brain concentration (0.5–1.5 mmol/L) and overlap by a total creatine (tCr) resonance at 3.0 ppm, which is up to 10-fold stronger (Fig. 3B, spectrum a), in vivo brain GSH is technically challenging to detect by 1H MRS. However, using the J-edited spin echo difference technique (23,24), we reliably achieved unobstructed detection of the GSH β-cysteinyl resonance at 2.98 ppm in this study (Fig. 3B, C). Briefly, a standard point-resolved spectroscopy sequence was turned into a volume-selective J-editing sequence by inserting a pair of frequency-selective inversion RF pulses, flanked by spoiler gradient pulses of opposite signs, before and after the second 180° RF pulse of the double spin echo. The application of these ‘editing’ RF pulses at the frequency (4.56 ppm) of the GSH α-cysteinyl resonance on alternate scans, with TE/TR = 68/1500 ms, alternately inverts the GSH β-cysteinyl resonance at 2.98 ppm by inhibiting and allowing its J modulation. The subtraction of the two sub-spectra acquired yields the desired GSH resonance at 2.98 ppm (Fig. 3B), with complete elimination of the overlapping tCr singlet, which is not J modulated. Figure 3C shows a stacked plot of edited spectra from five subjects, with the GSH peak clearly and consistently detected at 2.98 ppm. To obtain the area under this peak, a robust and highly optimized public-domain Levenberg–Marquardt nonlinear least-squares minimization routine (25), written in Interactive Data Language (IDL; Visual Information Solutions, Boulder, CO, USA), was implemented to selectively fit the GSH resonance in the frequency domain (Fig. 3B, spectrum d) to a pseudo-Voigt function, which enables more precise analysis of resonances that consist of mixtures of Lorentzian and Gaussian lineshapes (26), as is often the case for in vivo spectra. For group comparisons and reporting, the resulting GSH peak areas were expressed semi-quantitatively as ratios relative to the area of the simultaneously acquired unsuppressed water resonance (W) in the OCC voxel.
As metabolite level differences have been reported between gray matter (GM) and white matter (WM) (27,28), we implemented volumetric MRI-based tissue segmentation to assess brain matter content and heterogeneity of our relatively large OCC voxel. MEDx software (Medical Numerics, Sterling, VA, USA) was used to segment the brain tissue based on the signal intensity histogram of each subject's volumetric (spoiled gradient-recalled echo) MRI. From the histogram, a segmentation mask of the OCC voxel was generated and the proportions of GM, WM and CSF for the voxel were computed. These were then compared between groups for use as covariates in the case of significant differences.
rCBF was measured in each participant with the same eight-channel head coil as used for the GSH scans, employing a pseudo-continuous spiral gradient echo ASL pulse sequence (13). A field of view of 24 cm and an acquisition matrix of 512 × 8 resulted in a reconstructed matrix of 128 × 128 and a resolution of 0.53 mm × 0.53 mm × 3.8 mm over 48 slices. Acquisition parameters included TR/TE = 10.4/2.4 ms, a labeling delay of 1500 ms and a bandwidth of 62.5 kHz, resulting in a total scan time of 5.25 min. Water spins were labeled at the cervico-medullary junction, approximately 1 cm superior to the cerebellum, using velocity-driven adiabatic inversion with an amplitude-modulated control labeling pulse, which reproduces the frequency-dependent off-resonance effects of the labeling pulse, and therefore permits multislice imaging. A fast T1 map was generated using the same sequence to allow the quantification of rCBF. A reference T1-weighted anatomic imaging series, matching the ASL acquisition in slice placement and thickness, was acquired for anatomical reference. Finally, an axial, high-resolution, isotropic, three-dimensional BRAVO (BRAin VOlume imaging) sequence was acquired with a spatial resolution of 1.1 mm × 1.1 mm × 1.5 mm. Acquisition parameters included TR/TE/TI = 8.3/1.6/725 ms, a flip angle of 7°, field of view of 24 cm, matrix of 256 × 192 and 124 slices for overlay of the rCBF data.
rCBF maps were calculated in real time on the MRI scanner to solve the following equation pixel by pixel:
where λ represents the brain–blood partition coefficient reflecting the water distribution between intravascular and extravascular compartments, and has been determined to assume the value of 0.9 mL/g, T1app is the apparent T1 of tissue as altered by the labeled blood, and Mcontrol and Mtagged are the longitudinal magnetizations per gram of brain tissue for the control and tagged acquisitions derived from recorded ASL data.
The resulting rCBF images were normalized to the Montreal Neurological Institute positron emission tomography template using Statistical Parametric Mapping software version 5 (Wellcome Department of Imaging Neuroscience, London, UK). An 8-mm smoothing kernel was applied to the source images during normalization. A Bayesian framework was used to determine the optimal 12-parameter affine transformation needed to normalize the rCBF image set to the Montreal Neurological Institute template. Groupwise voxel-based analysis of the ASL data was performed, followed by between-group t-test comparisons in which age and gender were regressed out of the model as independent covariates. The threshold on voxel-level rCBF values was set at an uncorrected p ≤ 0.001.
Prior to the analyses, the normality of all data was assessed using the Shapiro–Wilks test. The Type I error rate (i.e. chance findings) was protected against by limiting the statistical tests only to our specified hypotheses. Differences in the means of the primary outcome measures (GSH and CSF lactate) between the three groups were assessed using one-way analysis of variance (ANOVA), with Bonferroni correction for multiple comparisons. For all normally distributed secondary outcome measures, three-group comparisons were performed using general linear models (ANOVA), followed by post-hoc comparisons with Tukey's honestly significant difference test. Additional factors and covariates were chosen based on their previously reported influence on MRS markers. In the case of inhomogeneity of variance, Welch tests were performed for global group comparisons and Games–Howell tests for post-hoc pairwise comparisons.
Group comparisons for non-normally distributed data used the Kruskal–Wallis nonparametric test, with Mann–Whitney U-tests for post-hoc analyses. Exploratory associations between continuous variables were conducted using Pearson's product–moment correlation (r) or, for non-normally distributed data, Spearman ρ. All analyses were two-tailed, with the level of significance set at α < 0.05, unless otherwise noted. All statistical analyses were conducted using SPSS version 18.0 (SPSS Inc., Chicago, IL, USA).
With a sample size of, at most, 15 subjects per group, this study did not have sufficient statistical power to enable the assessment of within-group correlations, which we estimate would require 40 or more subjects per group to be meaningful. However, because most of the clinical variables were dimensional or continuous, rather than categorical, and none of these variables had been specified as inclusion or exclusion criteria for any of the groups, the full range of clinical scores remained possible for any of the participants after enrollment (i.e. a subject who met the eligibility criteria for the HV group could still score poorly in a given clinical test). Therefore, to increase the statistical power to test for potential associations between the continuous MRS measures and the dimensional clinical scores, the full participant sample was treated as a continuum, with no a priori assumptions about group membership other than recognizing that there would be a ‘ceiling effect’ in clinical scores because of the absence of discrete levels beyond the cut-off point representing ‘normal’ test results. As a result of this caveat, we only report clinical correlations that reached significance at the very conservative threshold of Pearson r or Spearman ρ of at least 0.5 and p = 0.001, except where noted. Treating the three groups as a continuum also enabled the assessment of the overlap between the outcome measures to establish their potential utility in differentiating between the groups.
Tables 1 and and22 provide the demographic and clinical characteristics, respectively, for subjects with viable data, consisting of 15 patients with CFS, 15 with MDD and 13 HVs, after excluding five patients with MDD, one with CFS and one HV with poor quality data mainly because of severe head motion, as determined by the presence of lipid peaks in voxels outside the boundaries of the brain in the 1H MRSI scans. The three groups were comparable for age, sex, ethnicity, marital status, body mass index (BMI) and years of education (Table 1). However, they differed by race (p = 0.001), as the CFS group consisted of Caucasian patients only.
The CFS group showed more CFS-associated symptoms relative to both the HV (p < 0.001) and MDD (p = 0.019) groups. There was a main effect of group for depression severity (p < 0.001), with post-hoc tests showing greater QIDS-SR scores in the MDD group relative to both the CFS (p < 0.001) and HV (p < 0.001) groups. There were also more depressive symptoms in patients with CFS relative to the HV group (p = 0.003).
The two patient groups exhibited a high degree of functional impairment and disability relative to the HV group, as measured by the Sheehan Disability Scale, but did not differ from each other (p = 0.233) (Table 2). In general, health quality across multiple domains was markedly reduced in the CFS group relative to the HV and MDD groups. The CFS group reported lower physical functioning, greater physical limitations, greater fatigue, higher pain scores, lower social functioning and worse general health scores than the HV group (Table 2). The MDD group fared poorly relative to the CFS group only in emotional well-being and associated limitations, but the two groups did not differ in social functioning or overall fatigue.
Using the MFI, fatigue severity was assessed in all subjects for the following five dimensions: general fatigue (GF), physical fatigue (PF), mental fatigue (MF), reduced motivation (RM) and reduced activity (RA). There was a significant main effect of group for fatigue severity as measured along each dimension of the MFI (Table 2). Post-hoc testing revealed that the group differences were driven by greater fatigue in both the CFS and MDD groups relative to the HV group, without any significant difference between CFS and HV with regard to motivational level, or between CFS and MDD with regard to activity level, mental fatigue or general fatigue (Table 2).
Mean 1H MRSI-derived ventricular lactate levels differed significantly (Fig. 4A) between the CFS (2.23 ± 1.10 i.u.), MDD (1.35 ± 0.93 i.u.) and HV (0.00 0.86 i.u.) groups (F2,33 = 16.78; p < 0.001). Post-hoc testing revealed significantly higher ventricular lactate levels in both CFS (p < 0.001) and MDD (p = 0.009) groups relative to HV, with a weak trend towards higher levels in CFS relative to MDD (p = 0.114). When controlling lactate levels for age and other covariates (gender, smoking status and BMI), the significant main effect of group remained (F2,2.4 = 13.634, p = 0.050). The means of the ventricular volumes from which the observed lactate signal was measured did not differ between the CFS (3554 ± 457 mm3), MDD (3971 ± 643 mm3) and HV (3057 ± 503 mm3) groups (p = 0.48). As expected, there was an association between ventricular volume and age (ρ = 0.437, p = 0.004) across all participants, but ventricular lactate did not correlate with ventricular volume (r = 0.263, p = 0.127).
There was a main effect of group on occipital GSH levels normalized to the peak area of the unsuppressed voxel tissue water (W) (F2,40 = 15.93; p < 0.001) which, on post-hoc testing, was attributable to reductions in GSH/W (Fig. 4B) in both the CFS [(1.86 ± 0.47) × 10–3; p < 0.001] and MDD [(2.27 ± 0.42) 10–3; p = 0.004] groups relative to HVs [(2.92 ± 0.59) × 10–3]. There was a nonsignificant trend towards lower GSH/W in the CFS relative to MDD group (p = 0.086). Levels of the reference voxel water signal did not differ between the groups (p > 0.5), so that, henceforth, we use GSH to mean GSH/W. Age, gender, smoking status and BMI did not have an impact on cortical GSH differences, nor did tissue type, which did not differ between the groups for GM (CFS, 10 021 ± 3158 mm3; MDD, 11 280 ± 2307 mm3; HV, 10 716 ± 2282 mm3; p = 0.656) or WM (CFS, 5427 ± 1814 mm3; MDD, 5351 ± 1889 mm3; HV, 6197 ± 2004 mm3; p = 0.530).
We found a strong negative correlation between ventricular lac-tate and GSH across all participants (r = –0.545, p = 0.001; Fig. 5), as well as within the CFS group (r = –0.544, p = 0.044), although the latter correlation was of lower statistical significance. This inverse relationship between ventricular lactate and GSH was such that most of the clinical variables found to correlate positively or negatively with one MRS variable were inversely correlated with the other, as illustrated in Fig. 6.
Across the three groups, significant associations between GSH and the following clinical variables were found (Table 3): (i) RAND SF-36 scores showed GSH to be positively correlated with physical functioning (Spearman ρ = 0.506; p = 0.001) and physical energy (ρ = 0.606, p < 0.001); (ii) GSH was negatively correlated with general fatigue (ρ = 0.565, p < 0.001), physical fatigue (ρ = –0.525, p < 0.001) and reduced activity (ρ = 0.508, p < 0.001), as assessed by MFI; (iii) GSH was negatively correlated with physical impairment and disability (ρ = –0.615, p < 0.001), as measured by SDS; and (iv) GSH was negatively correlated with CFS symptomatology (ρ = –0.633, p < 0.001), as determined by the CDC CFS Symptom Inventory.
For ventricular lactate, likewise, multiple significant correlations, but opposite to those for GSH (Fig. 6), were found across the whole sample (Table 3). Health quality across nearly all domains of RAND SF-36 negatively correlated with ventricular lactate (ρ = –0.601, p < 0.001 for total RAND SF-36). In addition, there were significant positive correlations between ventricular lactate and general fatigue (ρ = 0.537, p = 0.001), physical fatigue (ρ = 0.596, p < 0.001) and reduced activity (ρ = 0.496, p = 0.002). Moreover, ventricular lactate positively correlated with the degree of physical impairment and disability (ρ = 0.581, p < 0.001); however, there was only a weak positive correlation between ventricular lactate and CFS symptoms (ρ = 0.374, p = 0.038).
Finally, sleep quality across all participants, as assessed by the Pittsburgh Sleep Quality Index, was negatively correlated with GSH (r = –0.606, p = 0.001) and positively correlated with ventricular lactate (r = 0.601, p = 0.001), with higher PSQI scores signifying greater sleep disturbance.
Following intensity and morphological normalization and statistical analysis by Statistical Parametric Mapping software version 5, we found significantly different absolute CBF values, at a voxel-level significance threshold of p < 0.001, uncorrected for multiple comparisons, in two brain regions. The CFS group had lower rCBF values in the left anterior cingulate cortex (p = 0.039) and right lingual (p = 0.016) regions relative to the HV group (Fig. 7). In addition, there was a trend towards lower rCBF in the left anterior cingulate cortex in the MDD group relative to HVs (p = 0.08). There were no significant differences in rCBF values between the CFS and MDD groups in any brain region.
Although only two relatively limited regions emerged as having statistically significant CBF, examination of the raw values was highly suggestive of potentially more widespread hypoperfusion in patients with CFS and MDD. Figure 8 shows the percentage changes in mean rCBF (%ΔrCBF) measured relative to the HV group for and between the CFS and MDD groups (Fig. 8A), and the percentage changes in mean global CBF (%ΔgCBF ) relative to HV for and between the CFS and MDD groups (Fig. 8B). For nearly all individual regions (Fig. 8A) and globally (Fig. 8B), the mean %ΔrCBF and %ΔgCBF values for both the CFS and MDD groups were numerically lower than those for the HV group (Fig. 8A, B), whereas those for the CFS group were numerically lower than those for the MDD group. However, these numerical differences failed to reach statistical significance.
None of the exploratory statistical analyses conducted on the levels of high-energy phosphates, including ATP, creatine phosphate (PCr) and inorganic phosphate (Pi), as well as on their ratios (e.g. ATP/PCr or Pi/Cr) and on the intracellular pH derived from Pi chemical shifts, revealed a main effect of group in any brain region. In addition, there were no correlations between any of the phosphate metabolites and demographic and clinical variables.
This report marks the third in a series of neuroimaging investigations of CFS, a disorder with an uncertain pathophysiology. In this independent cohort, we replicated our previously reported elevations of ventricular CSF lactate in patients with CFS relative to age- and sex-matched HVs. In addition, we report, for the first time, significant reductions in the levels of cortical GSH in patients with CFS and MDD relative to healthy controls. Ventricular lactate and cortical GSH were inversely correlated with each other, as well as with several key indices of physical health and disability. These findings further implicate ventricular lactate elevations in patients with CFS, whereas decreases in the antioxidant GSH in patients with CFS and MDD strongly suggest a role of increased oxidative stress in the two disorders.
Although, in the second study of this series (4), we did not find any ventricular lactate differences between HVs and patients with MDD, or between patients with CFS and MDD, despite numerically different mean values, in the present study, we found significantly higher ventricular lactate levels in patients with MDD relative to HVs, and weak trend elevations in the CFS relative to the MDD group. This may indicate that the CFS and MDD patient groups recruited for the present study were more representative of each disorder, as the CFS group exhibited more CFS symptoms (by CDC Symptom Inventory) and physical fatigue (by MFI) than did the MDD group, and patients with MDD were significantly more depressed than those with CFS (by QIDS-SR), whereas the two groups were well matched with respect to overall fatigue and functional disability (Table 2). The improved group matching might have decreased the overlap in key symptoms, thereby decreasing the heterogeneity within each group and the associated variance. These considerations suggest that higher ventricular lactate in patients with CFS than in those with MDD might be found in larger studies. However, the relatively small magnitude of the difference suggests that ventricular lactate alone is unlikely to be a biomarker for differentiating between individual patients with CFS and MDD.
In this study, we found significant GSH deficits in patients with CFS and MDD relative to HVs, and a nonsignificant trend towards lower GSH in patients with CFS relative to those with MDD. To our knowledge, this is the first study to document GSH deficits in patients with MDD, consistent with recent hypotheses regarding redox dysregulation and increased oxidative stress in the disorder (29,30). Our finding of significant GSH deficits in patients with CFS is in contrast with a recent 1H MRS study, which failed to detect cortical GSH differences between patients with CFS and healthy subjects (31). However, that study had a number of notable methodological limitations. First, a conventional single-voxel point-resolved spectroscopy 1H MRS technique, which does not yield a detectable GSH signal, was used to measure this weak signal in the presence of an overlapping tCr peak at 3.0 ppm (see Fig. 3B, spectrum a), which is up to 10-fold stronger, requiring spectral simulation of an invisible resonance. Second, a relatively long TE of 144 ms was used, which unavoidably led to significant losses in the signal intensity of the already weak GSH peak as a result of J-coupling evolution and T2 decay. Finally, although a spectral estimation method was implemented to derive computationally the GSH peak levels despite the overlap, no spectral data showing the quality of the model fitting, or quantitative estimates of the associated errors, were provided to enable the assessment of the reliability of the fitting procedure. It is therefore highly unlikely that the reported ‘GSH’ values of 2.703 ± 2.311 (standard deviation) for patients with CFS and 5.191 ± 8.984 for healthy controls (31) could have been distinguished from background noise or from mixtures of overlapping and co-resonant GABA and mobile macromolecules. This is in contrast with the present study, which used a spectral editing technique that yields an unobstructed GSH peak (Fig. 3B, C), and has been employed extensively by our group (32) and others (23,33–35) to measure similarly low concentrations of brain GABA under equally challenging in vivo conditions.
The present 1H MRS measures of brain GSH are difficult to compare with previous peripheral measures in patients with CFS, which did not differ from controls (6,36,37), as GSH has poor blood–brain barrier permeability and transport properties (38), such that blood levels may not reflect cortical levels. A consequence of the poor blood–brain barrier permeability of GSH is that dietary supplements of the antioxidant are unlikely to increase the cortical levels significantly, justifying the heightened interest in synthetic GSH precursors, such as N-acetylcysteine, that can cross the blood–brain barrier to spur in situ synthesis and elevations of cortical GSH (38). Thus, our finding of significant GSH deficits in patients with CFS relative to HVs, which correlated strongly with various indices of physical health and functional disability, suggests that the investigation of such indirect approaches for the elevation of cortical GSH for neuro-protection in CFS might be warranted.
The results of our clinical evaluations of the CFS and MDD groups were consistent with the overall manifestations of each disorder, whilst also revealing significant symptom overlap across a wide range of indices of physical and functional disability, which may account for the difficulty in differentiating between the two disorders. The two groups differed with regard to the more pronounced depressive symptoms, reduced motivation and more ‘limitations due to emotional difficulties’ and less ‘emotional well-being’ in MDD, and more pain, physical limitations and physical fatigue, and poorer general health and physical functioning, in CFS.
Although GSH correlated strongly with CFS symptoms, correlations between CFS symptoms and ventricular lactate were much weaker, potentially suggesting a tighter coupling of GSH deficits than lactate, whose elevations might be more distal to pathophysiology, with CFS symptoms. Importantly, examination of the relationships between the two MRS variables and clinical indices across the full sample (Fig. 6; Table 3) revealed such a high degree of overlap that these MRS variables, alone or in concert, cannot be considered as sensitive biomarkers of CFS or MDD, or a basis for differentiating between the two disorders. However, the fact that we found direct correlations between GSH and clinical indices of good health or function, and between ventricular lactate and indices of poor health and physical disability, across multiple domains is qualitatively consistent with the beneficial role of GSH in detoxification and in mitigating the effects of oxidative stress, and with the role of lactate as a marker of nonsalutary biological challenge, such as defective cell energetics or infarction.
In this study, we found limited regions of hypoperfusion in patients with CFS relative to HVs, a finding that is in agreement with most (8–11,39), but not all (40,41), previous measurements of rCBF using single-photon emission computed tomography (SPECT) (8,9,39–41), Xe computed tomography (10) or ASL (11). In nearly all studies that found CBF deficits in patients with CFS, multiple regions were reported to be involved (8,9,11,39), with some of the regions of hypoperfusion being as small as or less than 1 cm (33). These findings of extensive areas of hypoperfusion in patients with CFS suggest that our failure to find widespread regions with decreased rCBF in this study, despite numerically lower values in patients with CFS and MDD relative to HVs (Fig. 8A, B), was probably the result of a Type II statistical error, i.e. our study was insufficiently powered to detect more widespread rCBF deficits by ASL. As in the present study, one of the published studies (41) compared rCBF in CFS, MDD and HV groups, but measured with SPECT, and found decreases in the frontal lobe in the MDD group relative to both CFS and HV groups, with no differences between the last two groups. This is qualitatively in agreement with the present study, which found a trend towards decreased anterior cingu-late cortex CBF in the MDD group, a region in which we (32) and others (33–35) have previously found significant amino acid neurotransmitter abnormalities.
It is unclear whether the degree of hypoperfusion found in patients with CFS by these multiple techniques is sufficient to trigger increased glycolytic activity that could account for our observed elevations of ventricular lactate. However, at least three threshold levels of cell hypoxia have been described (42), one of which could plausibly trigger increased glycolytic rates in patients with CFS: (i) decreased oxygen tension, but the cell's ATP production remains able to meet energy demand through changes in phosphorylation state and increased glycolysis; (ii) an oxygenation level at which steady-state ATP levels can be maintained only by increased anaerobic glycolysis through the very inefficient Embden–Meyerhof pathway, which is not sustainable in a high-energy-consuming organ, such as the brain, leading to rapid ATP depletion; and (iii) an oxygenation level at which anaerobic glycolysis is no longer able to generate sufficient ATP to maintain cell function. As we did not find deficits in ATP in this study (see Results section), and the levels of ventricular lactate measured in patients with CFS are relatively small relative to those reported previously in fully symptomatic patients with primary mitochondrial disorders (43) or cerebral infarction (44), we postulate that the observed ventricular lactate elevations may be triggered by a decrease in oxygenation level corresponding to the first aforementioned threshold level (i.e. mild to moderate hypoxia), which may result from CBF decreases found in this and previous studies (8–11,39). However, direct evidence of altered tissue oxygenation in patients with CFS is necessary to establish a definite mechanistic link between rCBF decreases and the observed lactate elevations.
Taken in the context of recent studies reporting credible evidence of increased markers of oxidative stress in patients with CFS (6,7), the present study supports an emerging oxidative stress hypothesis (5), which implicates sustained oxidative damage as an integral aspect of CFS pathophysiology. Our findings of significant cortical GSH deficits in patients with CFS, which are associated with multiple clinical indices and physical disability, are consistent with studies showing significantly elevated levels of 8-iso-prostaglandin-F2α-isoprostanes (6,7), which are reliable blood markers of oxidative stress that are also correlated with clinical symptoms in CFS (45). Our explanatory model, depicted in Fig. 9, shows how increased oxidative stress might relate to and/or explain our findings of elevated ventricular lactate in three independent CFS samples.
The view that oxidative stress might play a key role in the pathophysiology of CFS, first coherently advanced by Pall (5), posits that an initial infectious process or associated immuno-logic response (‘cloud’ in model shown in Fig. 9) induces one or more cytokines (interleukin-1, interleukin-6, tumor necrosis factor-α, interferon-γ), which, in turn, induce nitric oxide synthase, leading to increased levels of nitric oxide (·NO) that subsequently react with superoxide (·OO–) to produce peroxynitrite (OONO), a potent pro-inflammatory and pro-apoptotic reactive oxygen/nitrogen species (ROS/RNS) and mediator of oxidative stress (45,46). In the absence of adequate antioxidant defense capacity or reserves, of which GSH is the most abundant and important component, peroxynitrite (Fig. 9, left-hand pathway) will react with arachidonic acid (AA) in cell membrane phospholipids, forming isoprostanes (46,47) and compromising membrane and cell function (5,45–47). Therefore, the results of recent studies (6,7) reporting significant elevations of blood levels of isoprostanes in patients with CFS, which correlated with clinical symptoms in one study (6), represent exceptionally strong evidence for the view that oxidative stress is implicated in the disorder. Importantly, our finding of significant cortical deficits in the primary intracellular antioxidant GSH in patients with CFS has added both mechanistic and face validity to this model by establishing that oxidative stress in the disorder is, at least in part, caused by decreased antioxidant capacity.
In addition to their role as markers of oxidative stress, isoprostanes have potent vasoconstrictor effects on peripheral vasculature, including cerebral arterioles (45,46). The finding of significantly elevated levels of these byproducts of membrane lipid peroxidation in patients with CFS may therefore also explain the results of this and multiple previous studies that have found reduced absolute (10) and relative (8,9,11,39) cortical and subcortical blood flow in patients with CFS. As cerebral hypoper-fusion is known to increase brain lactate (42), our finding of elevated CSF lactate, together with previous reports of increased markers of oxidative stress and evidence of hypoperfusion in patients with CFS, suggest a pathophysiological model of the disorder (Fig. 9) in which increased oxidative stress elevates blood levels of isoprostanes, which then constrict cerebral arterioles, leading to decreased cerebral blood flow and to attendant increases in anaerobic metabolism and lactate elevations.
Finally, as oxidative stress and the associated accumulation of free radicals can lead to mitochondrial dysfunction (48), the involvement of oxidative stress in CFS also raises the possibility that the observed cross-sectional elevations of ventricular lactate could be the result of secondary mitochondrial dysfunction (Fig. 9; right-hand branch of model), which would lead to ATP and energy production deficits that would be compensated by the upregulation of glycolytic activity, with attendant lactate production. However, our measurement of ATP and other high-energy phosphates that might be affected in the case of significant mitochondrial dysfunction found no abnormalities. Therefore, the results of this study do not support a significant involvement of mitochondrial energy metabolism in CFS, although our failure to find such evidence might be caused by the reported insensitivity of the 31P MRS method to subtle mitochondrial defects at rest (49,50).
A limitation of this pilot study is that our relatively small sample size did not provide sufficient power to enable reliable assessment of within-group correlations, or to confirm more widespread regions of hypoperfusion in patients with CFS. Therefore, these results require replication in larger samples. Moreover, we did not obtain circulating levels of isoprostanes in our study participants, inferring increases in the levels of this important component of the model from published findings, and our resting-state 31P MRSI assessment of mitochondrial dysfunction is inconclusive, because of the possibility that the unmasking of subtle mitochondrial metabolism defects by this technique might require the acquisition of spectra at rest and then during and following symptom provocation or a suitable task activation (49,50). Finally, the neuroimaging outcome measures assessed in this study, individually or in combination, do not seem likely to enable individual patients with CFS and MDD to be differentiated clinically. However, the results presented here have demonstrated the high overall potential of these methods to contribute to a new pathophysiological understanding by enabling in vivo testing of the consistency and validity of a promising model of CFS, and by providing a compelling rationale to undertake studies to investigate neuro-protective treatments (e.g. N-acetylcysteine) that can elevate cortical GSH to mitigate against the effects of oxidative stress-mediated damage in the disorder.
Contract/grant sponsors included the Chronic Fatigue and Immune Dysfunction Syndrome (CFIDS) Association of America, Inc. and the National Institutes of Health (contract/grant number: R01 MH-075895, K23 MH-069656).
We thank the patients who participated in this study and their families, and Mr Josefino Borja for his valuable contribution in operating the MRI scanner.