A Quantitative Electrophysiological Biomarker of Duplication 15q11.2-q13.1 Syndrome

Abstract

Background

Duplications of 15q11.2-q13.1 (Dup15q syndrome) are highly penetrant for autism spectrum disorder (ASD). A distinct electrophysiological (EEG) pattern characterized by excessive activity in the beta band has been noted in clinical reports. We asked whether EEG power in the beta band, as well as in other frequency bands, distinguished children with Dup15q syndrome from those with non-syndromic ASD and then examined the clinical correlates of this electrophysiological biomarker in Dup15q syndrome.

Methods

In the first study, we recorded spontaneous EEG from children with Dup15q syndrome (n = 11), age-and-IQ-matched children with ASD (n = 10) and age-matched typically developing (TD) children (n = 9) and computed relative power in 6 frequency bands for 9 regions of interest (ROIs). Group comparisons were made using a repeated measures analysis of variance. In the second study, we recorded spontaneous EEG from a larger cohort of individuals with Dup15q syndrome (n = 27) across two sites and examined age, epilepsy, and duplication type as predictors of beta power using simple linear regressions.

Results

In the first study, spontaneous beta1 (12–20 Hz) and beta2 (20–30 Hz) power were significantly higher in Dup15q syndrome compared with both comparison groups, while delta (1–4 Hz) was significantly lower than both comparison groups. Effect sizes in all three frequency bands were large (|d| > 1). In the second study, we found that beta2 power was significantly related to epilepsy diagnosis in Dup15q syndrome.

Conclusions

Here, we have identified an electrophysiological biomarker of Dup15q syndrome that may facilitate clinical stratification, treatment monitoring, and measurement of target engagement for future clinical trials. Future work will investigate the genetic and neural underpinnings of this electrophysiological signature as well as the functional consequences of excessive beta oscillations in Dup15q syndrome.

About the Authors:

Joel Frohlich

* E-mail: [email protected]

Affiliation: Center for Autism Research and Treatment University of California, Los Angeles Semel Institute for Neuroscience 760 Westwood Plaza, Suite A7-452 Los Angeles, California United States of America

Damla Senturk

Affiliation: Department of Biostatistics UCLA School of Public Health Room 21-254C CHS Los Angeles, CA United States of America

Vidya Saravanapandian

Affiliation: Center for Autism Research and Treatment University of California, Los Angeles Semel Institute for Neuroscience 760 Westwood Plaza, Suite A7-452 Los Angeles, California United States of America

Peyman Golshani

Affiliation: Department of Neurology and Psychiatry David Geffen School of Medicine 710 Westwood Plaza Los Angeles, CA United States of America

Lawrence T. Reiter

Affiliation: Departments of Neurology and Pediatrics The University of Tennessee Health Science Center 855 Monroe Ave., Link 415 Memphis, TN United States of America

Raman Sankar

Affiliation: Departments of Neurology and Pediatrics David Geffen School of Medicine at UCLA Los Angeles, CA United States of America

Ronald L. Thibert

Affiliation: Department of Neurology, Massachusetts General Hospital Harvard Medical School Boston, MA United States of America

Charlotte DiStefano

Affiliation: Center for Autism Research and Treatment University of California, Los Angeles Semel Institute for Neuroscience 760 Westwood Plaza, Suite A7-452 Los Angeles, California United States of America

Scott Huberty

Affiliation: Center for Autism Research and Treatment University of California, Los Angeles Semel Institute for Neuroscience 760 Westwood Plaza, Suite A7-452 Los Angeles, California United States of America

Edwin H. Cook

Affiliation: Department of Psychiatry University of Illinois at Chicago (MC 747) Chicago, IL United States of America

Shafali S. Jeste

Affiliation: Center for Autism Research and Treatment University of California, Los Angeles Semel Institute for Neuroscience 760 Westwood Plaza, Suite A7-452 Los Angeles, California United States of America

Introduction

Advances in genetic testing have accelerated the diagnosis of causative genetic syndromes in the context of neurodevelopmental disorders, in particular autism spectrum disorder (ASD), with the diagnostic yield of combined chromosomal microarray and whole exome sequencing of 3–20% [1–3]. In parallel, there has been an emerging interest in defining neurophysiological markers that may stratify children with ASD into biologically meaningful subgroups [4]. The convergence of discoveries in genetics and neurophysiology in ASD holds tremendous potential for the identification of biomarkers, grounded in genetic mechanisms, that can improve diagnosis, selection of treatment targets, and treatment monitoring in future clinical trials for genetically defined syndromes within the spectrum of neurodevelopmental disorders

This convergence of genetics with neurophysiology finds promise in the investigation of duplications of 15q11.2-q13.1 (or Dup15q syndrome), one of the most common copy number variants (CNVs) associated with ASD, accounting for 1–3% of cases [5]. Clinical features of Dup15q syndrome include hypotonia, global developmental delay, intellectual disability, social communication impairments, and often severe epilepsy, including infantile spasms [6–9]. Two duplication types exist: isodicentric duplications, which are characterized by two extra copies of the 15q11.2-q13.1 region of maternal origin on a supernumerary chromosome, and interstitial duplications, which can be one or more extra copies of this region on the q arm of chromosome 15 [8].

The 15q11.2-q13.1 region contains many genes critical to neural function, including UBE3A, a ubiquitin-protein ligase important in synaptic function [10–12], and three gamma aminobutyric acid type A receptor (GABA_AR) subunits genes, GABRA5, GABRB3, and GABRG3. Recently, two studies have qualitatively noted the presence of persistent beta frequency (12–30 Hz) activity in clinical electroencephalogram (EEG) recordings from children with Dup15q syndrome [6, 7]. The clinical EEGs are quite distinctive, with spontaneous beta oscillations (SBOs) evident with visual inspection alone. These SBOs resemble the EEG effects of allosteric modulation of GABA_ARs by benzodiazepine drugs, even though the children for whom this observation was reported were not taking these medications [6, 7]. Since the SBOs in these studies occurred in both maternal and paternal duplication cases, it is possible that these SBOs could result from aberrant GABAergic transmission and not primarily from the maternally expressed UBE3A gene.

Herein, we sought to quantify EEG beta power in Dup15q syndrome to determine if it distinguished Dup15q syndrome from non-syndromic ASD and typical development. To accomplish this goal, we obtained spontaneous EEG recordings from (1) children with Dup15q syndrome, (2) age-and-IQ matched children with non-syndromic ASD and (3) age matched typically developing (TD) children. We hypothesized that spontaneous beta power would differentiate children with Dup15q syndrome from these comparison groups. In a follow-up study we examined the variation in beta power within a larger Dup15q cohort by analyzing age, duplication type, and epilepsy status as predictors of SBO strength. Given the likelihood that SBOs are related to copy number variation and seizures in Dup15q syndrome, we hypothesized that both duplication type and epilepsy would relate to spontaneous beta power.

