A prediction model of working memory across health and psychiatric disease using whole-brain functional connectivity

  1. Masahiro Yamashita
  2. Yujiro Yoshihara
  3. Ryuichiro Hashimoto
  4. Noriaki Yahata
  5. Naho Ichikawa
  6. Yuki Sakai
  7. Takashi Yamada
  8. Noriko Matsukawa
  9. Go Okada
  10. Saori C Tanaka
  11. Kiyoto Kasai
  12. Nobumasa Kato
  13. Yasumasa Okamoto
  14. Ben Seymour  Is a corresponding author
  15. Hidehiko Takahashi  Is a corresponding author
  16. Mitsuo Kawato  Is a corresponding author
  17. Hiroshi Imamizu  Is a corresponding author
  1. Advanced Telecommunications Research Institute International, Japan
  2. Kyoto University, Japan
  3. Showa University, Japan
  4. The University of Tokyo, Japan
  5. National Institute of Radiological Sciences, Japan
  6. Hiroshima University, Japan
  7. Kyoto Prefectural University of Medicine, Japan
  8. University of Cambridge, United Kingdom
  9. National Institute of Information and Communications Technology, Japan
6 figures, 2 tables and 3 additional files

Figures

Schematic diagram of model construction and generalization tests using independent datasets.

(A) Model was developed using a whole-brain resting state FC and a learning plateau of a letter 3-back task within healthy individuals from ATR dataset. (B) We applied the model to resting state FC patterns and predicted individual participant’s working memory ability. We first examined the external validity using an independent USA healthy dataset (HCP dataset: the upper flow chart in (B)). The predicted working memory ability was compared to actual working memory performance (visual-object N-back task and the NIH toolbox list sorting test). Then we examined the generalizability to a clinical population using a schizophrenia dataset (the lower flow chart in (B)). The predicted working memory ability was compared to actual working memory score measured by Digit sequencing test. (C) Using the multiple psychiatric diagnoses dataset, degree of working memory impairment for each diagnosis was predicted as differences from corresponding controls. The predicted impairments were validated by previous meta-analysis studies on digit-span across multiple diagnoses. Note that the HCP dataset’s task stimuli images are just illustration purpose and different from the original stimuli.

https://doi.org/10.7554/eLife.38844.002
Figure 2 with 1 supplement
Normative model of working memory ability (WMA).

Circle plot of networks and their connections in the model. Individual letter 3-back learning plateaus are predicted by a linear weighted summation of 16 FC values at 16 connections selected by a sparse linear regression algorithm. Connection thicknesses indicate contribution ratios (weight x FC at each connection). Connections are labeled ‘Positive/Negative (P/N)’ based on correlation coefficient signs with letter 3-back learning performances, whereas numbers indicate descending orders of contribution ratio. Each network’s color indicates relevance with working memory function based on BrainMap ICA (Laird et al., 2011); warmer colors indicate closer relevance to working memory function. See Table 1 for the networks connected by the selected 16 connections, and precise values of contribution ratio of each connection. Each network’s label and regions included in it are summarized in Supplementary file 2.

https://doi.org/10.7554/eLife.38844.003
Figure 2—figure supplement 1
Letter 3-back learning curves for each participant.

Data points were five-session moving averaged (e.g. training session one plots average d-prime from sessions 1 to 5).

https://doi.org/10.7554/eLife.38844.004
Figure 3 with 3 supplements
Generalizability to HCP dataset and schizophrenia dataset.

(A) Significant Spearman’s rank partial correlation between predicted letter 3-back learning performance and measured visual-object N-back accuracy while factoring out general fluid intelligence and head motion (ρ = 0.110, p = 0.0072). (B) Significant Pearson partial correlation between predicted letter 3-back performances and measured digit-sequencing scores while factoring out the composite BACS score and age (ρ = 0.248, p = 0.033).

