Skip to main content

Network-level functional connectivity alterations in chemotherapy treated breast cancer patients: a longitudinal resting state functional MRI study

Abstract

Background

Previous studies have found abnormal structural and functional brain alterations in breast cancer survivors undergoing chemotherapy. However, the network-level brain changes following chemotherapy remain unknown. The purpose of this study was to investigate the dynamic changes of large-scale within- and between-network functional connectivity in chemotherapy-treated breast cancer patients.

Methods

Seventeen breast cancer patients were evaluated with resting state functional MRI (rs-fMRI), neuropsychological tests and blood examination before postoperative chemotherapy (t0), one week after completing chemotherapy (t1) and six months after completing chemotherapy (t2). Nineteen age- and education level-matched healthy controls (HC) were also recruited. Independent components analysis (ICA) was performed to assess network component using rs-fMRI data. The functional network changes were then correlated with cognitive assessment scores and blood biochemical indexes.

Results

One-way repeated measures ANOVA revealed significantly changed within-network functional connectivity in the anterior and posterior default mode network (ADMN and PDMN), left and right frontoparietal network (LFPN and RFPN), visual network and self-referential network. Post-hoc test showed that decreased within-network functional connectivity in ADMN, PDMN, LFPN, RFPN, SRN and central network one week after chemotherapy and increased six months after chemotherapy (all P < 0.05). As for the between-network functional connectivity, the PDMN- sensorimotor network connectivity showed the same tendency. Most of these within- and between-network functional connectivity changes were negatively associated with blood biochemical indexes and cognitive assessment scores (all P < 0.05).

Conclusions

These results indicated that chemotherapy may induce widespread abnormalities in resting state networks, which may serve as a potential biomarker of chemotherapy related cognitive impairment, providing insights for further functional recovery treatment.

Background

Chemotherapy-related cognitive impairments (CRCI) are common in non-central nervous system cancers which show multiple-domain cognitive deficits such as attention, executive function, learning ability, memory and information processing speed during or after chemotherapy [1, 2]. Although accumulating research had focused on these cognitive problems, the neural mechanism remains largely unclear. Cross-sectional [3,4,5] and longitudinal neuroimaging studies [6,7,8,9] have shown structural and functional alterations during or after chemotherapy, thereby providing neural substrates for CRCI. Brain activation and connectivity studies also showed that breast cancer survivors had abnormal brain activity and decreased neural network transfer efficiency from large-scale perspectives [10,11,12]. These damaged brain structural and functional areas partially overlapped, mostly located in the frontal and temporal lobes. This indicates that chemotherapy induced common cognitive deficits may be related to the large-scale abnormal brain activity/connectivity.

Human cognitive activity is the result of the synergetic action of multiple functional neural networks in the brain, not only from the activities within the network, but also from the synergetic activities between the networks [13, 14]. Resting-state functional MRI (rs-fMRI) evaluates the resting state brain functional connectivity and makes it possible to describe the interaction of sub-networks, which are spatially distinct regions [15], and facilitates the understanding of the cooperative relationship of different brain regions at a large scale. Studies have found about 7 common resting state networks (RSNs) in healthy adults [16]. Alterations in RSNs have been observed in many neuropsychiatric diseases such as Parkinson’s disease, mild cognitive impairment and dementia [17,18,19]. Default mode network (DMN) is the most commonly investigated RSN in breast cancer survivors for its preferential vulnerability and sensitivity to tumor chemotherapeutics [20]. Decreased functional connectivity of DMN has been reported particularly in the medial prefrontal cortex, posterior cingulate cortex and the medial temporal cortex [15, 21]. Our previous study found abnormal hippocampal connectivity in superior/middle temporal gyrus, insula and frontal gyrus, most are parts of the RSNs [22]. These altered functional connectivities were involved in the regulation of executive, memory and emotion, supporting chemotherapy induced widespread cognitive networks disruption.

Most of the previous CRCI-related rs-fMRI researches focused on a specific RSN alteration in breast cancer patients [23, 24]. To the best of our knowledge, no studies have reported the dynamic changes of RSN functional connectivity in breast cancer survivors from a large-scale network perspective. Therefore, the purpose of this study was to investigate the dynamic changes of within- and between- network functional connectivity in chemotherapy-treated breast cancer survivors. We hypothesized that chemotherapeutic agents not only induced extensive RSNs abnormality, but also disrupted the connections between these RSNs. The correlations between changes in RSNs and potential risk factors such as anxiety, depression, hemoglobin and blood lipids, as well as cognitive assessment scores were also analyzed.

Materials and methods

Participants

The study was approved by the Medical Ethics Committee of Jinling Hospital (2016NZGKJ-069). All participants signed the written informed consent prior to the study. Women newly diagnosed with primary non-metastatic breast cancer (stages I-III) were recruited from March 2017 to December 2018. Brain MRI scans, cognitive assessment and blood biochemical examinations were performed at baseline (after surgery but before the start of chemotherapy or before neoadjuvant chemotherapy, t0), one week after completing chemotherapy (t1, mean 6 days), and six months after completing chemotherapy (t2, mean 171 days). Healthy controls (HC) recruited from the local community were evaluated at baseline. The inclusion criteria were: right handedness; with no history of psychiatric, neurological diseases and smoking; with 9 years or higher educational level; no chronic illness, no intracranial radiotherapy and other diseases that might affect brain function or cognitive assessment. Study flowchart with included and excluded patient numbers is shown in Fig. 1. Seventeen female breast cancer patients (mean age: 45.6 ± 10 years old) completing six to eight courses of chemotherapy and 19 age (mean age: 45.8 ± 9.4 years old) and education level matched HCs were finally included in this study (Fig. 1). There was no concurrent radiotherapy or endocrine therapy during chemotherapy for the 17 breast cancer patients, but 6/17 (35.3%) patients received endocrine therapy after the end of chemotherapy.

