Next Article in Journal
Emerging Neurological and Psychobiological Aspects of COVID-19 Infection
Previous Article in Journal
Screening of Cognitive Changes in Adults with Intellectual Disabilities: A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Functional and Structural Connectome Features for Machine Learning Chemo-Brain Prediction in Women Treated for Breast Cancer with Chemotherapy

1
School of Medicine, Chang Gung University, Taoyuan 33302, Taiwan
2
Department of Psychiatry, Chang Gung Memorial Hospital, Chiayi 61363, Taiwan
3
Department of Medical Imaging and Radiological Sciences, Chang Gung University, Taoyuan 33302, Taiwan
4
Breast Medical Center, Cheng Ching Hospital Chung Kang Branch, Taichung 40764, Taiwan
5
Department of Radiology, Taichung Veterans General Hospital, Taichung 40705, Taiwan
6
College of Medicine, China Medical University, Taichung 406040, Taiwan
7
Medical Imaging Research Center, Institute for Radiological Research, Chang Gung University and Chang Gung Memorial Hospital at Linkou, Taoyuan 33302, Taiwan
*
Author to whom correspondence should be addressed.
Submission received: 22 October 2020 / Revised: 7 November 2020 / Accepted: 11 November 2020 / Published: 12 November 2020

Abstract

:
Breast cancer is the leading cancer among women worldwide, and a high number of breast cancer patients are struggling with psychological and cognitive disorders. In this study, we aim to use machine learning models to discriminate between chemo-brain participants and healthy controls (HCs) using connectomes (connectivity matrices) and topological coefficients. Nineteen female post-chemotherapy breast cancer (BC) survivors and 20 female HCs were recruited for this study. Participants in both groups received resting-state functional magnetic resonance imaging (rs-fMRI) and generalized q-sampling imaging (GQI). Logistic regression (LR), decision tree classifier (CART), and xgboost (XGB) were the models we adopted for classification. In connectome analysis, LR achieved an accuracy of 79.49% with the functional connectomes and an accuracy of 71.05% with the structural connectomes. In the topological coefficient analysis, accuracies of 87.18%, 82.05%, and 83.78% were obtained by the functional global efficiency with CART, the functional global efficiency with XGB, and the structural transitivity with CART, respectively. The areas under the curves (AUCs) were 0.93, 0.94, 0.87, 0.88, and 0.84, respectively. Our study showed the discriminating ability of functional connectomes, structural connectomes, and global efficiency. We hope our findings can contribute to an understanding of the chemo brain and the establishment of a clinical system for tracking chemo brain.

1. Introduction

Chemotherapy is a widely used, long-established therapy for cancers; however, it is often associated with long-term toxicity inducing multiple side effects. Cognitive impairment in post-chemotherapy patients has been a serious issue in recent years with the increase in breast cancer survival rate and the number of breast cancer survivors [1,2,3,4]. A term, chemo brain, had been proposed and used to refer to a brain damaged by chemotherapy. Chemotherapy-induced cognitive impairment indicated by impaired cognitive functions such as memory, concentration, multiple operation, executive functions, processing speed, reaction time abilities, and word retrieval, is common among breast, lung, prostate, and ovarian cancer survivors who have received chemotherapeutic agents [3,5,6]. The cognitive impairment in patients can affect work performance and social relationships [6]. A recent study [7], where the daily life of 42 post-chemotherapy participants was evaluated, indicated the requirement of independent diagnostics for cognitive impairment. Although neuropsychological testing is the current standard, neuroimaging has the potential for diagnosing chemo-brain [8]. The National Cancer Institute (NCI) has also called for increased use of neuroscience-based methods. To accomplish this purpose, we conducted a study [9] using several MRI features, including generalized q-sampling imaging (GQI) [10] and resting-state functional magnetic resonance imaging (rs-fMRI) indices, for voxelwise analysis and regional summation analysis to discriminate chemo-brain participants from healthy controls (HCs).
Magnetic resonance imaging (MRI) methods are commonly applied in human brain studies. Blood oxygenation level-dependent functional magnetic resonance imaging (fMRI) has been widely used to explore cerebral function. fMRI is an efficient tool that can demonstrate functional changes in the human brain, and fMRI can be separated into task-based fMRI and resting-state fMRI (rs-fMRI). Compared with task-based fMRI, rs-MRI simplifies the experimental design by instructing the patient to relax and clear the mind and does not require them to follow certain task instructions. Functional imaging can aid in portraying the impacts of breast cancer and treatment.
Diffusion MRI is increasingly applied in the study of many white matter disorders. For example diffusion tensor imaging (DTI) can provide a non-invasive method to investigate brain microstructure and to measure the white matter tract integrity of anatomical connectivity [11,12,13]. However, DTI is unable to detect the crossing or branching patterns of complex regions, and DTI which is based on the Gaussian diffusion and monoexponential b-value dependence reflects the weighted average of all compartments even though the partial volumes of different diffusion compartments may vary [14,15,16]. To better characterize the complicated fiber patterns and distinguish fiber orientations, generalized q-sampling imaging (GQI) have been proposed, providing an opportunity for more accurate, higher-order descriptions when compared to DTI [10]. GQI is a unique q-space reconstruction method derived from q-space imaging. GQI could be applied to a wide range of q-space datasets for a more accurate and sophisticated diffusion MR approach.
Graph theory is a mathematical analysis that uses nodes and edges to represent objects and connections between objects [17]. Thus, researchers can apply graph theoretical analysis for denoting the complex connections of the human brain, which are believed to have segregated as well as integrated functions. In such studies, nodes usually represent brain regions, while edges represent connections between the brain regions. Furthermore, the edges could represent either functional connections (functional connectome; the gray matter network), or structural connections (structural connectome; the white matter network). More specifically, the gray matter network reflects the functional connectivity among brain areas, whereas the white matter network reflects the physical connections of nerve fibers between brain regions. According to previous graph theoretical studies, it is believed that “small-worldness” is a common characteristic of human brains that refers to the high clustering of nodes and short path lengths between nodes that compose an efficient network [18]. The topological coefficients of graph theory, including the mean clustering coefficient (C), normalized clustering coefficient (gamma, γ), local efficiency, characteristic path length (L), normalized characteristic path length (lambda, λ), global efficiency, transitivity, modularity, and the small-worldness index (sigma, σ), can be acquired by calculation.
A previous similar study [19], which conducted a machine learning classification between post-chemotherapy breast cancer survivors (C+), nonchemotherapy-treated breast cancer survivors (C−), and HCs using the connectomes of 19 independent regions of interest (ROIs) and support vector machines (SVMs), achieved an accuracy of 91.23% in C+ vs. C− and an accuracy of 90.74% in C+ vs. HC. fMRI connectivity data were also adopted by researchers who evaluated the accuracy of several models for the prediction of early multiple sclerosis [20], resulting in an accuracy of 85.7% by random forest (RF) and SVM. There was an enormous machine learning study [21] that also used structural connectomes as well as functional connectomes for prediction tasks. In particular, a total of 168 structural connectomes and 1013 functional connectomes were used to predict the cognitive and neuromotor development of preterm neonates and the autism spectrum category of the Autism Brain Imaging Data Exchange (ABIDE) database.
Continuing our previous study [22,23], we further processed our data using graph theoretical analysis and conducted another machine learning classification. Hence, the functional and structural connectomes (connectivity matrices) and their topological coefficients of the participants were used as features in this study. Note that we have two different types of connectomes in the current study, namely, functional connectomes calculated from the rs-fMRI data and structural connectomes calculated from the GQI data.
Logistic regression (LR), decision tree classifier (CART), and xgboost (XGB) were the three best classifiers in our former multiclassifier analysis, which all achieved an accuracy of 84%. Therefore, we prioritized and adopted these models in the current study. In the medical field, various studies have utilized LR and indicated the prominent efficacy of LR, such as predicting clinical drug response [24], differentiating between transition zone cancers and benign prostatic hyperplasia [25], and predicting the risk of a repeat surgical intervention for participants having septic arthritis of the hip [26]. The utilization of CART model has included predicting membrane protein types [27], differentiating brain pilocytic astrocytoma from glioblastoma [28], and diagnosis of biliary atresia in infants with jaundice [29]. The XGB models have won several machine learning competitions, and they have also been applied in many studies, such as coronary plaque characterization [30], predicting glucose metabolism disorder risk [31], and predicting invasive disease-free survival for early-stage breast cancer patients [32]. Theoretically, linear models such as LR were expected to attain the best performance on our high-dimensional dataset compared to ensemble models, RF and XGB. Although SVM, which is compatible with small sample size, should be an appropriate choice, the outcomes of it in our previous studies did not surpass or equal to the best ones; thus, we decided to dispose SVM but focus on the algorithms we used. Among them, LR is the most remarkable algorithm in our point of view. In the aspect of time consumption, LR took the least time to perform predictions, while XGB required the most time to do it. CART and RF are relatively unstable algorithms since their results are sometimes slightly different with the same configuration. Our purpose was to classify the connectomes and topological coefficients of participants into an HC group or a chemo-brain group using the LR, CART, and XGB models.