Subjects and Methods

Study 1: Comparison of Dup15q syndrome with ASD and TD

Participants.

All data were acquired in accordance with the Institutional Review Board of the University of California, Los Angeles. This study was specifically approved by the Institutional Review Board. Parents of participants provided informed written consent prior to the start of study activities. EEG datasets analyzed for this study will be deposited in a public repository following publication of this manuscript. Participants were clinically referred through the Dup15q clinic at UCLA and the national Dup15q Alliance. Children were excluded from the study if treated with medications known to pharmacologically induce beta oscillations (benzodiazepines, benzodiazepine derivatives, or barbiturates). A total of 16 participants were recruited for the first study, 5 of whom were omitted due to treatment with exclusionary medications (n = 2), duplication type that did not include the canonical 15q11.2-q13.1 region (n = 1), or insufficient length or quality of EEG recordings (n = 2). The remaining sample included 11 participants (5 male), 16–144 months of age (median = 54 months). Details of the sample, including age, intelligence quotient (IQ), medication, and duplication type can be viewed in Table 1. A wide age range was included to ensure that a clinically representative sample was being studied, and age matching of the comparison groups ensured that the group level comparisons would not be confounded by age differences. Both isodicentric (n = 8) and interstitial (n = 3) duplications were represented in this cohort, and 2 participants with isodicentric duplications had a diagnosis of epilepsy. Data from an ongoing study of electrophysiological biomarkers in ASD were utilized for the two comparison groups: (1) an age and IQ-matched cohort of children with non-syndromic ASD (n = 10) and (2) an age-matched group of TD children (n = 9). Preschool age children with ASD were recruited as part of a larger study investigating predictors of treatment outcome in preschoolers enrolled in a UCLA early intervention program. All children enter the program with a prior clinical diagnosis of ASD, made through the California State Regional Center, independent clinical psychologists, child psychiatrist, and/or developmental pediatricians. Diagnoses were confirmed by UCLA psychologists based on DSM-IV criteria. Non-syndromic ASD was defined by normal clinical chromosomal microarray testing, but most children had not undergone whole exome sequencing. Details of both comparison groups are available in Table 1.

[Figure omitted. See PDF.]

Table 1. Dup15q syndrome participant characteristics.

Older participants tested in Orlando did not undergo cognitive testing, as age-appropriate cognitive tests were not available. N/A = not available.

https://doi.org/10.1371/journal.pone.0167179.t001

Clinical assessment.

Owing to the large range in age and developmental ability amongst participants in our study, several assessments were used to evaluate cognition, language, and motor skills. The following measures were used to match participants by cognitive function: the Mullen Scales of Early Learning (MSEL)[13], the Stanford Binet Intelligence Scales-Fifth Edition (SB5) [14], the Differential Ability Scales Second Edition (DAS-II)[15], Preschool Language Scales-Fifth Edition (PLS-5) [16], and the Leiter International Performance Scales–Revised (Leiter-R) [17].

EEG recording.

Spontaneous EEG was recorded at 500 Hz using high-density 129 channel geodesic nets with Ag/AgCl electrodes (Electrical Geodesics, Inc., Eugene, OR, USA) while participants watched nonsocial silent videos of bouncing soap bubbles and other abstract images on a computer monitor for 2 to 6 minutes, depending on the child’s level of compliance with the paradigm. No sedation was employed for electrode placement or EEG recording. EEG signals were amplified using a Net Amps 300 amplifier (Electrical Geodesics, Inc., Eugene, OR, USA) with a low-pass analog filter cutoff frequency of 6 KHz. EEG signals were vertex-referenced at the time of recording and later re-referenced to average after preprocessing and artifact reduction.

Data processing.

EEG recordings were band pass filtered at 1–50 Hz using a finite impulse response (FIR) filter with the EEGLAB toolbox [18]. Recordings were then segmented into 1000 ms segments for preprocessing. Noisy or loose channels were spherically interpolated using EEGLAB, and EEG recording segments with more than 11 interpolated channels were rejected. All remaining segments were manually inspected for non-stereotyped artifacts, e.g., electromyogram (EMG), and rejected based on qualitative inspection. Following manual artifact rejection, a combined principal component analysis (PCA) and independent component analysis (ICA) approach was used to eliminate stereotyped artifacts, e.g., ocular artifacts. PCA was performed prior to ICA to reduce each dataset to 24 dimensions, an important consideration for successful extraction of meaningful independent components (ICs) from EEG recordings [19]. ICs corresponding to physiological artifact were subtracted from EEG recordings. All EEG recordings were re-referenced to an average reference prior to power calculations.

EEG recordings from ASD and TD cohorts were acquired and processed according to the same protocol described above with the sole discrepancy that 8 ASD and 9 TD recordings were sampled at a lower frequency (250 Hz) than those with Dup15q syndrome. Accordingly, EEG signals from all groups were downsampled to 250 Hz prior to spectral analysis if originally sampled at 500 Hz to ensure that all power spectral densities (PSDs) are computed with the same frequency resolution.

Statistical analysis of spectral power.

To study spectral power, 9 regions of interest (ROIs) were defined corresponding to the locations of channels F3, Fz, F4, C3, Cz , C4, P3, Pz, P4 in the international 10–20 montage, as per prior studies by our group [20]. Each ROI consisted of 4 electrodes and was chosen to allow for maximum spatial coverage of the scalp (see Fig 1A). For each electrode, PSDs were computed according to Welch’s method [21] using the pwelch function in MATLAB. This method divides the signal into sections of equal length with 50% overlap and uses a Hamming window to estimate a modified periodogram for each segment. Segments are overlapped by 50% owing to the fact that the Hamming window weighs the center of the data segment more strongly than the sidelobes, which are attenuated by 42.5 dB. Periodograms for all segments are then averaged to provide a final spectral estimate. Relative power was then calculated for specific frequency bands, and these relative power values were then averaged across a particular ROI. 2000 ms of signal sampled at this rate were zero-padded to create 512-point Hamming windows to compute PSDs with approximately 0.5 Hz frequency resolution. The 2000 ms signal length was chosen as twice the period of the slowest oscillation examined (1 Hz delta activity). The following frequency bands were examined separately: delta (1–4 Hz), theta (4–8 Hz), alpha (8–12 Hz), beta1 (12–20 Hz), beta2 (20–30 Hz), and gamma (30–48 Hz). EEG power was calculated from PSDs as relative power, i.e., the proportion of total (1–48 Hz) spectral power accounted for by a given frequency band. Normalizing spectral power in this manner allows for meaningful comparisons between subjects with different overall levels of signal power.