https://doi.org/10.7554/eLife.38844.005
Figure 3—figure supplement 1
Spearman’s rank correlation matrix for HCP dataset.
https://doi.org/10.7554/eLife.38844.006
Figure 3—figure supplement 2
Pearson’s correlation matrices for schizophrenia samples.

(A) Correlation among schizophrenia patients (N = 58). (B) Correlation among schizophrenia patients and controls (N = 118).

https://doi.org/10.7554/eLife.38844.007
Figure 3—figure supplement 3
Distribution of BACS digit-sequencing score.

The count of healthy controls increased toward the maximum score (28), and the largest count was found for the maximum score. This suggests a ceiling effect in the healthy controls.

https://doi.org/10.7554/eLife.38844.008
Figure 4 with 1 supplement
Prediction of diagnosis-specific alterations of working memory ability.

(A) Predicted letter 3-back working memory ability for patients (N = 58, 77, 45, and 69 for SCZ, MDD, OCD, and ASD, respectively) and their age- and gender-matched healthy/typically developed controls (HC, N = 60, 62, 47, and 71) shown as kernel density. For illustration purposes, distribution of each control group was standardized to that of the ATR dataset, and the same linear transformation was applied to patients’ distributions. μ indicates mean value for each group. (B) Violin plots of Z-scores for predicted working memory ability alterations. White circles indicate medians. Box limits indicate 25th and 75th percentiles. Whiskers extend 1.5 times interquartile range from 25th and 75th percentiles. (C) Comparison of estimated effect sizes for working memory deficits. k indicates number of studies included in the meta-analyses (Forbes et al., 2009; Snyder, 2014; Snyder et al., 2015). Error bars indicate 95% confidence intervals.

https://doi.org/10.7554/eLife.38844.009
Figure 4—figure supplement 1
Effect sizes of IQ, digit-span, and predicted working memory ability.

k indicates number of studies included in the meta-analyses (Abramovitch et al., 2018; Ahern and Semkovska, 2017; Heinrichs and Zakzanis, 1998). Error bars indicate 95% confidence intervals.

https://doi.org/10.7554/eLife.38844.010
Figure 5 with 3 supplements
Accumulation of function connectivity differences exhibits diagnosis-specific working memory ability.

(A) Accumulation of averaged D-scores for all 16 connections. Bold black line indicates summation of contributions by all connections, corresponding to predicted working memory ability alteration. This figure shows how diagnosis-specific working memory impairment results from complex disturbances of multiple connections. Upper panel depicts two representative alteration patterns across diagnoses. While connection P1 commonly decreased working memory ability across diagnoses, connection N6 distinctly affected working memory ability (decrease in SCZ and MDD and increase in OCD and ASD). (B) Z-scores (normalized D-scores) for each diagnosis. Left asterisks and lines indicate significant differences in mean Z-scores between two diagnoses (p < 0.05, Bonferroni corrected). Vertical lines across horizontal bars indicate Z-scores averaged across connections. (C) Z-scores for connection that showed a significant effect of diagnosis. Connections were sorted by small p values of diagnosis effect (Kruskal-Wallis test, Q < 0.05, FDR corrected).

https://doi.org/10.7554/eLife.38844.013
Figure 5—figure supplement 1
Connections of non-significant effect of diagnosis.
https://doi.org/10.7554/eLife.38844.014
Figure 5—figure supplement 2
Distribution of magnitude of diagnosis effect on connectivity change between the connections in the model and other connections.

The horizontal axis shows chi-square values of a Kruskal-Wallis test for each connection. There was no significant difference between the two groups (two-sample Kolmogorov-Smirnov test, p = 0.30).

https://doi.org/10.7554/eLife.38844.015
Figure 5—figure supplement 3
D-scores related to network clusters.
https://doi.org/10.7554/eLife.38844.016
Author response image 1
Pearson’s correlation matrix of BACS sub-scores.

Tables

