Changes for page SGA3 D1.2 Showcase 1
Last modified by fousekjan on 2022/07/04 18:31
Summary
-
Page properties (1 modified, 0 added, 0 removed)
-
Attachments (0 modified, 4 added, 0 removed)
Details
- Page properties
-
- Content
-
... ... @@ -12,38 +12,47 @@ 12 12 ((( 13 13 (% class="col-xs-12 col-sm-8" %) 14 14 ((( 15 -(% class="wikigeneratedid" %) 16 -== Siibra - Python client for accessing the human atlas == 15 +The Showcase is implemented as a series of interactive Jupyter notebooks covering the individual logical steps and can be accessed in a dedicated public EBRAINS collab. 17 17 18 -The multilevel atlas interface called siibra (formerly known as Brainscapes) is designed to allow safe and convenient interaction with brain definitions from different parcellations, facilitating the implementation of reproducible neuroscience workflows on the basis of brain atlases. It allows to work with reference brain templates both at millimeter and micrometer resolutions and provides streamlined access to multimodal data features linked to brain regions. In particular in this demonstrator it is used to retrieve the dataset of receptor density spatial maps. 19 19 20 -The S howcase1a isbuilt onthe[[connectivitydatafromthe1000BRAINSstudy>>https://doi.org/10.25493%2F6640-3XH]]available in theKnowledge graph as adatasetwithprotectedaccessavailable toEBRAINS users. The accesstoprotecteddatasetsisprovidedinEBRAINS through the[[HumanDataGateway>>https://wiki.ebrains.eu/bin/view/Collabs/data-proxy/Human%20Data%20Gateway/]]whichallowsfullaccessoncetheuserhasvalidatedthe termsofuse.18 +The EBRAINS collab consists of interlinked Drive, Bucket, Wiki, and Lab. The Drive provides small file storage and contains the notebooks and all supporting code. The Bucket is a large file storage service and holds the pre-computed results of the extensive parameter sweeps and model optimizations to allow skipping the computationally demanding steps. The documentation of the showcase implementation is collected in the Wiki. The Lab service is an instance of JupyterLab—an interactive computing environment where the notebooks can be run and worked with. 21 21 22 -== Brain Network Modeling: TVB == 23 23 24 -The second major componentoftheshowcaseisTheVirtualBrain(TVB) simulationplatform. TVB isa simulationplatformthatusesempirical structural and functionaldata to buildwholebrainmodels ofindividual subjects.Forconvenientmodelconstruction,thesystemis basedonaprocessingpipelineforstructural,functional, and diffusion-weightedmagneticresonanceimaging(MRI)data.Thepipelinecombinesseveralstate-of-the-artneuroinformaticstoolstogeneratesubject-specific cortical andsubcorticalparcellations,surface-tessellations,structural andfunctionalconnectomes, andregion-wiseaggregatedblood oxygen level-dependent (BOLD) functionalMRI(fMRI) time-series.21 +The notebooks in this collab will load all required Python modules including Siibra and The Virtual Brain, and the interfaces for launching the computationally demanding parts in the HPC infrastructure. Running the notebooks requires an EBRANS account with permissions to access the Lab and the Knowledge Graph API. In addition, to be able to interact with the HPC infrastructure, the user has to have access to an active allocation on the corresponding FENIX site. 25 25 26 -== Virtual ageing trajectories: 1000BRAINS cohort == 23 +== 24 +Simulation of resting state == 27 27 28 - (% style="font-size:11pt;font-variant:normal; white-space: pre-wrap; font-family: Arial; color: rgb(0, 0, 0); font-weight: 400; font-style:normal; text-decoration:none" %)Beyondthegeneral common architecture, the humanbrain is characterized by a high interindividual variability withregard toitsstructureand functional abilities,includingmeasurable cognitiveoutcome. Disentanglingthe sourcesandconsequencesof thisvariability iskeytounderstandingthecapabilities ofourbrainsand any pathologies occurringoverthelifespan.Thisalsoincludesimplementingmeasuresofvariabilityinto simulationsandmodellingapproachestomakethembiologicallyplausibleandasrealistic aspossible.26 +The fist notebook in the inter-individual variability workflow explores the resting-state simulation for a subject of the 1000BRAINS dataset. Functional data are simulated by means of a brain network model implemented in TVB, which is an ensemble of neural mass models linked via the weights of the structural connectivity (SC) matrix. Following topics are covered: 29 29 30 -(% style="line-height:1.38; margin-top:13px" %) 31 -(% style="font-size: 11pt; font-variant: normal; white-space: pre-wrap; font-family: Arial; color: rgb(0, 0, 0); font-weight: 400; font-style: normal; text-decoration: none" %)The aging brain is a particularly remarkable use case in this respect as the interindividual variability of brains and cognitive abilities of older subjects from the general population is very pronounced. The relevant factors influencing this variability, be it genetics or environmental and lifestyle factors, are manifold and barely understood yet. Similarly, the different levels within the organization of the brain, from the molecular, cellular to the systems level, contribute to these effects to varying degrees. The aim of the first branch of Showcase 1 is thus to model brain aging with biologically plausible assumptions about the architecture of the brain, using the EBRAINS platform and integrating information from the different organizational levels of the brain, and to create a virtually aged cohort of brain models which allow testing of hypotheses and causal linkage of potential sources and effects of the empirically observed variability. 28 +1. The neural mass model by Montbrio, Pazo and Roxin 29 +1. Construction of the TVB model for a particular subject 30 +1. Dynamics of the model, and execution of a parameter study 31 +1. Summary of the simulation results for the whole cohort 32 32 33 - === Technical specification===33 +[[image:image-20220103100841-2.png]] 34 34 35 35 36 -== Brain region variance:incorporating receptordensities==36 +== Regional variability data == 37 37 38 - (% style="font-size:11pt;font-variant:normal; white-space:pre-wrap; font-family:Arial; color: rgb(0, 0, 0); font-weight: 400; font-style:normal; text-decoration: none" %)The second branchof the Showcase 1 focuses on the impact ofvariability of the properties ofbrain regions on the expressed dynamics. Over thelast twodecades,whole-brainnetwork models have proven useful instruments todescribethebrain's resting-stateactivity:fromexplainingtheintrinsic dynamicalbehavior of thebrain to theclassificationof brain states indisorders ofconsciousness.These whole-brain models are builtat the resolution of interconnected cortical regionsand subcortical brainareas,with each regionsimulated at thepopulation levels. So far, amajor limitation of these models lies in the assumption that all brain regionsare of identical characteristics whileit iswell-knownthat brain areas differin a varietyofproperties, e.g., neuronal densities, local cytoarchitecture,andtypesofneuroreceptors [Deco2021]. Accounting for these regionalvariancesis crucial fortranslational purposesasforexample differencesinneuroreceptor densitiesdeterminetheimpact of pharmacologicalagents on particular brain regionsandthustheirpotentialto modulatethe whole-braindynamics. Thepurpose ofShowcase #1(b) is,therefore,toillustrate how whole-brain network modelscan be adapted toinclude such regionalvariances and todevelop workflowsto run the simulations takingfull advantageoftheHBP'sEBRAINSinfrastructure, leadingto workflows that are accessibleandre-usable bythe scientific community.38 +The first step of the regional variability workflow loads the data from the Knowledge Graph and defines the regional bias on the model. In this case we require: 1) structural connectivity matrices, 2) GABA and AMPA receptor densities at each region, and 3) empirical resting-state fMRI data for fitting and validation of the simulations. The three datasets shall be characterised in the same parcellation. 39 39 40 +Link to the notebook: 40 40 41 - ===Technicalspecification===42 +* [[regional_variability/notebooks/1_data_setup.ipynb>>https://lab.ch.ebrains.eu/user-redirect/lab/tree/shared/SGA3%20D1.2%20Showcase%201/regional_variability/notebooks/1_data_setup.ipynb]] 42 42 43 -(% style=" line-height:1.38; margin-top:13px" %)44 - (% style="font-size: 11pt; font-variant: normal; white-space: pre-wrap; font-family: Arial; color: rgb(0, 0, 0); font-weight:400; font-style: normal; text-decoration: none" %)In implementation of Showcase 1b, the simulation is based on a network of nodes that represent brain ROI, each runningthe(% style="font-size: 11pt; font-variant: normal; white-space: pre-wrap; font-family: Arial; color: rgb(0,0, 0); font-weight: 700; font-style: normal; text-decoration: none" %)**BEI**(% style="font-size:11pt; font-variant: normal; white-space: pre-wrap; font-family: Arial; color: rgb(0,0, 0); font-weight:400; font-style: normal; text-decoration: none" %) (Balanced Excitation-Inhibition) model [Deco2014], which defines two neuronal subpopulations, an excitatory and an inhibitory one, both interacting through neurochemical (e.g., AMPA, GABA, NMDA, etc.) currents. The objective of the Showcase 1b is to demonstrate and find out the more precise mechanism through which regional neurotransmitter density affects this excitation-inhibition balance. Numerically, the goal is to minimize the difference in sliding window Functional Connectivity Dynamics ((% style="font-size:11pt; font-variant: normal; white-space: pre-wrap; font-family: Arial; color: rgb(0, 0, 0); font-weight: 700; font-style: normal; text-decoration: none" %)**swFCD**(% style="font-size: 11pt; font-variant: normal; white-space: pre-wrap; font-family: Arial; color: rgb(0, 0, 0); font-weight: 400; font-style: normal; text-decoration: none" %)) between a reference empirical signal and a simulated one.The swFCD subdivides the time-series into successive windows and, for each one, computes its correspondingFunctional Connectivity through the usual Pearson correlation, resulting in a series of NxN matrices. Two swFCD can be compared by means of the Kolmogorov-Smirnov statistic.44 +(% style="text-align:center" %) 45 +[[image:image-20220103095424-1.png]] 45 45 46 - 47 +== Fitting model parameters for models with regional bias == 48 + 49 +A series of simulations of the whole-brain network model is launched in the EBRAINS HPC facilities in order to identify the optimal model parameters leading to simulated resting-state brain activity that best resembles the empirically observed activity. 50 + 51 +Link to the notebook: 52 + 53 +* [[regional_variability/notebooks/2_parameter_swep_proto.ipynb>>https://lab.ch.ebrains.eu/user-redirect/lab/tree/shared/SGA3%20D1.2%20Showcase%201/regional_variability/notebooks/2_parameter_swep_proto.ipynb]] 54 + 55 +[[image:image-20220103095948-1.png]] 47 47 ))) 48 48 49 49
- image-20220103095424-1.png
-
- Author
-
... ... @@ -1,0 +1,1 @@ 1 +XWiki.fousekjan - Size
-
... ... @@ -1,0 +1,1 @@ 1 +33.7 KB - Content
- image-20220103095948-1.png
-
- Author
-
... ... @@ -1,0 +1,1 @@ 1 +XWiki.fousekjan - Size
-
... ... @@ -1,0 +1,1 @@ 1 +4.1 KB - Content
- image-20220103100841-2.png
-
- Author
-
... ... @@ -1,0 +1,1 @@ 1 +XWiki.fousekjan - Size
-
... ... @@ -1,0 +1,1 @@ 1 +116.4 KB - Content
- image-20220103101022-3.png
-
- Author
-
... ... @@ -1,0 +1,1 @@ 1 +XWiki.fousekjan - Size
-
... ... @@ -1,0 +1,1 @@ 1 +87.1 KB - Content