Abstract
BACKGROUND AND PURPOSE: Schizophrenia is well-known to be associated with hippocampal structural abnormalities. We used 1H-MR spectroscopy to test the hypothesis that these abnormalities are accompanied by NAA deficits, reflecting neuronal dysfunction, in patients compared with healthy controls.
MATERIALS AND METHODS: Nineteen patients with schizophrenia (11 men; mean age, 40.6 ± 10.1 years; mean disease duration, 19.5 ± 10.5 years) and 11 matched healthy controls (5 men; mean age, 33.7 ± 10.1 years) underwent MR imaging and multivoxel point-resolved spectroscopy (TE/TR, 35/1400 ms) 1H-MRS at 3T to obtain their hippocampal GM absolute NAA, Cr, Cho, and mIns concentrations. Unequal variance t tests and ANCOVA were used to compare patients with controls. Bilateral volumes from manually outlined hippocampal masks were compared by using unequal variance t tests.
RESULTS: Patients' average hippocampal GM Cr concentrations were 19% higher than that of controls, 8.7 ± 2.2 versus 7.4 ± 1.2 mmol/L (P < .05); showing no differences, concentrations in NAA were 8.8 ± 1.6 versus 8.7 ± 1.2 mmol/L; in Cho, 2.3 ± 0.7 versus 2.1 ± 0.3 mmol/L; and in mIns, 6.1 ± 1.5 versus 5.2 ± 0.9 (all P > .1). There was a positive correlation between mIns and Cr in patients (r = 0.57, P = .05) but not in controls. The mean bilateral hippocampal volume was ∼10% lower in patients: 7.5 ± 0.9 versus 8.4 ± 0.7 cm3 (P < .05).
CONCLUSIONS: These findings suggest that the hippocampal volume deficit in schizophrenia is not due to net loss of neurons, in agreement with histopathology studies but not with prior 1H-MR spectroscopy reports. Elevated Cr is consistent with hippocampal hypermetabolism, and its correlation with mIns may also suggest an inflammatory process affecting some cases; these findings may suggest treatment targets and markers to monitor them.
ABBREVIATIONS:
- CSI
- chemical shift imaging
- 1H-MRSI
- 3D multivoxel 1H-MRS imaging
- SZ
- schizophrenia
Schizophrenia (SZ) is a chronic psychiatric disorder that profoundly alters a person's perception, cognition, and behavior. Due to its high prevalence (∼1%), early onset, and the limited efficacy of existing treatments, SZ exacts enormous personal and economic tolls.1 The characteristic “positive” symptoms of SZ psychosis (delusions, hallucinations, disorganized speech, and behavior) are accompanied by cognitive decline and “negative” symptoms, including diminished emotional expression and avolition.2 Despite recognition as fundamental to the disorder as early as 1919, cognitive and negative symptoms were largely ignored once relatively effective pharmacologic treatments for the positive symptoms emerged in the 1950s.3 In recent years, the development of sophisticated investigative techniques has rekindled interest in the neurobiologic substrate underlying SZ so that therapy may be developed to target the full range of symptoms and alter the clinical course.
This surge of effort has resulted in the development of a plurality of sometimes competing, sometimes overlapping theories of SZ pathophysiology, from the classic neurotransmitter-based theories to genetic-, immune-, synapse-, and network-based theories. The hippocampus is implicated in many of these theories. SZ is known to be associated with reduced hippocampal volume, increased basal perfusion, decreased activation during certain memory tasks, decreased neurogenesis in the dentate gyrus, and reduced connectivity with cortical and subcortical regions.4⇓⇓⇓⇓–9 Recent findings show that conversion to psychosis in high-risk subjects is predicted by hypermetabolism in the hippocampal CA1 subregion.4 Psychosis may be conceptualized as a disruption in learning and memory involving impaired habituation and “runaway” pattern completion due to hippocampal hyperactivity.4,10,11 The mechanism underlying these hippocampal abnormalities remains unclear. The absence of gliosis on postmortem histopathology and reduced volume in prodromal and first-episode cases has fomented a shift from neurodegenerative to neurodevelopmental hypotheses.12 Lack of classic neurodegeneration with gliosis, however, does not imply the absence of subtler progressive damage and inflammation. Mounting genetic and epidemiologic data suggest a role for aberrant immune function and inflammation in SZ.13
1H-MR spectroscopy measures metabolites used as markers for underlying physiologic processes, most prominently: NAA (NAA and N-acetyl aspartylglutamate) for neuronal integrity, Cr (creatine and phosphocreatine) for energy metabolism, Cho (phosphocholine, choline, and glycerophosphocholine) for membrane turnover, and mIns (myo-inositol) for astroglial proliferation.14 Most previous hippocampal 1H-MR spectroscopy studies in SZ reported lower NAA15; one found increased Cr and Cho.15 None of 10 others that measured Cho, 7 that measured Cr, or 8 that measured mIns, found changes.16 Difficulty interpreting results may stem from diagnostic heterogeneity within and between studies, insufficiently powered samples, and methodologic variation.16
This study compares absolute hippocampal GM NAA, Cr, Cho, and mIns levels, obtained with 3D multivoxel 1H-MR spectroscopy imaging (1H-MRSI) at 3T, between patients with SZ and controls, to test the hypothesis that patients with SZ have decreased hippocampal GM NAA, reflecting neuronal damage. Higher field strength and 1H-MRSI (compared with single-voxel spectroscopy) yield better coverage of the irregular shape of the hippocampus and better SNR and spatial resolution.
Materials and Methods
Human Subjects
The patients with SZ or schizoaffective disorders were recruited from the outpatient clinics of New York University Langone Medical Center and Bellevue Hospital, diagnosed on the basis of the Diagnostic Interview for Genetic Studies conducted by clinicians trained for reliability, and had ongoing assessments.17 All were taking stable doses of medications and had no other psychiatric or neurologic disorders. Age-matched controls recruited from hospital postings met the criteria of no personal or family history of psychosis, no axis I disorder in the past 2 years, and no known neurologic disorder. Exclusion criteria for both patients and controls included uncontrolled medical illness, MR imaging contraindication or inability to tolerate an MR imaging examination, and substance (except tobacco) use in the past 6 months. Substance use history of the patients was the following: Six of 19 reported current or past tobacco use; 6/19, past cannabis use; 9/19, past cocaine use; and 7/19, past alcohol use. All participant demographics are compiled in Table 1, and all gave institutional review board–approved written consent.
MR Imaging Data Acquisition
All experiments were performed at 3T in a whole-body MR imaging scanner (Magnetom Trio; Siemens, Erlangen Germany) with a transmit-receive head coil (TEM3000; MR Instruments, Minneapolis, Minnesota). For anatomic reference, tissue segmentation, and 1H-MRSI VOI guidance, T1-weighted 3D MPRAGE images were obtained from each subject: TE/TI/TR, 2.6/800/1360 ms; matrix, 256 × 256; FOV, 256 × 256 mm2; 160 sections, 1-mm-thick each. These were reformatted into 192 axial, sagittal, and coronal sections at 1-mm3 isotropic resolution.
Our noniterative, B0 map–based, BOLERO (B0 Loop-Encoded ReadOut) in-house software adjusted the first- and second-order shims of the scanner to optimize the magnetic field homogeneity over the hippocampi in 3–5 minutes.18 A 6-cm anteroposterior × 9-cm left-right × 2-cm inferior-superior = 108-cm3 parallelepiped 1H-MRSI VOI was then image-guided over the bilateral hippocampus, as shown in Fig 1. This VOI was excited by using point-resolved spectroscopy (TE/TR, 35/1400 ms) with 2 second-order Hadamard-encoded slabs (4 sections) interleaved along the inferior-superior direction at every TR, as shown in Fig 1B, for optimal SNR and spatial coverage.19 Interleaving also enabled strong, 9 mT/m, Hadamard section-selection gradients, reducing the NAA→mIns chemical shift displacement to ∼0.13 mm, ∼3% of the section thickness.20 Thin sections were also chosen to reduce broadening from susceptibility gradients in the inferior-superior direction from the air-tissue interface with the maxillary sinuses just below, as seen in Fig 1A.
The 4 Hadamard section planes were encoded with 16 × 16 2D chemical shift imaging (CSI) over a 16 × 16 cm2 (left-right × anteroposterior) FOV to form 1.0 × 1.0 × 0.5 cm3 voxels. (Note that the actual voxel size [full width at half maximum of the point spread function] for such uniform 2D phase encoding is 1.12 × 1.12 × 0.5 = 0.63 cm3,21,22 because in the Hadamard direction, the nominal equals the actual voxel size23). The VOIs were defined in their planes by two 11.2-ms numerically optimized 180° radiofrequency pulses (4.8-kHz bandwidth) under 1.8 and 1.2 mT/m in the anteroposterior and left-right directions, to yield 9 × 6 × 4 = 216 voxels (Fig 1C). Such gradients lead the NAA VOI to experience a relative section in-plane chemical shift displacement of 2.3 mm in the anteroposterior and 3.5 mm in the left-right directions24 (ie, of the 6 × 9 × 2 cm3 anteroposterior × left-right × inferior-superior nominal VOIs, at least 5.7 × 8.6 × 2 cm2 [93.4%] is common for all metabolites).25 Note that Cho and Cr have even smaller relative displacements; the CSI localization grid does not experience this displacement.26 Because it is a relative shift, this error is encountered only at the VOI edges,25 and it is smaller than the 1 × 1 cm2 in-plane CSI resolution. Therefore, to avert it for all metabolites, we chose the in-plane VOI size large enough, 9 × 6 cm2, to have these displacement-error prone voxels at the VOI edges completely outside the hippocampus, as shown in Figs 1A and 2A. The MR spectroscopy signals were acquired for 256 ms at ±1-kHz bandwidth. At 2 averages, the 1H-MR spectroscopy was ∼25 minutes, and the examination took under an hour.
MR Spectroscopy Postprocessing
The MR spectroscopy data were processed off-line by using in-house software. Residual water was removed from the MR signals in the time domain27; the data was static-field drift-corrected,28 voxel-shifted to align the CSI grid with the NAA VOI, Fourier transformed in the temporal, anteroposterior, and left-right directions, and Hadamard reconstructed along the inferior-superior direction. Spectra were automatically corrected for frequency and zero-order phase shifts in reference to the NAA peak in each voxel.
Relative levels of the ith (i = NAA, Cr, Cho, mIns) metabolite in the jth (j = 1..216) voxel of the kth (k = 1..30) subject were estimated from their peak area, Sijk, by using parametric spectral modeling,29 with Glx, Cho, Cr, mIns, NAA, and taurine model functions, as shown in Fig 1D. The Sijk was scaled into absolute concentration, Cijk, relative to a 2 L reference sphere of Civitro = 12.5, 10.0, 3.0, and 7.5-mmol/L of NAA, Cr, Cho, and mIns in water at physiologic ionic strength: where SijR is the metabolite signals of the phantom; Vk180° and VR180° are subject and reference radiofrequency voltages for nonselective 1-ms 180° pulses; and Fjk is the tissue fraction of that voxel, estimated from the WM-, GM-, and CSF-segmented MPRAGE images, as described below. The fi corrects Cijk for in vivo (T1vivo, T2vivo assuming small patient-control differences) and in vitro (T1vitro, T2vitro) relaxation time differences:30 by using the T1vivo = 1.4, 1.3, 1.1, and 1.2 seconds and T2vivo = 343, 172, 248, and 200 ms reported for NAA, Cr, Cho, and mIns at 3T31⇓–33; T1vitro = 605, 336, 235, and 319 ms and T2vitro = 483, 288, 200, and 233 ms in the reference phantom.
Bilateral hippocampi masks were manually traced on the axial MPRAGE images on the basis of an MR imaging atlas,34 as shown in Fig 1A, and visually verified on sagittal and coronal planes. The axial MR images were segmented into CSF and gray and white matter masks by using SPM12 software (http://www.fil.ion.ucl.ac.uk/spm/software/spm12),35 as shown in Fig 2B–D. An in-house program (Matlab 14; MathWorks, Natick, Massachusetts) calculated the fraction of each tissue mask inside each voxel,36 from which Fjk = (GM volume + WM volume) / (voxel volume) for Equation 1 was obtained. We retained only voxels with at least 30% of their volume within the hippocampus mask. To minimize inclusion of degraded signals, we retained only voxels that also had the following: 1) Fjk > 70% (ie, contained <30% of CSF), 2) Cramer-Rao lower bounds of <20% for a given metabolite; and 3) 4 Hz < linewidths < 13 Hz, as shown in Fig 2. The software then estimated the global GM and WM concentrations of each metabolite in the retained n ≥ 2 hippocampal voxels by using linear regression, as described previously.36 This approach overcomes the GM/WM partial volume issue encountered by single-voxel MR spectroscopy. Note that although metabolic gradients are reported along the long axis of the hippocampus37,38 and because the entire structure is taken here as a whole and the intrasubject gradients of these metabolites are reported to be quite similar, our average concentration reported here is nevertheless a useful intersubject comparison metric.
Statistical Analyses
Exact Mann-Whitney and unequal variance t tests were used to compare the groups in terms of age, the volume measure, and each metabolite concentration without adjustment for age and sex. ANOVA was used to compare the groups in terms of the volume and each metabolite concentration adjusted for age and sex. ANOVA for each measure allowed the error variance to differ across subject groups to avoid an assumption of variance equality. F and Levene tests were used to compare the groups in terms of the variance of the volume measure and each metabolite concentration. All tests were conducted at the 5% significance level.
Results
Nineteen patients and 11 controls were included in the analyses, as shown in Table 1. The groups were not different with respect to age (P = .119). The BOLERO shim procedure yielded metabolite linewidths of 8.1 ± 3.0 Hz across all 216 VOI voxels for patients and controls, as shown in Fig 1. The mean number of voxels that passed the selection criteria described above and were used to estimate metabolite concentrations was 10.5 ± 3.4 in patients and 14.5 ± 5.4 in controls for NAA; 8.3 ± 3.4 and 11.3 ± 4.9 for Cho; 8.0 ± 3.5 and 12.8 ± 4.4 for Cr; and 4.8 ± 1.7 and 7.0 ± 3.8 for mIns. Average hippocampal GM metabolite concentrations and the volumes of the bilateral hippocampi, as shown in Fig 1A, are given in Table 2. Note that because the analysis was performed independently for each metabolite, different numbers of subjects are included for each on the basis of how many had at least 2 voxels that passed the exclusion criteria.
Patients had 19% (P <.05) higher hippocampal GM Cr levels than controls, even after adjusting for age and sex, as shown in Table 2 and Fig 3. NAA, Cho, and mIns were not significantly different. The variance of the Cho levels was larger in patients than in controls (P = .05) and trended that way for Cr (P = .06). There was a positive correlation between mIns and Cr levels in patients (r = 0.57, P = .05) but not in controls (r = 0.42, P = .26). Other metabolites were not significantly correlated. Bilateral hippocampal volume was ∼10% (P < .05) lower in patients.
Discussion
Most surprising, our findings do not support the hypothesis that patients with SZ exhibit decreased hippocampal GM NAA (due to neuronal dysfunction) compared with healthy controls, as generally reported.16 They are, however, concordant with neuropathology studies that consistently find reduced hippocampal volume without a net loss of neurons, but with neuronal architectural disarray and molecular alterations (eg, changes in receptor density and deficits in the inhibitory interneurons that compose 10% of hippocampal neurons).12
Our finding of elevated hippocampal GM Cr in SZ may suggest either altered hippocampal energy metabolism or an increased proportion of glial cells. The Cr 1H-MR spectroscopy peak is a composite of free creatine and phosphocreatine that interconvert to regenerate adenosine triphosphate, thus acting as a short-term energy “buffer” for the cell. Cr elevation may, therefore, reflect altered energy metabolism, likely hypermetabolism because a more metabolically active cell reasonably requires greater capacity to regenerate adenosine triphosphate (decreased Cr has similarly been interpreted as hypometabolism).10,39 Indeed, an increase in energy metabolism is consistent with the hippocampal hyperactivity model of SZ, supported by evidence from animal models, neuroimaging, and histopathology studies.40 Nuclear imaging studies show increased basal perfusion,41 which correlates with psychotic symptoms and is normalized by D2 antagonists,42 while resting-state fMRI shows hyperactivity correlated with cognitive deficits and negative symptoms.43 Hippocampal hyperactivity may be mediated by increased glutamatergic activity, which has been found in the hippocampus of unmedicated patients with SZ and is associated with reduced hippocampal volume and worse executive functioning and global clinical state.44 Furthermore, a reciprocal pathway between the hippocampus and dopaminergic neurons in the midbrain supports a link to the dopamine hypothesis of SZ.40
Cho and mIns are elevated alongside Cr in astrogliosis, reflecting their higher levels in glial cells than neurons.14 Although we did not find significant elevations in hippocampal GM Cho or mIns, the correlation between Cr and mIns in patients (but not controls) suggests that in at least some cases, elevation in Cr may reflect astrogliosis, a response to inflammation. Evidence supporting an immune or inflammatory component in SZ continues to accumulate. This includes the increased risk of SZ after exposure to certain maternal viruses in utero; the increased expression in the hippocampus of genes related to the immune response identified by messenger RNA sequencing; the significant associations with inflammatory-related genes identified in genome-wide association meta-analyses (including in the human leukocyte antigen complex, of which a recent study suggests that increased expression of a complement component gene may mediate excessive synaptic pruning); and the increase in proinflammatory cytokines, including interleukin-6 and tumor necrosis factor-α in patients with psychosis.45,46
That histopathology studies do not show glial cell proliferation in SZ does not rule out the possibility that a subset of patients might undergo astrogliosis.12 This notion is consistent with the emerging concept of the clinical syndrome of SZ as composed of a group of separate conditions differing in relevant pathophysiologies, termed “SZ and related psychoses.” Consistent with the idea of etiologic heterogeneity, patients exhibited greater coefficients of variation than controls for both Cho (28% versus 16%, P = .05) and Cr (25% versus 16%, P = .06). A subset of patients may undergo an inflammatory process, causing elevated Cho, Cr, and mIns through astrogliosis or subtler inflammation that damages glial cells and/or myelin.
The reduced hippocampal volume seen in patients is consistent with prior reports. Along with the prefrontal cortex, the hippocampus is consistently observed to undergo volume loss in SZ.4 A recent collaborative analysis of brain volumes by using imaging data from 2028 individuals with SZ and 2540 healthy controls found the largest patient control effect size for the hippocampus of any subcortical area, with a mean 4.10% volume decrease compared with controls.47 The lack of correlation between hippocampal volume and NAA levels in our study suggests that the lower volume in patients is not due to neuronal loss from progression of a pathogenic process and may, therefore, be suggestive of a developmental trait.
There may be several reasons for the discordance between our findings and most previous studies. First, clinical heterogeneity (disease duration, symptomatology, medications) both between and within studied populations is known to affect 1H-MR spectroscopy results.16 Second, methodology differences are well-known to affect results. Our 3D technique yields better spatial resolution (0.5 versus 3.5 cm3 or larger) and delineation of the irregular shape, with less inclusion of neighboring WM than single voxels, even accounting for partial WM and CSF volume. Third, studies that found lower NAA in patients had, on average, longer TEs, making them more sensitive to quantification errors due to T2-weighting variations, than our short TE = 35 ms acquisition.48 Fourth, some studies used Cr as an “internal reference” and reported its ratios with other metabolites. This approach magnifies errors and is further complicated by findings that Cr is altered in several brain regions in SZ.16 Finally, as others have pointed out, there is a publication bias in favor of NAA deficit reports as opposed to a finding of “no change.”16
Admittedly, our study also has several limitations. First, our sample comprised patients with varying disease durations, medication regimens, and histories of substance use at insufficient numbers to distinguish the effects of these variables. Second, a relatively small sample size may have limited the power to detect small differences (eg, Cho or mIns elevations in patients). Third, even our improved spatial resolution is insufficient to resolve hippocampal subfields that may be differentially affected in SZ.4 Fourth, the regressions used to obtain metabolite concentrations relied on relatively few, 8–14, voxels per individual, because a balance was struck among the 3D 1H-MRSI spatial resolution, measurement time, and hippocampal volume. While the number of voxels in the hippocampus may be improved at higher fields, at 3T, this is likely nearly as good as can be achieved. Fifth, because of the anatomic milieu (deep in the brain, irregularly shaped, near air-filled sinuses), it is difficult to obtain adequate SNR to quantify hippocampal glutamate and glutamine levels. Finally, our interpretation of elevated Cr is limited by the inability of 3T 1H-MR spectroscopy to distinguish Cr from phosphocreatine resonances, and our proposals of altered energy metabolism and astrogliosis are here supported by inference and existing literature rather than direct measurement.
Conclusions
Decreased volume without GM NAA decline in the hippocampus of patients with SZ suggests a process that preserves neuronal integrity and could represent a developmental rather than pathologic disruption. Increased hippocampal GM Cr suggests hypermetabolism and/or possibly astrogliosis in a subset of patients. These findings may have implications for both the pathogenesis and treatment of SZ by supporting the hippocampal hyperactivity model. In some patients, an inflammatory process (possibly involving astrogliosis) may be a treatment target in addition to the antidopaminergic treatments. Indeed, there is growing recognition that the efficacy of antipsychotics may be mediated, in part, through anti-inflammatory action, and there is ongoing effort with promising results for anti-inflammatory agents, including aspirin, COX-2 inhibitors, and N-acetylcysteine as adjuvant therapies.49 Therefore, we propose 1H-MRSI as both a probe to identify patient subgroups and a method to monitor therapies.
Footnotes
Disclosures: Emma J. Meyer—RELATED: National Institutes of Health.* Ivan I. Kirov—UNRELATED: Travel/Accommodations/Meeting Expenses Unrelated to Activities Listed: North American Brain Injury Society, Comments: for travel to their annual meeting to give a lecture on proton MR spectroscopy in traumatic brain injury. Mariana Lazar—RELATED: Grant: National Institute of Mental Health,* Comments: This work was funded, in part, by a grant from the National Institute of Mental Health. Dolores Malaspina—RELATED: Grant: National Institute of Mental Health*; UNRELATED: Expert Testimony: Provided “expert opinion” over 10,000 times, Comments: As a psychiatrist, I consult on issues such as psychiatric diagnoses of mass shooters and possible mental illness in transportation agents; Grants/Grants Pending: National Institute of Mental Health.* *Money paid to the institution.
This work was supported by National Institutes of Health grants: RC1MH088843, 2K24MH00169, MH108962, EB011639, EB01015, NS081772, NS090417, NS097494; and the Center for Advanced Imaging Innovation and Research (www.cai2r.net), a National Institute of Biomedical Imaging and Bioengineering Biomedical Technology Resource Center (NIH P41 EB017183) and CTSI UL1TR000038 (New York University). Assaf Tal acknowledges the support of the Monroy-Marks Career Development Fund, the Carolito Stiftung Fund, the Leona M. and Harry B. Helmsley Charitable Trust, and the historic generosity of the Harold Perlman Family. This research was supported by a research grant from the Sylvia Schaefer Alzheimer's Research Fund.
Indicates open access to non-subscribers at www.ajnr.org
References
- Received March 31, 2016.
- Accepted after revision June 3, 2016.
- © 2016 by American Journal of Neuroradiology