Aberrant amplitude low-frequency fluctuation (ALFF) and regional homogeneity (ReHo) in generalized anxiety disorder (GAD) and their roles in predicting treatment remission.

BACKGROUND: Generalized anxiety disorder (GAD) is a common affective disorder characterized by comprehensive anxiety with dysregulation of brain activity which can be reflected by functional magnetic resonance imaging (f-MRI). We aimed to examine abnormal aberrant amplitude low-frequency fluctuation (ALFF) and regional homogeneity (ReHo) in GAD and evaluate their ability to predict treatment remission.
METHODS: Using resting-state fMRI (Rs-fMRI), we examined ALFF and ReHo in 30 GAD patients and 30 healthy control (HC) participants. Using on DEPASF4.3 Advanced Edition, voxel-based two-sample t-test analysis was performed on the ALFF and ReHo maps to compare GAD to HC groups, and to compare remitters (n=9) and non-remitters (n=21). Pearson's correlation analysis was used to explore the relationship between baseline Hamilton Anxiety Rating Scale (HAM-A) scores/illness duration and mean ALFF/ReHo values. The severity of GAD symptoms was rated with HAM-A. Remission was defined as HAM-A ≤7 by week 8.
RESULTS: Compared to the HC group, GAD patients showed lower ALFF in the right postcentral and right precentral gyrus; lower ReHo in the right precentral, right postcentral, and left precentral gyrus; and higher ReHo in the left posterior cingulate cortex. ALFF values for left postcentral gyrus was negatively correlated with baseline HAM-A, while that of the middle frontal gyrus was positively correlated with baseline HAM-A scores. ReHo value of the left postcentral gyrus was negatively correlated with baseline HAM-A, while that of the right middle frontal gyrus was positively correlated with baseline HAM-A scores. ALFF of the right frontal_superior_orbital and right frontal-medial-orbital cortex was positively correlated with illness duration. ReHo of the left supplementary motor area cortex was negatively correlated with illness duration. Remitters showed higher ALFF in the left hippocampus and higher ReHo value in the right postcentral cortex compared to nonremitters.
CONCLUSIONS: These results suggest that altered regional brain activity and local synchronization may be related to the pathophysiology of GAD and have certain value in predicting remission in treatment. 2020 Annals of Translational Medicine. All rights reserved.

Introduction

Of the psychiatric disorders, anxiety disorders are the most common, with its prevalence being as high as 7.3% (4.8% to 10.9%) (1). Furthermore, the prevalence of anxiety disorders varies greatly across different countries, and within countries depending on the study conducted. This variability is particularly pronounced in generalized anxiety disorder (GAD); for example, the lifetime prevalence of GAD according to the Diagnostic and Statistical Manual of Mental Disorders (DSM)-IV criteria ranges from 0.1% in Nigeria to 6.2% in New Zealand (2), and the lifetime prevalence of GAD in China ranges from 0.3% (3) to 4.1–4.6% (4). GAD is a mental disorder characterized by comprehensive anxiety, unreasonable worries, nervous and high muscle tension, and autonomic hyperactivity (5). Excessive, uncontrollable worries concerning comprehensive areas with no respective stimuli or an inappropriate sensation of potential risk is the most important diagnostic item of GAD (6). Compared to other mental diseases such as mood disorders, the diagnostic criterion of GAD is considerably ambiguous which makes the diagnosis and prognosis of GAD challenging. Thus, the exploration for, and definition of objective criterion or biomarkers in identifying disease risk and treatment response in GAD has become an urgent need for psychiatrists.

In order resolve this issue, considerable work has been done to uncover the pathophysiological mechanisms of GAD, including gene studies, protein biomarker studies, as well as functional and structural neuroimaging studies (7-10). As a noninvasive method, functional magnetic resonance imaging (fMRI) utilizes intensity changes in a type of magnetic resonance (MR) signal to track hemodynamic changes in the brain. The hemodynamic changes can be reflected with blood oxygenation level-dependent (BOLD) signals, and BOLD can be utilized as a proxy measure for neuronal function (including the degree of these activations and the synchrony of brain function across two or more regions) (11). fMRI has been applied to reflect the functional brain abnormalities in GAD, while resting-state (Rs)-fMRI has been applied to assess functional abnormalities in individuals’ resting state (12).

