If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Prediction of Behavioral Improvement Through Resting-State Electroencephalography and Clinical Severity in a Randomized Controlled Trial Testing Bumetanide in Autism Spectrum Disorder
1 ELJ-M, JJS, and GC contributed equally to this work.
Erika L. Juarez-Martinez
Footnotes
1 ELJ-M, JJS, and GC contributed equally to this work.
Affiliations
Department of Integrative Neurophysiology, Center for Neurogenomics and Cognitive Research, Neuroscience, VU University Amsterdam, Amsterdam, The NetherlandsNBT Analytics BV, Amsterdam, The NetherlandsChild and Adolescent Psychiatry and Psychosocial Care, Emma Children’s Hospital, Amsterdam UMC, Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
1 ELJ-M, JJS, and GC contributed equally to this work.
Gianina Cristian
Footnotes
1 ELJ-M, JJS, and GC contributed equally to this work.
Affiliations
Child and Adolescent Psychiatry and Psychosocial Care, Emma Children’s Hospital, Amsterdam UMC, Vrije Universiteit Amsterdam, Amsterdam, The NetherlandsDepartment of Health Evidence, Radboud University Medical Center, Nijmegen, The NetherlandsDepartment of Health Evidence, Radboud University Medical Center, Nijmegen, The Netherlands
Department of Integrative Neurophysiology, Center for Neurogenomics and Cognitive Research, Neuroscience, VU University Amsterdam, Amsterdam, The Netherlands
Department of Integrative Neurophysiology, Center for Neurogenomics and Cognitive Research, Neuroscience, VU University Amsterdam, Amsterdam, The NetherlandsNBT Analytics BV, Amsterdam, The Netherlands
Department of Integrative Neurophysiology, Center for Neurogenomics and Cognitive Research, Neuroscience, VU University Amsterdam, Amsterdam, The Netherlands
Child and Adolescent Psychiatry and Psychosocial Care, Emma Children’s Hospital, Amsterdam UMC, Vrije Universiteit Amsterdam, Amsterdam, The NetherlandsDepartment of Psychiatry, UMC Utrecht Brain Center, University Medical Centre Utrecht, Utrecht, The NetherlandsN=You Neurodevelopmental Precision Center, Amsterdam Neuroscience, Amsterdam Reproduction and Development, Amsterdam UMC, Amsterdam, The Netherlands
Department of Integrative Neurophysiology, Center for Neurogenomics and Cognitive Research, Neuroscience, VU University Amsterdam, Amsterdam, The Netherlands
Department of Psychiatry, UMC Utrecht Brain Center, University Medical Centre Utrecht, Utrecht, The NetherlandsDepartment of Biostatistics & Research Support, Julius Center for Health Sciences and Primary Care, University Medical Centre Utrecht, Utrecht, The Netherlands
Department of Integrative Neurophysiology, Center for Neurogenomics and Cognitive Research, Neuroscience, VU University Amsterdam, Amsterdam, The Netherlands
Child and Adolescent Psychiatry and Psychosocial Care, Emma Children’s Hospital, Amsterdam UMC, Vrije Universiteit Amsterdam, Amsterdam, The NetherlandsDepartment of Psychiatry, UMC Utrecht Brain Center, University Medical Centre Utrecht, Utrecht, The NetherlandsN=You Neurodevelopmental Precision Center, Amsterdam Neuroscience, Amsterdam Reproduction and Development, Amsterdam UMC, Amsterdam, The NetherlandsLevvel, Center for Child and Adolescent Psychiatry, Amsterdam, The Netherlands
Mechanism-based treatments such as bumetanide are being repurposed for autism spectrum disorder. We recently reported beneficial effects on repetitive behavioral symptoms that might be related to regulating excitation-inhibition (E/I) balance in the brain. Here, we tested the neurophysiological effects of bumetanide and the relationship to clinical outcome variability and investigated the potential for machine learning–based predictions of meaningful clinical improvement.
Methods
Using modified linear mixed models applied to intention-to-treat population, we analyzed E/I-sensitive electroencephalography (EEG) measures before and after 91 days of treatment in the double-blind, randomized, placebo-controlled Bumetanide in Autism Medication and Biomarker study. Resting-state EEG of 82 subjects out of 92 participants (7–15 years) were available. Alpha frequency band absolute and relative power, central frequency, long-range temporal correlations, and functional E/I ratio treatment effects were related to the Repetitive Behavior Scale-Revised (RBS-R) and the Social Responsiveness Scale 2 as clinical outcomes.
Results
We observed superior bumetanide effects on EEG, reflected in increased absolute and relative alpha power and functional E/I ratio and in decreased central frequency. Associations between EEG and clinical outcome change were restricted to subgroups with medium to high RBS-R improvement. Using machine learning, medium and high RBS-R improvement could be predicted by baseline RBS-R score and EEG measures with 80% and 92% accuracy, respectively.
Conclusions
Bumetanide exerts neurophysiological effects related to clinical changes in more responsive subsets, in whom prediction of improvement was feasible through EEG and clinical measures.
). Novel treatments are being developed on the basis of specific pathophysiological mechanisms, which have been indicated in studies with animal models of genetic disorders (
). Here, we show the results of the first randomized controlled ASD medication trial to incorporate quantitative electroencephalography (EEG) before and after treatment.
A large body of literature exists on quantitative EEG in ASD (
), which may inform ASD treatment development. For instance, excitation/inhibition (E/I) ratio dysregulation has become an important pathophysiological target in ASD since both changes in inhibition and excitation have been implicated (
Electrophysiological biomarkers predict clinical improvement in an open-label trial assessing efficacy of autologous umbilical cord blood for treatment of autism.
Nonlinear measure of EEG complexity in the Eker Rat model of autism disorder: A pilot study.
in: In: ICoMS 2020: Proceedings of the 2020 3rd International Conference on Mathematics and Statistics. Association for Computing Machinery,
New York2020: 49-53
). To date, no pharmacological placebo-controlled clinical trials have used pre- and post-treatment EEG measurements, a strategy that has been tested for nonpharmacological interventions such as neurofeedback and repetitive transcranial magnetic stimulation (
Evaluation of a machine learning model based on pretreatment symptoms and electroencephalographic features to predict outcomes of antidepressant treatment in adults with depression: A prespecified secondary analysis of a randomized clinical trial [published correction appears in JAMA Netw Open 2020;3:e2016001].
), evaluation of EEG biomarkers in ASD together with machine learning analysis may show promise to develop treatment prediction models.
To test the potential of EEG analysis for treatment optimization for ASD, we developed Bumetanide in Autism Medication and Biomarker (BAMBI) (EudraCT 2014-001560-35), the first completed randomized, placebo-controlled ASD medication trial to incorporate EEG measurements to indicate brain effects of bumetanide and to develop stratification and prediction biomarkers of treatment response. Bumetanide is a diuretic drug currently repurposed for ASD treatment. It acts as a selective chloride importer NKCC1 antagonist, which is important to regulate intraneuronal chloride concentration and GABA (gamma-aminobutyric acid) polarity. Animal models have shown that the postnatal GABA shift through downregulation of NKCC1 activity is abolished, leading to excitatory actions of GABA and increased network activity (
). From these findings, it was postulated that bumetanide may have favorable clinical effects since elevated E/I has become an important theory in autistic development (
The potential neurologic effect of bumetanide has been disputed, given its poor brain bioavailability due to limited crossing of the blood-brain-barrier (
Symptom improvement in children with autism spectrum disorder following bumetanide administration is associated with decreased GABA/glutamate ratios [published correction appears in Transl Psychiatry 2020;10:63].
). We recently published the clinical results of the BAMBI study; although there was no superior effect on the primary Social Responsiveness Scale 2 (SRS-2) end point between placebo and bumetanide, we did find a superior effect of bumetanide on the Repetitive Behavior Scale-Revised (RBS-R), a scale of core symptomatology. In addition, a subgroup of most responsive individuals seemed delineated from that of placebo-treated individuals (
Bumetanide for core symptoms of autism spectrum disorder (BAMBI): A single center, double-blinded, participant-randomized, placebo-controlled, phase-2 superiority trial.
J Am Acad Child Adolesc Psychiatry.2021; 60: 865-876
). It is plausible that particularly these subgroups share a mechanistic element targeted by bumetanide. In this subset of participants with stronger improvement, neurophysiological measures may be more prominently related to clinical improvement.
In line with this reasoning, we first hypothesized that resting-state EEG measures would confirm that bumetanide has neurophysiological effects in the brain; second, that clinical improvement would be related to bumetanide’s neurophysiological effects; and third, that EEG measures could be used to develop prediction models for bumetanide in ASD. We used a set of alpha frequency band EEG measures that have proven sensitive to the ratio of excitation and inhibition in computational models, including a method to quantify a functional form of E/I at the network level (fE/I; see Methods and Materials) (
This study is a secondary analysis of the BAMBI trial (Eudra-CT 2014-001560-35), a single-center, double-blind, participant-randomized, placebo-controlled, phase 2 superiority trial testing bumetanide (twice-daily up to 1.0 mg) in otherwise medication-free children with ASD. The trial included children aged 7–15 years with expert-confirmed ASD diagnosis. Exclusion criteria were an IQ < 55 and use of psychoactive medication. Outcomes were assessed at pretreatment (D0), after 91 days of treatment (D91), and after 28-day washout (D119). Clinical outcome measures were the SRS-2, the RBS-R, and resting-state EEG. The study was approved by the medical ethical committee of the UMC Utrecht, Utrecht, the Netherlands, and all participants or their legal guardians signed informed consent. The CONSORT (Consolidated Standards of Reporting Trials) flow diagram with the included participants and measurements performed is given in Figure 1. Detailed description of the protocol and clinical effects have been published previously (
Bumetanide for core symptoms of autism spectrum disorder (BAMBI): A single center, double-blinded, participant-randomized, placebo-controlled, phase-2 superiority trial.
J Am Acad Child Adolesc Psychiatry.2021; 60: 865-876
Figure 1CONSORT (Consolidated Standards of Reporting Trials) flow diagram for the patients included for electroencephalography (EEG) analysis from the Bumetanide in Autism Medication and Biomarker trial. A total of 125 subjects were screened for eligibility between June 2016 to December 2018 at the UMC Utrecht, Utrecht, the Netherlands. Ninety-two participants were randomized to receive bumetanide (n = 47) or placebo (n = 45). EEG recordings of 82 participants were available for analysis (bumetanide n = 42; placebo n = 40). ASD, autism spectrum disorder; DFA, detrended fluctuation analysis; fE/I, functional excitation-inhibition ratio; KCL, potassium chloride; LRTC, long-range temporal correlations; RBS-R, Repetitive Behavior Scale-Revised; SRS-2, Social Responsiveness Scale 2.
EEGs were recorded during 5 minutes of eyes-closed rest with a 64-channel BioSemi system (2048 Hz). EEG analyses were done using the Neurophysiological Biomarker Toolbox and custom-made scripts (
). All recordings were manually cleaned for artifacts and re-referenced to the average reference. After preprocessing, on average 210 seconds per recording (80–327 s) were available for analysis.
EEG Analysis
Computational neuronal network models generating oscillations have shown that changes in excitation and inhibition affect the amplitude, frequency, temporal correlations, and fE/I in the alpha band (8–13 Hz) (
). Therefore, we applied these measures to the EEG data in the alpha band.
Spectral power was computed using the Welch method with an 8192-point Blackman window and a frequency resolution of 0.125 Hz. Relative alpha power is expressed in percent and was calculated by dividing the absolute power in the alpha band by the integrated power in the range of 1 to 45 Hz. Central frequency provides a measure of the frequency at which most of the power in the alpha band is concentrated. Central frequency has the advantage over peak frequency that it does not require the presence of a strong peak or—as it sometimes appears—that the peak is identified close to 8 Hz because of increasing power at lower frequencies. The detrended fluctuation analysis was used to quantify the long-range temporal correlations in the amplitude fluctuation of alpha oscillations (
). The detrended fluctuation analysis exponents in the interval of 0.5 to 1.0 indicate long-range temporal correlation and rich temporal structure fluctuations, whereas an exponent of 0.5 characterizes an uncorrelated signal. Finally, we used the recently introduced measure of functional E/I ratio (fE/I), which has been shown to track pharmacological intervention and reveal large heterogeneity in an ASD sample (
). Inhibition-dominated networks are characterized by an fE/I < 1, excitation-dominated networks fE/I > 1, and E/I-balanced networks will have fE/I = 1. A detailed description of the EEG measures is provided in the Supplement.
Statistical Analysis
The secondary analyses in this paper were conducted according to the predefined statistical analysis in the study protocol, which included the development of EEG biomarkers for stratification and functional assessment of treatment effects [see (
Bumetanide for core symptoms of autism spectrum disorder (BAMBI): A single center, double-blinded, participant-randomized, placebo-controlled, phase-2 superiority trial.
J Am Acad Child Adolesc Psychiatry.2021; 60: 865-876
To assess whole-brain treatment effect of bumetanide versus placebo on EEG measures, we analyzed the average biomarker value (mean) of the 64 channels using a linear mixed model (SAS, version 9.4; SAS Institute). Treatment and treatment-by-time interaction were included to assess the difference between bumetanide and placebo and washout effects. Sex, age, and pretreatment EEG measures were included to correct for potential confounding factors (
). We derived estimated means for each treatment group and a mean difference between treatment groups at 91 days with 95% confidence intervals and p values. This analysis and the machine learning parts were run as sensitivity analyses, excluding participants with unreliable clinical reports as motivated in the initial report of BAMBI (
Bumetanide for core symptoms of autism spectrum disorder (BAMBI): A single center, double-blinded, participant-randomized, placebo-controlled, phase-2 superiority trial.
J Am Acad Child Adolesc Psychiatry.2021; 60: 865-876
). For the EEG effects irrespective of clinical report, we also ran an intention-to-treat analysis. The topographical distribution of EEG treatment effects was analyzed per treatment group by comparing each electrode at different time points using a paired t test. False discovery rate was used to correct for multiple testing at the electrode level (see Supplemental Methods).
Clinical Outcome Analyses
The SRS-2 (primary end point; total raw score; range 0−195; higher score indicates more affected) and RBS-R (secondary end point; total score; range 0−129; higher score indicates more affected) were used to relate EEG parameters to clinical outcome in core symptomatology. The original study findings showed that a subset of participants with the greatest SRS-2 and RBS-R improvement were in the bumetanide group, which we took as a basis to define thresholds for treatment response stratification and prediction. We set a high threshold to indicate the level of clinical improvement outcome above which there were no or few (n = 3 for SRS-2) placebo-treated subjects present (16 points for RBS-R and 20 points for SRS-2) and may be more successful to identify true responders in terms of effects being related to bumetanide targeting mechanisms. In addition, we defined a medium response threshold, using mean clinical score improvement as a reference (7 and 10 points for RBS-R and SRS-2, respectively) to test stability of stratification and prediction under a medium treatment response threshold. We investigated the association between EEG and symptom change in the whole sample and for the different treatment thresholds to test the stability of correlations under different stratification scenarios.
Machine Learning Models for Clinical Outcome Predictions
The data preprocessing, feature selection, classifier training, and validation steps were performed in RStudio (version 1.3.959) (
Figure 2Machine learning workflow. The dataset comprised 5 electroencephalography biomarkers over 63 channels and 5 clinical features used to predict clinical improvement of n = 34 bumetanide-treated subjects above or below the thresholds Repetitive Behavior Scale-Revised: 7- and 16-point. To develop and validate the random forests classifier using a shuffle-split method, the data were divided into training and validation sets of 80% and 20%, respectively, with a ±5% margin to ensure a minimum number of above threshold improvement subjects in each validation set. To generalize results and aid replicability, these steps were repeated 100 times, and each time, a new random split was generated. For each of the 100 training splits, the most optimal features were selected using a recursive method to discriminate subjects with outcomes above and below the selected treatment thresholds. The selected feature sets served as input for the random forest classifiers trained on each of the 100 training partitions. This process generated 100 models with different selected feature constellations that were optimized using leave-one-out cross-validation. The resulting optimal models were validated on the out-of-sample partitions. The number of correctly and incorrectly classified datapoints for each of the 100 iterations rendered 100 accuracy results, which were then averaged to obtain overall results for each improvement threshold.
For D91 treatment response predictions from pretreatment EEG using machine learning analysis, we compiled the pool of potential features from 5 resting-state measures on 63 channels. We also incorporated pretreatment SRS-2 and RBS-R scores, total IQ, age, and a binary variable describing whether they had received treatment prior to the study, resulting in a total of 5 clinical features, thus 320 features in total. After data preprocessing, using a shuffle-split technique, which randomly and repeatedly splits the data (
), we created 100 random data partitions for training and 100 random data partitions for out-of-sample validation. In each partition, 80% (±5%) of the data was reserved for training and testing and 20% (±5%) for validation.
Feature selection was performed to identify from the 320 starting features the subsets that could classify subjects for being below or above medium and high response treatment thresholds. This resulted in 100 feature subsets, 1 per training partition. For clinical interpretation of which features were most important in the prediction models, we reported the 10 features with the highest frequency of appearance in 70%-or-higher raw accuracy runs as most informative.
Training, Hyperparameter Tuning, and Testing
A fast implementation of a random forest classifier (
) was trained and optimized 100 times on each of the training sets containing 80% (±5%) of the data, using the feature subsets obtained in the feature selection step using leave-one-out cross-validation.
Validation
The leave-one-out cross-validation results were validated on each of the 100 out-of-sample sets containing 20% (±5%) of the data. The performance measures expressed in terms of balanced and raw accuracy, sensitivity, and specificity were averaged over the out-of-sample runs. We report the raw accuracy showing all the correctly classified observations irrespective of class. Balanced accuracy was reported to even out imbalances between sensitivity and specificity resulting from the higher prevalence of low-clinical-improvement cases (Figure S2).
Results
Ninety-two participants were randomized, of whom 88 completed the trial (
Bumetanide for core symptoms of autism spectrum disorder (BAMBI): A single center, double-blinded, participant-randomized, placebo-controlled, phase-2 superiority trial.
J Am Acad Child Adolesc Psychiatry.2021; 60: 865-876
). EEG recordings of 82 of these trial completers were available for analysis: 42 from the bumetanide group (mean age 10.5 ± 2.5 years; 14 females) and 40 from the placebo group (mean age 10.5 ± 2.5 years; 12 females). EEG recordings from 10 participants (n = 12) were excluded after preprocessing because of poor quality of EEG data or no eyes-closed rest EEG could have been recorded (D0 n = 6, D91 n = 3, D119 n = 3; 1 participant had poor quality EEG at each time point). This left 76 EEGs for analysis at D0, 79 at D91, and 79 at D119 (for details, see Figure 1). Tables S1 and S2 describe the participants included, their demographics, clinical scores, and EEG biomarker values. To test EEG effects of bumetanide, we investigated 5 EEG measures in the alpha band (8–13 Hz) that have proven sensitive to pharmacological E/I intervention (see Methods and Materials).
Effect of Bumetanide Treatment on EEG Measures
After 91 days, superior EEG effects of bumetanide to placebo were found for absolute power, relative power, fE/I, and central frequency (Table 1 and Figure 3A–D). These effects were of medium size (Cohen’s d = 0.5), except for central frequency showing a small effect size (Cohen’s d = 0.3). No effects were found between D91 and the D119 washout period. Subanalysis of treatment interaction with age, IQ, and attention-deficit/hyperactivity disorder comorbidity did not show significant effects. These findings were not different in the intention-to-treat analysis (Table S3). Subsequently, we analyzed the topographical distribution of EEG changes between D0 and D91 in each treatment group. Significant effects were only present in the bumetanide group, where absolute power, relative power, and fE/I increased and central frequency decreased in a cluster of channels (Figure 3E–H).
Table 1Superior Effect of Bumetanide to Placebo in Power, Central Frequency, and fE/I
Data are shown for participants who completed D91. Linear mixed model analysis of whole-brain average electroencephalography (EEG) biomarkers reveal superior effect of bumetanide to placebo on absolute power, relative power, central frequency, and fE/I after 91 days of treatment. Treatment effect, 95% confidence intervals, and p values are reported.
D0, day 0 baseline recording; D91, day 91 of treatment; D119, day 119 (after 28-day washout period); DFA, detrended fluctuation analysis; fE/I, functional excitation-inhibition ratio.
Figure 3Bumetanide treatment—but not placebo—affects alpha oscillations in children with autism spectrum disorder. (A–D) Whole-brain average electroencephalography (EEG) measures at treatment time points. No baseline differences (D0) were found for EEG measures between treatment groups (A–D) (t test). Linear mixed model analysis of treatment effects (D0–D91) showed a superior effect of bumetanide to placebo in absolute power (A), relative power (B), central frequency (C), and functional excitation-inhibition ratios (fE/I) (D) (see also Table 1). Analysis of treatment-by-time interaction revealed no washout effect at D119 for any of the EEG measures. Bumetanide group (solid red lines). Placebo group (black dashed lines). Mean ± SEM are plotted. ∗p < .05; ∗∗p < .01. (E–H) Grand-average topographies of the treatment effects on EEG measures (D91−D0). Significant channels are shown in white circles (p < .05, false discovery rate corrected) indicating widespread effects in the bumetanide group with an increase in absolute power (E), relative power (F), and fE/I (H) and a cluster of channels with decreased central frequency (G). No significant channels were found in the placebo group. D0, day 0 baseline recording; D91, day 91 of treatment; D119, day 119 (after 28-day washout period).
Correlation Between EEG Measures and Clinical Outcome
Correlations of EEG biomarker D91-D0 changes with change in SRS-2 or RBS-R after bumetanide treatment did not reveal significant associations at the group level. Stratification by the medium treatment response threshold (RBS-R ≥ 7 points improvement [n = 8]) did reveal high, significant correlations between an increment in absolute and relative power and RBS-R improvement after bumetanide (Figure 4A, B), albeit just a trend for fE/I (Figure 4C). The subgroup with improvement < 7 points (n = 24) showed no significant correlations between EEG changes and clinical outcome (Figure 4D, E) or even an opposite trend for fE/I by contrast to the ≥7 points group (Figure 4C, F). No correlations were found for the RBS-16 threshold or both SRS-2 tested thresholds.
Figure 4Correlations between electroencephalography and clinical outcome changes in participants with below (D–F) and equal to or exceeding (A–C) 7 points improvement on the Repetitive Behavior Scale-Revised (RBS-R). (A, B) An increment in absolute and relative power after bumetanide was strongly correlated (rho ≤ −0.8) to repetitive behavior improvement, in participants with improvement of at least 7 points. (C) A medium-high (rho = −0.6), albeit nonsignificant, correlation was observed between an increment in functional excitation-inhibition ratio (fE/I) and repetitive behavior improvement. (D, E) No significant correlations were found between electroencephalography change and repetitive behavior improvement in participants with improvement below 7 points. (F) An increase in fE/I is moderately (rho = 0.4) and borderline significantly (p = .06) associated with worsening in repetitive behavior in the opposite direction when compared with the ≥7 points improvement group (C). Correlations were calculated using Spearman correlation coefficient. Channels with significant (p < .05, false discovery rate corrected) correlations are shown in white circles. Abs., absolute; Rel., relative.
Prediction of Clinical Improvement Through Machine Learning Analysis
To test the potential of clinical outcome predictions, we entered channel-level pretreatment EEG measures together with 5 clinical features (pretreatment SRS-2 and RBS-R scores, total IQ, age, and prior treatment) into a random forest machine learning classifier to predict clinical improvement of bumetanide-treated subjects above or below the thresholds RBS-R: 7- and 16-point, and SRS-2: 10- and 20-point (for details, see Methods and Materials) (Figures 5 and 6; Figure S1). The distributions of the individual change in SRS-2 and RBS-R clinical outcomes in relationship to each of the chosen treatment thresholds are illustrated in Figures 5A and 6A for those individuals included in the machine learning analysis (n = 34).
Figure 5Machine learning delineation between clinical improvement levels did not exceed chance when predicting improvement in social responsiveness. (A) Distribution of the individual change in Social Responsiveness Scale 2 (SRS-2) (D91−D0) in relation to treatment thresholds used for machine learning analysis (SRS 10-point orange line and SRS 20-point yellow line). Data presented are only for individuals with pretreatment electroencephalography measures and clinical outcome scores available (bumetanide n = 34, placebo n = 35). (B) Receiver operating characteristic (ROC) curves to illustrate the diagnostic ability of the binary classifiers for the different improvement thresholds across the SRS-2 clinical outcome scale. This ability is expressed by the relationship between true (sensitivity) and false-positive rates, as depicted by the ROC curves. The curve for the perfect model that correctly classifies 100% below and above threshold improvers contains the point (0, 1), while a random estimate curve is closer to the nondiscrimination line (gray). An area under the curve (AUC) closer to 1 indicates a classifier having both high sensitivity and specificity. Classifiers for 10- and 20-point SRS improvement demonstrate close to or worse than chance performance. (C) Classification performance of SRS 10-point and SRS 20-point treatment thresholds. As per out-of-sample validation, the classifiers predicted 10- and 20-point improvement on SRS-2 with an average sensitivity of 44% and 25% and 57% and 77% average specificity for each threshold, respectively, which corresponded to 56% and 51% average balanced accuracy (mean between specificity and sensitivity), respectively. (D, E) The 10 most frequently selected electrodes and electroencephalography measures in iterations exceeding 70% accuracy on out-of-sample validation partitions are shown for the outcome thresholds. Bal., balanced; C. I., confidence interval; D0, day 0 baseline recording; D91, day 91 of treatment; DFA, detrended fluctuation analysis; fE/I, functional excitation-inhibition ratio.
Figure 6Machine learning shows best delineation between clinical improvement levels when predicting improvement in repetitive behavior. (A) Distribution of the individual change in Repetitive Behavior Scale-Revised (RBS-R) (D91−D0) in relation to treatment thresholds used for machine learning analysis (RBS 7-point blue line and RBS 16-point purple line). The data presented are only for individuals with pretreatment electroencephalography measures and clinical outcome scores available (bumetanide n = 34, placebo n = 35). (B) Receiver operating characteristic (ROC) curves to illustrate the diagnostic ability of the binary classifiers for the different improvement thresholds across the RBS-R outcome scales. This ability is expressed by the relationship between true (sensitivity) and false-positive rates, as depicted by the ROC curves. The curve for the perfect model that correctly classifies 100% below and above threshold improvers contains the point (0, 1), while a random estimate curve is closer to the nondiscrimination line (gray). An area under the curve (AUC) closer to 1 indicates a classifier having both high sensitivity and specificity. The best performance is indicated for 16-point RBS-R improvement, followed by a 7-point reduction. (C) Classification performance of RBS-R 7-point and RBS-R 16-point treatment thresholds. As per out-of-sample validation, the classifiers predicted 7-point and 16-point improvement on RBS-R with an average sensitivity of 63% and 86%, and 96% and 98.8% average specificity for each threshold, respectively, which corresponded to 79.5% and 92.4% average balanced accuracy (mean between specificity and sensitivity), respectively. (D, E) The 10 most frequently selected electrodes and electroencephalography measures in iterations exceeding 70% accuracy on validation partitions are shown for the two outcome thresholds. The baseline RBS-R score was most frequently selected as an important feature to inform RBS-R improvement, followed by absolute power in centroparietal channels (CP5 and CP3) and functional excitation-inhibition ratio (fE/I) in central channels (C6) for a medium RBS-R improvement (RBS-R 7-point) (D), and absolute power in centroparietal channels (P3, CP5, CP3 and P5) for the high RBS-R improvement threshold (RBS-R 16-point) (E). Bal., balanced; C. I., confidence interval; D0, day 0 baseline recording; D91, day 91 of treatment; DFA, detrended fluctuation analysis.
To generalize results and aid replicability, the classifiers were trained in 100 runs using leave-one-out cross-validation and further validated on 100 out-of-sample partitions. For SRS-2, 10- and 20-point improvement predictions, the out-of-sample validation accuracy, sensitivity, and specificity did not exceed chance levels (Figure 5B, C; Figure S3A, B). In contrast, the classifiers predicted 7- and 16-point improvement on RBS-R with 80% and 92% average balanced accuracy, the equivalent of 63% and 86% average sensitivity and 96% and 99% average specificity for each threshold, respectively (Figure 6C; Figure S3C, D). The 16-point RBS-R improvement had a better trade-off between sensitivity and specificity than the 7-point improvement (Figure 6B).
The selection of features with the highest frequency of appearance (out of 100) in runs achieving a high validation accuracy (>70%) revealed that there was no consistent pattern for SRS-2 in the most predictive features (Figure 5D, E), apparently in line with poor prediction success for this scale. RBS-R baseline score was ranked as the most important predictor of RBS-R improvement, followed by EEG features of absolute power and fE/I in the parietal and central-parietal regions (Figure 6D, E).
Discussion
The presented analysis of resting-state placebo-controlled EEG measures has implications for the application of bumetanide and potentially for other similar ASD treatments. First, bumetanide and not placebo altered spectral and temporal characteristics of neuronal oscillations, indicating a neurophysiological treatment effect. Second, changes in brain activity after bumetanide were only related to improvement in repetitive behavior in more responsive subsets. Finally, predictions of improvement in repetitive behavior were feasible by implementing pretreatment EEG and clinical severity in machine learning analysis.
The observed EEG effects suggest that bumetanide enters the brain sufficiently to alter both power and network-level E/I of neuronal oscillations after 91 days of treatment, effects that did not significantly decrease after the washout period. Other studies have indicated functional brain effects of bumetanide in ASD using magnetic resonance spectroscopy and eye tracking, albeit in open-label trial designs (
Symptom improvement in children with autism spectrum disorder following bumetanide administration is associated with decreased GABA/glutamate ratios [published correction appears in Transl Psychiatry 2020;10:63].
Symptom improvement in children with autism spectrum disorder following bumetanide administration is associated with decreased GABA/glutamate ratios [published correction appears in Transl Psychiatry 2020;10:63].
). Another question is whether the observed EEG effects relate to a shift from depolarizing to hyperpolarizing GABAergic transmission. In healthy subjects, GABAA receptor positive allosteric modulators such as benzodiazepines have been shown to reduce alpha power (
). We recently described a decrease in fE/I following the administration of the positive allosteric modulator zolpidem in healthy adults, suggesting enforcement of inhibition, in line with the supposed drug effect (
). Bumetanide seems to show opposite alpha-oscillation effects in terms of power and fE/I in ASD when compared with the acute effect of positive allosteric modulators in healthy subjects. A similar opposite E/I effect in ASD has also been shown previously for the E/I-regulating drug, riluzole, measured by magnetic resonance spectroscopy (
). Importantly, our finding of an increase in fE/I through bumetanide in ASD is consistent with the open-label trial that found an increase in glutamate-to-GABA magnetic resonance spectroscopy concentrations in sensory cortices after 3 months of bumetanide treatment (
Symptom improvement in children with autism spectrum disorder following bumetanide administration is associated with decreased GABA/glutamate ratios [published correction appears in Transl Psychiatry 2020;10:63].
)—a shift toward a larger ratio between excitatory and inhibitory neurotransmission increases the power and fE/I also in the computational model used to develop the measure of E/I (
). Nonetheless, the direction of effects is counterintuitive when bumetanide is expected to reinstate hyperpolarizing GABA activity. Additional experiments are needed to gain understanding of the bumetanide effects on EEG, but we speculate that our network-level measure of E/I is sensitive to compensatory mechanisms, which might also be affected after 3 months of bumetanide treatment and may explain the seemingly paradoxical increase in network-level E/I. A homeostatic view has been postulated to gain understanding on how activity propagates in cortical circuits after E/I disruptions and how they may predict contradictory findings in ASD (
). In this scenario, we hypothesize from our findings that depolarizing GABA activity in neurodevelopmental disorders is balanced by a reduction in excitation. A bumetanide-induced shift in GABA polarity from depolarizing to hyperpolarizing effects may reduce the need for the compensatory excitatory downregulation and hence manifest as fE/I elevation in EEG signals. To test these ideas further, we plan follow-up EEG studies comparing immediate and sustained effects of treatment with bumetanide. Several other existing off-label drugs with a variety of GABAergic or glutamatergic actions have been tested in ASD using well-powered placebo-controlled trial designs (
). Most of these trials did not meet their primary end points, but most studies did find effects on secondary end points or subscales. Our findings encourage us to test these observations further and incorporate functional measures such as the EEG measures presented here with the collective aim to gain understanding on how they influence complex homeostatic E/I balance regulations and develop more stratified application strategies that acknowledge the physiological heterogeneity of ASD.
We did not find significant linear associations between EEG and outcome measures at the group level. We conducted subsequent correlation analyses testing different responsiveness thresholds under the hypothesis that participants with stronger improvement, and putatively less mixed with placebo-effects, might represent more true responsive subjects, i.e., where treatment effects are more likely related to the mechanism of action. The correlations between EEG and outcome measures observed in patients above the medium improvement threshold in repetitive behavior (RBS-R) may support this hypothesis and indicate utility of EEG for stratification of individuals likely to benefit from treatment. However, these subsample correlations need to be regarded as preliminary, requiring further validation as they lack power for generalization.
Through the machine learning analysis, we shifted from subgroup delineation to individualized predictions based on unique EEG and clinical values for each patient. The machine learning results indicated that baseline EEG measures can be used to predict improvement on the RBS-R scale in the bumetanide group with high accuracy, under both cross-validation and out-of-sample scenarios. The highest sensitivity (86%) was obtained when predicting a 16-point improvement on RBS-R, which suggests that pretreatment EEG measures and pretreatment RBS-R scores are more instrumental in identifying a high rather than a medium clinical improvement. SRS-2 predictions failed, which may be consistent with the fact that no superior treatment effect was found on this measure. We discussed in the first report of the BAMBI trial that repetitive behavior changes may be more readily observed after 3 months treatment than social behavioral change, which may be a more complex phenotype requiring longer treatment duration or additional behavioral therapy (
). Indeed, we found that RBS-R baseline score was selected most frequently in the predictions of both medium and high improvement thresholds. The results indicate that baseline RBS-R severity together with the most frequently selected absolute alpha power and fE/I in the parietal and central-parietal regions are most informative in predicting clinical improvement. The parietal cortex is relevant to ASD as it is a key area for integration of multiple stimuli (
). However, further prospective selected trials are needed to corroborate these feature combinations as potential ASD biotypes amenable to improvement through bumetanide treatment. Together, the correlations between the EEG and clinical changes and the successful prediction of clinical improvement indicate clinical importance of the neurophysiological effect induced by bumetanide.
It is possible that bumetanide EEG effects are associated with improvement in symptoms intuitively related to network excitability, such as lack of energy, affect problems, or sleeping difficulty, which are not covered by the SRS-2 or RBS-R scales. Evaluation of individual, patient-relevant symptoms (rather than summed severity scores) may be valuable in future research. We restricted the current analysis to a selected set of EEG measures in the alpha frequency band that we previously found sensitive to the ratio of excitation and inhibition in computational models of neuronal networks (
) and frequency bands can be considered for future studies. For the present study, however, this would have resulted in an excessive number of features relative to the sample size available for machine learning. The modest sample size coupled with a high dimensionality of the dataset may also have resulted in overfitting, albeit this was mitigated by out-of-sample validation. The validation in this study was conducted in such a manner that the optimistic bias resulting from a limited sample size was minimized. In addition, we reduced high dimensionality by selecting only the most informative features for prediction during feature selection. Because of the distribution of treatment responses in the sample, the high clinical improvement class is less likely to be identified correctly. This was reflected by sensitivity being generally lower than the specificity.
Conclusions
Quantitative analysis of alpha band resting-state EEG shows that bumetanide has neurophysiological effects in children with ASD. The out-of-sample treatment response predictions support the applicability of machine learning classification to future patient cohorts and encourage further development of EEG-assisted treatment decision-support systems.
Acknowledgments and Disclosures
This work was supported by the Netherlands Organization for Scientific Research (NWO) Physical Sciences (Grant No. 612.001.123 [to KL-H]), Netherlands Organization for Scientific Research (NWO) Social Sciences (Grant No. 406-15-256 [to A-EA and KL-H]), ZonMW Rational Pharmacotherapy Program, GGG project (Grant No. 836041015 [to HB and JJS]), and EU H2020 “Human Brain Project” (Grant No. 604102 [to HDM]).
Neither the funders of the study nor Neurochlore, who provided the study medication, had a role in study design, data collection, data analysis, data interpretation, or writing of the report.
Conceptualization: HB and KL-H; methodology: ELJ-M, JJS, GC, GJvdW, MJCE, KL-H, HB; investigation: ELJ-M, JJS, GC; visualization: ELJ-M, JJS, GC, SS, A-EA, SJH; acquisition, analysis, or interpretation of data: all authors; funding acquisition: KL-H, HB; administrative, technical, or material support: SS, A-EA, RH, SJH, BO; supervision: HDM, KL-H, HB; writing - original draft: ELJ-M, JJS, GC, KL-H, HB; and writing - reviewing and editing: all authors.
We thank Simon-Shlomo Poil for contributions to the Neurophysiological Biomarker Toolbox. We thank Neurochlore for providing the study medication.
Due to privacy regulations of human subjects, we cannot provide the electroencephalography files of the subjects included in our study. However, we provide the individual demographics and the mean of electroencephalography measures in Tables S1 and S2. Analysis scripts to reproduce the figures and statistics will be made available on figshare (https://figshare.com/). The code for the fE/I algorithm is publicly available at https://github.com/rhardstone/fEI.
KL-H and Simon-Shlomo Poil are shareholders of NBT Analytics BV, which provides EEG analysis services for clinical trials. HB, KL-H, and Simon-Shlomo Poil are shareholders of Aspect Neuroprofiles BV, which develops physiology-informed prognostic measures for neurodevelopmental disorders. RH and KL-H have filed the patent claim (PCT/NL2019/050167) “Method of determining brain activity”; with priority date March 16, 2018. All other authors report no biomedical financial interests or potential conflicts of interest.
Electrophysiological biomarkers predict clinical improvement in an open-label trial assessing efficacy of autologous umbilical cord blood for treatment of autism.
Nonlinear measure of EEG complexity in the Eker Rat model of autism disorder: A pilot study.
in: In: ICoMS 2020: Proceedings of the 2020 3rd International Conference on Mathematics and Statistics. Association for Computing Machinery,
New York2020: 49-53
Evaluation of a machine learning model based on pretreatment symptoms and electroencephalographic features to predict outcomes of antidepressant treatment in adults with depression: A prespecified secondary analysis of a randomized clinical trial [published correction appears in JAMA Netw Open 2020;3:e2016001].
Symptom improvement in children with autism spectrum disorder following bumetanide administration is associated with decreased GABA/glutamate ratios [published correction appears in Transl Psychiatry 2020;10:63].
Bumetanide for core symptoms of autism spectrum disorder (BAMBI): A single center, double-blinded, participant-randomized, placebo-controlled, phase-2 superiority trial.
J Am Acad Child Adolesc Psychiatry.2021; 60: 865-876