Trial PaperDepressive DisordersNeuroimaging & Brain MeasuresKetamine

High-order brain interactions during ketamine-induced state changes: A functional marker of response in late-life treatment-resistant depression?

This secondary analysis of a randomised, double-blind, midazolam-controlled trial (n=30) examined EEG brain interaction patterns in late-life veterans with treatment-resistant depression after ketamine infusion. Ketamine caused time-dependent changes in these higher-order interactions, and larger 24-hour increases in alpha-band redundancy were linked to greater improvement in depressive symptoms by day 7.

1 linked clinical trial·4 references indexed in Blossom

Authors

  • Sanjay Mathew
  • Brittany O'Brien
  • Alan Craig Swann

Published

Translational Psychiatry
individual Study

Abstract

Ketamine is a fast-acting intervention for treatment-resistant depression (TRD), yet only a subset of patients show robust clinical response, and the underlying neural mechanisms remain unclear. High-order interactions (HOI) derived from multivariate information theory provide a framework for examining nonlinear dependencies among brain regions beyond pairwise connectivity. One such metric, the O-information, captures the balance between synergistic and redundant interactions across three or more variables. In this secondary analysis of a randomized, double-blind, midazolam-controlled trial (NCT02556606), we examined EEG-derived HOI in 30 late-life veterans with TRD following a single 40-minute intravenous infusion of ketamine (0.1, 0.25, 0.5 mg/kg; n = 18) or midazolam (0.03 mg/kg; n = 12). Resting state and mismatch negativity data were analyzed at baseline, 1 h, 24 h, and 7 d post-infusion. Ketamine induced temporally dynamic alterations in redundancy-dominant O-info, with maximal effects in the alpha-band at 1 h (Cohen’s d = 2.57), attenuation at 24 h that shifted toward the theta-band, and partial resurgence in beta and gamma by Day 7. Linear mixed-effects modeling identified significant group effects across most band x metric families, with the strongest effects in alpha, beta, and gamma redundancy. Greater increases in 24-hour alpha-band redundancy were associated with greater improvement in depressive symptoms at Day 7 (β = 69.31, q = 0.05). HOI metrics also tracked acute dissociative states, with several 24-hour alpha and beta features remaining positively associated with symptom severity after correction. These findings extend prior HOI work in healthy samples to a controlled TRD cohort and suggest that ketamine induces temporally structured reorganization of higher-order brain interactions, with exploratory associations to clinical outcomes.

Unlocked with Blossom Pro

Research Summary of 'High-order brain interactions during ketamine-induced state changes: A functional marker of response in late-life treatment-resistant depression?'

Editorial

βBlossom's Take

This study matters because it pushes ketamine research beyond power spectra and pairwise connectivity, asking whether higher-order EEG structure tracks state change and later antidepressant response in a controlled late-life depression sample. The association with 24-hour alpha redundancy is intriguing, but the signal is narrow and exploratory, so the main value is as a biomarker hypothesis rather than a settled clinical marker.

MADRS over time

Sourced

Per-arm mean scores at each reported timepoint. The experimental arm is marked in orange; comparators in grey.

-1.212.826.840.8Day 0Day 7Mean scoreKetamine 0.10 mg/kg · Day 0 · 33.9 · n=18Ketamine 0.10 mg/kg · Day 7 · 15.9 · n=15Ketamine 0.25 mg/kg · Day 0 · 33.9 · n=18Ketamine 0.25 mg/kg · Day 7 · 15.9 · n=15Ketamine 0.50 mg/kg · Day 0 · 33.9 · n=18Ketamine 0.50 mg/kg · Day 7 · 15.9 · n=15Midazolam 0.03 mg/kg · Day 0 · 35 · n=12Midazolam 0.03 mg/kg · Day 7 · 15.6 · n=9
Ketamine 0.10 mg/kgn=18
Ketamine 0.25 mg/kgn=18
Ketamine 0.50 mg/kgn=18
Midazolam 0.03 mg/kgn=12

Extracted summary values for one outcome measure (MADRS); see the Results tab for all outcomes and comparisons.

Introduction

Ketamine is recognised as a rapid-acting option for treatment-resistant depression, but only a subset of patients show a robust clinical response, and the neural mechanisms that support benefit remain uncertain. The paper places this problem in the context of prior EEG findings suggesting that ketamine alters oscillatory activity and may shift the brain towards a more flexible functional state. The authors argue that conventional pairwise connectivity measures may miss more complex, nonlinear relations between brain regions, and that high-order interactions (HOI) derived from multivariate information theory could better capture these network changes. Shah and colleagues set out to test whether ketamine produces temporally evolving changes in HOI in late-life veterans with treatment-resistant depression, and whether those changes relate to antidepressant response and acute dissociation. They specifically hypothesised that ketamine-induced changes would be associated with increases in brain redundancy, as measured by the O-information, and that these effects would differ from those seen with midazolam. The study is presented as a secondary analysis intended both to extend earlier HOI work in healthy volunteers and to explore whether HOI might serve as a clinically useful biomarker in depression treatment.