[Figure omitted. See PDF.]

Fig 1. Broadband EEG recordings.

Qualitative analysis of spontaneous electroencephalogram (EEG) recordings from participants with Dup15q syndrome revealed overt beta frequency oscillations apparent upon visual inspection. (A) EEG recordings from 9 scalp regions of interest (ROIs) were analyzed: left frontal (orange), left central (yellow), left posterior (green), mid-frontal (aquamarine), mid-central (cyan), mid-posterior (blue), right frontal (purple), right central (pink), right posterior (red). (B) 3-dimenisonal head model showing ROI electrode locations. (C) 3 s of broadband EEG recordings from a representative 29-month-old TD child from 36 channels across 9 ROIs. (D) Same duration of EEG recorded from a 27-month-old child with nonsyndromic ASD. (E) EEG from a representative Dup15q syndrome participant (age 28 months) reveals spontaneous beta oscillations (SBOs) in virtually all channels and all ROIs. The overt quality of SBOs likely allows for their easy detection in clinical EEG recordings. By contrast, (F) a 43-month-old participant with both Dup15q syndrome and epilepsy does not show nearly such distinct SBOs. It is possible that beta activity is reduced in children with both Dup15q syndrome and epilepsy.

https://doi.org/10.1371/journal.pone.0167179.g001

Relative power values from 6 frequency bands and 9 ROIs were analyzed using repeated measures analysis of variance (ANOVA), with multiple testing corrected for using false discovery rates (FDR). The final model for each band included group (Dup15q, ASD, TD) and ROI main effects. Because no significant group-by-ROI interactions were found, interactions were removed from the final model. Since estimated variances for relative power in the Dup15q group were larger compared to ASD and TD groups, especially for beta power, separate group variances and within-subject correlations were allowed for in the fitted repeated measures ANOVA models. Comparisons of spectral power were performed both with and without outliers (defined by being more than three standard deviations from the mean) in order to assess the extent to which comparisons might be biased by outliers.

Study 2: Subgroup analyses of EEG power within Dup15q syndrome

Participants.

All data were acquired in accordance with the Institutional Review Board of the University of California, Los Angeles. This study was specifically approved by the Institutional Review Board of the University of California, Los Angeles. Parents of participants provided informed written consent prior to the start of study activities. After this initial study, we expanded the cohort in order to examine the variability in beta power within Dup15q syndrome. Spontaneous EEG data were acquired from individuals at the 2015 National Dup15q Alliance conference in Orlando, Florida, with consenting procedures as described in the first study. EEG datasets analyzed for this study will be deposited in a public repository following publication of this manuscript.

EEG recording and processing.

EEG recording followed the protocol described in the initial cohort study. A total of 24 participants underwent EEG recordings at the Orlando site, with 9 participants omitted due to treatment with exclusionary medications or insufficient length or quality of EEG recordings. Mean beta power did not differ between testing sites using Welch two-sample t-tests (beta1 power: t = 0.41, p = 0.69; beta2 power: t = 0.68, p = 0.50). Therefore, the data acquired from the Orlando site were combined with the 11 participants in the initial Dup15q cohort and 1 additional adult tested at UCLA, resulting in a final sample of 27 participants, 16–384 months of age (median = 81.2 months). As with the first study, no sedation was employed for electrode placement or EEG recording. Table 1 provides details about this cohort. EEG data were processed according to the protocol described in the initial study.

Statistical analysis of spectral power.

To model the effects of age, duplication type, and epilepsy on beta power, several simple linear regression models were implemented using two outcome measures, beta1 power and beta2 power, averaged across all ROIs. Simple linear regressions were used, with each variable as a separate predictor of beta1 and beta2 power. Age was treated as a continuous variable while duplication type and epilepsy were treated as binary variables. Prior to regression modeling, Welch two-sample t-tests were used to test for differences in mean beta power between testing sites. Having tested the null hypothesis that means of outcome measures do not differ between testing sites, three univariate regressions were performed for beta1 power and beta2 power using the aforementioned predictors.

Results

Study 1 –Comparison of Dup15q syndrome with ASD and TD

Results of behavioral testing, along with duplication type and epilepsy history, are summarized for all participants in Table 1. Qualitative analysis of bandpass filtered EEG traces from several participants with Dup15q syndrome revealed fast sinusoidal oscillations at beta frequencies apparent upon visual inspection. Fig 1 shows beta activity from 9 scalp regions of interest (ROIs) in a 29-month old TD child (Fig 1C), a 27-month-old child with nonsyndromic ASD (Fig 1D), a 28-month-old child with Dup15q syndrome (Fig 1E), and a 43-month old participant with both Dup15q syndrome and epilepsy (Fig 1F). ROIs were selected for maximum scalp coverage and their correspondence to 10–20 montage channels [20]. Visual inspection of topographic scalp plots for both beta1 (12–20 Hz) relative power and beta2 (20–30 Hz) relative power averaged across participants revealed a diffuse pattern of beta activity in Dup15q syndrome that appeared strongest over frontotemporal regions (Fig 2), in stark contrast to both comparison groups. Frontotemporal distributions of beta power in Dup15q syndrome are best visualized by projection onto 3-dimensional head models (Fig 2D).

[Figure omitted. See PDF.]

Fig 2. Group averaged EEG topoplots.

Topographic scalp plots of relative beta1 (12–20 Hz) power (row 1) and relative beta2 (20–30 Hz) power (row 2) averaged across all participants. Three groups are compared, Dup15q syndrome group (A), nonsyndromic ASD group (B), and TD group (C). Qualitative inspection of these scalp plots reveals profound differences in relative power between Dup15q syndrome and comparison groups. This contrast is most pronounced in the beta2 band. Beta power in both subbands is most pronounced over frontotemporal scalp regions and attenuated over the midline. This is best visualized in the Dup15q syndrome cohort by projecting averaged power onto a generic 3-dimensional head model (D).

https://doi.org/10.1371/journal.pone.0167179.g002

Because no significant group-by-ROI interactions were found in modeling power in each frequency band, group means were averaged across regions, and comparisons were made between (1) the Dup15q cohort and the ASD cohort and (2) the Dup15q cohort and the TD cohort, resulting in a total of 12 frequency comparisons (2 pairwise comparisons for 6 frequency bands, see Table 2). Group comparisons were robust to the removal of outliers greater than 3 standard deviations from the mean, and thus outliers were included in the final analysis.

[Figure omitted. See PDF.]

Table 2. Results from contrast analysis following repeated measures ANOVA.

12 pairwise comparisons were performed following the repeated measures analysis of variance (ANOVA), one comparison of duplication 15q11.2-q13.1 (Dup15q) syndrome to the nonsyndromic ASD group and one comparison of Dup15q syndrome to the TD group for each of 6 frequency bands. Mean estimate is the mean difference between groups. All comparisons were corrected for multiple testing using false discovery rates (FDR). Statistically significant comparisons appear in bold.

