NSC 17987

No. of Pages 10

30 August 2017 Please cite this article in press as: Wang X-H et al. Predicting clinical symptoms of attention deficit hyperactivity disorder based on temporal patterns between and within intrinsic connectivity networks. Neuroscience (2017), http://dx.doi.org/10.1016/j.neuroscience.2017.08.038 1

Neuroscience xxx (2017) xxx–xxx

4

PREDICTING CLINICAL SYMPTOMS OF ATTENTION DEFICIT HYPERACTIVITY DISORDER BASED ON TEMPORAL PATTERNS BETWEEN AND WITHIN INTRINSIC CONNECTIVITY NETWORKS

5

XUN-HENG WANG, a* YUN JIAO, b AND LIHUA LI a*

6 7

a College of Life Information Science and Instrument Engineering, Hangzhou Dianzi University, Hangzhou 310018, China

2 3

8 9 10

11

Key words: ADHD, clinical symptoms, intrinsic connectivity networks, temporal patterns, machine learning.

b

Jiangsu Key Laboratory of Molecular and Functional Imaging, Department of Radiology, Zhongda Hospital, Medical School of Southeast University, Nanjing 210009, China

12

Abstract—Attention deficit hyperactivity disorder (ADHD) is a common brain disorder with high prevalence in schoolage children. Previously developed machine learningbased methods have discriminated patients with ADHD from normal controls by providing label information of the disease for individuals. Inattention and impulsivity are the two most significant clinical symptoms of ADHD. However, predicting clinical symptoms (i.e., inattention and impulsivity) is a challenging task based on neuroimaging data. The goal of this study is twofold: to build predictive models for clinical symptoms of ADHD based on resting-state fMRI and to mine brain networks for predictive patterns of inattention and impulsivity. To achieve this goal, a cohort of 74 boys with ADHD and a cohort of 69 age-matched normal controls were recruited from the ADHD-200 Consortium. Both structural and resting-state fMRI images were obtained for each participant. Temporal patterns between and within intrinsic connectivity networks (ICNs) were applied as raw features in the predictive models. Specifically, sample entropy was taken as an intra-ICN feature, and phase synchronization (PS) was used as an inter-ICN feature. The predictive models were based on the least absolute shrinkage and selectionator operator (LASSO) algorithm. The performance of the predictive model for inattention is r = 0.79 (p < 108), and the performance of the predictive model for impulsivity is r = 0.48 (p < 108). The ICN-related predictive patterns may provide valuable information for investigating the brain network mechanisms of ADHD. In summary, the predictive models for clinical symptoms could be beneficial for personalizing ADHD medications. Ó 2017 IBRO. Published by Elsevier Ltd. All rights reserved.

*Corresponding authors. Fax: +86-0571-87713528. E-mail addresses: [email protected] (X.-H. Wang), [email protected]. cn (L. Li). Abbreviations: ADHD, attention deficit hyperactivity disorder; aDMN, anterior default mode network; AN, auditory network; DAN, dorsal attention network; FD, frame displacement; fDMN, fusiform-related default mode network; ICNs, intrinsic connectivity networks; LASSO, least absolute shrinkage and selectionator operator; LN, language network; LVN, lateral visual network; MRI, magnetic resonance imaging; MVN, medial visual network; OVN, occipital visual network; PS, phase synchronization; RFPN, right frontoparietal network; SCN, sub-cortical network; SMN, sensorimotor network; SN, salience network; SRN, self-referential network; TLN, temporal lobe network. http://dx.doi.org/10.1016/j.neuroscience.2017.08.038 0306-4522/Ó 2017 IBRO. Published by Elsevier Ltd. All rights reserved. 1

INTRODUCTION

13