Methods

This was a secondary analysis of EEG data from a randomised, double-blind, midazolam-controlled, multi-arm trial of intravenous ketamine in late-life treatment-resistant depression (ClinicalTrials.gov NCT02556606). The analysed sample comprised 30 veterans: 18 received ketamine and 12 received midazolam. The ketamine arm pooled three doses, 0.1, 0.25, and 0.5 mg/kg, on the basis of prior work suggesting dose-invariant EEG effects; the authors nevertheless ran dose-sensitivity analyses to check whether pooling was reasonable. Missing data were excluded. EEG was recorded from a 64-channel system. For resting-state analysis, alternating eyes-open and eyes-closed segments were collected, but the analysis was restricted to the first 100 seconds to focus on the eyes-closed condition and reduce visual variability. The data were filtered, downsampled, cleaned using artefact subspace reconstruction and independent component analysis, and then analysed on a reduced 15-electrode montage because the original HOI method required a comparable electrode set. For mismatch negativity (MMN), data were epoched around stimulus onset, cleaned similarly, and then concatenated into continuous time series for standard and deviant conditions. The main outcomes were high-order information-theoretic metrics: O-information, which reflects whether interactions are more redundancy-dominant or synergy-dominant, and S-information, which reflects overall higher-order interdependence. Analyses focused on baseline, 1 hour, 24 hours, and 7 days after infusion for resting-state data, and 2 hours and 7 days for MMN data. Post-infusion values were expressed as change from baseline. The authors first used effect-size-based feature selection to identify the strongest electrode combinations and then evaluated those features in linear mixed-effects models with subject-level random intercepts. Model selection used Bayesian Information Criterion and considered age, sex, BMI, baseline HOI, group, time, and group-by-time interactions. Associations with clinical change were then examined using Spearman correlations, permutation testing, false discovery rate correction, and ordinary least squares models adjusted for age, sex, BMI, and baseline HOI. Clinical outcomes were MADRS change from baseline to day 7 and CADSS change from baseline to 1 hour.

Results

