AustralieFrV1, Pmc, Corpus, bibRecord, 000E20

The dynamics of resting fluctuations in the brain: metastability and its dynamical cortical core

Identifieur interne : 000E20 ( Pmc/Corpus ); précédent : 000E19; suivant : 000E21

The dynamics of resting fluctuations in the brain: metastability and its dynamical cortical core

Auteurs : Gustavo Deco ; Morten L. Kringelbach ; Viktor K. Jirsa ; Petra Ritter

Source :

Scientific Reports [ 2045-2322 ] ; 2017.

RBID : PMC:5465179

Abstract

In the human brain, spontaneous activity during resting state consists of rapid transitions between functional network states over time but the underlying mechanisms are not understood. We use connectome based computational brain network modeling to reveal fundamental principles of how the human brain generates large-scale activity observable by noninvasive neuroimaging. We used structural and functional neuroimaging data to construct whole- brain models. With this novel approach, we reveal that the human brain during resting state operates at maximum metastability, i.e. in a state of maximum network switching. In addition, we investigate cortical heterogeneity across areas. Optimization of the spectral characteristics of each local brain region revealed the dynamical cortical core of the human brain, which is driving the activity of the rest of the whole brain. Brain network modelling goes beyond correlational neuroimaging analysis and reveals non-trivial network mechanisms underlying non-invasive observations. Our novel findings significantly pertain to the important role of computational connectomics in understanding principles of brain function.

Url:

http://www.ncbi.nlm.nih.gov/pmc/articles/PMC5465179

DOI: 10.1038/s41598-017-03073-5
PubMed: 28596608
PubMed Central: 5465179

Links to Exploration step

PMC:5465179

Le document en format XML

<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">The dynamics of resting fluctuations in the brain: metastability and its dynamical cortical core</title>
<author><name sortKey="Deco, Gustavo" sort="Deco, Gustavo" uniqKey="Deco G" first="Gustavo" last="Deco">Gustavo Deco</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2172 2676</institution-id>
<institution-id institution-id-type="GRID">grid.5612.0</institution-id>
<institution>Center for Brain and Cognition, Computational Neuroscience Group, Department of Information and Communication Technologies,</institution>
<institution>Universitat Pompeu Fabra,</institution>
</institution-wrap>
Roc Boronat 138, Barcelona, 08018 Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff2">Institució Catalana de la Recerca i Estudis Avançats (ICREA), Passeig Lluís Companys 23, Barcelona, 08010 Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0041 5028</institution-id>
<institution-id institution-id-type="GRID">grid.419524.f</institution-id>
<institution>Department of Neuropsychology,</institution>
<institution>Max Planck Institute for Human Cognitive and Brain Sciences,</institution>
</institution-wrap>
04103 Leipzig, Germany</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff4"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 7857</institution-id>
<institution-id institution-id-type="GRID">grid.1002.3</institution-id>
<institution>School of Psychological Sciences,</institution>
<institution>Monash University,</institution>
</institution-wrap>
Melbourne, Clayton VIC 3800, Australia</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kringelbach, Morten L" sort="Kringelbach, Morten L" uniqKey="Kringelbach M" first="Morten L." last="Kringelbach">Morten L. Kringelbach</name>
<affiliation><nlm:aff id="Aff5"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 8948</institution-id>
<institution-id institution-id-type="GRID">grid.4991.5</institution-id>
<institution>Department of Psychiatry,</institution>
<institution>University of Oxford,</institution>
</institution-wrap>
Oxford, UK</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff6"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1956 2722</institution-id>
<institution-id institution-id-type="GRID">grid.7048.b</institution-id>
<institution>Center of Music in the Brain (MIB), Clinical Medicine,</institution>
<institution>Aarhus University,</institution>
</institution-wrap>
Aarhus C, DK Denmark</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Jirsa, Viktor K" sort="Jirsa, Viktor K" uniqKey="Jirsa V" first="Viktor K." last="Jirsa">Viktor K. Jirsa</name>
<affiliation><nlm:aff id="Aff7"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0541 5643</institution-id>
<institution-id institution-id-type="GRID">grid.462494.9</institution-id>
<institution>Aix Marseille Univ, INSERM,</institution>
<institution>INS, Inst Neurosci Syst,</institution>
</institution-wrap>
Marseille, France</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ritter, Petra" sort="Ritter, Petra" uniqKey="Ritter P" first="Petra" last="Ritter">Petra Ritter</name>
<affiliation><nlm:aff id="Aff8"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2105 1091</institution-id>
<institution-id institution-id-type="GRID">grid.4372.2</institution-id>
<institution></institution>
<institution>Max-Planck Institute for Cognitive and Brain Sciences,</institution>
</institution-wrap>
Leipzig, Germany</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff9">Department of Neurology, Charité, Charitéplatz 1, 10117 Berlin, Germany</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">28596608</idno>
<idno type="pmc">5465179</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC5465179</idno>
<idno type="RBID">PMC:5465179</idno>
<idno type="doi">10.1038/s41598-017-03073-5</idno>
<date when="2017">2017</date>
<idno type="wicri:Area/Pmc/Corpus">000E20</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000E20</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">The dynamics of resting fluctuations in the brain: metastability and its dynamical cortical core</title>
<author><name sortKey="Deco, Gustavo" sort="Deco, Gustavo" uniqKey="Deco G" first="Gustavo" last="Deco">Gustavo Deco</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2172 2676</institution-id>
<institution-id institution-id-type="GRID">grid.5612.0</institution-id>
<institution>Center for Brain and Cognition, Computational Neuroscience Group, Department of Information and Communication Technologies,</institution>
<institution>Universitat Pompeu Fabra,</institution>
</institution-wrap>
Roc Boronat 138, Barcelona, 08018 Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff2">Institució Catalana de la Recerca i Estudis Avançats (ICREA), Passeig Lluís Companys 23, Barcelona, 08010 Spain</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0041 5028</institution-id>
<institution-id institution-id-type="GRID">grid.419524.f</institution-id>
<institution>Department of Neuropsychology,</institution>
<institution>Max Planck Institute for Human Cognitive and Brain Sciences,</institution>
</institution-wrap>
04103 Leipzig, Germany</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff4"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 7857</institution-id>
<institution-id institution-id-type="GRID">grid.1002.3</institution-id>
<institution>School of Psychological Sciences,</institution>
<institution>Monash University,</institution>
</institution-wrap>
Melbourne, Clayton VIC 3800, Australia</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kringelbach, Morten L" sort="Kringelbach, Morten L" uniqKey="Kringelbach M" first="Morten L." last="Kringelbach">Morten L. Kringelbach</name>
<affiliation><nlm:aff id="Aff5"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 8948</institution-id>
<institution-id institution-id-type="GRID">grid.4991.5</institution-id>
<institution>Department of Psychiatry,</institution>
<institution>University of Oxford,</institution>
</institution-wrap>
Oxford, UK</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff6"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1956 2722</institution-id>
<institution-id institution-id-type="GRID">grid.7048.b</institution-id>
<institution>Center of Music in the Brain (MIB), Clinical Medicine,</institution>
<institution>Aarhus University,</institution>
</institution-wrap>
Aarhus C, DK Denmark</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Jirsa, Viktor K" sort="Jirsa, Viktor K" uniqKey="Jirsa V" first="Viktor K." last="Jirsa">Viktor K. Jirsa</name>
<affiliation><nlm:aff id="Aff7"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0541 5643</institution-id>
<institution-id institution-id-type="GRID">grid.462494.9</institution-id>
<institution>Aix Marseille Univ, INSERM,</institution>
<institution>INS, Inst Neurosci Syst,</institution>
</institution-wrap>
Marseille, France</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Ritter, Petra" sort="Ritter, Petra" uniqKey="Ritter P" first="Petra" last="Ritter">Petra Ritter</name>
<affiliation><nlm:aff id="Aff8"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2105 1091</institution-id>
<institution-id institution-id-type="GRID">grid.4372.2</institution-id>
<institution></institution>
<institution>Max-Planck Institute for Cognitive and Brain Sciences,</institution>
</institution-wrap>
Leipzig, Germany</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff9">Department of Neurology, Charité, Charitéplatz 1, 10117 Berlin, Germany</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Scientific Reports</title>
<idno type="eISSN">2045-2322</idno>
<imprint><date when="2017">2017</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>In the human brain, spontaneous activity during resting state consists of rapid transitions between functional network states over time but the underlying mechanisms are not understood. We use connectome based computational brain network modeling to reveal fundamental principles of how the human brain generates large-scale activity observable by noninvasive neuroimaging. We used structural and functional neuroimaging data to construct whole- brain models. With this novel approach, we reveal that the human brain during resting state operates at maximum metastability, i.e. in a state of maximum network switching. In addition, we investigate cortical heterogeneity across areas. Optimization of the spectral characteristics of each local brain region revealed the dynamical cortical core of the human brain, which is driving the activity of the rest of the whole brain. Brain network modelling goes beyond correlational neuroimaging analysis and reveals non-trivial network mechanisms underlying non-invasive observations. Our novel findings significantly pertain to the important role of computational connectomics in understanding principles of brain function.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Deco, G" uniqKey="Deco G">G Deco</name>
</author>
<author><name sortKey="Tononi, G" uniqKey="Tononi G">G Tononi</name>
</author>
<author><name sortKey="Boly, M" uniqKey="Boly M">M Boly</name>
</author>
<author><name sortKey="Kringelbach, Ml" uniqKey="Kringelbach M">ML Kringelbach</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Berridge, Kc" uniqKey="Berridge K">KC Berridge</name>
</author>
<author><name sortKey="Kringelbach, Ml" uniqKey="Kringelbach M">ML Kringelbach</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Friston, K" uniqKey="Friston K">K Friston</name>
</author>
<author><name sortKey="Kiebel, S" uniqKey="Kiebel S">S Kiebel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kringelbach, Ml" uniqKey="Kringelbach M">ML Kringelbach</name>
</author>
<author><name sortKey="Mcintosh, Ar" uniqKey="Mcintosh A">AR McIntosh</name>
</author>
<author><name sortKey="Ritter, P" uniqKey="Ritter P">P Ritter</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
<author><name sortKey="Deco, G" uniqKey="Deco G">G Deco</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tognoli, E" uniqKey="Tognoli E">E Tognoli</name>
</author>
<author><name sortKey="Kelso, Ja" uniqKey="Kelso J">JA Kelso</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Friston, Kj" uniqKey="Friston K">KJ Friston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chaudhuri, R" uniqKey="Chaudhuri R">R Chaudhuri</name>
</author>
<author><name sortKey="Knoblauch, K" uniqKey="Knoblauch K">K Knoblauch</name>
</author>
<author><name sortKey="Gariel, Ma" uniqKey="Gariel M">MA Gariel</name>
</author>
<author><name sortKey="Kennedy, H" uniqKey="Kennedy H">H Kennedy</name>
</author>
<author><name sortKey="Wang, Xj" uniqKey="Wang X">XJ Wang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Simony, E" uniqKey="Simony E">E Simony</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Tversky, A" uniqKey="Tversky A">A Tversky</name>
</author>
<author><name sortKey="Kahneman, D" uniqKey="Kahneman D">D Kahneman</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Ghosh, A" uniqKey="Ghosh A">A Ghosh</name>
</author>
<author><name sortKey="Rho, Y" uniqKey="Rho Y">Y Rho</name>
</author>
<author><name sortKey="Mcintosh, Ar" uniqKey="Mcintosh A">AR McIntosh</name>
</author>
<author><name sortKey="Kotter, R" uniqKey="Kotter R">R Kotter</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Deco, G" uniqKey="Deco G">G Deco</name>
</author>
<author><name sortKey="Jirsa, V" uniqKey="Jirsa V">V Jirsa</name>
</author>
<author><name sortKey="Mcintosh, Ar" uniqKey="Mcintosh A">AR McIntosh</name>
</author>
<author><name sortKey="Sporns, O" uniqKey="Sporns O">O Sporns</name>
</author>
<author><name sortKey="Kotter, R" uniqKey="Kotter R">R Kotter</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cabral, J" uniqKey="Cabral J">J Cabral</name>
</author>
<author><name sortKey="Kringelbach, Ml" uniqKey="Kringelbach M">ML Kringelbach</name>
</author>
<author><name sortKey="Deco, G" uniqKey="Deco G">G Deco</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Schirner, M" uniqKey="Schirner M">M Schirner</name>
</author>
<author><name sortKey="Rothmeier, S" uniqKey="Rothmeier S">S Rothmeier</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
<author><name sortKey="Mcintosh, Ar" uniqKey="Mcintosh A">AR McIntosh</name>
</author>
<author><name sortKey="Ritter, P" uniqKey="Ritter P">P Ritter</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Deco, G" uniqKey="Deco G">G Deco</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
<author><name sortKey="Mcintosh, Ar" uniqKey="Mcintosh A">AR McIntosh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Deco, G" uniqKey="Deco G">G Deco</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
<author><name sortKey="Mcintosh, Ar" uniqKey="Mcintosh A">AR McIntosh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Damoiseaux, Js" uniqKey="Damoiseaux J">JS Damoiseaux</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cabral, J" uniqKey="Cabral J">J Cabral</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Breakspear, M" uniqKey="Breakspear M">M Breakspear</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wildie, M" uniqKey="Wildie M">M Wildie</name>
</author>
<author><name sortKey="Shanahan, M" uniqKey="Shanahan M">M Shanahan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cabral, J" uniqKey="Cabral J">J Cabral</name>
</author>
<author><name sortKey="Hugues, E" uniqKey="Hugues E">E Hugues</name>
</author>
<author><name sortKey="Sporns, O" uniqKey="Sporns O">O Sporns</name>
</author>
<author><name sortKey="Deco, G" uniqKey="Deco G">G Deco</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Freyer, F" uniqKey="Freyer F">F Freyer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Freyer, F" uniqKey="Freyer F">F Freyer</name>
</author>
<author><name sortKey="Roberts, Ja" uniqKey="Roberts J">JA Roberts</name>
</author>
<author><name sortKey="Ritter, P" uniqKey="Ritter P">P Ritter</name>
</author>
<author><name sortKey="Breakspear, M" uniqKey="Breakspear M">M Breakspear</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hansen, Ec" uniqKey="Hansen E">EC Hansen</name>
</author>
<author><name sortKey="Battaglia, D" uniqKey="Battaglia D">D Battaglia</name>
</author>
<author><name sortKey="Spiegler, A" uniqKey="Spiegler A">A Spiegler</name>
</author>
<author><name sortKey="Deco, G" uniqKey="Deco G">G Deco</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Deco, G" uniqKey="Deco G">G Deco</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Deco, G" uniqKey="Deco G">G Deco</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Llinas, Rr" uniqKey="Llinas R">RR Llinas</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Varela, F" uniqKey="Varela F">F Varela</name>
</author>
<author><name sortKey="Lachaux, Jp" uniqKey="Lachaux J">JP Lachaux</name>
</author>
<author><name sortKey="Rodriguez, E" uniqKey="Rodriguez E">E Rodriguez</name>
</author>
<author><name sortKey="Martinerie, J" uniqKey="Martinerie J">J Martinerie</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Rabinovich, M" uniqKey="Rabinovich M">M Rabinovich</name>
</author>
<author><name sortKey="Huerta, R" uniqKey="Huerta R">R Huerta</name>
</author>
<author><name sortKey="Varona, P" uniqKey="Varona P">P Varona</name>
</author>
<author><name sortKey="Afraimovich, V" uniqKey="Afraimovich V">V Afraimovich</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rabinovich, Mi" uniqKey="Rabinovich M">MI Rabinovich</name>
</author>
<author><name sortKey="Huerta, R" uniqKey="Huerta R">R Huerta</name>
</author>
<author><name sortKey="Laurent" uniqKey="Laurent">Laurent.</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huys, R" uniqKey="Huys R">R Huys</name>
</author>
<author><name sortKey="Perdikis, D" uniqKey="Perdikis D">D Perdikis</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Perdikis, D" uniqKey="Perdikis D">D Perdikis</name>
</author>
<author><name sortKey="Huys, R" uniqKey="Huys R">R Huys</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Beckmann, Cf" uniqKey="Beckmann C">CF Beckmann</name>
</author>
<author><name sortKey="Deluca, M" uniqKey="Deluca M">M DeLuca</name>
</author>
<author><name sortKey="Devlin, Jt" uniqKey="Devlin J">JT Devlin</name>
</author>
<author><name sortKey="Smith, Sm" uniqKey="Smith S">SM Smith</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Greicius, Md" uniqKey="Greicius M">MD Greicius</name>
</author>
<author><name sortKey="Krasnow, B" uniqKey="Krasnow B">B Krasnow</name>
</author>
<author><name sortKey="Reiss, Al" uniqKey="Reiss A">AL Reiss</name>
</author>
<author><name sortKey="Menon, V" uniqKey="Menon V">V Menon</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mueller, S" uniqKey="Mueller S">S Mueller</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Allen, Ea" uniqKey="Allen E">EA Allen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Baker, Ap" uniqKey="Baker A">AP Baker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ghosh, A" uniqKey="Ghosh A">A Ghosh</name>
</author>
<author><name sortKey="Rho, Y" uniqKey="Rho Y">Y Rho</name>
</author>
<author><name sortKey="Mcintosh, Ar" uniqKey="Mcintosh A">AR McIntosh</name>
</author>
<author><name sortKey="Kotter, R" uniqKey="Kotter R">R Kotter</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ponce Alvarez, A" uniqKey="Ponce Alvarez A">A Ponce-Alvarez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gusnard, Da" uniqKey="Gusnard D">DA Gusnard</name>
</author>
<author><name sortKey="Raichle, Me" uniqKey="Raichle M">ME Raichle</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Addis, Dr" uniqKey="Addis D">DR Addis</name>
</author>
<author><name sortKey="Wong, At" uniqKey="Wong A">AT Wong</name>
</author>
<author><name sortKey="Schacter, Dl" uniqKey="Schacter D">DL Schacter</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Botvinick, M" uniqKey="Botvinick M">M Botvinick</name>
</author>
<author><name sortKey="Nystrom, Le" uniqKey="Nystrom L">LE Nystrom</name>
</author>
<author><name sortKey="Fissell, K" uniqKey="Fissell K">K Fissell</name>
</author>
<author><name sortKey="Carter, Cs" uniqKey="Carter C">CS Carter</name>
</author>
<author><name sortKey="Cohen, Jd" uniqKey="Cohen J">JD Cohen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Goodkind, M" uniqKey="Goodkind M">M Goodkind</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Deco, G" uniqKey="Deco G">G Deco</name>
</author>
<author><name sortKey="Kringelbach, Ml" uniqKey="Kringelbach M">ML Kringelbach</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ritter, P" uniqKey="Ritter P">P Ritter</name>
</author>
<author><name sortKey="Schirner, M" uniqKey="Schirner M">M Schirner</name>
</author>
<author><name sortKey="Mcintosh, Ar" uniqKey="Mcintosh A">AR McIntosh</name>
</author>
<author><name sortKey="Jirsa, Vk" uniqKey="Jirsa V">VK Jirsa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sanz Leon, P" uniqKey="Sanz Leon P">P Sanz Leon</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ritter, P" uniqKey="Ritter P">P Ritter</name>
</author>
<author><name sortKey="Villringer, A" uniqKey="Villringer A">A Villringer</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Desikan, Rs" uniqKey="Desikan R">RS Desikan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tournier, Jd" uniqKey="Tournier J">JD Tournier</name>
</author>
<author><name sortKey="Calamante, F" uniqKey="Calamante F">F Calamante</name>
</author>
<author><name sortKey="Gadian, Dg" uniqKey="Gadian D">DG Gadian</name>
</author>
<author><name sortKey="Connelly, A" uniqKey="Connelly A">A Connelly</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Matthews, Pc" uniqKey="Matthews P">PC Matthews</name>
</author>
<author><name sortKey="Strogatz, Sh" uniqKey="Strogatz S">SH Strogatz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aronson, Dg" uniqKey="Aronson D">DG Aronson</name>
</author>
<author><name sortKey="Ermentrout, Gb" uniqKey="Ermentrout G">GB Ermentrout</name>
</author>
<author><name sortKey="Kopell, N" uniqKey="Kopell N">N Kopell</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Glerean, E" uniqKey="Glerean E">E Glerean</name>
</author>
<author><name sortKey="Salmi, J" uniqKey="Salmi J">J Salmi</name>
</author>
<author><name sortKey="Lahnakoski, Jm" uniqKey="Lahnakoski J">JM Lahnakoski</name>
</author>
<author><name sortKey="Jaaskelainen, Ip" uniqKey="Jaaskelainen I">IP Jaaskelainen</name>
</author>
<author><name sortKey="Sams, M" uniqKey="Sams M">M Sams</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Biswal, B" uniqKey="Biswal B">B Biswal</name>
</author>
<author><name sortKey="Yetkin, Fz" uniqKey="Yetkin F">FZ Yetkin</name>
</author>
<author><name sortKey="Haughton, Vm" uniqKey="Haughton V">VM Haughton</name>
</author>
<author><name sortKey="Hyde, Js" uniqKey="Hyde J">JS Hyde</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Buckner, Rl" uniqKey="Buckner R">RL Buckner</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Achard, S" uniqKey="Achard S">S Achard</name>
</author>
<author><name sortKey="Salvador, R" uniqKey="Salvador R">R Salvador</name>
</author>
<author><name sortKey="Whitcher, B" uniqKey="Whitcher B">B Whitcher</name>
</author>
<author><name sortKey="Suckling, J" uniqKey="Suckling J">J Suckling</name>
</author>
<author><name sortKey="Bullmore, E" uniqKey="Bullmore E">E Bullmore</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
  <front><journal-meta><journal-id journal-id-type="nlm-ta">Sci Rep</journal-id>