Table 1
Selected connections and their contribution to working memory ability.
https://doi.org/10.7554/eLife.38844.012
LabelConnection (rank)Contribution ratio [%]
Positive features
P1Left fronto-parietal network (1)(within-network)33.9%
P2Supplemental motor network (3)Primary sensorimotor network (hand) (11)15.4%
P3Middle frontal and parietal network (2)Lateral temporal network (6)1.9%
Negative features
N1Cingulo-opercular network (5)Midbrain (10)13.9%
N2Right fronto-parietal network (4)Midbrain (10)11.4%
N3Right fronto-parietal network (4)Superior parietal network (18)9.2%
N4Supplemental motor network (3)Orbitofrontal network (14)5.0%
N5Middle frontal and parietal network (2)Primary sensorimotor network (mouth) (9)3.0%
N6Lateral occipital network (7)Auditory (15)2.6%
N7Left fronto-parietal network (1)Midbrain (10)1.8%
N8Cerebellum (12)Auditory (15)1.3%
N9Left fronto-parietal network (1)Lateral occipital network (7)0.3%
N10Lateral occipital network (7)Primary sensorimotor network (hand) (11)0.3%
N11Lateral occipital network (7)Superior parietal network (18)0.1%
N12Primary sensorimotor network (mouth) (9)Cerebellum (12)−0.0%
N13Left fronto-parietal network (1)Basal ganglia (8)−0.2%
  1. Rank indicates relevance with working memory function according to the BrainMap ICA.

Table 2
Demographic data of the multiple psychiatric diagnoses dataset.
https://doi.org/10.7554/eLife.38844.011
DiagnosisSiteMeasurePatientsControlsTestP-value
SCZKYUN5860--
 Age37.935.2t116 = 1.70.1
(9.3)(8.4)
 Male %52%67%Fisher’s exact test0.13
MDDHRUN7763--
 Age41.639.3t138 = 1.30.21
(11.2)(12.0)
 Male %56%46%Fisher’s exact test0.31
OCDKPMN4647--
 Age32.230.3t91 = 1.10.28
(9.9)(8.7)
 Male %37%45%Fisher’s exact test0.53
ASDUTK (site1)N3333--
 Age32.834.7t64 = −1.00.3
(8.4)(7.0)
 Male %64%55%Fisher’s exact test0.62
 SHU (site2)N3638--
 Age29.932.5t72 = −1.50.14
(7.2)(7.4)
 Male %100%100%Fisher’s exact test1
 PooledN6971--
 Age31.333.5t138 = −1.80.08
(7.9)(7.2)
 Male %83%79%Fisher’s exact test0.67
  1. Site: KYU, Kyoto University; HRU, Hiroshima University; KPM, Kyoto Prefectural University of Medicine; UTK, University of Tokyo; SHU, Showa University. Measure:N’ indicates the number of subjects; ‘Age’ is shown as mean (SD); ‘Male %” is the fraction of male. The tests and p-values compare the patient and control groups within-site.

Additional files

Supplementary file 1

Resting state fMRI scan parameters and their values.

https://doi.org/10.7554/eLife.38844.017
Supplementary file 2

Whole-brain intrinsic functional network labels and their component regions.

https://doi.org/10.7554/eLife.38844.018
Transparent reporting form
https://doi.org/10.7554/eLife.38844.019

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Masahiro Yamashita
  2. Yujiro Yoshihara
  3. Ryuichiro Hashimoto
  4. Noriaki Yahata
  5. Naho Ichikawa
  6. Yuki Sakai
  7. Takashi Yamada
  8. Noriko Matsukawa
  9. Go Okada
  10. Saori C Tanaka
  11. Kiyoto Kasai
  12. Nobumasa Kato
  13. Yasumasa Okamoto
  14. Ben Seymour
  15. Hidehiko Takahashi
  16. Mitsuo Kawato
  17. Hiroshi Imamizu
(2018)
A prediction model of working memory across health and psychiatric disease using whole-brain functional connectivity
eLife 7:e38844.
https://doi.org/10.7554/eLife.38844