In resting-state EEG, ketamine produced temporally dynamic changes in HOI, with the largest single effect seen for alpha-band O-information at 1 hour after infusion (Cohen's d = 2.57). Across time, the effect pattern shifted: alpha effects were strongest early, theta effects were more evident at 24 hours, and gamma effects partially re-emerged by day 7. Across all 32,647 analysed n-plets, distribution-level summaries showed a similar early positive shift, while S-information tended to become negative by day 7, particularly in alpha and theta bands. In the mixed-effects models, the best-fitting structure included group, timepoint, baseline HOI, and a group-by-timepoint interaction. Significant group effects survived false discovery rate correction in 7 of 8 band-by-metric families; theta S-information was the only family that did not survive correction. The strongest effects were for alpha O-information (χ² = 48.73, q < .001), beta O-information (χ² = 24.70, q < .001), and gamma O-information (χ² = 23.41, q < .001). Baseline HOI was significant across all band-by-metric families. Significant group-by-time effects were present for alpha O-information (χ² = 25.50, q < .001) and gamma O-information (χ² = 10.27, q = .024). The maximal electrode configurations were not stable across timepoints, suggesting that different n-plets carried the peak effects at different stages. During MMN processing, ketamine also altered HOI, with early increases followed by attenuation. Mixed-effects modelling of the maximal n-plets supported a model including group, timepoint, and baseline HOI. FDR-corrected group effects remained significant for O-information in both deviant and standard responses (q = .001), whereas S-information did not survive correction. The strongest task-evoked effects shifted spatially from central-posterior sites early to more frontal-central configurations later. Clinical associations were selective. For depressive symptoms, only one feature remained significant after correction in the covariate-adjusted model: 24-hour alpha O-information predicted day 7 MADRS improvement (β = 69.31, SE = 25.77, p = .007, q = .050). This meant that greater increases in alpha-band redundancy at 24 hours were associated with greater symptom improvement at day 7. Other candidate associations were present but did not survive correction. Associations with dissociative symptoms were stronger. In the primary model, four 24-hour features survived FDR correction: alpha O-information (β = 6.10, q = .007), alpha S-information (β = 5.48, q = .019), beta O-information (β = 34.11, q < .001), and beta S-information (β = 6.00, q = .002). All were positive, indicating that greater increases in HOI were associated with greater CADSS scores. The authors note that some earlier candidate associations, including gamma-band features, did not survive correction in the final models.

Discussion

Shah and colleagues interpret the findings as evidence that ketamine induces large-scale, frequency-specific reorganisation of higher-order brain dynamics in late-life treatment-resistant depression. They emphasise that the effects were temporally structured rather than static: alpha-band redundancy rose earliest and most strongly, theta effects were more apparent later, and gamma-related effects persisted to day 7. In their view, this pattern supports the idea that ketamine changes not just the strength of brain activity, but the way information is distributed across the network over time. They also highlight that these effects were detectable with a relatively low-density electrode montage, which they present as encouraging for potential clinical scalability. The authors relate their results to earlier HOI work in healthy volunteers, stating that the present study extends those findings into a controlled clinical sample. They argue that the observed increase in alpha-band redundancy is notable because it does not simply mirror conventional oscillatory power changes; instead, it may reflect an aspect of network organisation that is not captured by spectral measures alone. They suggest that this could be relevant to models in which ketamine reduces hierarchical control and increases flexibility in information processing. For clinical outcomes, the authors take a cautious view. They report that the antidepressant association was limited to a single 24-hour alpha O-information feature, and they describe this as selective and exploratory rather than a broad HOI signal. By contrast, associations with acute dissociation were stronger and more consistent, with multiple 24-hour alpha and beta features surviving correction. They interpret this as suggesting that HOI-clinical relationships differ by outcome and are not reducible to one shared mechanism within this dataset. The authors acknowledge several limitations. The single-infusion design limits generalisability to repeated dosing, and the sample was small, late-life, predominantly male, and drawn from veterans, which reduces representativeness. The modest ketamine subgroup size limited power to detect dose heterogeneity, and the small midazolam arm together with adaptive randomisation may have created imbalance. The reduced electrode montage improves scalability but lowers spatial resolution and may miss finer network structure. They also note that HOI metrics are undirected, so they cannot determine the direction of information flow, and that some spatial differences from earlier studies suggest variability in exact electrode contributions. They further state that future work should assess temporal trajectories more explicitly, integrate HOI with spectral power, refine spatial modelling, and test predictive utility in larger longitudinal samples. The authors conclude that HOI may be a useful and potentially scalable framework for studying ketamine-related changes in brain organisation, while still requiring further validation.

Expert Research Summary

Go Pro to unlock section-by-section summaries, Blossom's editorial take, and the complete research toolkit.

See Pro Plans
View full paper sections

A R T I C L E I N P R E S S ARTICLE IN PRESS

One such approach, high-order interactions (HOI), uses generalizations of mutual information to characterize multivariate dependencies among electrode combinations, as well as the dominant quality of interaction within a network. For the former, the S-information measures the overall level of interdependencies, while for the latter, the O-information measures the dominance between synergy and redundancy. Synergy refers to relationships among variables in which the whole system conveys more information than the sum of its parts, while redundancy refers to relationships in which copies of the same information are distributed throughout the whole system. In previous work, HOI revealed increased alpha-band redundancy during rest, and increased redundancy in response to standard tones during the auditory mismatch negativity (MMN). Because redundancy reflects widespread broadcasting of shared information across the brain, this prior work suggested that ketamine might reduce top-down control, thereby increasing the influence of associative and sensory cortical systems on brain hierarchy. Notably, the effects of ketamine on resting alpha were also associated with an increased sense of derealization on the Clinician Administered Dissociative States Scale (CADSS) in healthy volunteers. These findings are relevant not only to models of consciousness, where altered redundancy may reflect changes in global information integration, but also to depression, which involves disruptions in how individuals process internal thoughts and external sensory information. To our knowledge, no prior study has evaluated HOI as a clinical biomarker of ketamine response in MDD. Here, we hypothesize that ketamine-induced changes in brain dynamics and depressive symptoms are associated with increases in brain redundancy. HOI provides a powerful framework for characterizing the spatiotemporal organization of functional brain networks critical to therapeutic response, positioning these metrics as promising targets for precision treatment. We present a secondary analysis of EEG data collected during the periinfusion period of a single ketamine or active placebo (midazolam) infusion in late-life veterans

S ARTICLE IN PRESS

with TRD. Specifically, we examined HOI among EEG channels during resting and cognitively engaged states to delineate the cortical information dynamics associated with clinical improvement following a single infusion of intravenous (IV) ketamine versus midazolam. Our study design builds upon previous researchto develop a dynamic framework aimed at advancing the precision treatment of MDD with ketamine, while also seeking to externally validate prior HOI findings in healthy volunteers in an independent clinical sample.

DESIGN AND OVERVIEW

This secondary analysis used EEG data from a randomized, midazolam-controlled, double-blind, multi-arm trial investigating dose-dependent ketamine effects in TRD [ClinicalTrials.gov: NCT02556606]. Here, we aimed to characterize temporally evolving changes in higher-order brain interactions following ketamine infusion. Detailed trial methods are described in. Based on prior work demonstrating dose-invariant EEG metrics, HOI metrics were examined using a dose-inclusive ketamine group (0.1, 0.25, 0.5 mg/kg) and compared with midazolam (0.03 mg/kg). The ketamine group included three participants at 0.1 mg/kg, five at 0.25 mg/kg, and ten at 0.5 mg/kg. Formal dose-sensitivity analyses were subsequently performed to evaluate whether pooling across ketamine doses was appropriate. Analyses focused on EEG data collected at baseline, 1 hour, 24 hours, and 7 days postinfusion, corresponding to timepoints of primary clinical interest. HOI metrics included Oinformation (O-info) and S-information (S-info). O-info quantifies the balance between redundancy-and synergy-dominated interactions, whereas S-info reflects the overall magnitude of higher-order interdependencies. Missing data were excluded from the analyses. To characterize treatment-related change, post-infusion HOI values were expressed as change scores relative to baseline (ΔHOI). infusion), 1h, 2h, 4h, 24h, and 7 days, while clinical scales were repeated at 40 min, 2h, 4h, 24h, 48h, 72h, and 7 days. Vital signs were monitored every 15 minutes for 4 hours. An overview of the randomization schema and assessment timepoints is provided in Fig.and presented in detail in previous literature.