2. Materials and Methods

2.1. Participants and Clinical Characteristics

Nineteen post-chemotherapy female BC survivors (age range: 32 to 55 years, mean age: 43.8 ± 6.4 years; years of education: 13.9 ± 2.2; BC stage: I (n = 2), II (n = 14), III (n = 3)) and 20 HCs (age range: 43 to 55 years, mean age: 50.1 ± 2.5 years; years of education: 13.3 ± 2.3) were recruited from Taichung Veterans General Hospital in Taichung, Taiwan (Table 1). All BC survivors had histological confirmation of primary BC and received scans within 6 months after chemotherapy with standard chemotherapeutic agents, i.e., docetaxel and epirubicin. The exclusion criteria for patients with BC were being in a terminal stage of the disease (defined as a life expectancy of less than 1 year), having a treatment history for any cancer other than breast cancer or double cancer, undergoing radiation therapy before the investigation, having an indication of brain metastasis on postcontrast MRI findings or any known brain lesion from previous brain image scans, having any contraindication of MRI scans, and having a history of psychiatric or neurological illnesses or substance-use disorders [22]. The exclusion criteria for healthy controls were a history of psychiatric or neurological illnesses or substance-use disorder, a family history of major psychiatric or neurological illnesses, a current prescription of any psychiatric or psychotropic medications, current pregnancy or breastfeeding, and any contraindication of MRI scans [22]. In addition to chemotherapy, four of the chemo-brain participants received radiation therapy, one received hormone therapy, and nine received surgical therapy. The number of menopausal women in the patients and controls were 5:5. All BC survivors were confirmed with no brain or other metastases screened by brain MRI, chest computed tomography, and bone scan. Regarding menopause status, among all participants, 5 HCs and 5 chemo-brain participants were postmenopausal. All data remained in the functional analysis, however, the data of 1 HC and 1 chemo-brain participant were not included in the structural analysis. As a result, the classification using structural connectomes and topological coefficients contained only 19 HCs and 18 chemo-brain participants, while all 20 HCs and 19 chemo-brain participants were retained in the functional classification.
This study was approved by the Institutional Review Board of Taichung Veterans General Hospital (No. SF14185A). The research assistant explained the research proposal to all participants so that they could understand the research purpose, process, and both their rights and interests, and all participants were compensated when they finished the study. The clinical physicians assessed the physiological status of all participants to ensure they could participate in the MRI examination. All participants participated in the study after providing written informed consent, and all research was performed in accordance with relevant guidelines and regulations. The neuropsychological assessment was performed by the clinical psychologist, and the subsequent MRI examination was performed by the radiologic technologist. The overall process took approximately 1.5 h in the same day. The datasets analyzed during the current study are available from the corresponding author on reasonable request.
We used objective psychological tests and subjective questionnaire to evaluate the cognitive function, mindfulness, and psychological trauma of the participants. The neuropsychological tests included the Mini-Mental State Examination (MMSE), Impact of Event Scale-Revised (IES-R), and Cognitive and Affective Mindfulness Scale-Revised (CAMS-R), and Functional Assessment of Cancer Therapy-Cognitive Function (FACT-Cog). The results of the neuropsychological assessments are shown in Table 1.

