CiftiStorm pipeline: facilitating reproducible EEG/MEG source connectomics

We present CiftiStorm, an electrophysiological source imaging (ESI) pipeline incorporating recently developed methods to improve forward and inverse solutions. The CiftiStorm pipeline produces Human Connectome Project (HCP) and megconnectome-compliant outputs from dataset inputs with varying degrees of spatial resolution. The input data can range from low-sensor-density electroencephalogram (EEG) or magnetoencephalogram (MEG) recordings without structural magnetic resonance imaging (sMRI) to high-density EEG/MEG recordings with an HCP multimodal sMRI compliant protocol. CiftiStorm introduces a numerical quality control of the lead field and geometrical corrections to the head and source models for forward modeling. For the inverse modeling, we present a Bayesian estimation of the cross-spectrum of sources based on multiple priors. We facilitate ESI in the T1w/FSAverage32k high-resolution space obtained from individual sMRI. We validate this feature by comparing CiftiStorm outputs for EEG and MRI data from the Cuban Human Brain Mapping Project (CHBMP) acquired with technologies a decade before the HCP MEG and MRI standardized dataset.

Spectral quantitative and semi-quantitative EEG provide complementary information on the life-long effects of early childhood malnutrition on cognitive decline

Objective

This study compares the complementary information from semi-quantitative EEG (sqEEG) and spectral quantitative EEG (spectral-qEEG) to detect the life-long effects of early childhood malnutrition on the brain.

Methods

Resting-state EEGs (N = 202) from the Barbados Nutrition Study (BNS) were used to examine the effects of protein-energy malnutrition (PEM) on childhood and middle adulthood outcomes. sqEEG analysis was performed on Grand Total EEG (GTE) protocol, and a single latent variable, the semi-quantitative Neurophysiological State (sqNPS) was extracted. A univariate linear mixed-effects (LME) model tested the dependence of sqNPS and nutritional group. sqEEG was compared with scores on the Montreal Cognitive Assessment (MoCA). Stable sparse classifiers (SSC) also measured the predictive power of sqEEG, spectral-qEEG, and a combination of both. Multivariate LME was applied to assess each EEG modality separately and combined under longitudinal settings.

Results

The univariate LME showed highly significant differences between previously malnourished and control groups (p < 0.001); age (p = 0.01) was also significant, with no interaction between group and age detected. Childhood sqNPS (p = 0.02) and adulthood sqNPS (p = 0.003) predicted MoCA scores in adulthood. The SSC demonstrated that spectral-qEEG combined with sqEEG had the highest predictive power (mean AUC 0.92 ± 0.005). Finally, multivariate LME showed that the combined spectral-qEEG+sqEEG models had the highest log-likelihood (-479.7).

Conclusion

This research has extended our prior work with spectral-qEEG and the long-term impact of early childhood malnutrition on the brain. Our findings showed that sqNPS was significantly linked to accelerated cognitive aging at 45-51 years of age. While sqNPS and spectral-qEEG produced comparable results, our study indicated that combining sqNPS and spectral-qEEG yielded better performance than either method alone, suggesting that a multimodal approach could be advantageous for future investigations.

Significance

Based on our findings, a semi-quantitative approach utilizing GTE could be a valuable diagnostic tool for detecting the lasting impacts of childhood malnutrition. Notably, sqEEG has not been previously explored or reported as a biomarker for assessing the longitudinal effects of malnutrition. Furthermore, our observations suggest that sqEEG offers unique features and information not captured by spectral quantitative EEG analysis and could lead to its improvement.

Short and long term prediction of seizure freedom in drug-resistant focal epilepsy surgery

BACKGROUND

The selection of candidates for drug-resistant focal epilepsy surgery is essential to achieve the best post-surgical outcomes.

OBJECTIVE

To develop two prediction models for seizure freedom in the short and long-term follow-up and from them to create a risk calculator in order to individualize the selection of candidates for surgery and future therapies in each patients.