EEG ACQUISITION AND PREPROCESSING

EEG was recorded using a 64-channel Curry 7 system (SynAmps2 amplifier) with an analog-to-digital conversion rate of 1000 Hz. Resting-state recordings consisted of alternating 2minute eyes-open and eyes-closed segments (total duration: 4 minutes), with eyes-closed always recorded first and no counterbalancing across visits or participants. Both conditions were retained during preprocessing to improve artifact identification during Independent Component Analysis (ICA). After artifact removal, analyses were restricted to the first 100 seconds of data to isolate the eyes-closed condition and minimize visual inputrelated variability. Continuous data were band-pass filtered (1-50 Hz, finite impulse response) and downsampled to 250 Hz. Line noise (59-61 Hz) was regressed out using the CleanLine plugin. We removed artifacts using a two-step approach to differentiate spontaneous mechanical artifacts from biological artifacts (electrocardiogram/electrooculogram/blinks). First, artifact subspace reconstruction (ASR 23 ) was applied to identify bad channels and recover cortical activity obscured by bursts of sporadic artifacts (channel criterion r = 0.75; burst threshold: 25 SD). Spherical interpolation was used to reconstruct excluded channels, and the dataset was rereferenced to the grand average. Second, ICA was used to identify and remove components with spatial, temporal, and spectral properties indicative of ocular and cardiac activity. The FastICA algorithm 24 was applied using the simultaneous decomposition method and dimension reduction to match the dataset rank. ICLabel 25 was used to remove components with a 90% probability of being artifacts. MMN data were epoched from -100 to 500 ms relative to stimulus onset and preprocessed similarly, with additional probabilistic rejection of trials exceeding ± 5 SD. Stimulus parameters and trial counts followed prior work. Epoched MMN data were concatenated into whole-band continuous time series for standard and deviant conditions. During preprocessing, an average of 6 seconds (±18 seconds) of unrecoverable burst activity was removed following ASR. ICLabel rejected an average of 8.9 components (±5.7). For MMN data, an average of 36.2 (±66.3) standard trials and 3.5 (±7.2) deviant trials were excluded. These preprocessing steps were implemented to balance artifact removal with retention of sufficient data for reliable HOI estimation.

HIGH-ORDER INTERACTIONS METRICS

We applied previously established high-order information-theoretic metricsto assess nonlinear, multivariate EEG interactions. The original methodology utilized 16 electrodes; however, our montage did not include FCz or a suitable substitute due to the placement of the ground and reference electrodes. To maintain methodological consistency with prior work, analyses were performed on a reduced 15-electrode subset derived from a 64-channel recording. From these 15 electrodes, we computed all the possible combinations of interacting electrodes (n-plets), ranging from 2-plets (pairwise) to a single 15-plet (global interaction). Each n-plet, representing a unique combination of n electrodes, was analyzed using four entropybased metrics. The calculus of HOI is described in detail in previous work. We present a paraphrased version of this description below for ease of access. Let 𝑋 𝑛 = (𝑋 1 , …, 𝑋 𝑛 ) be a set of n electrodes, where 𝑋 1 , 𝑋 2 , and 𝑋 𝑛 correspond to the time series of electrodes 1, 2, and n, where 𝐻(𝑋 1 , . . . , 𝑋 𝑛 ) is the joint Shannon's entropy of the n electrodes, 𝐻(𝑋 𝑖 ) the entropy of the i-th electrode, and 𝐻(𝑋 𝑖 |𝑋 𝑛 -𝑖 ) is the entropy of the i-th electrode conditioned by the activity of all the remaining electrodes, or "residual entropy"𝑋 𝑖 . Estimations were performed using the Gaussian copula approximationusing the THOI package 28 . This approach enabled comprehensive analysis of both pairwise and high-order connectivity changes across critical post-infusion timepoints relative to baseline (1h, 24h, 7d). To maintain statistical power, electrode and frequency-band selection were optimized using effect size-based feature selection (see below). The reduced montage ensured both methodological continuity and scalability for the longitudinal clinical dataset.

DOSE SENSITIVITY ANALYSIS

To assess whether HOI features varied as a function of ketamine dose (0.1, 0.25, 0.5 mg/kg), we conducted ketamine-only analyses at each timepoint (1h, 24h, Day 7) for resting state and mismatch negativity data. For each selected feature, linear models of the form: were fit, with dose treated as a continuous predictor. Given the modest sample size (n = 18), nonparametric bootstrap confidence intervals (5,000 resamples) were computed for each slope, and empirical p-values were derived. False Discovery Rate (FDR) correction (Benjamini-Hochberg) was applied within each timepoint across all tested features. In the absence of significant dose effects, analyses were conducted on pooled data. Because n-plets were selected in a data-driven manner to maximize observed effects within a modest sample, subsequent analyses were interpreted as exploratory.

RESTING STATE (RS):

