Altered small-world properties of gray matter networks in breast cancer
© Hosseini et al.; licensee BioMed Central Ltd. 2012
Received: 1 December 2011
Accepted: 28 May 2012
Published: 28 May 2012
Skip to main content
© Hosseini et al.; licensee BioMed Central Ltd. 2012
Received: 1 December 2011
Accepted: 28 May 2012
Published: 28 May 2012
Breast cancer survivors, particularly those treated with chemotherapy, are at significantly increased risk for long-term cognitive and neurobiologic impairments. These deficits tend to involve skills that are subserved by distributed brain networks. Additionally, neuroimaging studies have shown a diffuse pattern of brain structure changes in chemotherapy-treated breast cancer survivors that might impact large-scale brain networks.
We therefore applied graph theoretical analysis to compare the gray matter structural networks of female breast cancer survivors with a history of chemotherapy treatment and healthy age and education matched female controls.
Results revealed reduced clustering coefficient and small-world index in the brain network of the breast cancer patients across a range of network densities. In addition, the network of the breast cancer group had less highly interactive nodes and reduced degree/centrality in the frontotemporal regions compared to controls, which may help explain the common impairments of memory and executive functioning among these patients.
These results suggest that breast cancer and chemotherapy may decrease regional connectivity as well as global network organization and integration, reducing efficiency of the network. To our knowledge, this is the first report of altered large-scale brain networks associated with breast cancer and chemotherapy.
Breast cancer is one of the most frequently diagnosed cancers and the leading cause of cancer death among females, accounting for 23% of the total cancer cases . In the last decade, several neuropsychological studies have shown the negative influence of breast cancer and chemotherapy on various cognitive skills with executive function and memory impairments being the most common [2, 3]. These deficits have been reported both prior and following chemotherapy with evidence showing increased and/or more severe cognitive changes in breast cancer patients treated with chemotherapy [4–9]. Neuroimaging studies corroborate these findings by showing changes in both brain structure and function associated with breast cancer and chemotherapy . However, it is currently unknown whether breast cancer and chemotherapy affect large-scale brain networks.
There are several lines of evidence suggesting that breast cancer may negatively impact whole brain networks. First, while the mechanisms by which breast cancer and its treatments affect cognitive function are largely unknown, possible candidates include neurotoxic effects of chemotherapy, oxidative damage, cytokine dysregulation and individual variation in genes related to neural repair and/or plasticity [2, 3, 11, 12]. These candidate mechanisms are likely to have diffuse effects on brain structure. Second, neuroimaging studies indicate that breast cancer survivors show altered brain structure, which would disrupt large-scale networks . Specifically, these patients demonstrate reduced gray matter in bilateral frontal, temporal, cerebellar, thalamic, and cingulate regions as well as decreased white matter integrity in corpus callosum, frontal, and temporal white matter tracts [14–17]. Third, the specific cognitive domains that tend to be most commonly impaired in breast cancer involve executive functions and memory, as noted above. These skills are known to be subserved by distributed, integrated neural networks . Finally, cognitive impairment following breast cancer often tends to be quite subtle , potentially suggesting a more diffuse brain injury. In the current study, we investigated whether breast cancer and chemotherapy are associated with alterations in large-scale structural brain networks.
Recent graph-theoretical analyses have consistently shown that brain structural networks in healthy individuals have small-world characteristics ; an architecture that has dense local clustering of connections between neighboring nodes with short path length between any pair of nodes due to the existence of relatively few long-range connections . These characteristics, shared by various biological systems, reflect a network that is simultaneously highly segregated and integrated and allows for higher, more efficient rates of information processing and learning than random networks .
Since small-world characteristics were described quantitatively for brain structural networks, there have been multiple graph-theoretical studies seeking to assess the structural correlation networks constructed from regional gray matter volume, cortical thickness and surface area [12, 20, 22, 23]. The unique feature of graph-theoretical analysis, compared with the more traditional univariate neuroimaging approaches, is that it can directly test the differences in topological parameters of the brain network such as small-worldness, highly connected hubs and regional network parameters. Whereas univariate neuroimaging approaches have typically shown limited correlations with cognitive function and dysfunction, network parameters may provide a more robust model of cognitive status [24, 25]. Recent graph-theoretical studies have illustrated an alteration of arrangements in structural correlation networks associated with normal aging, multiple sclerosis, Alzheimer’s disease, schizophrenia and epilepsy [26–30].
In the present study, we applied graph theoretical analyses to compare magnetic resonance imaging (MRI)-based gray matter correlation networks of female breast cancer patients treated with chemotherapy and female healthy controls. Considering the lines of evidence regarding a diffuse pattern of gray matter atrophy in breast cancer patients, we hypothesized that such alterations should be reflected in small-world characteristics of the brain structural correlation network. We also examined the between group differences in highly connected hubs as well as in regional network measures such as node betweenness and degree.
Demographic data for the breast cancer and healthy control groups
BC (N = 37)
CON (N = 38)
t or Chi Sq.
premorbid cognitive status*
MRI scanning was performed on a GE Discovery MR750 3.0 Tesla whole body scanner (GE Medical Systems, Milwaukee, WI). High-resolution T1-weighted images were acquired with 3D spoiled gradient echo pulse sequence using the following parameters: TR = 8.5 ms, TE = 3.396, TI = 400 ms, flip angle = 15, FOV = 220 mm, number of excitation =1, acquisition matrix = 256 x 192. Totally, 124 contiguous coronal slices were with in-plane resolution of 0.859 mm x 0.859 mm.
Image preprocessing was performed using Statistical Parametric Mapping 8 (SPM8; Wellcome Department of Cognitive Neurology, London, UK; http://www.fil.ion.ucl.ac.uk/spm/). The images were initially segmented into gray matter (GM), white matter, and cerebrospinal fluid images based on the ICBM Tissue Probabilistic Maps (http://www.loni.ucla.edu/ICBM/ICBM TissueProb.html). A study-specific a priori probability map of GM was created from the modulated spatially normalized segmented GM images using the Template-O-Matic (TOM8) toolbox . Next, the custom priors were affine-registered to the standard Montreal Neurological Institute (MNI) space.
Voxel-based morphometry (VBM) preprocessing steps were performed using the VBM8 toolbox for SPM8 (http://dbm.neuro.uni-jena.de/vbm/) and this involved segmentation of MR images into GM tissue segments using the standard unified segmented and non-linearly warping of the tissue segments to the GM study-specific customized template model . Images were then modulated to ensure that relative volumes of GM were preserved following the spatial normalization procedure. Sample homogeneity was checked to identify any outliers in the study population. Data for one participant in the BC group were excluded from the analysis because of a covariance below 2 SD and confirmation of visual motion artifact.
There are different nodal definition methods in brain network analysis. While the results might be affected by the choice of parcellation scheme, recent evidence showed that the results of between-group comparison remain intact regardless of the applied parcellation scheme . We generated 90 cortical and subcortical regions of interest (ROIs), excluding the cerebellum, from the Automated Anatomical Labeling (AAL) atlas using the WFU PickAtlas Toolbox . The ROIs were identical to those used in a previous graph analysis study of typical brain development by Fan and colleagues . These AAL ROIs were resliced to the same dimension as that of tissue segmented images obtained from the VBM preprocessing step. The ROIs were subsequently used to mask the individual modulated, normalized GM images and extract the average volume within each ROI using the REX toolbox (http://web.mit.edu/swg/software.htm). A linear regression analysis was performed at every ROI to remove the effects of age and total brain volume. The residuals of this regression were then substituted for the raw ROI volume values [22, 27, 36].
The extracted residual volumes of all 90 anatomical ROIs were used for construction of structural correlation networks. For each group, a 90 × 90 association matrix R was generated with each entry r ij defined as the Pearson correlation coefficient between residual volumes of regions i and j, across subjects. From each association matrix, a binary adjacency matrix A was derived where a ij was considered 1 if r ij was greater than a specific threshold and zero otherwise. The diagonal elements of the constructed association matrix were also set to zero. The resultant adjacency matrix represented a binary undirected graph G in which regions i and j were connected if g ij was unity. Therefore, a graph was constructed with N = 90 nodes (anatomical ROIs), with a network degree of E equal to number of edges (links), and a network density (cost) of D = E/[(N x (N-1))/2] representing the fraction of present connections to all possible connections. Since thresholding the association matrices of different groups at an absolute threshold results in networks with a different number of nodes (and degrees) that might influence the network measures and reduce interpretation of between group results , two approaches were implemented for thresholding the constructed association matrices based on previous studies [26–28]: (1) Thresholding the constructed association matrices at a minimum network density in which all nodes become fully connected in the brain networks of both groups; (2) Thresholding the constructed association matrices at a range of network densities and comparing the network topologies across that range.
The small-worldness of a complex network, as described above, has two key metrics: the clustering coefficient C and the characteristic path length L of the network. The clustering coefficient of a node is a measure of the number of edges that exist between its nearest neighbors. The clustering coefficient of a network is thus the average of clustering coefficients across nodes and is a measure of network segregation. The characteristic path length of a network is the average shortest path length between all pairs of nodes in the network and is the most commonly used measure of network integration . To evaluate the topology of the brain network, these parameters must be compared to the corresponding mean values of a random graph with the same number of nodes, total edges, and degree distribution as the network of interest [39, 40]. Thus, we obtained the small-worldness index of a network as [C/Crand/[L/Lrand where C rand and L rand are the mean clustering coefficient and the characteristic path length of the random network . In a small-world network, the clustering coefficient is significantly higher than that of random networks (C/C rand ratio greater than 1) while the characteristic path length is comparable to random networks (L/L rand ratio close to 1).
We also investigated the nodal characteristics of the constructed structural networks to identify differences in regional network measures between groups. Nodal betweenness centrality and nodal degree were calculated for each of the anatomical ROIs for the networks thresholded at minimum density with full-connectivity. Nodal betweenness centrality is defined as the fraction of all shortest paths in the network that pass through a given node and is used to detect important anatomical or functional connections. Nodes that bridge disparate parts of the network have a high betweenness centrality . On the other hand, nodal degree is defined as the number of connections that a node has with the rest of the network and is considered a measure of interaction of a node, structurally or functionally, with the network.
Hubs are crucial components for efficient communication in a network. Hubs are not only considered as important regulators of information flow but also play a key role in network resilience to insult . We considered a node as a hub if its betweenness centrality was at least 2SD higher than mean network centrality .
To test the significance of differences in overall inter-regional correlation of gray matter volume between networks, the correlation coefficients were converted to z values using Fisher’s r-to-z transform. This transformation resulted in values that were approximately normally distributed . Then, a two sample t-test was used to test the significance of the difference in mean overall correlation between groups.
In order to test the statistical significance of the between-group differences in network topology and regional network measures, a non-parametric permutation test with 1000 repetitions was used [26, 28]. In each repetition, the calculated residual volumes of each participant were randomly reassigned to one of the two groups so that each randomized group had the same number of subjects as the original groups. Then, an association matrix was obtained for each randomized group. The binary adjacency matrices were then estimated by thresholding the association matrices at a range of network densities. The network measures were then calculated for all the networks at each density. The differences in network measures between randomized groups were then calculated resulting in a permutation distribution of difference under the null hypothesis. The actual between-group difference in network measures was then placed in the corresponding permutation distribution and a two-tailed p-value was calculated based on its percentile position .
We used the Brain Connectivity Toolbox  for quantification of network measures in addition to in-house software, Graph Analysis Toolbox (http://nnl.stanford.edu/Tools.html), for comparing the structural networks. Brain Net Viewer (http://www.nitrc.org/projects/bnv/) was used for visualization of the graphs.
It should be noted that we could not statistically explore the relationship between cognitive-behavioral measures and the extracted network properties. It is because the extracted network measures are group-specific and can not be quantified for individuals. Therefore, we could not perform statistical analysis (e.g. correlation analysis) to relate network properties and cognitive-behavioral outcomes.
Network measures at minimum density of 0.184 with full connectivity
Mean clustering coefficient
Characteristic path length
Normalized path length
Below the minimum density of full connectivity (0.184), the BC network fragments sooner that the CON network and thus the normalized path length of the BC network increases relatively faster compared with that of CON network. This results in a relatively large drop in the difference in normalized path length between groups in densities below 0.184.
We applied graph theoretical analyses to compare gray matter volume correlation networks of female primary breast cancer (BC) survivors treated with chemotherapy and female healthy controls (CON). We found an alteration in small-world characteristics of the brain structural network in the BC group; an observation that confirms our hypothesis suggesting changes in large-scale brain network properties in BC survivors treated with chemotherapy. To our knowledge, this is the first report of altered large-scale network properties associated with BC and chemotherapy.
While the structural correlation networks of both groups followed a small-world organization, the results revealed a decrease in the small-world index of the BC group compared to controls. The characteristic path lengths of the BC and CON group networks were not different whereas the clustering coefficient was significantly smaller in the BC network than in the CON network. Decreased clustering coefficient and small-world index in the BC group was observed not only at the minimum network density but also over a range of density thresholds.
These findings suggest that the structural correlation networks of BC patients tend to have a more randomized configuration compared to CON group associated with weaker regional connectivity and disrupted global organization . Given that BC tends to be a disease associated with aging, altered network topology could make women with BC more vulnerable to age-related cognitive decline. Recent studies suggest that older BC patients show increased chemotherapy-related cognitive deficits [4, 5]. Additionally, this same profile of network topology (i.e. decreased clustering and small-world index) has been associated with Alzheimer’s dementia . The interface between aging, breast cancer and cognitive decline will be a critical area of continued investigation.
The observed lack of densely interconnected neighbors in the network of BC group is also supported by the lower overall correlation strength in the association matrix of BC group compared to CON group. Together, the results of network-level analysis corroborates previous structural neuroimaging findings that have demonstrated a diffuse pattern of atrophy in gray matter volumes in BC patients [15, 16] and further suggests an alteration in the coordinated patterns of brain morphology in the structural network of these patients.
One potential mechanism underlying this network-level alteration is white matter damage. There is a large body of evidence suggesting that global gray matter atrophy is associated with focal and global white matter damage due in part to the transection of axons and subsequent retrograde neuronal loss [42, 43]. Animal studies show that chemotherapy suppresses neural progenitor cell proliferation responsible for white-matter tract integrity and cortico-cortical connections [44–48]. Thus, the observed network-level alteration in structural correlation network of BC patients might arise from neurotoxic effects of chemotherapy on cortico-cortical connections. This idea is also supported by a recent diffusion weighted imaging study involving chemotherapy-treated BC survivors that reported a diffuse pattern of microstructural white matter damage .
High degree/betweenness nodes were observed mainly in frontotemporal regions in the CON network and in parietal, cingulate, and occipital regions in the BC network. Nodes with high degree/centrality in the structural network identify regions that are highly interactive and have the potential to participate in a large number of functional interactions . The results indicate that the network of BC group has less highly-interactive nodes in the frontotemporal regions.
Reduction in gray and white matter volumes in frontotemporal regions of patients with BC before and after treatment has been shown in a recent study . A number of functional neuroimaging studies have also reported increased activation in the frontotemporal network in these patients suggesting necessity of recruiting a compensatory neural mechanism to complete certain tasks, or decreased efficiency of these systems [10, 50, 51]. The observed changes in the centrality of the frontotemporal nodes in the network of BC patients is consistent with these previous results and may explain the impairment of these patients in memory and executive functioning.
There are several mechanisms that could explain the observed change in the distribution of highly influential nodes in BC group with less highly-interactive nodes in the prefrontal and temporal cortices. Certain genetic and immunologic factors have been proposed as important contributors to cognitive decline in patients with BC . For example, apolipoprotein E (APOE) is a complex protein that has an important role in neuronal repair and plasticity after injury . The relationship between this genotype and adverse cognitive performance and reduce mesial temporal (hippocampal) volume has been shown in long-term survivors of BC . Also, cytokines that become dysregulated during cancer and/or chemotherapy treatment have been shown to affect cognitive function in healthy adults . Cytokines have important roles in central nervous system functioning including modulation of dopamine and other neurotransmitters  that affect the frontotemporal network . The effects of genotypic variations and cytokine regulation on large-scale networks have not been established and require further investigation.
The BC group had significantly more post-menopausal women than the healthy comparison group. This is expected given that chemotherapy frequently causes early menopause  and tamoxifen can increase menopausal symptoms , both by blocking or reducing estrogen. Estrogen has a marked effect on brain function, particularly in frontotemporal regions and their concomitant executive and memory functions . Previous studies suggest that tamoxifen, in addition to chemotherapy, can have detrimental effects on cognitive and neurobiologic outcomes in BC survivors [59, 60]. The present findings may indicate potential secondary neuroendocrine effects on large-scale brain networks.
There was a difference in the number and distribution of hubs between the BC and CON groups. We found six network hubs in the CON group mainly across frontal and temporal regions compared to three hubs in the BC group primarily in parietal and cingulate areas. The identified hubs in the CON group are consistent with the results of previous graph-theoretical analysis involving healthy adults [26, 28]. The present findings indicate that the network of BC patients does not show expected prefrontal and temporal hubs suggesting network alterations involving these regions, again which are critical for the executive function and memory skills commonly disrupted in BC.
Additionally, the number of identified hubs in the BC network was half the number of hubs in the CON network. Since hubs are crucial components for efficient communication across the network, this observation suggests a less efficient information transfer in the network of BC patients. These results are consistent with our previous findings showing that BC survivors can often perform similarly to controls on cognitive tasks but require significantly more neural resources to do so . Standardized neurocognitive tests may lack sensitivity to these subtle changes  as most were designed to assess specific cognitive skills. Based on the present findings, measures of reaction time or parallel processing should be included in evaluations of BC-related cognitive outcome. Quantitative process measures might also be more effective as these assessments examine the efficiency of the patient’s performance strategies in addition to performance accuracy .
This study has several limitations. First, this was a cross-sectional study and therefore, we could not directly test the dissociated effects of cancer and chemotherapy on network measures. Future studies can address this issue through longitudinal evaluations of network measures in patients with BC. In addition, future studies might investigate whether the alteration in gray matter network is due to any cancer (and not just breast cancer). It is also important to dissociate between cancer-related physiological changes and cancer effects on mental states. Second, the network measures were identified by calculating the correlations of gray matter volume across subjects. While this methodology is common for investigating between-group differences in brain structural network, it lacks the ability to explore the individual differences in network parameters. Specifically, we could not explore the effect of treatment period, regimen, age or other individual differences on network parameters or investigate the effects of altered network parameters on cognitive-behavioral outcome.
Despite these limitations, these findings increase our understanding of BC and chemotherapy-related cognitive impairment by demonstrating alterations in specific network properties. Survivors of BC treated with chemotherapy showed decreased regional connectivity (clustering coefficient) and global network organization (small-worldness) and integration (hubs) suggesting reduced robustness and efficiency of the network. A pattern of frontotemporal abnormality was also noted. These results contribute novel insights regarding the neurobiologic mechanisms underlying cognitive deficits in these patients and highlight critical areas for future research related to the potential vulnerability of BC patients to the neurologic effects of aging and other forms of structural damage.
This work was supported by the National Institutes of Health Director’s New Innovator Award [DP2 OD004445-01 to SK].