2.2. MRI Acquisition

A 1.5-T MRI scan (Magnetom Aera, Siemens Medical Systems, Erlangen, Germany) with a standard, 8-channel head coil was used for all the MRI examinations. A total of three examinations were performed for each participant in this study. 1. T1-weighted 3D volume images (magnetization-prepared rapid gradient-echo sequence; TR/TE/TI = 9.11 ms/1.77 ms/450 ms; flip angle = 7°; 124 axial slices; voxel size = 1.0 × 1.0 × 1.4 mm3); 2. rs-fMRI (256 volumes; gradient echo-echo planar imaging sequence; 33 axial slices per volume; TR/TE = 2000 ms/30 ms; flip angle = 90°; voxel size = 3.4 × 3.4 × 4.0 mm3; temporal resolution = 2 s); 3. GQI (single-shot, diffusion-weighted spin echo-planar imaging sequence; TR/TE = 7200 ms/107 ms; field of view = 256 mm; matrix = 128 × 128; slice thickness = 4 mm; voxel size = 2 × 2 × 4 mm3; b-values = 1000, and 2000 s/mm2 in 128 noncollinear directions and one null image (b0); number of excitations = 1).

2.3. MRI Image Analysis

We used Statistical Parametric Mapping 8 (SPM8, Wellcome Department of Cognitive Neurology, London, UK) based on MATLAB (Math-Work, Natick, MA, USA) for MRI image preprocessing. For the fMRI data, we first conducted slice-timing correction as well as motion correction to reduce the influence caused by slice-timing and motion. Second, we normalized the data to the Montreal Neurological Institute (MNI, Montreal, QC, Canada) standard space [33] and performed spatial smoothing to improve image quality. Finally, linear detrending and bandpass temporal filtering with the frequency range from 0.01 to 0.08 Hz were performed to remove physiological noise by using the functional connectivity toolbox (CONN, The Gabrieli Lab. McGovern Institute for Brain Research MIT, Cambridge, MA, USA). For GQI data, the FSL Eddy correction (FMRIB Software Library, Oxford, UK) was first applied, followed by motion correction, normalizing, spatial smoothing, and normalizing. Eventually, normalized quantitative anisotropy (NQA) was calculated using DSI Studio (National Taiwan University, Taipei, Taiwan) for the subsequent analysis [10].

2.4. Graph Theoretical Analysis

For the functional connectome, the functional connectivity toolbox (CONN, The Gabrieli Lab. McGovern Institute for Brain Research MIT, Cambridge, MA, USA) was used, the whole brain was divided into 90 brain regions (45 per hemisphere) with an automated anatomical labeling (AAL) template, each of which was considered a node [34,35]. The brain functional connectivity between two nodes could be represented as an edge. The degree of a node is the number of edges connecting it to the rest of the network, and the degree can be used to characterize the edge distribution of all nodes in the network [36]. The functional MR image was registered to the T1-weighted image and then to the MNI space. The transformation matrix from functional image to MNI space was calculated by the transformation matrices created in the above two register processing steps. We normalized the resting-state functional images to the AAL template in MNI native space spatially, and the functional connectivity matrix was obtained after the functional connectivity analysis. Finally, the graph theoretical analysis was performed by using the functional connectivity matrix.
For the structural connectome, GQI was used to detect the direction of water molecule diffusion in the white matter non-invasively. The pathways of nerve fibers in the brain were reconstructed using fiber assignment by continuous tracking (FACT) with DSI studio. Network edges were established using FACT and the AAL templates, which divided the brain into 90 brain regions in MNI space. The number of virtual fibers connecting each pair of AAL regions was determined, resulting in a 90 × 90 structural connectivity matrix for each participant [17]. The network edges were defined as follows [1]:
E = Fiber count Fiber length × NQA
Finally, both the functional and structural topological properties of the complex network measures across the ranges of network density were obtained by graph theoretical algorithm. The network segregation was evaluated with the mean clustering coefficient (C), normalized clustering coefficient (gamma, γ), and local efficiency, and the network integration was evaluated with the characteristic path length (L), normalized characteristic path length (lambda, λ), and global efficiency. The transitivity, modularity, and small-worldness index (σ) were also evaluated [37]. The density represents the ratio of existing connections to all possible connections from the random network. The differences between groups in network measures below the network density depend on the number of individual networks that fragment in each group, therefore group comparisons below the density are not meaningful [38]. The minimum density value was defined by the limit density of the individual network not to be fragmented, and the maximum density value was defined by the density when the topology indices remained unchanged [39]. The network density range was calculated from 0.16 to 0.45 in 0.01 increments for the functional connectome and from 0.05 to 0.26 in 0.01 increments for the structural connectome.

2.5. Machine Learning Algorithms

For machine learning analysis, we used LR, CART, and XGB to classify the connectomes and topological coefficients into a chemo-brain group or an HC group using scikit-learn [40] in python. In addition to the scikit-learn library, we imported the xgboost library for XGB, which is a function called XGBClassifier in the library. Scikit-learn, a convenient, formidable library for machine learning, contains several preprocessing, feature selection methods, algorithms, and so on, which makes the machine learning tasks easier for researchers. LR is a linear classification with a sigmoid activation function, learning the weights for each feature during the training. The essential algorithm of LR is relatively simple and results in the rapid training speed of LR, which is why LR is one of the most fascinating models in terms of high-dimensional and/or sparse data. Nonetheless, LR may acquire inferior performance with low-dimensional data. CART is another prevalent model that sets a series of discriminative questions in each node during the training phase to classify samples into different groups. XGB, one of the most formidable machine learning models, is an ensemble method integrating multiple trees in a single model. Unlike RF, another iconic kind of ensemble method, trees (weak learners) in XGB will attempt to correct the error made by the former tree to build a stronger model. The connectomes and topological coefficients were standardized before being sent into leave-one-out cross-validation (LOOCV). Receiver operating characteristic (ROC) curves and area under the curves (AUCs) were recorded in figures that were generated with the matplotlib library [41].

3. Results

3.1. Connectome Analysis

