Anger is a common and debilitating symptom of post-traumatic stress disorder (PTSD). Although studies have identified brain circuits underlying anger experience and expression in healthy individuals, how these circuits interact with trauma remains unclear. Here, we performed the first study examining the neural correlates of anger in patients with PTSD. Using a data-driven approach with resting-state fMRI, we identified two prefrontal regions whose overall functional connectivity was inversely associated with anger: the left anterior middle frontal gyrus (aMFG) and the right orbitofrontal cortex (OFC). We then used concurrent TMS-EEG to target the left aMFG parcel previously identified through fMRI, measuring its cortical excitability and causal connectivity to downstream areas. We found that low-anger PTSD patients exhibited enhanced excitability in the left aMFG and enhanced causal connectivity between this region and visual areas. Together, our results suggest that left aMFG activity may confer protection against the development of anger, and therefore may be an intriguing target for circuit-based interventions for anger in PTSD.

Repetitive transcranial magnetic stimulation (rTMS) is a commonly- used treatment for major depressive disorder (MDD). However, our understanding of the mechanism by which TMS exerts its antidepressant effect is minimal. Furthermore, we lack brain signals that can be used to predict and track clinical outcome. Such signals would allow for treatment stratification and optimization. Here, we performed a randomized, sham-controlled clinical trial and measured electrophysiological, neuroimaging, and clinical changes before and after rTMS. Patients (N = 36) were randomized to receive either active or sham rTMS to the left dorsolateral prefrontal cortex (dlPFC) for 20 consecutive weekdays. To capture the rTMS-driven changes in connectivity and causal excitability, resting fMRI and TMS/EEG were performed before and after the treatment. Baseline causal connectivity differences between depressed patients and healthy controls were also evaluated with concurrent TMS/fMRI. We found that active, but not sham rTMS elicited (1) an increase in dlPFC global connectivity, (2) induction of negative dlPFC-amygdala connectivity, and (3) local and distributed changes in TMS/EEG potentials. Global connectivity changes predicted clinical outcome, while both global connectivity and TMS/EEG changes tracked clinical outcome. In patients but not healthy participants, we observed a perturbed inhibitory effect of the dlPFC on the amygdala. Taken together, rTMS induced lasting connectivity and excitability changes from the site of stimulation, such that after active treatment, the dlPFC appeared better able to engage in top-down control of the amygdala. These measures of network functioning both predicted and tracked clinical outcome, potentially opening the door to treatment optimization.

Predictive models have succeeded in distinguishing between individuals with Alcohol use Disorder (AUD) and controls. However, predictive models identifying who is prone to develop AUD and the biomarkers indicating a predisposition to AUD are still unclear. Our sample (n = 656) included offspring and non-offspring of European American (EA) and African American (AA) ancestry from the Collaborative Study of the Genetics of Alcoholism (COGA) who were recruited as early as age 12 and were unaffected at first assessment and reassessed years later as AUD (DSM-5) (n = 328) or unaffected (n = 328). Machine learning analysis was performed for 220 EEG measures, 149 alcohol-related single nucleotide polymorphisms (SNPs) from a recent large Genome-wide Association Study (GWAS) of alcohol use/misuse and two family history (mother DSM-5 AUD and father DSM-5 AUD) features using supervised, Linear Support Vector Machine (SVM) classifier to test which features assessed before developing AUD predict those who go on to develop AUD. Age, gender, and ancestry stratified analyses were performed. Results indicate significant and higher accuracy rates for the AA compared with the EA prediction models and a higher model accuracy trend among females compared with males for both ancestries. Combined EEG and SNP features model outperformed models based on only EEG features or only SNP features for both EA and AA samples. This multidimensional superiority was confirmed in a follow-up analysis in the AA age groups (12-15, 16-19, 20-30) and EA age group (16-19). In both ancestry samples, the youngest age group achieved higher accuracy score than the two other older age groups. Maternal AUD increased the model's accuracy in both ancestries' samples. Several discriminative EEG measures and SNPs features were identified, including lower posterior gamma, higher slow wave connectivity (delta, theta, alpha), higher frontal gamma ratio, higher beta correlation in the parietal area, and 5 SNPs: rs4780836, rs2605140, rs11690265, rs692854, and rs13380649. Results highlight the significance of sampling uniformity followed by stratified (e.g., ancestry, gender, developmental period) analysis, and wider selection of features, to generate better prediction scores allowing a more accurate estimation of AUD development.

