The Relationship between Glucose Metabolism,Resting-State fMRI BOLD Signal,and GABA A -Binding Potential: A Preliminary Study in Healthy Subjects and Those with Temporal Lobe Epilepsy

Abstract

Glucose metabolism has been associated with magnitude of blood oxygen level-dependent (BOLD) signal and connectivity across subjects within the default mode and dorsal attention networks. Similar correlations within subjects across the entire brain remain unexplored. [¹⁸F]-fluorodeoxyglucose positron emission tomography ([¹⁸F]-FDG PET), [¹¹C]-flumazenil PET, and resting-state functional magnetic resonance imaging (fMRI) scans were acquired in eight healthy individuals and nine with temporal lobe epilepsy (TLE). Regional metabolic rate of glucose (rMRGlu) was correlated with amplitude of low frequency fluctuations (ALFFs) in the fMRI signal, global fMRI connectivity (GC), regional homogeneity (ReHo), and gamma-aminobutyric acid A—binding potential (GABA_A BP_ND) across the brain. Partial correlations for ALFFs, GC, and ReHo with GABA_A BP_ND were calculated, controlling for rMRGlu. In healthy subjects, significant positive correlations were observed across the brain between rMRGlu and ALFF, ReHo and GABA_A BP_ND, and between ALFFs and GABA_A BP_ND, controlling for rMRGlu. Brain-wide correlations between rMRGlu and ALFFs were significantly lower in TLE patients, and correlations between rMRGlu and GC were significantly greater in TLE than healthy subjects. These results indicate that the glutamatergic and GABAergic systems are coupled across the healthy human brain, and that ALFF is related to glutamate use throughout the healthy human brain. TLE may be a disorder of altered long-range connectivity in association with glutamate function.

Introduction

Temporal lobe epilepsy (TLE) is a chronic neurologic disorder characterized by abnormal synchronization of populations of neurons that lead to clinical seizures arising from mesial or neocortical structures. When uncontrolled, TLE has serious adverse medical, social, and psychological consequences.^{1, 2, 3, 4} It is estimated that even with pharmacologic therapy, seizures will remain uncontrolled in up to one-third of all TLE patients, underscoring the urgent need for more effective treatments for TLE. ⁵ The etiology of TLE is poorly understood, although a large body of literature implicates both gamma-aminobutyric acid (GABA), the primary inhibitory neurotransmitter, and glutamate, the primary excitatory neurotransmitter. Although most studies support a role for GABA in the inhibition of epileptiform discharges, evidence also suggests that GABA may support spontaneous network synchronization. ⁶ Glutamatergic system changes also appear to be involved. In the healthy brain, presynaptic metabotropic glutamate receptors regulate glutamate release via a negative feedback mechanism that results in neuronal inhibition; this inhibitory control mechanism—and thus the balance between excitation and inhibition—appears to be lost in TLE. ⁷

Positron emission tomography (PET) has been used extensively to examine the physiology of seizure foci. PET studies using the ligand [¹¹C]-flumazenil, a GABA_A receptor antagonist, have consistently found decreased binding in hippocampus ipsilateral to seizure foci,^{8, 9, 10, 11} particularly in patients exhibiting atrophy or mesial temporal sclerosis. ⁸ Interpreting these findings, however, is complicated. Some evidence suggests that greater binding of flumazenil is associated with higher concentrations of extracellular GABA; higher GABA levels initiate a conformational shift in the GABA_A receptor, increasing its affinity for flumazenil. ¹² In contrast, an overabundance of GABA may downregulate GABA_A receptor—binding sites, thus reducing binding potential. ¹³ In addition, reduced [¹¹C]-flumazenil affinity for GABA_A receptors has been observed in seizure foci in TLE. ¹⁴ Regional metabolic rate of glucose (rMRGlu), measured by [¹⁸F]-fluorodeoxyglucose (FDG) PET—which is dominated by glial uptake of glucose in response to neuronal glutamate release ¹⁵ —is also abnormal in TLE. In addition to remote areas of hypometabolism, there appears to be marked hypometabolism in the ipsilateral temporal cortex. ¹⁶ Thus, in the interictal period, net activity in excitatory glutamatergic pathways may be decreased. Concomitant abnormalities in the GABAergic system are complex, and the nature of the alterations in GABA levels is not entirely understood.

In addition to studies of neurotransmitter abnormalities in TLE, a growing body of evidence suggests that TLE may be a disorder of connectivity. The default-mode network (DMN), a collection of regions activated in the absence of any imposed task or stimulus, appears to show reduced activity as well as reduced connectivity in TLE (reviewed by Cataldi et al ¹⁷ ). Other studies have shown reduced connectivity between hippocampus and some DMN areas, but increased connectivity to temporal areas as well as insular regions, often identified as part of the salience network. ¹⁸ The DMN was originally identified as the network of regions that are most metabolically active at baseline. New research in healthy controls indicates that connectivity, or the magnitude of fluctuations of the functional magnetic resonance imaging (fMRI) signal, is related to glutamate use across healthy subjects in both the DMN and dorsal attention network, ¹⁹ and that these functional networks also display the greatest degree of aerobic glycolysis. ²⁰ It is unknown whether these relationships are preserved in subjects with TLE.

The glutamatergic and GABAergic systems appear to be intimately involved in the pathophysiology of TLE, and there also appear to be widespread alterations in connectivity that may themselves be driven by altered glutamatergic function. Based on these data, we sought to examine the relationship between glutamate use (as reflected by [¹⁸F]-FDG PET), GABAergic function (as measured with [¹¹C]-flumzenil PET), connectivity, and resting-state brain activity (as measured with resting-state fMRI (rsfMRI)). To our knowledge, this study—which investigates whether fundamental alterations in these relationships may potentially underlie neurologic diseases such as TLE—is the first to examine the relationships between these systems across the entire brain in healthy subjects.

Materials and methods

Subjects

Seventeen healthy subjects and 14 subjects with TLE were recruited into the study. Written informed consent was obtained from all subjects, and the study was approved by the Combined Neuroscience Institutional Review Board and the Radiation Safety Committee at the National Institutes of Health, in compliance with the federal policy for the protection of human subjects (‘Common Rule’, Department of Health and Human Services regulation 45 Code of Federal Regulations 46). Of the subjects who participated in imaging procedures, eight healthy subjects and nine subjects with TLE had complete datasets and met all inclusion/exclusion criteria. Most withdrawals were because of either failure to obtain blood samples for quantification of PET data, identification of additional seizure foci in TLE subjects, or not meeting all inclusion/exclusion criteria.

Healthy subjects were free of personal or first-degree family history of psychiatric disorders as assessed by the structured clinical interview for Diagnostic and Statistical Manual of Mental Disorders IV, text revision (SCID ²¹ ) and the family interview for genetic studies (FIGS ²² ). TLE subjects’ seizure foci were localized by clinical neuroimaging and ictal video electroencephalography (EEG). Three had right mesial TLE and three had left mesial TLE confirmed by successful surgery. Two had left anterior temporal ictal onset, and one right anterior temporal ictal onset on video EEG. All subjects with TLE were also assessed using the SCID; only subjects not currently meeting criteria for a major depressive episode were enrolled. All subjects were administered the Hamilton Depression Rating Scale (HAM-D) at screening, and at the time of the scanning sessions, to assess depressive symptoms. Mean HAM-D score across all time points was used in subsequent analyses. Subjects with TLE were maintained on antiepileptic drugs (AEDs) throughout the study, although all exhibited seizures were not adequately controlled by adequate doses of medication. In addition, only subjects taking AEDs that had no primary effects on the GABAergic system were recruited; no subject's medication was changed for the purpose of inclusion in the study. Psychotropic medications other than AEDs were not allowed.

All healthy subjects were free of major medical or neurologic illnesses, as assessed by medical history and physical examination, blood tests, echocardiogram, and urinalysis. Subjects with TLE were free of major medical or neurologic illnesses besides epilepsy. Additional exclusion criteria for all subjects included pregnancy, breastfeeding, past drug or alcohol dependence, drug or alcohol abuse within 1 year, and contraindications to MRI or PET.

Positron Emission Tomography Acquisition and Preprocessing

All subjects underwent PET scanning using both [¹⁸F]FDG and [¹¹C]-flumazenil. In some cases, both scans were performed the same day, however, in most cases scans were acquired in separate sessions, as close as logistically possible given scheduling constraints. All brain PET scans were acquired on a CPS Innovations ECAT high-resolution research tomograph (Knoxville, TN, USA). Scans had a resolution of ∼2.5 mm full width at half maximum. Before each scanning session, subjects received two intravenous catheters for radiotracer administration and blood sampling. In the [¹¹C]-flumazenil session, subjects received an injection of 0.25-mCi/kg [¹¹C]-flumazenil, with a maximum dosage of 20 mCi for all subjects. PET scans were then acquired for up to 60 minutes. During the [¹⁸F]FDG session, subjects received 10-mCi [¹⁸F]FDG. A 35-minute cardiac emission scan was then performed on a GE Advance PET scanner (Fairfield, CT, USA). Venous blood samples were taken in 5-minute increments to acquire a cardiac input function for subsequent modeling of rMRGlu. After the cardiac scan, subjects were moved to the high-resolution research tomograph scanner and images were acquired over 40 to 60 minutes. Frames from emission scans were corrected for subject motion, and all emission scans were corrected using transmission scans. rMRGlu images were calculated using the cardiac input function as previously validated. ²³ In brief, a left ventricular time activity curve was calculated using the cardiac emission scan, and this curve was extrapolated in time using the venous blood sample data. This was then used as the input function to calculate the rMRGlu images using a simplified three-compartment model. ²⁴ Frames from the [¹¹C]-flumazenil PET images were corrected for subject motion, and a two-step simplified reference region method was then used to calculate GABA_A-binding potential images as BP_ND (according to consensus nomenclature ²⁵ ). The pons, which has little to no specific binding to flumazenil, was used as a reference region to derive the input function. All PET scans were coregistered to the anatomic MRI scans, and then transformed to Talairach ²⁶ space. To facilitate comparison between subjects, all scans were reflected left-right for subjects with a left-lateralized seizure focus so that all TLE subjects had their seizure focus in the right hemisphere.

To compare images across modalities within a single subject, a suitable mask had to be created for each subject. First, a mask was created from the spatially normalized rMRGlu images that included only voxels with a magnitude >0.8 times the mean within a simple binary brain mask derived from the Talairach ²⁶ template. This had the effect of primarily restricting the analysis to the gray matter, and is identical to the ‘threshold masking’ option implemented in statistical parametric mapping (Functional Imaging Laboratory, UCL, London, UK) widely used in PET analysis. This mask was combined with the mask created during fMRI preprocessing (described below). In addition, voxels with a value of zero in the GABA_A BP_ND image (i.e., voxels where GABA_A BP_ND was not calculated because of low signal in the raw [¹¹C]-flumazenil image) were excluded to create the final mask. We adopted this procedure in lieu of using masks derived from the MRI images to avoid problems with misclassification of voxels, and to optimize exclusion of voxels where one of our primary measures (rMRGlu) was very low and likely reflected noise rather than true uptake.

MRI Acquisition and Preprocessing

All MRI scans were acquired on a GE 3T system (Fairfield, CT, USA). Anatomic scans were acquired for coregistration using a fast spoiled gradient—echo sequence. For the resting-state scans, subjects were instructed to close their eyes, relax, and try to stay awake. Scans lasted 8 minutes, and were acquired using a gradient-echo planar imaging (EPI) sequence (repetition time (TR)=2 seconds, echo time (TE)=30 ms, resolution=64 × 64 and 34 4-mm thick sagittal slices). A respiratory belt and pulse oximeter were used during scanning to record physiologic signals. Resting-state scans were preprocessed using AFNI (Analysis of Functional NeuroImages, SSCC/NIMH/NIH, Bethesda, MD, USA). The 3dDespike routine was used to eliminate any temporal spikes, and cardiac and respiratory signals were removed using a slightly modified version of RETROICOR, as implemented in AFNI, which models heart rate, respiration, and respiration volume per time. ²⁷ This correction was performed on a slice-by-slice basis to account for changes in physiologic state during the interslice time intervals. Next, image slice timing was corrected, and EPI images were realigned to the first image. The anatomic images were warped to stereotaxic (Talairach ²⁶ ) space, and the EPI images were aligned to the fast spoiled gradient—echo anatomic images and then again transformed to Talairach space. Images were smoothed to a full width at half maximum of 10 mm, masked, and regression was performed to remove baseline effects, motion effects, and any remaining physiologic effects. All further analyses and indices were performed for the residuals from this regression analysis.

The FATCAT toolbox within AFNI was used to calculate maps of the amplitude of low frequency fluctuations (ALFFs) in the bandpass 0.01 to 0.1 Hz, as well as regional homogeneity (ReHo, a measure of local connectivity). ²⁸ The method for deriving ALFFs is described in detail by Zang et al. ²⁹ Briefly, the time series was bandpass filtered, the square root of the power spectrum was calculated, averaged over the band pass, and normalized by the global mean ALFF. The ReHo parameter is identical to Kendall's W (or Kendall's coefficient of concordance). The parameter reflects the concordance of each voxel to its neighboring 27 voxels (i.e., voxels that share a face, edge, or corner). To derive a measure of global connectivity, the mean fMRI series was normalized over time, then averaged across all voxels in the mask derived during the fMRI analysis. The time series of each voxel in the fMRI image was then correlated with this mean time series to obtain a global connectivity (GC) map for each subject. This is mathematically equivalent in calculating the correlation between each voxel and every other voxel in the brain, as shown by Saad et al. ³⁰ This measure was previously used as a measure of brain-wide connectivity (see, for example Driesen et al ³¹ ). Although weak correlations are present across the entire brain, there is significant regional heterogeneity in GC.

To determine if TLE measures differed significantly between the hippocampus contralateral and ipsilateral to the seizure focus, we used MedX (Medical Numerics, Germantown, MD, USA) to register individual MRI structural images to a template on which we had previously placed regions of interest (ROIs).^{32, 33} The ROI was an ellipsoid that fit well within the anterior hippocampus and spanned five 2-mm axial slices. For each subject, the ROI in each slice was repositioned such that only hippocampus was included within the ROI. Slices of the ROI were removed if they did not fit adequately inside the hippocampus, so that values in subjects with mesial temporal lobe sclerosis were not artefactually decreased.

Statistical Analysis

Five masked maps were obtained for each subject: rMRGlu, GABA_A BP_ND, ALFFs, GC, and ReHo. Mean values within the subject masks were compared between subject groups by independent samples t-tests. To determine the relationship between rMRGlu and rsfMRI parameters, ALFFs, GC, and ReHo were correlated with rMRGlu across all voxels (within the mask) for each subject. To determine a measure of the coupling between glutamatergic and GABAergic systems, rMRGlu and GABA_A BP_ND were similarly correlated. To assess the potential effects of GABAergic inhibition on rsfMRI parameters, partial correlations were performed between ALFFs and GABA_A BP_ND, between GC and GABA_A BP_ND, and between ReHo and GABA_A BP_ND, controlling for rMRGlu. Resulting r values were transformed to z scores using the Fisher r-to-z transformation, and assessed within groups using one sample t-tests. Owing to the large number of comparisons (seven), we applied a Bonferroni correction for multiple comparisons (requiring P<0.007 for P_corr<0.05). If a correlation differed significantly from zero, or showed a trend toward significance (P<0.014 or P_corr<0.1) after Bonferroni correction, we compared correlation coefficients between healthy and TLE groups using t-tests. Bonferroni correction would be overly conservative given that the outcome measures are correlated; thus, between-group testing on correlations showing trend-level differences from zero within groups are warranted. Finally, because we consider this study to be preliminary in nature and therefore hypothesis generating, our between-group t-test results are presented at P<0.05.

To ensure that our TLE subjects were representative of TLE subjects in the general population, we report the mean values for rMRGlu and GABA_A BP_ND in the ipsilateral and contralateral hippocampi, corrected for partial volume effects. We assessed differences between hemispheres and between healthy subjects and those with TLE using paired t-tests. Correlations are also given between rMRGlu and GABA_A BP_ND within the hippocampi.

Voxel-Wise Between-Subject Analysis

To show regional variation in correlations between PET and rsfMRI measures, and potential differences between groups, we used AFNI's 3dttest++ routine to construct a two-group linear model with rMRGlu as the outcome measure, diagnostic group as a factor, and ALFF as a covariate. The rMRGlu images were smoothed with a 10-mm Gaussian smoothing kernel for these analyses. Separate models were used to compare rMRGlu and GC, rMRGlu and ReHo, and rMRGlu and GABA BP_ND. From these models, we examined the relationship between rMRGlu and the other variables across all subjects within each diagnostic group, as well as any differences between groups. All images were masked using the intersection of the individual subject masks, dilated by one voxel in all directions. We initially set a voxel threshold of P=0.01 uncorrected, then used Monte Carlo simulations, implemented in AFNI's 3dClustSim, to correct clusters for multiple comparisons. We present results only for clusters significant at P_corr<0.05.

Results

Subjects

The healthy (N=8, mean age=34±10.6, 3 women) and TLE (N=9, age=36±9.5, 4 women) groups did not differ on the basis of age or gender. Of the nine subjects with TLE, four had right-lateralized foci, and five had left-lateralized foci. Eight subjects with TLE had depressive symptomatology rated at the time of screening and at the time of imaging; mean HAM-D score in these subjects was 5.8. All subjects with TLE were medicated at the time of the scan; five were taking levetiracetam, six lamotrigine, two lacosamide, and one subject each was taking Trileptal, trazodone, and phenytoin (note that most subjects were taking more than one AED). Mean rMRGlu, GABA_A BP_ND, ALFFs, GC, and ReHo did not differ between groups (Table 1). Mean rMRGlu, and mean rMRGlu normalized by gray matter rMRGlu, differed between ipsilateral and contralateral hippocampi (P=0.044 and P=0.038, respectively). Although GABA_A BP_ND and normalized GABA_A BP_ND were nominally reduced in the ipsilateral compared with contralateral hippocampus, these differences were not significantly different (Table 2). A difference was observed between the normalized rMRGlu and GABA_A BP_ND correlation in the healthy controls between hemispheres; in addition, the correlation differed between healthy subjects and those with TLE in the contralateral/right hippocampus. These results, however, would not survive correction for multiple comparisons, and should be interpreted with caution given the small sample size and the fact that not all correlations were themselves significant within groups.

Table 1.

Mean values for each parameter of interest

A	Group means (std)
Group	rMRGlu (mg/m/mL)	ALFFs	GC		ReHo	GABA A BP ND
Healthy controls	0.0769 (0.0278)	1482 (255.4)		0.5432 (0.0656)	0.7275 (0.0399)	3.7863 (0.5237)
TLE patients	0.0852 (0.0126)	1396 (458.4)		0.5617 (0.0854)	0.6993 (0.03774)	4.0416 (0.5312)

B	Correlations Z; mean (std)				Partial correlations Z, controlling for rMRGlu; mean (std)
Group	rMRGlu versus ALFFs a	rMRGlu versus GC a	rMRGlu versus ReHo	rMRGlu versus GABA A BP ND	GABA A BP ND versus ALFFs b	GABA A BP ND versus GC	GABA A BP ND versus ReHo
Healthy controls	0.274 (0.049) c	0.016 (0.054)	0.227 (0.053) c	0.784 (0.070) c	0.205 (0.050) c	−0.022 (0.079)	0.006 (0.056)
TLE patients	0.212 (0.061) c	0.102 (0.090) d	0.248 (0.062) c	0.760 (0.096) c	0.152 (0.068) c	−0.005 (0.064)	0.017 (0.076)

ALFFs, amplitude of low frequency functional magnetic resonance imaging fluctuations; GABA_A BP_ND, gamma-aminobutyric acid A—binding potential over all voxels in the mask for each subject group; GC, global connectivity; ReHo, regional homogeneity; rMRGlu, regional metabolic rate of glucose; TLE, temporal lobe epilepsy; std, standard deviation. Groups did not significantly differ on any measure. Mean z-transformed correlations and partial correlations between measures are also given for each group.

a

Significantly different between groups at P<0.05.

b

Trend toward a group difference at P<0.1.

c

Significantly different from 0 by one sample t-test across subjects at P<0.001.

d

Significantly different from 0 by one sample t-test across subjects at P<0.01.

Table 2.

Mean values for rMRGlu and GABA_A BP_ND in hippocampus ipsilateral and contralateral to the seizure focus, and correlations between measures

	rMRGlu mean (std)	GABA A BP ND mean (std)	Normalized rMRGlu mean (std)	Normalized GABA A BP ND mean (std)	Correlation r(P): rMRGlu versus GABA A BP ND	Correlation r(P): normalized rMRGlu versus GABA A BP ND
TLE: Ipsilateral hippocampus	0.051 (0.010) a	3.689 (2.084) b	0.529 (0.116) a	0.670 (0.309) b	0.448 (0.226)	0.704 (0.034)
TLE: Contralateral hippocampus	0.065 (0.020)	4.451 (1.450)	0.657 (0.122)	0.820 (0.205)	0.735 (0.032)	0.879 (0.002) b
HC: Left hippocampus	0.060 (0.024)	4.38 (0.868)	0.639 (0.102)	0.848 (0.090)	0.549 (0.159)	0.878 (0.004) a
HC: Right hippocampus	0.064 (0.027)	1.94 (0.701)	0.688 (0.098)	0.964 (0.080)	0.277 (0.506)	0.163 (0.700)

GABA_A BP_ND, gamma-aminobutyric acid A—binding potential; HC, healthy control; rMRGlu, regional metabolic rate of glucose; TLE, temporal lobe epilepsy; std, standard deviation.

a

Significantly different between ipsilateral and contralateral (or left and right) hippocampi at P<0.05, not corrected for multiple comparisons.

b

Significantly different between patients and controls at P<0.05, TLE ipsilateral compared with HC left, not corrected for multiple comparisons.

Correlation Results

Mean correlation coefficients, transformed to z scores using the Fisher r-to-z transformation, are given in Table 1; scatterplots for a representative subject are shown in Figure 1. Healthy subjects exhibited significant positive correlations for rMRGlu versus ALFFs, rMRGlu versus ReHo, rMRGlu versus GABA_A BP_ND, and for GABA_A BP_ND versus ALFFs controlling for rMRGlu. Subjects with TLE exhibited positive correlations for rMRGlu versus ALFFs, rMRGlu versus GC, rMRGlu versus ReHo, rMRGlu versus GABA_A BP_ND, and GABA_A BP_ND versus ALFFs controlling for rMRGlu. After Bonferroni correction for multiple comparisons, correlations between rMRGlu versus ALFFs, rMRGlu versus ReHo, rMRGlu versus GABA_A BP_ND, and GABA_A BP_ND versus ALFFs controlling for rMRGlu remained significant within both groups. A trend toward significance was also noted for rMRGlu versus GC in TLE subjects.

Figure 1.

(A–E) Regional metabolic rate of glucose (rMRGlu), amplitude of low frequency functional magnetic resonance imaging fluctuation (ALFF), global connectivity (GC), regional homogeneity (ReHo), and gamma-aminobutyric acid A—binding potential (GABA_A BP_ND), respectively, for a single subject. (F–I) Scatterplots of ALFF, GC, and GABA_A BP_ND plotted versus rMRGlu for the same subject.

Only those correlations that differed significantly from zero, or showed a trend toward significance in at least one group, were assessed between groups. Correlations between rMRGlu and ALFFs, and between rMRGlu and GC differed significantly between groups (P=0.034, Figure 2A and P=0.035, Figure 2B, respectively). In addition, a trend toward differing correlations was noted between GABA_A BP_ND and ALFFs when controlling for rMRGlu (P=0.089, Figure 2C). Although these results did not survive additional Bonferroni correction over the five correlations tested between groups, we also include effect size as additional support for the results. The between-group results for rMRGlu versus ALFFs, rMRGlu versus GC, and GABA_A BP_ND versus ALFFs had Cohen's ds of 1.12, 1.16, and 0.89, respectively, indicating large effects for all three measures.

Figure 2.

Scatterplots illustrating group differences in the z-transformed correlation coefficients for (A) global connectivity (GC) versus regional metabolic rate of glucose (rMRGlu); (B) amplitude of low frequncy functional magnetic resonance imaging fluctuation (ALFF) versus rMRGlu; and (C) partial correlation coefficients for ALFF versus gamma-aminobutyric acid A—binding potential (GABA_A BP_ND), controlling for rMRGlu. Within-group correlations significantly different from zero by one sample t-test are indicated by ∗. P values are given for group differences.

Finally, we compared the correlation between rMRGlu and GABA_A BP_ND in the hippocampi ipsilateral and contralateral to the epileptic focus. Although correlations were nominally higher in the contralateral hippocampi, this difference was not statistically significant (Table 2). As a final post-hoc analysis, we correlated mean rMRGlu, GABA_A BP_ND, ALFFs, GC, and ReHo across the entire brain with HAM-D scores and found no significant correlations; it should be noted, however, that our study was underpowered to detect such effects. In addition, no correlation was noted between contra- and ipsilateral rMRGlu and GABA_A BP_ND and HAM-D score.

Voxel-Wise Between-Subject Analyses

Across subjects, no significant correlations between rMRGlu and ALFFs, rMRGlu and GC, or rMRGlu and ReHo were observed within either group. However, in both healthy and TLE subjects, a relationship between rMRGlu and GABA_A BP_ND was observed in widespread regions. To display only the most closely related regions, the voxel-wise threshold was reduced to P=0.005, and only clusters significant at P_corr<0.05 are shown (Figures 3A and 3B; Table 3). The regions in both groups were largely similar, and no significant differences were found between groups. Partial overlap was seen in these regions with the executive control network and DMN regions.

Figure 3.

Voxelwise t-maps for the regression of glucose metaboolism (rMRGlu) and gamma-aminobutyric acid A—binding potential (GABA_A BP_ND), in (A) healthy controls and (B) subjects with temporal lobe epilepsy (TLE). Voxel-wise threshold is P<0.005 uncorrected, and only clusters surviving a P_corr<0.05 threshold are shown.

Table 3.

Regions where GABA_A BP_ND and rMRGlu were significantly correlated across subjects within each group ^a

TLE	x	y	z	Extent	t
Dorsal ACC	−1.8	−23.5	20.5	644	14.037
R Precentral	−26.2	11.5	52	361	9.84542
R Precuneus	−22.8	60.5	41.5	228	23.11272
L Precuneus	19.2	64	45	178	17.8649
L Middle temporal	47.2	39.5	3	130	9.459752
L Middle frontal	36.8	−16.5	24	104	11.83465
L Lingual gyrus	8.8	74.5	−11	98	10.64685
R Inferior frontal	−33.2	−27	6.5	92	14.09222

HC	x	y	z	Extent	t
L Postcentral	26.2	39.5	48.5	1095	20.95914
L Inferior frontal	54.2	−16.5	17	190	13.81607
Superior frontal	−1.8	−13	48.5	163	25.75655
L Middle frontal	19.2	−20	38	127	13.88949
R Superior temporal	−57.8	50	20.5	88	13.39557
R Middle frontal	−26.2	−34	24	86	10.11646

ACC, anterior cingulate cortex; GABA_A BP_ND, gamma-aminobutyric acid A—binding potential; HC, healthy control; L, left; R, right; rMRGlu, regional metabolic rate of glucose; TLE, temporal lobe epilepsy.

a

Talairach coordinates of cluster peak, cluster extent (in number of voxels), and peak t value are given.

Discussion

This preliminary study is the first to show significant group differences in the relationship between rsfMRI measures of functional connectivity and rMRGlu between healthy subjects and those with TLE. In contrast to our expectations, no differences were observed between groups in the correlation between rMRGlu and GABA_A BP_ND.

Subjects with TLE consistently showed weak but significant positive correlations between GC and rMRGlu, while healthy controls showed no such significant correlation. The finding in healthy controls would seem to contradict prior research regarding correlations between glutamatergic function and functional connectivity. For example, a magnetic resonance spectroscopy and rsfMRI study showed correlations between glutamate concentration in a medial prefrontal area and both functional and structural connectivity to subcortical structures. ³⁴ Similarly, another magnetic resonance spectroscopy and rsfMRI study showed correlations between glutamate in a posterior medial cortex ROI and connectivity within the DMN. Perhaps most convincingly, Tomasi et al ¹⁹ found positive correlations between global degree (a measure of overall connectivity) and rMRGlu in regions largely overlapping the posterior DMN in a large study of 54 healthy individuals. The most important point separating our study from these results is that prior studies examined the relationship between metabolism and connectivity across subjects, generally in focal areas, while we examined the same parameter within subjects across all gray matter. In addition to the across-subject analysis, Tomasi et al ¹⁹ also examined the relationship between connectivity and rMRGlu within subjects, but only in the DMN, dorsal attention network, and cerebellar networks, rather than examining the relationship across all gray matter, as we did. The study by Tomasi et al ¹⁹ also found no correlation between overall connectivity (global degree), and rMRGlu averaged over the entire brain; however, they did observe a correlation with local degree (a measure of local connectivity), consistent with our finding of positive correlations between rMRGlu and ReHo in healthy subjects. These results provide ample support for the notion that subjects with higher metabolic utilization in focal brain regions will also exhibit higher connectivity; however, our study also supports the idea that this relationship is confined to specific regions. Within any healthy individual's brain, our study is consistent with the idea that there is very little correlation between metabolism and global connectivity across all gray matter. This may lend support to the idea that brain connectivity is organized around ‘hub’ regions; outside the hub regions, metabolic demands may reflect more localized connectivity and distance to the nearest hub rather than long-range connectivity to distant parts of the brain.

In contrast to the healthy control subjects, subjects with TLE showed significant correlations between GC and rMRGlu. Even in the absence of clinical seizures, subjects with TLE exhibit abnormal electrophysiologic activity known as interictal discharges, which depend on both glutamatergic and GABAergic activity that may induce hypersynchronization. ³⁵ Although patients with epilepsy exhibit hypometabolism in the extended epileptogenic zone, evidence suggests that even interictal discharges cause transient focal increases in both blood flow and metabolism.^{36, 37, 38} Although purely speculative, it is conceivable that increases in metabolism after these abnormal discharges may pathologically stimulate long-range connectivity, enabling abnormal interictal discharges to propagate across the brain, resulting in seizures. Although the partial correlations we observed between GC and GABA_A BP_ND controlling for rMRGlu did not differ significantly between groups, this parameter only reflects correlations above and beyond those because of correlations between GC and rMRGlu or GABA_A BP_ND and rMRGlu. It is also worth noting that TLE subjects exhibit metabolic abnormalities beyond the seizure focus, in regions throughout the ipsilateral hemisphere,^{39, 40} consistent with the idea that long-range connectivity may be related to pathologic glucose metabolism in epileptic subjects.

We also found that healthy subjects and those with TLE differed in the coupling between ALFFs and rMRGlu. A correlation between rMRGlu and ALFFs is generally interpreted as reflecting the metabolic demands of synaptic currents and action potentials, which are reflected in the magnitude of the blood oxygen level-dependent signal. ¹⁹ Thus, the absence of such a relationship across the brains of subjects with TLE is notable. The reason for this, however, is unclear. If we assume a relatively intact neurovascular coupling, we can infer that the blood oxygen level-dependent activity measured by ALFFs accurately reflects underlying neural activity, and that the mechanisms underlying the glutamate-glutamine cycle in neurons and astrocytes are disturbed. Some evidence from animal models suggests that impaired glutamate cycling is widespread throughout the cortex, particularly implicating neuronal mitochondrial dysfunction,^{41, 42} though inadequate evidence is available in humans to firmly support this hypothesis. An alternate explanation for our finding is that neurovascular coupling is impaired in TLE. In this case, neuronal activity may not properly stimulate blood flow, resulting in decoupling from the measured blood oxygen level-dependent signal and underlying neuronal function. It appears, however, that in general neurovascular coupling is intact, although there may be a variable period of time after ictal or interictal spike events where cerebral blood flow is inadequate to meet metabolic demands (reviewed by Schwartz ⁴³ ).

A pattern in the regions of abnormality does emerge from the literature on resting-state activity and connectivity in TLE. One study that directly examined ALFFs in TLE subjects found increased ALFFs in limbic structures, particularly ipsilateral to the seizure focus, but decreased ALFFs in regions overlapping the DMN, ⁴⁴ areas that typically show the greatest rMRGlu. In addition, resting-state connectivity within the DMN, dorsal attention network, and executive control networks appears to be impaired, and regions involved in the DMN are smaller (reviewed by Cataldi et al ¹⁷ ). Strikingly, these areas also overlap regions where we found significant correlations between rMRGlu and GABA_A BP_ND in both subject groups. Thus, those areas where we showed the tightest coupling between glutamate utilization and GABA levels overlapped with those areas consistently showing abnormalities in resting-state activity and connectivity in TLE. Although we saw no differences between healthy subjects and those with TLE within these regions, our preliminary study was not adequately powered to detect such differences.

The fact that no significant differences were observed between rMRGlu and GABA_A BP_ND correlations runs counter to our original hypothesis. It may be that the GABA and glutamate balance in the human brain represents a tight window of homeostasis, deviation from which is incompatible with normal human cognition. Other studies have investigated the balance between glutamate and GABA using single-voxel magnetic resonance spectroscopy in healthy humans, and found correlations of similar magnitude; to our knowledge, however, no studies have examined this coupling across the entire brain. ⁴⁵ The positive correlation suggests that both extracellular GABA and glutamate turnover are coupled in lock step. Although it is possible that altered glutamate/GABA balance may be involved in epileptogenesis, the normal relationship appears to be largely preserved, at least during the interictal period; any potential differences are subtle in nature.

The primary limitation of this preliminary study was the small sample size. As a result, the study was underpowered to detect focal regions showing correlations across subjects or altered correlations between groups. We were therefore limited to correlative measures calculated within subjects, and our between-group comparisons did not survive correction for multiple comparisons. Notably, however, our between-group differences did show large effect sizes (D>0.8). Furthermore, we would note that a strict Bonferroni correction would not be appropriate in this study, given that the outcome measures are themselves correlated. Despite this limitation, the study proved to be a fruitful approach for detecting alterations in brain-wide function in TLE, and revealed unique aspects of normal brain function. Another important limitation was that rMRGlu, derived from [18F]-FDG PET, does not directly measure glutamate function; results should thus be interpreted with some caution. A secondary limitation was the fact that all subjects with TLE were currently medicated with AEDs. Although a medication-free study in this population group is not possible because of ethical concerns, the use of these drugs in our sample nevertheless confounds the result. It is important to note, however, that in all cases patients’ seizures were not adequately controlled by their medications. It is also possible that our results were at least partially confounded by partial volume effects, although the partial volume effects should be similar for each measure, given the high resolution of the PET scans acquired. It is also evident from the scatterplots in Figure 1 that the correlations would be preserved even if the threshold on rMRGlu for inclusion inside the analysis mask were raised. Finally, it would have been desirable to acquire all scans in a smaller time period, perhaps even the same day. The fact that we did find correlations between several of the measures seems to indicate, however, that these processes have a reasonably high degree of temporal stability.

Despite the limitations, our study provides at least preliminary evidence suggesting that metabolic-functional correlations are disrupted in TLE. We also found evidence that glutamatergic function was abnormally coupled to long-range connectivity, potentially indicating a mechanism by which discharges in the seizure focus could result in abnormal discharges in widespread regions of the brain. Further elucidation of this process could potentially lead to new targets for the treatment of TLE and other epileptic disorders.