To quantify treatment-related group differences, we computed standardized effect sizes (Cohen's d for independent groups) between the ketamine and midazolam groups at each postinfusion timepoint (1h, 24h, Day 7). For each HOI metric and frequency band (theta, alpha, beta, gamma), per n-plet effect sizes were calculated as: where 𝑋 ̅ 𝑘𝑒𝑡𝑎𝑚𝑖𝑛𝑒 and 𝑋 ̅ 𝑚𝑖𝑑𝑎𝑧𝑜𝑙𝑎𝑚 represent the mean HOI value of a given n-plet across participants within each group, and 𝑆𝐷 𝑝𝑜𝑜𝑙𝑒𝑑 reflects pooled variance across groups. Missing, zero, and infinite values were excluded. Effect sizes were ranked within each band x metric x timepoint, and the n-plet with the maximal absolute effect size was selected. Corresponding participant-level HOI values and electrode configurations were extracted for downstream analyses.

MISMATCH NEGATIVITY (MMN):

An analogous procedure was applied to MMN-derived time series. HOI metrics were computed separately for deviant and standard continuous data at 2 hours and Day 7. After exclusion of invalid values, maximal effect size n-plets were identified and extracted at the participant level. This approach isolated condition-specific effects while maintaining consistency with resting-state analyses.

FEATURE SELECTION AND MODEL-BASED EVALUATION

Selected maximal n-plets from the effect size analysis were carried forward for modelbased evaluation. For each timepoint, the n-plet exhibiting the largest absolute effect size within each band x metric combination was identified. To evaluate robustness and account for covariates, we implemented a two-stage modeling approach. First, candidate fixed-effects structures were compared using backward elimination and exhaustive model selection based on the Bayesian Information Criterion (BIC), considering age, sex, BMI, baseline HOI, group, time, and group x time interactions. The optimal fixed-effects structure was then incorporated into a linear mixed-effects model, with subject-level random intercepts to account for repeated measures. To assess temporal stability, n-plets with the highest effect size at each timepoint were tracked across all three timepoints. We additionally characterized distributions across bands and metrics using summary statistics (mean, median, proportion of positive vs negative effects) and visualized temporal stability using violin plots.

ASSOCIATIONS BETWEEN HOI METRICS AND CLINICAL OUTCOMES

To examine associations between HOI metrics and clinical outcomes (MADRS and CADSS), we applied a data-driven correlation framework designed to identify candidate HOI features associated with clinical change independent of the effect-size-based feature selection procedure. For each band x metric x timepoint, Spearman correlations were computed between ΔHOI and changes in clinical scores (ΔMADRS: baseline to Day 7; ΔCADSS: baseline to 1h). Statistical significance thresholds were derived using permutation testing (10,000 iterations), generating null distributions via random shuffling of clinical outcomes. Two-tailed significant thresholds were defined at the 2.5th and 97.5th percentiles of the permutation distribution, preserving both positive and negative effects. Surviving correlations were subjected to false discovery rate (FDR) correction using the Benjamini-Hochberg procedure (p < 0.05). For each timepoint, n-plets exhibiting the strongest positive and negative FDR-surviving associations were identified and extracted at the participant level.

MODEL-BASED INFERENCE OF CLINICAL ASSOCIATIONS

To evaluate whether identified HOI features independently predicted clinical outcomes, we fit ordinary least squares (OLS) models with ΔHOI as the primary predictor and age, sex, BMI, and baseline HOI as covariates. Although ΔHOI captures within-subject changes, individual variability in baseline network architecture, such as baseline synergy or redundancy, may influence its relationship with clinical outcomes. Prior work suggests that baseline electrophysiological features, including gamma, theta, and alpha activity, may moderate treatment response to ketamine, supporting adjustment for baseline HOI. Sensitivity analyses excluding baseline HOI were conducted to assess potential collinearity effects.

KETAMINE ADMINISTRATION INDUCED TEMPORALLY DYNAMIC ALTERATIONS IN RESTING-STATE HIGHORDER INTERACTIONS.

We first examined ketamine-related changes in resting-state HOI using effect-size-based feature selection followed by covariate-adjusted mixed-effects modeling. The largest single effect was observed for alpha-band O-info at 1 hour (d = 2.57) (Fig.). Across timepoints, maximal effects followed a structured temporal progression, with alpha O-info peaking at 1h, theta O-info at 24h, and gamma O-info at Day 7, corresponding to peak separation at 1h, attenuation at 24h, and partial re-emergence by Day 7 (Fig.). Distribution-level analyses across all 32,647 n-plets supported this pattern (Fig.positive). S-info showed a similar early positive shift but was predominantly negative by Day 7, particularly in alpha and theta (median ds = -0.45 and -0.48, respectively) (Supplement Fig.

S2).

Linear mixed-effects modeling showed that the best-fitting model (BIC-based) included group, timepoint, baseline HOI, and group x timepoint interaction. Significant group effects survived FDR correction in 7 of 8 band × metric families, with theta S-info not surviving correction (q = .078). Effects were strongest for alpha O-info (χ² = 48.73, q < .001), beta O-info (χ² = 24.70, q < .001), and gamma O-info (χ² = 23.41, q < .001). Baseline HOI was a significant predictor across all band × metric families after FDR correction. Significant group × timepoint effects were present for alpha O-info (χ² = 25.50, q < .001) and gamma O-info (χ² = 10.27, q = .024), indicating temporally evolving treatment effects (Supplement Table). Temporal tracking of the maximal n-plets showed limited stability across timepoints, with peak configurations often attenuating or changing signs outside their defining timepoint (Fig.). Spatially, peak effects shifted from midline regions at 1h to more distributed frontoparietal and posterior configurations at later timepoints. Together, these findings suggest a dynamic reorganization of resting-state HOI after ketamine, with the clearest and most time-sensitive effects in O-info.

KETAMINE INDUCES CHANGES IN REDUNDANCY DURING MMN.

Second, we examined ketamine-related changes in HOI during mismatch negativity processing using effect-size-based feature selection followed by mixed-effects validation. Distribution-level analyses across all 32,647 n-plets showed a clear temporal shift (Supplement Linear mixed-effects modeling of maximal n-plet features supported a BIC-selected model including group, timepoint, and baseline HOI. Within this framework, FDR-corrected group effects remained significant for O-info in both deviant and standard responses (q = .001), whereas S-info did not survive correction (Supplement Table). Spatially, maximal early Oinfo effects were centered over central-posterior electrodes, with later configurations shifting toward more frontal-central regions (Fig.). Together, these findings indicate early increases followed by attenuation of task-evoked higher-order interactions, with O-info showing the most robust model-supported effects.

HOI METRICS SHOWED SELECTIVE ASSOCIATIONS WITH CLINICAL OUTCOMES FOLLOWING KETAMINE.

Having established ketamine-related changes in HOI, we next tested whether ΔHOI tracked changes in clinical outcomes. Candidate features were identified using permutationbased thresholding and FDR correction, then evaluated in covariate-adjusted outcome models. Analyses were restricted to the ketamine group, as midazolam does not induce a comparable neurophysiological perturbation and associations in that condition would be expected to reflect non-specific variability rather than mechanistically meaningful coupling.

MADRS:

For antidepressant response, model-supported associations were limited. In the baseline-adjusted primary model, only 24h alpha O-info remained significant after FDR correction (β = 69.31, SE = 25.77, p = .007, q = .050), indicating that greater increases in redundancy at 24h were associated with greater Day 7 symptom improvement (Fig.). Several additional features showed directional effects but did not survive correction, including 24h). Thus, the MADRS findings support a selective and exploratory association, centered primarily on 24h alpha-band redundancy, rather than a broad HOI-antidepressant signal.

CADSS

Associations with dissociative symptoms were stronger and more consistent. In the primary model, four 24h features survived FDR correction: alpha O-info (β = 6.10, q = .007), alpha S-info (β = 5.48, q = .019), beta O-info (β = 34.11, q < .001), and beta S-info (β = 6.00, q = .002) (Fig.). All were positive, indicating that greater increases in HOI were associated with greater dissociative symptoms (Fig.). Earlier 1h and 24h features, including negative gammaband O-info associations, were observed at the candidate-feature stage but did not survive FDR in the outcome model (Supplement Table). Together, these analyses indicate that HOI-clinical associations were outcome-specific: antidepressant effects were limited and centered on a single 24h alpha O-info feature, whereas dissociative effects showed a more robust 24h pattern involving both O-info and S-info. Although HOI features were associated with both dissociative and antidepressant outcomes, the resulting patterns differed in timing and feature composition, suggesting that these associations are not reducible to a single shared clinical process within this dataset.

DISCUSSION

Our primary objective was to identify temporally evolving changes in higher-order brain interactions following ketamine and evaluate their relationship to clinical outcomes in late-life TRD. Building on previous work on HOI, we hypothesized that ketamine would induce temporally dynamic changes in network-level information structure, distinct from those induced by midazolam, and that these changes would relate to clinical improvement. We aimed not only

S ARTICLE IN PRESS

to characterize the neurophysiological statistical organization associated with ketamine exposure but also to evaluate the translational potential of HOI as a clinically viable, scalable biomarker. Notably, we demonstrate these effects using a relatively low-density electrode montage, suggesting that HOI metrics may be deployable in affordable, wireless EEG systems suitable for outpatient settings. Ketamine induced large-scale, frequency-and region-specific shifts in high-order brain dynamics. At 1-hour post-infusion, we observed widespread increases in O-info across frequency bands, peaking in the alpha range, consistent with previous results (Fig.). This carried more overlapping, rather than complementary information, consistent with models proposing reduced hierarchical constraint during auditory prediction, and reinforcing previous literature. By Day 7, the negative effect size shifts suggest that ketamine's longer-term effects may selectively shift the balance between redundancy and synergy without affecting overall interaction strength (Supplementary Fig.). Importantly, the same electrodes contributed to peak O-info for both standard and deviant conditions at the 2h timepoint, replicating previous findings that link these scalp sites to the superior temporal gyrus, a known source of auditory processing 32 . Alpha effects warrant closer examination because they emerged as the earliest and most robust HOI changes following ketamine administration and showed the clearest association with clinical response. Prior studies have reported reductions in alpha power postketamine. Although we did not assess power directly, our findings suggest that higher-order information structure may change independently of conventional spectral interpretations. One might expect reduced alpha power following ketamine to coincide with lower alpha redundancyfewer oscillations and fewer overlapping signals. However, we observed the opposite: alpha redundancy increased. This dissociation between oscillatory power and redundancy is potentially important, as it suggests that HOI captures an aspect of network organization not reducible to band-limited amplitude alone. Direct comparison with spectral power will be an important next step. The evolving effect-size trajectories (Fig.), with alpha peaking at 1 hour, theta at 24 hours, and gamma persisting most clearly by Day 7, suggest a temporally ordered progression, with early alpha effects, intermediate theta effects, and later gamma persistence. Ketamineassociated network influence appeared to abate over time, consistent with Herzog et al.. In our data, alpha-band redundancy may suggest diminished hierarchical control, potentially permitting greater local micro-circuitry flexibility (gamma band) as downstream effects increase. Importantly, however, temporal tracking of maximal n-plets showed limited stability across visits,

S ARTICLE IN PRESS

indicating that these effects were not carried by a single fixed configuration, but rather by changing subsets of n-plets over time (Fig.). Altogether, these results highlight how HOI metrics allow us to detect not just whether brain activity is connected, but how that coordination may be structured across time and frequency: a richer characterization of ketamine-related network dynamics than that provided by traditional spectral or pairwise measures. Associations with clinical outcomes were more selective after permutation testing than those suggested by the initial candidate-feature stage. For MADRS, only 24-hour alpha O-info remained significant after FDR correction in the covariate-adjusted model, including baseline HOI, age, sex, and BMI (Fig.). This finding suggests that increases in alpha band redundancy at 24 hours may track greater Day 7 antidepressant improvement. Other candidate associations, including 24-hour gamma O-info, were directionally informative, but did not survive correction and should therefore be interpreted as exploratory. In sensitivity models excluding baseline HOI, the same 24-hour alpha O-info remained associated with clinical improvement, while one Day 7 theta O-info feature became significant; however, the latter effect was not retained in the baseline-adjusted model and should be interpreted cautiously given marked collinearity between baseline HOI and ΔHOI. These findings align with group-level effect size results, which showed early and prominent increases in alpha-band redundancy following ketamine. Greater channel-level redundancy, especially at 24 hours, could reflect a network state more conducive to clinical recovery. Accordingly, our data support a cautious interpretation: alpha-band redundancy may be a candidate marker of antidepressant-related neural reorganization. More broadly, this pattern is compatible with frameworks such as the relaxed beliefs under psychedelics (REBUS) framework, which proposes that psychedelic-like states increase the importance of sensory evidence and reduce top-down influence on bottom-up signaling. In doing so, they lower the precision weighting of prior beliefs and increase sensitivity to incoming evidence, forcing internal Associations with dissociative symptoms (CADSS change from baseline to 1h) were stronger and more consistent. In the primary covariate-adjusted model, four 24-hour features survived FDR correction: alpha O-info, alpha S-info, beta O-info, and beta S-info, all positively associated with CADSS (Fig.). Thus, greater dissociation was linked to increases in network redundancy and greater higher-order interaction strength, particularly within alpha-and betaband configurations. Notably, the strongest surviving HOI predictors of acute dissociation were observed at 24 hours rather than 1 hour, suggesting that post-acute network changes may retain information about earlier dissociative response. Earlier negative candidate-stage associations, including gamma-band O-info features, did not survive FDR in the outcome model. Still, important limitations remain. Our use of a single-infusion protocol limits generalizability to repeated dosing paradigms, and the late-life, predominantly male veteran sample reduces population representativeness. Although we evaluated dose effects prior to pooling and did not observe dose-dependent relationships, the modest sample size (n=18) limits power to detect more subtle dose heterogeneity. The relatively small midazolam arm and use of adaptive randomizations may have introduced imbalances affecting group comparisons. In addition, the reduced electrode montage, while intentionally aligned with a low-density, clinically scalable setup, limits spatial resolution and may omit combinations that contain finer-grained network structure. MMN parameters are different from prior studies, potentially influencing task-based HOI outcomes. Methodologically, HOI metrics are undirected and therefore cannot resolve the directionality of information flow (e.g., top-down vs. bottom-up processes) but instead characterize how information is distributed across networks. Moreover, while we observed spatial convergence with earlier work, some topographical discrepancies suggest variability in precise electrode contributions. Finally, although our modeling framework accounted for key Future work should explicitly assess the temporal trajectories of HOI metrics under ketamine, as delineating time-dependent dynamics will be essential to understanding their clinical utility and informing dosing strategies. It will also be important to focus on refining spatial modeling of HOI dynamics, integrating HOI with complementary measures like spectral power, and assessing their predictive utility in larger, longitudinal samples. In sum, these findings support HOI as a useful and potentially scalable framework for characterizing ketamine-related changes in higher-order brain organization. Rather than indexing only whether regions co-vary, HOI captures how information is distributed across O-info effects peaked in the alpha band at 1 h, shifted toward theta-band effects at 24 h, and showed partial gamma-band re-emergence by Day 7. Maximal n-plets showed limited temporal stability across visits. Electrode plots for all bands, metrics, and timepoints are depicted in Supplement Figure.

Full Paper Sections

Go Pro to read the extracted full-text sections.

Full Text PDF

Full Paper PDF

Create a free account to open full-text PDFs.

Study Details

References (4)

References cited by this study and indexed in Blossom.

Ketamine versus ECT for Nonpsychotic Treatment-Resistant Major Depression

Anand, A., Mathew, S. J., Sanacora, G. et al. · New England Journal of Medicine (2023)

248 cited
Distinct trajectories of antidepressant response to intravenous ketamine

O'Brien, B., Lijffijt, M., Lee, J. et al. · Journal of Affective Disorders (2021)

17 cited

Your Personal Research Library

Go Pro to save papers, add notes, rate research, and organize custom shelves.