The understanding and treatment of psychiatric disorders, which are known to be neurobiologically and clinically heterogeneous, could benefit from the data-driven identification of disease subtypes. Here, we report the identification of two clinically relevant subtypes of post-traumatic stress disorder (PTSD) and major depressive disorder (MDD) on the basis of robust and distinct functional connectivity patterns, prominently within the frontoparietal control network and the default mode network. We identified the disease subtypes by analysing, via unsupervised and supervised machine learning, the power-envelope-based connectivity of signals reconstructed from high-density resting-state electroencephalography in four datasets of patients with PTSD and MDD, and show that the subtypes are transferable across independent datasets recorded under different conditions. The subtype whose functional connectivity differed most from those of healthy controls was less responsive to psychotherapy treatment for PTSD and failed to respond to an antidepressant medication for MDD. By contrast, both subtypes responded equally well to two different forms of repetitive transcranial magnetic stimulation therapy for MDD. Our data-driven approach may constitute a generalizable solution for connectome-based diagnosis.

A mechanistic understanding of the pathology of psychiatric disorders has been hampered by extensive heterogeneity in biology, symptoms, and behavior within diagnostic categories that are defined subjectively. We investigated whether leveraging individual differences in information-processing impairments in patients with post-traumatic stress disorder (PTSD) could reveal phenotypes within the disorder. We found that a subgroup of patients with PTSD from two independent cohorts displayed both aberrant functional connectivity within the ventral attention network (VAN) as revealed by functional magnetic resonance imaging (fMRI) neuroimaging and impaired verbal memory on a word list learning task. This combined phenotype was not associated with differences in symptoms or comorbidities, but nonetheless could be used to predict a poor response to psychotherapy, the best-validated treatment for PTSD. Using concurrent focal noninvasive transcranial magnetic stimulation and electroencephalography, we then identified alterations in neural signal flow in the VAN that were evoked by direct stimulation of that network. These alterations were associated with individual differences in functional fMRI connectivity within the VAN. Our findings define specific neurobiological mechanisms in a subgroup of patients with PTSD that could contribute to the poor response to psychotherapy.

Results of neuroimaging datasets aggregated from multiple sites may be biased by site-specific profiles in participants' demographic and clinical characteristics, as well as MRI acquisition protocols and scanning platforms. We compared the impact of four different harmonization methods on results obtained from analyses of cortical thickness data: (1) linear mixed-effects model (LME) that models site-specific random intercepts (LME_INT), (2) LME that models both site-specific random intercepts and age-related random slopes (LME_INT+SLP), (3) ComBat, and (4) ComBat with a generalized additive model (ComBat-GAM). Our test case for comparing harmonization methods was cortical thickness data aggregated from 29 sites, which included 1,340 cases with posttraumatic stress disorder (PTSD) (6.2-81.8 years old) and 2,057 trauma-exposed controls without PTSD (6.3-85.2 years old). We found that, compared to the other data harmonization methods, data processed with ComBat-GAM was more sensitive to the detection of significant case-control differences (Χ²(3) = 63.704, p < 0.001) as well as case-control differences in age-related cortical thinning (Χ²(3) = 12.082, p = 0.007). Both ComBat and ComBat-GAM outperformed LME methods in detecting sex differences (Χ²(3) = 9.114, p = 0.028) in regional cortical thickness. ComBat-GAM also led to stronger estimates of age-related declines in cortical thickness (corrected p-values < 0.001), stronger estimates of case-related cortical thickness reduction (corrected p-values < 0.001), weaker estimates of age-related declines in cortical thickness in cases than controls (corrected p-values < 0.001), stronger estimates of cortical thickness reduction in females than males (corrected p-values < 0.001), and stronger estimates of cortical thickness reduction in females relative to males in cases than controls (corrected p-values < 0.001). Our results support the use of ComBat-GAM to minimize confounds and increase statistical power when harmonizing data with non-linear effects, and the use of either ComBat or ComBat-GAM for harmonizing data with linear effects.

BACKGROUND: Recent advances in data-driven computational approaches have been helpful in devising tools to objectively diagnose psychiatric disorders. However, current machine learning studies limited to small homogeneous samples, different methodologies, and different imaging collection protocols, limit the ability to directly compare and generalize their results. Here we aimed to classify individuals with PTSD versus controls and assess the generalizability using a large heterogeneous brain datasets from the ENIGMA-PGC PTSD Working group. METHODS: We analyzed brain MRI data from 3,477 structural-MRI; 2,495 resting state-fMRI; and 1,952 diffusion-MRI. First, we identified the brain features that best distinguish individuals with PTSD from controls using traditional machine learning methods. Second, we assessed the utility of the denoising variational autoencoder (DVAE) and evaluated its classification performance. Third, we assessed the generalizability and reproducibility of both models using leave-one-site-out cross-validation procedure for each modality. RESULTS: We found lower performance in classifying PTSD vs. controls with data from over 20 sites (60 % test AUC for s-MRI, 59 % for rs-fMRI and 56 % for d-MRI), as compared to other studies run on single-site data. The performance increased when classifying PTSD from HC without trauma history in each modality (75 % AUC). The classification performance remained intact when applying the DVAE framework, which reduced the number of features. Finally, we found that the DVAE framework achieved better generalization to unseen datasets compared with the traditional machine learning frameworks, albeit performance was slightly above chance. CONCLUSION: These results have the potential to provide a baseline classification performance for PTSD when using large scale neuroimaging datasets. Our findings show that the control group used can heavily affect classification performance. The DVAE framework provided better generalizability for the multi-site data. This may be more significant in clinical practice since the neuroimaging-based diagnostic DVAE classification models are much less site-specific, rendering them more generalizable.

Although the cerebellum contributes to higher-order cognitive and emotional functions relevant to posttraumatic stress disorder (PTSD), prior research on cerebellar volume in PTSD is scant, particularly when considering subregions that differentially map on to motor, cognitive, and affective functions. In a sample of 4215 adults (PTSD n = 1642; Control n = 2573) across 40 sites from the ENIGMA-PGC PTSD working group, we employed a new state-of-the-art deep-learning based approach for automatic cerebellar parcellation to obtain volumetric estimates for the total cerebellum and 28 subregions. Linear mixed effects models controlling for age, gender, intracranial volume, and site were used to compare cerebellum volumes in PTSD compared to healthy controls (88% trauma-exposed). PTSD was associated with significant grey and white matter reductions of the cerebellum. Compared to controls, people with PTSD demonstrated smaller total cerebellum volume, as well as reduced volume in subregions primarily within the posterior lobe (lobule VIIB, crus II), vermis (VI, VIII), flocculonodular lobe (lobule X), and corpus medullare (all p_-FDR < 0.05). Effects of PTSD on volume were consistent, and generally more robust, when examining symptom severity rather than diagnostic status. These findings implicate regionally specific cerebellar volumetric differences in the pathophysiology of PTSD. The cerebellum appears to play an important role in higher-order cognitive and emotional processes, far beyond its historical association with vestibulomotor function. Further examination of the cerebellum in trauma-related psychopathology will help to clarify how cerebellar structure and function may disrupt cognitive and affective processes at the center of translational models for PTSD.