Abstract
Whole-brain network modeling of epilepsy is a data-driven approach that combines personalized anatomical information with dynamical models of abnormal brain activity to generate spatio-temporal seizure patterns as observed in brain imaging signals. Such a parametric simulator is equipped with a stochastic generative process, which itself provides the basis for inference and prediction of the local and global brain dynamics affected by disorders. However, the calculation of likelihood function at whole-brain scale is often intractable. Thus, likelihood-free inference algorithms are required to efficiently estimate the parameters pertaining to the hypothetical areas in the brain, ideally including the uncertainty. In this detailed study, we present simulation-based inference for the virtual epileptic patient (SBI-VEP) model, which only requires forward simulations, enabling us to amortize posterior inference on parameters from low-dimensional data features representing whole-brain epileptic patterns. We use state-of-the-art deep learning algorithms for conditional density estimation to retrieve the statistical relationships between parameters and observations through a sequence of invertible transformations. This approach enables us to readily predict seizure dynamics from new input data. We show that the SBI-VEP is able to accurately estimate the posterior distribution of parameters linked to the extent of the epileptogenic and propagation zones in the brain from the sparse observations of intracranial EEG signals. The presented Bayesian methodology can deal with non-linear latent dynamics and parameter degeneracy, paving the way for reliable prediction of neurological disorders from neuroimaging modalities, which can be crucial for planning intervention strategies.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
This work was funded by the French National Research Agency (ANR; https://anr.fr) as part of the second "Investissements d'Avenir" program, ANR-17-RHUS-0004, EPINOV, to VJ and FB, by European Union's Horizon 2020 Framework Programme for Research and Innovation under the Specific Grant Agreement No. 785907 and 945539, Human Brain Project SGA2 and SGA3 (https://ec.europa.eu/programmes/horizon2020) to VJ, by European Union's Horizon 2020 Framework Programme for Research and Innovation under the Specific Grant Agreement and No. 826421, VirtualBrainCloud (https://ec.europa.eu/programmes/horizon2020) to VJ, PHRC-I 2013 EPISODIUM (grant number 2014-27), the Fondation pour la Recherche Medicale (DIC20161236442), and by SATT Sud-Est, 827-SA-16-UAM (https://www.sattse.com) to VJ and FB. The funders had no role in study design, data collection, and analysis, decision to publish, or preparation of the manuscript.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
Informed written consent was obtained for all patients in compliance with the ethical requirements of the Declaration of Helsinki and the study protocol was approved by the local Ethics Committee (Comite de Protection des Personnes sud Mediterranee 1).
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Data Availability
The patient data sets cannot be made publicly available due to the data protection concerns regarding potentially identifying and sensitive patient information. Interested researchers may access the data sets by contacting Clinical Data Manager Aurelie Ponz (aurelie.ponz@univ-amu.fr) at the Institut de Neurosciences des Systemes, Aix-Marseille Universite. The main source codes needed to reproduce the presented results are available on GitHub (https://github.com/ins-amu/BVEP).
Abbreviations
- VEP
- virtual epileptic patient
- TVB
- the virtual brain
- EZ
- epileptogenic zone
- PZ
- propagation zone
- HZ
- healthy zone
- SEEG
- stereoelectroencephalography
- MCMC
- Markov chain Monte Carlo
- HMC
- Hamiltonian Monte Carlo
- VI
- variational inference
- SBI
- simulation-based-inference
- ANNs
- artificial neural networks
- NFs
- normalizing flows
- MDN
- mixture density network
- NSFs
- neural spline flows
- MAF
- masked autoregressive flow
- SNPE
- sequential neural posterior estimation
- SNLE
- sequential neural likelihood estimation
- SNRE
- sequential neural ratio estimation