https://doi.org/10.1371/journal.pone.0167179.t002

Effect sizes were computed from raw data averaged across ROIs using Cohen’s d. Large effect sizes were found in both beta1 power (Dup15q –ASD, d = 1.08; Dup15q –TD, d = 1.12) and beta2 power (Dup15q –ASD, d = 1.73; Dup15q –TD, d = 1.63). These are especially large effect sizes for studies of EEG power in neurodevelopmental disorder. The Dup15q cohort was best distinguished from comparison groups in the beta2 band, in which participants exhibited significantly stronger spontaneous power than both the ASD (p = 6 × 10⁻⁴, FDR corrected) and TD (p = 6 × 10⁻⁴, FDR corrected) groups (see Table 2 and Fig 3). Children with Dup15q syndrome also exhibited stronger spontaneous beta1 power (Fig 3) than both the ASD group (p = 0.0252, FDR corrected) and the TD group (p = 0.02, FDR corrected). Two other significant findings included weaker spontaneous delta power in Dup15q syndrome compared to both the ASD group (p = 0.03, FDR corrected; d = -1.07) and the TD group (p = 0.02, FDR corrected; d = -1.13) and higher spontaneous gamma power in Dup15q syndrome compared to the ASD group (p = 0.04, FDR corrected; d = 0.97). Spontaneous gamma power was not significantly different between children with Dup15q syndrome and the TD group.

[Figure omitted. See PDF.]

Fig 3. Distributions of relative beta1 and beta2 power for all groups.

Dot plots of relative beta1 (12–20 Hz) power (A) and beta2 (20–30 Hz) power (B) averaged across all ROIs in the Dup15q syndrome group (left), nonsyndromic ASD group (center), and TD group (right). Participants with interstitial duplications [int Dup(15)] are colored red, while participants with isodicentric duplications [idic Dup(15)] are colored blue. Comparison group participants with nonsyndromic ASD are colored green, while children in the TD comparison group are colored purple. The Dup15q cohort features greater mean and standard deviation (mean ± S.D.) in beta1 power (0.060 ± 0.031) and beta2 power (0.092 ± 0.051) in relation to comparison groups [beta1 power: 0.034 ± 0.010 (ASD), 0.034 ± 0.0052 (TD); beta2 power: 0.027 ± 0.012 (ASD), 0.030 ± 0.0070 (TD)].

https://doi.org/10.1371/journal.pone.0167179.g003

These data suggest that a pattern of high beta power and low delta power distinguishes Dup15q syndrome from both children with nonsyndromic ASD and TD children, with the largest effect size found in the beta2 band. High gamma power may also distinguish Dup15q syndrome from nonsyndromic ASD. We next asked if gamma power is related to beta2 power, possibly as an extension of the same broadband electrophysiological activity, by adding EEG recordings from Study 2 to examine a larger Dup15q syndrome cohort (n = 27). In fact, beta2 power and gamma power were significantly correlated in every cohort (Dup15q syndrome, r = 0.43, p = 0.024; non-syndromic ASD, r = 0.78, p = 0.0080; TD, r = 0.88, p = 0.0017), suggesting that gamma power and beta2 power reflect the same electrophysiological process. Given the reciprocal pattern of high beta power and low delta power in Dup15q syndrome, we then asked a similar question about how these two frequency bands might be related. We observed that delta power and beta power are significantly related only in the Dup15q cohort (beta1, r = -0.63, p = 5 × 10⁻³; beta2, r = -0.47, p = 0.013) and are thus possibly codependent features of the same EEG signature. This relationship between delta power and beta power was not significant in either comparison group, but appeared as a trend in the ASD cohort (beta1, r = -0.62, p = 0.056; beta2, r = -0.62, p = 0.056).

While spectral power measures are integrated across frequency bins, they do not capture peaks and valleys in the frequency domain. To capture such spectral features, grand averaged PSDs were visualized in Fig 4A by computing the mean across ROIs and all participants. A clear spectral peak in the beta band can be seen for the Dup15q syndrome cohort studied at UCLA in Fig 4A. Fig 4B shows the presence of beta spectral peaks in PSDs averaged across ROIs for individual participants with Dup15q syndrome. After removing the 1/f trends (dotted lines) that account for most variance in EEG PSDs (Fig 4C), the peak frequency for participants with Dup15q syndrome (~23 Hz) was higher than that of both the ASD (~8 Hz) and TD (~9 Hz) comparison groups (Fig 4D).

[Figure omitted. See PDF.]

Fig 4. Grand averaged power spectral densities from all groups.

(A) Power spectral densities (PSDs) averaged across all regions of interest (ROIs) and participants for the Dup15q syndrome group (red), nonsyndromic ASD group (blue), and TD group (green). Before averaging, participant PSDs are normalized such that the area under the curve equals 1 to emphasize relative power. Translucent highlights represent standard error of the mean (SEM) computed across participants. An enormous peak from 12–30 Hz reveals the presence of powerful spontaneous beta oscillations (SBOs) in the Dup15q cohort. PSDs are normalized to represent relative power. (B) Individual PSDs, averaged across ROIs, from participants with Dup15q syndrome. (C) Group averaged linear trends (dotted lines) fitted from log-log transformed PSDs. Linear trends represent the 1/f distribution inherent in the EEG. (D) Group averaged PSDs with linear trends removed to emphasize deviations from the 1/f trend. Dup15q syndrome shows the largest deviation, with a peak frequency (~23 Hz) in the beta band. Both comparison groups feature peak frequencies in the alpha band.

https://doi.org/10.1371/journal.pone.0167179.g004

Study 2 –Within-group analysis

Of the regression models tested, only epilepsy diagnosis statistically predicted beta2 power, with stronger beta2 power in participants with Dup15q syndrome who did not have epilepsy (R² = 0.17, p = 0.03; Fig 5). Qualitative evidence for this finding can be seen in averaged scalp plots of beta power in individuals with Dup15q with and without epilepsy (Fig 5B and 5D). Neither age nor duplication type significantly predicted beta power within the Dup15q cohort. Because between group differences were also found in delta power and gamma power, we asked if these variables could also predict epilepsy status. Neither delta power nor gamma power significantly predicted epilepsy status in Dup15q syndrome, although a trend level finding is observed for gamma power (R² = 0.12, p = 0.076).

[Figure omitted. See PDF.]

Fig 5. Age, duplication type, and epilepsy as predictors of beta1 and beta2 power.