<journal-id journal-id-type="iso-abbrev">Sci Rep</journal-id>
<journal-title-group><journal-title>Scientific Reports</journal-title>
</journal-title-group>
<issn pub-type="epub">2045-2322</issn>
<publisher><publisher-name>Nature Publishing Group UK</publisher-name>
<publisher-loc>London</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">28596608</article-id>
<article-id pub-id-type="pmc">5465179</article-id>
<article-id pub-id-type="publisher-id">3073</article-id>
<article-id pub-id-type="doi">10.1038/s41598-017-03073-5</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>The dynamics of resting fluctuations in the brain: metastability and its dynamical cortical core</article-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Deco</surname>
<given-names>Gustavo</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
<xref ref-type="aff" rid="Aff2">2</xref>
<xref ref-type="aff" rid="Aff3">3</xref>
<xref ref-type="aff" rid="Aff4">4</xref>
</contrib>
<contrib contrib-type="author" corresp="yes"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0002-3908-6898</contrib-id>
<name><surname>Kringelbach</surname>
<given-names>Morten L.</given-names>
</name>
<address><email>morten.kringelbach@psych.ox.ac.uk</email>
</address>
<xref ref-type="aff" rid="Aff5">5</xref>
<xref ref-type="aff" rid="Aff6">6</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Jirsa</surname>
<given-names>Viktor K.</given-names>
</name>
<xref ref-type="aff" rid="Aff7">7</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Ritter</surname>
<given-names>Petra</given-names>
</name>
<xref ref-type="aff" rid="Aff8">8</xref>
<xref ref-type="aff" rid="Aff9">9</xref>
</contrib>
<aff id="Aff1"><label>1</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2172 2676</institution-id>
<institution-id institution-id-type="GRID">grid.5612.0</institution-id>
<institution>Center for Brain and Cognition, Computational Neuroscience Group, Department of Information and Communication Technologies,</institution>
<institution>Universitat Pompeu Fabra,</institution>
</institution-wrap>
Roc Boronat 138, Barcelona, 08018 Spain</aff>
<aff id="Aff2"><label>2</label>
Institució Catalana de la Recerca i Estudis Avançats (ICREA), Passeig Lluís Companys 23, Barcelona, 08010 Spain</aff>
<aff id="Aff3"><label>3</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0041 5028</institution-id>
<institution-id institution-id-type="GRID">grid.419524.f</institution-id>
<institution>Department of Neuropsychology,</institution>
<institution>Max Planck Institute for Human Cognitive and Brain Sciences,</institution>
</institution-wrap>
04103 Leipzig, Germany</aff>
<aff id="Aff4"><label>4</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 7857</institution-id>
<institution-id institution-id-type="GRID">grid.1002.3</institution-id>
<institution>School of Psychological Sciences,</institution>
<institution>Monash University,</institution>
</institution-wrap>
Melbourne, Clayton VIC 3800, Australia</aff>
<aff id="Aff5"><label>5</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 8948</institution-id>
<institution-id institution-id-type="GRID">grid.4991.5</institution-id>
<institution>Department of Psychiatry,</institution>
<institution>University of Oxford,</institution>
</institution-wrap>
Oxford, UK</aff>
<aff id="Aff6"><label>6</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1956 2722</institution-id>
<institution-id institution-id-type="GRID">grid.7048.b</institution-id>
<institution>Center of Music in the Brain (MIB), Clinical Medicine,</institution>
<institution>Aarhus University,</institution>
</institution-wrap>
Aarhus C, DK Denmark</aff>
<aff id="Aff7"><label>7</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0541 5643</institution-id>
<institution-id institution-id-type="GRID">grid.462494.9</institution-id>
<institution>Aix Marseille Univ, INSERM,</institution>
<institution>INS, Inst Neurosci Syst,</institution>
</institution-wrap>
Marseille, France</aff>
<aff id="Aff8"><label>8</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0001 2105 1091</institution-id>
<institution-id institution-id-type="GRID">grid.4372.2</institution-id>
<institution></institution>
<institution>Max-Planck Institute for Cognitive and Brain Sciences,</institution>
</institution-wrap>
Leipzig, Germany</aff>
<aff id="Aff9"><label>9</label>
Department of Neurology, Charité, Charitéplatz 1, 10117 Berlin, Germany</aff>
</contrib-group>
<pub-date pub-type="epub"><day>8</day>
<month>6</month>
<year>2017</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>8</day>
<month>6</month>
<year>2017</year>
</pub-date>
<pub-date pub-type="collection"><year>2017</year>
</pub-date>
<volume>7</volume>
<elocation-id>3095</elocation-id>
<history><date date-type="received"><day>16</day>
<month>10</month>
<year>2015</year>
</date>
<date date-type="accepted"><day>6</day>
<month>4</month>
<year>2017</year>
</date>
</history>
<permissions><copyright-statement>© The Author(s) 2017</copyright-statement>
<license license-type="OpenAccess"><license-p><bold>Open Access</bold>
 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 <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>
.</license-p>
</license>
</permissions>
<abstract id="Abs1"><p>In the human brain, spontaneous activity during resting state consists of rapid transitions between functional network states over time but the underlying mechanisms are not understood. We use connectome based computational brain network modeling to reveal fundamental principles of how the human brain generates large-scale activity observable by noninvasive neuroimaging. We used structural and functional neuroimaging data to construct whole- brain models. With this novel approach, we reveal that the human brain during resting state operates at maximum metastability, i.e. in a state of maximum network switching. In addition, we investigate cortical heterogeneity across areas. Optimization of the spectral characteristics of each local brain region revealed the dynamical cortical core of the human brain, which is driving the activity of the rest of the whole brain. Brain network modelling goes beyond correlational neuroimaging analysis and reveals non-trivial network mechanisms underlying non-invasive observations. Our novel findings significantly pertain to the important role of computational connectomics in understanding principles of brain function.</p>
</abstract>
<custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name>
<meta-value>© The Author(s) 2017</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body><sec id="Sec1" sec-type="introduction"><title>Introduction</title>
<p><disp-quote><p><italic>“When we take a general view of the wonderful stream of our consciousness, what strikes us first is the different pace of its parts. Like a bird’s life, it seems to be made of an alternation of flights and perchings.” William James</italic>
<sup><xref ref-type="bibr" rid="CR1">1</xref>
</sup>
</p>
</disp-quote>
</p>
<p>Survival remains the perhaps most important problem faced by brains and a key challenge is how to segregate and integrate relevant information over different timescales when faced with hostile, often constantly changing environments<sup><xref ref-type="bibr" rid="CR2">2</xref>
</sup>
. Reconciling different speeds of information processing, from fast to slow, is especially important, and could be key to the relative evolutionary success of mammals whose sophisticated brains are able to combine prior information from past memories with current stimuli to predict the future and to adapt behaviour accordingly<sup><xref ref-type="bibr" rid="CR3">3</xref>
–<xref ref-type="bibr" rid="CR5">5</xref>
</sup>
.</p>
<p>This was recognized well over a century ago by William James, generally acknowledged as one of the fathers of modern cognitive psychology<sup><xref ref-type="bibr" rid="CR1">1</xref>
</sup>
. Speaking of this problem using the apt metaphor of the stream of consciousness, James noted that there is a different pace to its parts, comparing it to the life of a bird whose journey consists of an “alternation of flights and perchings”. In the language of today’s dynamical systems, the flights are akin to fast, segregative tendencies and the perchings to slower, integrative tendencies of the dynamic brain in action<sup><xref ref-type="bibr" rid="CR2">2</xref>
, <xref ref-type="bibr" rid="CR6">6</xref>
, <xref ref-type="bibr" rid="CR7">7</xref>
</sup>
. In addition, motivated by recent experimental and modelling work of other labs<sup><xref ref-type="bibr" rid="CR8">8</xref>
, <xref ref-type="bibr" rid="CR9">9</xref>
</sup>
, we investigate cortical heterogeneity across areas. By optimizing the spectral characteristics of each local brain node (in the coupled network), this allowed us to discover a dynamical core of the brain, i.e. the set of brain regions, which through their oscillations are driving the rest of the brain. Furthermore, with regards to balancing the different speeds of processing, a large body of psychological research has focused on what is known as dual process theories<sup><xref ref-type="bibr" rid="CR10">10</xref>
, <xref ref-type="bibr" rid="CR11">11</xref>
</sup>
, identifying competing fast and slow systems which have to co-exist and function on multiple time-scales in order for the brain to efficiently allocate the resources necessary for survival<sup><xref ref-type="bibr" rid="CR12">12</xref>
, <xref ref-type="bibr" rid="CR13">13</xref>
</sup>
.</p>
<p>Yet, the temporal dynamics and underlying neural mechanisms of this temporal processing on multiple timescales are poorly understood. Here we aim to provide a better understanding of the dynamics using computational brain network modelling which has emerged as a powerful tool for investigating the <italic>causal</italic>
 dynamics of the human brain, when carefully constrained by functional (FC) and structural connectivity (SC) obtained from empirical neuroimaging data<sup><xref ref-type="bibr" rid="CR14">14</xref>
–<xref ref-type="bibr" rid="CR18">18</xref>
</sup>
. This theoretical framework has been largely successful in explaining the highly structured dynamics arising from spontaneous brain activity in the so-called resting-state-networks (RSN)<sup><xref ref-type="bibr" rid="CR19">19</xref>
–<xref ref-type="bibr" rid="CR21">21</xref>
</sup>
, even if the resting brain never truly rests<sup><xref ref-type="bibr" rid="CR20">20</xref>
</sup>
. Efficient task-related brain activity has been shown to rely on <italic>metastability</italic>
 of spontaneous brain activity allowing for optimal exploration of the dynamical repertoire<sup><xref ref-type="bibr" rid="CR22">22</xref>
</sup>
 but it is not known if this metastability is maximally metastable<sup><xref ref-type="bibr" rid="CR6">6</xref>
