Predicting Age From Brain EEG Signals—A Machine Learning Approach

Objective: The brain age gap estimate (BrainAGE) is the difference between the estimated age and the individual chronological age. BrainAGE was studied primarily using MRI techniques. EEG signals in combination with machine learning (ML) approaches were not commonly used for the human age prediction, and BrainAGE. We investigated whether age-related changes are affecting brain EEG signals, and whether we can predict the chronological age and obtain BrainAGE estimates using a rigorous ML framework with a novel and extensive EEG features extraction.Methods: EEG data were obtained from 468 healthy, mood/anxiety, eating and substance use disorder participants (297 females) from the Tulsa-1000, a naturalistic longitudinal study based on Research Domain Criteria framework. Five sets of preprocessed EEG features across channels and frequency bands were used with different ML methods to predict age. Using a nested-cross-validation (NCV) approach and stack-ensemble learning from EEG features, the predicted age was estimated. The important features and their spatial distributions were deduced.Results: The stack-ensemble age prediction model achieved R2 = 0.37 (0.06), Mean Absolute Error (MAE) = 6.87(0.69) and RMSE = 8.46(0.59) in years. The age and predicted age correlation was r = 0.6. The feature importance revealed that age predictors are spread out across different feature types. The NCV approach produced a reliable age estimation, with features consistent behavior across different folds.Conclusion: Our rigorous ML framework and extensive EEG signal features allow a reliable estimation of chronological age, and BrainAGE. This general framework can be extended to test EEG association with and to predict/study other physiological relevant responses

The unique face of comorbid anxiety and depression: increased interoceptive fearfulness and reactivity

Anxiety and depression commonly co-occur, yet the underlying brain and behavioral processes are poorly understood. Here we examined the hypothesis that individuals with comorbid anxiety and depression would show increased fearful reactivity to an aversive interoceptive perturbation relative to depressed-only individuals. One-hundred and eighty anxious and/or depressed participants from the Tulsa 1000 study completed multi-level behavioral or functional magnetic resonance imaging assessments of interoception and nociception including breath-hold and cold-pressor challenges, and heartbeat perception and interoceptive attention tasks. One-hundred and four individuals with comorbid depression and anxiety disorders (Dep+Anx) were propensity matched with 52 individuals with depression-only (Dep). Data were analyzed using mixed-effects linear regression. The Dep+Anx group showed significantly greater self-reported fear of suffocation during breath holding (Wilcoxon r = 0.23) and reduced cold pain tolerance (R2 = 0.027) signified by hand removal during immersion. However, these groups did not differ with respect to neutrally-valenced behavioral indices of heartbeat perception or neural indices of interoceptive attention. Individuals with comorbid depression and anxiety, vs. those with only depression, show increased respiratory fearfulness and nociceptive reactivity during perturbations of these signals, whilst showing similar interoceptive awareness in the absence of perturbation. Our findings suggest that individuals with comorbid anxiety and depression process aversive interoceptive and nociceptive signals differently than those with depression alone, providing support for a process model of increased threat sensitivity and hyperarousal in anxious depression

A Nonlinear Simulation Framework Supports Adjusting for Age When Analyzing BrainAGE

Several imaging modalities, including T1-weighted structural imaging, diffusion tensor imaging, and functional MRI can show chronological age related changes. Employing machine learning algorithms, an individual's imaging data can predict their age with reasonable accuracy. While details vary according to modality, the general strategy is to: (1) extract image-related features, (2) build a model on a training set that uses those features to predict an individual's age, (3) validate the model on a test dataset, producing a predicted age for each individual, (4) define the “Brain Age Gap Estimate” (BrainAGE) as the difference between an individual's predicted age and his/her chronological age, (5) estimate the relationship between BrainAGE and other variables of interest, and (6) make inferences about those variables and accelerated or delayed brain aging. For example, a group of individuals with overall positive BrainAGE may show signs of accelerated aging in other variables as well. There is inevitably an overestimation of the age of younger individuals and an underestimation of the age of older individuals due to “regression to the mean.” The correlation between chronological age and BrainAGE may significantly impact the relationship between BrainAGE and other variables of interest when they are also related to age. In this study, we examine the detectability of variable effects under different assumptions. We use empirical results from two separate datasets [training = 475 healthy volunteers, aged 18–60 years (259 female); testing = 489 participants including people with mood/anxiety, substance use, eating disorders and healthy controls, aged 18–56 years (312 female)] to inform simulation parameter selection. Outcomes in simulated and empirical data strongly support the proposal that models incorporating BrainAGE should include chronological age as a covariate. We propose either including age as a covariate in step 5 of the above framework, or employing a multistep procedure where age is regressed on BrainAGE prior to step 5, producing BrainAGE Residualized (BrainAGER) scores

EEG Microstates Temporal Dynamics Differentiate Individuals with Mood and Anxiety Disorders From Healthy Subjects

Electroencephalography (EEG) measures the brain’s electrophysiological spatio-temporal activities with high temporal resolution. Multichannel and broadband analysis of EEG signals is referred to as EEG microstates (EEG-ms) and can characterize such dynamic neuronal activity. EEG-ms have gained much attention due to the increasing evidence of their association with mental activities and large-scale brain networks identified by functional magnetic resonance imaging (fMRI). Spatially independent EEG-ms are quasi-stationary topographies (e.g., stable, lasting a few dozen milliseconds) typically classified into four canonical classes (microstates A through D). They can be identified by clustering EEG signals around EEG global field power (GFP) maxima points. We examined the EEG-ms properties and the dynamics of cohorts of mood and anxiety (MA) disorders subjects (n = 61) and healthy controls (HCs; n = 52). In both groups, we found four distinct classes of EEG-ms (A through D), which did not differ among cohorts. This suggests a lack of significant structural cortical abnormalities among cohorts, which would otherwise affect the EEG-ms topographies. However, both cohorts’ brain network dynamics significantly varied, as reflected in EEG-ms properties. Compared to HC, the MA cohort features a lower transition probability between EEG-ms B and D and higher transition probability from A to D and from B to C, with a trend towards significance in the average duration of microstate C. Furthermore, we harnessed a recently introduced theoretical approach to analyze the temporal dependencies in EEG-ms. The results revealed that the transition matrices of MA group exhibit higher symmetrical and stationarity properties as compared to HC ones. In addition, we found an elevation in the temporal dependencies among microstates, especially in microstate B for the MA group. The determined alteration in EEG-ms temporal dependencies among the cohorts suggests that brain abnormalities in mood and anxiety disorders reflect aberrant neural dynamics and a temporal dwelling among ceratin brain states (i.e., mood and anxiety disorders subjects have a less dynamicity in switching between different brain states)

Parameter Space and Potential for Biomarker Development in 25 Years of fMRI Drug Cue Reactivity: A Systematic Review.

IMPORTANCE: In the last 25 years, functional magnetic resonance imaging drug cue reactivity (FDCR) studies have characterized some core aspects in the neurobiology of drug addiction. However, no FDCR-derived biomarkers have been approved for treatment development or clinical adoption. Traversing this translational gap requires a systematic assessment of the FDCR literature evidence, its heterogeneity, and an evaluation of possible clinical uses of FDCR-derived biomarkers. OBJECTIVE: To summarize the state of the field of FDCR, assess their potential for biomarker development, and outline a clear process for biomarker qualification to guide future research and validation efforts. EVIDENCE REVIEW: The PubMed and Medline databases were searched for every original FDCR investigation published from database inception until December 2022. Collected data covered study design, participant characteristics, FDCR task design, and whether each study provided evidence that might potentially help develop susceptibility, diagnostic, response, prognostic, predictive, or severity biomarkers for 1 or more addictive disorders. FINDINGS: There were 415 FDCR studies published between 1998 and 2022. Most focused on nicotine (122 [29.6%]), alcohol (120 [29.2%]), or cocaine (46 [11.1%]), and most used visual cues (354 [85.3%]). Together, these studies recruited 19 311 participants, including 13 812 individuals with past or current substance use disorders. Most studies could potentially support biomarker development, including diagnostic (143 [32.7%]), treatment response (141 [32.3%]), severity (84 [19.2%]), prognostic (30 [6.9%]), predictive (25 [5.7%]), monitoring (12 [2.7%]), and susceptibility (2 [0.5%]) biomarkers. A total of 155 interventional studies used FDCR, mostly to investigate pharmacological (67 [43.2%]) or cognitive/behavioral (51 [32.9%]) interventions; 141 studies used FDCR as a response measure, of which 125 (88.7%) reported significant interventional FDCR alterations; and 25 studies used FDCR as an intervention outcome predictor, with 24 (96%) finding significant associations between FDCR markers and treatment outcomes. CONCLUSIONS AND RELEVANCE: Based on this systematic review and the proposed biomarker development framework, there is a pathway for the development and regulatory qualification of FDCR-based biomarkers of addiction and recovery. Further validation could support the use of FDCR-derived measures, potentially accelerating treatment development and improving diagnostic, prognostic, and predictive clinical judgments

