1. Kantor, E. D., Rehm, C. D., Haas, J. S., Chan, A. T. & Giovannucci, E. L. Trends in Prescription Drug Use Among Adults in the
United States From 1999–2012. JAMA 314, 1818–1831 (2015).
2. Abbing-Karahagopian, V. et al. Antidepressant prescribing in five European countries: application of common definitions to assess
the prevalence, clinical observations, and methodological implications. Eur. J. Clin. Pharmacol. 70, 849–857 (2014).
3. Kjosavik, S. R., Hunskaar, S., Aarsland, D. & Ruths, S. Initial prescription of antipsychotics and antidepressants in general practice
and specialist care in Norway. Acta Psychiatr. Scand. 123, 459–465 (2011).
4. Uchida, N. et al. International study on antidepressant prescription pattern at 20 teaching hospitals and major psychiatric
institutions in East Asia: Analysis of 1898 cases from China, Japan, Korea, Singapore and Taiwan. Psychiatry Clin. Neurosci. 61,
522–528 (2007).
5. Trivedi, M. H. et al. Evaluation of outcomes with citalopram for depression using measurement-based care in STAR*D: implications
for clinical practice. Am. J. Psychiatry 163, 28–40 (2006).
6. Vai, B. et al. Fronto-limbic effective connectivity as possible predictor of antidepressant response to SSRI administration. Eur.
Neuropsychopharmacol. 26, 2000–2010 (2016).
7. Cao, H. et al. Test-retest reliability of fMRI-based graph theoretical properties during working memory, emotion processing, and
resting state. Neuroimage 84, 888–900 (2014).
8. Drysdale, A. T. et al. Resting-state connectivity biomarkers define neurophysiological subtypes of depression. Nat. Med. 23, 28–38
(2017).
9. Sundermann, B. et al. Diagnostic classification of unipolar depression based on resting-state functional connectivity MRI: effects of
generalization to a diverse sample. J Neural Transm, https://doi.org/10.1007/s00702-016-1673-8 (2016).
10. Zeng, L. L. et al. Identifying major depression using whole-brain functional connectivity: a multivariate pattern analysis. Brain 135,
1498–1507 (2012).
Scientific Reports |
(2020) 10:3542 | https://doi.org/10.1038/s41598-020-60527-z
10
www.nature.com/scientificreports/
www.nature.com/scientificreports
11. Kupfer, D. J., Frank, E. & Phillips, M. L. Major depressive disorder: new clinical, neurobiological, and treatment perspectives. Lancet
379, 1045–1055 (2012).
12. Singh, M. K. & Gotlib, I. H. The neuroscience of depression: implications for assessment and intervention. Behav. Res. Ther. 62,
60–73 (2014).
13. Dichter, G. S., Gibbs, D. & Smoski, M. J. A systematic review of relations between resting-state functional-MRI and treatment
response in major depressive disorder. J. Affect. Disord. 172, 8–17 (2015).
14. Hamilton, J. P., Chen, M. C. & Gotlib, I. H. Neural systems approaches to understanding major depressive disorder: an intrinsic
functional organization perspective. Neurobiol. Dis. 52, 4–11 (2013).
15. Arbabshirani, M. R., Plis, S., Sui, J. & Calhoun, V. D. Single subject prediction of brain disorders in neuroimaging: Promises and
pitfalls. Neuroimage 145, 137–165 (2017).
16. Kendler, K. S. The diagnostic validity of melancholic major depression in a population-based sample of female twins. Arch. Gen.
Psychiatry 54, 299–304 (1997).
17. Sun, N. et al. A comparison of melancholic and nonmelancholic recurrent major depression in Han Chinese women. Depress.
Anxiety 29, 4–9 (2012).
18. CONVERGE consortium. Sparse whole-genome sequencing identifies two loci for major depressive disorder. Nat. 523, 588–591
(2015).
19. Hyett, M. P., Breakspear, M. J., Friston, K. J., Guo, C. C. & Parker, G. B. Disrupted effective connectivity of cortical systems
supporting attention and interoception in melancholia. JAMA Psychiatry 72, 350–358 (2015).
20. Parker, G., Hadzi-Pavlovic, D. Melancholia: A Disorder of Movement and Mood: A Phenomenological and Neurobiological Review.
Cambridge and New York: Cambridge University Press (1996).
21. Yan, C. G. et al. Reduced default mode network functional connectivity in patients with recurrent major depressive disorder. Proc.
Natl Acad. Sci. USA 116, 9078–9083 (2019).
22. Di, X. & Biswal, B. B. Dynamic brain functional connectivity modulated by resting-state networks. Brain Struct. Funct. 220, 37–46
(2015).
23. Borsboom, D., Cramer, A. O. J. & Kalis, A. Brain disorders? Not really: Why network structures block reductionism in
psychopathology research. Behav. Brain Sci. 42(e2), 1–63 (2019).
24. Seeley, W. W. et al. Dissociable intrinsic connectivity networks for salience processing and executive control. J. Neurosci. 27,
2349–2356 (2007).
25. Woo, C. W., Chang, L. J., Lindquist, M. A. & Wager, T. D. Building better biomarkers: brain models in translational neuroimaging.
Nat. Neurosci. 20, 365–377 (2017).
26. Sheehan, D. V. et al. The Mini-International Neuropsychiatric Interview (M.I.N.I.): the development and validation of a structured
diagnostic psychiatric interview for DSM-IV and ICD-10. J Clin Psychiatry 59 (Suppl 20), 22–33; quiz 34–57. Review (1998).
27. Otsubo, T. et al. Reliability and validity of Japanese version of the Mini-International Neuropsychiatric Interview. Psychiatry Clin.
Neurosci. 59, 517–526 (2005).
28. Beck, A. T., Steer, R. A., Brown, G. K. Manual for the Beck Depression Inventory-II. San Antonio, TX: Psychological Corporation
(1996).
29. Matsuoka, K., Uno, M., Kasai, K., Koyama, K. & Kim, Y. Estimation of premorbid IQ in individuals with Alzheimer’s disease using
Japanese ideographic script (Kanji) compound words: Japanese version of National Adult Reading Test. Psychiatry Clin. Neurosci.
60, 332–339 (2006).
30. Hamilton, M. A rating scale for depression. J. Neurology, Neurosurg. Psychiatry 23, 56–62 (1960).
31. Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L. & Petersen, S. E. Spurious but systematic correlations in functional
connectivity MRI networks arise from subject motion. Neuroimage 59, 2142–2154 (2012).
32. Perrot, M., Rivière, D. & Mangin, J. F. Cortical sulci recognition and spatial normalization. Med. Image Anal. 15, 529–550 (2011).
33. Rivière, D. et al. Automatic recognition of cortical sulci of the human brain using a congregation of neural networks. Med. Image
Anal. 6, 77–92 (2002).
34. Witten, D. M., Tibshirani, R. & Hastie, T. A penalized matrix decomposition, with applications to sparse principal components and
canonical correlation analysis. Biostatistics 10, 515–534 (2009).
35. Yamashita, O., Sato, M. A., Yoshioka, T., Tong, F. & Kamitani, Y. Sparse estimation automatically selects voxels relevant for the
decoding of fMRI activity patterns. Neuroimage 42, 1414–1429 (2008).
36. Tipping, M. E., Faul, A. Fast marginal likelihood maximisation for sparse Bayesian models. Proceedings of the Ninth International
Workshop on Artificial Intelligence and Statistics (2003).
37. Yahata, N. et al. A small number of abnormal brain connections predicts adult autism spectrum disorder. Nat. Commun. 7, 11254
(2016).
38. McIntyre, R., Kennedy, S., Bagby, R. M. & Bakish, D. Assessing full remission. J. Psychiatry Neurosci. 27, 235–239 (2002).
39. Yoshihara, Y. et al. Overlapping but asymmetrical relationships between schizophrenia and autism revealed by brain connectivity.
bioRxiv 403212; https://doi.org/10.1101/403212 (2018).
40. Greening, S. G., Finger, E. C. & Mitchell, D. G. Parsing decision making processes in prefrontal cortex: response inhibition,
overcoming learned avoidance, and reversal learning. Neuroimage 54, 1432–1441 (2011).
41. Dombrovski, A. Y. et al. Corticostriatothalamic reward prediction error signals and executive control in late-life depression. Psychol.
Med. 45, 1413–1424 (2015).
42. Mascaro, J. S., Rilling, J. K., Negi. L. T., Raison, C. L. Compassion meditation enhances empathic accuracy and related neural activity.
Soc Cogn Affect Neurosci. 48–55 (2013).
43. LeWinn, K. Z. et al. An exploratory examination of reappraisal success in depressed adolescents: Preliminary evidence of functional
differences in cognitive control brain regions. J. Affect. Disord. 240, 155–164 (2018).
44. Ferrari, C. et al. The Dorsomedial Prefrontal Cortex Plays a Causal Role in Integrating Social Impressions from Faces and Verbal
Descriptions. Cereb. Cortex. 26, 156–165 (2016).
45. Mattavelli, G., Cattaneo, Z. & Papagno, C. Transcranial magnetic stimulation of medial prefrontal cortex modulates face expressions
processing in a priming task. Neuropsychologia. 49, 992–998 (2011).
46. Jonides, J. & Nee, D. E. Brain mechanisms of proactive interference in working memory. Neurosci. 139, 181–193 (2006).
47. Engels, A. S. et al. Co-occurring anxiety influences patterns of brain activity in depression. Cogn. Affect. Behav. Neurosci. 10, 141–156
(2010).
48. Andrews, P. W., Bharwani, A., Lee, K. R., Fox, M. & Thomson, J. A. Jr. Is serotonin an upper or a downer? The evolution of the
serotonergic system and its role in depression and the antidepressant response. Neurosci. Biobehav. Rev. 51, 164–188 (2015).
49. Baldinger, P. et al. Regional differences in SERT occupancy after acute and prolonged SSRI intake investigated by brain PET.
Neuroimage 88, 252–262 (2014).
50. Ardila, A., Bernal, B. & Rosselli, M. How Localized are Language Brain Areas? A Review of Brodmann Areas Involvement in Oral
Language. Arch. Clin. Neuropsychol. 31, 112–122 (2016).
51. Pascual-Leone, A., Rubio, B., Pallardó, F. & Catalá, M. D. Rapid-rate transcranial magnetic stimulation of left dorsolateral prefrontal
cortex in drug-resistant depression. Lancet 348, 233–237 (1996).
52. Gotlib, I. H. Frontal EEG Alpha Asymmetry, Depression, and Cognitive Functioning. Cognition Emot. 12, 449–478 (1998).
Scientific Reports |
(2020) 10:3542 | https://doi.org/10.1038/s41598-020-60527-z
11
www.nature.com/scientificreports/
www.nature.com/scientificreports
53. Brzezicka, A., Kamiński, J., Kamińska, O. K., Wołyńczyk-Gmaj, D. & Sedek, G. Frontal EEG alpha band asymmetry as a predictor of
reasoning deficiency in depressed people. Cogn. Emot. 31, 868–878 (2017).
54. Mennella, R., Patron, E. & Palomba, D. Frontal alpha asymmetry neurofeedback for the reduction of negative affect and anxiety.
Behav. Res. Ther. 92, 32–40 (2017).
55. Zotev, V. et al. Correlation between amygdala BOLD activity and frontal EEG asymmetry during real-time fMRI neurofeedback
training in patients with depression. Neuroimage Clin. 11, 224–238 (2016).
56. Okada, G., Okamoto, Y., Morinobu, S., Yamawaki, S. & Yokota, N. Attenuated left prefrontal activation during a verbal fluency task
in patients with depression. Neuropsychobiology 47, 21–26 (2003).
57. Takamura, M. et al. Disrupted Brain Activation and Deactivation Pattern during Semantic Verbal Fluency Task in Patients with
Major Depression. Neuropsychobiology 74, 69–77 (2016).
58. Chen, A. C. et al. Causal interactions between fronto-parietal central-executive and default-mode networks in humans. Proc. Natl
Acad. Sci. USA 110, 19944–19949 (2013).
59. Cho, S. S. & Strafella, A. P. rTMS of the Left Dorsolateral Prefrontal Cortex Modulates Dopamine Release in the Ipsilateral Anterior
Cingulate Cortex and Orbitofrontal Cortex. PLoS ONE 4, e6725, https://doi.org/10.1371/journal.pone.0006725 (2009).
60. Fox, M. D., Buckner, R. L., White, M. P., Greicius, M. D. & Pascual-Leone, A. Efficacy of transcranial magnetic stimulation targets
for depression is related to intrinsic functional connectivity with the subgenual cingulate. Biol. Psychiatry 72, 595–603 (2012).
61. Bassett, D. S. & Khambhati, A. N. A network engineering perspective on probing and perturbing cognition with neurofeedback.
Ann. N. Y. Acad. Sci. 1396, 126–143 (2017).
62. Orndorff-Plunkett, F., Singh, F., Aragón, O. R. & Pineda, J. A. Assessing the Effectiveness of Neurofeedback Training in the Context
of Clinical and Social Neuroscience. Brain Sci. 7, 95, https://doi.org/10.3390/brainsci7080095 (2017).
63. Yamada, T. et al. Resting-State Functional Connectivity-Based Biomarkers and Functional MRI-Based Neurofeedback for
Psychiatric Disorders: A Challenge for Developing Theranostic Biomarkers. Int. J. Neuropsychopharmacol. 20, 769–781 (2017).
64. Dinga, R. et al. Evaluating the evidence for biotypes of depression: attempted replication of Drysdale 2017. bioRxiv, https://doi.
org/10.1101/416321 (2018).
65. Oathes, D. J., Patenaude, B., Schatzberg, A. F. & Etkin, A. Neurobiological signatures of anxiety and depression in resting-state
functional magnetic resonance imaging. Biol. Psychiatry 77, 385–393 (2015).
Acknowledgements
This research was conducted as a contract project supported by AMED under Grant Number JP18dm0307008,
and under the “Application of DecNef for development of diagnostic and cure system for mental disorders and
construction of clinical application bases” of the Strategic Research Program for Brain Sciences from Japan
Agency for Medical Research and development, AMED. N.I., G.O., M.T., S.Y., Y.O. are also partially supported
by “Integrated research on neuropsychiatric disorders”, “the Integrated Research on Depression, Dementia and
Development Disorders”, and Brain/MINDS Beyond (19dm0307002 h0002) by AMED. B.S. is partially supported
by Versus Arthritis (21537).
Author contributions
Y.O. had full access to all of the data in the study and takes responsibility for the integrity of the data and the
accuracy of the data analysis. N.I. and G.L. contributed equally as first authors. G.O., M.T, N.I., N.Y., R.H., T.Y.,
M.Y., T.S., S.M., M.M., Y.Y., H.T., K.K., N.K., S.Y. and Y.O. involved in collecting the data. N.I. and Y.O. discussed
the design of clinical data and ran the analyses. G.L. and N.Y. designed the preprocessing pipeline and prepared
the scripts for the Biomarker Toolbox. N.I. and G.L. created tables and figures. G.L., J.M., N.Y. and M.K. consulted
on all statistical analyses. N.I., G.L., N.Y., B. S., M.K. and Y.O. wrote the manuscript. All authors discussed the
results and conclusions for editing the manuscript.
Competing interests
The authors declare no competing interests.
Additional information
Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-60527-z.
Correspondence and requests for materials should be addressed to B.S. or Y.O.
Reprints and permissions information is available at www.nature.com/reprints.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and
institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International
License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
© The Author(s) 2020
Scientific Reports |
(2020) 10:3542 | https://doi.org/10.1038/s41598-020-60527-z
12
Supplementary information for “Primary functional brain connections
associated with melancholic major depressive disorder and modulation by
antidepressants”
Naho Ichikawa1, Giuseppe Lisi2, Noriaki Yahata3, Go Okada1, Masahiro Takamura1, Ryu-ichiro
Hashimoto4, Takashi Yamada2, Makiko Yamada3,5, Tetsuya Suhara3, Sho Moriguchi6, Masaru
Mimura6, Yujiro Yoshihara7, Hidehiko Takahashi7,8, Kiyoto Kasai9, Nobumasa Kato4, Shigeto
Yamawaki1, Ben Seymour2,10*, Mitsuo Kawato2, Jun Morimoto2, & Yasumasa Okamoto1*
*oy@hiroshima-u.ac.jp
Department of Psychiatry and Neurosciences, Graduate School of Biomedical Sciences,
Hiroshima University, Hiroshima, JAPAN; 2ATR Brain Information Communication Research
Laboratory Group, Kyoto, JAPAN; 3Institute for Quantum Life Science, National Institute of
Radiological Sciences, National Institutes for Quantum and Radiological Science and Technology,
Chiba, JAPAN; 4Medical Institute of Developmental Disabilities Research, Showa University,
Tokyo, JAPAN; 5Department of Functional Brain Imaging Research, National Institutes for
Quantum and Radiological Science and Technology, Chiba, JAPAN;
Department of
Neuropsychiatry, Keio University School of Medicine, Tokyo, JAPAN; 7Department of
Psychiatry, Kyoto University Graduate School of Medicine, Kyoto, JAPAN; 8Graduate School of
Medical and Dental Sciences Tokyo Medical and Dental University, Tokyo, JAPAN 9Department
of Youth Mental Health, Graduate School of Medicine, The University of Tokyo, Tokyo, JAPAN;
10
Computational and Biological Learning Lab, Cambridge University, Cambridge, UK.
Supplementary Information
Supplementary Methods.
Feature selection with reduction of nuisance variable effects
Application to non-melancholic MDD, treatment-resistant MDD.
Application to other psychiatric disorders.
Application to euthymic MDD.
Changes of each FC: post minus pre antidepressant treatments.
Supplementary Results.
Tables S1: Demographic information of the participants
Tables S2: Scanner information and resting-state fMRI protocols
Tables S3: Head motion during resting-state
Figure S1: Generalization to an independent cohort assessed by permutation tests
Supplementary Methods.
Feature selection with reduction of nuisance variable effects
As specific steps, in the first data matrix X1 of demographic information, the column of diagnosis
contained either 1 or 0 (i.e., MDD or HC). The next four columns of scanner/site contained either
[1 0 0 0] for the Site1, [0 1 0 0] for the Site2, [0 0 1 0] for the Site3, or [0 0 0 1] for the Site4. The
sixth column was for age, and the seventh column was for sex (i.e., 1 for male and 0 for female).
The second data matrix X2 pools the off-diagonal lower-triangular elements of correlation matrix
as the FC of a single subject. By applying L1-SCCA to these data, the sparse projection matrices
of V1 and V2 were derived from the pair of data matrices (i.e., X1 and X2). We here defined the
diagnostic canonical variables and the diagnostic canonical constraint, and which should be
associated only with the diagnostic label. We used the columns of V1 that have non-zero elements
only in the row corresponding to the diagnostic label, and only the columns of V2 corresponding
to the diagnostic canonical variables based on the sum of absolute value across columns as union
of features across diagnostic canonical variables. This means, the columns which had non-zero
values in the rows corresponding to the scanner/site, age, sex were disregarded. In this way,
identifying the FCs corresponding to the diagnostic canonical variables enabled us to select only
essential FCs for MDD/HC classification, and simultaneously reduced undesirable effects of
nuisance variables (i.e., scanner/site, age, sex).
Application to non-melancholic MDD, treatment-resistant MDD.
In order to test if the classifier was specific to the characteristics of melancholic MDD, we
applied the same classifier to the datasets of non-melancholic and treatment-resistant MDD
(Demographics are shown in Supplementary Table S1b, with two MDD patients lacking MINI
scores of previous episodes, melancholia, and comorbidity). Non-melancholic MDDs are from
the all MDD dataset with more than mild depressive state (BDI score > = 17).
Application to other psychiatric disorders.
Autism spectrum disorder (ASD) dataset was adopted from our previous investigation. In
order to minimize any effect from comorbidity of depression, the autism spectrum disorder dataset
was limited to the data with no active antidepressant medication. A total of 110 participants
included 74 ASD patients (No. of male/female = 58/16, Age (year): mean (SD) = 31.5 (8.5)) and
36 healthy controls (No. of male/female = 30/6, Age (year): mean (SD) = 30.9 (6.9)).
For the schizophrenia spectrum disorder (SSD) data (Yoshihara et al., in submission), all
participants provided written informed consent that was approved by the Committee on Medical
Ethics of Kyoto University. A total of 170 participants included 68 SSD patients (No. of
male/female = 33/35, Age (year): mean (SD) = 38.4 (9.1), 64 Schizophrenia and 4 Schizoaffective
disorder, Duration of illness (year): mean (SD) = 12.8 (7.8)) and 102 healthy controls (No. of
male/female = 62/40, Age (year): mean (SD) = 31.1 (9.6)). See the Supplementary Table S1 for
details of MRI experimental settings for each site.
Application to euthymic MDD.
Thirty-six patients who are in euthymic states for more than two months were recruited, and
their demographic information was shown in Supplementary Table S1c. Their resting state
fMRI and T1 structural data were scanned using the identical scan protocols at the Site 4, and
went through the identical preprocessing steps.
Changes of each FC: post minus pre antidepressant treatments.
For the analysis of pre-post treatment effects, we computed contribution scores of each FC
using the following equations:
Equation (1)
where N is the number of patients who underwent the treatment and
is the classifier weight.
In order to have a point of reference, we computed the same contribution score, for the MDD and
healthy control populations:
Equation (2)
where M is the total number of MDD patients, and N is the number of healthy control subjects.
The
represents the difference in average
between healthy controls and MDD, weighted
by the classifier’s weight. It should be noted that the
the positive class in the classifier, and
subtraction) in Formula 2. The
is the subtrahend (i.e. negative sign in the
represents the difference in average
pre- antidepressant treatments. Negative value of
the post-treatment
treatment
should always be negative since MDD is
between post- and
(i.e. post- minus pre- treatment) means that
is closer to healthy controls, while positive value means that post-
is closer to MDD. For each
, we compared
and
using the Welch’s t-
test, due to the different sample sizes between HCs, MDDs, and the MDD patients who underwent
treatment:
Equation (3)
where
is the sample variance of the
the sum of the variances of
the sum of variances of
, and it is obtained using the variance sum law as
and
. Similarly,
and
is computed as
. The sample sizes are computes as
and
. The P value of each test was corrected for multiple
comparisons by the Benjamini–Hochberg procedure.
Supplementary Results.
Table S1: Demographic information of the melancholic MDD classifier.
Table S2: Scanner information and resting-state fMRI protocols of melancholic MDD, healthy
control, other subtypes of MDD, autism spectrum disorder, and schizophrenia. All the euthymic
MDD data shown in Table S1c was collected at the Site 4 with identical settings.
Table S3: Head motion of MDD and healthy control in the training dataset and test dataset.
10
Figure S1. Classification results of the melancholic MDD biomarker and its generalization
performance to other subtypes and psychiatric disorders. Density of the weighted linear sum
computed based on the melancholic MDD biomarker for a.) training dataset (melancholic MDD
n = 65, healthy control: n = 65) and b.) test dataset from an independent site (melancholic MDD:
n = 11, healthy control: n = 40). Permutation tests show the histogram of the permutation test
(1,000 repetitions) for c) the training dataset LOOCV, and show d.) a completely independent test
dataset accuracies, and the binomial distribution is shown as a green curve. The accuracies of the
melancholic MDD classifier trained and tested without permutation were shown as red vertical
lines. The results of permutation test were significant for LOOCV (p = .002) and for the
independent test dataset (p = .040). * p < .05, *** p < .005.
11
...