Fig. 1
figure1

Flowchart of this study. BC, breast cancer; HC, healthy controls

Neuropsychological assessment

All participants completed the following neuropsychological assessments before each MR examination: 1) Mini-mental State Examination (MMSE) and Montreal Cognitive Assessment (MoCA), which were used to assess the participants’ common cognitive function; 2) Self-Rating Anxiety Scale (SAS) and Self-Rating Depression Scale (SDS), which were used to evaluate anxiety and depression state; 3) Number Connection Test type A (NCT-A) and Digital Symbol Test (DST) for the assessment of attention and executive ability; 4) Auditory verbal learning memory (WDT) for intelligence and memory assessment; 5) Line Tracing Test (LTT) for the visual ability; 6) Serial Dotting Test (SDT) for assessing the fine motor skills; 7) The Chinese version of the Stroop Color- Word Test (SCWT) for measuring reaction and conversion [21, 25]. All participants were scored according to the standard of each scale. The specific operation and instructions were completed by the same skilled researcher.

Blood biochemical examinations

Fasting venous blood samples of all participants were obtained on the day of their MRI examinations to assess basic endocrine, immune and digestive function. Blood tests included Estradiol (E2), prolactin (PRL), follicle stimulating hormone (FSH), luteinizing hormone (LH), total cholesterol, triglyceride, high sensitivity C-reactive protein and hemoglobin (Hb) concentration. The fasting venous blood samples were obtained in the morning after sitting quietly for 15 min. For HCs with normal menstrual cycles, blood samples were collected 3–10 days after menstruation. All blood samples were collected, stored and tested in the same way.

MRI data acquisition

Breast cancer patients underwent brain MRI scans at three time points on the same GE 3.0-T MRI scanner (Discovery MR 750, GE Healthcare, Milwaukee, USA) equipped with a standard 32-channel head coil. Our structural and functional MRI scan paradigms had been described in our previous study in detail [22]. The sagittal three-dimensional T1-weighted (3D-T1) images for anatomic reference were collected by using 3D inversion recovery prepared fast spoiled gradient recalled sequence. Repetition time (TR) = 8.14 ms, echo time (TE) = 3.17 ms, flip angle (FA) = 12°, matrix = 256 × 256, number of slices = 176 and slice thickness = 1.0 mm. Rs-fMRI data were acquired with a gradient echo EPI sequence during the motionless and unintentional thinking condition: total volume = 250, slices number = 10,750, TR = 2000 ms, TE = 30 ms, FA = 90°, matrix size = 64 × 64, fields of view (FOV) = 240 × 240 mm2. Routine T2- fluid attenuated inversion recovery (FLAIR) sequence was used to exclude primary brain lesion. The excluded primary brain lesions were brain infarction (n = 1) and multiple sclerosis (n = 1) at t0, brain metastases (n = 1) at t1 (Fig. 1). The HC group completed the MRI examinations at the baseline, follow-up examinations were not performed.

Image preprocessing