METHODS

A sample of 64 consecutive patients who underwent epilepsy surgery at two Cuban tertiary health institutions between 2012 and 2020 constituted the basis for the prediction models. Two models were obtained through the novel methodology, based on biomarker selection reached by resampling methods, cross-validation and high-accuracy index measured through the area under the receiving operating curve (ROC) procedure.

RESULTS

The first, to pre-operative model included five predictors: epilepsy type, seizures per month, ictal pattern, interictal EEG topography and normal or abnormal magnetic resonance imaging,. it's precision was 0.77 at one year, and with four years and more 0.63. The second model including variables from the trans-surgical and post-surgical stages: the interictal discharges in the post-surgical EEG, incomplete or complete resection of the epileptogenic zone, the surgical techniques employed and disappearance of the discharge in post-resection electrocorticography; the precision of this model was 0.82 at one year, and with four years and more 0.97.

CONCLUSIONS

The introduction of trans-surgical and post-surgical variables increase the prediction of the pre-surgical model. A risk calculator was developed using these prediction models, which could be useful as an accurate tool to improve the prediction in epilepsy surgery.

Harmonized-Multinational qEEG norms (HarMNqEEG)

This paper extends frequency domain quantitative electroencephalography (qEEG) methods pursuing higher sensitivity to detect Brain Developmental Disorders. Prior qEEG work lacked integration of cross-spectral information omitting important functional connectivity descriptors. Lack of geographical diversity precluded accounting for site-specific variance, increasing qEEG nuisance variance. We ameliorate these weaknesses. (i) Create lifespan Riemannian multinational qEEG norms for cross-spectral tensors. These norms result from the HarMNqEEG project fostered by the Global Brain Consortium. We calculate the norms with data from 9 countries, 12 devices, and 14 studies, including 1564 subjects. Instead of raw data, only anonymized metadata and EEG cross-spectral tensors were shared. After visual and automatic quality control, developmental equations for the mean and standard deviation of qEEG traditional and Riemannian DPs were calculated using additive mixed-effects models. We demonstrate qEEG "batch effects" and provide methods to calculate harmonized z-scores. (ii) We also show that harmonized Riemannian norms produce z-scores with increased diagnostic accuracy predicting brain dysfunction produced by malnutrition in the first year of life and detecting COVID induced brain dysfunction. (iii) We offer open code and data to calculate different individual z-scores from the HarMNqEEG dataset. These results contribute to developing bias-free, low-cost neuroimaging technologies applicable in various health settings.

Brain charts for the human lifespan