In the connectome analysis, LR was the most discriminative model and attained a higher accuracy with functional connectomes (Figure 1) than with structural connectomes (Figure 2). An accuracy of 79.49% with an AUC of 0.93, a sensitivity of 0.80, and a specificity of 0.95 was achieved using functional connectomes, and an accuracy of 71.05% with an AUC of 0.94, a sensitivity of 0.82 and a specificity of 0.91 was achieved using structural connectomes. The figure was drawn using BrainNet [42].

3.2. Topological Coefficients Analysis

In the topological coefficient analysis, machine learning models could not distinguish participants with most features. Apart from most of the topological coefficients, functional global efficiency with CART and XGB and structural transitivity with CART reached accuracies of 87.18%, 82.05%, and 83.78%, with AUCs of 0.87, 0.88, and 0.84, respectively. The ROC curves and AUCs of the connectome analysis and topological analysis are shown in Figure 3.

4. Discussion

In the Figure 1, the coefficients derived from the attributes of LR, represent the independent coefficient of the features in the decision function. We found higher absolute values of the functional connection weighted coefficients in the LR model over the frontoparietal and occipital lobes in the cancer patients compared with the healthy controls. The altered brain regions may be associated with the dorsal attention network, and the results were consistent with previous functional neuroimaging studies [43,44]. In Figure 2, higher absolute values of the structural connection weighted coefficients in the LR model over the parietal and occipital lobes were also found in the cancer patients compared with the healthy controls, and the results were consistent with previous structural neuroimaging studies [45].
According to our results, we validated the discriminant ability of functional connectomes in differentiating HCs and chemo-brain participants as a previous study [19] had achieved. Comparing our results with two categories of connectomes, a higher accuracy of 79.49% was obtained with the functional connectomes, and a lower accuracy of 71.05% was obtained with the structural connectomes. However, utilization of structural connectomes could be found in a DL MRI study [46], where researchers used a DL network to predict the survival time of amyotrophic lateral sclerosis patients. They constructed four DL models based on: 1. eight clinical characteristics, 2. structural connectomes from MRI data, 3. cortical thickness and subcortical volume measurements from morphology MRI data, and 4. a combination of the previous three models. A highest accuracy of 84.4% was attained by the fourth model.
For the topological coefficient analysis, the global efficiency of the functional connectome and the transitivity of the structural connectome were the only two topological coefficients that showed an accuracy higher than 80%. The former reached 87.18% and 82.05% with CART and XGB, and the latter reached 83.78% with CART. According to our results, a subtle alteration in global efficiency and transitivity within chemo-brain participants might exist and be detectable with machine learning models. These results might plausibly suggest that there are disturbances in the optimal balance between functional/structural local segregation and global integration in chemo brains. In another study [47] using electroencephalogram (EEG), global efficiency was found to be reduced in young smokers compared with healthy nonsmokers, and the transitivity of major depressive disorder (MDD) participants appeared to decline when compared to HCs [48].
As we mentioned in the Introduction, the choice of machine learning algorithms was based on the outcome of our previous study (6), in which these three algorithms outperformed other algorithms. Additionally, the reason to use standardization remained the same since standardization had increased the classification accuracy. We also attempted to process feature sets to evaluate the performance of the models. For example, we once used a feature set that included all topological coefficients. However, the classification accuracy eventually decreased. From our point of view, this result may be attributed to the fact that only a few topological coefficients vary in chemo-brain participants, and the addition of low discriminative features detrimentally impacted the final performance of the model. Models using combined feature sets, were also considered, e.g., a model using functional and structural connectome. However, in most of the trials, a combined model did not acquire a result surpassing the model using either of them. Finally, we decided to only present results using feature separately.
By advanced computer technology and artificial intelligence, we were able to address numerous problems that may be difficult for humans to perform in an efficient way. For instance, drawing ROIs in images can be time-consuming and tedious, but now we can rely on AI to do the task for us. A recent study [49] provided an overview of the utilization of deep learning (DL)-based segmentation approaches on brain MRI, indicating several models in the past years and the advantages of using automated segmentation methods, such as alleviating intra-expert variabilities and interexpert variabilities. Apart from image segmentation, image reconstruction is another great example of benefits AI has brought to us. Researchers had proposed a DL model [50] that could reconstruct a contrast-enhanced brain MRI image with an image that was enhanced by merely one-tenth of the dosage allowing a reduced gadolinium dose. Using DL to reconstruct fully sampled MR images from undersampled images [29,51] is another popular topic since the undersampling of k-space usually represents a shorter scan time. Furthermore, DL could even generate pseudo-CT scans for attenuation correction in PET/MR imaging [52].
The interpretability of the machine learning model may be limited by the black-box nature of the classifiers. The results from machine learning depend on many a priori decisions regarding the different parameters chosen, and the conclusions may differ somewhat if a different set of parameters were used. It is recommended that future studies should aim to design a new variant to visualize the important brain source patterns and to realize the goal of interpretable machine learning. It is possible that the functional/structural connectome may be affected by behavior, aging, and other clinical issues [53,54], and this could be explored in future studies.
The first limitation was the lack of samples, which has been a common limitation in most MRI studies. However, in the connectome analysis it can be neglected since different values in the connectivity matrix and topological parameters in different density can be treated as one patient (one row in the final feature matrix). To isolate effects of chemotherapy from effects of cancer diagnosis (regardless of type of treatment) [55], we thought that a group of prechemotherapy BC survivors would be most appropriate for comparing pre- and post-chemotherapy participants, which was not achieved in the current study. The control group in the study is older than the patient group. Age may be the confounder and may result in bias close to null. Nonetheless, we can still use machine learning to distinguish patients from control group in the present study but the scale of the findings may underestimate. For our future work, we aim to increase the number of participants, including breast cancer patients and age-matched healthy controls, in our dataset and add prechemotherapy group and post treatment-but no-chemo group to address the disadvantages above.

5. Conclusions

Our study validated the alterations of the connectome in chemo-brain participants and identified the discriminating ability of connectome to machine learning. We are hopeful that this study could help the establishment of a clinically available model to track chemo brain in the future.

Author Contributions

