|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Neuroradiology |
1 From the Biophysics Research Institute (S.J.L., Z.L., G.W.), Division of BioStatistics (M.J.Z.), and Department of Neurology (M.F., P.G.A.), Medical College of Wisconsin, 8701 Watertown Plank Rd, Milwaukee, WI 53226. Received July 30, 2001; revision requested September 25; final revision received March 21, 2002; accepted April 12. Supported in part by the Extendicare Foundation, the Dana Foundation, and National Institutes of Health research grants DA10214, MH51358, and RR00058. Address correspondence to S.J.L. (e-mail: sjli@mcw.edu).
| ABSTRACT |
|---|
|
|
|---|
MATERIALS AND METHODS: Three subject groups (nine cognitively healthy elderly control subjects, 10 patients with probable AD, and five subjects with MCI) underwent resting-state functional magnetic resonance (MR) imaging for measurement of functional synchrony in the hippocampus. Functional synchrony was defined and quantified as the mean of the cross-correlation coefficients of spontaneous low frequency (COSLOF) components between possible pairs of voxel time courses in a brain region, or the COSLOF index. The two-tailed Student t test was used to determine differences in the COSLOF index between the control group, the probable AD group, and the MCI group. An operating characteristic curve was calculated to graphically depict the tradeoff between sensitivity and specificity of the COSLOF index.
RESULTS: Functional synchrony quantified with the COSLOF index was obtained in AD, MCI, and control subjects. COSLOF index values were significantly lower in AD patients than in control subjects (t = 4.32, P < .0012). For MCI subjects, COSLOF index values were significantly higher than those of AD patients (t = -2.4052, P < .047) but significantly lower than those of control subjects (t = 2.257, P < .043). The exponential-class curve significantly fits the relationship between the COSLOF index and the Mini-Mental Status Examination score (
2 = 20.4), indicating the rapid decrease in cognitive capacity below a threshold of the COSLOF index.
CONCLUSION: Our results suggest that the COSLOF index could be used as a noninvasive quantitative marker for the preclinical stage of AD.
© RSNA, 2002
Index terms: Alzheimer disease, 10.83 Brain, function, 10.83 Brain, MR, 10.121412, 10.121416, 10.121419 Magnetic resonance (MR), functional imaging, 10.121419
| INTRODUCTION |
|---|
|
|
|---|
Several methods have been proposed to identify signs and establish biologic markers of the preclinical phase. Measurement of plasma concentration of the 42-residue ß-amyloid has been proposed as a preclinical marker (5). Genetic markers such as the PS1 and PS2 genes are predictive for a small number of cases of early-onset AD. The ApoE4 allele, present in 50%75% of AD patients, can be used to predict when, but not if, a person is predisposed to develop AD (6,7). Hippocampal atrophy, measured with high-spatial-resolution magnetic resonance (MR) imaging of the brain, has been shown to have important diagnostic value (8,9). One disadvantage of these methods is their limited sensitivity and specificity (10).
The results of several functional neuroimaging studies with positron emission tomography (PET) and functional MR imaging suggest that compensatory functional responses may be present in asymptomatic subjects and in early AD (11,12). With disease progression, however, such a compensatory response may no longer be observable and, consequently, cannot provide reliable preclinical detection of AD because of its overlap with findings in control subjects.
Developments in functional MR imaging technology with high spatial and temporal resolutions made it possible to detect the spontaneous low-frequency fluctuation (SLF) (13,14), which has been studied for more than 40 years (15,16) and intensively reviewed (17). The SLF has also been used to map brain function (18) and to study brain lesions (19). Unlike conventional functional MR imaging, which typically requires stimulus paradigms to activate neuronal circuitry, the SLFs are obtained while a subject is at rest (rest here is relative to stimuli paradigms). We hypothesize that the gradual accumulation of AD lesions in the brain will affect the SLF signals and reduce their functional synchrony within brain regions. To determine functional synchrony, the cross-correlation coefficients between SLF of voxel time courses in a brain region are used. Since the hippocampus is one of the earliest loci affected by the accumulation of AD lesions, the functional synchrony in the hippocampus of subjects with MCI and early AD should be affected. Thus, the purpose of our study was to measure changes in functional synchrony in the hippocampus region in patients with MCI and in those with AD.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The diagnosis of probable AD was made on the basis of nationally accepted criteria (20). In addition, a score of 4 or less on the modified Hachinski Cerebral Ischemia Scale was used to exclude patients with major risk factors for vascular disease (21). Other exclusionary criteria included the presence of infarction or focal lesion on imaging studies acquired before the functional MR imaging study for all probable AD patients. Mini-Mental Status Examination (MMSE) scores were measured but were not considered as a factor for exclusion (22).
Criteria for MCI (4) consisted of the following: (a) subjective symptoms of memory loss, (b) normal activities of daily living, (c) normal general cognitive function, (d) impaired memory function for age as measured with a Rey Auditory Verbal Learning Test (<2 SDs), (e) absence of dementia, and (f) MMSE score of 24/30 or higher (22). The Clinical Dementia Rating Scale was used as an inclusion tool (23). The requirement for scale-removed sum of box scores was 0.5 or less. Subjects were required to have a modified Hachinski Cerebral Ischemia Scale score of 4 or less (21). These criteria were applied consistently to all subjects in this category.
All cognitively healthy subjects underwent a cognitive examination; none reported subjective symptoms of cognitive impairment, MMSE score of 27/30 or higher, or modified Hachinski score of 4 or less. Functional MR imaging was performed within a maximum of 2 months of cognitive testing.
MR Imaging and Functional MR Imaging
All experiments were conducted with a 1.5-T MR imager (Signa; GE Medical Systems, Milwaukee, Wis) with a three-axis balanced-torque local gradient coil and a shielded quadrature elliptic end-capped transmit or receive birdcage radio-frequency coil. Foam padding was used to limit head motion within the head coil. In all MR imaging sessions, localized T1-weighted transverse and sagittal images were acquired to provide structural information and to define the number of sections and location for functional MR imaging experiments. All experiments were performed in the sagittal plane with a gradient-recalled-echo (GRE) echo-planar MR imaging pulse sequence (2,000/40 [repetition time msec/echo time msec], in-plane resolution of 3.75 mm with 64 x 64 image matrix, field of view of 24 cm, section thickness of 7 mm). Each functional MR imaging session required 6 minutes to obtain 15 sagittal sections; each section contained 180 images. To reduce the aliasing effect with repetition time of less than 2 seconds, a crusher gradient was applied (24). Determination of the hippocampus region was conducted with a spoiled GRE sequence (24/5, flip angle of 45°, one signal acquired, section thickness of 1.1 mm). A total imaging time of 30 minutes, including functional MR imaging and spoiled GRE MR imaging, was required for each imaging session. A repeated 6-minute functional MR image was acquired in two control subjects and two AD patients for the test-retest experiments.
Functional MR Imaging and Postprocessing
All functional MR imaging data sets were subjected to a head-motion correction routine. The functional MR imaging data sets were registered among different time points with an auxiliary program in a software package (3dvolreg in Analysis of Functional Neuro-Imaging [25]). The movement parameters were checked. Subjects with excessive movement (>1 pixel) that was not correctable were excluded from further analysis. On the basis of findings in this procedure, the following subjects were excluded: four AD patients, one MCI subject, and four control subjects. Occasionally, a sharp spike in voxel time courses can be observed. Such spikes were detected and removed as follows: Assuming the mean signal intensity and its SD as µ and
, respectively, if a signal intensity is greater than µ + 4
or less than µ - 4
, and the length of these data points is two time points or less (four time points with repetition time of 2,000 msec), the intensity of the corresponding time points was then replaced with a weighted-average value of the adjacent time points of this period. All image postprocessing was performed by one author (G.W.).
Selection of Voxels in the Hippocampus and the Primary Visual Cortex
The T1-weighted MR images acquired with the spoiled GRE pulse sequence were processed with software (25). First, the software masks the entire hippocampus region, according to the brain atlas, on three-dimensional spoiled GRE images of the Talairach space and then excludes the white matter and CSF voxels (and therefore atrophic voxels), according to the intensity of the hippocampus histogram. To determine which voxel time courses should be included in the calculation of the coefficients of spontaneous low frequency (COSLOF) index, the first step is to determine which voxels are considered within the hippocampus. Since the spoiled GRE images had 1.1-mm spatial resolution and the functional MR images had 3.75-mm spatial resolution, the voxels included in the hippocampus region in functional MR images were determined according to the masked volume in the spoiled GRE images by using a deresolution program. (Note, functional MR images were also transformed into the Talairach space.) Only those voxels in functional MR images that contained more than 50% of masked spoiled GRE voxels were included for the voxel time course analysis. The same method was applied to the determination of voxels in the primary visual cortex. The voxels were selected by one author (Z.L.).
Calculation of the COSLOF Index
After the voxels in the hippocampus were selected, the corresponding voxel time courses were processed to obtain the SLF components. First, the linear trend that may be present in voxel time courses was removed by means of the linear regression method; then, the voxel time course was convolved with a Hamming bandpass filter. The expression of a Hamming filter (26) with a length of N and a passband between frequencies from fl to fh is the following:
|
To quantify the functional synchrony within the hippocampus, we conducted cross correlation between these SLF components. The cross-correlation coefficient (ccij) between any two voxel time courses si(n) and sj(n) was calculated as follows:
|
The COSLOF index defines the mean components from all possible pairs of filtered voxel time courses in a brain region. The COSLOF index is mathematically calculated from a pairwise correlation coefficient matrix:
|
j) in this matrix are symmetric, and their histogram distribution can be plotted as the cross-correlation coefficient versus its relative frequencies. A similar approach can be repeated for other brain regions such as the primary visual cortex.
Calculation of Receiver Operating Characteristic Curves
We calculated the receiver operating characteristic curves and SDs according to an asymptotic normality distribution and a simulation method (27,28), to graphically depict the trade-off between sensitivity and specificity of the COSLOF index test. When we calculated the receiver operating characteristic curves, the five MCI subjects and nine control subjects were grouped together and assumed as test negative; the 10 AD patients were assumed as test positive. The receiver operating characteristic curve is a standard summary method for displaying the accuracy of a test; it is a plot of the sensitivity versus 1 - specificity for all possible threshold values. In tests, the diseased and nondiseased states are often denoted as D+ and D-, respectively. If T denotes the diagnostic test with lower values more indicative of disease, then for possible threshold values of t, the sensitivity (true-positive fraction) and 1 - specificity (false-positive fraction) are P(T
t | D+) and P(T
t | D-), respectively. The higher the receiver operating characteristic curve, the better the test is for distinguishing subjects with disease from those without disease.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
2 = 20.4 (Fig 4). This exponential relationship between the MMSE scores and the COSLOF indexes indicates a rapid decrease in the cognitive capacity when the COSLOF index is below a threshold of 0.19. The sensitivity and specificity of the COSLOF index test are shown with a receiver operating characteristic curve analysis in Figure 5. The area under the curve of the COSLOF index test is greater than that under the 45° diagonal line; the result is consistent with the significant difference between the AD group and control group with the t test.
|
| DISCUSSION |
|---|
|
|
|---|
Developments in functional MR imaging technology allow us to observe functional synchrony in the human brain through the analysis of temporally correlated SLF components present in the functional MR imaging signal. The sensitivity and temporal resolution of functional MR imaging make it an ideal modality for assessing functional synchrony. Biswal and colleagues observed that the low-frequency resting-state fluctuations from single-section functional MR imaging time series are correlated to the left and right regions of the motor cortex (13). Other research groups have repeated these results and demonstrated consistent findings (14,37). Furthermore, SLF components are sensitive to capillary blood flow and oxygenation (38). This observation of spatial-temporal correlation provides a basis for the COSLOF index to be used to assess functional synchrony in the human brain without a focal task.
Similar to SLF components observed with blood oxygen leveldependent functional MR imaging, early reports of the measurement of brain tissue PO2 with microelectrodes (39,40) prompted Cooper et al (39) to hypothesize that these oxygenation fluctuations play an important role in maintaining an optimum balance between cerebral blood flow and cerebral metabolic rate. This was considered the earliest effort to establish the neurophysiologic foundation of SLF. Findings in the Weliky and Katz study (41), with multielectrode recordings of neuronal action potential in the lateral geniculate nucleus of awake baby ferrets, revealed patterns of spontaneous activity. It has been demonstrated that each region of the nervous system can generate its own cyclical patterns that interact with those of the other regions to which it is interconnected (42,43). The most important finding by Weliky and Katz (41) is that although there is no direct connection that links both eyes, a significant cross-correlation coefficient exists between their spontaneous bursts of neuronal activity when corticothalamic feedback is intact. It has been directly established that the spatial and temporal properties of SLF have a neuronal foundation (4446) and that synchrony of neuron firing may provide important information (47,48).
Neuropathologic changes are thought to begin in the hippocampal formation and become severe with disease progression (ie, the floor effect); thus, changes in the COSLOF index could reflect the earliest sign of AD in the region of the hippocampus. In this hypothesis, the hippocampus in the preclinical stage of AD is one of the first areas affected by neurofibrillary tangles and neuritic plaques, which in turn could affect functional synchrony, which can be measured and quantified with the COSLOF index. This characteristic of AD progression in the hippocampus region, together with the neuronal origin of SLF components and the high spatial and temporal resolutions of functional MR imaging detection, provide the foundation for the COSLOF index as a preclinical marker for the disease.
Unlike the compensatory response to functional stimuli in functional MR imaging or PET measurements, the COSLOF index has a unidirectional decrease with disease progression. In the present study, a wide range of individual COSLOF index values (from 0.05 to 0.46) and small intrasubject variation provide us with a powerful discrimination tool to characterize the different degrees of changes in functional synchrony among control subjects, MCI subjects, and AD patients. It is possible that the decrease in the COSLOF index reflects the degree of disruption in SLF synchrony in the brain of AD patients (49).
This property of the COSLOF index may fit the threshold theory for AD progression. Close inspection of Figure 4 reveals that the threshold of the COSLOF index, which separates the control group and the probable AD group, can be set around 0.19. We hypothesize that if a subject has a COSLOF index value in the hippocampus region that is less than 0.19, the subject may be at high risk of developing AD. This threshold of 0.19 may provide a quantitative evaluation for individuals with different amounts of brain reserve and with different degrees of an AD-type neuropathologic condition. Findings in numerous studies support a threshold model in which AD dementia is a result of the interplay between brain reserve and AD-type lesions (50). When the brain reserve is below a certain threshold, AD dementia will occur, and when the brain reserve remains above a certain threshold, AD dementia will not occur, even if many AD lesions are present. The COSLOF index may reflect a persons remaining brain reserve and provide a possible explanation for why a large number of individuals who meet the neuropathologic criteria for AD after death are not cognitively impaired during life (51). This concept has been further supported and demonstrated by the results of a structural MR imaging study in which patients with smaller hippocampi and accelerated loss of temporal lobe volume developed dementia (52).
By applying the concept of threshold and the disconnection models of AD, the COSLOF index provides an opportunity to objectively study MCI subjects. By measuring the COSLOF index and its rate of change in MCI subjects, we may be able to identify those indexes that occur closer to the preclinical stage of AD. In the present study, there were two control subjects and four MCI subjects with COSLOF index values below the cut-off value of 0.19 who were judged to have high risk of developing AD. findings reveal If this judgment is found to be true at clinical follow-up, then the COSLOF index could be used as a preclinical marker for AD. With such a hypothesis, future prospective studies are needed in which the COSLOF index can be used to predict progression to AD in MCI subjects.
Our methods have technical advantages over typical functional MR brain mapping techniques. Because our functional MR procedures did not require patient participation in task-driven activities for cognitive paradigms, there was no performance-related head motion. Also, because the functional MR acquisition lasted 6 minutes or less, subject tolerance (a particular concern with the elderly) was substantially improved. Furthermore, because there were no cognitive demands for probable AD patients, compliance and tolerance for the experimental procedures were much greater than we have observed in cognitive functional MR imaging studies. Nevertheless, the high failure rate (nine of 33 [27%] subjects [four AD patients, one MCI subject, four control subjects]) due to head motion with this elderly population needs to be improved in future studies by either further constraining head movement or using real-time functional MR imaging techniques and immediately repeating the functional MR acquisition if head movement occurs.
A diagnosis of definite AD requires histopathologic confirmation. In our study, we could not provide histopathologic diagnosis for AD, and therefore we could not confirm a direct correlation between the COSLOF index and neuropathologic lesions in the hippocampus region. Nevertheless, the lack of postmortem verification in our study is mitigated by the high correlation of clinical and pathologic diagnosis (92%) at our institution (53); this percentage is similar to data reported from other national specialized dementia centers.
Finally, there are technical and biologic confounding factors that may affect the determination of the COSLOF index. For instance, the magnetic susceptibility artifacts near the air-tissue interfaces may affect the SLF measurement in the hippocampus region. Extensive data processing strategies must be implemented to reduce these artifacts (54). In addition, the temporal and spatial resolutions of voxel time courses obtained with different acquisition parameters may affect the COSLOF index determination. The reduction of cardiac aliasing effect with a longer repetition time must be performed to avoid steady-state free precession disturbance (24). Other factors that may alter the COSLOF index measurements include head motion and hippocampus atrophy. In addition, the lower COSLOF index measured in the region of the hippocampus will not absolutely indicate the progression of AD processes, although AD patients had lower COSLOF index values in our study. In other words, the lower COSLOF index may be necessary but not sufficient to mark AD, because other types of pathophysiology present in the hippocampus could also affect the COSLOF index value.
In conclusion, functional MR imaging technology can be applied to AD research and generate the COSLOF index to objectively quantify functional synchrony in the human hippocampus. The lower COSLOF index in AD patients may reflect the disruption of functional synchrony by AD lesions. Therefore, the COSLOF index may be a useful tool to distinguish between probable AD patients, MCI subjects, and cognitively healthy elderly subjects. The COSLOF index measured in the human hippocampus could be used as a noninvasive marker of the cumulative effects of AD and brain reserve.
| FOOTNOTES |
|---|
Author contributions: Guarantor of integrity of entire study, S.J.L.; study concepts, S.J.L., Z.L., M.J.Z., P.G.A.; study design, S.J.L., Z.L., G.W., M.J.Z.; literature research, S.J.L., Z.L., G.W., M.J.Z.; clinical and experimental studies, Z.L.; data acquisition, S.J.L., Z.L., G.W.; data analysis/interpretation, S.J.L., Z.L., G.W., M.J.Z., P.G.A.; statistical analysis, S.J.L. G.W., M.J.Z.; manuscript preparation, S.J.L., Z.L., G.W., M.J.Z., P.G.A.; manuscript definition of intellectual content, S.J.L., M.F.; manuscript editing, S.J.L., Z.L., M.J.Z., M.F.; manuscript revision/review and final version approval, S.J.L., Z.L., G.W., M.J.Z., M.F., P.G.A.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
H. Lu, Y. Zuo, H. Gu, J. A. Waltz, W. Zhan, C. A. Scholl, W. Rea, Y. Yang, and E. A. Stein Synchronized delta oscillations correlate with the resting-state functional MRI signal PNAS, November 13, 2007; 104(46): 18265 - 18269. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Allen, H. Barnard, R. McColl, A. L. Hester, J. A. Fields, M. F. Weiner, W. K. Ringe, A. M. Lipton, M. Brooker, E. McDonald, et al. Reduced Hippocampal Functional Connectivity in Alzheimer Disease Arch Neurol, October 1, 2007; 64(10): 1482 - 1487. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Anand, Y. Li, Y. Wang, K. Gardner, and M. J. Lowe Reciprocal Effects of Antidepressant Treatment on Activity and Connectivity of the Mood Regulating Circuit: An fMRI Study J Neuropsychiatry Clin Neurosci, August 1, 2007; 19(3): 274 - 282. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |