Spatial–temporal signature of resting-state BOLD signals in classic trigeminal neuralgia

Abstract

Resting-state functional magnetic resonance imaging (R-fMRI) signals are spatiotemporally organized. R-fMRI studies in patients with classic trigeminal neuralgia (CTN) have suggested alterations in functional connectivity. However, far less attention has been given to investigations of the local oscillations and their frequency-specific changes in these patients. The objective of this study was to address this issue in patients with CTN. R-fMRI data from 17 patients with CTN and 19 age- and gender-matched healthy controls (HCs) were analyzed using amplitude of low-frequency fluctuation (ALFF). The ALFF was computed across different frequencies (slow-4: 0.027–0.073 Hz; slow-5: 0.01–0.027 Hz; and typical band: 0.01–0.08 Hz) in patients with CTN compared to HCs. In the typical band, patients with CTN showed increases of ALFF in bilateral temporal, occipital, and left middle frontal regions and in the left middle cingulate gyrus, as well as decreases of ALFF in the right inferior temporal region and in regions (medial prefrontal regions) of default mode network. These significant group differences were identified in different sub-bands, with greater brainstem findings in higher frequencies (slow-4) and extensive default mode network and right postparietal results in lower frequencies (slow-5). Furthermore, significant relationships were found between subjective pain ratings and both amplitudes of higher frequency (slow-4) blood oxygen level-dependent (BOLD) signals in pain localization brain regions and lower frequencies (slow-5) in pain signaling/modulating brain regions in the patients, and decreased ALFF within the prefrontal regions was significantly correlated with pain duration in the patients. This result supports our hypothesis that trigeminal pain has a characteristic spatiotemporal distribution of low-frequency BOLD signals. These findings might contribute to a better understanding of the impact of CTN on the brain’s intrinsic architecture. Future studies should take the frequencies into account when measuring brain resting BOLD signals of patients with CTN.

Keywords:

• frequency dependent
• resting fMRI
• slow-4
• slow-5

Introduction

Classic trigeminal neuralgia (CTN) is a unique pain disorder that is characterized by paroxysmal, electric-lancinating-like pain. It is estimated that the annual incidence is 4–5 in 100,000. It affects >360,000 people worldwide and most frequently affects individuals >50 years of age.^1,2 With the aggravating aging problem, the proportion of the affected persons with CTN is gradually increasing. Despite its prevalence and potential public burden, the changes in brain underlying CTN remain poorly understood.

Neuroimaging research has greatly improved our understanding of the impact of CTN on brain. An important proportion of these studies stands for task functional magnetic resonance imaging (fMRI), which measures stimuli-evoked brain activation.³ Previous studies of task fMRI in CTN have examined whether this syndrome stems from the central mechanism of a “pain matrix”.^3,4 However, variabilities in the activation results obtained by these task-related fMRI studies are quite large. One important reason may be that patients with CTN continuously process dull pain, therefore trigeminal neuralgia is easily concealed by these background activities, leading to greater variations in the functional activation.

During the past two decades, resting-state fMRI (R-fMRI), or task-free fMRI, has become an increasingly popular research paradigm in the human brain research community and an important modality in the human brain project.^5–7 R-fMRI does not require specific behavioral inputs and can be obtained in the task fMRI. During the R-fMRI scan, a subject is generally informed to stay awake, to “not think about special things” or “engage in day dreaming/mind wandering”, and to keep his/her eyes open or closed. Studies from R-fMRI have revealed some basic features of the human brain’s intrinsic functional architecture, of which one important facet is the spatiotemporal structure.⁸ These spatial distributions are characterized by coherent oscillations of the low-frequency band of about <0.1 Hz in resting state blood oxygen level-dependent (BOLD) data. Zang et al^9,10 have developed an index termed as amplitude of low-frequency fluctuation (ALFF) to measure the extent, or amplitude, of the low-frequency oscillations.

So far, many studies have shown that ALFF can serve as a potential biomarker in different neurologic diseases,^11–13 physiologic phenotypes,^14,15 and so on. Moreover, the brain’s oscillatory activity measured by resting BOLD signals is suitable for the investigation of alterations in spontaneous brain activity in patients with CTN. Our recent study² also found that patients with CTN showed a disrupted local homogeneity of the resting BOLD signals across the distributed brain regions. However, due to the limitations of the local coherence of resting BOLD signals (i.e., regional homogeneity [ReHo]), the measure only reflects the similarity of local adjacent voxels in time series and cannot reflect the amplitude or intensity of local spontaneous activity.¹⁶ ALFF needs no a priori and offers a valuable tool to investigate the local features of brain oscillatory activity.

In this study, we used the ALFF to investigate the amplitude of spontaneous activity and its frequency-dependent characteristics in CTN patients. We examined changes in ALFF across different frequencies (slow-4: 0.027–0.073 Hz; slow-5: 0.01–0.027 Hz; and typical band: 0.01–0.08 Hz). We hypothesized that the patterns of low-frequency brain activities would be altered in the patients, and these alterations may be frequency-specific and related to pain behaviors.

Materials and methods

Participants

Between February 2014 and March 2015, 17 preoperative outpatients with CTN (10 females and 7 males, mean age 63.41±7.25 years) from our hospital were included for study entry. The diagnoses of CTN were based on the International Classification of Headache Disorders-3 and confirmed by experienced neurologists (Yanping Wang and Xiaoling Zhang). The inclusion criteria as defined by Headache Classification Subcommittee, 2004 were: 1) unilateral pain, within the distribution of one or more branches of the trigeminal nerve; 2) stereotyped attacks, and characteristic (intense, sharp, superficial, or stabbing paroxysmal) pains precipitated from trigger areas or by trigger factors; 3) no additional neurologic or sensory deficits; and 4) no previous surgical or other invasive procedures for CTN. Additional neurologic examination was unremarkable in all patients. All patients were given carbamazepine 300–900 mg daily. In addition, 19 age- and gender-matched healthy subjects (11 females and 8 males, mean age 61.75±6.02 years) were also recruited. All of the participants were right-handed according to their self-reports.

The exclusion criteria for both groups were: 1) related neurologic conditions or chronic pains other than CTN (for the patients with CTN); 2) current chronic pains, or family history of neuralgia (for the healthy controls [HCs]); 3) history of substance abuse, systemic, psychiatric, or neurologic disorders; 4) hypertension, diabetes, a history of brain surgery; and 5) participants unable to undergo magnetic resonance imaging (MRI) scanning. The Institutional Review Board of the Second Hospital of Jiaxing City approved the study, and all participants provided written informed consent before participation. Table 1 lists the demographics of the participants.

Table 1 Demographics

	Patients with CTN (n=17)	Controls (n=19)	p-value
Age, years	62.53±7.41	61.75±6.02	0.72
Gender, M/F	7/10	8/11	>0.99
Handedness	17R	19R	>0.99
MMSE	28.9±1.7	28.1±1.3	0.105
BDI-II	2.4±1.8	2.2±1.4	0.67
HAMD	2.9±2.1	1.9±1.1	0.147
HAMA	3.5±2.9	1.7±0.9	0.035*
Duration	6.98±5.64 years, median 5.00 years	NA	NA

Notes: Data are presented as mean ± SD unless otherwise specified.

* Statistically significant differences (p<0.05).

Abbreviations: BDI-II, Beck’s Depression Inventory-II; CTN, classic trigeminal neuralgia; HAMD/HAMA, Hamilton Depression/Anxiety Rating Scale; M/F, male/female; MMSE, Mini-Mental Status Examination; NA, not applicable; R, right.

Questionnaires

Before brain imaging, we asked all patients of CTN to complete the following questionnaires based on his or her current pain. Pain severity was measured by using visual analog pain scales. This scale, widely used in neuropsychology studies, provides ratings of subjective pain on 10 grades. Depression, anxiety, and general cognitive abilities were assessed by the Beck Depression Inventory-II,¹⁷ the Hamilton Depression/Anxiety Rating Scale (HAMD/A),¹⁸ and the Mini-Mental State Examination.¹⁹ Details are shown in Table 1.

MRI image acquisition

Details of MRI data acquisition are as described in our previous work.² Briefly, T1-weighted (field of view=25.6 cm, 140 axial slices with slices thickness=1 mm, no gap) anatomic and resting-state functional (30 axial slices, TR=2 s, field of view=24 cm, 200 volumes) datasets were acquired on a 1.5 T MRI scanner (GE Medical Systems, Milwaukee, WI, USA) using a standard GE whole-head coil. Each subject was laid supine, with the head in a neutral position and fixed comfortably by a belt and foam pads during the test.

Data processing

R-fMRI data were preprocessed and analyzed using DPARSFA (http://www.restfmri.net), SPM8 (http://www.fil.ion.ucl.ac.uk/spm), and MATLAB.²⁰ The preprocessing of R-fMRI data included exclusion of the first 10 dummy scans, slice-timing correction, within-subject realignment, co-registration to the T1 structural image with diffeomorphic anatomic registration through exponentiated lie algebra,²¹ nuisance regression (white matter, cerebrospinal fluid, and the 24 Friston head motion parameters²² as covariates), spatial normalization of the functional images to the Montreal Neurological Institute space, smoothing using an isotropic 6 mm³ full-width at half maximum Gaussian kernel, and low-frequency filtering (0.01–0.08 Hz). No subject of this study showed head motion with >2.0 mm maximum displacement or 2.0° of any angular motion.

ALFF calculation

On the basis of the aforementioned preprocessing, the low-frequency filtered data were then calculated for the ALFF values of all voxels in the whole brain at each individual level. Briefly, for the four-dimensional fMRI data filtered by different low frequencies (i.e., the typical, slow-4, and slow-5), the time series were first converted to the frequency domain using Fast Fourier Transform, and then the square root of the power spectrum was calculated. The average root mean square is termed as ALFF; ALFF can reflect the absolute strength of the low-frequency oscillation signal.⁹ To improve normality of the ALFF, we also performed Fisher Z-transform.

Statistical analyses

Clinical variables were analyzed using SPSS 16 (SPSS Inc., Chicago, IL, USA). Statistical significance was set at p<0.05.

For brain metrics of ALFF, statistics were adjusted for HAMA score and completed with SPM8. First, one-sample t-test for each group and independent two-sample t-tests for intergroups in the typical frequency band (0.01–0.08 Hz) were performed. Then, for the main effects of groups (CTN and HC) and frequencies (slow-4 and slow-5) and their interactions, a 2×2 analysis of variance model was used via SPM8’s second full factorial analysis. The regions that showed significant main effects for groups were then used for a post hoc analysis. All results were corrected by Gaussian random field cluster level p<0.05, which corresponded to a voxel p<0.001 and cluster level with a p<0.05.

Results

Demographic and clinical data

There was no significant between-group difference in gen-der, age, handedness, Mini-Mental State Examination, Beck Depression Inventory or HDMA (p>0.05), except for HAMA (p<0.05). Table 1 shows the demographics of the two groups. Table 2 shows the clinical variables of the CTN patients.

Table 2 Clinical variables of the patients

Patient	Gender	Age, years	Duration, years	Location	Pain days per month	Pain intensity (VAS)	Vessel–nerve contact on MRI
1	M	46	20	V2	15	6	R
2	F	69	10	V2	10	7	L
3	F	56	3	V1+V2	10	5	R
4	M	58	6	V2+V3	8	8	L
5	M	60	1	V2	6	5	R
6	F	70	3	V2	7	7	R
7	F	68	10	V3	5	7	L
8	M	69	5	V1+V2	8	4	R
9	M	64	1	V2	11	8	L
10	F	68	2.5	V3	15	7	R
11	F	63	3	V2+V3	9	8	L
12	F	63	20	V2+V3	8	6	R
13	F	75	3	V3	6	6	R
14	M	58	10	V2	4	3	R
15	F	55	10	V3	3	7	R
16	F	65	3	V2	6	6	R
17	M	71	5	V2	7	4	R

Abbreviations: F, female; L, left; M, male; MRI, magnetic resonance imaging; R, right; V2, the second branch of the trigeminal nerve; V3, the third branch of the trigeminal nerve; VAS, visual analog scale.

ALFF in the typical band (0.01–0.08 Hz)

We started by calculating within-group ALFF results from the typical frequency band (0.01–0.08 Hz) with two-tailed one-sample t-tests. Both the CTN and HC groups exhibited a similar spatial topography, with significantly higher ALFF intensity than the whole brain mean in the default mode systems, bilateral visual cortices, bilateral middle temporal gyri, and dorsolateral prefrontal cortices. These results are compatible with the published findings and empirical data (Figure 1).

Figure 1 Within-group effects of the ALFF.

Notes: Regions of significant ALFF in the CTN and HC groups in the typical frequency band (0.01–0.08 Hz), slow-4 band (0.027–0.073 Hz), and slow-5 band (0.01–0.027 Hz) are shown separately. The effects are significant at p<0.05, FDR corrected, ≥100 contiguous voxels for one-sample t-tests. A cool color indicates that the CTN group had decreased ALFF compared with the controls and a hot color indicates the opposite. This result is rendered by using the BrainNet Viewer (https://www.nitrc.org/projects/bnv/).
Abbreviations: ALFF, amplitude of low-frequency fluctuation; CTN, classic trigeminal neuralgia; FDR, false discovery rate; HC, healthy control; L, left; R, right.

We then contrasted these ALFF patterns between the two groups. Patients with CTN showed significant ALFF increases in the left middle occipital gyrus, left middle frontal gyrus, right middle cingulate gyrus, right cerebellum (Crus 2), and bilateral temporal cortices. We also found significant ALFF decreases in the right inferior temporal and left superior medial gyri in the CTN group compared with the HC group (Figure 2; Table 3).

Figure 2 Between-group differences of the ALFF in the typical frequency band (0.01–0.08 Hz).

Notes: The effects are significant at voxel p<0.001, and cluster level p<0.05, GRF corrected for two-sample t-tests. A cool color indicates that the CTN group had decreased ALFF compared with the controls and a hot color indicates the opposite. This result is rendered by using the BrainNet Viewer (https://www.nitrc.org/projects/bnv/).
Abbreviations: ALFF, amplitude of low-frequency fluctuation; CTN, classic trigeminal neuralgia; GRF, Gaussian random field; HC, healthy control; L, left; R, right.

Table 3 In the typical frequency band, detailed information for clusters showing group ALFF differences at the given threshold (p<0.05, GRF cluster corrected)

Contrast items	Region label	Extent	z value	MNI coordinates
				x	y	z
Positive	L middle occipital gyrus	282	3.424	−21	−63	45
	L temporal pole	160	4.685	−36	18	−27
	L middle frontal gyrus	183	4.584	−42	30	33
	R MCC	176	3.285	3	−3	36
	R inferior temporal gyrus	114	4.280	54	−48	−21
	R cerebellum (Crus 2)	114	2.593	42	−42	−36
	R superior temporal gyrus	111	3.999	63	−33	12
	Cerebellar vermis (3)	94	3.167	3	−45	−6
	L middle temporal gyrus	109	2.813	−54	−33	0
Negative	R inferior temporal gyrus	113	−4.270	60	−21	−27
	L medial prefrontal/ACC	178	−4.069	−6	45	24

Notes: Regions were automatically reported and labeled using bspmview (http://www.bobspunt.com/software/bspmview) and AnatomyToolbox atlas (http://www.fz-juelich.de/inm/inm1/DE/Forschung/_docs/SPMAnatomyToolbox/SPMAnatomyToolbox_node.html). x, y, and z=MNI coordinates in the left–right, anterior–posterior, and inferior–superior dimensions, respectively.

Abbreviations: ACC, anterior cingulate cortex ALFF, amplitude of low-frequency fluctuation; GRF, Gaussian random field; L, left; MCC middle cingulate cortex; MNI, Montreal Neurological Institute; R, right.

Recent studies utilizing a R-fMRI scheme have demonstrated frequency-dependent changes of intrinsic architecture of BOLD signals in chronic back pain, which indicate different low frequencies may be distinct across different brain parcellations.^23,24 To explore the frequency-specific alterations of CTN, a special type of pain syndrome, we also divided the low-frequency range into sub-bands as previously defined: slow-5 (0.01–0.027 Hz) and slow-4 (0.027–0.073 Hz). To test for the presence of main effects and interactions, we carried out a 2×2 analysis of variance.

Brain areas displaying a significant main effect for group were found in the bilateral cerebellum, middle frontal gyri, lateral temporal regions, and the left insula (slow-4 > slow-5), as shown in Figure 3A and Table 4, as well as in the bilateral precuneus, the left superior/medial frontal gyrus, and the right fusiform (slow-4 < slow-5), as shown in Figure 3A and Table 4. Brain areas displaying a significant main effect for frequencies were found in the bilateral brainstem (slow-4 > slow-5), as shown in Figure 3B and Table 5, as well as in the bilateral superior orbital gyri, right postparietal region, and right precuneus (slow-4 < slow-5), as shown in Figure 3B and Table 5. There was no significant interaction between frequencies and groups.

Figure 3 Main effects for group (A) and frequency band (B) from the 2×2 ANOVA analysis.

Notes: Hot color represents higher ALFF in the CTN group (or slow-4 frequency band) than in the control group (or slow-5 frequency band), respectively, whereas cool color represents the opposite. This result is rendered by using the BrainNet Viewer (https://www.nitrc.org/projects/bnv/).
Abbreviations: ALFF, amplitude of low-frequency fluctuation; ANOVA, analysis of variance; CTN, classic trigeminal neuralgia; HC, healthy control; L, left; R, right.

Table 4 Main effects for groupa

Contrast items	Region label	Extent	z Value	MNI coordinates
				x	y	z
Positive (CTN>HC)	L middle frontal gyrus	246	6.861	−42	30	33
	R inferior temporal gyrus	268	6.576	54	−48	−21
	L temporal pole	154	5.762	−33	18	−24
	R superior temporal gyrus	348	4.971	63	−33	12
	R supramarginal gyrus	348	3.679	54	−33	30
	R cerebellum (Crus 2)	243	4.813	21	−84	−36
	L cerebellum (Crus 2)	243	2.672	0	−81	−30
	L superior parietal lobule	248	4.734	−24	−63	48
	L cerebellum (Crus 2)	122	4.724	−33	−84	−39
	L middle orbital gyrus	99	4.057	−33	57	−9
	L middle temporal gyrus	229	3.951	−45	−21	−9
	L insula lobe	229	3.071	−39	−12	9
	R cerebellum (VII)	108	3.624	42	−69	−51
	Cerebellar vermis (3)	107	3.345	3	−45	−9
Negative (CTN<HC)	L cerebellum (Crus 1)	424	−6.376	−9	−87	−18
	R fusiform gyrus	424	−2.714	36	−75	−12
	L superior medial gyrus	441	−5.809	−9	42	27
	L superior frontal gyrus	441	−2.572	−24	36	42
	R inferior temporal gyrus	222	−5.708	60	−24	−27
	R precuneus	263	−4.946	6	−75	48
	L precuneus	263	−4.885	0	−51	48
	R cuneus	88	−4.804	6	−93	21
	L middle occipital gyrus	123	−4.431	−30	−78	36

Notes: Regions were automatically reported and labeled using bspmview (http://www.bobspunt.com/software/bspmview) and AnatomyToolbox atlas (http://www.fz-juelich.de/inm/inm-1/EN/Forschung/_docs/SPMAnatomyToolbox/SPMAnatomyToolbox_node.html). x, y, and z=MNI coordinates in the left–right, anterior–posterior, and inferior–superior dimensions, respectively.

a Positive and negative effects at a given threshold (cluster p<0.05, GRF corrected) are shown.

Abbreviations: CTN, classic trigeminal neuralgia; GRF, Gaussian random field; HC, healthy control; L, left; MNI, Montreal Neurological Institute; R, right.

Table 5 Main effects for frequency bandsa

Contrast items	Region label	Extent	z Value	MNI coordinates
				x	y	z
Positive (S4>S5)	Brain stem	355	3.719	−6	−36	−9
	R brain stem	355	2.560	12	−30	−21
Negative (S4<S5)	R superior orbital gyrus	660	−4.498	18	27	−15
	L superior orbital gyrus	660	−4.251	−18	36	−18
	R superior parietal lobule	116	−2.791	36	−48	45
	R precuneus	116	−2.604	12	−63	54

Notes: Regions were automatically reported and labeled using bspmview (http://www.bobspunt.com/software/bspmview) and AnatomyToolbox atlas (http://www.fz-juelich.de/inm/inm-1/EN/Forschung/_docs/SPMAnatomyToolbox/SPMAnatomyToolbox_node.html). x, y, and z=MNI coordinates in the left–right, anterior–posterior, and inferior–superior dimensions, respectively.

a Positive and negative effects at a given threshold (cluster p<0.05, GRF corrected) are shown.

Abbreviations: GRF, Gaussian random field; L, left; MNI, Montreal Neurological Institute; R, right; S4, slow-4; S5, slow-5.

ALFF changes in different frequency bands (slow-4 and slow-5)

Further post hoc t-tests revealed that the group differences in ALFF in the slow-4 band were greater than those in the slow-5 band (Figure 4; Tables 6 and 7).

Figure 4 Two-sample t-tests for the slow-4 (0.027–0.073 Hz) and slow-5 (0.01–0.027 Hz) frequency bands at voxel-level p<0.001 and cluster-level p<0.05, GRF corrected.

Notes: Cool color indicates that the CTN group had decreased ALFF compared with the controls and hot color indicates the opposite. This result is rendered by using the BrainNet Viewer (https://www.nitrc.org/projects/bnv/).
Abbreviations: ALFF, amplitude of low-frequency fluctuation; CTN, classic trigeminal neuralgia; GRF, Gaussian random field; HC, healthy control; L, left; R, right.

Table 6 In the slow-4 (0.027–0.073 Hz) frequency band, detailed information for clusters showing group ALFF differences at a given threshold (cluster p<0.05, GRF corrected)

Contrast items	Region label	Extent	z Value	MNI coordinates
				x	y	z
Positive	L medial temporal pole	99	4.3745	−36	15	−30
	L middle frontal gyrus	125	4.0331	−42	30	33
	L postcentral gyrus	94	4.0083	−39	−33	39
Negative	L inferior temporal gyrus	114	−4.8141	−60	−27	−24
	R inferior temporal gyrus	134	−4.3597	60	−24	−27
	L superior medial gyrus	195	−4.1154	−6	45	24
	R superior frontal gyrus	81	−3.9903	24	66	3

Notes: Regions were automatically reported and labeled using bspmview (http://www.bobspunt.com/software/bspmview) and AnatomyToolbox atlas (http://www.fz-juelich.de/inm/inm-1/EN/Forschung/_docs/SPMAnatomyToolbox/SPMAnatomyToolbox_node.html). x, y, and z=MNI coordinates in the left–right, anterior–posterior, and inferior–superior dimensions, respectively.

Abbreviations: ALFF, amplitude of low-frequency fluctuation; GRF, Gaussian random field; L, left; MNI, Montreal Neurological Institute; R, right.

Table 7 In the slow-5 (0.01–0.027 Hz) frequency band, detailed information for clusters showing group ALFF differences at a given threshold (cluster p<0.05, GRF corrected)

Contrast items	Region label	Extent	t Value	MNI coordinates
				x	y	z
Positive (CTN>HC)	L middle frontal gyrus	199	5.489	−42	33	36
	R inferior temporal gyrus	270	5.425	54	−48	−21
	R ACC	216	4.126	3	15	27
	R cerebellum (Crus 2)	91	3.714	24	−84	−36
	R superior temporal gyrus	137	3.598	51	−36	27
Negative (CTN<HC)	R linual gyrus	146	−2.500	15	−84	−6
	L medial prefrontal cortex	210	−4.359	−9	45	15

Notes: Regions were automatically reported and labeled using bspmview (http://www.bobspunt.com/software/bspmview) and AnatomyToolbox atlas (http://www.fz-juelich.de/inm/inm-1/EN/Forschung/_docs/SPMAnatomyToolbox/SPMAnatomyToolbox_node.html). x, y, and z=MNI coordinates in the left–right, anterior–posterior, and inferior–superior dimensions, respectively.

Abbreviations: ACC, anterior cingulate cortex; ALFF, amplitude of low-frequency fluctuation; CTN, classic trigeminal neuralgia; GRF, Gaussian random field; HC, healthy control; L, left; MNI, Montreal Neurological Institute; R, right.

Correlations between ALFF and neurobehavior variables

On evaluating the relationship between ALFF and behaviors in the CTN patients, we found significantly positive relationships between subjective pain ratings and amplitudes of higher frequency (slow-4) BOLD signals in pain localization brain regions and significantly negative relationships between subjective pain ratings and amplitudes of lower frequencies (slow-5) in pain signaling/modulating brain regions in the patients. We also found that decreased ALFF within the medial/orbital prefrontal regions was negatively correlated with pain duration in the patients.

Discussion

We investigated the ALFF, a widely used index of resting fMRI, and its spatiotemporal (i.e., location- and frequency-dependent) changes in patients with CTN. In the typical band (0.01–0.08 Hz), we found that patients with CTN showed increases in ALFF in the bilateral temporal, occipital, and left middle frontal regions and in the left middle cingulate gyrus, as well as decreases of ALFF in the right inferior temporal gyrus and the medial prefrontal cortex. These significant group differences were identified in different sub-bands, and even more differences were found in sub-bands, with greater brainstem findings in slow-4 (0.027–0.073 Hz) and extensive bilateral preorbital frontal, right post-parietal, and precuneus results in slow-5 (0.01–0.027 Hz). Furthermore, significant relationships were found among subjective pain ratings and amplitudes of higher frequency (slow-4) BOLD signals in pain localization brain regions and lower frequencies (slow-5) in pain signaling/modulating brain regions in the patients. We also found that decreased ALFF within the medial/orbital prefrontal regions was significantly correlated with pain duration in the patients. This result supports our hypothesis that trigeminal pain has a characteristic spatiotemporal low-frequency BOLD signal property. These findings might contribute to a better understanding of the impact of CTN on the brain’s intrinsic architecture.

Most R-fMRI studies address functional connectivity or functional interactions between distant brain regions, while the local features that describe these intrinsic signals are also important. ALFF and ReHo belong to two analyses that describe local spontaneous brain activity. These two measures are relatively straightforward and complement each other. Many physiologic and neurologic processes can affect the low-frequency oscillatory signals, and ALFF is higher in gray matter than white matter.^6,8,9 In the resting state, the ALFF is higher across the DMN and visual cortices and has a certain hierarchical distribution in different anatomic functional partitioning units.⁸

Our results of increased ALFF in bilateral temporal, occipital, and the left middle cingulate regions, as well as decreased ALFF in the medial prefrontal cortex (also noted as the hub region of the DMN) in patients with CTN are generally consistent with our previous report² and other published findings.^25–31 One interesting result is the decreases in ALFF in the medial prefrontal cortex, which constitutes a hub region of the DMN.³² The DMN³² is the most stable and extensively studied so far and is generally believed to address spontaneous, unconstrained, and self-referential cognitive processes. Intrinsic activity in the DMN is highest during rest, and stimuli-evoked or goal-directed tasks that engage in the attentional and control resources usually reduce the activity.³² Similarly, a variety of pain processes are associated with decreases in the DMN activity, including chronic low back pain,^23,24 fibromyalgia,³³ and migraine,³⁴ and the default activity is increasingly recognized as a characteristic change of chronic pain. The reduction in DMN activity may represent the maintenance of attention and vigilance of pain and reflect CTN as a background pain sensation that interferes with the resting state of the brain that is dominated by the default network and continues to be nuisance stimulating. One illustrative example is that a more recent study found that patients with fibromyalgia who listened to their preferred relaxing music and an auditory control can increase ALFF in the DMN.³⁵

In addition to the DMN changes, we also found increased ALFF in the bilateral temporal, occipital, and left middle cingulate regions. Similar to other chronic pain syndromes, previous studies have found that morphologic and functional changes in the temporal lobe and fusiform gyrus are more prominent. For example, Obermann et al²⁹ used voxel-based morphometry analysis and found that temporal gray matter size of CTN patients was significantly reduced and correlated with the course of the disease. Desouza et al^25,26 utilized a finer morphometric analysis and found that cortical thickness of patients with CTN has changed significantly. Cortical thickness changes occur across the areas that process noxious stimulation, such as orbital frontal cortex, insula, sensorimotor, (anterior) cingulate, and temporal regions. May³⁶ summarized the voxel-based morphometry results of 15 chronic pains from 30 published papers using meta-analysis and showed that changes in the cortices (increased or decreased gray matter) occur in different pain syndromes but in some common regions such as the anterior cingulate, temporal lobe, prefrontal cortex, thalamus, motor cortex, brainstem, dorsal lateral prefrontal cortex, and insula. In another study that included both patients with trigeminal neuralgia and trigeminal neuropathy, Gustin et al³⁷ reported a reduction of gray matter size in the sensorimotor and basal areas. In a report of migraine headache, it was also found that the volume of the fusiform gyrus was reduced.³⁸ The temporal and occipital regions are involved in multidimensional sensory processing, which may reflect the CTN pathophysiology related to including continuous processing of nociceptive input information and large-scale functional brain network communication.

Like many other chronic pains, CTN is a unique maladaptive process and affects the brain function and behavior through adjusting the cortical oscillatory rhythms both locally and globally. Baliki et al²³ established the link between chronic back pain and its frequency-dependent architecture of resting BOLD signals; they found the signals in patients showed a shift to high-frequency dynamics, with increases focusing on the medial prefrontal cortex and the DMN.

Previous studies have also suggested the role of precuneus in chronic pain and the pathologic changes of CTN. Parise et al³⁰ have reported cortical thinning of the precuneus and inferior temporal/fusiform gyrus in CTN patients. Although the cuneus and inferior temporal gyrus are often associated with visual information processing, they are limited by the area of the superior colliculus – a multidimensional information integration. One of the presumed functions of the cuneus seems to be the integration of somatosensory information and other sensory stimulation and cognitive processes such as attention, learning, and memory. These regions have also been reported to be associated with sensory motor zones and other pain-related areas, such as the trigeminal nerve distribution region and the distribution area that are induced by thermal pain stimulation and Aδ fibrosis selective stimulation. In healthy subjects, electric shock pain also activates a series of cortical regions, including the left fusiform, hippocampus, primary and secondary visual cortex, extending to the precuneus and a more dopamine pathway network. The dorsal attention network, activated in the shock stimulus, may involve the process of mental imagination related to the extraction of electric shock perception.³⁹ The pain of CTN is essentially similar to the nature of the electric shock and, therefore, may also lead to similar activation location. Desouza et al^25,26 also found that the lateral frontal lobe and the orbital frontal cortex of CTN patients were also thinner, and the loss of the cortical gray matter could also explain the reduction in the spontaneous activity of the medial frontal lobe in this study. In addition, the changes in the brain activity of the fusiform gyrus and precuneus can be activated by trigeminal nerve Aδ fibers and may be involved in the interpretation of the triggering of trigeminal mental activity.⁴⁰

There are several limitations in this study. The role of carbamazepine on intrinsic brain activity is largely unknown, but we think this effect would be limited. Our sample size was relatively small, which limited the statistical power. Future studies using a larger sample size to increase the test–retest reliability are required.

Conclusion

Results from resting-state oscillations in typical frequency bands showed increases in ALFF in sensorimotor, multidimensional sensory processing, and interoceptive systems, as well as decreases in the default system, suggesting that the pathophysiologic basis of CTN involves the inclusion of sustained processing of noxious inputs and disruptions in modular segregation. The significant group differences identified in different sub-bands with greater brainstem findings in slow-4 (0.027–0.073 Hz) and extensive bilateral preorbital frontal, right posterior parietal, and precuneus results in slow-5 (0.01–0.027 Hz) support our hypothesis that trigeminal pain has a characteristic spatiotemporal resting state in the fMRI low-frequency BOLD signal properties. These findings might contribute to a better understanding of the impact of CTN on the brain’s intrinsic architecture. Future research should take the frequencies into account when investigating the brain resting BOLD signals of patients with CTN.

Acknowledgments

This work was supported by grants from the Science and Technology Plan Projects of Jiaxing (grant no. 2016AY23049).

Disclosure

The authors report no conflicts of interest in this work.