Conceptualization, V.C.-H.C. and J.-C.W.; methodology, J.-C.W. and V.C.-H.C.; software, T.-Y.L.; validation, J.-C.W.; formal analysis, T.-Y.L.; investigation, V.C.-H.C. and J.-C.W.; resources, V.C.-H.C. and D.-C.Y.; data curation, V.C.-H.C., D.-C.Y. and J.-W.C.; writing—original draft preparation, T.-Y.L. and J.-C.W.; writing—review and editing, V.C.-H.C., D.-C.Y. and J.-W.C.; visualization, T.-Y.L. and J.-W.C.; supervision, J.-C.W.; project administration, J.-C.W. and V.C.-H.C.; funding acquisition, J.-C.W. and V.C.-H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the research programs MOST107-2221-E-182-054-MY3 and MOST106-2221-E-182-079, which were sponsored by the Ministry of Science and Technology, Taipei, Taiwan. This study was also supported by grants NMRPD1H0101~3 and NMRPD1G1321 of Chang Gung University, Taoyuan, Taiwan, and grants CORPG6G0101~3 and CORPG6G0121~3 of Chang Gung Memorial Hospital, Chiayi, Taiwan.

Conflicts of Interest

The authors report no conflict of interest in this work.

Abbreviations

AALautomated anatomical labeling
ABIDEAutism Brain Imaging Data Exchange
AUCareas under curve
BCbreast cancer
Cmean clustering coefficient
C+breast cancer survivors
C−nonchemotherapy-treated breast cancer survivors
CAMS-RCognitive and Affective Mindfulness Scale-Revised
CARTdecision tree classifier
DLdeep learning
EEGelectroencephalogram
FACTfiber assignment by continuous tracking
FACT-CogFunctional Assessment of Cancer Therapy-Cognitive Function
GQIgeneralized q-sampling imaging
HChealthy control
IES-RImpact of Event Scale–Revised
Lcharacteristic path length
LOOCVleave-one-out cross-validation
LRlogistic regression
MDDmajor depressive disorder
MMSEMini-Mental State Examination
MNIMontreal Neurological Institute
NQAnormalized quantitative anisotropy
RFrandom forest
ROIregion of interest
rs-fMRIresting-state functional magnetic resonance imaging
SPMStatistical Parametric Mapping
SVMsupport vector machine
XGBxgboost
γnormalized clustering coefficient
λnormalized characteristic path length
σsmall-worldness index