Diagnosis of bipolar disorders and body mass index predict clustering based on similarities in cortical thickness-ENIGMA study in 2436 individuals

AIMS: Rates of obesity have reached epidemic proportions, especially among people with psychiatric disorders. While the effects of obesity on the brain are of major interest in medicine, they remain markedly under-researched in psychiatry. METHODS: We obtained body mass index (BMI) and magnetic resonance imaging-derived regional cortical thickness, surface area from 836 bipolar disorders (BD) and 1600 control individuals from 14 sites within the ENIGMA-BD Working Group. We identified regionally specific profiles of cortical thickness using K-means clustering and studied clinical characteristics associated with individual cortical profiles. RESULTS: We detected two clusters based on similarities among participants in cortical thickness. The lower thickness cluster (46.8% of the sample) showed thinner cortex, especially in the frontal and temporal lobes and was associated with diagnosis of BD, higher BMI, and older age. BD individuals in the low thickness cluster were more likely to have the diagnosis of bipolar disorder I and less likely to be treated with lithium. In contrast, clustering based on similarities in the cortical surface area was unrelated to BD or BMI and only tracked age and sex. CONCLUSIONS: We provide evidence that both BD and obesity are associated with similar alterations in cortical thickness, but not surface area. The fact that obesity increased the chance of having low cortical thickness could explain differences in cortical measures among people with BD. The thinner cortex in individuals with higher BMI, which was additive and similar to the BD-associated alterations, may suggest that treating obesity could lower the extent of cortical thinning in BD

Mega-analysis of association between obesity and cortical morphology in bipolar disorders:ENIGMA study in 2832 participants

Background: Obesity is highly prevalent and disabling, especially in individuals with severe mental illness including bipolar disorders (BD). The brain is a target organ for both obesity and BD. Yet, we do not understand how cortical brain alterations in BD and obesity interact. Methods: We obtained body mass index (BMI) and MRI-derived regional cortical thickness, surface area from 1231 BD and 1601 control individuals from 13 countries within the ENIGMA-BD Working Group. We jointly modeled the statistical effects of BD and BMI on brain structure using mixed effects and tested for interaction and mediation. We also investigated the impact of medications on the BMI-related associations. Results: BMI and BD additively impacted the structure of many of the same brain regions. Both BMI and BD were negatively associated with cortical thickness, but not surface area. In most regions the number of jointly used psychiatric medication classes remained associated with lower cortical thickness when controlling for BMI. In a single region, fusiform gyrus, about a third of the negative association between number of jointly used psychiatric medications and cortical thickness was mediated by association between the number of medications and higher BMI. Conclusions: We confirmed consistent associations between higher BMI and lower cortical thickness, but not surface area, across the cerebral mantle, in regions which were also associated with BD. Higher BMI in people with BD indicated more pronounced brain alterations. BMI is important for understanding the neuroanatomical changes in BD and the effects of psychiatric medications on the brain.</p

Principal component analysis as an efficient method for capturing multivariate brain signatures of complex disorders—ENIGMA study in people with bipolar disorders and obesity

Multivariate techniques better fit the anatomy of complex neuropsychiatric disorders which are characterized not by alterations in a single region, but rather by variations across distributed brain networks. Here, we used principal component analysis (PCA) to identify patterns of covariance across brain regions and relate them to clinical and demographic variables in a large generalizable dataset of individuals with bipolar disorders and controls. We then compared performance of PCA and clustering on identical sample to identify which methodology was better in capturing links between brain and clinical measures. Using data from the ENIGMA-BD working group, we investigated T1-weighted structural MRI data from 2436 participants with BD and healthy controls, and applied PCA to cortical thickness and surface area measures. We then studied the association of principal components with clinical and demographic variables using mixed regression models. We compared the PCA model with our prior clustering analyses of the same data and also tested it in a replication sample of 327 participants with BD or schizophrenia and healthy controls. The first principal component, which indexed a greater cortical thickness across all 68 cortical regions, was negatively associated with BD, BMI, antipsychotic medications, and age and was positively associated with Li treatment. PCA demonstrated superior goodness of fit to clustering when predicting diagnosis and BMI. Moreover, applying the PCA model to the replication sample yielded significant differences in cortical thickness between healthy controls and individuals with BD or schizophrenia. Cortical thickness in the same widespread regional network as determined by PCA was negatively associated with different clinical and demographic variables, including diagnosis, age, BMI, and treatment with antipsychotic medications or lithium. PCA outperformed clustering and provided an easy-to-use and interpret method to study multivariate associations between brain structure and system-level variables. Practitioner Points: In this study of 2770 Individuals, we confirmed that cortical thickness in widespread regional networks as determined by principal component analysis (PCA) was negatively associated with relevant clinical and demographic variables, including diagnosis, age, BMI, and treatment with antipsychotic medications or lithium. Significant associations of many different system-level variables with the same brain network suggest a lack of one-to-one mapping of individual clinical and demographic factors to specific patterns of brain changes. PCA outperformed clustering analysis in the same data set when predicting group or BMI, providing a superior method for studying multivariate associations between brain structure and system-level variables.</p

Principal component analysis as an efficient method for capturing multivariate brain signatures of complex disorders—ENIGMA study in people with bipolar disorders and obesity

Multivariate techniques better fit the anatomy of complex neuropsychiatric disorders which are characterized not by alterations in a single region, but rather by variations across distributed brain networks. Here, we used principal component analysis (PCA) to identify patterns of covariance across brain regions and relate them to clinical and demographic variables in a large generalizable dataset of individuals with bipolar disorders and controls. We then compared performance of PCA and clustering on identical sample to identify which methodology was better in capturing links between brain and clinical measures. Using data from the ENIGMA-BD working group, we investigated T1-weighted structural MRI data from 2436 participants with BD and healthy controls, and applied PCA to cortical thickness and surface area measures. We then studied the association of principal components with clinical and demographic variables using mixed regression models. We compared the PCA model with our prior clustering analyses of the same data and also tested it in a replication sample of 327 participants with BD or schizophrenia and healthy controls. The first principal component, which indexed a greater cortical thickness across all 68 cortical regions, was negatively associated with BD, BMI, antipsychotic medications, and age and was positively associated with Li treatment. PCA demonstrated superior goodness of fit to clustering when predicting diagnosis and BMI. Moreover, applying the PCA model to the replication sample yielded significant differences in cortical thickness between healthy controls and individuals with BD or schizophrenia. Cortical thickness in the same widespread regional network as determined by PCA was negatively associated with different clinical and demographic variables, including diagnosis, age, BMI, and treatment with antipsychotic medications or lithium. PCA outperformed clustering and provided an easy-to-use and interpret method to study multivariate associations between brain structure and system-level variables. Practitioner Points: In this study of 2770 Individuals, we confirmed that cortical thickness in widespread regional networks as determined by principal component analysis (PCA) was negatively associated with relevant clinical and demographic variables, including diagnosis, age, BMI, and treatment with antipsychotic medications or lithium. Significant associations of many different system-level variables with the same brain network suggest a lack of one-to-one mapping of individual clinical and demographic factors to specific patterns of brain changes. PCA outperformed clustering analysis in the same data set when predicting group or BMI, providing a superior method for studying multivariate associations between brain structure and system-level variables.</p