Front. Neurosci.Frontiers in NeuroscienceFront. Neurosci.1662-453XFrontiers Media S.A.10.3389/fnins.2023.1131247NeuroscienceOriginal ResearchAltered dynamic functional connectivity in the primary visual cortex in patients with primary angle-closure glaucomaLiDong-Jin1HuangBing-Lin2*PengYuan3LiangLing-Yan4LiuHui51Health Management Center, The First People’s Hospital of Jiujiang City, Jiujiang, Jiangxi, China2College of Clinical Medicine, Jiangxi University of Traditional Chinese Medicine, Nanchang, Jiangxi, China3Department of Ophthalmology, Affiliated Hospital of Jiangxi University of Traditional Chinese Medicine, Nanchang, Jiangxi, China4College of Graduate, Jiangxi University of Traditional Chinese Medicine, Nanchang, Jiangxi, China5Department of Ophthalmology, Jiangxi Provincial People’s Hospital, Nanchang, Jiangxi, China
Edited by: Xin Huang, Renmin Hospital of Wuhan University, China
Reviewed by: Yu Ji, The First Affiliated Hospital of Nanchang University, China; Chen-Xing Qi, Renmin Hospital of Wuhan University, China; Tianming Huo, Wuhan University, China
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Purpose
Glaucoma is the main blindness-causing disease in the world. Previous neuroimaging studies demonstrated that glaucoma not only causes the loss of optic ganglion cells but also leads to the abnormal function of the optic nerve pathway and the visual cortex. However, previous studies also reported that patients with glaucoma have dysfunction in the visual cortex in a static state. Whether or not patients with primary angle-closure glaucoma (PACG) were accompanied by dynamic functional connectivity (FC) changes in the primary visual cortex (V1) remains unknown.
Methods
A total of 34 patients with PACG (23 men and 11 women) and 34 well-matched healthy controls (HCs) were enrolled in the study. The dynamic functional connectivity (dFC) with the sliding window method was applied to investigate the dynamic functional connectivity changes in the V1.
Results
Compared with HCs, patients with PACG showed increased dFC values between left V1 and bilateral calcarine (CAL). Meanwhile, patients with PACG showed increased dFC values between right V1 and bilateral CAL.
Conclusion
Our study demonstrated that patients with PACG showed increased dFC within the visual network, which might indicate the increased variability FC in the V1 in patients with PACG.
Primary angle-closure glaucoma (PACG) is a serious irreversible blinding eye disease. PACG leads to increased intraocular pressure, which in turn leads to the apoptosis of optic ganglion cells. PACG is characterized by ocular pain and visual field defects, which might lead to blindness in advanced cases. PACG is more common in Asians than Europeans and Africans, with >80% of PACG cases worldwide occurring in Asia. At present, the prevalence of PACG is 0.75% in adult Asians worldwide (Cheng et al., 2014). Glaucoma not only leads to the loss of retinal ganglion cells but also leads to abnormal changes in the optic nerve pathway and visual center function. Retinal ganglion cell death could leads to anterograde and retrograde in visual pathway (Lawlor et al., 2018).
With the development of functional magnetic resonance imaging (fMRI) technology, multimodal MRI techniques, including diffusion tensor imaging (DTI) and fMRI, have been used widely to investigate the functional changes of the visual cortex in glaucoma. Numerous studies have shown that patients with glaucoma have primary visual cortex (V1) and V1-related dysfunction. Li et al. (2017) demonstrated that patients with PACG have decreased functional connectivity (FC) between the left V1 and right V2 and increased FC between the left V1 and higher cognitive cortex.
Meanwhile, Jiang et al. (2020) reported that patients with PACG showed abnormal effective connectivity between V1 and the higher visual area, motor cortices, somatosensory cortices, and frontal lobe. Graham et al. (2021) found that dogs with unilateral PACG showed reduced fractional anisotropy (FA) in the visual pathway. Nuzzi et al. (2018) demonstrated that patients with glaucoma showed abnormal FC changes in the visual pathway. In addition, Pankowska et al. (2022) reported that the patients with advanced glaucoma showed an increase in the gray matter thickness in the V1 region. Meanwhile, early glaucoma was associated with reduced thickness in the right lateral occipital gyrus and the left lingual gyrus (Pankowska et al., 2022). Haykal et al. (2022) reported that patients with primary open-angle glaucoma (POAG) showed a loss of both axonal coherence and density in the pregeniculate visual pathways, while the postgeniculate pathways exhibited a loss of axonal coherence. Meanwhile, patients with POAG showed a decrease in fiber density and fiber-bundle cross-section in the pregeniculate optic tracts, whereas the postgeniculate optic radiation showed a decrease in fiber density (Haykal et al., 2019).
Thus, combined with our findings, previous studies demonstrated that patients with glaucoma were along with visual pathway and visual cortex dysfunction. However, the abovementioned studies focused mainly on static brain activity changes in the visual cortex. It is recognized increasingly that the functional activity changes in the brain not only are constant over time but also show dynamic changes over time.
The human brain is a complex dynamic system capable of non-stationary neural activity. The human brain showed inherently dynamic activity, which is related to the functional ability of neural networks. The dynamics of neural activity in the brain is closely related to neurophysiological activity. Growing evidence demonstrated that the fluctuations in FC are of neural origin and are also even more prominent during the resting period when mental activity is unconstrained. Temporal variability of the dynamic brain activity is related closely to various neurophysiological activities, including working memory (Shunkai et al., 2021), vision (Di and Biswal, 2020), and cognition (Patil et al., 2021). Therefore, recent studies have shown that the dynamics of brain activity can better reflect neural activity. Dynamic FC with the sliding window method is applied widely to investigate the dynamic neural activity changes in neuroimaging studies (Yuan et al., 2019; Li et al., 2020; Tsurugizawa and Yoshimaru, 2021). Previous neuroimaging studies demonstrated that patients with PACG showed abnormal neural activity changes in the visual cortex and visual-related cortex. However, there have been no studies on the dynamic FC changes in the V1 in patients with PACG. Thus, we hypothesized that patients with PACG might be associated with dynamic FC changes in the V1.
Based on these assumptions, our study is the first to determine whether or not patients with PACG are associated with the dynamic FC changes in the V1. Moreover, we chose the dFC maps as a feature of machine-learning classification. Machine learning for MRI-related neuroimaging has been widely applied in the diagnosis of nervous system diseases. In addition, the support vector machine (SVM) method is the most commonly used supervised machine learning algorithm for MRI classification related to random forests, decision trees, and convolutional neural networks.
Materials and methodsParticipants
In total, 34 patients with PACG and 34 healthy controls (matched for sex and age) were recruited in this study. The inclusion criteria for individuals with PACG were as follows: (1) the intraocular pressure is higher than 21 mmHg; (2) without other eye diseases (cataract, optic neuritis, high myopia, etc.); (3) no medical treatment of PACG; (4) eliminate patients with primary open-angle glaucoma (POAG) and normal tension glaucoma (NTG); and (5) without brain lesions including cerebral hemorrhage and cerebral infarction.
The inclusion criteria for healthy controls were as follows: (1) normal visual acuity [>1.0] in both eyes; (2) no ophthalmic diseases (i.e., optic neuritis, cataract, keratitis, etc.). This study was performed in accordance with the tenets of the Declaration of Helsinki. Each participant provided written informed consent before inclusion in the study.
MRI acquisition and experimental procedure
Magnetic resonance imaging scanning was performed on a 3-tesla magnetic resonance scanner (Discovery MR 750W system; GE Healthcare, Milwaukee, WI, USA) with the eight-channel head coil. Functional images were obtained using a gradient-echo-planar imaging sequence. fMRI scanning parameter:repetition time = 2,000 ms, echo time = 25 ms, thickness = 3.0 mm, gap = 1.2 mm, acquisition matrix = 64 × 64, field of view = 240 mm2 × 240 mm2, flip angle = 90°, voxel size = 3.6 mm3 × 3.6 mm3 × 3.6 mm3, and 35 axial slices. Before fMRI scanning, all study subjects were informed of all the caveats of the experiment and were asked to remove any metal objects. Then, during the fMRI scanning, all subjects were asked to keep their eyes closed and relaxed without falling asleep.
fMRI data analysis
All preprocessing was performed using the toolbox for Data Processing and Analysis of Brain Imaging (DPABI)1 (Yan et al., 2016), which is based on the Statistical Parametric Mapping (SPM12)2 implemented in MATLAB 2013a (MathWorks, Natick, MA, USA) and the details of steps according to a previous study (Chen et al., 2022).
dFC analysis
The dFC method was performed using DPABI software. Based on previous studies (Wen et al., 2018; Qi et al., 2021; Tong et al., 2021), left V1(x = –8, y = –76, z = 10) and right V1 (x = 7, y = –76, z = 10) with 6 mm radius. Specifically, the R-fMRI indices mentioned above were computed with hamming windows (window length = 30 TR, window step = 1 TR). To avoid the introduction of spurious fluctuations during the sliding window select, the minimum window length should be larger than 1/fmin, where fmin is the minimum frequency of the time series.
Support vector machine (SVM) analysis
The support vector machine (SVM) algorithm for binary classification is implemented on the Pattern Recognition for Neuroimaging Toolbox (PRoNTo) software3 (Schrouff et al., 2013). The details of the steps for the SVM method were as follows: (1) The dFC maps were selected as a classification feature. (2) The leave-one-out cross-validation (LOOCV) was applied to perform the SVM classifier validation. (3) The permutation test was applied to assess the statistical significance of the total accuracy of this classification. (4) The total accuracy, specificity, sensitivity, and area under the receiver operating characteristic curve (AUC) were determined to assess and classify the patients with PACG and the HCs.
Statistical analysis
The independent sample t-test was used to investigate the clinical features between the two groups. In this study, a one-sample t-test was applied to assess the group the mean of dFC maps within two groups, and the two-sample t-test was used to compare the two group differences in the dFC maps between the two groups using the Gaussian random field (GRF) method (two-tailed, voxel-level p < 0.01, GRF correction, cluster-level p < 0.05).
ResultsAnalysis of ophthalmic clinical data
There were no statistically significant differences between the PACG and HC groups in gender or age. However, there were significant differences between the two groups in visual acuity. More details of the results are shown in Table 1.
Clinical characteristics for health controls (HCs) and patients with primary angle-closure glaucoma (PACG).
PACG group
HC group
T-Value
P-Value
Sex (male/female)
23/11
23/11
N/A
0.780
Age (years)
45.15 ± 14.95
45.30 ± 13.87
−0.038
0.970
BCVA-OD
0.44 ± 0.27
1.16 ± 0.16
−11.474
<0.001*
BCVA-OS
0.43 ± 0.37
1.19 ± 0.16
−9.352
<0.001*
Chi-square test for sex. Independent t-test was used for other normally distributed continuous data.
The group means of dFC maps of the PACG and HC are shown in Figures 1A, B, 2A, B. Compared with HCs, patients with PACG showed increased dFC values between left V1 and bilateral calcarine (CAL) (Figure 1C and Table 2). The mean values of altered dFC in left V1 values were shown with a histogram (Figure 1D). Meanwhile, patients with PACG showed increased dFC values between right V1 and bilateral calcarine (CAL) (Figure 2C and Table 3). The mean values of dFC in right V1 values were shown with a histogram (Figure 2D).
One-sample t-test results of dynamic functional connectivity (dFC) in left V1 maps within primary angle-closure glaucoma (PACG) group (A) and HC group (B); significant dFC differences in left V1 between two groups (C); the mean of altered dFC values of left V1 between the patients with PACG and HCs (D). PACG, primary angle-closure glaucoma; HC, health control; V1, primary visual cortex; dFC, dynamic functional connectivity; CAL, calcarine; R, right; L, left.
Different dynamic functional connectivity (dFC) values in the left V1 between two groups.
Condition
Brain regions
BA
Peak T-scores
MNI coordinates (x, y, z)
Cluster size (voxels)
PACG > HC
Calcarine_R
–
3.5784
24, −57, 3
10
PACG > HC
Calcarine_L
–
4.0202
−6, −69, 12
51
Different dFC values in left V1 between two groups (voxel-level p < 0.01, GRF correction, cluster-level p < 0.05).
PACG, primary angle-closure glaucoma; HC, health control; MNI, Montreal Neurological Institute; dFC, dynamic functional connectivity; R, right; L, left.
One-sample t-test results of dFC in right V1 maps within PACG group (A) and HC group (B); significant dFC differences in right V1 between two groups (C); the mean of altered dFC values of right V1 between the patients with PACG and HCs (D). PACG, primary angle-closure glaucoma; HC, health control; V1, primary visual cortex; dFC, dynamic functional connectivity; CAL, calcarine; B, bilateral; R, right; L, left.
Different dFC values in the right V1 between two groups.
Condition
Brain regions
BA
Peak T-scores
MNI coordinates (x, y, z)
Cluster size (voxels)
PACG > HC
Calcarine_B
–
4.5846
21, −75, 3
237
Different dFC values in right V1 between two groups (voxel-level p < 0.01, GRF correction, cluster-level p < 0.05).
PACG, primary angle-closure glaucoma; HC, health control; MNI, Montreal Neurological Institute; dFC, dynamic functional connectivity; B, bilateral.
Support vector machine (SVM) result
For the left V1 map, the total classification accuracy was 55.88%. Function value of SVM algorithm for two groups (in scatter diagram. class 1: PACG group; class 2: HC group) (Figure 3A). A receive operating characteristic curve of the binary classifiers was generated to evaluate the system’s performance in distinguishing individuals with PACG from HCs and the AUC was 0.40 (Figure 3B).
Support vector machine classification results based on the mean dFC maps. Function dFC in left V1 value of SVM algorithm for two groups (class 1: PACG group; class 2: HC group) (A); the receive operating characteristic curve of the SVM classifier (B); function dFC in left V1 value of SVM algorithm for two groups (class 1: PACG group; class 2: HC group) (C); the receive operating characteristic curve of the SVM classifier (D); weight maps for SVM models for left V1 (E); weight maps for SVM models for right V1 (F). PACG, primary angle-closure glaucoma; HC, health control; V1, primary visual cortex; dFC, dynamic functional connectivity; SVM, support vector machine.
For the right V1 map, total classification accuracy was 51.47%. Function value of SVM algorithm for two groups (in scatter diagram. class 1: PACG group; class 2: HC group) (Figure 3C). A receive operating characteristic curve of the binary classifiers was generated to evaluate the system’s performance in distinguishing individuals with PACG from HCs and the AUC was 0.44 (Figure 3D). Weight maps for SVM models for left V1 (Figure 3E), and weight maps for SVM models for right V1 (Figure 3F).
Discussion
Our study is the first to investigate the dFC alterations of the V1 in patients with PACG. Patients with PACG showed increased dFC values between the left V1 and bilateral CAL. Meanwhile, patients with PACG showed increased dFC values between the right V1 and bilateral CAL related to the HC group.
Our most interesting finding is that the patients with PACG showed an increased dFC within the visual network. As we all know, the most important pathological mechanism of glaucoma is the loss of retinal ganglion cells. Meanwhile, retinal ganglion cell death can lead to anterograde and retrograde retinal ganglion cell degeneration (Lawlor et al., 2018). Thus, the damage to the optic nerve leads to abnormal visual signaling in the visual cortex. Thus, we speculated that the reduced afferent visual signals may lead to increased variability and decreased stability of primary visual cortex neural activity. Hernowo et al. (2011) reported that reduced volume in the visual pathway was found in patients with glaucoma. Chen et al. (2013) found that patients with glaucoma showed decreased gray matter (GM) volume in the visual cortex. Frezzotti et al. (2014) also demonstrated that patients with POAG had brain atrophy in both the visual cortex and other distant GM regions (the frontoparietal cortex, the hippocampi, and the cerebellar cortex). Wang et al. (2016) found that patients with glaucoma had significantly reduced cortical thickness in the right frontal pole and decreased GM volume in LGN, the right V1, and the left amygdala. Thus, brain structural atrophy in VI may lead to increased dFC values of V1 in patients with glaucoma. We speculated that the abnormal brain structure may lead to local neural activity flexibility changes. Consistent with these findings, our study found that patients with PACG showed increased dFC values within the visual network. Increased dFC might reflect increased neural activity flexibility in the visual cortex, which might indicate visual function compensation in patients with PACG.
In addition, in our study, the SVM method was used to investigate the predictive values of classifying individual patient populations with PACG versus HCs. For the left V1 map, the total classification accuracy was 55.88%. A receiver operating characteristic (ROC) curve of the binary classifiers was generated to evaluate the system’s performance in distinguishing individuals with PACG from HCs and the area under the ROC curve (AUC) was 0.40. For the right V1 map, the total classification accuracy was 51.47%. A ROC curve of the binary classifiers was generated to evaluate the system’s performance in distinguishing individuals with PACG from HCs and the AUC was 0.44. The SVM algorithm was trained with input data labeled previously by PACG or HCs to predict the desired outcome. This powerful multivariate analysis enabled us to make clinical predictions at the individual subject level.
However, there were some limitations to the study. First, the sample size of our study was small. Second, blood oxygen level-dependent (BOLD) signals would be influenced by a variety of physiological noises such as breathing and heartbeat. Third, the machine noise of MRI can also have an effect on the BOLD signal. Finally, a different time window length might affect the results. We would select a different window length in the future study for the stability of the results.
In conclusion, our study demonstrated that patients with PACG showed an increased dFC within the visual network, which might be indicative of an increased variability of FC in the primary visual cortex. Our study also provides an important imaging reference to aid the understanding of the mechanism of nerve damage in the visual center of patients with glaucoma.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by tenets of the Declaration of Helsinki of Jiangxi Provincial People’s Hospital. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
D-JL, B-LH, YP, L-YL, and HL contributed to the data collection, statistical analyses, and wrote the manuscript. B-LH, YP, L-YL, and HL designed the protocol and contributed to the MRI analysis and designed the study, oversaw all clinical aspects of study conduct, and manuscript preparation. All authors contributed to the article and approved the submitted version.
Funding
This study was funded by the National Natural Science Foundation of China Youth Science Fund Project (No. 821049352), the Natural Science Foundation of Jiangxi Province (No. 20202BAB2060743), and the Science and Technology Research Project of Jiangxi Provincial Department of Education (No. GJJ190647).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
ReferencesChenJ.ChenW. Y.HuangX. (2022). More than static: Altered dynamic functional connectivity of primary visual cortex in the comitant exotropia patients.Neuroreport33812–818. 10.1097/WNR.000000000000184936367793ChenW. W.WangN.CaiS.FangZ.YuM.WuQ. (2013). Structural brain abnormalities in patients with primary open-angle glaucoma: A study with 3T MR imaging.Invest. Ophthalmol. Vis. Sci.54545–554. 10.1167/iovs.12-989323258150ChengJ.ZongY.ZengY.WeiR. (2014). The prevalence of primary angle closure glaucoma in adult Asians: A systematic review and meta-analysis.PLoS One9:e103222. 10.1371/journal.pone.010322225057993DiX.BiswalB. B. (2020). Intersubject consistent dynamic connectivity during natural vision revealed by functional MRI.Neuroimage216:116698. 10.1016/j.neuroimage.2020.11669832130972FrezzottiP.GiorgioA.MotoleseI.LeucioA. D.IesterM.MotoleseE. (2014). Structural and functional brain changes beyond visual system in patients with advanced glaucoma.PLoS One9:e105931. 10.1371/journal.pone.010593125162716GrahamK. L.JohnsonP. J.BarryE. F.OrricoM. P.SoligoD. J.LawlorM. (2021). Diffusion tensor imaging of the visual pathway in dogs with primary angle-closure glaucoma.Vet. Ophthalmol.2463–74. 10.1111/vop.1282432990378HaykalS.Curcic-BlakeB.JansoniusN.CornelissenF. (2019). Fixel-Based analysis of visual pathway white matter in primary open-angle glaucoma.Invest. Ophthalmol. Vis. Sci.603803–3812. 10.1167/iovs.19-2744731504081HaykalS.InvernizziA.CarvalhoJ.JansoniusN.CornelissenF. (2022). Microstructural visual pathway white matter alterations in primary open-angle glaucoma: A neurite orientation dispersion and density imaging study.AJNR Am. J. Neuroradiol.43756–763. 10.3174/ajnr.A749535450857HernowoA. T.BoucardC. C.JansoniusN. M.HooymansJ. M.CornelissenF. W. (2011). Automated morphometry of the visual pathway in primary open-angle glaucoma.Invest. Ophthalmol. Vis. Sci.522758–2766. 10.1167/iovs.10-568221398286JiangF.FangJ.YeY.TianY.ZengX.ZhongY. (2020). Altered effective connectivity of primary visual cortex in primary angle closure glaucoma using Granger causality analysis.Acta Radiol.61508–519. 10.1177/028418511986764431390872LawlorM.Danesh-MeyerH.LevinL. A.DavagnanamI.VitaE. D.PlantG. T. (2018). Glaucoma and the brain: Trans-synaptic degeneration, structural change, and implications for neuroprotection.Surv. Ophthalmol.63296–306. 10.1016/j.survophthal.2017.09.01028986311LiS.LiP.GongH.JiangF.LiuD.CaiF. (2017). Intrinsic functional connectivity alterations of the primary visual cortex in primary angle-closure glaucoma patients before and after surgery: A resting-state fMRI study.PLoS One12:e170598. 10.1371/journal.pone.017059828122025LiY.ZhuY.NguchuB. A.WangY.WangH.QiuB. (2020). Dynamic functional connectivity reveals abnormal variability and hyper-connected pattern in autism spectrum disorder.Autism Res.13230–243. 10.1002/aur.221231614075NuzziR.DallortoL.RolleT. (2018). Changes of visual pathway and brain connectivity in glaucoma: A systematic review.Front. Neurosci.12:363. 10.3389/fnins.2018.0036329896087PankowskaA.MatwiejczukS.KoziolP.ŻarnowskiT.PieturaR.Kosior-JareckaE. (2022). Visual tract degradation in bilateral normal-tension glaucoma-cortical thickness maps and volumetric study of visual pathway areas.J. Clin. Med.11:1907. 10.3390/jcm1107190735407515PatilA. U.GhateS.MadathilD.TzengO. J.HuangH.HuangC. (2021). Static and dynamic functional connectivity supports the configuration of brain networks associated with creative cognition.Sci. Rep.11:165. 10.1038/s41598-020-80293-233420212QiC.HuangX.TongY.ShenY. (2021). Altered functional connectivity strength of primary visual cortex in subjects with diabetic retinopathy.Diabetes Metab. Syndr. Obes.143209–3219. 10.2147/DMSO.S31100934285528SchrouffJ.RosaM. J.RondinaJ. M.MarquandA. F.ChuC.AshburnerJ. (2013). PRoNTo: Pattern recognition for neuroimaging toolbox.Neuroinformatics11319–337. 10.1007/s12021-013-9178-123417655ShunkaiL.SuT.ZhongS.ChenG.ZhangY.ZhaoH. (2021). Abnormal dynamic functional connectivity of hippocampal subregions associated with working memory impairment in melancholic depression.Psychol. Med.1–13. 10.1017/S003329172100490634870570TongY.HuangX.QiC.ShenY. (2021). Altered functional connectivity of the primary visual cortex in patients with iridocyclitis and assessment of its predictive value using machine learning.Front. Immunol.12:660554. 10.3389/fimmu.2021.66055434025659TsurugizawaT.YoshimaruD. (2021). Impact of anesthesia on static and dynamic functional connectivity in mice.Neuroimage241:118413. 10.1016/j.neuroimage.2021.11841334293463WangJ.LiT.SabelB. A.ChenZ.WenH.LiJ. (2016). Structural brain alterations in primary open angle glaucoma: A 3T MRI study.Sci. Rep.6:18969. 10.1038/srep1896926743811WenZ.ZhouF.HuangX.DanH. D.XieB.ShenY. (2018). Altered functional connectivity of primary visual cortex in late blindness.Neuropsychiatr. Dis. Treat.143317–3327. 10.2147/NDT.S18375130584305YanC. G.WangX. D.ZuoX. N.ZangY. F. (2016). DPABI: Data processing & analysis for (resting-state) brain imaging.Neuroinformatics14339–351. 10.1007/s12021-016-9299-427075850YuanY.ZhangL.LiL.HuangG.AnterA.LiangZ. (2019). Distinct dynamic functional connectivity patterns of pain and touch thresholds: A resting-state fMRI study.Behav. Brain Res.375:112142. 10.1016/j.bbr.2019.11214231394144