(A) Scatter plots of age, duplication type, and epilepsy against relative beta1 (12–20 Hz) power. Duplication type and epilepsy are treated as binary variables. Interstitial duplications and no epilepsy are represented as 0; Isodicentric duplications and epilepsy are represented as 1. Highlighted area around regression line represents the 95% confidence region. (B) Topographic scalp plots of beta1 power averaged across participants with epilepsy (left) and without epilepsy (right). (C) Scatter plots of age, duplication type, and epilepsy against relative beta2 (20–30 Hz) power. (D) Topographic scalp plots of beta2 power averaged across participants with epilepsy (left) and without epilepsy (right). The relationship between epilepsy and beta2 power is statistically significant (R² = 0.17, p = 0.032).

https://doi.org/10.1371/journal.pone.0167179.g005

Discussion

Dup15q syndrome is highly penetrant for intellectual disability, epilepsy, and ASD. Several clinical reports have described a distinctive feature on clinical EEG that may represent an electrophysiological biomarker of this syndrome in the form of increased beta oscillations. We found that EEG beta power (SBOs) most strongly distinguished children with Dup15q syndrome from both (1) age/IQ-matched children with nonsyndromic ASD and (2) age-matched TD children. Changes in the EEG beta band were qualitatively obvious upon visual inspection of data as SBOs. These SBOs, as measured by power in the beta1 and beta2 bands, correlate with epilepsy diagnosis but not age or duplication type. Although our clinically referred sample included a small cohort with a relatively large age range, the robustness of this EEG signature in both this cohort and the larger sample examined in the second study is evident by the large effect sizes (|d| > 1).

The promise of biomarkers

There has been a tremendous interest in the field of neurodevelopmental disorders in the identification of quantitative measures of brain function that may relate to specific genetic etiologies, as these “biomarkers” can help provide clues into the neurobiological sequelae of a genetic variation (linking genes to brain) and shed light on the impact of aberrant brain function on behavior (linking brain to behavior). A quantitative measure of neural function in a genetically defined subgroup may provide a more refined assay of subtle individual differences that can inform predictors of outcome, particularly in the context of interventions that target the specific mechanism underlying the measure. In other words, one may see a change in a biomarker with treatment that precedes any overt behavioral change but that suggests engagement of the biological target and, therefore, hope for clinical improvement.

For both practical and scientific reasons, EEG is a particularly robust method to measure neural function in developmental disorders. Not only does it have excellent motion tolerance, but its temporal resolution also allows it to resolve neurophysiological oscillations and dynamics on a millisecond scale. Power Spectral Densities (PSDs) from EEG recordings follow a characteristic 1/f distribution, so named because spectral power is inversely proportional to frequency. 1/f distributions are ubiquitous in the brain and are a likely signature of balance between neural excitation and inhibition (E/I balance) [22, 23].

Presumed mechanisms of spontaneous oscillations in Dup15q syndrome

Enhanced SBOs and diminished delta oscillations observed in Dup15q syndrome represent deviations from the 1/f distribution (Fig 4C). Because E/I balance is believed to be necessary for varied and complex electrophysiological signals [22, 24], deviations from the 1/f distribution likely represent a disruption of balanced neurotransmission. In Dup15q syndrome, a disruption of E/I balance could be created by GABA_AR subunit gene overexpression. SBOs observed in Dup15q syndrome [6, 7] strongly resemble those induced by positive allosteric modulators (PAMs) of GABA_ARs such as benzodiazepines and barbiturates [25–27]. Benzodiazepines, barbiturates, and other GABA_AR PAMs increase the net chloride flux through the GABA_AR’s ion pore [28, 29]. Barbiturates [28, 30] and at least one benzodiazepine compound, zolpidem [31, 32], have been shown to increase the time constant of GABA_ARs by lengthening the duration of hyperpolarizing chloride currents through the receptor’s ion pore. Beta activity resulting from GABA_AR modulation or other dysfunction may actually represent slowed gamma activity, as the GABA_AR time constant is known to control the frequency of gamma oscillations [24]. In the healthy brain, inhibitory interneurons with reciprocal connections fire synchronously to inhibit pyramidal cells, silencing themselves in the process. Pyramidal cells recover from inhibition to enjoy a period of excitability before being silenced again by interneurons, beginning the gamma cycle anew [24, 33]. Lengthening the time constant of GABA_ARs through altered GABA_AR gene expression would lengthen the period of the gamma cycle to that characteristic of beta oscillations. Together, high-frequency beta and gamma oscillations are hypothesized to play a critical role in temporal binding of local circuits during cognitive tasks.

In addition to elevated beta power, spectral anomalies were also observed in the delta and gamma bands of EEG recordings from participants with Dup15q syndrome. It is possible that reduced delta power observed in Dup15q syndrome may be linked to enhanced beta power in a reciprocal manner. A significant negative correlation was observed between delta power and both beta1 power and beta2 power only in the Dup15q syndrome cohort. Angelman syndrome, most commonly caused by maternal deletion of the 15q11-q13 region, including UBE3A, features enhanced delta oscillations in clinical EEG recordings [34, 35], suggesting a reciprocal relationship between deletion and duplication of the GABA_AR subunits. Furthermore, loss of UBE3A has been associated with enhanced delta oscillations [36] and suppression of ventral striatal GABA co-release in mouse models of Angelman syndrome [37], underscoring the relationship between the ubiquitin ligase and GABAergic transmission.

Evidence from pharmacological studies also suggests that reduced delta power in Dup15q syndrome may be directly related to GABA_AR subunit gene overexpression rather than UBE3A per se. For instance, the benzodiazepine compounds diazepam and zolpidem decrease cortical EEG delta power in awake, behaving rats [27, 38]. Similarly, midazolam has been shown to reduce EEG delta power during sleep in rats [39]. In humans with generalized anxiety disorder, the benzodiazepine clorazepate has been shown to reduce delta power in scalp EEG recordings [40]. All studies considered so far also associated increased beta power with benzodiazepine challenge. The foregoing evidence from both rodents and humans suggests that GABA_AR potentiation and EEG delta power are inversely related.

Finally, we suggest that our finding of stronger gamma oscillations in Dup15q syndrome as compared with nonsyndromic ASD may reflect a common mechanism for beta2 and gamma oscillations involving feedback inhibition between pyramidal cells and interneurons. This hypothesis is supported by the fact that beta2 power and gamma power were strongly correlated in all three cohorts. Nonetheless, the effect size observed in beta2 (d = 1.73) was greater than that of gamma (d = 0.97), suggesting that SBOs are a clearer biomarker of Dup15q syndrome than gamma oscillations.

SBOs as markers of gene expression?