</sup>
. In dynamic systems, metastability refers to a state that falls outside the natural equilibrium state of the system but persists for an extended period of time. One example of a metastable dynamical system is a winnerless competition<sup><xref ref-type="bibr" rid="CR23">23</xref>
</sup>
, however metastability can arise from a number of underlying mechanisms and it is in this broader sense that we use the term metastability. Here, we use metastability to denote the variability of the global synchronization as measured by the Kuramoto order parameter, following the work of Wildie and Shanahan<sup><xref ref-type="bibr" rid="CR24">24</xref>
</sup>
. This is motivated by previous research which has shown that this variability can be linked to the underlying metastable cluster synchronization<sup><xref ref-type="bibr" rid="CR25">25</xref>
</sup>
.</p>
<p>We investigated the dynamics of the brain network system through a local node neural mass description based on the most general form of expressing both noisy asynchronous dynamics and oscillations, namely a normal form of a Hopf bifurcation<sup><xref ref-type="bibr" rid="CR26">26</xref>
–<xref ref-type="bibr" rid="CR28">28</xref>
</sup>
. Previous research has shown the usefulness, richness and generality of this type of model for describing EEG dynamics at the local node level<sup><xref ref-type="bibr" rid="CR27">27</xref>
, <xref ref-type="bibr" rid="CR28">28</xref>
</sup>
. This normal form allowed us to fit the model to neuroimaging data <italic>over time</italic>
, i.e. not only by fitting the grand average FC but also by fitting the temporal structure of the fluctuations, functional connectivity dynamics<sup><xref ref-type="bibr" rid="CR29">29</xref>
</sup>
 (FCD, Fig. <xref rid="Fig1" ref-type="fig">1A,B</xref>
).<fig id="Fig1"><label>Figure 1</label>
<caption><p>Methods for measuring fit between simulated and empirical data. (<bold>A</bold>
) The fitting of the FC is measured by the Pearson correlation coefficient between corresponding elements of the upper triangular part of the matrices. (<bold>B</bold>
) For comparing the FCD statistics, we collected the upper triangular elements of the matrices (over all participants or sessions) and compared the simulated and empirical distribution by means of the Kolmogorov-Smirnov distance between them. The Kolmogorov–Smirnov distance quantifies the maximal difference between the cumulative distribution functions of the two samples. (<bold>C</bold>
) We measure the metastability as the standard deviation of the Kuramoto order parameter across time. The Kuramoto order parameter measures the global level of synchronization of the n oscillating signals. Under complete independence, the n phases are uniformly distributed and thus R is nearly zero, whereas R = 1 if all phases are equal (full synchronization). For calculating the metastability of the empirical and simulated BOLD signals, we first band-pass filtered within the narrowband 0.04–0.07 Hz and computed the instantaneous phase φ<sub>k</sub>
(t) of each narrowband signal k using the Hilbert transform. The Hilbert transform yields the associated analytical signals. The analytic signal represents a narrowband signal, s(t), in the time domain as a rotating vector with an instantaneous phase, φ(t), and an instantaneous amplitude, A(t). Bottom panel visualizes a single example scenario (of many possible others) where the model system’s metastability increases as a function of G. We also indicate the metastability measured in empirical data. Part of figure B is based on the work of Allen and colleagues<sup><xref ref-type="bibr" rid="CR44">44</xref>
</sup>
.</p>
</caption>
<graphic xlink:href="41598_2017_3073_Fig1_HTML" id="d29e476"></graphic>
</fig>
</p>
<p>We further explored if the optimal working point where FC and FCD are fitted corresponds to a dynamical region where the global metastability of the whole brain is maximized<sup><xref ref-type="bibr" rid="CR6">6</xref>
</sup>
. In addition, motivated by recent experimental and modelling work of other labs<sup><xref ref-type="bibr" rid="CR8">8</xref>
, <xref ref-type="bibr" rid="CR9">9</xref>
</sup>
, we investigate cortical heterogeneity across areas. By optimizing the spectral characteristics of each local brain node (in the coupled network), this allowed us discovering a dynamical core of the brain, i.e. the set of brain regions, which through their oscillations is driving the rest of the brain. As such this investigation was designed to provide an empirical, scientific footing for James’ metaphorical speculations of the flights and perchings of human brain dynamics, and to demonstrate the potential of sophisticated brain network computational modelling to provide new insights into the causal mechanisms of neuroimaging results.</p>
</sec>
<sec id="Sec2" sec-type="results"><title>Results</title>
<p>The results arose from using brain network computational models for the analysis of empirical neuroimaging data characterising the functional and structural connectivity of 24 healthy human participants acquired using standard MRI techniques<sup><xref ref-type="bibr" rid="CR17">17</xref>
</sup>
 (see <italic>Methods</italic>
). In particular, we were able to gain new insights on the emergence of transiently spatiotemporal structured networks among segregated brain regions by examining a whole-brain network model using a very general neural mass model known as the normal form of a Hopf bifurcation (also known as Landau-Stuart Oscillators), which is the canonical model for studying the transition from noisy to oscillatory dynamics<sup><xref ref-type="bibr" rid="CR26">26</xref>
</sup>
 (Fig. <xref rid="Fig2" ref-type="fig">2</xref>
). Here, we extended previous research on local node dynamics<sup><xref ref-type="bibr" rid="CR27">27</xref>
, <xref ref-type="bibr" rid="CR28">28</xref>
</sup>
 by studying the whole-brain network dynamics, i.e. by investigating how those local noisy oscillators interact, and how the emerging whole-brain network activity relates to fMRI resting state dynamics. Within this model, each node of the network is modeled by a normal Hopf bifurcation, with an intrinsic frequency <italic>ω</italic>
<sub><italic>i</italic>
</sub>
 in the 0.04–0.07 Hz band (<italic>i</italic>
 = 1, …,<italic>n</italic>
). The intrinsic frequencies were estimated directly from the data, as given by the averaged peak frequency of the narrowband BOLD signals of each brain region (see <italic>Methods</italic>
). The state of each node <italic>i</italic>
 is determined by its phase, <italic>φ</italic>
<sub><italic>i</italic>
</sub>
(<italic>t</italic>
), and the interaction between nodes depends both on the structural couplings and the phase difference between the nodes. The model has only two types of control parameters, namely: one single global parameter, <italic>G</italic>
, that represents the global scaling of the anatomical connectivity matrix, and the bifurcation parameters <italic>a</italic>
<sub><italic>j</italic>
</sub>
 for each node (see Fig. <xref rid="Fig3" ref-type="fig">3</xref>
 and methods for the general structure and strategy of the brain network model).<fig id="Fig2"><label>Figure 2</label>
<caption><p>Construction of individual brain network models. (<bold>A</bold>
) The brain network model was based on individual structural connectivity (SC) matrices from 24 participants derived from tractography of DTI (left) between the 68 regions of the Desikan-Kahilly parcellation (middle). The control parameters of the models were tuned using the grand average FC and FCD derived from fMRI BOLD data (right). (<bold>B</bold>
) For modelling local neural masses we used the normal form of a Hopf bifurcation, where depending on the bifurcation parameter, the local model generates a noisy signal (left), a mixed noisy and oscillatory signal (middle) or an oscillatory signal (right). It is at the border between noisy and oscillatory behaviour (middle), where the simulated signal looks like the empirical data, i.e. like noise with an oscillatory component around 0.05 Hz.</p>
</caption>
<graphic xlink:href="41598_2017_3073_Fig2_HTML" id="d29e576"></graphic>
</fig>
<fig id="Fig3"><label>Figure 3</label>
<caption><p>Fitting of the empirical data by the brain network Hopf model for different working points. (<bold>A</bold>
) Level of fitting of the FC, FCD and metastability as a function of the global scaling parameter G for three different bifurcation parameters a = [−0.2 0 0.2], namely at the noisy oscillatory region, at the edge of the bifurcation and at the oscillatory regime. (<bold>B</bold>
) The three measures for assessing fitting between simulated and empirical data are shown color-coded as a function of bifurcation parameter a and global scaling parameter, G. The best fitting of the three measures is obtained for a region at the brink of the Hopf bifurcation, i.e. for bifurcation parameter a, at the edge of zero on the negative side. In this region not only the correlation between the empirical and simulated FC is maximized (upper panel), but also the statistics of the rapid switching between FC(t) across time (FCD) is minimized in Kolmogorov-Smirnov sense (middle panel), and the level of metastability of the data is perfectly reproduced (bottom panel).</p>
</caption>
<graphic xlink:href="41598_2017_3073_Fig3_HTML" id="d29e592"></graphic>
</fig>
</p>
<sec id="Sec3"><title>Maximal metastability at the optimal working point of model</title>
<p>Using the Hopf model, we were able discern the dynamical properties of the optimal working point of the system that is able to fit the characteristics of the empirical fMRI data. We were able to distinguish the origin of resting activity between the two hypothesized scenarios, namely: 1) noisy excursions at the edge of a critical bifurcation<sup><xref ref-type="bibr" rid="CR19">19</xref>
, <xref ref-type="bibr" rid="CR20">20</xref>
, <xref ref-type="bibr" rid="CR30">30</xref>
, <xref ref-type="bibr" rid="CR31">31</xref>
</sup>
 or 2) metastable oscillations<sup><xref ref-type="bibr" rid="CR16">16</xref>
</sup>
. The first scenario refers to the entrainment of noisy dynamics through the underlying anatomical connectivity matrix, i.e. inducing correlations of the local noise because of the underlying SC connections. The second scenario refers to the structuring of metastable cluster synchronizations of the underlying local oscillatory dynamics through the underlying anatomical SC connections. We define metastability as the standard deviation of synchrony at the network level described by order parameter R(t), where R(t) measures the phase uniformity and varies between 0 for a fully desynchronized network and 1 for a fully synchronized network (see methods and Fig. <xref rid="Fig1" ref-type="fig">1C</xref>
)<sup><xref ref-type="bibr" rid="CR24">24</xref>
</sup>
. The present model is able to describe both types of dynamics, and the smooth transitions from one to the other, i.e. the transition from noisy to oscillatory dynamics (Fig. <xref rid="Fig2" ref-type="fig">2</xref>
). In order to distinguish the dynamical scenario, we investigated the capabilities of the model for fitting the grand average FC and also the time dependent characteristics of the RSN as reflected in the FCD in the different dynamical working regions (i.e. as a function of the control parameters). The grand average FC describes the mean spatial structure of the resting activity, whereas the FCD captures the statistical characteristic of the temporal structure of those spatial correlations (see <italic>Methods</italic>
 and ref. <xref ref-type="bibr" rid="CR29">29</xref>
).</p>
<p>Figure <xref rid="Fig3" ref-type="fig">3</xref>
 shows that the best fit to the empirical data of Hopf model is found at the brink of the Hopf bifurcation. We equalized all local bifurcation parameters to a common value i.e. <italic>a</italic>
<sub><italic>j</italic>
</sub>
 = <italic>a</italic>
, in order to reduce the investigations to just two parameters, namely global bifurcation parameter (<italic>a</italic>
) and global coupling strength (<italic>G</italic>
). Figure <xref rid="Fig3" ref-type="fig">3</xref>
 shows how the empirical data are fitted in the Hopf model for different working points. The right column of Fig. <xref rid="Fig3" ref-type="fig">3</xref>
 shows the level of fitting of the FC, FCD and metastability. As can be seen, the best fitting of the three measures is obtained at the region on the brink of the Hopf bifurcation, i.e. for bifurcation parameter <italic>a</italic>
, at the edge of zero on the negative side, such that the oscillators remain damped still. In this region not only the correlation between the empirical and simulated FC is maximized, but also the statistics of the rapid switching between FC(t) across time (FCD) is minimized in Kolmogorov-Smirnov sense, and the level of metastability of the data is reproduced. The fitting of the FC was measured by the Pearson correlation coefficient between corresponding elements of the upper triangular part of the matrices (see Fig. <xref rid="Fig1" ref-type="fig">1</xref>
 and <italic>Methods</italic>
). For comparing the FCD statistics, we collected the upper triangular elements of the matrices (over all participants or sessions) and compared the simulated and empirical distribution by means of the Kolmogorov-Smirnov distance between them (see <italic>Methods</italic>
).</p>
<p>Furthermore, the results showed that only in the region at the border between noisy and oscillatory behaviour, is where the signals resembles the data, i.e. like noise with an oscillatory component around 0.05 Hz (Fig. <xref rid="Fig2" ref-type="fig">2</xref>
). The first three columns of Fig. <xref rid="Fig3" ref-type="fig">3</xref>
 show the dependence of those measurements as a function of the global scaling parameter <italic>G</italic>
 for three specific values of the bifurcation parameter <italic>a</italic>
, namely at the noisy region, at the edge of the bifurcation and at the oscillatory regime. Clearly, the best results are obtained for the second column (at the edge of the bifurcation). The same panel shows that the FCD is the best constraining measure. There is a broad range of G where the FC and the metastability is well fitted, but only a relative narrow range where the FCD statistics is minimal, i.e. maximally fitted. In other words, the spatiotemporal structure of the FC is more informative than the grand average of the FC (i.e. the “classical” RSN). This is important, because until now, brain network models have always been fitted with the grand average FC - but see also the work of Hansen and colleagues<sup><xref ref-type="bibr" rid="CR29">29</xref>
</sup>
.</p>
<p>We would like to remark that Fig. <xref rid="Fig3" ref-type="fig">3</xref>
 characterizes some of the bifurcation behaviour of the whole system. Indeed, the metastability for example serves as a network metric and characterizes the variability of this global synchronization as a function of those two control parameters. All three parameter spaces in Fig. <xref rid="Fig3" ref-type="fig">3B</xref>
, in conjunction, present a full picture of the spatiotemporal organization of the system. The three metrics characterize computationally the bifurcation properties of the full network dynamics.</p>
<p>Perhaps most importantly, as shown in Fig. <xref rid="Fig3" ref-type="fig">3</xref>
, the brain network model shows maximal <italic>metastability</italic>
 at the optimal working point of the model (<italic>a</italic>
 = 0 and <italic>G</italic>
 = 2.85), where the metastability is reflecting the variability of the synchronization between different nodes, i.e. the fluctuations of the states of phase configurations as a function of time<sup><xref ref-type="bibr" rid="CR24">24</xref>
</sup>
. Further characterisation of these results is shown in Fig. <xref rid="Fig4" ref-type="fig">4</xref>
 which shows the optimal working point at the edge of the Hopf bifurcation (i.e. bifurcation parameter <italic>a</italic>
 = 0), the FC, FCD and FCD statistics for three levels of global coupling <italic>G</italic>
 namely low, optimal and large. For comparison, the same matrices and distributions are plotted on the rightmost column for the empirical data (Fig. <xref rid="Fig4" ref-type="fig">4B</xref>
). Only the FCD and its statistics (bottom row) are constraining enough for optimizing the working point. Please note that for low <italic>G</italic>
 the FCD statistics does not show any switching between states in the RSN and that for very large <italic>G</italic>
 there are too much switching between states.<fig id="Fig4"><label>Figure 4</label>
<caption><p>Fitting to the grand average FC is a necessary but not sufficient condition for best empirical fitting. (<bold>A</bold>
) The figure shows the result of fitting the model to the empirical as a function of the global coupling parameter, G, at the optimal working point at the edge of the Hopf bifurcation (i.e. bifurcation parameter a = 0). Three different coupling points were selected (low, optimal and large in the three columns) and we show the resulting FC correlation, FCD correlations and FCD histogram. Note that for low G the FCD statistics does not show any switching between RSN and that for very large G there are too much switching between states. (<bold>B</bold>
) For comparison, the same matrices and distributions are plotted for the empirical data. Note how only the FCD (row 2) and its statistics (row 3) are constraining enough for optimizing the working point the model to fit the empirical data (compare the distributions in row 3 and compare plots for FCD and FC fitting in row 1).</p>
</caption>
<graphic xlink:href="41598_2017_3073_Fig4_HTML" id="d29e752"></graphic>
</fig>
</p>
</sec>
<sec id="Sec4"><title>Dynamical core: contribution of individual brain regions to dynamics</title>
<p>In order to obtain information about the dynamical characteristics of each single brain area and to generate a heterogeneous brain network model (i.e. with different dynamics at each node), we optimized each single bifurcation parameter <italic>a</italic>
<sub><italic>j</italic>
</sub>
 independently by fitting for each value of global coupling <italic>G</italic>
 the spectral characteristics of the simulated and empirical BOLD signals at each brain area (see <italic>Methods</italic>
). The main results are plotted in Fig. <xref rid="Fig5" ref-type="fig">5</xref>
, where Fig. <xref rid="Fig5" ref-type="fig">5A</xref>
 shows the evolution of the fitting of the FC and FCD statistics as a function of <italic>G</italic>
. For large enough value of the global coupling a good fitting of both is obtained, i.e. large correlation between the empirical and simulated grand average FC and low difference in the statistics of the empirical and simulated FCD (Kolmogorov-Smirnov distance). Please note that Fig. <xref rid="Fig5" ref-type="fig">5A</xref>
 is generated in a different way than Fig. <xref rid="Fig4" ref-type="fig">4A</xref>
 (which uses only the optimum fit a = 0 for all regions and G = 2.85). Instead, in Fig. <xref rid="Fig5" ref-type="fig">5A</xref>
, for each G we optimize the bifurcation value, <italic>a</italic>
, for <italic>each</italic>
 region (shown in Fig. <xref rid="Fig5" ref-type="fig">5B</xref>
). As can be seen at a critical value of G, the bifurcation values remain the same, only scaled. Thus the FCD fit in Fig. <xref rid="Fig5" ref-type="fig">5A</xref>
 will asymptote as G increases.<fig id="Fig5"><label>Figure 5</label>
<caption><p>Spectral characteristics of the dynamical core of the human brain. To generate a heterogeneous brain network model (i.e. with different dynamics at each node), we optimized each single bifurcation parameter independently by fitting for each value of global coupling G the spectral characteristics of the simulated and empirical BOLD signals at each brain area. (<bold>A</bold>
) The evolution of the fitting of the FC and FCD statistics as a function of G. For large enough value of the global coupling a good fitting of both is obtained, i.e. large correlation between the empirical and simulated grand average FC and low difference in the statistics of the empirical and simulated FCD (Kolmogorov-Smirnov distance). (<bold>B</bold>
) The evolution of the single values of the local bifurcations parameters as a function of the global coupling G. For low values of G homogeneous local bifurcation parameters around zero are obtained. When the level of fitting improves for larger values of G a more heterogeneous distribution of is obtained. (<bold>C</bold>
) The local bifurcation parameters for each region for the uncoupled network (i.e. G = 0) and for the optimal coupling (G = 5.4). If the network is uncoupled, each single brain area fitted the spectral characteristics of the empirical BOLD signals in a very homogeneous way by local bifurcations parameters at the edge of the local Hopf bifurcation, i.e. at zero. (<bold>D</bold>
) When the whole-brain network is coupled, we can discover the “true” intrinsic local dynamics that fits the local empirical BOLD characteristics and the global quantities FC, FCD and metastability.</p>
</caption>
<graphic xlink:href="41598_2017_3073_Fig5_HTML" id="d29e824"></graphic>
</fig>
</p>
<p>For optimizing <italic>a</italic>
<sub><italic>j</italic>
</sub>
 values, we use a greedy optimisation strategy, where we iteratively increase or decrease the <italic>a</italic>
<sub><italic>j</italic>
</sub>
 value according to the local power of the signal in a given region <italic>j</italic>
. Greedy algorithms exploit local optima, but often approximate optimal solutions well in reasonable time and produce good results as shown in Fig. <xref rid="Fig5" ref-type="fig">5</xref>
. The local bifurcation parameters for each region for the uncoupled network (i.e. <italic>G</italic>
 = 0) and for the optimal coupling (<italic>G</italic>
 = 5.4) can be seen in Fig. <xref rid="Fig5" ref-type="fig">5C</xref>
. If the network is uncoupled, each single brain area fitted the spectral characteristics of the empirical BOLD signals in a very homogeneous way by local bifurcations parameters at the edge of the local Hopf bifurcation, i.e. at zero. When the brain network is coupled, the “true” intrinsic local dynamics for the profile of optimal local bifurcation parameters <italic>a</italic>
<sub><italic>j</italic>
</sub>
 observed at that point that fit the local empirical BOLD characteristics and the global quantities FC, FCD and metastability (Fig. <xref rid="Fig5" ref-type="fig">5D</xref>
).</p>
<p>Brain regions, for which best predictions were achieved in an oscillatory mode, i.e. with bifurcation parameters <italic>a</italic>
 > 0.1 are visualised in Fig. <xref rid="Fig6" ref-type="fig">6</xref>
. We found that the dynamical core within this parcellation consisted of eight lateralised brain regions: medial orbitofrontal cortex, posterior cingulate cortex and transverse temporal gyrus in the right hemisphere, and caudal middle frontal gyrus, precentral gyrus, precuneus cortex, rostral anterior cingulate cortex and transverse temporal gyrus in the left hemisphere. Those nodes working at the edge of the bifurcation are highlighted as a “dynamical core” whose perturbations can propagate in an optimal way to the rest of the network.<fig id="Fig6"><label>Figure 6</label>
<caption><p>Dynamical core in the human brain. The figure shows the dynamical core regions on the edge of bifurcation (location of neural masses shown in light blue and transparent blue for the full region). These are the nodes with the ability to react immediately to changes in the predicted input and thus likely to drive the rest of the brain networks. The eight regions are clearly lateralised; and in the right hemisphere encompass medial orbitofrontal cortex, posterior cingulate cortex and transverse temporal gyrus, while in the left hemisphere include caudal middle frontal gyrus, precentral gyrus, precuneus cortex, rostral anterior cingulate cortex and transverse temporal gyrus. Interestingly, some of these regions are part of the default mode network (medial orbitofrontal cortex, posterior cingulate cortex and precuneus cortex) while others have been implicated in memory processing (parahippocampal and transverse temporal gyrus), auditory processing (transverse temporal gyrus), selection for action (rostral anterior cingulate cortex and caudal middle frontal gyrus) and motor execution (precentral gyrus).</p>
</caption>
<graphic xlink:href="41598_2017_3073_Fig6_HTML" id="d29e885"></graphic>
</fig>
</p>
</sec>
</sec>
<sec id="Sec5" sec-type="discussion"><title>Discussion</title>
<p>We provide mechanistic explanations of the complex spatiotemporal dynamics of brain function arising from James’ early speculations<sup><xref ref-type="bibr" rid="CR1">1</xref>
</sup>
 to much more detailed scientific enquiry<sup><xref ref-type="bibr" rid="CR2">2</xref>
, <xref ref-type="bibr" rid="CR32">32</xref>
–<xref ref-type="bibr" rid="CR34">34</xref>
</sup>
. This confirms that brain results from complex interactions in a system of non-linearly coupled, non-linear oscillatory processes which display dynamical system phenomena such as multiple stable states, instability, state transitions and metastability, of which the latter has been proposed to form a core dynamical description of coordinated brain and behavioral activity<sup><xref ref-type="bibr" rid="CR6">6</xref>
</sup>
.</p>
<p>In the 1980s the physicist Hermann Haken suggested to mechanistically interpret brain processes of segregation and integration as a sequence of semistable states, so-called saddle states<sup><xref ref-type="bibr" rid="CR35">35</xref>
</sup>
. He proposed to view the complex integrative and segregative tendencies as expressions of emergent lower-dimensional behavior of collective variables, which he termed ‘order parameters’. Scott Kelso popularized this concept using the term ‘metastability’ based on his brain-behaviour experiments and drawing inspiration from other researchers including Rodolfo Llinás and Francisco Varela<sup><xref ref-type="bibr" rid="CR33">33</xref>
, <xref ref-type="bibr" rid="CR34">34</xref>
</sup>
. He generalized metastability to include the oscillatory states of brain processes found in between complete synchronization and independence<sup><xref ref-type="bibr" rid="CR32">32</xref>
, <xref ref-type="bibr" rid="CR36">36</xref>
</sup>
. Later research has formalized these concepts more rigorously, e.g. via the heteroclinic channel<sup><xref ref-type="bibr" rid="CR37">37</xref>
, <xref ref-type="bibr" rid="CR38">38</xref>
</sup>
 and Structured Flows on Manifolds (SFM)<sup><xref ref-type="bibr" rid="CR39">39</xref>
, <xref ref-type="bibr" rid="CR40">40</xref>
</sup>
.</p>
<p>We shed new causal light on the mechanisms underlying RSNs by extending previous research which has demonstrated the existence of RSNs, i.e. brain networks correlated within the grand average FC during resting state<sup><xref ref-type="bibr" rid="CR21">21</xref>
, <xref ref-type="bibr" rid="CR41">41</xref>
, <xref ref-type="bibr" rid="CR42">42</xref>
</sup>
. FC has become routinely used as a biomarker in various clinical applications, even though its predictive value holds only for group analyses, and not currently for the individual<sup><xref ref-type="bibr" rid="CR43">43</xref>
</sup>
. This problem arises most likely from the lack of taking time into account, i.e. the non-stationary nature of the resting state dynamics<sup><xref ref-type="bibr" rid="CR44">44</xref>
, <xref ref-type="bibr" rid="CR45">45</xref>
</sup>
. Hansen and colleagues demonstrated that the grand average FC is more closely linked to the SC and linear models of FC<sup><xref ref-type="bibr" rid="CR29">29</xref>
</sup>
. When non-linearities are considered in the network models, the spatiotemporally dynamic repertoire of the network is significantly enhanced and the resting state dynamics shows the non-stationary FCD, which expresses itself as the switching dynamics of the FC. While Hansen and colleagues proposed FCD as a novel biomarker and demonstrated that all known RSNs can be derived from the non-linear network dynamics of FCD, they did not fit the model to the empirical functional time series data. The patterns in the FCD matrix arise from what is essentially a random process and thus different for different measurements. This renders the fitting process for brain network models more complex than fitting with the grand average FC, for which a Pearson correlation across empirical and simulated FC matrices is sufficient.</p>
<p>We have addressed this issue through a systematic fitting approach of the random process in FCD to the empirical data. The conjunction of using sophisticated fitting and systematic parameter analysis allowed us to test the mechanistic hypotheses underlying the resting state, i.e. whether the brain at rest operates close to the edge of a bifurcation and/or occupies a metastable state. Both scenarios can be mechanistically realized by non-linearly coupling Hopf bifurcators<sup><xref ref-type="bibr" rid="CR26">26</xref>
</sup>
. Hopf oscillators have been used previously in connectome-based modelling of resting state dynamics in EEG/MEG and fMRI<sup><xref ref-type="bibr" rid="CR14">14</xref>
</sup>
, as well as for the modelling of the detailed temporal dynamics in EEG/MEG<sup><xref ref-type="bibr" rid="CR27">27</xref>
, <xref ref-type="bibr" rid="CR28">28</xref>
</sup>
. The usage here though is different from the previous research, since the Hopf oscillators act as the sources of BOLD signal in the connectome based network model. Ghosh and colleagues used the Hopf oscillators as the sources of the electrophysiological signal and employed the Balloon Windkessel to derive the BOLD signal<sup><xref ref-type="bibr" rid="CR46">46</xref>
</sup>
. Given this interpretation, they needed to include all the signal transmission delays. In our present approach, the oscillation frequencies are significantly slower and thus permit the neglect of the time delays, which simplifies the computational effort of the simulation and thus the computational fitting of the models against empirical data.</p>
<p>Our key finding is the demonstration that the optimal operating regime is at the edge of the local Hopf bifurcation, i.e. a balance of noisy excursions in the oscillatory state. We not only were able to demonstrate that previous findings on the optimal operating point based on grand-average FC hold true if we take into account the temporal dynamics of FC, i.e. FCD. We also demonstrated that a better way of constraining brain network models is by not only fitting the grand average FC, but by also fitting the temporal structure of the fluctuations using the FCD.</p>
<p>Another remarkable and important finding is that high metastability is only present in a narrow range of bifurcation parameter when <italic>a</italic>
 is close to the edge of the bifurcation. In other words, the FCD of the spontaneous resting state, in conjunction with brain network modelling provide evidence that the brain at rest is maximally metastable, refining and demonstrating the hypothesis of Tognoli and Kelso<sup><xref ref-type="bibr" rid="CR6">6</xref>
</sup>
. Note that there is also a region for very small <italic>G</italic>
 and positive <italic>a</italic>
 (oscillatory regime) where a relatively good fitting is obtained. This dynamic regime was previously observed with a pure oscillatory Kuramoto model of the BOLD signals at the mesoscopic level<sup><xref ref-type="bibr" rid="CR47">47</xref>
</sup>
. Nevertheless, the level of fitting for the FC, metastability and even FCD is not as good as the one obtained in the region at the edge of the Hopf bifurcation. On the other hand, besides the extreme sensitivity of that working point (ultra-narrow regime of optimality) which means that the result is not so robust, the qualitative description of the BOLD signals is not realistic in the pure oscillatory regime in comparison with the noisy/oscillatory excursions evidenced in the regime of the bifurcation parameter <italic>a</italic>
 near zero.</p>
<p>For constructing a heterogeneous brain network model with different local parameter values, we took into account the spectral information of the BOLD data. We addressed the question if the oscillations at the individual nodes play a mechanistic role for the emergence of FC/FCD. In particular, we identified a cortical core of eight brain regions with the optimal fit of bifurcation parameter <italic>a</italic>
 close to the edge of bifurcation. We propose to call this the <italic>dynamical cortical core</italic>
 of the brain. Interestingly, three of these regions (the medial orbitofrontal cortex, posterior cingulate cortex and precuneus cortex) are part of the default mode network and thus re-experience past events and pre-experience possible future events<sup><xref ref-type="bibr" rid="CR48">48</xref>
, <xref ref-type="bibr" rid="CR49">49</xref>
</sup>
. In this vein other regions (parahippocampal and transverse temporal gyrus) have also been implicated in memory processing and may thus perhaps be helping integrate information over different timescales, binding fast and slow processes over time<sup><xref ref-type="bibr" rid="CR2">2</xref>
</sup>
. This information is always contextual and in the noisy, unpredictable scanner it is perhaps not surprising that the brain is attending to the auditory signals (transverse temporal gyrus). As such this information processing is available for conflict monitoring and selection for action (rostral anterior cingulate cortex and caudal middle frontal gyrus) and motor execution (precentral gyrus)<sup><xref ref-type="bibr" rid="CR50">50</xref>
</sup>
. Equally, the involvement of the cingulate cortex is interesting given that this region recently has been shown to be part of the common neurobiological substrate for mental illness across across six diverse diagnostic groups (schizophrenia, bipolar disorder, depression, addiction, obsessive-compulsive disorder, and anxiety) based on a meta-analysis of grey matter loss in 193 neuroimaging studies of 15892 individuals<sup><xref ref-type="bibr" rid="CR51">51</xref>
</sup>
. This reinforces the potential use of brain network computational modelling for understanding the underlying mechanisms of neuropsychiatric disorders<sup><xref ref-type="bibr" rid="CR52">52</xref>
</sup>
. The right-handed quality of Fig. <xref rid="Fig6" ref-type="fig">6</xref>
, presumably arises from the specifics of the data used to fit the model and we will be exploring its biological validity in subsequent studies with larger group sizes.</p>
<p>Although the bifurcation parameter does not have a direct biophysical correlate, it seems to be involved in mediating biophysical effects. Another note of caution: we have presented a mesoscopic phenomenological model, i.e. the dynamical equation corresponds directly to the measured BOLD signals and not to the neural signals. It is a <italic>phenomenological</italic>
 model since the real coupling between regions does not, of course, occur between the hemodynamic signals, but between the underlying neural activity. Still, by using this simplified approach, we do not need to convolve modelled neural signals with the haemodynamic response function<sup><xref ref-type="bibr" rid="CR47">47</xref>
</sup>
.</p>
<p>We therefore propose that in future both the global bifurcation parameter as well as the individual parameters could potentially serve as biomarkers for disease. It will be important to explore the changes for different brain diseases, e.g. within a standardized framework for connectome-based modelling such as The Virtual Brain (TVB)<sup><xref ref-type="bibr" rid="CR53">53</xref>
, <xref ref-type="bibr" rid="CR54">54</xref>
</sup>
, and applications such as fitting of TVB’s dynamic regime and TVB Processing pipeline<sup><xref ref-type="bibr" rid="CR17">17</xref>
</sup>
.</p>
<p>Overall, we have shown that neuroimaging data can be causally analysed by constructing a relatively simple brain network computational model using a Hopf bifurcation. This model was shown to be maximally metastable at the optimal fitting with the spatiotemporal dynamics of spontaneous brain activity. This dynamical regime may well allow for the optimal integration and segregation of fast and slow information over different time-scales, the “flights” and “perchings” of the stream of consciousness alluded to by William James over 100 years ago.</p>
</sec>
<sec id="Sec6" sec-type="materials|methods"><title>Methods</title>
<sec id="Sec7"><title>Ethics Statement</title>
<p>All participants of this study gave written informed consent before the study, which was performed in compliance with the relevant laws and institutional guidelines and approved by the ethics committee of the Charité University Berlin.</p>
</sec>
<sec id="Sec8"><title>Empirical MRI Data Collection</title>
<p>Structural data from DTI and resting-state BOLD signal time series were acquired for 24 healthy participants (age between 18 and 33 years old, mean 25.7, 12 females, 12 males). A full description of the generation of SC and FC matrices from those data can be found in ref. <xref ref-type="bibr" rid="CR17">17</xref>
. Here, we provide a quick overview of the employed methods. Empirical data were acquired at Berlin Center for Advanced Imaging, Charité University Medicine, Berlin, Germany. For simultaneous EEG-fMRI<sup><xref ref-type="bibr" rid="CR55">55</xref>
, <xref ref-type="bibr" rid="CR56">56</xref>
</sup>
, participants were asked to stay awake and keep their eyes closed. No other controlled task had to be performed. In addition, a localizer, DTI and T2 sequence were recorded for each participant. MRI was performed using a 3 Tesla Siemens Trim Trio MR scanner and a 12-channel Siemens head coil. Specifications for the employed sequences can be found in ref. <xref ref-type="bibr" rid="CR56">56</xref>
. For each participant anatomical T1-weighted scans were acquired. DTI and GRE field mapping were measured directly after the anatomical scans. Next, functional MRI (BOLD-sensitive, T2*-weighted, TR 1940 ms, TE 30 ms, FA 78°, 32 transversal slices (3 mm), voxel size 3 × 3 × 3 mm, FoV 192 mm, 64 matrix) was recorded simultaneously to the EEG recording.</p>
</sec>
<sec id="Sec9"><title>MRI Data Analysis</title>
<p>Processing steps executed by the public Berlin automatized processing pipeline<sup><xref ref-type="bibr" rid="CR56">56</xref>
</sup>
 comprised 1) preprocessing of T1-weighted scans, cortical reconstruction, tessellation and parcellation, 2) transformation of anatomical masks to diffusion space, 3) processing of diffusion data, 4) transformation of anatomical masks to fMRI space, 5) Processing of fMRI data.</p>
</sec>
<sec id="Sec10"><title>Anatomical MRI Data Analysis</title>
<p>The highly resolved anatomical images are important to create a precise parcellation of the brain. For each of those parcellated units, empirical functional data time series are spatially aggregated. T1-weighted images are pre-processed using FREESURFER including probabilistic atlas based cortical parcellation, here using Desikan-Killany (DK) atlas<sup><xref ref-type="bibr" rid="CR57">57</xref>
</sup>
 (Table <xref rid="Tab1" ref-type="table">1</xref>
). This generates volumes that contain all cortical and subcortical parcellated regions with corresponding region labels used for fiber-tracking and BOLD time-series extraction.<table-wrap id="Tab1"><label>Table 1</label>
<caption><p>Anatomical labels for the 68 regions in the Desikan-Kahilly parcellation.</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th>Region number</th>
<th>Region name</th>
</tr>
</thead>
<tbody><tr><td>1;35</td>
<td>Superior temporal sulcus, banks of</td>
</tr>
<tr><td>2;36</td>
<td>Caudal anterior cingulate cortex</td>
</tr>
<tr><td>3;37</td>
<td>Caudal middle frontal gyrus</td>
</tr>
<tr><td>4;38</td>
<td>Cuneus cortex</td>
</tr>
<tr><td>5;39</td>
<td>Entorhinal cortex</td>
</tr>
<tr><td>6;40</td>
<td>Fusiform gyrus</td>
</tr>
<tr><td>7;41</td>
<td>Inferior parietal cortex</td>
</tr>
<tr><td>8;42</td>
<td>Inferior temporal gyrus</td>
</tr>
<tr><td>9;43</td>
<td>Isthmus of cingulate cortex</td>
</tr>
<tr><td>10;44</td>
<td>Lateral occipital cortex</td>
</tr>
<tr><td>11;45</td>
<td>Lateral orbitofrontal cortex</td>
</tr>
<tr><td>12;46</td>
<td>Lingual gyrus</td>
</tr>
<tr><td>13;47</td>
<td>Medial orbitofrontal cortex</td>
</tr>
<tr><td>14;48</td>
<td>Middle temporal gyrus</td>
</tr>
<tr><td>15;49</td>
<td>Parahippocampal gyrus</td>
</tr>
<tr><td>16;50</td>
<td>Paracentral lobule</td>
</tr>
<tr><td>17;51</td>
<td>Pars opercularis</td>
</tr>
<tr><td>18;52</td>
<td>Pars orbitalis</td>
</tr>
<tr><td>19;53</td>
<td>Pars triangularis</td>
</tr>
<tr><td>20;54</td>
<td>Pericalcarine cortex</td>
</tr>
<tr><td>21;55</td>
<td>Postcentral gyrus</td>
</tr>
<tr><td>22;56</td>
<td>Posterior cingulate cortex</td>
</tr>
<tr><td>23;57</td>
<td>Precentral gyrus</td>
</tr>
<tr><td>24;58</td>
<td>Precuneus cortex</td>
</tr>
<tr><td>25;59</td>
<td>Rostral anterior cingulate cortex</td>
</tr>
<tr><td>26;60</td>
<td>Rostral middle frontal gyrus</td>
</tr>
<tr><td>27;61</td>
<td>Superior frontal cortex</td>
</tr>
<tr><td>28;62</td>
<td>Superior parietal cortex</td>
</tr>
<tr><td>29;63</td>
<td>Superior temporal gyrus</td>
</tr>
<tr><td>30;64</td>
<td>Supramarginal gyrus</td>
</tr>
<tr><td>31;65</td>
<td>Frontal pole</td>
</tr>
<tr><td>32;66</td>
<td>Temporal pole</td>
</tr>
<tr><td>33;67</td>
<td>Transverse temporal cortex (primary auditory cortex)</td>
</tr>
<tr><td>34;68</td>
<td>Insula</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p>The two region numbers per line refer to right and left hemisphere respectively.</p>
</table-wrap-foot>
</table-wrap>
</p>
</sec>
<sec id="Sec11"><title>Empirical DTI Data Analysis and Tractography</title>
<p>Tractography requires binary WM masks to restrict tracking to WM voxels. Upon extraction of gradient vectors and values (known as b-table) using MRTrix, dw-MRI data are pre-processed using FREESURFER. Besides motion correction and eddy current correction (ECC) the b0 image is linearly registered (6 degrees of freedom, DOF) to the participant’s anatomical T1-weighted image and the resulting registration rule is stored for later use. We transformed the high-resolution mask volumes from the anatomical space to the participant’s diffusion space, to further use it for fiber tracking. The cortical and subcortical parcellations are resampled into diffusion space, one time using the original 1 mm isotropic voxel size (for subvoxel seeding) and one time matching that of our dw-MRI data, i.e., 2.3 mm isotropic voxel size. During MRTrix pre-processing diffusion tensor images that store the diffusion tensor (i.e., the diffusion ellipsoid) for each voxel location are computed. Based on that, a fractional anisotropy (FA) and an eigenvector map are computed and masked by the binary WM mask created previously. For subsequent fiber-response function estimation, a mask containing high-anisotropy voxels is computed. Fibre orientation distributions are estimated using constrained spherical deconvolution<sup><xref ref-type="bibr" rid="CR58">58</xref>
</sup>
 based on a response function estimated in voxels that are expected to contain a single, coherently-oriented fibre bundle (commands dwi2response tournier and dwi2fod; see MRTrix Documentation: <ext-link ext-link-type="uri" xlink:href="http://mrtrix.readthedocs.io/en/latest/">http://mrtrix.readthedocs.io/en/latest/</ext-link>
). In order to resolve crossing pathways, fibers are prolonged by employing a probabilistic tracking approach as provided by MRTrix. In order to exclude spurious tracks, three types of masks are used to constrain tracking: seeding-, target- and stop-masks. In order to restrict track-prolongation to WM, a WM-mask that contains the union of GM-WM-interface and cortical WM voxels is defined as a global stop mask for tracking. To address several confounds in the estimation of connection strengths (information transmission capacities), a new seeding and fiber aggregation strategy was employed developed for this pipeline and described in detail in ref. <xref ref-type="bibr" rid="CR17">17</xref>
. In combination with a new aggregation scheme, it is based on an appropriate selection of seed voxels and controlling for the number of generated tracks in each seed voxel. Instead of using every WM voxel, tracks are initiated from GM-WM-interface voxels and a fixed number of tracks are generated for each seed-voxel. Since a GM parcellation-based aggregation is performed, each seed-mask is associated with a ROI of the GM atlas. Along with seeding-masks complementary target-masks are defined specifying valid terminal regions for each track that was initiated in a specific seed voxel. The capacity measures that we derive between each pair of regions are intended to estimate the strength of the influence that one region exerts over another, i.e., their SC. In order to improve existing methods for capacities estimation the approach makes use of several assumptions with regard to seed-ROI selection, tracking and aggregation of generated tracks<sup><xref ref-type="bibr" rid="CR17">17</xref>
</sup>
. Upon tractography the pipeline aggregates generated tracks to structural connectome matrices. The weighted distinct connection count used here divides each distinct connection by the number of distinct connections leaving the seed-voxel (yielding asymmetric capacities matrix). Values have been normalized by the total surface area of the GWI of a participant.</p>
</sec>
<sec id="Sec12"><title>Empirical fMRI Data Analysis</title>
<p>In order to generate the FC matrices, FSL’s FEAT pipeline is used to perform the following operations: deleting the first five images of the series to exclude possible saturation effects in the images, high-pass temporal filtering (100 seconds high-pass filter), motion correction, brain extraction and a 6 DOF linear registration to the MNI space. Functional data is registered to the participant’s T1-weighted images and parcellated according to FREESURFER’s cortical segmentation. By inverting the mapping rule found by registration, anatomical segmentations are mapped onto the functional space. Finally, average BOLD signal time series for each region are generated by computing the mean over all voxel time-series for each region. From the region wise aggregated BOLD data, FC matrices are computed within MATLAB using and Pearson’s linear correlation coefficient as FC metrics. We did not perform global signal regression on data.</p>
</sec>
<sec id="Sec13"><title>Brain Network Model</title>
<p>The brain network model consists of 68 coupled brain areas (nodes) derived from the parcellation explained above. The global dynamics of the brain network model used here results from the mutual interactions of local node dynamics coupled through the underlying empirical anatomical structural connectivity matrix <italic>C</italic>
<sub><italic>ij</italic>
</sub>
 (see Fig. <xref rid="Fig2" ref-type="fig">2</xref>
). The structural matrix <italic>C</italic>
<sub><italic>i</italic>
</sub>
 denotes the density of fibres between cortical area i and j as extracted from the DTI based tractography (scaled to a maximum value of 0.2). The local dynamics of each individual node is described by the normal form of a supercritical Hopf bifurcation, which is able to describe the transition from asynchronous noisy behavior to full oscillations. Thus, in complex coordinates, each node <italic>j</italic>
 is described by following equation:<disp-formula id="Equ1"><label>1</label>
<alternatives><tex-math id="M1">\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym} 
				\usepackage{amsfonts} 
				\usepackage{amssymb} 
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\frac{d{z}_{j}}{dt}=z[{a}_{j}+i{\omega }_{j}-|{z}_{j}^{2}|]+\beta {\eta }_{j}(t)$$\end{document}</tex-math>
<mml:math id="M2" display="block"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi>
<mml:msub><mml:mrow><mml:mi>z</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow><mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mi>z</mml:mi>
<mml:mo stretchy="false">[</mml:mo>
<mml:msub><mml:mrow><mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mi>i</mml:mi>
<mml:msub><mml:mrow><mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>−</mml:mo>
<mml:mo stretchy="false">|</mml:mo>
<mml:msubsup><mml:mrow><mml:mi>z</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo stretchy="false">|</mml:mo>
<mml:mo stretchy="false">]</mml:mo>
<mml:mo>+</mml:mo>
<mml:mi>β</mml:mi>
<mml:msub><mml:mrow><mml:mi>η</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
<graphic xlink:href="41598_2017_3073_Article_Equ1.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
where<disp-formula id="Equ2"><label>2</label>
<alternatives><tex-math id="M3">\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym} 
				\usepackage{amsfonts} 
				\usepackage{amssymb} 
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${z}_{j}={\rho }_{j}{e}^{i{\theta }_{j}}={x}_{j}+i{y}_{j}$$\end{document}</tex-math>
<mml:math id="M4" display="block"><mml:msub><mml:mrow><mml:mi>z</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub><mml:mrow><mml:mi>ρ</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msup><mml:mrow><mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:msub><mml:mrow><mml:mi>θ</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msup>
<mml:mo>=</mml:mo>
<mml:msub><mml:mrow><mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mi>i</mml:mi>
<mml:msub><mml:mrow><mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
<graphic xlink:href="41598_2017_3073_Article_Equ2.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
and <italic>η</italic>
<sub><italic>i</italic>
</sub>
(<italic>t</italic>
) is additive Gaussian noise with standard deviation β=0.02. This normal form has a supercritical bifurcation at <italic>a</italic>
<sub><italic>j</italic>
</sub>
 = 0, so that for <italic>a</italic>
<sub><italic>j</italic>
</sub>
 < 0 the local dynamics has a stable fixed point at <italic>z</italic>
<sub><italic>j</italic>
</sub>
 = 0 (which because of the additive noise corresponds to a low activity asynchronous state) and for <italic>a</italic>
<sub><italic>j</italic>
</sub>
 > 0 there exists a stable limit cycle oscillation with frequency <inline-formula id="IEq1"><alternatives><tex-math id="M5">\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym} 
				\usepackage{amsfonts} 
				\usepackage{amssymb} 
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${f}_{j}={\omega }_{j}/2\pi $$\end{document}</tex-math>
<mml:math id="M6"><mml:msub><mml:mrow><mml:mi>f</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub><mml:mrow><mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>π</mml:mi>
</mml:math>
<inline-graphic xlink:href="41598_2017_3073_Article_IEq1.gif"></inline-graphic>
</alternatives>
</inline-formula>
. We insert equation <xref rid="Equ2" ref-type="">2</xref>
 in equation <xref rid="Equ1" ref-type="">1</xref>
 and separate real part in equation <xref rid="Equ3" ref-type="">3</xref>
 and imaginary part in equation <xref rid="Equ4" ref-type="">4</xref>
.</p>
<p>Thus, the whole-brain dynamics is defined by following set of coupled equations:<disp-formula id="Equ3"><label>3</label>
<alternatives><tex-math id="M7">\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym} 
				\usepackage{amsfonts} 
				\usepackage{amssymb} 
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\frac{d{x}_{j}}{dt}=[{a}_{j}-{x}_{j}^{2}-{y}_{j}^{2}]{x}_{j}-{\omega }_{j}{y}_{j}+G\,\sum _{i}{C}_{ij}({x}_{i}-{x}_{j})\beta {\eta }_{j}(t)$$\end{document}</tex-math>
<mml:math id="M8" display="block"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi>
<mml:msub><mml:mrow><mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow><mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mo stretchy="false">[</mml:mo>
<mml:msub><mml:mrow><mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>−</mml:mo>
<mml:msubsup><mml:mrow><mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>−</mml:mo>
<mml:msubsup><mml:mrow><mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo stretchy="false">]</mml:mo>
<mml:msub><mml:mrow><mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>−</mml:mo>
<mml:msub><mml:mrow><mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mrow><mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mi>G</mml:mi>
<mml:mspace width=".25em"></mml:mspace>
<mml:munder><mml:mo>∑</mml:mo>
<mml:mi>i</mml:mi>
</mml:munder>
<mml:msub><mml:mrow><mml:mi>C</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:mo stretchy="true">(</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>−</mml:mo>
<mml:msub><mml:mrow><mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="true">)</mml:mo>
</mml:mrow>
<mml:mi>β</mml:mi>
<mml:msub><mml:mrow><mml:mi>η</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:mo stretchy="true">(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="true">)</mml:mo>
</mml:mrow>
</mml:math>
<graphic xlink:href="41598_2017_3073_Article_Equ3.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ4"><label>4</label>
<alternatives><tex-math id="M9">\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym} 
				\usepackage{amsfonts} 
				\usepackage{amssymb} 
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$\frac{d{y}_{j}}{dt}=[{a}_{j}-{x}_{j}^{2}-{y}_{j}^{2}]{y}_{j}+{\omega }_{j}{x}_{j}+G\,\sum _{i}{C}_{ij}({y}_{i}-{y}_{j})\beta {\eta }_{j}(t)$$\end{document}</tex-math>
<mml:math id="M10" display="block"><mml:mfrac><mml:mrow><mml:mi>d</mml:mi>
<mml:msub><mml:mrow><mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow><mml:mi>d</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mo stretchy="false">[</mml:mo>
<mml:msub><mml:mrow><mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>−</mml:mo>
<mml:msubsup><mml:mrow><mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>−</mml:mo>
<mml:msubsup><mml:mrow><mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo stretchy="false">]</mml:mo>
<mml:msub><mml:mrow><mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub><mml:mrow><mml:mi>ω</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub><mml:mrow><mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mi>G</mml:mi>
<mml:mspace width=".25em"></mml:mspace>
<mml:munder><mml:mo>∑</mml:mo>
<mml:mi>i</mml:mi>
</mml:munder>
<mml:msub><mml:mrow><mml:mi>C</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:mo stretchy="false">(</mml:mo>
<mml:mrow><mml:msub><mml:mrow><mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>−</mml:mo>
<mml:msub><mml:mrow><mml:mi>y</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>β</mml:mi>
<mml:msub><mml:mrow><mml:mi>η</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:mo stretchy="true">(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="true">)</mml:mo>
</mml:mrow>
</mml:math>
<graphic xlink:href="41598_2017_3073_Article_Equ4.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
</p>
<p>Please note that following the literature from physics, the equations are written in Cartesian, rather than polar coordinates<sup><xref ref-type="bibr" rid="CR59">59</xref>
–<xref ref-type="bibr" rid="CR62">62</xref>
</sup>
. We couple the equations using the common difference coupling, which approximates the simplest (linear) part of a general coupling function. These equations are valid in the weakly coupled oscillator limit, in which the coupling preserves the periodic orbit of the uncoupled oscillators. If the linear coupling (following a Taylor expansion of the full coupling) does not exist, the next non-vanishing higher order term should be considered, which is a case we do not address here (please see Kuramoto<sup><xref ref-type="bibr" rid="CR59">59</xref>
</sup>
 (see Eq 5.3.1) and Pikovsky, Arkady and Kurths<sup><xref ref-type="bibr" rid="CR60">60</xref>
</sup>
 (see Eq. 8.12) for more detailed analytic treatments of the equations).</p>
<p>In the latter equations, G is a global scaling factor (global conductivity parameter scaling equally all synaptic connections). The global scaling factor G and the bifurcation parameters <italic>a</italic>
<sub><italic>j</italic>
</sub>
 are the control parameters with which we study the optimal dynamical working region where the simulations maximally fit the empirical FC and the FCD. We model with the variables <italic>x</italic>
<sub><italic>j</italic>
</sub>
 the BOLD signal of each node <italic>j</italic>
. The empirical BOLD signals were band-pass filtered within the narrowband 0.04–0.07 Hz. This frequency band has been mapped to the gray matter and it has been shown to be more reliable and functionally relevant than other frequency bands<sup><xref ref-type="bibr" rid="CR63">63</xref>
–<xref ref-type="bibr" rid="CR66">66</xref>
</sup>
. Within this model, the intrinsic frequency <italic>ω</italic>
<sub><italic>j</italic>
</sub>
 of each node is in the 0.04–0.07 Hz band (<italic>i</italic>
 = 1, …, <italic>n</italic>
). The intrinsic frequencies were estimated from the data, as given by the averaged peak frequency of the narrowband BOLD signals of each brain region.</p>
</sec>
<sec id="Sec14"><title>Grand average FC and FCD matrices</title>
<p>The grand average FC is defined as the matrix of correlations of the BOLD signals between two brain areas over the whole time window of acquisition. In order to characterize the time dependent structure of the resting fluctuations, we estimate the FCD matrix<sup><xref ref-type="bibr" rid="CR29">29</xref>
</sup>
 (see Fig. <xref rid="Fig1" ref-type="fig">1</xref>
). Each full-length BOLD signal of 22 min is split up into M=61 sliding windows of 60 sec, overlapping by 40 sec. For each sliding window, centered at time t, we calculated a separate FC matrix, FC(t). The FCD is a MxM symmetric matrix whose (t1, t2) entry is defined by the Pearson correlation between the upper triangular parts of the two matrices FC(t1) and FC(t2). Epochs of stable FC(t) configurations are reflected around the FCD diagonal in blocks of elevated inter-FC(t) correlations.</p>
<p>The grand average FC and the FCD matrices were estimated for the recordings of each of the 24 participants as well as for 24 simulations of 22 minutes of the computational model. We compared the FC matrices of the model (averaged Fisher’s z-transformed over the 24 sessions) and the empirical data (averaged Fisher’s z-transformed over the 24 participants), adopting as a measure of similarity between the two matrices the Pearson correlation coefficient between corresponding elements of the upper triangular part of the matrices. For comparing the FCD statistics, we collected the upper triangular elements of the matrices (over all participants or sessions) and generated the distribution of them. Then, we compared the simulated and empirical distribution by means of the Kolmogorov-Smirnov distance between them. The Kolmogorov–Smirnov distance quantifies the maximal difference between the cumulative distribution functions of the two samples.</p>
</sec>
<sec id="Sec15"><title>Metastability</title>
<p>Here, we refer to metastability as a measure of how variable are the states of phase configurations as a function of time, i.e. how the synchronization between the different nodes fluctuates across time<sup><xref ref-type="bibr" rid="CR24">24</xref>
</sup>
. Thus, we measure the metastability as the standard deviation of the Kuramoto order parameter across time. The Kuramoto order parameter is defined by following equation:<disp-formula id="Equ5"><label>5</label>
<alternatives><tex-math id="M11">\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym} 
				\usepackage{amsfonts} 
				\usepackage{amssymb} 
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$R(t)=|\sum _{k=1}^{n}{e}^{i{\phi }_{k}(t)}|/n$$\end{document}</tex-math>
<mml:math id="M12" display="block"><mml:mi>R</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>=</mml:mo>
<mml:mrow><mml:mo stretchy="true">|</mml:mo>
<mml:mrow><mml:munderover><mml:mo>∑</mml:mo>
<mml:mrow><mml:mi>k</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:msup><mml:mrow><mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>i</mml:mi>
<mml:msub><mml:mrow><mml:mi>ϕ</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="true">|</mml:mo>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mi>n</mml:mi>
</mml:math>
<graphic xlink:href="41598_2017_3073_Article_Equ5.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
where <italic>φ</italic>
<sub><italic>k</italic>
</sub>
(<italic>t</italic>
) is the instantaneous phase of each narrowband BOLD signal at node <italic>k</italic>
. The Kuramoto order parameter measures the global level of synchronization of the <italic>n</italic>
 oscillating signals. Under complete independence, the <italic>n</italic>
 phases are uniformly distributed and thus <italic>R</italic>
 is nearly zero, whereas <italic>R</italic>
=1 if all phases are equal (full synchronization). Thus, for calculating the metastability of the empirical and simulated BOLD signals, we first band-pass filtered within the narrowband 0.04–0.07 Hz (as previously explained) and computed the instantaneous phase <italic>φ</italic>
<sub><italic>k</italic>
</sub>
(<italic>t</italic>
) of each narrowband signal <italic>k</italic>
 using the Hilbert transform. The Hilbert transform yields the associated analytical signals. The analytic signal represents a narrowband signal, <italic>s</italic>
(<italic>t</italic>
), in the time domain as a rotating vector with an instantaneous phase, <italic>φ</italic>
<sub><italic>k</italic>
</sub>
(<italic>t</italic>
), and an instantaneous amplitude, <italic>A</italic>
(<italic>t</italic>
), i.e., <inline-formula id="IEq2"><alternatives><tex-math id="M13">\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym} 
				\usepackage{amsfonts} 
				\usepackage{amssymb} 
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$s(t)=A(t)\cos (\phi (t))$$\end{document}</tex-math>
<mml:math id="M14"><mml:mi>s</mml:mi>
<mml:mrow><mml:mo stretchy="true">(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="true">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>A</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mrow><mml:mo stretchy="true">(</mml:mo>
<mml:mrow><mml:mi>ϕ</mml:mi>
<mml:mrow><mml:mo stretchy="true">(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="true">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="true">)</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="41598_2017_3073_Article_IEq2.gif"></inline-graphic>
</alternatives>
</inline-formula>
. The phase and the amplitude are given by the argument and the modulus, respectively, of the complex signal <italic>z</italic>
(<italic>t</italic>
), given by <inline-formula id="IEq3"><alternatives><tex-math id="M15">\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym} 
				\usepackage{amsfonts} 
				\usepackage{amssymb} 
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$$z(t)=s(t)+i.H[s(t)]$$\end{document}</tex-math>
<mml:math id="M16"><mml:mi>z</mml:mi>
<mml:mrow><mml:mo stretchy="true">(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="true">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>s</mml:mi>
<mml:mrow><mml:mo stretchy="true">(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="true">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mi>i</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>H</mml:mi>
<mml:mrow><mml:mo stretchy="true">[</mml:mo>
<mml:mrow><mml:mi>s</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo stretchy="true">]</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="41598_2017_3073_Article_IEq3.gif"></inline-graphic>
</alternatives>
</inline-formula>
, where <italic>i</italic>
 is the imaginary unit and H[<italic>s</italic>
(<italic>t</italic>
)] is the Hilbert transform of <italic>s</italic>
(t).</p>
</sec>
<sec id="Sec16"><title>Local Optimization of Brain Nodes</title>
<p>The local optimization of each single bifurcation parameter <italic>a</italic>
<sub><italic>j</italic>
</sub>
 is based on the fitting of the spectral information of the empirical BOLD signals in each node. In particular, we aim to fit the proportion of power in the 0.04–0.07 Hz band with respect to the 0.04–0.25 Hz band (i.e. we remove the smallest frequencies below 0.04 Hz and consider the whole spectra until the Nyquist frequency which is 0.25 Hz)<sup><xref ref-type="bibr" rid="CR47">47</xref>
</sup>
. For this, we filtered the BOLD signals in the 0.04–0.25 Hz band, and calculated the power spectrum <italic>P</italic>
<sub><italic>j</italic>
</sub>
(<italic>f</italic>
) for each node <italic>j</italic>
. We define the proportion,<disp-formula id="Equ6"><label>6</label>
<alternatives><tex-math id="M17">\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym} 
				\usepackage{amsfonts} 
				\usepackage{amssymb} 
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${p}_{j}=\frac{{\int }_{0.04}^{0.07}{P}_{j}(f){\rm{d}}f}{{\int }_{0.04}^{0.25}{P}_{j}(f){\rm{d}}f}$$\end{document}</tex-math>
<mml:math id="M18" display="block"><mml:msub><mml:mrow><mml:mi>p</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac><mml:mrow><mml:msubsup><mml:mo>∫</mml:mo>
<mml:mrow><mml:mn>0.04</mml:mn>
</mml:mrow>
<mml:mrow><mml:mn>0.07</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:msub><mml:mrow><mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:mo stretchy="true">(</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo stretchy="true">)</mml:mo>
</mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
<mml:mrow><mml:msubsup><mml:mo>∫</mml:mo>
<mml:mrow><mml:mn>0.04</mml:mn>
</mml:mrow>
<mml:mrow><mml:mn>0.25</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:msub><mml:mrow><mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:mo stretchy="true">(</mml:mo>
<mml:mi>f</mml:mi>
<mml:mo stretchy="true">)</mml:mo>
</mml:mrow>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi>f</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:math>
<graphic xlink:href="41598_2017_3073_Article_Equ6.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
and update the local bifurcation parameters by a gradient descendent strategy, i.e.:<disp-formula id="Equ7"><label>7</label>
<alternatives><tex-math id="M19">\documentclass[12pt]{minimal}
				\usepackage{amsmath}
				\usepackage{wasysym} 
				\usepackage{amsfonts} 
				\usepackage{amssymb} 
				\usepackage{amsbsy}
				\usepackage{mathrsfs}
				\usepackage{upgreek}
				\setlength{\oddsidemargin}{-69pt}
				\begin{document}$${a}_{j}={a}_{j}+\eta ({p}_{j}^{empirical}-{p}_{j}^{simulated})$$\end{document}</tex-math>
<mml:math id="M20" display="block"><mml:msub><mml:mrow><mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub><mml:mrow><mml:mi>a</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mi>η</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:msubsup><mml:mrow><mml:mi>p</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>e</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>l</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>−</mml:mo>
<mml:msubsup><mml:mrow><mml:mi>p</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>u</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo stretchy="false">)</mml:mo>
</mml:math>
<graphic xlink:href="41598_2017_3073_Article_Equ7.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
until convergence. We used here <italic>η</italic>
 = 0.1. The updates of the <italic>a</italic>
<sub><italic>j</italic>
</sub>
 values are done in each optimization step in parallel.</p>
</sec>
</sec>
</body>
<back><fn-group><fn><p><bold>Publisher's note:</bold>
 Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p>
</fn>
</fn-group>
<ack><title>Acknowledgements</title>
<p>GD is supported by the ERC Advanced Grant DYSTRUCTURE (n. 295129), by the Spanish Research Project PSI2016-75688-P. MLK is supported by the ERC Consolidator Grant: CAREGIVING (n. 615539) and Center for Music in the Brain, funded by the Danish National Research Foundation (DNRF117). VJ and GD are supported by the European Union’s Horizon 2020 research and innovation programme under grant agreement n. 720270 (HBP SGA1). VJ and PR are supported by the James S. McDonnell Foundation (Brain Network Recovery Group JSMF22002082). VJ is supported by FHU EPINEXT [A<sup>*</sup>
MIDEX project (ANR-11-IDEX-0001-02) funded by the ‘Investissements d’Avenir’ French Government]. PR is supported the German Ministry of Education and Research (US-German Collaboration in Computational Neuroscience 100258846 and Bernstein Focus State Dependencies of Learning 01GQ0971-5), the Max-Planck Society and funding from the European Union Horizon 2020 (ERC Consolidator grant BrainModes 683049).</p>
</ack>
<notes notes-type="author-contribution"><title>Author Contributions</title>
<p>G.D., M.K., V.J. & P.R. designed research, contributed data, carried out the analyses, prepared figures, wrote the manuscript.</p>
</notes>
<notes notes-type="COI-statement"><sec id="FPar1"><title>Competing Interests</title>
<p>The authors declare that they have no competing interests.</p>
</sec>
</notes>
<ref-list id="Bib1"><title>References</title>
<ref id="CR1"><label>1.</label>
<mixed-citation publication-type="other">James, W. <italic>The Principles of Psychology</italic>
. (Henry Holt, 1890).</mixed-citation>
</ref>
<ref id="CR2"><label>2.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Deco</surname>
<given-names>G</given-names>
</name>
<name><surname>Tononi</surname>
<given-names>G</given-names>
</name>
<name><surname>Boly</surname>
<given-names>M</given-names>
</name>
<name><surname>Kringelbach</surname>
<given-names>ML</given-names>
</name>
</person-group>
<article-title>Rethinking segregation and integration: contributions of whole-brain modelling</article-title>
<source>Nat Rev Neurosci</source>
<year>2015</year>
<volume>16</volume>
<fpage>430</fpage>
<lpage>439</lpage>
<pub-id pub-id-type="doi">10.1038/nrn3963</pub-id>
<pub-id pub-id-type="pmid">26081790</pub-id>
</element-citation>
</ref>
<ref id="CR3"><label>3.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Berridge</surname>
<given-names>KC</given-names>
</name>
<name><surname>Kringelbach</surname>
<given-names>ML</given-names>
</name>
</person-group>
<article-title>Pleasure systems in the brain</article-title>
<source>Neuron</source>
<year>2015</year>
<volume>86</volume>
<fpage>646</fpage>
<lpage>664</lpage>
<pub-id pub-id-type="doi">10.1016/j.neuron.2015.02.018</pub-id>
<pub-id pub-id-type="pmid">25950633</pub-id>
</element-citation>
</ref>
<ref id="CR4"><label>4.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Friston</surname>
<given-names>K</given-names>
</name>
<name><surname>Kiebel</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Predictive coding under the free-energy principle</article-title>
<source>Philosophical transactions of the Royal Society of London. Series B, Biological sciences</source>
<year>2009</year>
<volume>364</volume>
<fpage>1211</fpage>
<lpage>1221</lpage>
<pub-id pub-id-type="doi">10.1098/rstb.2008.0300</pub-id>
<pub-id pub-id-type="pmid">19528002</pub-id>
</element-citation>
</ref>
<ref id="CR5"><label>5.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kringelbach</surname>
<given-names>ML</given-names>
</name>
<name><surname>McIntosh</surname>
<given-names>AR</given-names>
</name>
<name><surname>Ritter</surname>
<given-names>P</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
<name><surname>Deco</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>The rediscovery of slowness: exploring the timing of cognition</article-title>
<source>TICS</source>
<year>2015</year>
<volume>19</volume>
<fpage>616</fpage>
<lpage>628</lpage>
</element-citation>
</ref>
<ref id="CR6"><label>6.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tognoli</surname>
<given-names>E</given-names>
</name>
<name><surname>Kelso</surname>
<given-names>JA</given-names>
</name>
</person-group>
<article-title>The metastable brain</article-title>
<source>Neuron</source>
<year>2014</year>
<volume>81</volume>
<fpage>35</fpage>
<lpage>48</lpage>
<pub-id pub-id-type="doi">10.1016/j.neuron.2013.12.022</pub-id>
<pub-id pub-id-type="pmid">24411730</pub-id>
</element-citation>
</ref>
<ref id="CR7"><label>7.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Friston</surname>
<given-names>KJ</given-names>
</name>
</person-group>
<article-title>The labile brain. I. Neuronal transients and nonlinear coupling</article-title>
<source>Philosophical transactions of the Royal Society of London. Series B, Biological sciences</source>
<year>2000</year>
<volume>355</volume>
<fpage>215</fpage>
<lpage>236</lpage>
<pub-id pub-id-type="doi">10.1098/rstb.2000.0560</pub-id>
<pub-id pub-id-type="pmid">10724457</pub-id>
</element-citation>
</ref>
<ref id="CR8"><label>8.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chaudhuri</surname>
<given-names>R</given-names>
</name>
<name><surname>Knoblauch</surname>
<given-names>K</given-names>
</name>
<name><surname>Gariel</surname>
<given-names>MA</given-names>
</name>
<name><surname>Kennedy</surname>
<given-names>H</given-names>
</name>
<name><surname>Wang</surname>
<given-names>XJ</given-names>
</name>
</person-group>
<article-title>A Large-Scale Circuit Mechanism for Hierarchical Dynamical Processing in the Primate Cortex</article-title>
<source>Neuron</source>
<year>2015</year>
<volume>88</volume>
<fpage>419</fpage>
<lpage>431</lpage>
<pub-id pub-id-type="doi">10.1016/j.neuron.2015.09.008</pub-id>
<pub-id pub-id-type="pmid">26439530</pub-id>
</element-citation>
</ref>
<ref id="CR9"><label>9.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Simony</surname>
<given-names>E</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Dynamic reconfiguration of the default mode network during narrative comprehension</article-title>
<source>Nat Commun</source>
<year>2016</year>
<volume>7</volume>
<fpage>12141</fpage>
<pub-id pub-id-type="doi">10.1038/ncomms12141</pub-id>
<pub-id pub-id-type="pmid">27424918</pub-id>
</element-citation>
</ref>
<ref id="CR10"><label>10.</label>
<mixed-citation publication-type="other">Stanovich, K. E. & West, R. F. Individual differences in reasoning: implications for the rationality debate? <italic>The Behavioral and brain sciences</italic>
<bold>23</bold>
, 645-665; discussion 665–726 (2000).</mixed-citation>
</ref>
<ref id="CR11"><label>11.</label>
<mixed-citation publication-type="other">Posner, M. I. & Snyder, C. R. R. <italic>Information processing and cognition: The Loyola Symposium</italic>
 (ed R. L. Solso) 55–85 (Wiley, 1975).</mixed-citation>
</ref>
<ref id="CR12"><label>12.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tversky</surname>
<given-names>A</given-names>
</name>
<name><surname>Kahneman</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Judgment under Uncertainty: Heuristics and Biases</article-title>
<source>Science</source>
<year>1974</year>
<volume>185</volume>
<fpage>1124</fpage>
<lpage>1131</lpage>
<pub-id pub-id-type="doi">10.1126/science.185.4157.1124</pub-id>
<pub-id pub-id-type="pmid">17835457</pub-id>
</element-citation>
</ref>
<ref id="CR13"><label>13.</label>
<mixed-citation publication-type="other">Kahneman, D. <italic>Thinking, fast and slow</italic>
. (Farrar, Straus & Giroux, 2011).</mixed-citation>
</ref>
<ref id="CR14"><label>14.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ghosh</surname>
<given-names>A</given-names>
</name>
<name><surname>Rho</surname>
<given-names>Y</given-names>
</name>
<name><surname>McIntosh</surname>
<given-names>AR</given-names>
</name>
<name><surname>Kotter</surname>
<given-names>R</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
</person-group>
<article-title>Cortical network dynamics with time delays reveals functional connectivity in the resting brain</article-title>
<source>Cogn Neurodyn</source>
<year>2008</year>
<volume>2</volume>
<fpage>115</fpage>
<lpage>120</lpage>
<pub-id pub-id-type="doi">10.1007/s11571-008-9044-2</pub-id>
<pub-id pub-id-type="pmid">19003478</pub-id>
</element-citation>
</ref>
<ref id="CR15"><label>15.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Deco</surname>
<given-names>G</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>V</given-names>
</name>
<name><surname>McIntosh</surname>
<given-names>AR</given-names>
</name>
<name><surname>Sporns</surname>
<given-names>O</given-names>
</name>
<name><surname>Kotter</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Key role of coupling, delay, and noise in resting brain fluctuations</article-title>
<source>Proceedings of the National Academy of Sciences of the United States of America</source>
<year>2009</year>
<volume>106</volume>
<fpage>10302</fpage>
<lpage>10307</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0901831106</pub-id>
<pub-id pub-id-type="pmid">19497858</pub-id>
</element-citation>
</ref>
<ref id="CR16"><label>16.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cabral</surname>
<given-names>J</given-names>
</name>
<name><surname>Kringelbach</surname>
<given-names>ML</given-names>
</name>
<name><surname>Deco</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Exploring the network dynamics underlying brain activity during rest</article-title>
<source>Prog Neurobiol</source>
<year>2014</year>
<volume>114</volume>
<fpage>102</fpage>
<lpage>131</lpage>
<pub-id pub-id-type="doi">10.1016/j.pneurobio.2013.12.005</pub-id>
<pub-id pub-id-type="pmid">24389385</pub-id>
</element-citation>
</ref>
<ref id="CR17"><label>17.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schirner</surname>
<given-names>M</given-names>
</name>
<name><surname>Rothmeier</surname>
<given-names>S</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
<name><surname>McIntosh</surname>
<given-names>AR</given-names>
</name>
<name><surname>Ritter</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>An automated pipeline for constructing personalized virtual brains from multimodal neuroimaging data</article-title>
<source>NeuroImage</source>
<year>2015</year>
<pub-id pub-id-type="pmid">25837600</pub-id>
</element-citation>
</ref>
<ref id="CR18"><label>18.</label>
<mixed-citation publication-type="other">Schirner, M., McIntosh, A. R., Jirsa, V., Deco, G. & Ritter, P. Bridging multiple scales in the human brain using computational modeling. <italic>bioRxiv</italic>
, <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/">http://dx.doi.org/</ext-link>
10.1101/085548 (2016).</mixed-citation>
</ref>
<ref id="CR19"><label>19.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Deco</surname>
<given-names>G</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
<name><surname>McIntosh</surname>
<given-names>AR</given-names>
</name>
</person-group>
<article-title>Emerging concepts for the dynamical organization of resting-state activity in the brain</article-title>
<source>Nature reviews. Neuroscience</source>
<year>2011</year>
<volume>12</volume>
<fpage>43</fpage>
<lpage>56</lpage>
<pub-id pub-id-type="doi">10.1038/nrn2961</pub-id>
<pub-id pub-id-type="pmid">21170073</pub-id>
</element-citation>
</ref>
<ref id="CR20"><label>20.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Deco</surname>
<given-names>G</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
<name><surname>McIntosh</surname>
<given-names>AR</given-names>
</name>
</person-group>
<article-title>Resting brains never rest: computational insights into potential cognitive architectures</article-title>
<source>Trends in neurosciences</source>
<year>2013</year>
<volume>36</volume>
<fpage>268</fpage>
<lpage>274</lpage>
<pub-id pub-id-type="doi">10.1016/j.tins.2013.03.001</pub-id>
<pub-id pub-id-type="pmid">23561718</pub-id>
</element-citation>
</ref>
<ref id="CR21"><label>21.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Damoiseaux</surname>
<given-names>JS</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Consistent resting-state networks across healthy subjects</article-title>
<source>Proceedings of the National Academy of Sciences of the United States of America</source>
<year>2006</year>
<volume>103</volume>
<fpage>13848</fpage>
<lpage>13853</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0601417103</pub-id>
<pub-id pub-id-type="pmid">16945915</pub-id>
</element-citation>
</ref>
<ref id="CR22"><label>22.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cabral</surname>
<given-names>J</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Exploring mechanisms of spontaneous functional connectivity in MEG: How delayed network interactions lead to structured amplitude envelopes of band-pass filtered oscillations</article-title>
<source>Neuroimage</source>
<year>2014</year>
<volume>90</volume>
<fpage>423</fpage>
<lpage>435</lpage>
<pub-id pub-id-type="doi">10.1016/j.neuroimage.2013.11.047</pub-id>
<pub-id pub-id-type="pmid">24321555</pub-id>
</element-citation>
</ref>
<ref id="CR23"><label>23.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Breakspear</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Dynamic models of large-scale brain activity</article-title>
<source>Nature neuroscience</source>
<year>2017</year>
<volume>20</volume>
<fpage>340</fpage>
<lpage>352</lpage>
<pub-id pub-id-type="doi">10.1038/nn.4497</pub-id>
<pub-id pub-id-type="pmid">28230845</pub-id>
</element-citation>
</ref>
<ref id="CR24"><label>24.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wildie</surname>
<given-names>M</given-names>
</name>
<name><surname>Shanahan</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Metastability and chimera states in modular delay and pulse-coupled oscillator networks</article-title>
<source>Chaos</source>
<year>2012</year>
<volume>22</volume>
<fpage>043131</fpage>
<pub-id pub-id-type="doi">10.1063/1.4766592</pub-id>
<pub-id pub-id-type="pmid">23278066</pub-id>
</element-citation>
</ref>
<ref id="CR25"><label>25.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cabral</surname>
<given-names>J</given-names>
</name>
<name><surname>Hugues</surname>
<given-names>E</given-names>
</name>
<name><surname>Sporns</surname>
<given-names>O</given-names>
</name>
<name><surname>Deco</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Role of local network oscillations in resting-state functional connectivity</article-title>
<source>NeuroImage</source>
<year>2011</year>
<volume>57</volume>
<fpage>130</fpage>
<lpage>139</lpage>
<pub-id pub-id-type="doi">10.1016/j.neuroimage.2011.04.010</pub-id>
<pub-id pub-id-type="pmid">21511044</pub-id>
</element-citation>
</ref>
<ref id="CR26"><label>26.</label>
<mixed-citation publication-type="other">Kuznetsov, Y. A. <italic>Elements of applied bifurcation theory</italic>
. (Springer, 1998).</mixed-citation>
</ref>
<ref id="CR27"><label>27.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Freyer</surname>
<given-names>F</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Biophysical mechanisms of multistability in resting-state cortical rhythms</article-title>
<source>The Journal of neuroscience: the official journal of the Society for Neuroscience</source>
<year>2011</year>
<volume>31</volume>
<fpage>6353</fpage>
<lpage>6361</lpage>
<pub-id pub-id-type="doi">10.1523/JNEUROSCI.6693-10.2011</pub-id>
<pub-id pub-id-type="pmid">21525275</pub-id>
</element-citation>
</ref>
<ref id="CR28"><label>28.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Freyer</surname>
<given-names>F</given-names>
</name>
<name><surname>Roberts</surname>
<given-names>JA</given-names>
</name>
<name><surname>Ritter</surname>
<given-names>P</given-names>
</name>
<name><surname>Breakspear</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>A canonical model of multistability and scale-invariance in biological systems</article-title>
<source>PLoS Comput Biol</source>
<year>2012</year>
<volume>8</volume>
<fpage>e1002634</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1002634</pub-id>
<pub-id pub-id-type="pmid">22912567</pub-id>
</element-citation>
</ref>
<ref id="CR29"><label>29.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hansen</surname>
<given-names>EC</given-names>
</name>
<name><surname>Battaglia</surname>
<given-names>D</given-names>
</name>
<name><surname>Spiegler</surname>
<given-names>A</given-names>
</name>
<name><surname>Deco</surname>
<given-names>G</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
</person-group>
<article-title>Functional connectivity dynamics: modeling the switching behavior of the resting state</article-title>
<source>NeuroImage</source>
<year>2015</year>
<volume>105</volume>
<fpage>525</fpage>
<lpage>535</lpage>
<pub-id pub-id-type="doi">10.1016/j.neuroimage.2014.11.001</pub-id>
<pub-id pub-id-type="pmid">25462790</pub-id>
</element-citation>
</ref>
<ref id="CR30"><label>30.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Deco</surname>
<given-names>G</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
</person-group>
<article-title>Ongoing cortical activity at rest: criticality, multistability, and ghost attractors</article-title>
<source>The Journal of neuroscience: the official journal of the Society for Neuroscience</source>
<year>2012</year>
<volume>32</volume>
<fpage>3366</fpage>
<lpage>3375</lpage>
<pub-id pub-id-type="doi">10.1523/JNEUROSCI.2523-11.2012</pub-id>
<pub-id pub-id-type="pmid">22399758</pub-id>
</element-citation>
</ref>
<ref id="CR31"><label>31.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Deco</surname>
<given-names>G</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Resting-state functional connectivity emerges from structurally and dynamically shaped slow linear fluctuations</article-title>
<source>The Journal of neuroscience: the official journal of the Society for Neuroscience</source>
<year>2013</year>
<volume>33</volume>
<fpage>11239</fpage>
<lpage>11252</lpage>
<pub-id pub-id-type="doi">10.1523/JNEUROSCI.1091-13.2013</pub-id>
<pub-id pub-id-type="pmid">23825427</pub-id>
</element-citation>
</ref>
<ref id="CR32"><label>32.</label>
<mixed-citation publication-type="other">Kelso, J. A. S. Dynamic Patterns: The Self-Organization of Brain and Behavior., (MIT Press, 1995).</mixed-citation>
</ref>
<ref id="CR33"><label>33.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Llinas</surname>
<given-names>RR</given-names>
</name>
</person-group>
<article-title>The intrinsic electrophysiological properties of mammalian neurons: insights into central nervous system function</article-title>
<source>Science</source>
<year>1988</year>
<volume>242</volume>
<fpage>1654</fpage>
<lpage>1664</lpage>
<pub-id pub-id-type="doi">10.1126/science.3059497</pub-id>
<pub-id pub-id-type="pmid">3059497</pub-id>
</element-citation>
</ref>
<ref id="CR34"><label>34.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Varela</surname>
<given-names>F</given-names>
</name>
<name><surname>Lachaux</surname>
<given-names>JP</given-names>
</name>
<name><surname>Rodriguez</surname>
<given-names>E</given-names>
</name>
<name><surname>Martinerie</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>The brainweb: phase synchronization and large-scale integration</article-title>
<source>Nature reviews. Neuroscience</source>
<year>2001</year>
<volume>2</volume>
<fpage>229</fpage>
<lpage>239</lpage>
<pub-id pub-id-type="doi">10.1038/35067550</pub-id>
<pub-id pub-id-type="pmid">11283746</pub-id>
</element-citation>
</ref>
<ref id="CR35"><label>35.</label>
<mixed-citation publication-type="other">Haken, H. Information and Self-Organization. A macroscopic approach to Complex Systems. (Springer 1988).</mixed-citation>
</ref>
<ref id="CR36"><label>36.</label>
<mixed-citation publication-type="other">Kelso, J. A. S. & Tognoli, E. <italic>Neurodynamics of Cognition and Consciousness</italic>
 (eds R. Kozma & L. Perlovsky) 39–60 (Springer, 2007).</mixed-citation>
</ref>
<ref id="CR37"><label>37.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rabinovich</surname>
<given-names>M</given-names>
</name>
<name><surname>Huerta</surname>
<given-names>R</given-names>
</name>
<name><surname>Varona</surname>
<given-names>P</given-names>
</name>
<name><surname>Afraimovich</surname>
<given-names>V</given-names>
</name>
</person-group>
<article-title>Transient Cognitive Dynamics, Metastability, and Decision Making</article-title>
<source>PLoS computational biology</source>
<year>2008</year>
<volume>4</volume>
<fpage>e1000072</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1000072</pub-id>
<pub-id pub-id-type="pmid">18452000</pub-id>
</element-citation>
</ref>
<ref id="CR38"><label>38.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rabinovich</surname>
<given-names>MI</given-names>
</name>
<name><surname>Huerta</surname>
<given-names>R</given-names>
</name>
<name><surname>Laurent.</surname>
</name>
</person-group>
<article-title>Transient Dynamics for Neural Processing</article-title>
<source>Science</source>
<year>2008</year>
<volume>321</volume>
<fpage>48</fpage>
<lpage>50</lpage>
<pub-id pub-id-type="doi">10.1126/science.1155564</pub-id>
<pub-id pub-id-type="pmid">18599763</pub-id>
</element-citation>
</ref>
<ref id="CR39"><label>39.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huys</surname>
<given-names>R</given-names>
</name>
<name><surname>Perdikis</surname>
<given-names>D</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
</person-group>
<article-title>Functional architectures and structured flows on manifolds: a dynamical framework for motor behavior</article-title>
<source>Psychol Rev</source>
<year>2014</year>
<volume>121</volume>
<fpage>302</fpage>
<lpage>336</lpage>
<pub-id pub-id-type="doi">10.1037/a0037014</pub-id>
<pub-id pub-id-type="pmid">25090422</pub-id>
</element-citation>
</ref>
<ref id="CR40"><label>40.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Perdikis</surname>
<given-names>D</given-names>
</name>
<name><surname>Huys</surname>
<given-names>R</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
</person-group>
<article-title>Time scale hierarchies in the functional organization of complex behaviors</article-title>
<source>PLoS computational biology</source>
<year>2011</year>
<volume>7</volume>
<fpage>e1002198</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1002198</pub-id>
<pub-id pub-id-type="pmid">21980278</pub-id>
</element-citation>
</ref>
<ref id="CR41"><label>41.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Beckmann</surname>
<given-names>CF</given-names>
</name>
<name><surname>DeLuca</surname>
<given-names>M</given-names>
</name>
<name><surname>Devlin</surname>
<given-names>JT</given-names>
</name>
<name><surname>Smith</surname>
<given-names>SM</given-names>
</name>
</person-group>
<article-title>Investigations into resting-state connectivity using independent component analysis</article-title>
<source>Philos. Trans. R. Soc. Lond. B Biol. Sci</source>
<year>2005</year>
<volume>360</volume>
<fpage>1001</fpage>
<lpage>1013</lpage>
<pub-id pub-id-type="doi">10.1098/rstb.2005.1634</pub-id>
<pub-id pub-id-type="pmid">16087444</pub-id>
</element-citation>
</ref>
<ref id="CR42"><label>42.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Greicius</surname>
<given-names>MD</given-names>
</name>
<name><surname>Krasnow</surname>
<given-names>B</given-names>
</name>
<name><surname>Reiss</surname>
<given-names>AL</given-names>
</name>
<name><surname>Menon</surname>
<given-names>V</given-names>
</name>
</person-group>
<article-title>Functional connectivity in the resting brain: a network analysis of the default mode hypothesis</article-title>
<source>Proc. Natl Acad. Sci. USA</source>
<year>2003</year>
<volume>100</volume>
<fpage>253</fpage>
<lpage>258</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0135058100</pub-id>
<pub-id pub-id-type="pmid">12506194</pub-id>
</element-citation>
</ref>
<ref id="CR43"><label>43.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mueller</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Individual variability in functional connectivity architecture of the human brain</article-title>
<source>Neuron</source>
<year>2013</year>
<volume>77</volume>
<fpage>586</fpage>
<lpage>595</lpage>
<pub-id pub-id-type="doi">10.1016/j.neuron.2012.12.028</pub-id>
<pub-id pub-id-type="pmid">23395382</pub-id>
</element-citation>
</ref>
<ref id="CR44"><label>44.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Allen</surname>
<given-names>EA</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Tracking whole-brain connectivity dynamics in the resting state</article-title>
<source>Cerebral cortex</source>
<year>2014</year>
<volume>24</volume>
<fpage>663</fpage>
<lpage>676</lpage>
<pub-id pub-id-type="doi">10.1093/cercor/bhs352</pub-id>
<pub-id pub-id-type="pmid">23146964</pub-id>
</element-citation>
</ref>
<ref id="CR45"><label>45.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Baker</surname>
<given-names>AP</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Fast transient networks in spontaneous human brain activity</article-title>
<source>eLife</source>
<year>2014</year>
<volume>3</volume>
<fpage>e01867</fpage>
<pub-id pub-id-type="doi">10.7554/eLife.01867</pub-id>
<pub-id pub-id-type="pmid">24668169</pub-id>
</element-citation>
</ref>
<ref id="CR46"><label>46.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ghosh</surname>
<given-names>A</given-names>
</name>
<name><surname>Rho</surname>
<given-names>Y</given-names>
</name>
<name><surname>McIntosh</surname>
<given-names>AR</given-names>
</name>
<name><surname>Kotter</surname>
<given-names>R</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
</person-group>
<article-title>Noise during rest enables the exploration of the brain’s dynamic repertoire</article-title>
<source>PLoS computational biology</source>
<year>2008</year>
<volume>4</volume>
<fpage>e1000196</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1000196</pub-id>
<pub-id pub-id-type="pmid">18846206</pub-id>
</element-citation>
</ref>
<ref id="CR47"><label>47.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ponce-Alvarez</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Resting-state temporal synchronization networks emerge from connectivity topology and heterogeneity</article-title>
<source>PLoS computational biology</source>
<year>2015</year>
<volume>11</volume>
<fpage>e1004100</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pcbi.1004100</pub-id>
<pub-id pub-id-type="pmid">25692996</pub-id>
</element-citation>
</ref>
<ref id="CR48"><label>48.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gusnard</surname>
<given-names>DA</given-names>
</name>
<name><surname>Raichle</surname>
<given-names>ME</given-names>
</name>
</person-group>
<article-title>Searching for a baseline: functional imaging and the resting human brain</article-title>
<source>Nature Reviews Neuroscience</source>
<year>2001</year>
<volume>2</volume>
<fpage>685</fpage>
<lpage>694</lpage>
<pub-id pub-id-type="doi">10.1038/35094500</pub-id>
<pub-id pub-id-type="pmid">11584306</pub-id>
</element-citation>
</ref>
<ref id="CR49"><label>49.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Addis</surname>
<given-names>DR</given-names>
</name>
<name><surname>Wong</surname>
<given-names>AT</given-names>
</name>
<name><surname>Schacter</surname>
<given-names>DL</given-names>
</name>
</person-group>
<article-title>Remembering the past and imagining the future: common and distinct neural substrates during event construction and elaboration</article-title>
<source>Neuropsychologia</source>
<year>2007</year>
<volume>45</volume>
<fpage>1363</fpage>
<lpage>1377</lpage>
<pub-id pub-id-type="doi">10.1016/j.neuropsychologia.2006.10.016</pub-id>
<pub-id pub-id-type="pmid">17126370</pub-id>
</element-citation>
</ref>
<ref id="CR50"><label>50.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Botvinick</surname>
<given-names>M</given-names>
</name>
<name><surname>Nystrom</surname>
<given-names>LE</given-names>
</name>
<name><surname>Fissell</surname>
<given-names>K</given-names>
</name>
<name><surname>Carter</surname>
<given-names>CS</given-names>
</name>
<name><surname>Cohen</surname>
<given-names>JD</given-names>
</name>
</person-group>
<article-title>Conflict monitoring versus selection-for-action in anterior cingulate cortex</article-title>
<source>Nature</source>
<year>1999</year>
<volume>402</volume>
<fpage>179</fpage>
<lpage>181</lpage>
<pub-id pub-id-type="doi">10.1038/46035</pub-id>
<pub-id pub-id-type="pmid">10647008</pub-id>
</element-citation>
</ref>
<ref id="CR51"><label>51.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Goodkind</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Identification of a Common Neurobiological Substrate for Mental Illness</article-title>
<source>JAMA psychiatry</source>
<year>2015</year>
<volume>72</volume>
<fpage>305</fpage>
<lpage>315</lpage>
<pub-id pub-id-type="doi">10.1001/jamapsychiatry.2014.2206</pub-id>
<pub-id pub-id-type="pmid">25651064</pub-id>
</element-citation>
</ref>
<ref id="CR52"><label>52.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Deco</surname>
<given-names>G</given-names>
</name>
<name><surname>Kringelbach</surname>
<given-names>ML</given-names>
</name>
</person-group>
<article-title>Great Expectations: Using Whole-Brain Computational Connectomics for Understanding Neuropsychiatric Disorders</article-title>
<source>Neuron</source>
<year>2014</year>
<volume>84</volume>
<fpage>892</fpage>
<lpage>905</lpage>
<pub-id pub-id-type="doi">10.1016/j.neuron.2014.08.034</pub-id>
<pub-id pub-id-type="pmid">25475184</pub-id>
</element-citation>
</ref>
<ref id="CR53"><label>53.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ritter</surname>
<given-names>P</given-names>
</name>
<name><surname>Schirner</surname>
<given-names>M</given-names>
</name>
<name><surname>McIntosh</surname>
<given-names>AR</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>VK</given-names>
</name>
</person-group>
<article-title>The virtual brain integrates computational modeling and multimodal neuroimaging</article-title>
<source>Brain connectivity</source>
<year>2013</year>
<volume>3</volume>
<fpage>121</fpage>
<lpage>145</lpage>
<pub-id pub-id-type="doi">10.1089/brain.2012.0120</pub-id>
<pub-id pub-id-type="pmid">23442172</pub-id>
</element-citation>
</ref>
<ref id="CR54"><label>54.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sanz Leon</surname>
<given-names>P</given-names>
</name>
<etal></etal>
</person-group>
<article-title>The Virtual Brain: a simulator of primate brain network dynamics</article-title>
<source>Frontiers in neuroinformatics</source>
<year>2013</year>
<volume>7</volume>
<fpage>10</fpage>
<pub-id pub-id-type="doi">10.3389/fninf.2013.00010</pub-id>
<pub-id pub-id-type="pmid">23781198</pub-id>
</element-citation>
</ref>
<ref id="CR55"><label>55.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ritter</surname>
<given-names>P</given-names>
</name>
<name><surname>Villringer</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Simultaneous EEG-fMRI</article-title>
<source>Neurosci.Biobehav.Rev.</source>
<year>2006</year>
<volume>30</volume>
<fpage>823</fpage>
<lpage>838</lpage>
<pub-id pub-id-type="doi">10.1016/j.neubiorev.2006.06.008</pub-id>
<pub-id pub-id-type="pmid">16911826</pub-id>
</element-citation>
</ref>
<ref id="CR56"><label>56.</label>
<mixed-citation publication-type="other">Ritter, P., Becker, R., Freyer, F. & Villringer, A. <italic>EEG-fMRI Physiology, Technique and Application</italic>
 (eds C. Mulert & L. Lemieux) Ch. 9, 153–171 (Springer, 2010).</mixed-citation>
</ref>
<ref id="CR57"><label>57.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Desikan</surname>
<given-names>RS</given-names>
</name>
<etal></etal>
</person-group>
<article-title>An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest</article-title>
<source>NeuroImage</source>
<year>2006</year>
<volume>31</volume>
<fpage>968</fpage>
<lpage>980</lpage>
<pub-id pub-id-type="doi">10.1016/j.neuroimage.2006.01.021</pub-id>
<pub-id pub-id-type="pmid">16530430</pub-id>
</element-citation>
</ref>
<ref id="CR58"><label>58.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tournier</surname>
<given-names>JD</given-names>
</name>
<name><surname>Calamante</surname>
<given-names>F</given-names>
</name>
<name><surname>Gadian</surname>
<given-names>DG</given-names>
</name>
<name><surname>Connelly</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Direct estimation of the fiber orientation density function from diffusion-weighted MRI data using spherical deconvolution</article-title>
<source>Neuroimage</source>
<year>2004</year>
<volume>23</volume>
<fpage>1176</fpage>
<lpage>1185</lpage>
<pub-id pub-id-type="doi">10.1016/j.neuroimage.2004.07.037</pub-id>
<pub-id pub-id-type="pmid">15528117</pub-id>
</element-citation>
</ref>
<ref id="CR59"><label>59.</label>
<mixed-citation publication-type="other">Kuramoto, Y. Chemical Oscillations,Waves, and Turbulence. Springer-Verlag, Berlin (1984).</mixed-citation>
</ref>
<ref id="CR60"><label>60.</label>
<mixed-citation publication-type="other">Pikovsky, A., Rosenblum, M. & Kurths, J. <italic>Synchronization</italic>
: <italic>a universal concept in nonlinear sciences</italic>
. (Cambridge University Press, 2003).</mixed-citation>
</ref>
<ref id="CR61"><label>61.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Matthews</surname>
<given-names>PC</given-names>
</name>
<name><surname>Strogatz</surname>
<given-names>SH</given-names>
</name>
</person-group>
<article-title>Phase diagram for the collective behavior of limit-cycle oscillators</article-title>
<source>Phys. Rev. Lett.</source>
<year>1990</year>
<volume>65</volume>
<fpage>1701</fpage>
<lpage>1704</lpage>
<pub-id pub-id-type="doi">10.1103/PhysRevLett.65.1701</pub-id>
<pub-id pub-id-type="pmid">10042341</pub-id>
</element-citation>
</ref>
<ref id="CR62"><label>62.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aronson</surname>
<given-names>DG</given-names>
</name>
<name><surname>Ermentrout</surname>
<given-names>GB</given-names>
</name>
<name><surname>Kopell</surname>
<given-names>N</given-names>
</name>
</person-group>
<article-title>Amplitude response of coupled oscillators</article-title>
<source>Physica D: Nonlinear Phenomena</source>
<year>1990</year>
<volume>41</volume>
<fpage>403</fpage>
<lpage>449</lpage>
<pub-id pub-id-type="doi">10.1016/0167-2789(90)90007-C</pub-id>
</element-citation>
</ref>
<ref id="CR63"><label>63.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Glerean</surname>
<given-names>E</given-names>
</name>
<name><surname>Salmi</surname>
<given-names>J</given-names>
</name>
<name><surname>Lahnakoski</surname>
<given-names>JM</given-names>
</name>
<name><surname>Jaaskelainen</surname>
<given-names>IP</given-names>
</name>
<name><surname>Sams</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Functional magnetic resonance imaging phase synchronization as a measure of dynamic functional connectivity</article-title>
<source>Brain connectivity</source>
<year>2012</year>
<volume>2</volume>
<fpage>91</fpage>
<lpage>101</lpage>
<pub-id pub-id-type="doi">10.1089/brain.2011.0068</pub-id>
<pub-id pub-id-type="pmid">22559794</pub-id>
</element-citation>
</ref>
<ref id="CR64"><label>64.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Biswal</surname>
<given-names>B</given-names>
</name>
<name><surname>Yetkin</surname>
<given-names>FZ</given-names>
</name>
<name><surname>Haughton</surname>
<given-names>VM</given-names>
</name>
<name><surname>Hyde</surname>
<given-names>JS</given-names>
</name>
</person-group>
<article-title>Functional connectivity in the motor cortex of resting human brain using echo-planar MRI</article-title>
<source>Magn.Reson.Med.</source>
<year>1995</year>
<volume>34</volume>
<fpage>537</fpage>
<lpage>541</lpage>
<pub-id pub-id-type="doi">10.1002/mrm.1910340409</pub-id>
<pub-id pub-id-type="pmid">8524021</pub-id>
</element-citation>
</ref>
<ref id="CR65"><label>65.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Buckner</surname>
<given-names>RL</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Cortical hubs revealed by intrinsic functional connectivity: mapping, assessment of stability, and relation to Alzheimer’s disease</article-title>
<source>J Neurosci</source>
<year>2009</year>
<volume>29</volume>
<fpage>1860</fpage>
<lpage>1873</lpage>
<pub-id pub-id-type="doi">10.1523/JNEUROSCI.5062-08.2009</pub-id>
<pub-id pub-id-type="pmid">19211893</pub-id>
</element-citation>
</ref>
<ref id="CR66"><label>66.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Achard</surname>
<given-names>S</given-names>
</name>
<name><surname>Salvador</surname>
<given-names>R</given-names>
</name>
<name><surname>Whitcher</surname>
<given-names>B</given-names>
</name>
<name><surname>Suckling</surname>
<given-names>J</given-names>
</name>
<name><surname>Bullmore</surname>
<given-names>E</given-names>
</name>
</person-group>
<article-title>A resilient, low-frequency, small-world human brain functional network with highly connected association cortical hubs</article-title>
<source>The Journal of neuroscience: the official journal of the Society for Neuroscience</source>
<year>2006</year>
<volume>26</volume>
<fpage>63</fpage>
<lpage>72</lpage>
<pub-id pub-id-type="doi">10.1523/JNEUROSCI.3874-05.2006</pub-id>
<pub-id pub-id-type="pmid">16399673</pub-id>
</element-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>

Pour manipuler ce document sous Unix (Dilib)

EXPLOR_STEP=$WICRI_ROOT/Wicri/Asie/explor/AustralieFrV1/Data/Pmc/Corpus

HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000E20 | SxmlIndent | more

HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000E20 | SxmlIndent | more

Pour mettre un lien sur cette page dans le réseau Wicri

{{Explor lien
   |wiki=    Wicri/Asie
   |area=    AustralieFrV1
   |flux=    Pmc
   |étape=   Corpus
   |type=    RBID
   |clé=     PMC:5465179
   |texte=   The dynamics of resting fluctuations in the brain: metastability and its dynamical cortical core
}}

Pour générer des pages wiki

HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i   -Sk "pubmed:28596608" \
       | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd   \
       | NlmPubMed2Wicri -a AustralieFrV1

This area was generated with Dilib version V0.6.33.
Data generation: Tue Dec 5 10:43:12 2017. Site generation: Tue Mar 5 14:07:20 2024

	Serveur d'exploration sur les relations entre la France et l'Australie
	Attention, ce site est en cours de développement ! Attention, site généré par des moyens informatiques à partir de corpus bruts. Les informations ne sont donc pas validées.

Serveur d'exploration sur les relations entre la France et l'Australie

The dynamics of resting fluctuations in the brain: metastability and its dynamical cortical core

The dynamics of resting fluctuations in the brain: metastability and its dynamical cortical core

Source :

Abstract

Links to Exploration step

Le document en format XML

Pour manipuler ce document sous Unix (Dilib)

Pour mettre un lien sur cette page dans le réseau Wicri

Pour générer des pages wiki