The fMRI data analysis was done by using Data Processing Assistant for Resting-State fMRI (DPARSF) (http://rfmri.org/DPARSF). After removing the first 10 time points for the signal stabilization, a total of 240 time points were left for the following processing: slice timing correction was used to eliminate time errors between slices; realignment and head motion correction. Participants with head motion > 1.5 mm translation or rotation > 1.5° would be removed. Following realignment, there were 3 steps: 1) co-registration of functional images to the participants’ structural images; 2) T1 structural images were bias-corrected, segmented into grey matter, white matter and cerebrospinal fluid, and initially normalized to standard Montreal Neurological Institute (MNI) space; 3) the normalization parameters obtained were then applied to the functional images to normalize them into MNI space, and resampled with a voxel size of 3 × 3 × 3 mm3, the final images were smoothed using a 6 × 6 × 6 mm Gaussian kernel at full-width at half-maximum.

Independent component analysis (ICA)

Analyses were conducted using ICA with GIFT software (http://icatb.sourceforge.net/). ICA is a data-driven, multivariate approach without priori assumptions which can separate unknown mixed fMRI signal sources into maximum spatial activation maps or independent temporal components. Independent components were estimated by GIFT software from pre-processed data of t0, t1 and t2 time points in the breast cancer group and HC group at baseline. Principal component analysis (PCA) was used to reduce the data dimension. Using the minimum description length criteria, 39 independent components were separated from four groups of smooth data dimensions [26]. We selected the ICASSO algorithm as stability analysis type by running ICA estimation 100 times. Subsequently, the Z value of each independent component was obtained by back-reconstruction using spatial-temporal regress. The z-value measured the correlation between the time series of each voxel and each independent component, which was used for further statistical analyses. Finally, the time series of 39 independent components reflecting the spontaneous brain activity and the spatial activation distribution map reflecting the intensity of brain activity were obtained. According to the maximum spatial correlation, the normalized average spatial activation map and the interested spatial template were calculated using automatic template matching process in the GIFT software. The component with the best “goodness of fit” of standard RSNs template was taken as the final resting state sub-network map [27]. In this study, 9 sub-network templates of RSN were used: default mode network (DMN), frontoparietal network (FPN), dorsal attention network, (DAN), sensorimotor network (SMN), central executive network (CEN), self-referential network (SRN), visual network (VN), auditory network (AN), and central network (CN).

Statistical analyses

Demographic and clinical data

SPSS17.0 software (SPSS Inc. Chicago, IL) was used for the evaluation of the demographic and clinical data of all participants. Continuous variables were calculated the medians and the 25–75% interquartile ranges (IQRs). At baseline, the two-sample t-test was used to examine the differences between the breast cancer group and HC group. Non-normal distribution data were analyzed using two independent sample nonparametric tests (Mann-Whitney U method). One-way repeated measures ANOVA was used to test the within participants’ effects in breast cancer group, time as repeated effect, clinical data as the dependent variables. The P value was corrected by the Greenhouse-Gesisser correction method when it failed the spherical test. Bonferroni post hoc test was used to conduct pairwise comparisons between groups (t0, t1, t2).

Within–network functional connectivity

SPM8 software was used to obtain the group-level spatial distribution map of each RSN subnetwork by using the one sample t-test, family-wise error (FWE) corrected (P < 0.05). In order to restrict the within–network functional connectivity differences, a group mask was also generated. (1) Two-sample t-test was used to compare between the patients and controls at baseline, the Gaussian Random Field theory correction (GRF) (voxel P value < 0.001, cluster P value < 0.05) was used for multiple comparisons. Parameters such as SAS, SDS, age, education level and gray matter volume that may affect functional connectivity were used as covariates. (2) One-way repeated measures ANOVA was used to assess the dynamic changes of within-network functional connectivity in breast cancer patients at three time points (GRF correction, voxel P value < 0.001, cluster P value < 0.05). In the ANOVA model, we used REST software (http://www.restfmri.net/) to extract the average time series of significantly different brain regions within the subject, and evaluated the changes at t0, t1, and t2. Pairwise comparisons based on three time points were performed in the breast cancer group (post hoc test, P < 0.05).

Between-network functional connectivity

ICA time series was used to analyze between-network connectivity and the functional connectivity coefficients between the pairs of RSNs were obtained using the GIFT toolkit. The functional connectivity coefficients of all pairs of sub-networks were analyzed and the correlation coefficients results were converted to z values. In the breast cancer group, one-way repeated measures ANOVA was used to test the between-network functional connectivity coefficients differences at three time points (Bonferroni multiple comparison correction, P < 0.05). Post hoc test was used to conduct pairwise comparisons between groups (t0, t1, t2), and P < 0.05 was considered to be statistically significant.

Correlation analysis

According to the formula (Δ value = t2-t1), the z value changes of significantly different within-network connectivity and between-network connectivity of each participant were calculated. Pearson correlation analysis was used to analyze the correlation between brain function connectivity parameter changes and clinical index changes. P < 0.05 was considered to be statistically significant.

Results

Demographic and clinical data

A total of 17 breast cancer patients receiving chemotherapy completed all tests at three time points (t0, t1 and t2). Nineteen HCs with matched age and education level only completed baseline assessments (t0). In breast cancer group, 11 women had normal menstrual status at baseline, 10 of the 11 premenopausal women had amenorrhea (90.9%, 10/11) at t1, and one had irregular menstrual cycles (9.1%, 1/11). In the breast cancer group, 6 patients received endocrine therapy after the end of chemotherapy. Six patients (35.3%, 6/17) recovered their menstruation at t2.

At baseline, there was no significant difference in menstrual status, E2, LH, and FSH levels between the breast cancer patients group and HC group (all P > 0.05). The baseline anxiety and depression scores in the breast cancer group were significantly higher than those in the HC group, however, these scores were in the normal range (P = 0.012, P = 0.002; respectively). The total cholesterol, triglyceride and fasting blood glucose in the breast cancer group were higher than those in the HC group (all P < 0.05), and the scores of SDT, WDT and Stroop test were higher than those of HC group (all P < 0.05). Demographic data and clinical data at baseline are shown in Table 1.

Table 1 Demographic, Clinical, and Neuropsychological Data at Baseline

Demographic and clinical data changes

One-way repeated measures ANOVA showed that E2 values were significantly different among three time points (F = 5.089; P = 0.012). Post-hoc tests showed significant decrease from t0 to t1 (P = 0.041), and a rising trend 6 months after chemotherapy. Pairwise comparisons showed significantly increased triglycerides from t0 to t1 (P = 0.041) and decreased triglycerides levels at t2 without statistical significance. The Hb levels decreased slightly from t0 to t1, and significantly increased from t1 to t2 (P = 0.003) (Table 2).

Table 2 Blood Examination and Cognitive Assessment Results in Breast Cancer Group at Three Time Points

Repeated measures ANOVA showed significantly different WDT and LTT scores (P = 0.008, P = 0.013, respectively), the breast cancer patients showed slightly decreased WDT scores from t0 to t1, and significantly decreased scores from t0 to t2 (P = 0.002); LTT scores decreased significantly at t1 (P = 0.010) and slightly lower from t1 to t2, while the difference was statistically significant from t0 to t2 (P = 0.012) (Table 2). The NST, SDT and Stroop-test scores fluctuated after chemotherapy, but the differences were not statistically significant. Anxiety and depression scores were increased from t0 to t1 and reduced from t1 to t2, but the difference was not statistically significant (Table 2).

Spatial distribution of RSNs

Thirty nine independent components were separated from the temporally concatenated 4D population data. After automatic template matching and visual discrimination confirmation, networks with the best matching were selected as the interested network. They were anterior and posterior DMN (ADMN, PDMN), DAN, left and right FPN (LFPN, RFPN), SMN, CEN, SRN, VN, AN and CN. One sample t-test showed that both breast cancer group and HC group demonstrated a typical spatial distribution pattern of RSNs.

Within-network functional connectivity

Independent sample t-test showed no significant difference in within-network functional connectivity between the breast cancer group and HC group at baseline (GRF corrected).

The component of ADMN, PDMN, LFPN, RFPN, CN, SRN and VN revealed significantly different functional connectivity between time points in the breast cancer group (GRF corrected) (Tables 3, 4). The locations of the peak value cluster are shown in Table 3. Post-hoc tests showed significantly increased within-network functional connectivity from t0 to t1 in ADMN, PDMN, LFPN, RFPN, SRN and CN (P = 0.015, 0.001, 0.037, 0.006, 0.001 and 0.011, respectively) and significantly deceased connectivity in these above-mentioned networks from t1 to t2. Only the VN functional connectivity decreased slightly from t0 to t1 and significantly increased from t1 to t2 (P = 0.002) (Fig. 2).

Table 3 Statistically Significant Differences of Within-network Functional Connectivity for Breast Cancer Group
Table 4 Summary of Within-network Functional Connectivity for Breast Cancer Group at three time-points
Fig. 2
figure2

Seven RSNs show significant within-network connectivity changes in breast cancer patients among three time points. The histograms display the longitudinal evaluation of the within-network functional connectivity changes in the peak cluster regions for each RSN at three time points. t0, baseline assessment; t1, one week after chemotherapy; t2, six months after chemotherapy. SFG. R, right superior frontal gyrus, PCC, posterior cingulate gyrus, SFG. L, left superior frontal gyrus, Tri-IFG.R, Pars triangularis of right inferior frontal gyrus, SMA. R, right supplementary motor area, Orb-IFG.L, Left orbital inferior frontal gyrus, Calc. R, right calcarine sulcus cortex

Between-network functional connectivity

The between-network connectivity did not show significant difference between the patient and HC group at baseline. One-way repeated measure ANOVA analysis showed that the connectivity between ADMN and CN, PDMN and SMN, SMN and VN were different at three time points. Post-hoc test showed the functional connectivity between PDMN and SMN decreased from t0 to t1 (P = 0.051), while increased from t1 to t2 (P = 0.023). The functional connectivity between ADMN and CN, SMN and VN showed continuous decrease (all P < 0.05) (Fig. 3).

Fig. 3
figure3

The between-network connectivity changes in breast cancer patients among three time points. t0, baseline assessment; t1, one week after chemotherapy; t2, six months after chemotherapy. ADMN, anterior default mode network; PDMN, posterior default mode network; SMN, sensorimotor network; VN, visual network; CN, central network

Correlation analysis

The decreased between-network functional connectivity changes of ADMN and CN were negatively correlated with the decreased fasting glucose and increased DST score changes (r = − 0.497, P = 0.043; r = − 0.547, P = 0.035, respectively). The decreased between-network functional connectivity changes of SMN and VN were negatively correlated with increased changes in blood estrogen levels and decreased SDS scores (r = − 0.655, P = 0.039; r = − 0.498, P = 0.041, respectively) (Fig. 4). The decreased within-network connectivity changes in the supplementary motor area for CN were positively correlated with decreased fasting blood glucose changes (r = 0.561, P = 0.019) and negatively correlated with increased total cholesterol changes (r = − 0.484, P = 0.049). The increased within-network connectivity changes in the right calcarine sulcus cortex for VN were negatively correlated with the decreased changes of Stroop-W (r = − 0.563, P = 0.019).

Fig. 4
figure4

The correlation between functional connectivity changes and clinical variables changes in breast cancer survivors from t1 to t2. a-b: Correlation analysis of changes in the functional connectivity between ADMN and CN with changes in fasting glucose and DST scores. c-d: Correlation analysis of changes in the functional connectivity between SMN and VN with changes in blood estrogen levels and SDS scores. t1, one week after chemotherapy; t2, six months after chemotherapy. ADMN, anterior default mode network; SMN, sensorimotor network; VN, visual network; CN, central network; SDS, Self-Rating Depression Scale; DST, digit symbol test

Discussion

In this study, we found increased within-network connectivity in most RSNs in breast cancer patients one week after chemotherapy, and all the within-network connectivity of RSNs tended to recover towards the baseline level six months after chemotherapy. Chemotherapeutics showed selective damage to the between-network connectivity. The connectivity between PDMN and SMN showed a recovery trend six months after chemotherapy, while the connectivity between ADMN and CN, SMN and VN continued to decline. Furthermore, both within- and between-network connectivity changes were significantly correlated with blood biomarkers and cognitive function alterations. This prospective longitudinal study provided the evidence that chemotherapy may induce widespread connectivity abnormalities in RSNs, which might serve as potential biomarkers of chemotherapy related cognitive deficits in breast cancer patients.

Our results suggested that chemotherapeutics can damage the large-scale resting state brain networks, including ADMN, PDMN, LFPN, RFPN, SRN, CN and VN. These findings provided the evidence that chemotherapy induced widespread cognitive deficits across various domains. DMN is one of the most common RSNs, which is considered to support processes such as active episodic memory and introspection, and would be deactivated during specific goal-directed tasks [28]. Previous studies have suggested that the ADMN is associated with emotional regulation, while the PDMN is involved in consciousness formation and memory [29, 30]. In the study, patients mainly exhibited memory impairment which may be associated with abnormal PDMN functional connectivity. CN is related to activity inhibition, playing an important role in adaptive cognitive control [16]. In this study, changes of ADMN-CN functional connectivity were negatively correlated with DST scores. This result indicated that dysregulated ADMN-CN may induce executive functioning and processing speed impairment. Interestingly, mind tracking and visual mobility showed improved in breast cancer survivors. We hypothesized that chemotherapy induced decreased VN and increased SMN functional connectivity may contribute to effectively coordinate this complex visual movement. However, we also found some slight but not significant cognitive impairment in attention, executive and reaction ability. The abnormal cognitive function may reflect brain mechanisms for attention and executive regulation involving DMN, CN and FPN [31].

In our study, chemotherapy induced abnormal RSNs functional connectivity was mainly located in the superior frontal gyrus, PCC, supplementary motor area, orbital inferior frontal gyrus and calcarine sulcus, which was consistent with the most previous studies [12, 20, 32]. We speculated that these regions may be more vulnerable to chemotherapeutic agent attack. PCC and superior frontal gyrus mainly located in the task-negative network and deactivated for self-referential processing tasks, while the supplementary motor area and orbital inferior frontal gyrus located primarily within or near the task-positive network and activated for demanding cognitive tasks. These different task-evoked neuronal responses were considered to support the process of target activity by regulating the allocation of neural resources [33]. We speculated that the increased functional connectivity may represent a compensation for chemotherapy induced cognitive dysfunction which required the recruitment of more neural regions and increased the strength of functional connectivity to maintain normal neural activity. Decreased functional connectivity in these regions may be attributed to the reallocation of neural resources in information processing [34, 35]. The functional connectivity strength of RSNs might predict chemotherapy induced higher or lower neural activity during task [30]. Previous multitask-based fMRI of CRCI showed increased [36,37,38] and decreased activity [39,40,41] located mainly in several frontal and parietal regions, anterior cingulate cortex and supplementary motor area. The functional connectivity regions of altered RSNs in our study were matched well with previous task-related studies.

Chemotherapeutics also showed selective damage to the between-network connectivity in this study. The DMN can integrate primary perception and advanced cognitive functions. It accepts more information from other RSNs than it outputs. Liao et al. showed that SRN, CEN, CN and AN had significant effect on connections to the DMN [42]. Sridharan et al. found that CN may regulate DMN, CEN and dorsal attention network [43]. The collaborative work among these RSNs completes the integration and transformation of information. In this study, the decreased functional connectivity between networks after chemotherapy suggested that chemotherapeutic drugs may alter function central brain regions through multiple inter-network connectives. The chemotherapy related cognitive impairments were the results of reduced coordination between multiple networks. The RSNs competed for more processing resources from the “central cognitive operator” to compensate for their impaired cognition [34]. This may explain the impairment of multiple cognitive domains after chemotherapy, but no significant difference was observed between baseline and post-chemotherapy.

Interestingly, most of the within-network connectives tended to recover to the baseline levels six months after chemotherapy, while the between-network connectives showed partial recovery. Recovery of some brain regions to baseline levels has also been reported in other CRCI-related functional and structural MRI studies [8, 44,45,46], and the acute damage mostly appeared 1 month after chemotherapy and (partially) recovered one year later. These studies suggested that chemotherapy-induced brain structural and functional damage may be temporary, the frontal regions (such as the frontal and temporal gyrus) may recover over time, and brain abnormalities in the posterior region may persist for a long time [47]. Brain function and structural recovery may be attributed to neuroplasticity mechanisms. In this study, most of the subjects were young women, cognitive challenges in daily work and social activities promoted early rehabilitation through neuroplasticity mechanisms. Additionally, the frontal regions began to show signs of recovery six months after chemotherapy. We hypothesized that the frontal lobe might be one of the first sensitive brain region to experience functional recovery. However, our finding was inconsistent with some previous studies as for recovery time, which may be related to the age, chemotherapy dose, chemotherapeutic drug, treatment stage, cognitive function type and the different control groups [1, 8, 48]. The recovery to the baseline level may reflect compensation ability improvements over time to some extent, it does not mean a return to normal level, and some brain dysfunction may persist for a long time. Sustained networks alterations may further affect patient’s cognitive function, such as sustained memory deterioration and executive ability impairment, which showed poor WDT and NST scores in this study. We are not sure whether the brain abnormality in the posterior brain region (such as visual function areas) alteration was a sustained change in acute effects after chemotherapy or a delayed brain injury occurring at a certain time after chemotherapy. Further long-term follow-up studies are thus needed.

Breast cancer patients also showed decreased visual ability and memory at baseline compared to HC group. This pre-treatment mild cognitive decline may be associated with tumor-related physical and psychological stress [41, 49]. Additionally, the improved test scores in LTT, SDT and Stroop test after chemotherapy may be related to the compensatory effect of the functional network. Furthermore, partial networks alterations were associated with estrogen, fasting blood glucose and blood lipids changes, though the improved blood indicator was still in the abnormal range compared to HC group. We speculated that the increased blood estrogen levels, decreased blood sugar and lipid levels six months later may play a positive role in the cognition improvement. Estrogen can alter the metabolic level of the frontotemporal lobe, affecting the structure and function of specific brain regions [50]. Animal experiments and clinical practice have shown that estrogen can be used to treat attention and memory loss. Abnormalities in blood lipids, fasting blood glucose, and hemoglobin have also been found in neuropsychiatric diseases such as Alzheimer’s disease and diabetic encephalopathy [51, 52]. The improved cholesterol metabolism would reduce the risk of vascular related cognitive impairment. However, the use of chemotherapeutic agents and endocrine therapy confused the assessment of underlying metabolic changes. Although there was no direct evidence that baseline metabolism changes were significantly associated with cognitive improvement six months after chemotherapy, we speculated that the improvement of baseline metabolism in patients can play an irreplaceable role in the recovery mechanism of CRCI.

Anxiety and depression scores were significantly higher in the breast cancer group than HC group at baseline, and continued to increase one week after chemotherapy, but decreased over time after the end of chemotherapy. The dynamic changes of depression and anxiety scores were similar as that of some cognitive functions. Studies had suggested that there were competitive interactions between emotion and cognition [7, 35], we speculated that the decreased negative emotion may benefit cognitive function recovery of patients in this study.

The study had some limitations. Although HC was included in this study, the control participants had no serial follow-up, making it impossible to evaluate the interaction between the groups and time. However, all breast cancer patients performed a series of neuropsychological tests, blood examination and MRI scans at three time points, which enabled the longitudinal analysis of RSNs in these patients. Secondly, the breast cancer patients were heterogeneous, some confounding factors such as different severities of the disease and different treatment strategies may induce some impacts on the results of our study. Thirdly, the basic metabolism and hormone level may fluctuate greatly in a relatively short time, and untangling these effects was difficult. Fourthly, the sample size is relatively small due to some patients refusing follow-up. Finally, memory effect can be present because of repeated neuropsychological testing at short time intervals, which may lead to no significant change and even improved in the cognitive test.

Conclusion

This prospective longitudinal study found that the breast cancer survivors showed a large-scale functional connectivity damage after chemotherapy, and these impaired RSNs partly tended to return to the baseline level 6 months after chemotherapy. The altered functional connectivity was related to the patient’s cognitive function and hematology changes. Therefore, RSNs could be promising markers for evaluating potential chemotherapy-related brain damage, which provides an important basis for the observation of brain function impairment process and the subsequent cognitive rehabilitation interventions.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

CRCI:

Chemotherapy related cognitive impairment

fMRI:

Functional magnetic resonance imaging

HC:

Healthy controls

RSNs:

Resting state networks

E2:

Estradiol

DMN:

Default mode network

FPN:

Frontoparietal network

DAN:

Dorsal attention network

SMN:

Sensorimotor network

CEN:

Central executive network

SRN:

Self-referential network

VN:

Visual network

AN:

Auditory network

CN:

Central network

PCC:

Posterior cingulate cortex

References

  1. 1.

    Wefel JS, Kesler SR, Noll KR, Schagen SB. Clinical characteristics, pathophysiology, and management of noncentral nervous system cancer-related cognitive impairment in adults. CA Cancer J Clin. 2015;65(2):123–38.

    PubMed  Article  Google Scholar 

  2. 2.

    Rick O, Reuß-Borst M, Dauelsberg T, Hass HG, König V, Caspari R, Götz-Keil G, Pfitzner J, Kerschgens C, Fliessbach K, Hoppe C. NeuroCog FX study: a multicenter cohort study on cognitive dysfunction in patients with early breast cancer. Psychooncology. 2018;27(8):2016–22.

    PubMed  Article  Google Scholar 

  3. 3.

    Li X, Chen H, Lv Y, Chao HH, Gong L, Li CSR, Cheng H. Diminished gray matter density mediates chemotherapy dosage-related cognitive impairment in breast cancer patients. Sci Rep. 2018;8(1):13801.

  4. 4.

    Apple AC, Ryals AJ, Alpert KI, Wagner LI, Shih PA, Dokucu M, Cella D, Penedo FJ, Voss JL, Wang L. Subtle hippocampal deformities in breast cancer survivors with reduced episodic memory and self-reported cognitive concerns. Neuroimage Clin. 2017;14(C):685–91.

    PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Koppelmans V, de Groot M, de Ruiter MB, Boogerd W, Seynaeve C, Vernooij MW, Niessen WJ, Schagen SB, Breteler MM. Global and focal white matter integrity in breast cancer survivors 20 years after adjuvant chemotherapy. Hum Brain Mapp. 2014;35(3):889–99.

    PubMed  Article  Google Scholar 

  6. 6.

    Deprez S, Vandenbulcke M, Peeters R, Emsell L, Smeets A, Christiaens MR, Amant F, Sunaert S. Longitudinal assessment of chemotherapy-induced alterations in brain activation during multitasking and its relation with cognitive complaints. J Clin Oncol. 2014;32(19):2031–8.

    PubMed  Article  Google Scholar 

  7. 7.

    Jung MS, Zhang M, Askren MK, Berman MG, Peltier S, Hayes DF, Therrien B, Reuter-Lorenz PA, Cimprich B. Cognitive dysfunction and symptom burden in women treated for breast cancer: a prospective behavioral and fMRI analysis. Brain Imaging Behav. 2017;11(1):86–97.

    PubMed  Article  PubMed Central  Google Scholar 

  8. 8.

    Billiet T, Emsell L, Vandenbulcke M, Peeters R, Christiaens D, Leemans A, Hecke WV, Smeets A, Amant F, Sunaert S. Recovery from chemotherapy-induced white matter changes in young breast cancer survivors. Brain Imaging Behav. 2017;12(1):1–14.

    Google Scholar 

  9. 9.

    Chen H, Ding K, Zhao J, Chao HH, Li CR, Cheng H. The dorsolateral prefrontal cortex is selectively involved in chemotherapy-related cognitive impairment in breast cancer patients with different hormone receptor expression. Am J Cancer Res. 2019;9(8):1776–85.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Piccirillo JF, Hardin FM, Nicklaus J, Kallogjeri D, Wilson M, Ma CX, Coalson RS, Shimony J, Schlaggar BL. Cognitive impairment after chemotherapy related to atypical network architecture for executive control. Oncology. 2015;88(6):360–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Kesler SR, Adams M, Packer M, Rao V, Henneghan AM, Blayney DW, Palesh O. Disrupted brain network functional dynamics and hyper-correlation of structural and functional connectome topology in patients with breast cancer prior to treatment. Brain & Behavior. 2017;7(3):e00643.

    Article  Google Scholar 

  12. 12.

    Kesler SR, Blayney DW. Neurotoxic effects of Anthracycline- vs Nonanthracycline-based chemotherapy on cognition in breast Cancer survivors. JAMA Oncol. 2016;2(2):185–92.

    PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Haken H. Self-organization of brain function. Scholarpedia. 2008;3(5):2555.

    Article  Google Scholar 

  14. 14.

    Segall JM, Allen EA, Jung RE, Erhardt EB, Arja SK, Kiehl K, Calhoun VD. Correspondence between structure and function in the human brain at rest. Front Neuroinform. 2012;6(6):10.

    PubMed  PubMed Central  Google Scholar 

  15. 15.

    Biswal B, Yetkin FZ, Haughton VM, Hyde JS. Functional connectivity in the motor cortex of resting human brain using echo-planar mri. Magn Reson Med. 1995;34(4):537–41.

  16. 16.

    Barkhof F, Haller S, Rombouts SA. Resting-state functional MR imaging: a new window to the brain. Radiology. 2014;272(1):29–49.

    PubMed  Article  Google Scholar 

  17. 17.

    Bharti K, Suppa A, Pietracupa S, Upadhyay N, Giannì C, Leodori G, Di Biasio F, Modugno N, Petsas N, Grillea G. Aberrant functional connectivity in patients with Parkinson's disease and freezing of gait: a within- and between-network analysis. Brain Imaging Behav. 2019. https://doi.org/10.1007/s11682-019-00085-9.

  18. 18.

    Petrella JR, Sheldon FC, Prince SE, Calhoun VD, Doraiswamy PM. Default mode network connectivity in stable vs progressive mild cognitive impairment. Neurology. 2011;76(6):511–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Schumacher J, Firbank M, Gallagher P. Functional connectivity in dementia with Lewy bodies: a within- and between-network analysis. Hum Brain Mapp. 2018;39(3):1118–29.

    PubMed  Article  Google Scholar 

  20. 20.

    Kesler SR. Default mode network as a potential biomarker of chemotherapy-related brain injury. Neurobiol Aging. 2014;35(Suppl 2):S11–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Miao H, Chen X, Yan Y, He X, Hu S, Kong J, Wu M, Wei Y, Zhou Y, Wang L, Wang K, Qiu B. Functional connectivity change of brain default mode network in breast cancer patients after chemotherapy. Neuroradiology. 2016;58(9):921–8.

    PubMed  Article  PubMed Central  Google Scholar 

  22. 22.

    Feng Y, Tuluhong D, Shi Z, Zheng LJ, Chen T, Lu GM, Wang S, Zhang LJ. Postchemotherapy hippocampal functional connectivity patterns in patients with breast cancer: a longitudinal resting state functional MR imaging study. Brain Imaging Behav. 2019. https://doi.org/10.1007/s11682-019-00067-x.

  23. 23.

    Chen BT, Jin T, Patel SK, Ye N, Ma H, Wong CW, Rockne RC, Root JC, Saykin AJ, Ahles TA. Intrinsic brain activity changes associated with adjuvant chemotherapy in older women with breast cancer: a pilot longitudinal study. Breast Cancer Res Treat. 2019;176(1):181–9.

    PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Feng Y, Zhang XD, Zheng G, Zhang LJ. Chemotherapy-induced brain changes in breast cancer survivors: evaluation with multimodality magnetic resonance imaging. Brain Imaging Behav. 2019;13(6):1799–814.

    PubMed  Article  PubMed Central  Google Scholar 

  25. 25.

    Wang RY, Zhou JH, Huang YC, Yang YR. Reliability of the Chinese version of the trail making test and Stroop color and word test among older adults. Int J Gerontol. 2018;12(4):336–9.

    Article  Google Scholar 

  26. 26.

    Li YO, Adali T, Calhoun VD. Estimating the number of independent components for functional magnetic resonance imaging data. Hum Brain Mapp. 2007;28(11):1251–66.

    PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Greicius MD. Decoding subject-driven cognitive states with whole-brain connectivity patterns. Cereb Cortex. 2012;22(1):158–65.

    PubMed  Article  PubMed Central  Google Scholar 

  28. 28.

    Damoiseaux JS, Rombouts SARB, Barkhof F, Scheltens P, Stam CJ, Smith SM, Beckmann CF. Consistent resting-state networks across healthy subjects. Proc Natl Acad Sci U S A. 2006;103(37):13848–53.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Andrews-Hanna JR, Reidler JS, Sepulcre J, Poulin R, Buckner RL. Functional-anatomic fractionation of the Brain's default network. Neuron. 2010;65(4):550–62.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    Andrews-Hanna JR, Smallwood J, Spreng RN. The default network and self-generated thought: component processes, dynamic control, and clinical relevance. Ann N Y Acad Sci. 2014;1316:29–52.

  31. 31.

    Ji L, Pearlson GD, Hawkins KA, Steffens DC, Hua G, Wang L. A new measure for neural compensation is positively correlated with working memory and gait speed. Front Aging Neurosci. 2018;10:71.

    PubMed  PubMed Central  Article  Google Scholar 

  32. 32.

    Hampson JP, Zick SM, Khabir T, Wright BD, Harris RE. Altered resting brain connectivity in persistent cancer related fatigue. Neuroimage Clin. 2015;8:305–13.

    PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Buzsáki G, Draguhn A. Neuronal oscillations in cortical networks. Science. 2004;304(5679):1926–9.

  34. 34.

    Fox MD, Snyder AZ, Vincent JL, Corbetta M, Essen DCV, Raichle ME. The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc Natl Acad Sci U S A. 2005;102(27):9673–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Mennes M, Kelly C, Zuo XN, Martino AD, Biswal BB, Castellanos FX, Milham MP. Inter-individual differences in resting-state functional connectivity predict task-induced BOLD activity. Neuroimage. 2010;50(4):1690–701.

    PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Ferguson RJ, Mcdonald BC, Saykin AJ, Ahles TA. Brain structure and function differences in monozygotic twins: possible effects of breast Cancer chemotherapy. J Clin Oncol. 2007;25(25):3866–70.

    PubMed  PubMed Central  Article  Google Scholar 

  37. 37.

    Kesler SR, Bennett FC, Mahaffey ML, Spiegel D. Regional brain activation during verbal declarative memory in metastatic breast cancer. Clin Cancer Res. 2009;15(21):6665–73.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    Menning S, Ruiter MBD, Veltman DJ, Koppelmans V, Kirschbaum C, Boogerd W, Reneman L, Schagen SB. Multimodal MRI and cognitive function in patients with breast cancer prior to adjuvant treatment — the role of fatigue. Neuroimage Clin. 2015;7:547–54. https://doi.org/10.1016/j.nicl.2015.02.005.

    Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Pomykala KL, Ruiter MB, Deprez S, McDonald BC, Silverman DH. Integrating imaging findings in evaluating the post-chemotherapy brain. Brain Imaging Behav. 2013;7:436–52. https://doi.org/10.1007/s11682-013-9239-y.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    López Zunini RA, Scherling C, Wallis N, Collins B, MacKenzie J, Bielajew C, Smith AM. Differences in verbal memory retrieval in breast cancer chemotherapy patients compared to healthy controls: a prospective fMRI study. Brain Imaging Behav. 2013;7(4):460–77.

    PubMed  Article  PubMed Central  Google Scholar 

  41. 41.

    Churchill NW, Cimprich B, Askren MK, Reuter-Lorenz PA, Jung MS, Peltier S, Berman MG. Scale-free brain dynamics under physical and psychological distress: pre-treatment effects in women diagnosed with breast cancer. Hum Brain Mapp. 2015;36(3):1077–92.

    PubMed  Article  PubMed Central  Google Scholar 

  42. 42.

    Liao W, Mantini D, Zhang Z, Pan Z, Ding J, Gong Q, Yang Y, Chen H. Evaluating the effective connectivity of resting state networks using conditional granger causality. Biol Cybern. 2010;102(1):57–69.

    PubMed  Article  PubMed Central  Google Scholar 

  43. 43.

    Sridharan D, Levitin DJ, Menon V. A critical role for the right fronto-insular cortex in switching between central-executive and default-mode networks. Proceedings Nat Acad Ences. 2008;105(34):12569–74.

    CAS  Article  Google Scholar 

  44. 44.

    Dumas JA, Makarewicz J, Schaubhut GJ, Devins R, Albert K, Dittus K, Newhouse PA. Chemotherapy altered brain functional connectivity in women with breast cancer: a pilot study. Brain Imaging Behav. 2013;7(4):524–32.

    PubMed  Article  PubMed Central  Google Scholar 

  45. 45.

    Lepage C, Smith AM, Moreau J, Barlow-Krelina E, Wallis N, Collins B, MacKenzie J, Scherling C. A prospective study of grey matter and cognitive function alterations in chemotherapy-treated breast cancer patients. Springerplus. 2014;3:444.

    PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    McDonald BC, Conroy SK, Ahles TA, West JD, Saykin AJ. Alterations in brain activation during working memory processing associated with breast cancer and treatment: a prospective functional magnetic resonance imaging study. J Clin Oncol. 2012;30(20):2500–8.

    PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Ruiter MBD, Reneman L, Boogerd W, Veltman DJ, Caan M, Douaud G, Lavini C, Linn SC, Boven E, FSAMV D. Late effects of high-dose adjuvant chemotherapy on white and gray matter in breast cancer survivors: Converging results from multimodal magnetic resonance imaging. Hum Brain Mapp. 2012;33(aop):2971–83.

    PubMed  Article  PubMed Central  Google Scholar 

  48. 48.

    Ahles TA, Root JC, Ryan EL. Cancer- and cancer treatment-associated cognitive change: an update on the state of the science. J Clin Oncol. 2012;30(30):3675–86.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  49. 49.

    Berman MG, Askren MK, Mi SJ, Therrien B, Cimprich B. Pretreatment worry and neurocognitive responses in women with breast cancer. Health Psychol. 2014;33(3):222–31.

    PubMed  Article  PubMed Central  Google Scholar 

  50. 50.

    Rasgon NL. Differences in regional brain metabolism associated with specific formulations of hormone therapy in postmenopausal women at risk for AD. Psychoneuroendocrinol. 2011;36(4):502–13.

    Article  CAS  Google Scholar 

  51. 51.

    Reynolds CA, Gatz M, Prince JA, Berg S, Pedersen NL. Serum lipid levels and cognitive change in late life. J Am Geriatr Soc. 2010;58(3):501–9.

  52. 52.

    Mayeux R. Early Alzheimer's disease. New Engl J Med. 2010;363(18):1773.

    CAS  Google Scholar 

Download references

Acknowledgements

We would like to thank the patients and volunteers who participated in the studies. The authors gratefully acknowledge the altruism of all participants and contributions to this research.

Funding

This study was funded by the grants from the Natural Scientific Foundation of China (81322020 and 81230032 to L.J.Z., 8180071661 to L.J.Z.).

Author information

Affiliations

Authors

Contributions

YF conducted the subject’s MRI examination, participated in the interpretation and analysis of data and drafted the manuscript. YFW participated in the interpretation of the data and the revision of the manuscript. ZS participated in the MRI examination and the collection of clinical data. LJZ recruited chemotherapy patients and normal control population. WH contributed to the statistical analysis. LJZ conceived of the study, participated in its design, and was responsible for the preparation of the manuscript. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Long Jiang Zhang.

Ethics declarations

Ethics approval

Ethical approval was obtained by the Medical Ethics Committee of Jinling Hospital (2016NZGKJ-069). Written informed consent was obtained from all subjects.

Consent for publication

All authors have read the manuscript and agreed with the submission.

Competing interests

All authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Feng, Y., Wang, Y.F., Zheng, L.J. et al. Network-level functional connectivity alterations in chemotherapy treated breast cancer patients: a longitudinal resting state functional MRI study. Cancer Imaging 20, 73 (2020). https://doi.org/10.1186/s40644-020-00355-6

Download citation

Keywords

  • Breast cancer
  • Chemotherapy
  • Cognitive impairment
  • Functional connectivity
  • Resting-state functional magnetic resonance imaging