Attention deficit hyperactivity disorder (ADHD) is a highly prevalent brain disorder among school-age children. Approximately 5–10% of children and approximately 4% of adults suffer from ADHD worldwide (Biederman, 2005). Recent neuroimaging techniques identified altered structural and functional brain regions in ADHD patients. Smaller basal ganglia volumes were noted in children with ADHD than in normal controls based on structural magnetic resonance imaging (MRI) (Qiu et al., 2009). Decreased cortical thicknesses in attention- and executive-related cortical areas have also been reported in adults with ADHD based on structural MRI (Makris et al., 2007). Abnormal frontal- and cerebellar-related fractional anisotropy of white matter integrity was identified in ADHD-affected brains using diffusion tensor MRI (Ashtari et al., 2005). Altered anterior cingulate cortexrelated resting-state functional connectivity patterns were observed in adolescents with ADHD (Tian et al., 2006). Moreover, neuroimaging data may be able to classify patients with ADHD from healthy controls based on machine learning techniques (Zhu et al., 2008;Wang et al., 2013a). The ADHD-200 competition has released a large dataset for research on disease classification and mechanisms (Milham et al., 2012). Several methods have been applied to automatically diagnosis patients with ADHD based on this dataset (Eloyan et al., 2012). Support vector machine-based feature selection methods have been applied to diagnose ADHD through structural and functional MRI (Colby et al., 2012). Non-negative matrix factorization of multimodal datasets has also revealed different changes in ADHD-related default mode networks (Anderson et al., 2014). The above machine learning-based methods could discriminate patients with ADHD from normal controls by providing label information of the disease for individuals. Inattention and impulsivity are the two most prominent clinical symptoms of ADHD. However, predicting individual clinical symptoms (i.e., inattention and impulsivity) based on neuroimaging data is a challenging task. Investigating the relationships between neuroimaging data and clinical symptoms might be beneficial for

14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54

NSC 17987

No. of Pages 10

30 August 2017

2 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115

X.-H. Wang et al. / Neuroscience xxx (2017) xxx–xxx

personalized medication of ADHD patients. Several significant relationships have been identified between default mode network-related functional connectivity and ADHD-related behavior scores (Chabernaud et al., 2012). Significant correlations were observed between regional homogeneity and clinical symptoms (Wang et al., 2013a). Moreover, inter-network functional connectivity could predict clinical symptoms based on multisite ADHD datasets (Cai et al., 2015). Whole-brain functional connectivity may also serve as a potential biomarker for inattention according to a study of resting-state fMRI across a large number of patients with ADHD (Rosenberg et al., 2016). However, prediction of clinical symptoms based on neuroimaging data remains a challenging task for ADHD. Conventional voxel-wise or interregional measures based on MRI or fMRI typically possess tens of thousands of attributes, which are unfit for machine learning tasks. Thus, novel resting-state measures with an appropriate number of attributes are required to build a behavioral predictive model for ADHD. Intrinsic connectivity networks (ICNs) are spatiotemporally coherent patterns of resting-state functional connectivity and are consistent across healthy subjects (Damoiseaux et al., 2006). The spatial patterns of ICNs exhibited moderate-to-high test–retest reliability based on independent component analysis and dual regression (Zuo et al., 2010). The temporal patterns of ICNs exhibited good test–retest reliability based on intra-ICN complexity (Wang et al., 2013b). There are two types of ICN-related temporal patterns: the temporal measure within ICNs and the temporal measure between ICNs. Moreover, the amplitude of low-frequency fluctuations within the posterior default mode network was significantly related to inattentive scores (Wang and Li, 2015). The phase synchronization (PS) between ICNs could reflect different attention-related eyes-open/closed resting-states (Wang et al., 2015). Therefore, the temporal patterns within and between ICNs might be beneficial for building predictive models of clinical behaviors for ADHD based on neuroimaging data. In this paper, we hypothesized that the clinical symptoms of patients with ADHD could be predicted by the temporal patterns within and between ICNs, and the predictive powers of the temporal features within and between ICNs could be discovered by machine learning techniques. To test this hypothesis, a cohort of schoolage children with and without ADHD were recruited from a publicly available database. In the methods employed here, sample entropy was used as a temporal pattern within ICNs, whereas PS was used as a temporal pattern between ICNs. The combined features of sample entropy and PS were taken as inputs for the predictive model. The performance of the predictive model was evaluated using cross-validation in the results section. Moreover, the predictive weights of the combined features are reported in various illustrations. In the discussion section, we describe the performance of the predictive model and the relationships between temporal patterns of ICNs and behavior scores (i.e., inattention and impulsivity). In addition, we discuss the advantages and limitations of this study.

EXPERIMENTAL PROCEDURES

116

Participants and MRI protocols

117