Although we did not measure UBE3A and GABA receptor gene expression from participants in our study, these data support the need for future investigations that can directly examine the relationship between EEG power and mRNA transcript levels from GABRA5, GABRB3, GABRG3, and UBE3A, all genes which are duplicated in Dup15q syndrome. In particular, there already exists evidence from pharmacological studies of GABA_AR PAMs (i.e., benzodiazepines) [25–27], as well as correlations between motor evoked beta power and resting GABA levels [41], suggesting an important relationship between beta activity and GABAergic transmission. As with beta power in our study, Scoles et al. found greater mean and variance in neural GABRB3 expression in a small cohort (n = 8) of postmortem tissue samples from individuals with Dup15q syndrome (isodicentric duplications) compared to nonsyndromic ASD and TD tissue samples [42]. The close resemblance between the distribution of GABRB3 expression in the Scoles et al. study and the distribution of beta power in our study can be visualized in Fig 4B of Scoles et al. (2011) (cf. Fig 3, this paper). Slight overlap in distributions of beta power between the Dup15q syndrome cohort and the ASD comparison group could potentially reflect point mutations of GABA_AR subunit genes in some children in the ASD comparison group. Finally, many studies of benzodiazepine GABA_AR PAMs have shown reductions of delta power in a variety of contexts [27, 38–40, 43, 44], suggesting that reduced delta power in the wakeful spontaneous EEG of Dup15q syndrome participants could also be explained by GABA_AR abnormalities such as altered subunit expression. Further work in humans and animal models will be necessary to test this hypothesis.

Another possible cause of reduced EEG delta power in Dup15q syndrome is duplications of UBE3A, the causative gene of Angelman syndrome. Patients with Angelman syndrome do not express UBE3A in neural tissue and show elevated EEG delta power [34, 35, 37], the opposite electrophysiological phenotype as Dup15q syndrome. For this reason, it is plausible that an inverse relationship exists between EEG delta power and UBE3A expression levels. It is also possible that UBE3A overexpression influences beta power, perhaps an indirect effect mediated through the GABAergic system. For instance, a recent study in a UBE3A-null mouse model found that loss of UBE3A can impair co-release of subcortical GABA, thus demonstrating their intrinsic functional relationship [37]. However, considering that SBOs have been reported in individuals with paternal Dup15q syndrome [6, 7], the paternally imprinted UBE3A alone cannot explain SBOs in Dup15q syndrome.

SBOs and epilepsy risk

We identified an inverse relationship between SBOs, particularly beta2 power, and epilepsy in Dup15q syndrome in the larger cohort analysis. At first glance, this result could be interpreted as SBOs being markers of enhanced GABAergic tone (as found in GABAergic medications such as benzodiazepenes), thus providing neural protection against seizures in this population. However, interpretation of this relationship requires caution, as likely there are modifying factors such as background EEG, antiepileptics, and developmental level of individuals with epilepsy that may influence this relationship [45]. In particular, with our small sample, we could not fully disentangle the relationship between age and epilepsy. To better understand the relationship between age, duplication type, epilepsy status, and beta2 power, we visualized all 4 variables in Fig 6. Log-transformed age (abscissa) is plotted against beta2 power (ordinate), with interstitial duplications represented by diamonds and isodicentric duplications represented by circles. Children are color coded by epilepsy status (pink for epilepsy, blue for no epilepsy). As seen in this figure, those participants with the highest beta2 power included our youngest children without epilepsy. Of note, while it is difficult to separate the role of epilepsy per se from the impact of medications on EEG, anticonvulsant medications typically increase beta activity rather than decrease such activity [46]. Although no subjects in this study were treated with benzodiazepine or barbiturate medication, most children with Dup15q syndrome with epilepsy were treated with levetiracetam, which has been shown to increase (rather than decrease) relative beta power in epilepsy patients [46].

[Figure omitted. See PDF.]

Fig 6. Dup15q syndrome participants by age, beta2 power, duplication type, and epilepsy status.

A scatter plot of age versus beta2 power in Dup15q syndrome reveals a cluster of very young participants with very high beta2 power (top left corner). None of these participants have epilepsy. All participants with epilepsy feature beta2 power < 0.08. Note that the abscissa has been log-transformed to accommodate a large number of young participants and a much smaller number of older participants.

https://doi.org/10.1371/journal.pone.0167179.g006

We must emphasize that the relationship between GABAergic activity and epilepsy is complex as is, most likely, the relationship between beta oscillations and epilepsy. It may seem paradoxical that a syndrome associated with overexpression of GABA receptor genes could also confer such a high risk for epilepsy. In general, synaptic, i.e. phasic, inhibition is mediated by α1 or α2 containing GABA receptors that also have a γ2 subunit producing an increase in chloride conductance favoring hyperpolarization. The extrasynaptic GABA receptors, usually α4 or α6 combined with δ (no gamma) are sensitive to small changes in ambient GABA, and they produce sustained lowering of the membrane potential [47]. Excessive GABAergic activity can produce two effects that are both potentially epileptogenic. First, the spike-wave bursts in the thalamo-cortical network are initiated by hyperpolarization, since the lowered membrane potential is the trigger for the low-threshold calcium currents [48]. Secondly, prolongation of inhibition can promote synchronization of such networks, which can occur with enhancement of tonic inhibition [49]. Thus spike-wave stupor can be treated with intravenous benzodiazepines which act only at synaptic receptors. On the other hand, generalized spike-wave paroxysms become more frequent and longer in duration when ambient GABA is increased by other medications.

Limitations and future directions

The relatively small sample size, albeit representative of the full clinical spectrum of this rare disorder, does limit further subgroup analyses and examination of clinical correlates of this biomarker. The sample size also precludes the development of a multi-variable prediction model of these elevated beta oscillations. The association with epilepsy, in particular, warrants further investigation through two approaches. First, through longitudinal studies we can examine changes in EEG power after the onset of epilepsy and, by doing so, elucidate whether beta oscillations represent a protective biomarker for the development of seizures.

However, these findings have laid the foundation for a larger scale study of the functional and clinical implications of electrophysiological biomarkers in this syndrome. Through a multi-site, coordinated effort with the National Dup15q Alliance, we will expand our sample size to ask the following questions: First, do changes in state modulate beta power, in particular during cognitive or perceptual tasks, or during sleep? One might hypothesize that persistent beta power in sleep could disrupt sleep architecture enough to impact cognition and behavior in these children with neurodevelopmental disabilities. Moreover, a lack of modulation of EEG oscillations during cognitive tasks could directly hinder learning. We will directly examine the relationship between beta power and changes in beta power with more quantitative measures of cognition and autism severity. Second, what are the exact genetic underpinnings of this biomarker? We will examine electrophysiological markers in several pre-clinical models of Dup15q syndrome (full genetic duplication compared to UBE3A overexpression mouse). Translational studies linking gene expression and SBOs in individuals with Dup15q syndrome will elucidate the specific role of UBE3A, GABRA5, GABRB3, and GABRG3 gene expression on this biomarker and help further understand the genetic effects on transcript levels that lead to the pathogenesis of Dup15q syndrome.