References

  1. Hislop, J.O. Yes, Virginia, Chemo Brain is Real. Clin. Breast Cancer 2015, 15, 87–89. [Google Scholar] [CrossRef]
  2. Moore, H.C. An overview of chemotherapy-related cognitive dysfunction, or ‘chemobrain’. Oncology 2014, 28, 797–804. [Google Scholar] [PubMed]
  3. Hermelink, K. Chemotherapy and Cognitive Function in Breast Cancer Patients: The So-Called Chemo Brain. J. Natl. Cancer Inst. Monogr. 2015, 2015, 67–69. [Google Scholar] [CrossRef] [Green Version]
  4. Bernstein, L.J.; McCreath, G.A.; Komeylian, Z.; Rich, J.B. Cognitive impairment in breast cancer survivors treated with chemotherapy depends on control group type and cognitive domains assessed: A multilevel meta-analysis. Neurosci. Biobehav. Rev. 2017, 83, 417–428. [Google Scholar] [CrossRef]
  5. Asher, A. Cognitive dysfunction among cancer survivors. Am. J. Phys. Med. Rehabil. 2011, 90 (Suppl. S1), S16–S26. [Google Scholar] [CrossRef]
  6. Reuter-Lorenz, P.A.; Cimprich, B. Cognitive function and breast cancer: Promise and potential insights from functional brain imaging. Breast Cancer Res. Treat. 2013, 137, 33–43. [Google Scholar] [CrossRef]
  7. Słowik, A.J.; Jabłoński, M.J.; Michałowska-Kaczmarczyk, A.M.; Jach, R. Evaluation of quality of life in women with breast cancer, with particular emphasis on sexual satisfaction, future perspectives and body image, depending on the method of surgery. Psychiatr. Pol. 2017, 51, 871–888. [Google Scholar] [CrossRef]
  8. Horowitz, T.S.; Suls, J.; Trevino, M. A Call for a Neuroscience Approach to Cancer-Related Cognitive Impairment. Trends Neurosci. 2018, 41, 493–496. [Google Scholar] [CrossRef]
  9. Chen, V.C.H.; Lin, T.Y.; Yeh, D.C.; Chai, J.W.; Weng, J.C. Predicting chemo-brain in breast cancer survivors using multiple MRI features and machine-learning. Magn. Reson. Med. 2019, 81, 3304–3313. [Google Scholar] [CrossRef]
  10. Yeh, F.C.; Wedeen, V.J.; Tseng, W.Y. Generalized q-sampling imaging. IEEE Trans. Med. Imaging 2010, 29, 1626–1635. [Google Scholar]
  11. Wise, T.; Radua, J.; Nortje, G.; Cleare, A.J.; Young, A.H.; Arnone, D. Voxel-Based Meta-Analytical Evidence of Structural Disconnectivity in Major Depression and Bipolar Disorder. Biol. Psychiatry 2016, 79, 293–302. [Google Scholar] [CrossRef] [Green Version]
  12. Srivastava, S.; Bhatia, M.S.; Bhargava, S.K.; Kumari, R.; Chandra, S. A Diffusion Tensor Imaging Study Using a Voxel-Based Analysis, Region-of-Interest Method to Analyze White Matter Abnormalities in First-Episode, Treatment-Naive Major Depressive Disorder. J. Neuropsychiatry Clin. Neurosci. 2016, 28, 131–137. [Google Scholar] [CrossRef]
  13. Beaulieu, C. The basis of anisotropic water diffusion in the nervous system—A technical review. NMR Biomed. 2002, 15, 435–455. [Google Scholar] [CrossRef]
  14. Roine, T.; Jeurissen, B.; Perrone, D.; Aelterman, J.; Leemans, A.; Philips, W.; Sijbers, J. Isotropic non-white matter partial volume effects in constrained spherical deconvolution. Front. Neuroinform. 2014, 8, 28. [Google Scholar] [CrossRef]
  15. Roine, T.; Jeurissen, B.; Perrone, D.; Aelterman, J.; Philips, W.; Leemans, A.; Sijbers, J. Informed constrained spherical deconvolution (iCSD). Med. Image Anal. 2015, 24, 269–281. [Google Scholar] [CrossRef]
  16. Vos, S.B.; Jones, D.K.; Viergever, M.A.; Leemans, A. Partial volume effect as a hidden covariate in DTI analyses. NeuroImage 2011, 55, 1566–1576. [Google Scholar] [CrossRef]
  17. Kesler, S.R.; Watson, C.L.; Blayney, D.W. Brain network alterations and vulnerability to simulated neurodegeneration in breast cancer. Neurobiol. Aging 2015, 36, 2429–2442. [Google Scholar] [CrossRef] [Green Version]
  18. Bassett, D.S.; Bullmore, E. Small-World Brain Networks. Neuroscientist 2006, 12, 512–523. [Google Scholar] [CrossRef]
  19. Kesler, S.R.; Wefel, J.S.; Hosseini, S.M.H.; Cheung, M.; Watson, C.L.; Hoeft, F. Default mode network connectivity distinguishes chemotherapy-treated breast cancer survivors from controls. Proc. Natl. Acad. Sci. USA 2013, 110, 11600–11605. [Google Scholar] [CrossRef] [Green Version]
  20. Saccà, V.; Sarica, A.; Novellino, F.; Barone, S.; Tallarico, T.; Filippelli, E.; Granata, A.; Chiriaco, C.; Bossio, R.B.; Valentino, P.; et al. Evaluation of machine learning algorithms performance for the prediction of early multiple sclerosis from resting-state FMRI connectivity data. Brain Imaging Behav. 2019, 13, 1103–1114. [Google Scholar] [CrossRef]
  21. Brown, C.J.; Miller, S.P.; Booth, B.G.; Zwicker, J.G.; Grunau, R.E.; Synnes, A.R.; Chau, V.; Hamarneh, G. Predictive connectome subnetwork extraction with anatomical and connectivity priors. Comput. Med. Imaging Graph. 2019, 71, 67–78. [Google Scholar] [CrossRef]
  22. Shen, C.-Y.; Chen, V.C.-H.; Yeh, D.-C.; Huang, S.-L.; Zhang, X.-R.; Chai, J.-W.; Huang, Y.-H.; Chou, M.-C.; Weng, J.-C. Association of functional dorsal attention network alterations with breast cancer and chemotherapy. Sci. Rep. 2019, 9, 104. [Google Scholar] [CrossRef]
  23. Li, T.-Y.; Chen, V.C.-H.; Yeh, D.-C.; Huang, S.-L.; Chen, C.-N.; Chai, J.-W.; Chen, C.C.-C.; Weng, J.-C. Investigation of chemotherapy-induced brain structural alterations in breast cancer patients with generalized q-sampling MRI and graph theoretical analysis. BMC Cancer 2018, 18, 1211. [Google Scholar] [CrossRef]
  24. Huang, H.-H.; Dai, J.-G.; Liang, Y. Clinical Drug Response Prediction by Using a Lq Penalized Network-Constrained Logistic Regression Method. Cell. Physiol. Biochem. 2018, 51, 2073–2084. [Google Scholar] [CrossRef]
  25. Iyama, Y.; Nakaura, T.; Katahira, K.; Nagayama, Y.; Oda, S.; Utsunomiya, D.; Yamashita, Y. Development and validation of a logistic regression model to distinguish transition zone cancers from benign prostatic hyperplasia on multi-parametric prostate MRI. Eur. Radiol. 2017, 27, 3600–3608. [Google Scholar] [CrossRef]
  26. Murphy, R.F.; Plumblee, L.; Barfield, W.B.; Murphy, J.S.; Fuerstenau, N.; Spence, D.D.; Kelly, D.M.; Dow, M.A.; Mooney, J.F. Septic Arthritis of the Hip—Risk Factors Associated With Secondary Surgery. J. Am. Acad. Orthop. Surg. 2019, 27, 321–326. [Google Scholar] [CrossRef]
  27. Sankari, E.S.; Manimegalai, D. Predicting membrane protein types using various decision tree classifiers based on various modes of general PseAAC for imbalanced datasets. J. Theor. Biol. 2017, 435, 208–217. [Google Scholar] [CrossRef]
  28. Dong, F.; Li, Q.; Xu, D.; Xiu, W.; Zeng, Q.; Zhu, X.; Xu, F.; Jiang, B.; Zhang, M. Differentiation between pilocytic astrocytoma and glioblastoma: A decision tree model using contrast-enhanced magnetic resonance imaging-derived quantitative radiomic features. Eur. Radiol. 2018, 29, 3968–3975. [Google Scholar] [CrossRef]
  29. Hyun, C.M.; Kim, H.P.; Lee, S.M.; Lee, S.; Seo, J.K. Deep learning for undersampled MRI reconstruction. Phys. Med. Biol. 2018, 63, 135007. [Google Scholar] [CrossRef]
  30. Masuda, T.; Nakaura, T.; Funama, Y.; Okimoto, T.; Sato, T.; Higaki, T.; Noda, N.; Imada, N.; Baba, Y.; Awai, K. Machine-learning integration of CT histogram analysis to evaluate the composition of atherosclerotic plaques: Validation with IB-IVUS. J. Cardiovasc. Comput. Tomogr. 2019, 13, 163–169. [Google Scholar] [CrossRef]
  31. Maeta, K.; Nishiyama, Y.; Fujibayashi, K.; Gunji, T.; Sasabe, N.; Iijima, K.; Naito, T.; Lopez, B.; Bian, J.; Billis, A.; et al. Prediction of Glucose Metabolism Disorder Risk Using a Machine Learning Algorithm: Pilot Study. JMIR Diabetes 2018, 3, e10212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Fu, B.; Liu, P.; Lin, J.; Deng, L.; Hu, K.; Zheng, H. Predicting Invasive Disease-Free Survival for Early Stage Breast Cancer Patients Using Follow-Up Clinical Data. IEEE Trans. Biomed. Eng. 2019, 66, 2053–2064. [Google Scholar] [CrossRef] [PubMed]
  33. Collins, D.L.; Neelin, P.; Peters, T.M.; Evans, A.C. Automatic 3D Intersubject Registration of MR Volumetric Data in Standardized Talairach Space. J. Comput. Assist. Tomogr. 1994, 18, 192–205. [Google Scholar] [CrossRef] [PubMed]
  34. Whitfield-Gabrieli, S.; Nieto-Castanon, A. Conn: A Functional Connectivity Toolbox for Correlated and Anticorrelated Brain Networks. Brain Connect. 2012, 2, 125–141. [Google Scholar] [CrossRef] [Green Version]
  35. Behzadi, Y.; Restom, K.; Liau, J.; Liu, T.T. A component based noise correction method (CompCor) for BOLD and perfusion based fMRI. Neuroimage 2007, 37, 90–101. [Google Scholar] [CrossRef] [Green Version]
  36. Bullmore, E.T.; Bassett, D.S. Brain graphs: Graphical models of the human brain connectome. Annu. Rev. Clin. Psychol. 2011, 7, 113–140. [Google Scholar] [CrossRef] [Green Version]
  37. Rubinov, M.; Sporns, O. Complex network measures of brain connectivity: Uses and interpretations. NeuroImage 2010, 52, 1059–1069. [Google Scholar] [CrossRef]
  38. Van Wijk, B.C.M.; Stam, C.J.; Daffertshofer, A. Comparing Brain Networks of Different Size and Connectivity Density Using Graph Theory. PLoS ONE 2010, 5, e13701. [Google Scholar] [CrossRef]
  39. Hosseini, S.M.H.; Hoeft, F.; Kesler, S.R. GAT: A Graph-Theoretical Analysis Toolbox for Analyzing Between-Group Differences in Large-Scale Structural and Functional Brain Networks. PLoS ONE 2012, 7, e40709. [Google Scholar] [CrossRef] [Green Version]
  40. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine Learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  41. Hunter, J.D. Matplotlib: A 2D Graphics Environment. Comput. Sci. Eng. 2007, 9, 90–95. [Google Scholar] [CrossRef]
  42. Xia, M.; Wang, J.; He, Y. BrainNet Viewer: A Network Visualization Tool for Human Brain Connectomics. PLoS ONE 2013, 8, e68910. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Chen, X.; He, X.; Tao, L.; Cheng, H.; Li, J.; Zhang, J.; Qiu, B.; Yu, Y.; Wang, K. The attention network changes in breast cancer patients receiving neoadjuvant chemotherapy: Evidence from an arterial spin labeling perfusion study. Sci. Rep. 2017, 7, 42684. [Google Scholar] [PubMed] [Green Version]
  44. McDonald, B.C.; Conroy, S.K.; Ahles, T.A.; West, J.D.; Saykin, A.J. 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, 2500–2508. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Deprez, S.; Amant, F.; Yigit, R.; Porke, K.; Verhoeven, J.; Stock, J.V.D.; Smeets, A.; Christiaens, M.-R.; Leemans, A.; Van Hecke, W.; et al. Chemotherapy-induced structural changes in cerebral white matter and its correlation with impaired cognitive functioning in breast cancer patients. Hum. Brain Mapp. 2011, 32, 480–493. [Google Scholar] [CrossRef] [PubMed]
  46. Van Der Burgh, H.K.; Schmidt, R.; Westeneng, H.-J.; De Reus, M.A.; Berg, L.H.V.D.; Heuvel, M.P.V.D. Deep learning predictions of survival based on MRI in amyotrophic lateral sclerosis. NeuroImage Clin. 2017, 13, 361–369. [Google Scholar]
  47. Su, S.; Yu, D.; Cheng, J.; Chen, Y.; Zhang, X.; Guan, Y.; Li, Y.; Bi, Y.; Xue, T.; Lu, X.; et al. Decreased Global Network Efficiency in Young Male Smoker: An EEG Study during the Resting State. Front. Psychol. 2017, 8, 1605. [Google Scholar] [CrossRef] [Green Version]
  48. Chen, V.C.-H.; Shen, C.-Y.; Liang, S.H.-Y.; Li, Z.-H.; Hsieh, M.-H.; Tyan, Y.-S.; Lu, M.-L.; Lee, Y.; McIntyre, R.S.; Weng, J.-C. Assessment of brain functional connectome alternations and correlation with depression and anxiety in major depressive disorders. PeerJ 2017, 5, 3147. [Google Scholar] [CrossRef]
  49. Akkus, Z.; Galimzianova, A.; Hoogi, A.; Rubin, D.L.; Erickson, B.J. Deep Learning for Brain MRI Segmentation: State of the Art and Future Directions. J. Digit. Imaging 2017, 30, 449–459. [Google Scholar] [CrossRef] [Green Version]
  50. Gong, E.; Pauly, J.M.; Wintermark, M.; Zaharchuk, G. Deep learning enables reduced gadolinium dose for contrast-enhanced brain MRI. J. Magn. Reson. Imaging 2018, 48, 330–340. [Google Scholar] [CrossRef]
  51. Tezcan, K.C.; Baumgartner, C.F.; Luechinger, R.; Pruessmann, K.P.; Konukoglu, E. MR Image Reconstruction Using Deep Density Priors. IEEE Trans. Med Imaging 2018, 38, 1633–1642. [Google Scholar] [CrossRef] [PubMed]
  52. Liu, F.; Jang, H.; Kijowski, R.; Bradshaw, T.; McMillan, A.B. Deep Learning MR Imaging–based Attenuation Correction for PET/MR Imaging. Radiology 2018, 286, 676–684. [Google Scholar] [CrossRef] [PubMed]
  53. Hsu, Y.-H.; Chen, V.C.-H.; Hsieh, C.-C.; Weng, Y.-P.; Hsu, Y.-T.; Hsiao, H.-P.; Wang, W.-K.; Chen, H.-M.; Weng, J.-C.; Wu, S.-I.; et al. Subjective and objective cognitive functioning among patients with breast cancer: Effects of chemotherapy and mood symptoms. Breast Cancer 2020. [Google Scholar] [CrossRef]
  54. Nugent, B.D.; Ren, D.; Bender, C.M.; Rosenzweig, M. The Impact of Age and Adjuvant Chemotherapy Modifications on Survival Among Black Women With Breast Cancer. Clin. Breast Cancer 2019, 19, 254–258. [Google Scholar] [CrossRef] [PubMed]
  55. Scherling, C.S.; Collins, B.; MacKenzie, J.; Bielajew, C.; Smith, A.M. Pre-Chemotherapy Differences in Visuospatial Working Memory in Breast Cancer Patients Compared to Controls: An fMRI Study. Front. Hum. Neurosci. 2011, 5, 122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The weighted coefficients of the functional connections in the logistic regression (LR) model. Warmer colors represent higher absolute values of the LR coefficients. A threshold of 0.02706 was set to reduce the connections for visual presentation. Coefficients derived from the attributes of LR, representing the independent coefficient of the features in the decision function. R: right; L: left; M: middle; A: anterior; P: posterior; S: superior; I: inferior; F: frontal; T: temporal P: parietal; O: occipital; G: gyrus.