A cohort of school-age boys with ADHD [74 subjects, mean age (11.98 ± 1.88)] and a cohort of age-matched male healthy controls [69 subjects, mean age (11.72 ± 1.8)] were recruited from the Institute of Mental Health, Peking University. All participants were initially diagnosed using the Computerized Diagnostic Interview Schedule IV (C-DIS-IV). In addition, one parent of each participant was asked to complete the test of the Schedule of Affective Disorders and Schizophrenia for Children-Present and Lifetime Version (KSADS-PL). The clinical symptoms of ADHD and normal controls were evaluated based on the ADHD Rating Scale (ADHD-RS) IV. Additional intelligence quotient (IQ) scores were obtained using the Wechsler Intelligence Scale for Chinese Children-Revised (WISCC-R). All participants included in this study had an IQ greater than 80. Each subject agreed to participate in this study and provided informed consent. The Research Ethics Review Board of Institute of Mental Health, Peking University, approved this study. The data usage agreement was provided by the functional 1000 connectome project. All data can be publicly obtained from the ADHD-200 Consortium (http://fcon_1000.projects.nitrc.org/indi/ adhd200/) (Milham et al., 2012). The demographic variables are listed in Table 1. For each participant, structural and resting-state functional brain images were collected from a SIEMENS 3T MRI scanner at Peking University. High-resolution structural brain images were collected using T1weighted sequences (magnetization prepared rapid acquisition gradient echo, MPRAGE). Resting-state brain images were based on standard functional sequences (echo-planar imaging, EPI). The detailed scan parameters were provided by the ADHD-200 Consortium website. All participants were requested to remain still and relaxed during the resting-state fMRI scan sessions. Each participant watched a black screen with a white cross during the scan.

118

Data preprocessing

156

The structural MRI images were preprocessed using the following steps: (1) deobliqued and reoriented to the fslfriendly space; (2) skull stripped; (3) segmented into different brain tissues (i.e., white matter, gray matter,

157

Table 1. Subjects’ demographic variables

Number of subjects Gender (male: female) Handless (R:L) Age (year) Full IQ Inattentive scores Impulsive scores

ADHD

Normal

p-value

74 74:0

69 69:0

1

74:0 11.98 ± 1.88 107.35 ± 12.31 28.12 ± 3.66 22.54 ± 6.04

69:0 11.72 ± 1.8 118.81 ± 14.26 16.04 ± 4.04 13.61 ± 3.55

1 0.398 n < =xmld el > inal y  predicted y; (4) increase (wj, wk)

236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251

255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295

in their joint least squares direction until a third predictor xm is identified that is most related to r; 5) traverse all xi to create the linear model. To achieve the LASSO solution, xi should be removed when wi reaches zero, and the least angle regression procedure should be repeated. The LASSO algorithm was performed using the glmnet package (https://cran.r-project.org/web/packages/glmnet/ index.html). The performances of the predictive models were evaluated by the Pearson correlation coefficients, the mean absolute error (MAE), and the root mean square error (RMSE) (Chai and Draxler, 2014) between the original values and the predicted values. The predictive models were trained and tested in a leave-one-out cross validation (LOOCV) (Zhu et al., 2008). Given n participants, the LOOCV procedure divided all samples into n pairs of training and testing samples. First, one sample was randomly selected as testing sample from all samples, leaving 142 samples as training samples. Then, the above step was repeated 143 times to traversal all samples. Finally, LOOCV produced 143 pairs of training subsets (142 samples) and a testing subset (1 sample).

296

RESULTS

318

297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317

Performance of the predictive models

319

Table 3 presents the performance of the predictive models for inattention based on different features. Table 4 presents the performance of the predictive models for impulsivity based on different features. Fig. 1 depicts the performance of the predictive models for clinical symptoms based on intra-ICN entropy and interICN PS. The predictive models with combined features of intra-ICN and inter-ICN perform better than models with inter-ICN features. The predictive model with the best performance for inattention has r = 0.81, MAE = 3.3, RMSE = 4.17, and p < 108 based on intra-ICN sample entropy. The predictive model with the best performance for impulsivity has r = 0.64, MAE = 4.08, RMSE = 5.19, and p < 108 based on intra-ICN sample entropy. Moreover, the PS-related

320

Table 3. Performances of the predictive models for inattention features

r

p-value

MAE

RMSE

Entropy FC PS Entropy and FC Entropy and PS

0.81 0.67 0.76 0.66 0.79

0 0 0 0 0

3.3 4.42 3.97 4.61 3.66

4.17 5.69 4.8 5.83 4.47

Table 4. Performances of the predictive models for impulsivity features

r

p-value

MAE

RMSE

Entropy FC PS Entropy and FC Entropy and PS

0.64 0.44 0.49 0.41 0.48

0

Predicting clinical symptoms of attention deficit hyperactivity disorder based on temporal patterns between and within intrinsic connectivity networks.

Attention deficit hyperactivity disorder (ADHD) is a common brain disorder with high prevalence in school-age children. Previously developed machine l...
2MB Sizes 3 Downloads 7 Views