Conclusions

The field of ASD research has sorely lacked quantifiable biomarkers that may help parse the neurobiological heterogeneity of this spectrum of disorders. Moreover, as genetic testing in ASD and related neurodevelopmental disorders becomes clinical gold standard, an increasing number of children are diagnosed with genetic variants that not only will elucidate causal mechanisms but also therapeutic targets [50]. The identification of these targets necessitates quantifiable biomarkers that relate directly to genetic mechanisms. Studies in Dup15q syndrome provide a promising path towards this mission, as the elucidation and quantification of an electrophysiological biomarker could improve diagnosis and prognostication, as well as measurement of target engagement and outcomes in clinical trials, a model that can inform similar investigations in the quickly expanding number of high-risk genetic syndromes associated with neurodevelopmental disorders.

Acknowledgments

We warmly thank the children and families who participated in this research for making this work possible. We especially thank Andrei Irimia and Iman Mohammad-Rezazadeh for their assistance with EEG signal processing and Amaya Miquelajauregui for giving her thoughts and insights on this manuscript. Joel Frohlich was supported by the National Institute of Mental Health of the National Institutes of Health under Award Number T32MH073526. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health This work was also funded by grant 20144605 from the Dup15q Alliance, to whom we extend our sincere gratitude, as well as R01 GM111378-01A1 from the National Institute of General Medical Sciences NIGMS (NIH), NIMH K23MH094517-01, and IDDRC U54HD087101.

Author Contributions

1. Conceptualization: SSJ JF.

2. Data curation: JF SSJ SH.

3. Formal analysis: DS JF SSJ PG RS VS.

4. Funding acquisition: SSJ PG.

5. Investigation: JF CD SSJ SH VS.

6. Methodology: SSJ DS JF.

7. Project administration: SSJ.

8. Resources: SSJ DS.

9. Software: JF.

10. Supervision: SSJ EHC.

11. Validation: SSJ JF SH.

12. Visualization: JF SH.

13. Writing – original draft: JF SSJ DS LTR PG.

14. Writing – review & editing: PG LTR RS RLT CD EHC SSJ JF DS VS.

Citation: Frohlich J, Senturk D, Saravanapandian V, Golshani P, Reiter LT, Sankar R, et al. (2016) A Quantitative Electrophysiological Biomarker of Duplication 15q11.2-q13.1 Syndrome. PLoS ONE 11(12): e0167179. https://doi.org/10.1371/journal.pone.0167179

References

1. Schaefer GB, Mendelsohn NJ, Professional P, Guidelines C. Clinical genetics evaluation in identifying the etiology of autism spectrum disorders: 2013 guideline revisions. Genetics in Medicine. 2013;15(5):399–407. pmid:23519317

2. Gaugler T, Klei L, Sanders SJ, Bodea CA, Goldberg AP, Lee AB, et al. Most genetic risk for autism resides with common variation. Nature genetics. 2014;46(8):881–5. pmid:25038753

3. Tammimies K, Marshall CR, Walker S, Kaur G, Thiruvahindrapuram B, Lionel AC, et al. Molecular diagnostic yield of chromosomal microarray analysis and whole-exome sequencing in children with autism spectrum disorder. JAMA. 2015;314(9):895–903. pmid:26325558

4. Jeste SS, Frohlich J, Loo SK. Electrophysiological biomarkers of diagnosis and outcome in neurodevelopmental disorders. Curr Opin Neurol. 2015.

5. Cook EH Jr, Lindgren V, Leventhal BL, Courchesne R, Lincoln A, Shulman C, et al. Autism or atypical autism in maternally but not paternally derived proximal 15q duplication. American journal of human genetics. 1997;60(4):928. pmid:9106540

6. Urraca N, Cleary J, Brewer V, Pivnick EK, McVicar K, Thibert RL, et al. The interstitial duplication 15q11.2-q13 syndrome includes autism, mild facial anomalies and a characteristic EEG signature. Autism Res. 2013;6(4):268–79. pmid:23495136

7. Al Ageeli E, Drunat S, Delanoë C, Perrin L, Baumann C, Capri Y, et al. Duplication of the 15q11-q13 region: clinical and genetic study of 30 new cases. European journal of medical genetics. 2014;57(1):5–14. pmid:24239951

8. DiStefano C, Gulsrud A, Huberty S, Kasari C, Cook E, Reiter LT, et al. Identification of a distinct developmental and behavioral profile in children with Dup15q syndrome. Journal of neurodevelopmental disorders. 2016;8(1):1.

9. Hogart A, Wu D, LaSalle JM, Schanen NC. The comorbidity of autism with the genomic disorders of chromosome 15q11.2-q13. Neurobiology of developmental disorders involving human higher cognition. 2010;38(2):181–91.

10. Greer PL, Hanayama R, Bloodgood BL, Mardinly AR, Lipton DM, Flavell SW, et al. The Angelman Syndrome protein Ube3A regulates synapse development by ubiquitinating arc. Cell. 2010;140(5):704–16. pmid:20211139

11. Smith SEP, Zhou Y-D, Zhang G, Jin Z, Stoppel DC, Anderson MP. Increased gene dosage of Ube3a results in autism traits and decreased glutamate synaptic transmission in mice. Science translational medicine. 2011;3(103):103ra97–ra97. pmid:21974935

12. Dindot SV, Antalffy BA, Bhattacharjee MB, Beaudet AL. The Angelman syndrome ubiquitin ligase localizes to the synapse and nucleus, and maternal deficiency results in abnormal dendritic spine morphology. Human molecular genetics. 2008;17(1):111–8. pmid:17940072

13. Mullen EM. Infant Mullen Scales of Early Learning (Infant MSEL): TOTAL Child, Incorporated; 1989 1989.

14. Roid GH. Stanford-binet intelligence scales (; SB5). Rolling Meadows, IL: Riverside. 2003.

15. Elliot CD. Differential Abilitity Scales (2nd Ed). San Antonio: Harcourt Assessment; 2007 2007.

16. Zimmerman IL, Steiner VG, Pond RE. Preschool language scale, (PLS-5). Age. 2011.

17. Roid G. Miller L. Leiter International Performance Scale–Revised. 1997.

18. Delorme A, Makeig S. EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J Neurosci Methods. 2004;134(1):9–21. pmid:15102499

19. Onton J, Westerfield M, Townsend J, Makeig S. Imaging human EEG dynamics using independent component analysis. Neurosci Biobehav Rev. 2006;30(6):808–22. pmid:16904745

20. McEvoy K, Hasenstab K, Senturk D, Sanders A, Jeste SS. Physiologic artifacts in resting state oscillations in young children: methodological considerations for noisy data. Brain Imaging Behav. 2015;9(1):104–14. pmid:25563227

21. Welch PD. The use of fast Fourier transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms. IEEE Transactions on audio and electroacoustics. 1967;15(2):70–3.

22. Poil S-S, Hardstone R, Mansvelder HD, Linkenkaer-Hansen K. Critical-state dynamics of avalanches and oscillations jointly emerge from balanced excitation/inhibition in neuronal networks. The Journal of Neuroscience. 2012;32(29):9817–23. pmid:22815496

23. Voytek B, Kramer MA, Case J, Lepage KQ, Tempesta ZR, Knight RT, et al. Age-Related Changes in 1/f Neural Electrophysiological Noise. The Journal of neuroscience. 2015;35(38):13257–65. pmid:26400953

24. Buzsaki G. Rhythms of the Brain. New York: Oxford University Press; 2006 2006.

25. Domino EF, French J, Pohorecki R, Galus CF, Pandit SK. Further observations on the effects of subhypnotic doses of midazolam in normal volunteers. Psychopharmacol Bull. 1989;25(3):460–5. pmid:2626519

26. Feshchenko VA, Veselis RA, Reinsel RA. Comparison of the EEG effects of midazolam, thiopental, and propofol: the role of underlying oscillatory systems. Neuropsychobiology. 1997;35(4):211–20. pmid:9246224

27. van Lier H, Drinkenburg WHIM, van Eeten YJW, Coenen AML. Effects of diazepam and zolpidem on EEG beta frequencies are behavior-specific in rats. Neuropharmacology. 2004;47(2):163–74. pmid:15223295

28. Löscher W, Rogawski MA. How theories evolved concerning the mechanism of action of barbiturates. Epilepsia. 2012;53(s8):12–25.

29. Nutt DJ, Malizia AL. New insights into the role of the GABAA—benzodiazepine receptor in psychiatric disorder. The British Journal of Psychiatry. 2001;179(5):390–6.

30. Macdonald RL, Twyman RE, editors. Biophysical properties and regulation of GABA A receptor channels. Seminars in Neuroscience; 1991 1991: Elsevier.

31. Fain G. Molecualr and Cellular Physiology of Neurons. Eastern Economy Edition ed. New Delhi: Prentice-Hall of India; 1999 1999.

32. Mody I, De Koninck Y, Otis TS, Soltesz I. Bridging the cleft at GABA synapses in the brain. Trends Neurosci. 1994;17(12):517–25. pmid:7532336

33. Buzsáki G, Wang X-J. Mechanisms of gamma oscillations. Annu Rev Neurosci. 2012;35:203–25. pmid:22443509

34. Robinson AA, Goldman S, Barnes G, Goodpaster L, Malow BA. Electroencephalogram (EEG) Duration Needed to Detect Abnormalities in Angelman Syndrome Is 1 Hour of Overnight Recording Sufficient? J Child Neurol. 2015;30(1):58–62. pmid:24820335

35. Muthaffar O, Minassian B, editors. Angelman Syndrome: review of clinical and EEG findings. Annual Meeting of American Epilepsy Society; 2015 2015. Philadelphia, PA.

36. Judson MC, Wallace ML, Sidorov MS, Burette AC, Gu B, van Woerden GM, et al. GABAergic Neuron-Specific Loss of Ube3a Causes Angelman Syndrome-Like EEG Abnormalities and Enhances Seizure Susceptibility. Neuron. 2016;90(1):56–69. pmid:27021170

37. Berrios J, Stamatakis AM, Kantak PA, McElligott ZA, Judson MC, Aita M, et al. Loss of UBE3A from TH-expressing neurons suppresses GABA co-release and enhances VTA-NAc optical self-stimulation. Nature Communications. 2016;7.

38. Nickel B, Szelenyi I. Comparison of changes in the EEG of freely moving rats induced by enciprazine, buspirone and diazepam. Neuropharmacology. 1989;28(8):799–803. pmid:2571102

39. Lancel M, Cro TA, Mu P, Holsboer F. Pregnenolone enhances EEG delta activity during non-rapid eye movement sleep in the rat, in contrast to midazolam. Brain research. 1994;646(1):85–94. pmid:8055343

40. Buchsbaum MS, Hazlett E, Sicotte N, Stein M, Wu J, Zetin M. Topographic EEG changes with benzodiazepine administration in generalized anxiety disorder. Biological psychiatry. 1985;20(8):832–42. pmid:2862924

41. Gaetz W, Edgar JC, Wang D, Roberts TP. Relating MEG measured motor cortical oscillations to resting γ-aminobutyric acid (GABA) concentration. Neuroimage. 2011;55(2):616–21. pmid:21215806

42. Scoles HA, Urraca N, Chadwick SW, Reiter LT, LaSalle JM. Increased copy number for methylated maternal 15q duplications leads to changes in gene and protein expression in human cortical samples. Molecular Autism. 2011;2:19–. pmid:22152151

43. Poyares D, Guilleminault C, Ohayon MM, Tufik S. Chronic benzodiazepine usage and withdrawal in insomnia patients. Journal of psychiatric research. 2004;38(3):327–34. pmid:15003439

44. Wisor JP, DeLorey TM, Homanics GE, Edgar DM. Sleep states and sleep electroencephalographic spectral power in mice lacking the β 3 subunit of the GABA A receptor. Brain research. 2002;955(1):221–8.

45. Pillai J, Sperling MR. Interictal EEG and the diagnosis of epilepsy. Epilepsia. 2006;47(s1):14–22.

46. Veauthier J, Haettig H, Meencke H-J. Impact of levetiracetam add-on therapy on different EEG occipital frequencies in epileptic patients. Seizure. 2009;18(6):392–5. pmid:19269864

47. Mody I. Distinguishing between GABA(A) receptors responsible for tonic and phasic conductances. Neurochemical research. 2001;26(8–9):907–13. Epub 2001/11/09. pmid:11699942

48. Crunelli V, Leresche N. Childhood absence epilepsy: genes, channels, neurons and networks. Nature reviews Neuroscience. 2002;3(5):371–82. Epub 2002/05/04. pmid:11988776

49. Errington AC, Di Giovanni G, Crunelli V, Cope DW. mGluR control of interneuron output regulates feedforward tonic GABAA inhibition in the visual thalamus. The Journal of neuroscience: the official journal of the Society for Neuroscience. 2011;31(23):8669–80. Epub 2011/06/10. PubMed Central PMCID: PMCPMC3130900.

50. Jeste SS, Geschwind DH. Clinical trials for neurodevelopmental disorders: At a therapeutic frontier. Science translational medicine. 2016;8(321):321fs1–fs1. pmid:26764154