Figure 1. The weighted coefficients of the functional connections in the logistic regression (LR) model. Warmer colors represent higher absolute values of the LR coefficients. A threshold of 0.02706 was set to reduce the connections for visual presentation. Coefficients derived from the attributes of LR, representing the independent coefficient of the features in the decision function. R: right; L: left; M: middle; A: anterior; P: posterior; S: superior; I: inferior; F: frontal; T: temporal P: parietal; O: occipital; G: gyrus.
Brainsci 10 00851 g001
Figure 2. The weighted coefficients of the structural connections in the LR model. Warmer colors represent higher absolute values of the LR coefficients. A threshold of 0.03105 was set to reduce the connections for visual presentation. R: right; L: left; M: middle; A: anterior; P: posterior; S: superior; I: inferior; F: frontal; T: temporal P: parietal; O: occipital; G: gyrus.
Figure 2. The weighted coefficients of the structural connections in the LR model. Warmer colors represent higher absolute values of the LR coefficients. A threshold of 0.03105 was set to reduce the connections for visual presentation. R: right; L: left; M: middle; A: anterior; P: posterior; S: superior; I: inferior; F: frontal; T: temporal P: parietal; O: occipital; G: gyrus.
Brainsci 10 00851 g002
Figure 3. The receiver operating characteristic (ROC) curves and areas under the curves (AUCs) of the predictions. (a) Functional connectome analysis and topological coefficient analysis (AUC: connectivity matrix-LR = 0.93; global efficiency-CART = 0.87; global efficiency = 0.88); (b) structural connectome analysis and topological coefficient analysis (AUC: connectivity matrix-LR = 0.94; transitivity-CART = 0.84).
Figure 3. The receiver operating characteristic (ROC) curves and areas under the curves (AUCs) of the predictions. (a) Functional connectome analysis and topological coefficient analysis (AUC: connectivity matrix-LR = 0.93; global efficiency-CART = 0.87; global efficiency = 0.88); (b) structural connectome analysis and topological coefficient analysis (AUC: connectivity matrix-LR = 0.94; transitivity-CART = 0.84).
Brainsci 10 00851 g003
Table 1. Demographic characteristics and summary of neuropsychological testing and questionnaire evaluation of the all participants in the study.
Table 1. Demographic characteristics and summary of neuropsychological testing and questionnaire evaluation of the all participants in the study.
CharacteristicBreast Cancer Patients after Chemotherapy (n = 19)Healthy Controls (n = 20)
MeanSDMinMaxMeanSDMinMaxp-ValueEffect Sizes
Age (years)44.846.56325549.952.6143550.004−1.03
Education (years)14.002.6191613.22.269160.2740.33
MMSE28.111.29263028.351.4625300.583−0.17
CAMS-R33.634.91244233.853.9427400.879−0.05
IES-R14.2623.910648.6012.570380.3660.30
FACT-Cog
Perceived cognitive impairments53.2610.49346759.456.8847700.035−0.70
Comments from others14.002.4771614.301.8411160.669−0.14
Perceived cognitive abilities17.215.7822618.005.490240.664−0.14
Impact on quality of life11.633.5241614.102.0510160.010−0.86
Total score102.3719.2453131112.5012.31901340.056−0.63
Breast cancer stagen% n%
000 N/AN/A
I210.53 N/AN/A
II1473.68 N/AN/A
III315.79 N/AN/A
IV00 N/AN/A
Chemotherapeutic drugs (Taxotere and Epirubicin)19100.00 N/AN/A
Radiation therapy421.05 N/AN/A
Hormonal treatment15.26 N/AN/A
Menopausal526.32 525.00
Abbreviations: MMSE: Mini-Mental State Examination; CAMS-R: Cognitive and Affective Mindfulness Scale-Revised; IES-R: Impact of Event Scale—Revised; FACT-Cog: Functional Assessment of Cancer Therapy-Cognitive Function; SD: standard deviation; N/A: not applicable.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, V.C.-H.; Lin, T.-Y.; Yeh, D.-C.; Chai, J.-W.; Weng, J.-C. Functional and Structural Connectome Features for Machine Learning Chemo-Brain Prediction in Women Treated for Breast Cancer with Chemotherapy. Brain Sci. 2020, 10, 851. https://0-doi-org.brum.beds.ac.uk/10.3390/brainsci10110851

AMA Style

Chen VC-H, Lin T-Y, Yeh D-C, Chai J-W, Weng J-C. Functional and Structural Connectome Features for Machine Learning Chemo-Brain Prediction in Women Treated for Breast Cancer with Chemotherapy. Brain Sciences. 2020; 10(11):851. https://0-doi-org.brum.beds.ac.uk/10.3390/brainsci10110851

Chicago/Turabian Style

Chen, Vincent Chin-Hung, Tung-Yeh Lin, Dah-Cherng Yeh, Jyh-Wen Chai, and Jun-Cheng Weng. 2020. "Functional and Structural Connectome Features for Machine Learning Chemo-Brain Prediction in Women Treated for Breast Cancer with Chemotherapy" Brain Sciences 10, no. 11: 851. https://0-doi-org.brum.beds.ac.uk/10.3390/brainsci10110851

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop