A Study of the Brain Functional Network of Post-Stroke Depression in Three Different Lesion Locations

Research on the mechanism of post stroke depression (PSD) is the key way to improve the treatment of PSD. However, the functional brain network of PSD has not been entirely supported by the results of functional magnetic resonance imaging (fMRI) studies. The aims of this study are to investigate the brain response of PSD in three different lesions. The brain responses of the three PSD subgroups were similar. However, each subgroup had its own characteristics of the brain network. In the temporal lobe subgroup, the right thalamus had increased degree centrality (DC) values which were different from the other two subgroups. In the frontal lobe subgroup, the left dorsolateral prefrontal cortex, caudate, and postcentral gyrus had increased DC values which were different from the other two subgroups. The hemodynamic response of PSD indicates that PSD has activities of similar emotional networks, of which the negative network realizes its function through the limbic system and default mode network. The brain network has unique characteristics for different lesion locations. The neurological function of the lesion location, the compensatory mechanism of the brain, and the mechanism of integrity and locality of the brain are the important factors in the individual emotional network.

In recent years, resting-state fMRI (rs-fMRI) has been extensively used to understand the mechanism of brain function 12 . Among the methods of rs-fMRI data analysis, functional connection (FC) and independent component analysis (ICA) are the most commonly used for examining connectivity patterns in distinct brain regions 13,14 . In the FC method, researchers generally need to select an a priori defined region of interest (ROI) before analysis. However, for poorly understood brain networks, it is very difficult to select an a priori defined ROI, and this affects the reliability of the conclusions. Compared with the FC method, the ICA method may be performed without a priori selection 15 . although it is difficult to distinguish the results as physiological noise or real brain response; this also affects the reliability of the results. At present, an increasing number of researchers are focusing on a new method termed degree centrality (DC) 16 . DC can count the number of direct connections for a given voxel in a brain network and reflect its FC within that network, without requiring an a priori defined ROI. Based on voxel level, DC technology takes each element as a node, and then calculates the amount of each node connecting with other nodes. This indirectly reflects the position and importance of the node or brain regions in the whole-brain networks 17,18 . Converging evidence indicates that DC is one of the main topological properties measured in graph theory analysis, and provides as an effective index. It has been widely used to identify changes in resting-state functional networks in mental illness, including depression 19 , schizophrenia 20 and autism. However, no literature is available on PSD research based on the DC method. This method provides a promising tool for the elucidation of the neural basis of PSD because it does not require an a priori defined ROI.
The frontal, temporal and parietal lobes serve as the main donor brain areas of the internal carotid artery system, the three lobes are closely related to advanced brain function (i.e. emotion, cognition and memory) and is the predilection site of stroke 21 . Stroke in these brain regions is also the focus of current research.
To reduce heterogeneity, we chose patients with three lesion locations (temporal, frontal and parietal lobes) with similar stroke type, duration, size of lesion, and other basic factors. Moreover, healthy volunteers with approximate conditions were selected as the blank control group (baseline). Our aim in the present investigation was to compare the brain response of PSD and non-PSD (non-depression after stoke) subjects by grouping fMRI datasets into three different lesion locations. We hypothesized that each different type of PSD has a unique feature of the brain network, but there are also some similarities between them. Through this research, we will contribute to the further understanding of the mechanism of PSD, and provide a bridge for future studies.

Method and Subjects
Participants. We reviewed the charts of 133 patients who were admitted for ischemic stroke between December 2012 and June 2017 to the Zhujiang Hospital of Southern Medical University. Patients were included if they met the following criteria. Firstly, the subjects met the WHO criteria for the diagnosis of cerebral infarction, which is based on both the presence of neurological symptoms and a compatible lesion, as demonstrated by magnetic resonance imaging (MRI); secondly, the subjects were in the recovery period (3 months <disease duration <1 year) with stable symptoms; thirdly, the subjects had a single infarcted brain area (3-5 cm) in the right temporal, frontal, or parietal lobe; fourthly, the subjects had a National Institute of Health Stroke Scale (NIHSS) score of six or lower; fifthly, the subjects were conscious and able to cooperate with the interview, and provide informed consent, and complete the scale evaluations and a clinical interview for the diagnosis of depression; sixthly, the subjects had basic self-care ability in their daily lives (Barthel Index score ≥ 60); seventhly; the subjects were aged between 60 and 70 years; eighthly, the subjects did not have a history of hemorrhagic or ischemic stroke; ninthly, the subjects did not have obvious cognitive dysfunction disorder and language understanding disorder; tenthly, the subjects did not have a history of schizophrenia, major depression, anxiety, dementia, drug abuse, or antidepressant use at stroke onset, or a family history of mental disorders; eleventhly, the subjects were not alcoholics or drug abusers. twelfthly, all subjects were right-handed. After a detailed evaluation of inclusion and exclusion criteria, 72 patients who suffered from ischemic stroke participated in this study (Fig. 1).
We also recruited healthy volunteers of similar age, educational background and lifestyle to the stroke patients, with no history of neuropsychiatric disorders or drug abuse. After screening, 15 healthy volunteers were included in the blank control group (baseline). The screened subjects had one or more conditions that did not meet the inclusion criteria.
The following information was collected for each subject: demographics (i.e., age, gender, education level, and whether they lived alone) and stroke severity, as measured by the NIHSS at the time of admission to the hospital. Simultaneously, we obtained scores on the Mini Mental State Examination (MMSE) and the Barthel Index (BI).
All the experiments and protocols were approved by the Ethics Committee of Zhujiang Hospital which is affiliated with the Southern Medical University, China 22 . According to the dictates of the State Council of China, each subject provided written informed consent after receiving detailed instructions and full explanations on the experimental procedures. All methods were performed in accordance with the relevant guidelines and regulations.
The stroke subjects were divided into three subgroups based on their different lesion locations, 1) temporal lobe group, 2) frontal lobe group and 3) parietal lobe group. An experienced neuropsychologist performed the clinical interview to diagnose depression, according to DSM-IV criteria. The severity of depression was assessed using the 24-item Hamilton Rating Scale for Depression . For inclusion in the depression group (PSD group) in our final analysis, participants had to meet DSM-IV criteria for depressive disorder, and score at least 17 on the HAMD. The other participants were assigned to the non-PSD group.
Brain imaging. The experiment was performed in the Department of Radiology of Zhujiang Hospital, Southern Medical University, China. Anatomical scans of the brain were collected prior to stimulation imaging. Then, all subjects were subjected to an rs-fMRI scan for six minutes.
Structural and functional scans were acquired with a 3.0 T Philips Achieva MRI System (Royal Philips Electronics, Eindhoven, The Netherlands) with an eight-channel head array coil equipped for echo planar imaging. The images were axial and parallel to the anterior commissure-posterior commissure line, which covered the whole brain. Structural images were collected prior to functional imaging using a T1-weighted fast spin echo sequence (repetition time/echo time = 500/14 ms, flip angle = 90°, 0.859 mm × 0.859 mm in-plane resolution, slice thickness = 1 mm). Blood oxygenation level-dependent functional imaging was acquired using a T2*-weighted, single-shot, gradient-recalled echo planar imaging sequence (repetition time/echo time = 2000/40 ms, flip angle = 90°, 3.4 mm × 3.4 mm in-plane resolution, 180 time points for a total of 360 seconds). In addition, fMRI image collection was preceded by five dummy scans to minimize gradient distortion.
Preprocessing of experimental functional MRI data. The fMRI image data were preprocessed and analyzed using the Data Processing Assistant for Resting-State fMRI (DPARSF, http://www.restfmri.net) by routines in MATLAB R2010a. The blood oxygen level-dependent (BOLD) time series preprocessing steps included removal of the first 10 volumes, slice-time correction, motion correction, intensity normalization, spatial smoothing, and linear high-pass temporal filtering. The first 10 volumes of each scan were discarded in order to eliminate any non-equilibrium effects of magnetization and to allow subjects to become familiar with the scanning environment. The motion time courses were used to select subjects' head movements of <2 mm in translation and 2° in rotation (no subjects were excluded). Each individual's functional images were normalized using the symmetric echo-planar imaging templates and resampled at a resolution of 3 mm × 3 mm × 3 mm. The normalized functional images were smoothed spatially using a 6 mm full width at half maximum (FWHM) Gaussian kernel 23 . Finally, voxel-wise linear trend removal and temporal high-pass filtering (0.01 Hz < f < 0.08 Hz) were applied. DC Calculation. Weighted DC measures were calculated using the "REST-DC" toolkit in the REST V1.8 package 18 , as previously described 24 . To obtain each participant's graph, Pearson correlation coefficients were computed between the time series of all pairs of brain voxels. Each voxel represented a node in the graph, and each significant functional connection (i.e., Pearson correlation) between any pair of voxels was an edge. As a result, we obtained an n × n matrix of Pearson correlation coefficients between any pair of voxels to construct the whole-brain FC matrix for each participant. Then, individual correlation matrices were transformed into a Z-score matrix using Fisher's r-to-z transformation to improve normality. The weighted DC strength of a voxel as the sum of the connections (Z-values) between a given brain voxel and all other voxels was then computed. To eliminate possible spurious connectivity, we used the Pearson correlation coefficient at r > 0.25 by thresholding each correlation at P ≤ 0.001 18 . Furthermore, standardized weighted DC maps were acquired by subtracting the mean value, and then dividing by the standard deviation within the whole gray matter mask 23 . Statistical analysis. SPSS 18.0 software (SPSS, Chicago, IL, USA) was used to calculate descriptive statistics (mean ± SD) for psychophysical data. All statistical assessments were two-tailed, and we considered results to be significant at p < 0.05, consistent with the preliminary status of the trial.
The DC value differences between PSD and non-PSD, and differences between non-PSD and baseline were calculated using two-tailed, paired t-tests (P < 0.05) and corrected for multiple comparisons [false discovery rate [FDR] (P < 0.05) in rest (http://restfmri.net/forum/rest). The resulting images were shown by rest.

Results
Patient characteristics. 87 subjects (female = 44) were recruited into the study, of whom 72 were stroke patients and 15 were healthy subjects. The mean age of the study sample was 64.60 ± 3.66 (range 60-70) years. Thirty-three of the stroke patients of (45.8%, 33/72) were diagnosed with PSD. There were 21 subjects in the temporal lobe subgroup, 30 subjects in the frontal lobe subgroup, and 21 subjects in the parietal lobe subgroup. A statistical difference was found in the HAMD score between the PSD and non-PSD groups for all three subgroups (p < 0.05). No significant differences were observed in the basic data (i.e., age, sex, education, duration, and whether they lived alone) and functional assessment scores (i.e., MMSE score, BI score and NIHSS score) amongst the three groups (p > 0.05). (See Tables 1 & 2). Hemodynamic responses. Compared with the baseline (healthy volunteers), the DC values of three non-PSD subgroups varied in several brain areas such as prefrontal cortex, limbic lobe, motor cortex, sensory cortex, temporal cortex, parietal cortex and cingulate cortex (See Table 3; Fig. 2). Compared with the non-PSD group, several brain regions of PSD showed decreased DC values, e.g., the ventromedial prefrontal cortex (VMPFC), dorsolateral prefrontal cortex (DLPFC), hippocampal gyrus (HP), thalamus, posterior cingulate cortex (PCC), angular gyrus, and cerebellum lobe. However, DC values were greater in brain areas such as the temporal gyrus, parahippocampal gyrus (PHP), insula, anterior cingulate cortex (ACC), caudate, and supplementary motor area (SMA) (See Table 4; Fig. 2). The brain network of each subgroup showed distinct characteristics. For example, our research demonstrated that the right thalamus in the temporal lobe subgroup showed higher DC values than the other two subgroups (See Table 4 & Fig. 3). In the frontal lobe subgroup, the DC values of the left DLPFC, caudate, and postcentral gyrus were higher than the other two subgroups (See Table 4 & Fig. 4). In addition, for the parietal lobe subgroup, the frontal lobe area of the brain showed decreased DC values compared to the other subgroups (See Table 4 & Fig. 5).

Discussion
Post-stroke depression (PSD) is a common complication of stroke that has a negative impact on the rehabilitation process, and always brings a great economic burden for patients 25 . clarification of the mechanism of PSD is very important, and may be achieved through fMRI, which provides an effective tool for the study of the brain network of the condition.
In the present study, no statistical differences were observed in the basic data (i.e., age and gender) amongst the three groups. This result suggests that, under the strict inclusion criteria, differences in age and gender amongst the three subgroups had minor influence. This led to reduced heterogeneity of the experiment and improvements in the reliability of the conclusion. Previous studies have shown that PSD results from brain injury and nerve damage caused by stroke, the severity of which is an important factor in PSD 26 . In our study, there were no significant differences in NIHSS score, MMSE score and BI score between the PSD and non-PSD groups. Although the stroke subjects showed similar degrees of nerve damage, some individuals did not appear to present depressive symptoms. This suggests that nerve damage may not be a key determining factor for PSD, and other factors are more influential in the onset of depression.

Non-PSD brain response.
Our results indicate that non-PSD can extensively alter the excitability of the brain network, including stroke lesions and remote sites. These altered brain networks are closely associated with emotion, cognition, memory, sensation, movement, etc. We believe that strokes can affect these networks and cause emotional, sensory, motor and other dysfunctions in patients. In addition, the brain responses of different subgroups were not identical and showed different characteristics. For instance, wide inactivation in the limbic system (i.e., insular, ACC, thalamus, and PHP) in the frontal lobe subgroup was apparent, whilst the parietal lobe subgroup showed limbic system over-activation. The main reason for this difference may be that the parietal lobe is predominantly involved in sensory, taste and logical functions 27 , and is not closely related to the limbic system. In contrast, the prefrontal cortex is closely related to the limbic system, with numerous axons connecting each other. Prefrontal cortex and limbic systems are involved in higher brain functions such as emotion and cognition 28 . Frontal stroke causes disruption of nerve connections and decreased transmission of information, resulting in decreased activity of the limbic system. Different brain network reactions can cause different prognosis of stroke patients, the nature of which deserves further attention.

Similar emotional network of PSD. Using structural and functional MRI data from the subjects, we
showed that brain areas located in the prefrontal lobe, thalamus, PCC, angular gyrus, and cerebellum had decreased DC values in the three subgroups. Additionally, the temporal lobe, PHP, insular, ACC, and SMA had increased DC values in these subgroups. These results suggest that strokes of different lesion locations have similar brain networks. Negative brain network. The DLPFC, thalamus, HP, and PCC belong to the limbic system. The limbic system supports a variety of functions including emotion, behavior, motivation and memory 29 . Emotional life is largely housed in the limbic system which is the site of generation and conduction of emotional information 30 . The decreased DC values of the limbic system may affect emotional information processing, which plays a role in PSD. As the functional connection of the limbic system reduces, the whole network is left in a restrained state, which affects other functions of the brain network.
The prefrontal lobe is well known to have a wide range of neural connections and complex structural schemas, as well as rich and complex axonal linkages, which are involved in the processing of emotional information 31 .
Decreased DC values of the prefrontal lobe would lead to a reduction in the neural association between the prefrontal lobe and other brain regions, thereby reducing the transmission of emotional information. In addition, the left prefrontal cortex is closely related to the reward system, which can promote the secretion of dopamine 32 . The decreased DC values of the left prefrontal cortex may reduce the frequency of the reward mechanism function, and aggravate the negative state of stroke patients.
The thalamus is an important part of the limbic system, which is associated with changes in emotional reactivity 33 . In emotional conduction, the thalamus shares a nerve connection with multiple brain areas. For example, the medial dorsal nucleus makes connections with cortical zones of the prefrontal lobe 34 . Also, the anterior nuclei connect with the mammillary bodies, and through them (via fornix), with the HP and the CC 35 . Based on these phenomena, some researchers believe that the thalamus is the relay station of emotional conduction. In our study, the thalamus showed a decreased DC value, which would result in reduced nerve conduction among the prefrontal lobe, HP, CC and thalamus. Moreover, some studies have shown that the frontal lobe and thalamus Continued together constitute the awareness system, which is the main center of spiritual activity (32). Decreased functional connectivity in both brain regions may affect the awareness system, and is a possible reason for the slow-thinking symptom of PSD. Most researchers believe that the HP is associated with recent memory, and involved in emotional reactions and control 36 . Decreased DC values of the HP would affect these emotional functions, and may be the reason that PSD patients are restless and easily angered. Likewise, the decreased memory function of HP also causes memory impairment in PSD patients. The PCC has been strongly linked to emotional salience. Our study showed that the PCC of both sides had decreased DC values; PCC abnormalities may be associated with the negative memory of depressed patients, which is thought to be related to emotional burden.
The prefrontal lobe, PCC, HP and angular gyrus belong to the default mode network (DMN). The DMN has two apparently opposite functions: a spontaneous cognition function and the function of monitoring the environment (sentinel hypothesis) 37 . Our results suggest that PSD may significantly inhibit the functional activities of the DMN, and thereby weaken both of these functions. In addition, PSD may reduce self-evaluation and alerts to environmental awareness, eventually leading to lower episodic memory and alertness, all of which play an important role in the PSD network.       thought to be responsible for emotional feelings, and processes a person's sense of disgust in society. The increased DC values of the insular may aggravate the disgust symptoms of PSD patients. Moreover, as the insular and ACC are important nodes of the pain network 40 , their excessive activation may also explain why many PSD patients suffer pain in the recovery period. The ACC is involved in higher-level functions such as decision-making, impulse control, and emotion 41 . Excessive activation of the ACC has been linked to the transmission of negative emotions. Other studies suggest that the temporal lobe controls the emotion response; if the temporal lobe appears negative during activation, the mood of the patient will stabilize 42 . In our study, the temporal lobe showed increased DC values, which may explain the irritability observed in PSD patients.
Briefly, the similar emotional network consists of a positive and a negative network. The positive network is associated with the generation of the emotion of disgust, negative emotional transmission, and stimulation of a pain sensation, and includes the temporal lobe, PHP, insular, ACC, etc. The negative network is mainly linked to decreased functional connection, emotional network and alertness inhibition, the generation of negative emotions, and awareness system inhibition, and includes brain areas in the limbic system and the DMN. The similar emotional network may be the common neurophysiological basis of PSD.

Individual emotional network of PSD.
Our experimental results indicate that each subgroup has its own characteristics. In the temporal lobe subgroup, the right thalamus had increased DC values. The temporal lobe is involved in emotional processes, and is responsible for recognizing familiar facial emotions and interpreting emotions through someone's body posture. Evidence suggests that the temporal lobe may also be involved in precipitating emotional empathy 43 . Damage to this area can result in problems with memory, understanding language, and maintaining emotional control. Extensive damage to the temporal lobe caused by stroke leads to emotional control reduction of the ipsilateral cerebral hemisphere 42 , and an increase in unstable negative nerve impulses. The thalamus is the relay station of emotional conduction, and is directly connected with the temporal lobe. Stroke injury causes the ipsilateral thalamus to transmit more negative emotions, leading to enhanced thalamus activation.
In the frontal lobe subgroup, the left DLPFC and caudate showed increased DC values which were different from the other subgroups. Based on the mechanism of plasticity 44 , the brain will form compensatory areas in the contralateral hemisphere. As a consequence of damage to ipsilateral brain function, DLPFC (as an important node of the limbic system and DMN) needs to transmit a high level of negative emotional information. Thus, more neurons are required for transmission in contralateral compensatory brain areas, which eventually leads to enhanced compensatory brain area activation. In addition, the compensatory activation of the DLPFC may cause activation of dopamine-related brain areas such as caudate, probably due to the close association between LPFC and the dopamine secretion mechanism 45 .
In the parietal lobe subgroup, the frontal lobe showed substantially larger areas of decreased DC values than the other two subgroups. The parietal lobe integrates sensory information among various modalities, including spatial sense and navigation, and is responsible for brain information integration 46 . Extensive activation of the prefrontal cortex may be due to the confusion within the integrated system of the damaged parietal lobe.
In brief, each subgroup possessed its own characteristics in terms of lesion locations. Four factors contributive to this observation: the neurological function of the stroke site, the compensatory mechanism of the brain, the mechanism of integrity and the locality of the brain.

Study limitations.
Although the subjects were divided into three subgroups and heterogeneity reduced, more lesion locations need to be assessed to improve the reliability of the conclusions. Furthermore, the brain network of PSD of different generations, durations, and lesion sizes require analysis. The sample size in this study was relatively small, and future experiments should recruit more participants to achieve robust conclusions. Moreover, DC measurments require a threshold (r > 0.25), which represents a shortcoming of our experimental method. Thus, we plan to use more advanced methodology in future studies, including participation coefficient and eigenvector centrality. Furthermore, multimodal brain imaging (i.e., diffusion tensor imaging and voxel-based morphometry) is useful for studying gray matter density and anatomical connectivity of the PSD brain, and will be employed in future experiments.

Conclusion
The hemodynamic response of PSD indicates that PSD has activities of similar emotional networks, of which the negative network realizes its function through the limbic system and DMN. The brain network has unique characteristics for different lesion locations. The neurological function of the lesion location, the compensatory mechanism of the brain, and the mechanism of integrity and locality of the brain are the important factors in the individual emotional network.