Regional aberrant amplitude low-frequency fluctuation (ALFF) is calculated as the averaged square root at each frequency (0.01–0.08 Hz) of the power spectrum at each voxel. It represents the spontaneous neuronal activity of the regional brain area (11). Compared to healthy controls (HC), GAD patients have been found to show higher ALFF in the bilateral dorsomedial prefrontal and dorsolateral prefrontal cortex, along with higher ALFF in the left precuneus and posterior cingulate cortex (13).

Finally, regional homogeneity (ReHo) can reflect the local synchronization of the time courses of the nearest neighboring voxels (usually 27 voxels) based on a voxel-wise measure. ReHo can thus indicate the changes of temporal neuronal activity in specific regional areas. The ReHo method has been used to investigate neural activity change in depression and GAD (11,14).

Thus far, however, few studies have combined ALFF and ReHo methods to examine GAD, and whether or not ALFF combined with ReHo can predict the treatment remission, is still unknown. We therefore conducted this study to examine the abnormal ALFF and ReHo brain regions between a GAD patient group and an HC group in order to identify the potential brain regions related to the pathophysiology of GAD and evaluate their ability to predict treatment remission. We present the following article in accordance with the MDAR reporting checklist (available at http://dx.doi.org/10.21037/atm-20-6598).

Methods

Participants

All participants were patients who met the diagnostic criteria for GAD according to DSM-IV (15), and all participants were recruited from Huzhou Third Hospital between June 2015 and June 2017. The inclusion criteria for patients were the following: (I) first episode and no treatment history; (II) age range, 18–50 years; (III) Han nationality and right-handed; (IV) education level over 9 years; (V) Hamilton Anxiety Rating Scale (HAM-A) score ≥17 and 17-item Hamilton Depression Rating Scale (HAM-D17) ≤14 (5). The exclusion criteria were as follows: (I) any other mental illness requiring medical intervention including schizophrenia, bipolar disorder, major depression disorder, etc.; (II) MR scanning contraindications; (III) structural abnormalities after routine brain MR scanning; (IV) received any psychotropic drug or psychotherapy within 2 weeks before admission. We screened 120 cases, 78 cases met the criteria of this study, and 40 GAD cases participated in this study. We recruited 32 HCs who were well matched for age and gender from the local community, and all subjects were assessed by a single psychiatrist using the Structural Clinical Interview for DSM-IV Disorders (SCID). The HAM-A scores of all HC subjects were <7. Those subjects with any psychiatric disorder history were excluded. The protocol was permitted by the Ethics Committee of Huzhou Third Hospital. The protocol was permitted by the Ethics Committee of Huzhou Third Hospital (Ethical approval number: 2018-010,036). After a complete description of the research scheme, written informed consent was obtained from all participants before the study. All procedures performed in this study involving human participants were in accordance with the Declaration of Helsinki (as revised in 2013).

Treatment

GAD patients received anti-anxiety treatment with escitalopram or venlafaxine for 8 weeks. The drug choice and change was determined by doctors based on the symptomatic characteristics and the side effects in the treatment. In all 30 GAD cases, 23 cases were treated by escitalopram while 7cases were treated by venlafaxine. The dosing of drugs was flexible and conducted according to the anxiety symptom severity and the side effects assessed by the clinician. Other anxiolytics and sedative-hypnotics were permitted.

The average dosage of escitalopram was 16.2±3.6 mg/d, and the average dosage of venlafaxine was 186.2±23.7 mg/d after 2 weeks of treatment and lasted for 8 weeks.

Clinical assessment

The severity and efficacy of the GAD symptoms was rated with HAM-A by the same trained rater throughout the treatment. HAM-A score on week 8 was applied as the primary efficacy endpoint. Remission was defined as HAM-A ≤7 after 8 weeks of treatment. GAD patients were classified as remitters if their HAM-A score ≤7; otherwise, they were classified as nonremitters. The raters were blind to the purpose of the rating scales. To maintain the blindness of this study, an experienced research coordinator was trained to manage the schedules and all data in this study.

MRI data acquisition

Rs-fMRI scan was performed using a GE Signa HDx 1.5T MRI scanner (GE Healthcare, Chicago, IL, USA) with an eight-channel head coil by a trained operator. All participants were asked to keep their eyes opened and to relax without falling asleep or thinking about specific things. To minimize the head motion in the course, a special pad was used. Resting-state MR images (BOLD) were acquired using the follow parameters: repetition time (TR)/echo time (TE) =2,000/40 ms; flip angle =83°; resolution =64×64; field of view (FOV) =240×240 mm2; voxel size =3×3×3 mm3; slice thickness =3.0 mm; interspace =1.0 mm; slices =28; duration =8 min 04 s. All subjects received a conventional high-resolution structural T1 scan for spatial normalization and localization using the follow parameters: flip angle =20°; resolution =250×250; thickness =1.2 mm; no interspace; slices =128; duration =4 min 04 s.

MRI data preprocessing

DEPASF 4.3 Advanced Edition software (http://rfmri.org/dpabi) based on Matlab 2013b was used to conduct the MRI data preprocessing and statistical analyses. First, Digital Imaging and Communications in Medicine (DICOM MRI) data files were converted into neuroimaging informatics technology initiative (NIFTI) images, and then the first 10 volumes were discarded as being the subjects’ adaptation to the scan. For slice timing and head motion correction, the remaining MRI images were slice-time–corrected and realigned to the first volume. For spatial normalization, realigned MRI images were spatially normalized according to each subject’s T1 structural images, and then all images were resampled to 3-mm isotropic voxels. Smoothing was completed using a 6-mm full-width-half-maximum (FWHM) Gaussian kernel, and the resampled images were spatially smoothed in order decrease spatial noise. The head motion of all subjects in the scanning was checked to identify whether the translation exceeded 2 mm or if the rotation exceeded 2°, if either occurred, the data were discarded for quality control. We found no significant differences of head motion and rotation between GAD patients and HC participants.

ALFF and ReHo value calculation

ALFF and ReHo analysis was performed using DEPASF4.3 Advanced Edition (http://rfmri.org/DPARSF). The ALFF value was calculated as follows: each voxel of the time series was in the power spectrum range of 0.01 to 0.08 Hz, and the average value of the square root was calculated by fast Fourier transform to calculate the ALFF; after subtracting the average value and dividing by the whole brain voxel deviation, the ALFF value was converted to z-distribution to achieve standardization; finally, the result obtained was filtered through a 0.01–0.08 Hz bandpass filter for band-pass filtering.

The ReHo value was calculated as follows: the spatially standardized data were passed through a 0.01–0.08 Hz bandpass filter to reduce low-frequency drift and high-frequency noise; Kendall’s Coefficient of Concordance (KCC) was then applied to calculate the similarity between a single voxel and the surrounding 27 voxels based on ReHo to measure the similarity of the time series in the functional cluster. Then, the individual’s ReHo image was divided by the average ReHo value of the whole brain of all the participants in this group. Finally, the ReHo brain map was spatially smoothed with a smooth core of 6 mm.

Statistical analysis

The study endpoint for all analyses was week 8. The demographic and clinical data were processed with SPSS 19.0 for Windows. Data were recorded as mean ± standard deviation (SD). All variables showed normal and near-normal distribution as checked by the Kolmogorov-Smirnov test. Student’s t-tests were used to compare the difference in demographic and clinical data. An analysis of two-sample t-test based on DEPASF 4.3 Advanced Edition Statistical Analysis (http://rfmri.org/DPARSF) was performed to identify brain areas which showed significant differences between the GAD and HC groups, and between remitters and nonremitters, with age and gender as covariates. The results were corrected with a permutation correction (voxel-level, P<0.01; cluster size, >20 voxels; number, 5,000; no acceleration; two tailed). The two-sample t-test on the voxel-based analysis were used to compare ALFF and ReHo maps between GAD and HC as well as between remitters and nonremitters (see details above). Pearson’s correlation analysis was performed to explore the relationship between baseline HAMA scores/illness duration and mean ALFF/ReHo values (see details above). A significant difference was defined as a two-tailed P value <0.05.

Results

Demographic information of the GAD vs. HC groups and remitters vs. nonremitters

In all 40 GAD patients and 32 HC who participated in the study, 10 GAD patients and 2 HC were excluded from the statistical analysis: 4 GAD patients and 2 HC for excessive head motion in quality control (translation >2 mm or rotation >2°) and 6 GAD patients for exceeding the age limit (age not matched). The final analysis consisted of 30 GAD patients (age, 18–50) and 30 HC (age, 20–48), and after 8 weeks of treatment, there were 9 remitters and 21 non-remitters in the GAD group.

There was no significant difference in the male-to-female ratio or in age between the GAD and HC groups. There was no obvious difference in gender (male-to-female ratio), age, illness duration, and education status between remitters and non-remitters. The baseline HAM-A/HAMD scores showed no significant difference, but the 8-week HAM-A/HAMD scores in the remitters were lower than those in the nonremitters. The demographic data and main clinical characteristics of all subjects are listed in .

Table 1

Comparison of demographic and clinical characteristics between GAD and HC

Characteristics	GAD (n=30)	HC (n=30)	Statistics
			χ2/t	P value
Sex (male/female)	7/23	8/22	0.089	0.766
Age (years)	40.9±10.0	37.1±6.3	1.195	0.238

GAD, generalized anxiety disorder; HC, healthy control.

Table 2

Comparison of demographic and clinical characteristics between remitters and nonremitters

Characteristics	Remitters (n=9)	Nonremitters (n=21)	Statistics
			χ2/t	P value
Sex (male/female)	3/6	4/17	1.216	0.70
Age (years)	43.92±7.37	39.40±10.93	1.512	0.140
Illness duration (months)	19.38±38.24	46.40±85.24	1.082	0.286
Education status (years)	16.8±4.2	15.9±4.5	0.5979	0.5537
Baseline HAMA	21.38±3.69	23.76±3.36	2.002	0.053
Baseline HAMD	14.15±2.64	15.48±2.12	1.68	0.102
8-week HAMA	3.38±1.45	15.16±5.60	9.906	<0.001
8-week HAMD	1.62±0.87	8.60±4.59	7.356	<0.001

ALFF and ReHo values between the GAD and HC groups

ALFF values of the GAD and HC groups

The ALFF values in GAD patients in the right postcentral gyrus, right supramarginal cortex, and right precentral gyrus were lower than those of the HC group (P<0.05, permutation corrected) ( and ).

Figure 1

ALFF comparison between GAD and HC. ALFF, amplitude low-frequency fluctuation; GAD, generalized anxiety disorder; HC, healthy control.

Table 3

Between-group differences of ALFF

Variable	MNI coordinate			Number of voxels	Peak t value
	x	y	Z
Postcentral_R	60	−3	−30	28	−4.2625
Supramarginal_R	48	−33	33	24	−3.2337
Precentral_R	48	−15	57	81	−4.4555

ALFF, amplitude low-frequency fluctuation.

ReHo values of the GAD and HC groups

The ReHo values in the GAD group in the right precentral gyrus, right postcentral gyrus, and left precentral gyrus were lower than those in the HC group. The ReHo values in the GAD group for the left posterior cingulate cortex, left angular gyrus, and left precuneus were higher than those in the HC group (P<0.05, permutation corrected) ( and ).

Figure 2

Comparison of ReHo values between the GAD and HC groups. ReHo, regional homogeneity; GAD, generalized anxiety disorder; HC, healthy control.

Table 4

Between-group differences of ReHo

Brain region	MNI coordinate			Number of voxels	Peak t value
	x	y	Z
Precentral_R	48	−15	57	793	−5.2905
Posterior Cingulate_L	−6	−30	9	56	3.3698
Angular_L	−42	−66	30	168	4.4242
Postcentral_R	30	−48	63	49	−3.1219
Precentral_L	−60	0	30	164	−5.3575
Precuneus_L	−9	−60	18	72	3.7177

Correlations between ALFF and ReHo values and clinical features (baseline HAMA and illness duration) in brain regions exhibiting significant difference

ALFF was associated with baseline HAMA

The ALFF value in the left postcentral gyrus was negatively correlated with baseline HAMA, while that of the right frontal cortex was positively correlated with baseline HAMA (P<0.05, permutation corrected) ().

Table 5

Correlations between baseline HAMA and ALFF

Brain region	MNI coordinate			Number of voxels	Peak t value
	x	y	Z
Postcentral_L	−45	−39	51	41	−0.16444
Frontal_Mid_R	27	9	48	23	0.65518

ALFF, amplitude low-frequency fluctuation.

ReHo was associated with baseline HAMA

The ReHo value in right fusiform cortex, left postcentral gyrus, and right rolandic operculum cortex was negatively correlated with baseline HAMA, while that of the right angular gyrus, right frontal-mid cortex was positively correlated with baseline HAMA (P<0.05, permutation corrected) ().

Table 6

Correlations between baseline HAMA and ReHo

Brain region	MNI coordinate			Number of voxels	Peak t value
	x	y	Z
Fusiform_R	39	−45	−21	40	−0.50909
Postcentral L	−27	−45	54	2,288	−0.7165
Rolandic Oper R	50	−15	15	50	−0.60993
Angular R	42	−72	30	62	0.61347
Frontal Mid_R	27	60	21	37	0.64469

ALFF was associated with illness duration

The ALFF value in the right frontal-superior-orbital cortex and right frontal-medial-orbital cortex was positively correlated with illness duration (P<0.05, permutation corrected) ().

Table 7

Correlations between baseline HAMA and ALFF

Brain region	MNI coordinate			Number of voxels	Peak t value
	x	y	Z
Frontal Sup Orb R	21	33	−21	24	0.62206
Frontal Med Orb R	6	30	3	69	0.83681

ALFF, amplitude low-frequency fluctuation.

ReHo was associated with illness duration

The ReHo value of the left supplementary motor area cortex was negatively correlated with illness duration, while those in the right frontal-superior-orbital cortex, left hippocampus, and frontal lobe were positively correlated with illness duration (P<0.05, permutation corrected) ().

Table 8

Correlations between baseline HAMA and ReHo

Brain region	MNI coordinate			Number of voxels	Peak t value
	x	y	Z
Frontal Sup Orb R	27	30	3	79	0.75787
Supp Motor Area L	−12	−6	66	33	−0.59128
Hippocampus L	−36	−18	−12	39	0.79572
Frontal lobe	−15	36	3	75	0.71101

ALFF and ReHo values between remitters and nonremitters

ALFF value between remitters and non-remitters

ALFF values in GAD patients who achieved remission were higher in the left fusiform cortex, left hippocampus, right thalamus, and left rolandic operculum cortex than in those who did not achieve remission (P<0.05, permutation corrected) ( and ).

Figure 3

ALFF value comparison between remitters and non-remitters. ALFF, amplitude low-frequency fluctuation.

Table 9

Between-group differences of ALFF

Brain region	MNI coordinate			Number of voxels	Peak t value
	x	y	Z
Fusiform L	−24	−78	9	320	5.0861
Hippocampus L	−27	−21	−15	964	5.6091
Thalamus R	6	−12	9	41	4.0983
Rolandic Oper L	−39	−15	18	252	4.8811

ALFF, amplitude low-frequency fluctuation.

ReHo value between remitters and non-remitters

ReHo values in GAD patients who achieved remission were higher in the left angular gyrus and right postcentral gyrus than those who did not achieve remission (P<0.05, permutation corrected) ( and ).

Figure 4

ReHo value comparison between remitters and nonremitters.

Table 10

Between-group differences of ReHo

Brain region	MNI coordinate			Number of voxels	Peak t value
	x	y	Z
Angular_L	30	−54	36	103	5.4062
Postcentral_R	15	−36	69	97	4.8186

Discussion

It is commonly accepted that ALFF can reflect regional brain activity (11,13) while ReHo can reflect local synchronization with nearby brain areas (11,14). Both ALFF and ReHo values have been used to reveal the neuropathology of psychiatric disease, including that of GAD (13,16). This study combined ALFF and ReHo values to investigate the characteristics of ALFF and ReHo in GAD patients in order to examine ALFF and ReHo as possible biomarkers for GAD, and to evaluate the ability of ALFF and ReHo values to predict treatment outcome in GAD patients. In this study, GAD patients showed several functional abnormalities in ALFF and ReHo. Compared to the HC, the GAD patients had lower ALFF values in the right postcentral gyrus and right precentral gyrus; lower ReHo values in the right precentral gyrus, right postcentral gyrus, and left precentral gyrus; and higher ReHo values in the left posterior cingulate cortex. ALFF of the left postcentral gyrus was negatively correlated with baseline HAMA, while that of the right frontal middle cortex was positively correlated with HAMA. ReHo values in the left postcentral gyrus was negatively correlated with baseline HAMA, while that of the right frontal middle cortex was positively correlated with HAMA. ALFF values in the right frontal-superior-orbital cortex and right frontal med-orbital-cortex was positively correlated with illness duration. ReHo value of the right superior-orbital cortex, left hippocampus, and frontal lobe was positively correlated with illness duration respectively. The remitters showed higher ALFF values in the left hippocampus and higher ReHo values in the right postcentral cortex compared to nonremitters.

The prevalence of GAD varies greatly depending on the country (3) and the study conducted (4,5). This is largely due to the pathological mechanism underlying GAD remaining obscure and the lack of objective diagnostic criteria. And there are no reliable markers which can predict the treatment outcome effectively, so it takes at least 2 weeks for a GAD patient to respond to anti-anxiety drug treatment while remission may take even longer (17). It is thus urgent that reliable biomarkers be found in order to guide clinicians in the diagnosis and the treatment of GAD. Accumulated evidence has proven the value of Rs-fMRI indices as biomarkers in GAD, as anxiety is the result of brain dysfunction in some brain regions (17). A core symptom of anxiety disorder, fear, (including panic) is regulated by an amygdala-centered neural circuit. The other core symptom of anxiety disorder, worry, (including anxiety pain, anticipated anxiety, and compulsion) is regulated by the cortic-striatal-thalamic-cortical (CSTC) neural circuit (18). However, anxiety disorder comprises many subtypes, with GAD being one among many, and thus it remains unclear whether these changes are specific to GAD or are instead related to anxiety in general.

In the present study, the GAD group exhibited lower ALFF and ReHo values in the right postcentral and precentral gyrus compared to the HC group, and higher ReHo values in the left posterior cingulate cortex, left angular gyrus (13), and left precuneus compared to the HC group. The GAD group also had higher ALFF values in the dorsomedial prefrontal cortex (DMPFC) and left posterior cingulate cortex compared to the HC group (19). These results indicate a hyperactivation of left posterior cingulate cortex in GAD patients. The precentral gyrus is a part of the prefrontal cortex and postcentral cortex, and is closely linked to the prefrontal cortex. The dysfunction of the prefrontal cortex in GAD has been explained in many studies (20), and emotional dysfunction due to hyperactivity of the amygdala and damage to the prefrontal cortex and hippocampal structures caused by chronic stress plays a key part in the pathological mechanism of anxiety and fear (21). Besides the abnormalities in the prefrontal cortex and left posterior cingulate cortex, we also found abnormal ReHo values in the left angular gyrus and left precuneus. Findings of studies that use rs-fMRI on GAD patients rarely report congruent results because GAD involves many brain regions and may consist of different subtypes; we thus need further large sample studies to confirm the pathological mechanism of GAD.

GAD is a chronic disease with recurrent episodes (6), and there are no reliable markers to evaluate its illness duration and the severity. The present study also found that ALFF and ReHo values in the left postcentral gyrus were negatively correlated with baseline HAMA, while those in the right frontal middle cortex were positively correlated with HAMA. This suggests that the presence of hypoactivity of the left postcentral gyrus and hyperactivity of the right frontal mid cortex may indicate more serious symptoms in GAD patients. ALFF and ReHo values in the right frontal-superior-orbital cortex were positively correlated with illness duration, suggesting that the hyperactivity of the frontal-superior-orbital cortex may indicate a longer illness duration in GAD patients. ReHo value in the left hippocampus and frontal lobe was also positively correlated with disease duration which means the dysfunction of these brain regions may also have a close correlation with the course of GAD. Andreescu et al. (18) also found that longer illness duration was positively correlated with functional connectivity (FC) between the posterior cingulate cortex and the insula in anxiety patients. Meanwhile, worry severity was negatively correlated with FC between the posterior cingulate cortex and the medial prefrontal cortex. These results stress the key importance of the postcentral cortex and prefrontal cortex in GAD. We used ALFF and ReHo values in our study while Andreescu et al. (18) used FC to evaluate brain function, which may explain the discrepancy in results between the studies. The results of rs-fMRI studies in adolescents and elderly GAD patients are usually inconsistent and difficult to duplicate (21,22), which suggests that future studies need to improve the sampling homogeneity, including that of age, course of disease, and disease severity, and this once again suggests that GAD may have different subtypes.

Antidepressants including escitalopram are the first-line drugs in the treatment of GAD (23). However, antidepressants often take weeks to take effect, and, as about half GAD patients cannot achieve remission (2), impairs patient compliance, markers which can predict treatment outcome before the treatment would be highly valuable. According to our findings, compared to nonremitters, remitters showed higher ALFF values in the left fusiform cortex, left hippocampus, right thalamus cortex, and left rolandic operculum cortex; and higher ReHo values in the left angular gyrus and right postcentral cortex. Previous research reported that elevated neural activation in the anterior cingulate cortex, hippocampus, insula, dorsolateral prefrontal cortex, amygdala, and inferior frontal gyrus predicted response to psychotherapy in panic disorder and GAD (24); naturally, as this study was focused on treatment via psychotherapy, the results are different to ours. Meanwhile, Andreescu et al. (18) found that after 12 weeks of treatment, GAD patients showed increased connectivity between the dorsolateral prefrontal cortex and prefrontal regions. Andreescu et al.’s work further highlighted the importance of the postcentral cortex in GAD, and suggested that the hyperactivity of the right postcentral cortex indicates a better outcome for patients treated with escitalopram. However, whether the baseline function of the fusiform cortex, left hippocampus, and right thalamus can predict treatment outcome requires further study.

Previous studies indicate spatiotemporal mechanisms play a key role in the psychopathological symptoms in anxiety disorders, such as dysregulation of synchronization in the neural activity of networks, and based on clinical finding in GAD and rsfMRI observations in healthy subjects, these network abnormalities in the resting state underlie core psychopathological symptoms of GAD such as unstable self, increased anxiety, increased interoceptive and somatic symptoms (25,26). So, the abnormalities of neural network activities and synchronization can be applied to assist the diagnosis of GAD. In our study, we only investigated reginal brain area but not the network of brain so the value of clinical application is limited, and in the future study, we will pay more attention on these areas.

In conclusion, the results of this study indicate that altered regional brain activity and local synchronization may be related to the pathophysiology of GAD and have certain value in predicting remission in treatment. Special attention should be paid to the precentral and postcentral gyrus, as aberrant ALFF and ReHo values in these two regions could be used as a biomarker for GAD.

There are a few limitations to our study. First, our 1.5 T scanner did not duplicate the results of abnormal amygdala function in GAD patients. Sample and scanning imprecision may account for this result, and we should improve the quality of our MRI data acquisition. Second, although we have considered the effects of depression on GAD pathophysiology and brain function, the comorbidity of GAD with major depressive disorder (MDD) is so high that we could not eliminate the effect of depression in this study, and stricter inclusion/exclusion criteria for depression, or perhaps the addition of an MDD control group, will be important for future studies. We also did not eliminate the effects of anti-anxiety drugs on brain function; in fact, brain image of f-MRI has been used to select anti-anxiety drugs in precision medicine of GAD. A larger sample size and longer study duration could also potentially reveal more significant characteristics of rs-fMRI in GAD.

The article’s supplementary files as