Over the past few decades, neuroimaging has become a ubiquitous tool in basic research and clinical studies of the human brain. However, no reference standards currently exist to quantify individual differences in neuroimaging metrics over time, in contrast to growth charts for anthropometric traits such as height and weight1. Here we assemble an interactive open resource to benchmark brain morphology derived from any current or future sample of MRI data ( http://www.brainchart.io/ ). With the goal of basing these reference charts on the largest and most inclusive dataset available, acknowledging limitations due to known biases of MRI studies relative to the diversity of the global population, we aggregated 123,984 MRI scans, across more than 100 primary studies, from 101,457 human participants between 115 days post-conception to 100 years of age. MRI metrics were quantified by centile scores, relative to non-linear trajectories2 of brain structural changes, and rates of change, over the lifespan. Brain charts identified previously unreported neurodevelopmental milestones3, showed high stability of individuals across longitudinal assessments, and demonstrated robustness to technical and methodological differences between primary studies. Centile scores showed increased heritability compared with non-centiled MRI phenotypes, and provided a standardized measure of atypical brain structure that revealed patterns of neuroanatomical variation across neurological and psychiatric disorders. In summary, brain charts are an essential step towards robust quantification of individual variation benchmarked to normative trajectories in multiple, commonly used neuroimaging phenotypes.

Early protein energy malnutrition impacts life-long developmental trajectories of the sources of EEG rhythmic activity

Protein Energy Malnutrition (PEM) has lifelong consequences on brain development and cognitive function. We studied the lifelong developmental trajectories of resting-state EEG source activity in 66 individuals with histories of Protein Energy Malnutrition (PEM) limited to the first year of life and in 83 matched classmate controls (CON) who are all participants of the 49 years longitudinal Barbados Nutrition Study (BNS). qEEGt source z-spectra measured deviation from normative values of EEG rhythmic activity sources at 5-11 years of age and 40 years later at 45-51 years of age. The PEM group showed qEEGt abnormalities in childhood, including a developmental delay in alpha rhythm maturation and an insufficient decrease in beta activity. These profiles may be correlated with accelerated cognitive decline.

The Cuban Human Brain Mapping Project, a young and middle age population-based EEG, MRI, and cognition dataset

The Cuban Human Brain Mapping Project (CHBMP) repository is an open multimodal neuroimaging and cognitive dataset from 282 young and middle age healthy participants (31.9 ± 9.3 years, age range 18–68 years). This dataset was acquired from 2004 to 2008 as a subset of a larger stratified random sample of 2,019 participants from La Lisa municipality in La Habana, Cuba. The exclusion criteria included the presence of disease or brain dysfunctions. Participant data that is being shared comprises i) high-density (64–120 channels) resting-state electroencephalograms (EEG), ii) magnetic resonance images (MRI), iii) psychological tests (MMSE, WAIS-III, computerized go-no go reaction time), as well as iv,) demographic information (age, gender, education, ethnicity, handedness, and weight). The EEG data contains recordings with at least 30 minutes in duration including the following conditions: eyes closed, eyes open, hyperventilation, and subsequent recovery. The MRI consists of anatomical T1 as well as diffusion-weighted (DWI) images acquired on a 1.5 Tesla system. The dataset presented here is hosted by Synapse.org and available at https://chbmp-open.loris.ca.

Measurement(s)	functional brain measurement
Technology Type(s)	electroencephalography (EEG) • magnetic resonance imaging (MRI) • neuropsychological testing
Factor Type(s)	age of participants • gender of participants • handedness of participants • educational level of participants
Sample Characteristic - Organism	Homo sapiens
Sample Characteristic - Location	Cuba

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.13277348

A Quantitative EEG Toolbox for the MNI Neuroinformatics Ecosystem: Normative SPM of EEG Source Spectra

The Tomographic Quantitative Electroencephalography (qEEGt) toolbox is integrated with the Montreal Neurological Institute (MNI) Neuroinformatics Ecosystem as a docker into the Canadian Brain Imaging Research Platform (CBRAIN). qEEGt produces age-corrected normative Statistical Parametric Maps of EEG log source spectra testing compliance to a normative database. This toolbox was developed at the Cuban Neuroscience Center as part of the first wave of the Cuban Human Brain Mapping Project (CHBMP) and has been validated and used in different health systems for several decades. Incorporation into the MNI ecosystem now provides CBRAIN registered users access to its full functionality and is accompanied by a public release of the source code on GitHub and Zenodo repositories. Among other features are the calculation of EEG scalp spectra, and the estimation of their source spectra using the Variable Resolution Electrical Tomography (VARETA) source imaging. Crucially, this is completed by the evaluation of z spectra by means of the built-in age regression equations obtained from the CHBMP database (ages 5-87) to provide normative Statistical Parametric Mapping of EEG log source spectra. Different scalp and source visualization tools are also provided for evaluation of individual subjects prior to further post-processing. Openly releasing this software in the CBRAIN platform will facilitate the use of standardized qEEGt methods in different research and clinical settings. An updated precis of the methods is provided in Appendix I as a reference for the toolbox. qEEGt/CBRAIN is the first installment of instruments developed by the neuroinformatic platform of the Cuba-Canada-China (CCC) project.