ORIGINAL RESEARCH article
Reduced Global-Brain Functional Connectivity of the Cerebello-Thalamo-Cortical Network in Patients With Dry Eye Disease
- 1National Clinical Research Center for Mental Disorders, and Department of Psychiatry, The Second Xiangya Hospital of Central South University, Changsha, China
- 2Department of Neurology, The First Affiliated Hospital of Guangxi Medical University, Nanning, China
- 3Department of Radiology, Tianjin Medical University General Hospital, Tianjin, China
- 4Department of Radiology, The Second Xiangya Hospital of Central South University, Changsha, China
- 5Department of Radiology, The First Affiliated Hospital of Guangxi Medical University, Nanning, China
- 6The Third People's Hospital of Foshan, Foshan, China
Background: The pathophysiology of patients with dry eye disease (DED) is associated with abnormal functional connectivity (FC). The present study aims to probe alterations of voxel-wise brain-wide FC in patient with DED at rest in an unbiased way.
Method: A total of 20 patients with DED and 23 controls matched by age, sex, and years of education underwent resting-state functional magnetic resonance imaging scans. Global-brain FC (GFC) was adopted to analyze the images. Support vector machine (SVM) was utilized to differentiate the patients from the controls.
Results: Compared with the controls, patients with DED exhibited decreased GFC in the right cerebellum lobule VIII/inferior semi-lunar lobule and left thalamus that belonged to the cerebello-thalamo-cortical network. The GFC values in the left thalamus were positively correlated to the illness duration (r = 0.589, p = 0.006) in the patients. Decreased GFC values in the left thalamus could be used to discriminate the patients from the controls with optimal accuracy, sensitivity and specificity (88.37, 85.00, and 91.30%).
Conclusions: Our findings indicate that decreased GFC in the brain regions associated with cerebello-thalamo-cortical network may provide a new insight for understanding the pathological changes of FC in DED. GFC values in the left thalamus may be utilized as a potential biomarker to identify the patients from the controls.
The surface of the eyes is covered with a layer of tears. A stable tear film provides a comfortable environment to the eyes and serves as a good refractive media of the ocular surface (cornea, conjunctiva, accessory lacrimal gland and meibomian gland). Lack of tear membrane components will break the stability of the tear film resulting in ocular surface drying and the conjunctive epithelium damaged by dehydration, known as dry eye syndrome (DED) (Pflugfelder et al., 1999, 2018; Brewitt and Sistani, 2001). DED, one of the most common ocular surface diseases, caused by a variety of factors, is characterized by the absence of mucous components in tears, abnormal distribution of tears on the ocular surface, and increased evaporation of tears (Tabbara and Sharara, 1998). The main clinical manifestations of DED include painful ocular symptoms often described as “burning and painful” (Lemp, 2007; Kalangara et al., 2017; Galor et al., 2018), visual impairment and unstable tear film with potential ocular surface damage, accompanied by increased tear osmotic pressure and ocular surface inflammation (Lemp, 2007). The aggravation of symptoms will lead to a range of serious eye disease including corneal abrasion, filamentous keratitis, and corneal ulcers eventually leading to nubecula and vision loss. However, the pathophysiology of DED remains unclear at present.
In recent years, “neurosensory abnormalities” have been added to the definition of DED (Craig et al., 2017). A growing literature suggested that dry eye symptoms in patients with DED should be better conceptualized as neuropathic eye pain, a manifestation of a central pain processing disorder (Galor et al., 2015; Kalangara et al., 2017; Levitt et al., 2017). Changes in surface microenvironment have taken place to control the secretion of ocular surface glands and blink activity through the ocular sensory pathway. Persistent changes will lead to ocular surface dyskinesia and neuropathic pain caused by functional dysfunction (Belmonte et al., 2017). Dysfunction is a medical condition that normal function of the body is impaired by autonomic dysfunction caused by cortical dysfunction (Liang et al., 2019). Previous studies showed that the occurrence of ocular surface pain is associated with tears deficiency, and repeated ocular sensory nerve injury is involved in the persistence of pain leading to transition from acute pain to chronic pain (Kalangara et al., 2017; Levitt et al., 2017). Sustained structural and functional changes in ocular sensory pathway could cause neuropathic pain and hypotonia of the ocular surface in patients with DED (Belmonte et al., 2017). Functional magnetic resonance imaging (fMRI) has been gradually applied in the clinical researches of DED. Research suggests that chronic peripheral nerve damage will lead to pathological neuroplasticity in the central nervous system (CNS) resulting in reduced nociceptor excitatory thresholds (Levitt et al., 2017). Previous researches showed that patients with DED had extensive alterations in brain function, and abnormal brain function played an important role in the maintenance and development of dry eye symptoms (Rahman et al., 2015; Levitt et al., 2017). fMRI has been used to explore abnormal functional activity in eye diseases such as glaucoma, amblyopia and corneal ulcers (Lin et al., 2012; Li et al., 2014; Rahman et al., 2015; Chen et al., 2017; Wang et al., 2017; Xu et al., 2019). However, it remains unclear whether patients with DED have abnormal GFC in certain brain regions.
To further clarify the pathophysiology of DED, a voxel-wise global-brain functional connectivity (GFC) approach was used to investigate the differences of functional organization between patients with DED and healthy controls. GFC is helpful to obtain voxel-wise whole-brain functional connectivity (FC) in an unbiased way (Li et al., 2016; Cui et al., 2018). Previous researches have demonstrated that GFC is a powerful and replicable data-driven analysis capable of identifying major intrinsic networks (Cole et al., 2009; Murrough et al., 2016). GFC provides an approach to measure the connectivity of all voxels in the brain relative to all other voxels by using a metric that does not require prior selection of seeds or networks (Meier et al., 2016; Zhurakovskaya et al., 2016; Cui et al., 2018; Pan et al., 2019a). Several researches of brain mechanisms have focused on the functional connectivity (FC) between preselected brain regions by using a region of interest (ROI) method (Pires et al., 2012; Meijer and Goraj, 2014; Giorgio and De Stefano, 2016; Lee et al., 2018). This approach is insufficient although the results from ROI are informative. Different ROI selection may yield different results due to potentially biased results based on preselected ROIs. In addition, it may not cover important brain regions associated with the core pathological changes in DED. By contrast, the novel aspect of the present study is that FC abnormalities in patients with DED were examined in a voxel-wise brain-wide way. GFC was considered as a suitable method to examine the differences in a large-scale functional organization in the brain for providing an unbiased way to measure brain function. Therefore, the purpose of the GFC method adopted in the study was to observe the mechanism of brain from the perspective of FC alterations across the whole brain. The potential FC alterations in patients with DED may be detected by this method. Finally, we used the receiver operating characteristic curve (ROC) and support vector machine (SVM) methods to investigate whether abnormal GFC in relevant brain areas could be considered as potential image biomarkers to discriminate patients from controls with good sensitivity and specificity.
Materials and Methods
A total of 20 right-handed patients with DED were recruited from the Frist Affiliated Hospital of Guangxi Medical University. DED was diagnosed according to DED diagnostic guidelines published by the Dry Eye Workshop in 2007 (Lemp, 2007). A total of 23 right-handed healthy controls without symptoms of neurological and ophthalmic disease were recruited from the local community at the same time. All participants aged from 18 to 65 years old. Most symptoms did not differ across different ages. However, burning sensation was common in patients aged 18~45 years old, whereas photophobia was common in patients aged 46~65 years old (Gao et al., 2011). Healthy controls were group-matched with the patients in terms of age, sex ratio, and years of education.
The exclusion criteria for patients with DED were as follows: (1) any history of connective tissue disease, such as rheumatoid arthritis and systemic lupus erythematosus; (2) any history hypertensive encephalopathy, metabolic encephalopathy, and lesions of the CNS caused by infection or other reasons. Healthy controls shared the following exclusion criteria: (1) any history of severe neurological diseases or ophthalmic diseases; (2) any history of serious surgery of internal medicine diseases; and (3) any family history of serious neuropsychiatric disorders or ophthalmic diseases in their first-degree relatives. The participants who did not meet the standard for MRI or showed alterations under conventional MRI scans were also excluded.
The study was approved by the Local Ethics Committee of the First Affiliated Hospital of Guangxi Medical University. All procedures performed in the studies involving human participants were in accordance with the ethical standards of the institutional or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. All participants were provided with a written informed consent before the experiment.
MRI images were captured by a Siemens 3.0T scanner. All participants were required to remain motionless and awake with their eyes closed. Participants were given soft earplugs and foam pads to reduce scanning noise and head motion. Resting-state functional images of slice-order type were obtained by an ascend type using the following parameters: repetition time/echo time = 2,000 ms/30 ms, inversion time = 900 ms, 30 slices, 64 × 64 matrix, 90° flip angle, 240 mm field of view, 4 mm slice thickness, 0.4 mm gap, and 250 volumes lasting for 500 s. After the scan, participants were asked some questions to confirm the wakefulness during the scan. For example, whether the participant fell asleep or the sound frequency of the machine changed during the scanning.
fMRI Data Analysis
Functional images were preprocessed automatically by using the DPABI software (Yan et al., 2016). Because there was time for participants to adapt to the scan environment, the first 10 volumes were removed to ensure a steady-state condition. We excluded the participants whose head movement was more than 2 mm of translation or 2° of rotation in any directions. The time series of fMRI were first corrected for within-scan acquisition time difference between slices and head motion. All realigned images were spatially normalized to the Montreal Neurological Institute (MNI) EPI space in SPM8 and resampled to 3 × 3 × 3 mm3 (Liu et al., 2015). After normalization, the images were smoothed with 4 mm full width on the half maximum Gaussian kernel. Due to the limitation of imaging methods, the lesions of fMRI images were often hid in noise signals which might cause missed diagnosis and even misdiagnosis. The features of the lesion site would be more obvious after proper noise reduction pretreatment. The most common method of noise reduction is smooth. The smooth method is also used in previous studies on ophthalmic diseases (Chen et al., 2017; Wang et al., 2017; Xu et al., 2019). The time series were further linearly detrended and temporally band-pass filtered (0.01–0.08 Hz). Finally, several covariates were removed including Friston-24 head motion parameters obtained by rigid body correction, signal from white-matter centered region, and signal from cerebrospinal fluid. As described in a previous study, the global signal was not removed (Hahamy et al., 2014). The frame displacement (FD) value of each subject was calculated based on a previous study (Power et al., 2012). Scrubbing that removing time point with FD > 0.2 mm was used to control the effect of head motion.
The GFC method used in the study is similar to that used in our previous studies (Cui et al., 2018, 2020; Ding et al., 2019; Pan et al., 2019a,b; Li et al., 2020). GFC is a measure of functional connectivity of all voxels in the brain with respect to other voxels. For each participant, we calculated the average values of correlations between each voxel's time series and every other voxel in the gray matter of the whole brain in MATLAB, which was defined as GFC of this voxel. The threshold setting classified voxel with probability of >0.2 as gray matter, and the gray matter mask would be produced by the gray matter probability map in SPM8 (Liu et al., 2015). The GFC evaluation was defined as,
The process of GFC involves calculating the Pearson correlation coefficients (r) between a given voxel's time series and all other voxels' time series, converting all correlations into Fisher z-scores matrix (Wang et al., 2015; Cui et al., 2018; Pan et al., 2019a,b; Li et al., 2020). The process generates a map for each participant, where each voxel value represents the average connectivity of the voxel with the rest of the brain. The GFC maps were generated by combining GFC of all voxels.
Two-sample t-tests were conducted on the GFC maps between patients with DED and healthy controls. The mean FD and age were used as covariates of no interest. The significance level was set as p < 0.05 by using the family wise error (FWE) correction method.
Mean z values of the brain clusters with abnormal GFC were extracted. Pearson correlation analyses were performed to explore the relationship between abnormal GFC values and illness duration in the DED group after the normality of the data being checked. In the present study, we tested the independent hypotheses on the same data at the significance level of p < 0.05 through an approach of Bonferroni which was the strictest multiple test correction method.
Classification Analysis by Using SVM
As a good tool of classification especially suitable for the case of small samples (Vapnik, 1995), SVM was utilized to examine whether abnormal GFC values in brain areas could be applied to distinguish patients from controls (Chih-Chung and Chih-Jen, 2011) using the LIBSVM software package in the MATLAB. The type of kernel was the default Gaussian kernel in the LIBSVM. The classification performance of the unobserved data was evaluated by dividing the sample set into a training set and a test set. A random SVM cluster was constructed for classification and feature selection based on the brain fMRI data of the subjects. The LIBSVM adopted a “leave-one-out” method that was cross-validated to optimize parameters and obtain satisfactory sensitivity and specificity (Liu et al., 2015; Yan et al., 2016).
The receiver operating characteristic curve (ROC) analysis was used to validate the SVM results, once significant differences in GFC values were observed between the patients and the controls.
Characteristics of the Subjects
The sample of the present study included 20 patients with DED and 23 healthy controls. Continuous variables including age, years of education and FD were analyzed with two-sample t-tests. A Chi-squared test was utilized for sex distribution. Compared with the control group, age (p = 0.23), sex ration (p = 0.19), FD (p = 0.23), and years of education (p = 0.08) of the patient group were not statistically significant. The detailed information of demographic and clinical characteristics of the participants were showed in Table 1.
Group Differences in GFC
Compared with the controls, patients with DED exhibited decreased GFC in the right cerebellum lobule VIII/inferior semi-lunar lobule and left thalamus (Figure 1 and Table 2). No brain region exhibited increased GFC in the patients relative to the controls.
Figure 1. Decreased GFC in patients with DED relative to healthy controls. GFC, global-brain functional connectivity; DED, dry eye disease.
Correlations Between GFC and Clinical Variable
As shown in Figure 2, a positive correlation was observed between GFC values in the left thalamus and illness duration in the patients (r = 0.589, p = 0.006).
Figure 2. A positive correlation (r = 0.589, p = 0.006) between GFC values in the left thalamus and illness duration in patients with dry eye disease. GFC, global-brain functional connectivity.
Distinguishing Patients With DED From Controls
SVM analysis was utilized to determine whether abnormal GFC of these brain regions could distinguish patients with DED from healthy controls. Decreased GFC in the left thalamus exhibited the highest accuracy (88.37%), sensitivity (85.00%) and specificity (91.30%) that could be applied to identify the patients from the controls (Figure 3). As shown in Figure 3 and Table 3, the accuracy of another brain region was unsatisfactory.
Figure 3. 3D view of classified accuracy with best parameters using the GFC values in the left thalamus to differentiate the patients from the controls. The result was obtained in the LIBSVM using a “leave-one-out” approach with default Gaussian kernel. GFC, global-brain functional connectivity.
Table 3. Differentiate the patients from the controls by GFC values in each brain region with the SVM method.
The SVM results were further validated by the ROC method. The results exhibited that the GFC values in the left thalamus could be applied to identify patients from controls with optimal specificity (95.00%) and sensitivity (78.26%) (Figure 4 and Table 4).
Figure 4. Differentiating the patients from the controls with the receiver operating characteristic curve method by using the GFC values in the left thalamus and right cerebellum lobule VIII/inferior semi-lunar.
Compared with healthy controls, patients with DED exhibited significantly decreased GFC values in the cerebello-thalamo-cortical network including the thalamus and cerebellar lobule. Furthermore, GFC values in the left thalamus were positivity correlated with illness duration in the patients. GFC values in the left thalamus could correctly distinguish patients from healthy controls with optimal accuracy, sensitivity, and specificity.
Previous animal neurophysiological studies (Optican and Robinson, 1981; Barash et al., 1999) and human neurotrauma psychology studies (Waespe and Baumgartner, 1992; Panouillères et al., 2013) provided clear evidence that the cerebellum acted as an important role in the saccade adaptation. Alteration in the cerebellum lobule activity was associated with saccade adaptation (Guillaume et al., 2018). Spontaneous nystagmus might affect visual function in patients with unilateral cerebellar lobules and inferior semi-lunar lobule injury (Lee and Kim, 2020). The cerebellar lobule plays an important role in different coding of visual spatial cognition and visual working memory. The cerebellum lobule VIII was the most significantly activated area during above process (Sobczak-Edmans et al., 2016). A previous study observed functional separation between spatial coding and visual working memory processing in the cerebellum lobule (Brissenden et al., 2018). To sum up, cerebellum lobule may participate in the formation of visual working memory. The occurrence of visual impairment in patient with DED may be related to decreased GFC in the right cerebellum lobule VIII/inferior semi-lunar lobule.
The thalamus, located bilaterally around the third ventricle, is the largest elliptic gray matter mixed nucleus mass in the diencephalon. The main function of the thalamus as a relay station for the cerebral cortex to transmit damage information is to synthesize and distribute the sensory to different brain areas (Ab Aziz and Ahmad, 2006; Yen and Lu, 2013). Previous studies showed that the thalamus plays an important role in acute or chronic aching, anesthesia, and analgesia. The thalamus was often activated in fMRI experiments with pain (Heinke and Schwarzbauer, 2002; Martuzzi et al., 2010; Mhuircheartaigh et al., 2010). The descending inhibitory system composed of the thalamus and dorsal root ganglion transmitted aching information to the sensory cortex, and responded to aching during the transmission process and then gave feedback to external traumatic stimuli (Ab Aziz and Ahmad, 2006). Significantly decreased GFC in the left thalamus was observed in patients with DED in this study. Decreased GFC in the thalamus might lead to repeated nerve injury, which changed the stress activated state into the continuously activate state, and then made the pain threshold from enhancement to reduced feedback during the transformation of ocular pain from acute phase to persistent pain in patient with DED. Decreased GFC in the thalamus may be one cause of persistent eye aching in patients with DED. Notably, decreased GFC in the left thalamus was positivity correlated with illness duration in the patients. Abnormal GFC in the left thalamus may be affected the duration of the disease in patients with DED and provides theoretical basis for the pathophysiological of course of the disease in DED.
Specific interconnections exist widely between the thalamus and the cerebral cortex. The thalamic nuclei were connected with the corresponding cerebral cortex to form the thalamus-sensory projection system through these connections, and received and transmitted sensory information in this way (Klingner et al., 2013; Cheng et al., 2015; Penner et al., 2016). The cerebello-thalamo-cortical network was a complete functional network that could carry out visual information and pain information transmission, processing and response in both internal and external of the visual brain region. GFC alteration in the cerebello-thalamo-cortical network may be associated with visual impairment and persistent pain symptoms in patients with DED.
Abnormal GFC values might be utilized as potential biomarkers to identify patients from healthy controls. SVM analyses were conducted to determine whether the GFC values in the right cerebellum lobule VIII/inferior semi-lunar lobule and left thalamus could differentiate the patients from the controls with good sensitivities, specificities, and accuracies. Sensitivity or specificity higher than 0.75 indicated that decreased GFC in the left thalamus to go an accurate indicators (Gong et al., 2011). However, specificity < 0.6 for decreased GFC in right cerebellum lobule VIII/inferior semi-lunar lobule seems to be an indicator with poor accuracy. The accuracy, sensitivity and specificity of GFC values in the left thalamus were more than 0.8 (Table 3), which were optimal for the established diagnostic indicators (Swets, 1988). The ROC results were used to validate the SVM results, which showed that the area under the curve of the left thalamus was more than 0.7, an acceptable accuracy for established diagnostic indicators. Hence, we suggested that decreased GFC in the left thalamus could serve as a potential biomarker to discriminate the patients with DED from the controls.
Several limitations should be considered in the present study. First, the data size is small, the methods and results need to be validated on a larger data set. Second, clinical characteristics such as the data on age of onset were collected retrospectively which might have limited the accuracy of the information. Further assessment of other non-sensory manifestation could be used to provide a comprehensive explanation in patients with DED. Finally, the MNI template used in the study came from a Caucasian population which might bias the present findings from Chinese subjects.
Despite the limitation, the present study indicates that decreased GFC exists in brain regions of the cerebello-thalamo-cortical network in patients with DED. Decreased GFC values in the left thalamus may be utilized as a potential biomarker to differentiate the patients from the controls. Thus, this study provides new insights into the pathological changes of GFC in DED.
Data Availability Statement
All datasets generated for this study are included in the article/supplementary material.
The studies involving human participants were reviewed and approved by Local Ethics Committee of the First Affiliated Hospital of Guangxi Medical University. The patients/participants provided their written informed consent to participate in this study.
WG and SL provided the conception of the work. SW, PP, YO, and WJ collected the data. FL, HL, WL, and YL were responsible for data analysis and interpretation. The manuscript was drafted by PP and critically revised by WG. All authors have given approval to final version of the manuscript. The manuscript was written through contributions of all authors.
This study was supported by grants from the National Key R&D Program of China (2016YFC1307100), the National Natural Science Foundation of China (Grant No. 81771447), and the Guangxi Appropriate Technology for Medical and Health Research and Development Project (Grant No. S201415-05).
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.
The authors thank all individuals with served as the research participants.
Barash, S., Melikyan, A., Sivakov, A., Zhang, M., Glickstein, M., and Thier, P. (1999). Saccadic dysmetria and adaptation after lesions of the cerebellar cortex. J. Neurosci. 19, 10931–10939. doi: 10.1523/JNEUROSCI.19-24-10931.1999
Brissenden, J. A., Tobyne, S. M., Osher, D. E., Levin, E. J., Halko, M. A., and Somers, D. C. (2018). Topographic cortico-cerebellar networks revealed by visual attention and working memory. Curr. Biol. 28, 3364–3372.e65. doi: 10.1016/j.cub.2018.08.059
Chen, W., Zhang, L., Xu, Y., Zhu, K., and Luo, M. (2017). Primary angle-closure glaucomas disturb regional spontaneous brain activity in the visual pathway: an fMRI study. Neuropsychiatr. Dis. Treat 13, 1409–1417. doi: 10.2147/NDT.S134258
Cheng, W., Palaniyappan, L., Li, M., Kendrick, K., Zhang, J., Luo, Q., et al. (2015). Voxel-Based, brain-wide association study of aberrant functional connectivity in schizophrenia implicates thalamocortical circuitry. NPJ Schizophr. 1:15016. doi: 10.1038/npjschz.2015.16
Craig, J. P., Nichols, K. K., Akpek, E. K., Caffery, B., Dua, H. S., Joo, C. K., et al. (2017). TFOS DEWS II definition and classification report. Ocul. Surf. 15, 276–283. doi: 10.1016/j.jtos.2017.05.008
Cui, G., Ou, Y., Chen, Y., Lv, D., Jia, C., Zhong, Z., et al. (2020). Altered global brain functional connectivity in drug-naive patients with obsessive-compulsive disorder. Front. Psychiatr. 11:98. doi: 10.3389/fpsyt.2020.00098
Cui, X., Liu, F., Chen, J., Xie, G., Wu, R., Zhang, Z., et al. (2018). Voxel-wise brain-wide functional connectivity abnormalities in first-episode, drug-naive patients with major depressive disorder. Am. J. Med. Genet. B Neuropsychiatr. Genet. 177, 447–453. doi: 10.1002/ajmg.b.32633
Ding, Y., Ou, Y., Su, Q., Pan, P., Shan, X., Chen, J., et al. (2019). Enhanced global-brain functional connectivity in the left superior frontal gyrus as a possible endophenotype for schizophrenia. Front. Neurosci. 13:145. doi: 10.3389/fnins.2019.00145
Galor, A., Felix, E., Feuer, W., Shalabi, N., Martin, E., Margolis, T., et al. (2015). Dry eye symptoms align more closely to non-ocular conditions than to tear film parameters. Br. J. Ophthalmol. 99, 1126–1129. doi: 10.1136/bjophthalmol-2014-306481
Gong, Q., Wu, Q., Scarpazza, C., Lui, S., Jia, Z., Marquand, A., et al. (2011). Prognostic prediction of therapeutic response in depression using high-field MR imaging. Neuroimage 55, 1497–1503. doi: 10.1016/j.neuroimage.2010.11.079
Hahamy, A., Calhoun, V., Pearlson, G., Harel, M., Stern, N., Attar, F., et al. (2014). Save the global: global signal connectivity as a tool for studying clinical populations with functional magnetic resonance imaging. Brain Connect. 4, 395–403. doi: 10.1089/brain.2014.0244
Heinke, W., and Schwarzbauer, C. (2002). In vivo imaging of anaesthetic action in humans: approaches with positron emission tomography (PET) and functional magnetic resonance imaging (fMRI). Br. J. Anaesth. 89, 112–122. doi: 10.1093/bja/aef155
Kalangara, J., Galor, A., Levitt, R., Covington, D., McManus, K., Sarantopoulos, C., et al. (2017). Characteristics of ocular pain complaints in patients with idiopathic dry eye symptoms. Eye Contact Lens 43, 192–198. doi: 10.1097/ICL.0000000000000249
Klingner, C., Langbein, K., Dietzek, M., Smesny, S., Witte, O., Sauer, H., et al. (2013). Thalamocortical connectivity during resting state in schizophrenia. Eur. Arch. Psychiatr. Clin. Neurosci. 264, 111–119. doi: 10.1007/s00406-013-0417-0
Lee, H., and Kim, H. A. (2020). Anatomical structure responsible for direction changing bilateral gaze-evoked nystagmus in patients with unilateral cerebellar infarction. Medicine 99:e19866. doi: 10.1097/MD.0000000000019866
Lemp, M. A. (2007). The definition and classification of dry eye disease: report of the definition and classification subcommittee of the international dry eye workshop 2007. Ocul. Surf. 5, 75–92. doi: 10.1016/S1542-0124(12)70081-2
Levitt, A., Galor, A., Chowdhury, A., Felix, E., Sarantopoulos, C., Zhuang, G., et al. (2017). Evidence that dry eye represents a chronic overlapping pain condition. Mol. Pain 13:1744806917729306. doi: 10.1177/1744806917729306
Li, H., Ou, Y., Liu, F., Chen, J., Zhao, J., Guo, W., et al. (2020). Reduced connectivity in anterior cingulate cortex as an early predictor for treatment response in drug-naive, first-episode schizophrenia: a global-brain functional connectivity analysis. Schizophr. Res. 215, 337–343. doi: 10.1016/j.schres.2019.09.003
Li, T., Liu, Z., Li, J., Liu, Z., Tang, Z., Xie, X., et al. (2014). Altered amplitude of low-frequency fluctuation in primary open-angle glaucoma: a resting-state FMRI study. Invest. Ophthalmol. Vis. Sci. 56, 322–329. doi: 10.1167/iovs.14-14974
Li, T., Wang, Q., Zhang, J., Rolls, E., Yang, W., Palaniyappan, L., et al. (2016). Brain-wide analysis of functional connectivity in first-episode and chronic stages of schizophrenia. Schizophr. Bul.l 43, 436–448. doi: 10.1093/schbul/sbw099
Lin, X., Ding, K., Liu, Y., Yan, X., Song, S., and Jiang, T. (2012). Altered spontaneous activity in anisometropic amblyopia subjects: revealed by resting-state FMRI. PLoS ONE 7:e43373. doi: 10.1371/journal.pone.0043373
Liu, F., Guo, W., Fouche, J. P., Wang, Y., Wang, W., Ding, J., et al. (2015). Multivariate classification of social anxiety disorder using whole brain functional connectivity. Brain Struct. Funct. 220, 101–115. doi: 10.1007/s00429-013-0641-4
Martuzzi, R., Ramani, R., Qiu, M., Rajeevan, N., and Constable, R. T. (2010). Functional connectivity and alterations in baseline brain state in humans. Neuroimage 49, 823–834. doi: 10.1016/j.neuroimage.2009.07.028
Meier, T., Bellgowan, P., and Mayer, A. (2016). Longitudinal assessment of local and global functional connectivity following sports-related concussion. Brain Imaging Behav. 11, 129–140. doi: 10.1007/s11682-016-9520-y
Mhuircheartaigh, R. N., Rosenorn-Lanng, D., Wise, R., Jbabdi, S., Rogers, R., and Tracey, I. (2010). Cortical and subcortical connectivity changes during decreasing levels of consciousness in humans: a functional magnetic resonance imaging study using propofol. J. Neurosci. 30, 9095–9102. doi: 10.1523/JNEUROSCI.5516-09.2010
Murrough, J., Abdallah, C., Anticevic, A., Collins, K., Geha, P., Averill, L., et al. (2016). Reduced global functional connectivity of the medial prefrontal cortex in major depressive disorder. Hum. Brain Mapp. 37, 3214–3223. doi: 10.1002/hbm.23235
Pan, P., Ou, Y., Su, Q., Liu, F., Chen, J., Zhao, J., et al. (2019a). Voxel-based global-brain functional connectivity alterations in first-episode drug-naive patients with somatization disorder. J. Affect. Disord. 254, 82–89. doi: 10.1016/j.jad.2019.04.099
Pan, P., Wei, S., Ou, Y., Jiang, W., Li, W., Lei, Y., et al. (2019b). Reduced global-brain functional connectivity and its relationship with symptomatic severity in cervical dystonia. Front. Neurol. 10:1358. doi: 10.3389/fneur.2019.01358
Panouillères, M., Alahyane, N., Urquizar, C., Salemme, R., Nighoghossian, N., Gaymard, B., et al. (2013). Effects of structural and functional cerebellar lesions on sensorimotor adaptation of saccades. Exp. Brain Res. 231, 1–11. doi: 10.1007/s00221-013-3662-6
Penner, J., Ford, K., Taylor, R., Schaefer, B., Théberge, J., Neufeld, R., et al. (2016). Medial prefrontal and anterior insular connectivity in early schizophrenia and major depressive disorder: a resting functional MRI evaluation of large-scale brain network models. Front. Hum. Neurosci. 10:132. doi: 10.3389/fnhum.2016.00132
Pflugfelder, S. C., Bian, F., Gumus, K., Farley, W., Stern, M. E., and De Paiva, C. S. (2018). Severity of Sjögren's syndrome keratoconjunctivitis sicca increases with increased percentage of conjunctival antigen-presenting cells. Int. J. Mol. Sci. 19:2760. doi: 10.3390/ijms19092760
Pflugfelder, S. C., Jones, D., Ji, Z., Afonso, A., and Monroy, D. (1999). Altered cytokine balance in the tear fluid and conjunctiva of patients with Sjögren's syndrome keratoconjunctivitis sicca. Curr. Eye Res. 19, 201–211. doi: 10.1076/ceyr.184.108.40.20609
Pires, C., Bahia, C., Rueda Lopes, F., Malfetano, F., Pereira, V., Kubo, T., et al. (2012). Brain MRI abnormalities in Brazilian patients with neuromyelitis optica. J. Clin. Neurosci. 19, 969–974. doi: 10.1016/j.jocn.2011.10.009
Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L., and Petersen, S. E. (2012). Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59, 2142–2154. doi: 10.1016/j.neuroimage.2011.10.018
Rahman, M., Okamoto, K., Thompson, R., Katagiri, A., and Bereiter, D. A. (2015). Sensitization of trigeminal brainstem pathways in a model for tear deficient dry eye. Pain 156:942. doi: 10.1097/j.pain.0000000000000135
Sobczak-Edmans, M., Ng, T. H. B., Chan, Y. C., Chew, E., Chuang, K. H., and Chen, S. H. A. (2016). Temporal dynamics of visual working memory. Neuroimage 124, 1021–1030. doi: 10.1016/j.neuroimage.2015.09.038
Waespe, W., and Baumgartner, R. (1992). Enduring dysmetria and impaired gain adaptivity of saccadic eye movements in Wallenberg's lateral medullary syndrome. Brain J. Neurol. 115, 1123–1146. doi: 10.1093/brain/115.4.1125
Wang, J., Li, T., Zhou, P., Wang, N., Xian, J., and He, H. (2017). Altered functional connectivity within and between the default model network and the visual network in primary open-angle glaucoma: a resting-state fMRI study. Brain Imag. Behav. 11, 1154–1163. doi: 10.1007/s11682-016-9597-3
Wang, L., Xia, M., Li, K., Zeng, Y., Su, Y., Dai, W., et al. (2015). The effects of antidepressant treatment on resting-state functional brain networks in patients with major depressive disorder. Hum. Brain Mapp. 36, 768–778. doi: 10.1002/hbm.22663
Xu, M., Liu, H., Tan, G., Su, T., Xiang, C., Wu, W., et al. (2019). Altered regional homogeneity in patients with corneal ulcer: a resting-state functional MRI study. Front. Neurosci. 13:743. doi: 10.3389/fnins.2019.00743
Zhurakovskaya, E., Paasonen, J., Shatillo, A., Lipponen, A., Salo, R., Aliev, R., et al. (2016). Global functional connectivity differences between sleep-like states in urethane anesthetized rats measured by fMRI. PLoS ONE 11:e0155343. doi: 10.1371/journal.pone.0155343
Keywords: dry eye disease, global-brain functional connectivity, network, functional magnetic resonance imaging, support vector machine
Citation: Pan P, Wei S, Ou Y, Liu F, Li H, Jiang W, Li W, Lei Y, Guo W and Luo S (2020) Reduced Global-Brain Functional Connectivity of the Cerebello-Thalamo-Cortical Network in Patients With Dry Eye Disease. Front. Hum. Neurosci. 14:572693. doi: 10.3389/fnhum.2020.572693
Received: 15 June 2020; Accepted: 24 August 2020;
Published: 25 September 2020.
Edited by:Dajiang Zhu, University of Texas at Arlington, United States
Reviewed by:Lu Zhang, University of Texas at Arlington, United States
Xi Jiang, University of Electronic Science and Technology of China, China
Shijie Zhao, Northwestern Polytechnical University, China
Copyright © 2020 Pan, Wei, Ou, Liu, Li, Jiang, Li, Lei, Guo and Luo. 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.