- 1Department of Neurology and Developmental Medicine, Kennedy Krieger Institute, Baltimore, MD, USA
- 2Department of Neurology, Johns Hopkins University School of Medicine, Baltimore, MD, USA
- 3Department of Psychological and Brain Sciences, Johns Hopkins University, Baltimore, MD, USA
- 4Human Motor Control Section, Medical Neurology Branch, National Institute of Neurological Disorders and Stroke, National Institutes of Health, Bethesda, MD, USA
- 5Center for Neurodevelopmental and Imaging Research, Kennedy Krieger Institute, Baltimore, MD, USA
- 6Department of Psychiatry and Behavioral Sciences, Johns Hopkins University School of Medicine, Baltimore, MD, USA
Autism spectrum disorders (ASD) are thought to result in part from altered cortical excitatory-inhibitory balance; this pathophysiology may impact the generation of oscillations on electroencephalogram (EEG). We investigated premotor-parietal cortical physiology associated with praxis, which has strong theoretical and empirical associations with ASD symptomatology. Twenty five children with high-functioning ASD (HFA) and 33 controls performed a praxis task involving the pantomiming of tool use, while EEG was recorded. We assessed task-related modulation of signal power in alpha and beta frequency bands. Compared with controls, subjects with HFA showed 27% less left central (motor/premotor) beta (18–22 Hz) event-related desynchronization (ERD; p = 0.030), as well as 24% less left parietal alpha (7–13 Hz) ERD (p = 0.046). Within the HFA group, blunting of central ERD attenuation was associated with impairments in clinical measures of praxis imitation (r = −0.4; p = 0.04) and increased autism severity (r = 0.48; p = 0.016). The modulation of central beta activity is associated, among other things, with motor imagery, which may be necessary for imitation. Impaired imitation has been associated with core features of ASD. Altered modulation of oscillatory activity may be mechanistically involved in those aspects of motor network function that relate to the core symptoms of ASD.
Introduction
A key hypothesis of the neurobiology of autism spectrum disorders (ASD) is that cortical excitatory function is not sufficiently balanced by inhibitory forces (Rubenstein and Merzenich, 2003; Rubenstein, 2010). This proposition is supported by converging evidence from mouse model histology (Rubenstein, 2010) as well as the clinical observation of an increased prevalence of epilepsy in ASD (Viscidi et al., 2013). As GABAergic activity in the brain is partly responsible for coordination of neural activity, inhibitory alterations in ASD may lead to altered synchronization, which is postulated by some to be a fundamental cause of the ASD phenotype (Uhlhaas and Singer, 2006), possibly by reducing the “signal-to-noise ratio” of task-relevant information (Markram and Markram, 2010).
Oscillatory electrical activity, as captured on the electroencephalogram (EEG), is a phenomenon molded by inhibitory activity. Specifically, the oscillatory activity seen on the EEG reflects, at least in part, the synchronization of neuronal activity through the action of GABAergic interneurons (Uhlhaas and Singer, 2010; Gaetz et al., 2011), whose firing patterns create rhythm in principal cells (Buzsáki and Chrobak, 1995). Alterations of electrical brain oscillatory activity may therefore be a useful window into the pathophysiology of the inhibitory interneuron system in conditions such as ASD (Wilson et al., 2007; Uhlhaas et al., 2010).
Indeed, there is evidence—though limited—for altered cortical oscillatory behavior in ASD. Much of the evidence to date comes from resting state (spontaneous, task-free) EEG and consists of somewhat disparate findings (Cantor et al., 1986; Dawson et al., 1995; Daoust et al., 2004; Chan and Leung, 2006; Murias et al., 2007; Orekhova et al., 2007; Cornew et al., 2012). The oscillatory activity seen on the EEG is not merely an epiphenominological marker of interneuron function, but is rather directly involved in cortical computation (Başar et al., 2001; Kelly et al., 2006; Fries et al., 2007; Engel and Fries, 2010). We therefore intended to examine the role of potential oscillatory changes in ASD in altered cortical processing. This goal is best achieved when oscillatory measures are directly associated with a disorder-relevant task.
Complex motor control is a useful model for studying cortical network dynamics and behavioral performance in ASD. Motor deficits were noted in Kanner’s (1943) original description of autism, and their study is important not only for their own sake, but also because motor control deficits relate theoretically (Mostofsky and Ewen, 2011) and empirically (Dziuk et al., 2007) to the social and communicative impairments in ASD. The benefits of studying the motor system are practical. Neural substrates of the motor system are relatively well characterized (Uttal, 2012), and tasks involving motor output lend themselves to experimental manipulation and reliable quantification. Our laboratory focuses on praxis, the performance of complex, skilled gestures used in functional skills and communication (Heilman and Valenstein, 2003). Neuropsychological and neurophysiological studies of acquired apraxia and developmental dyspraxia have allowed investigators to interrogate the praxis network, which includes inferior parietal and premotor (frontal) areas (Wheaton and Hallett, 2007). Convergent lesion and physiological studies have shown lateralization to the left (Moll et al., 2000; Wheaton and Hallett, 2007).
The link between praxis and the cardinal features of ASD is both theoretical and empirical. On the theoretical side, current frameworks suggest that the motor system plays an important role in the internal simulation of others’ actions, which allows for understanding of others’ intentions and affective states (Klin et al., 2003; Gallese, 2007; Mostofsky and Ewen, 2011). Empirically, we have shown a relationship between deficits in motor control on the one hand and social-communicative impairment in children with ASD on the other (Dziuk et al., 2007; Dowell et al., 2009).
In order to assess oscillatory modulation during a task relevant to the ASD phenotype, we recorded EEG while subjects performed a praxis task. Specifically, we looked at task-related modulation, relative to baseline, of power in the EEG signal (event-related spectral perturbations; ERSP) and hypothesized that children with ASD would show altered ERSP responses in both posterior and central sensors on the left, reflecting the activity of left parietal and frontal premotor praxis-related brain regions. To assess relevance to the ASD phenotype, we examined correlations between physiological EEG measurements and behavioral testing of ASD severity and praxis function.
Materials and Methods
Subjects
Analyses included data from 25 subjects with high-functioning ASD (HFA; defined as ASD with at least low-average IQ; (Ghaziuddin and Mountain-Kimchi, 2004) and 33 controls (or typically developing participants; TD). As described below, 40 additional subjects’ data were removed from analysis due to poor behavioral performance and/or excessive artifact, as set out in section “Signal Pre-Processing”. The bar for data quality was set high to reduce the likelihood of type I errors. Inclusion criteria relevant to both groups included age between the eighth and thirteenth birthday, a Wechsler Intelligence Scale for Children—IV (Wechesler, 2003) full-scale IQ > 80 (except in the presence of a split of greater than 12 points between Verbal Comprehension Index and Perceptual Reasoning Index, in which case the higher score was >80 and the lower >65) and right handedness [Edinburgh (Oldfield, 1971) and PANESS (Denckla, 1985)]. Edinburgh results were to be between 0.5 and 1. If the Edinburgh was not available, children were asked to do 11 gestures from the PANESS; at least 10 must be have been done spontaneously with the right hand for the subject to qualify for the study. General exclusion criteria were diagnosed neurological or chronic medical condition, severe visual impairment (corrected worse than 20/40) and pregnancy.
ASD diagnosis was conservative and based on DSM-IV (American Psychiatric Association, 2000) criteria, with positive scores on both ADOS (Lord et al., 2000, 2012) and ADI-R (Lord et al., 1994), and clinical confirmation by child neurologist experienced in autism diagnosis (SHM). Children with a known neurological, genetic, metabolic or other etiology associated with the ASD were excluded. All children were assessed with the Diagnostic Interview for Children and Adolescents-IV (DICA-IV; Reich et al., 1997) or Schedule for Affective Disorders and Schizophrenia for School-Age Children (K-SADS; Kaufman et al., 2013); children were excluded from the ASD group if they had any diagnosis other than: (1) ASD; (2) an anxiety disorder (given the extensive comorbidity with ASD as well as the large overlap between obsessive-compulsive symptomatology and ASD symptomatology in particular); and (3) attention-deficit/hyperactivity disorder (ADHD; given the high levels of comorbidity with ASD).
Seventeen participants with HFA also had a diagnosis of ADHD. Co-occurring ADHD was defined by convergence of performance on three measures: (1) meeting diagnostic criteria on either the DICA-IV or K-SADS; (2) obtaining a t-score ≥65 on the Connors Parent Rating Scale (Conners et al., 1998) for either the hyperactive/impulsive scale, inattentive scale, or both; and (3) having at least 6 symptoms rated as 2 (often) or 3 (very often) on the DuPaul rating scale (fourth edition; Dupaul et al., 1998) for either the hyperactive/impulsive scale, the inattentive scale, or both. In cases in which one or two measures, but not all three, were positive, the final diagnosis of ADHD was based on DICA/K-SADS and clinical impression from an experienced autism/ADHD clinician and researcher (SHM). Eight subjects with ASD were diagnosed with an anxiety disorder (including obsessive-compulsive disorder, generalized anxiety disorder and simple phobia). Six had no secondary psychiatric diagnosis. Eight participants in the HFA group were on medication (2 dextroamphetamine salts, 5 methylphenadate/dexmethylphenadate, 1 “homeopathic lithium carbonate” + clonidine); stimulant medication was stopped 24 h prior to testing.
Entry into the control group was only permitted for children who did not meet published criteria for ASD on the Social Responsiveness Scale (Constantino and Gruber, 2005) and absence of first-degree family members with ASD. Children were also ineligible for entry into the TD group if they met diagnostic criteria for any disorder on the DICA-IV (Reich et al., 1997), except for simple phobia. No subject in the control group had a psychiatric diagnosis. Seven control subjects were on second-generation histamine blockers for seasonal allergies, inhaled steroids and leukotriene inhibitors and as-needed decongestants.
The protocol was approved by Johns Hopkins Medicine Institutional Review Board. Written consent was obtained from parents and, when possible, from subjects. Otherwise, verbal assent was obtained from subjects.
Task and Behavioral Monitoring
The task was based on that used by Wheaton et al. (2009), as described in a previous study (Ewen et al., 2014). Modifications to Wheaton’s original paradigm were designed to mitigate confounds due to known working memory and linguistic deficits in ASD. The paradigm consisted of pantomiming the use of 10 common tools (scissors, spoon [to stir hot chocolate], ice cream scoop, doorknob, pencil, screwdriver, hammer, paintbrush, key, chalkboard eraser). Prior to the recording, the participants demonstrated the correct use of each of the tools to the experimenter, who remained next to the subject for the duration of the task.
The stimuli were presented using eevoke software (Advanced Neuro Technologies [ANT], Netherlands). During the pre-stimulus portion of each trial, subjects fixated on a cross at the center of the computer monitor for 4 s (Figure 1). Next, the Prepare stimulus appeared: this was a photograph of one of the ten tools; participants were instructed not to make any movements during this time. After the Prepare stimulus remained on the screen for 3 s, a green box appeared around the photograph of the tool; this was the “Go” stimulus. The Go stimulus remained on the screen for 3.5 s, during which time the subjects were to pantomime the use of the tool with their right hand until the word “Rest” appeared. “Rest” lasted 2 s and was replaced by the fixation cross for the next trial. Six blocks of 20 trials each (120 trials total) were recorded in each subject.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/187244/fnhum-10-00198-HTML/image_m/fnhum-10-00198-g001.jpg)
Figure 1. Praxis task. The tool photograph that appears immediately after the fixation cross is the Prepare stimulus, during the presentation of which subjects do not move. When the green frame appears around the photograph (Go stimulus), subjects pantomime the use of the tool.
To ensure adequate behavioral performance, a psychological associate, trained in the conduct of psychological and motor research in children with developmental disabilities and determined to be research-reliable in the offline coding of the pediatric adaptation of the Florida Apraxia Battery (FAB; Mostofsky et al., 2006), monitored the performance of the subject and used a button to mark the electronic record whenever a subject started the movement early (before the onset of Go) or late (>1000 ms following the onset of Go), finished early or made errors in the movement/no response. Movement errors were coded by the same stringent criteria as in Mostofsky et al. (2006) and other publications examining dyspraxia in ASD (Dziuk et al., 2007; Dowell et al., 2009; MacNeil and Mostofsky, 2012). All sessions were video recorded, and another experimenter examined the video of each trial by the same criteria. Only movements that were initiated within the allowable time-frame and were error-free by FAB criteria, by both raters, were included in the EEG analysis.
Recordings
Using an Advanced Neuro Technologies asa-lab system (Netherlands) and elastic mesh Waveguard cap (Figure 2), EEG data were recorded from 47 sites covering the whole scalp with approximately uniform density. The recording was performed in DC mode, with a 512 Hz sampling rate and an anti-aliasing filter with a 138 Hz cut off. Each channel was referenced to an average of all channels during recording. The electrode cap used active cable-shielding technology. Electrode impedance was kept below 5 kΩ in all channels.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/187244/fnhum-10-00198-HTML/image_m/fnhum-10-00198-g002.jpg)
Figure 2. Electrode layout. The 47 electrodes used in recording are represented by open circles. Electrodes within the left central regions of interest (ROI) are indicated in green, and electrodes within the left posterior ROI are indicated in red.
Signal Pre-Processing
To minimize the spatial blurring due to volume conduction, signals were converted from average reference to current source density (CSD) estimates, computed using the spherical spline algorithm (Kayser and Tenke, 2006). In the CSD montage, vertical electro-oculograms (VEOG) were recorded from frontal channels whose locations were designed specifically to capture eye blinks. Eye blink correction was then performed using a principal component analysis (PCA) method that models the brain signal and artifact subspaces (Ille et al., 2002).
CSD data from each trial were epoched by time-locking to the onset of the Prepare stimulus. Epochs were created with a pre-stimulus (baseline) length of 2 s and post-Prepare-stimulus-onset (trial) length of 6 s, for a total length of 8 s. In the remainder of the article, event-related timing is referred to in reference to the onset of the Prepare stimulus.
Epochs with horizontal eye movement artifacts, recognized as step-like potentials of opposite polarity in channels RE1 and LE1 (electrodes located 1–2 cm lateral to external epicanthi), were removed if the amplitude exceeded a threshold defined individually for each subject (based on analysis of definite lateral eye movements). Muscular and movement artifacts (such as jaw clenches) were identified by visually inspecting the signal and also by video verification. Epochs with these artifacts were removed from analysis.
We excluded from all analyses all subjects who did not finish all six runs of the experiment. We further excluded any subject for whom more than 50 percent of trials had to be excluded due to errors in behavioral performance, as described in section “Recordings”. In total, 18 subjects were removed from the HFA group and 14 from the TD group for insufficient behavioral performance. Following pre-processing of the EEG, we next excluded any subject who failed to have more than 50 percent of trials that were both artifact-free and behaviorally correct. Four subjects were removed from each group for insufficient trials after EEG pre-processing. All excluded subjects were removed prior to ERSP analysis.
Time-Frequency Decomposition
Time-frequency decomposition of the CSD data was accomplished using Morlet wavelet decomposition in EEGLAB (Delorme and Makeig, 2004), with a window size of 512 points and frequencies of 0.5–30 Hz. Each time-frequency point of the active portion of the ERSP analysis (i.e., following the onset of the Prepare stimulus) was calculated as a z-score, based on the baseline period.
For each trial, at each frequency, all measurements of signal energy during the 2-s baseline interval were aggregated into a distribution. Next, the ERSP value at each time-frequency point during the active phase (0–6 s, relative to the onset of Prepare) was calculated as a z-score value relative to the baseline distribution for the same trial, at the relevant frequency. We did not remove phase-locked activity prior to ERSP calculation.
Channel Selection
We next selected relevant channels based on the known anatomy of the praxis network and on precedent from previous EEG studies of praxis (Wheaton et al., 2005b; Ewen et al., 2014). We selected electrodes in scalp regions superficial to the central (premotor/primary motor) and parietal areas that constitute the principal nodes of the praxis network. Although the article from Hallett’s laboratory (Wheaton et al., 2005b) used a different electrode layout, and the present study used a more parsed sensor array compared with our prior study (Ewen et al., 2014), the electrodes selected in the current study represent sensor locations that overlap with those in the previous works. As described below (see “ERD Quantification” Section), we used regions of interest (ROI) to characterize the amount of ERSP in each individual’s central and parietal areas. The rationale for this approach was that using multiple channels in a ROI, rather than a single channel, would mitigate the variance in correspondence between channels and brain generators between patients. Further, the effects of volume conduction, though already limited through the use of CSD spatial filtering, would be further decreased by averaging out activity peripheral to the ROI.
Frequency Selection
Using the grand average spectrogram for each ROI, for each group (i.e., four in total), we defined the frequency ranges of interest (Figure 3) by examining delimited regions in the frequency axis of ERSP. There was great consistency in relevant frequency bands between central and posterior spectrograms, and these frequency ranges were consistent with traditionally described frequency bands: theta (2–8 Hz), alpha (7–13 Hz) and beta (18–22 Hz). Of note, the beta band selected comports with that used by Wheaton et al. (2005a, 2009) in a series of EEG articles examining praxis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/187244/fnhum-10-00198-HTML/image_m/fnhum-10-00198-g003.jpg)
Figure 3. Grand average spectrograms from children with high-functioning ASD (HFA; A,B) and typically developing (TD) controls (C,D). Cool colors represent event-related desynchronization (ERD), and warm colors represent event-related synchronization (ERS). ERD is indicated in unitless z-scores. Time is indicated along the x-axis (in seconds), and frequency is indicated along the y-axis (in Hz). The top spectrograms (A,C) represents the left central ROI, and the bottom, the left parietal ROI (B,D). Dotted boxes indicate time-frequency regions that were averaged and compared between groups; regions that were statistically significantly different between groups are in white boxes; those which were not different are in black boxes.
ERD Quantification
For each subject, we quantified mean event-related desynchronization (ERD) magnitude separately in eight time-frequency-topography windows: two time ranges (Prepare = 0–2 s, and Go = 3–5 s) × the two frequency ranges (alpha and beta) × two ROIs (left central and left posterior ROI, averaging over the three channels in each ROI; Figure 2). For alpha and beta separately, we performed a Group (ASD, TD) × Task (Prepare, Go) × ROI (LP, LC) repeated measures ANOVA, with Task and ROI as the within-subjects factors and Group as the between-subjects factor, using SPSS. Post hoc, we examined for group differences by using Welch’s t-test (no assumption of equal variances; two unmatched samples, two-sided) in each of the eight time-frequency-topography windows.
We additionally measured mean baseline EEG power in each ROI, in each frequency band (alpha and beta), using the time window of −2 to 0 s relative to the onset of the Prepare stimulus. We compared the samples from each group using Welch’s t-test (no assumption of equal variances; two unmatched samples, two-sided) as well as JZS Bayesian Factor and Information Bayesian Factor.
Additional Behavioral Measures and Comparisons with Physiological Measures
A priori, we defined a set of clinical and behavioral tests we believed could have a theoretical relationship to the physiological effects being measured. These included the ADOS, SRS and a pediatric modification of the FAB (Gonzalez Rothi et al., 1997; Mostofsky et al., 2006), which is a more comprehensive praxis competence assessment than the task used during EEG and better able to discriminate levels of ability amongst subjects. The modified FAB examines a wide range of praxis-relevant motor behavior, including tool use gestures, communicative gestures and meaningless gestures, under the conditions of gesture-to-command, gesture-to-imitation and performance of tool use gestures while holding the actual tool.
Subjects’ results on the praxis battery were compared between groups (using Student’s t-test). We then assessed correlations, separately within the two groups, between those ERSP measures found to be different between groups on the one hand and praxis, ADOS and SRS measures on the other. Correlations were assessed using Pearson r-tests. Some subjects were tested with the ADOS-G and others with the ADOS-2, so we adjusted ADOS-2 severity scores to be comparable to ADOS-G scores (Gotham et al., 2009).
Results
Participants
Average age for children with HFA was 10.7 ± 1.4 years; for control children, 10.52 ± 1.3 years; two-sided Student’s t-test p = 0.62). Twenty four percent of HFA subjects were female; 24% of TD subjects were female (Fisher’s exact test p = 1). Full-scale IQ was not different between groups (mean FSIQ for TD = 115.6, and for HFA = 109.3; p = 0.095). To further examine for group differences in IQ, we used two Bayesian analyses, each with two priors. The Scaled JZS Bayes Factor = 1.15, and Scaled-Information Bayes Factor = 0.825 (Rouder et al., 2009), both failing to provide evidence of group differences in IQ. Within the HFA group, the mean ADOS score was 12.5 (range: 8–20).
Event-Related Spectral Perturbations
Subjects from both groups showed a characteristic pattern of ERSP (Figures 3, 4). Specifically, brief event-related synchronization (ERS; task-related increase in power) was seen in the delta/theta band (2–8 Hz), while sustained event-related desynchronization (ERD; task-related decrease in power) was seen in both alpha (7–13 Hz) and beta (18–22 Hz). We examined the magnitude of alpha and beta ERD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/187244/fnhum-10-00198-HTML/image_m/fnhum-10-00198-g004.jpg)
Figure 4. Topographical distribution of alpha and beta ERD in Prepare (A) and Go (B). Alpha ERD (in unitless z-scores) was seen in occipital and parietal scalp regions following Prepare, and beta ERD was seen principally in parietal regions during this phase. Following Go, both alpha and beta ERD were seen in parietal and central (motor/premotor) regions. Signals were analyzed from left central and left posterior ROIs. Grand average topographical plots from children with HFA on the left and typically developing (TD) children on the right.
In the alpha band, the RM-ANOVA showed a significant effect of task (F(1,56) = 89.686, p < 0.0001, = 0.616) and ROI (F(1,56) = 5.137, p = 0.027, = 0.084). There were also significant interactions between ROI and Group (F(1,56) = 11.466, p < 0.017, = 0.097), and between ROI and Task (F(1,56) = 6.221, p = 0.016, = 0.1). In the beta band, we found a significant effect of task (F(1,56) = 74.797, p < 0.0001, = 0.572) and a significant interaction between ROI and task (F(1,56) = 5.598, p = 0.021, = 0.91). The full ANOVA tables are presented in the Supplementary Material (see “Supplementary Data 2”).
Post hoc t-tests showed children with HFA had, on average, a decreased magnitude of alpha ERD in the posterior ROI (29% difference in sample means, p = 0.0045) during Prepare. During Prepare, beta ERD showed near-significant differences (25% difference in sample means; p = 0.06).
During the Go phase, a decreased magnitude of left parietal alpha ERD was again seen in HFA (24%; p = 0.030), and group differences in left central beta ERD were statistically significant (27%; p = 0.046).
The ERD showed similar time courses in both groups (Figure 5), suggesting that group differences in timing of praxis initiation or temporal dispersion (group differences in variability of timing across trials) are unlikely to explain measured group differences in mean ERD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/187244/fnhum-10-00198-HTML/image_m/fnhum-10-00198-g005.jpg)
Figure 5. ERD-over-time, in topographies and bands found to be different between groups. The x-axis reflects time (in seconds), and the y-axis represents ERD magnitude (in unitless z-scores). The grand average for the TD group is plotted in blue, and the grand average for the HFA group is plotted in red. Shaded areas around the mean (heavy line) represent standard errors. The Prepare phase ERD measurement in statistical calculations is the average from 0 to 2 s on the x-axis, and the Go phase, from 3 to 5 s. Plot A represents left posterior alpha ERD. Group differences were seen during both the Prepare and Go phases (gray shaded areas). In left central beta ERD (Plot B), the Go phase ERD was significantly different between groups (gray shaded area). The plots suggest that a delay in onset of ERD or temporal dispersion is unlikely to explain group differences in average ERD.
Because of significant group differences in the number of valid trials per subject (HFA = 86.12 ± 14.9; TD = 102.5 ± 11.2; p < 0.001), the above t-tests were conducted without the assumption of equal variances, as described in section “ERD Quantification”.
There were no statistical differences between groups in the baseline spectral power, in either ROI × either frequency band (see “Supplementary Data 3”).
Behavioral Measures and Correlations with EEG
Children with HFA had significantly worse scores on the praxis performance battery (mean total percent correct = 74% for TD and 53% for HFA, p < 0.001), consistent with several previously published studies that included samples that did not overlap with that from the current study (Mostofsky et al., 2006; Dziuk et al., 2007; Dowell et al., 2009; MacNeil and Mostofsky, 2012). SRS data also, expectedly, differed between groups (HFA mean: 97; TD mean: 17.6; p < 0.001).
We correlated behavioral results (ADOS, SRS and pediatric FAB) with those ERSP measures that were different between groups. No statistically significant correlations were seen in the TD group. Within the HFA group, left central ROI beta ERD (during Go) correlated with praxis imitation performance (r = −0.4; p = 0.04). This left central beta ERD also correlated with ADOS total score (r = 0.48; p = 0.016). (See “Supplementary Data 4”). Post hoc analysis of ADOS sub-scores showed correlations between the same ERD measurement on one hand, and the Communications sub-score and Repetitive and Restricted Behaviors subscores of the ADOS on the other.
Discussion
The primary finding of this research is decreased magnitude of task-related EEG power modulation in HFA during the performance of skilled movements on a praxis task—which is known to be impaired in ASD and has been shown to be associated with the core features of the disorder (Dziuk et al., 2007). These findings suggest that autism and its associated dyspraxia may be associated with decreased activity in the frontal-parietal praxis network.
Alpha activity was originally described as representing an idling state of the underlying cortical region (Pfurtscheller et al., 1996), however more recent accounts have demonstrated the role of alpha activity in inhibiting sensory input and/or maintenance of sensory-linked representations in primary sensory cortex (Kelly et al., 2006; Tuladhar et al., 2007; Ikkai et al., 2014). Suppression (ERD) of alpha activity therefore correlates with increased task-related activity of a cerebral region and differentiated the two groups. Decreased modulation of alpha in the parietal cortex suggests decreased task-related parietal activation in HFA.
Beta is a critical frequency in the motor system, having been documented not only in motor cortex, but in basal ganglia and muscle activity, as well as coherently among the parts of the motor system (Baker et al., 1999). Beta activity has been linked with the resting state of the motor system (Neuper et al., 2006) or with the maintenance of the current state of the motor system (Engel and Fries, 2010), while beta ERD has been linked with top-down motor control and with motor imagery (Neuper et al., 2006; Zhang et al., 2008; Engel and Fries, 2010). The current study showed that decreased central beta ERD during the praxis task correlates within the HFA group with praxis imitation ability and with severity of autism symptoms. Imitation has long been known to be impaired in ASD (DeMeyer et al., 1972; Williams et al., 2004), and it may be that deficiencies in the ability of the imagery and production system to activate in complex motor tasks underlies both impaired imitation and impaired complex gesture production.
In this experiment, we have demonstrated correlations between behavior and physiology. One challenge endemic in this type of research is determining which level of analysis drives the other. On one hand, it may be cognitive differences—in attention, motor imagery or event task performance—that drive physiological markers of task participation—ERSP. On the other hand, it may be altered neurobiological mechanisms—the ability to modulate cortical rhythms to respond to task demands—that produce behavioral group differences. We took pains to limit the influence of behavioral group differences on the EEG measures, by creating an easy task and by removing all behaviorally bad trials from EEG analysis, but it is impossible to rule out that there could be performance-related differences that were not quantified yet influenced the EEG responses. Unlike in research on lower-dimensional motor behaviors, it is difficult to find one or a few kinematic variables that adequately describe the movement, so as to control completely for behavioral differences.
The alternative view is that biological differences in the ability to modulate physiological rhythms in the brain are causally responsible for group differences in performance on the praxis task and other measures of motor and social/communicative differences. While the current data are not conclusive, converging data may bolster the latter view. First, a collection of results demonstrating an association between specific deficits in ASD and well-established behavioral effects of oscillatory modulation would provide evidence supportive of the proposition that the physiological effects are primary. The literature is in early stages, but a recent study of multisensory attention selection, an impaired ability in ASD, has demonstrated deficient alpha ERS (Murphy et al., 2014); given the putative role of alpha in sensory suppression (Kelly et al., 2006; Tuladhar et al., 2007), there is a compelling argument that deficits in alpha modulation cause alterations of multisensory attentional behavior. Further, future work with animal models, cellular physiology and computer modeling of pathological differences may establish the presence of physiological differences build up from knowledge of genetic mechanisms, without resorting to behavioral tasks. Finally, and perhaps most compellingly, interventional approaches with non-invasive brain stimulation may be able to perturb oscillatory activity sufficiently to draw direct inferences about the role of oscillations in behavioral differences in ASD (Helfrich et al., 2014). Each of these approaches has its limitations, but together, they would shed the most conclusive light on the purported role of altered oscillatory activity and neural synchrony in ASD (Uhlhaas and Singer, 2006; Uhlhaas et al., 2009), due to altered excitatory/inhibitory balance (Rubenstein and Merzenich, 2003) and the functional consequences of altered cytoarchitecture (Casanova, 2007).
Because this study is a novel examination of multidimensional physiological data (i.e., modulation of oscillations within the praxis network, in ASD), the analysis is necessarily somewhat exploratory. A further limitation is that we only examined children with high functioning autism. Lower-functioning children are worthy of study, and, indeed electrophysiological techniques hold the promise of providing mechanistic insights and biomarkers in children who are cannot be tested behaviorally because they are unable to make complex behavioral responses. In order to develop these biomarkers, however, we need to validate the link between EEG measurement and behavior in higher functioning children who are indeed able to make behavioral responses so as to validate the relevance of the physiological markers.
A few notes should be made about the extensive efforts made to reduce the likelihood of spurious group differences in the EEG data (Webb et al., 2015). First, the task contained both a motor execution phase (to monitor compliance with the task) and a motor planning phase (which would not be affected by differences in motion artifact). That group differences were seen during the Prepare phase limits the ability to ascribe group differences to motion artifact. Further, experimental effects were specific in topography and frequency, whereas many types of EEG artifact are symmetrical and broadband. Absence of group differences in the baseline spectra also suggests that artifact does not play a role in the ERD differences. Absence of baseline spectral differences also minimizes the likelihood that differences in overall cognitive state (e.g., attention) impacted ERD results. Further, as described above, we minimized the potential impact of performance-related group differences by creating a relatively easy task and removing trials (and subjects) with excessive errors. Since children with ASD may have differences in the timing of cortical events, we constructed both the task and analysis such that they are relatively insensitive to changes in timing of a relevant magnitude; additionally, we removed trials where responses were delayed or anticipatory. As described above, we constructed the stimuli to minimize the effects of potential group differences in reading and working memory ability. Finally, we controlled for potentially spurious alterations of SNR by equating the number of trials in the two groups.
In summary, children with HFA show a decreased ability to generate task-related changes in ongoing oscillatory activity during a praxis task. These group differences occurred in specific cerebral regions and frequency bands that are known to be relevant to the network involved. Electrophysiological differences correlate with autism symptom severity and with diagnosis-relevant differences in praxis function. These findings suggest that alterations in the mechanisms that generate cortical oscillations may contribute to known differences in praxis function and in the core aspects of autism.
Author Contributions
JBE, BML, MH, NEC, and SHM contributed to the design of the experiments. ASP, DM, BML, and CN contributed to the acquisition of data. JBE, ASP, DMcA, and BML contributed to the analysis of data. JBE, ASP, DMcA, BML, MH, NEC, and SHM contributed to the interpretation of the data. All authors reviewed the manuscript for important intellectual content.
Funding
This work was funded by the National Institute of Neurological Disorders and Stroke/National Institutes of Health (Intramural Program—MH; K23NS073626 to JBE; and R01NS048527 to SHM) and Autism Speaks (to SHM).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We are grateful to Nancy Kopell, Michael Johnston, Gary Goldstein and Nitish Thakor for their helpful suggestions on a previous version of this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fnhum.2016.00198/abstract
References
American Psychiatric Association. (2000). Diagnostic and Statistical Manual of Mental Disorders: DSM-IV-TR. Washington, DC: American Psychiatric Association.
Baker, S. N., Kilner, J. M., Pinches, E. M., and Lemon, R. N. (1999). The role of synchrony and oscillations in the motor output. Exp. Brain Res. 128, 109–117. doi: 10.1007/s002210050825
Başar, E., Başar-Eroglu, C., Karakaş, S., and Schürmann, M. (2001). Gamma, alpha, delta and theta oscillations govern cognitive processes. Int. J. Psychophysiol. 39, 241–248. doi: 10.1016/s0167-8760(00)00145-8
Buzsáki, G., and Chrobak, J. J. (1995). Temporal structure in spatially organized neuronal ensembles: a role for interneuronal networks. Curr. Opin. Neurobiol. 5, 504–510. doi: 10.1016/0959-4388(95)80012-3
Cantor, D. S., Thatcher, R. W., Hrybyk, M., and Kaye, H. (1986). Computerized EEG analyses of autistic children. J. Autism Dev. Disord. 16, 169–187. doi: 10.1007/bf01531728
Casanova, M. F. (2007). The neuropathology of autism. Brain Pathol. 17, 422–433. doi: 10.1111/j.1750-3639.2007.00100.x
Chan, A. S., and Leung, W. W. (2006). Differentiating autistic children with quantitative encephalography: a 3-month longitudinal study. J. Child Neurol. 21, 391–399. doi: 10.2310/7010.2006.00094
Conners, C. K., Sitarenios, G., Parker, J. D., and Epstein, J. N. (1998). The revised conners’ parent rating scale (CPRS-R): factor structure, reliability, and criterion validity. J. Abnorm. Child Psychol. 26, 257–268. doi: 10.1023/A:1022602400621
Constantino, J., and Gruber, C. (2005). Social Responsiveness Scale. Los Angeles, CA: Western Psychological Services.
Cornew, L., Roberts, T. P., Blaskey, L., and Edgar, J. C. (2012). Resting-state oscillatory activity in autism spectrum disorders. J. Autism Dev. Disord. 42, 1884–1894. doi: 10.1007/s10803-011-1431-6
Daoust, A. M., Limoges, E., Bolduc, C., Mottron, L., and Godbout, R. (2004). EEG spectral analysis of wakefulness and REM sleep in high functioning autistic spectrum disorders. Clin. Neurophysiol. 115, 1368–1373. doi: 10.1016/j.clinph.2004.01.011
Dawson, G., Klinger, L. G., Panagiotides, H., Lewy, A., and Castelloe, P. (1995). Subgroups of autistic children based on social behavior display distinct patterns of brain activity. J. Abnorm. Child Psychol. 23, 569–583. doi: 10.1007/bf01447662
Delorme, A., and Makeig, S. (2004). EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J. Neurosci. Methods 134, 9–21. doi: 10.1016/j.jneumeth.2003.10.009
DeMeyer, M. K., Alpern, G. D., Barton, S., Demyer, W. E., Churchill, D. W., Hingtgen, J. N., et al. (1972). Imitation in autistic, early schizophrenic and non-psychotic subnormal children. J. Autism Child Schizophr. 2, 264–287. doi: 10.1007/bf01537618
Denckla, M. B. (1985). Revised neurological examination for subtle signs. Psychopharmacol. Bull. 21, 773–779.
Dowell, L. R., Mahone, E. M., and Mostofsky, S. H. (2009). Associations of postural knowledge and basic motor skill with dyspraxia in autism: implication for abnormalities in distributed connectivity and motor learning. Neuropsychology 23, 563–570. doi: 10.1037/a0015640
Dupaul, G., Power, T., Anastopoulos, A., and Reid, R. (1998). ADHD Rating Scale-IV. New York, NY: Guilford Press.
Dziuk, M. A., Gidley Larson, J. C., Apostu, A., Mahone, E. M., Denckla, M. B., and Mostofsky, S. H. (2007). Dyspraxia in autism: association with motor, social and communicative deficits. Dev. Med. Child Neurol. 49, 734–739. doi: 10.1111/j.1469-8749.2007.00734.x
Engel, A. K., and Fries, P. (2010). Beta-band oscillations—signalling the status quo? Curr. Opin. Neurobiol. 20, 156–165. doi: 10.1016/j.conb.2010.02.015
Ewen, J. B., Lakshmanan, B. M., Hallett, M., Mostofsky, S. H., Crone, N. E., and Korzeniewska, A. (2014). Dynamics of functional and effective connectivity within human cortical motor control networks. Clin. Neurophysiol. 126, 987–996. doi: 10.1016/j.clinph.2014.09.006
Fries, P., Nikolić, D., and Singer, W. (2007). The gamma cycle. Trends Neurosci. 30, 309–316. doi: 10.1016/j.tins.2007.05.005
Gaetz, W., Edgar, J. C., Wang, D. J., and Roberts, T. P. (2011). Relating MEG measured motor cortical oscillations to resting gamma-aminobutyric acid (GABA) concentration. Neuroimage 55, 616–621. doi: 10.1016/j.neuroimage.2010.12.077
Gallese, V. (2007). Before and below ‘theory of mind’: embodied simulation and the neural correlates of social cognition. Philos. Trans. R. Soc. Lond. B Biol. Sci. 362, 659–669. doi: 10.1098/rstb.2006.2002
Ghaziuddin, M., and Mountain-Kimchi, K. (2004). Defining the intellectual profile of Asperger syndrome: comparison with high-functioning autism. J. Autism Dev. Disord. 34, 279–284. doi: 10.1023/b:jadd.0000029550.19098.77
Gonzalez Rothi, L., Raymer, A., and Heilman, K. (1997). “Limb praxis assessment,” in Apraxia: The Cognitive Neuropsychology of Praxis, eds L. J. G. Rothi and K. M. Heilman (Hove, UK: Psychology Press), 61–73.
Gotham, K., Pickles, A., and Lord, C. (2009). Standardizing ADOS scores for a measure of severity in autism spectrum disorders. J. Autism Dev. Disord. 39, 693–705. doi: 10.1007/s10803-008-0674-3
Heilman, K. M., and Valenstein, E. (2003). Clinical Neuropsychology. Oxford, NY: Oxford University Press.
Helfrich, R. F., Schneider, T. R., Rach, S., Trautmann-Lengsfeld, S. A., Engel, A. K., and Herrmann, C. S. (2014). Entrainment of brain oscillations by transcranial alternating current stimulation. Curr. Biol. 24, 333–339. doi: 10.1016/j.cub.2013.12.041
Ikkai, A., Blacker, K. J., Lakshmanan, B. M., Ewen, J. B., and Courtney, S. M. (2014). Maintenance of relational information in working memory leads to suppression of the sensory cortex. J. Neurophysiol. 112, 1903–1915. doi: 10.1152/jn.00134.2014
Ille, N., Berg, P., and Scherg, M. (2002). Artifact correction of the ongoing EEG using spatial filters based on artifact and brain signal topographies. J. Clin. Neurophysiol. 19, 113–124. doi: 10.1097/00004691-200203000-00002
Kaufman, J., Birmaher, B., Axelson, D., Perepletchikova, F., Brent, D., and Ryan, N. (2013). Kiddie Schedule for Affective Disorders and Schizophrenia for School-Aged Children–Lifetime Version (Kiddie-SADS-PL 2013 Working Draft). Pittsburgh, PA: Western Psychiatric Institute and Clinic and Yale University.
Kayser, J., and Tenke, C. E. (2006). Principal components analysis of Laplacian waveforms as a generic method for identifying ERP generator patterns: I. Evaluation with auditory oddball tasks. Clin. Neurophysiol. 117, 348–368. doi: 10.1016/j.clinph.2005.08.034
Kelly, S. P., Lalor, E. C., Reilly, R. B., and Foxe, J. J. (2006). Increases in alpha oscillatory power reflect an active retinotopic mechanism for distracter suppression during sustained visuospatial attention. J. Neurophysiol. 95, 3844–3851. doi: 10.1152/jn.01234.2005
Klin, A., Jones, W., Schultz, R., and Volkmar, F. (2003). The enactive mind, or from actions to cognition: lessons from autism. Philos. Trans. R. Soc. Lond. B Biol. Sci. 358, 345–360. doi: 10.1098/rstb.2002.1202
Lord, C., Risi, S., Lambrecht, L., Cook, E. H., Jr., Leventhal, B. L., DiLavore, P. C., et al. (2000). The autism diagnostic observation schedule-generic: a standard measure of social and communication deficits associated with the spectrum of autism. J. Autism Dev. Disord. 30, 205–223. doi: 10.1023/A:1005592401947
Lord, C., Rutter, M., Dilavore, P., Risi, S., Gotham, K., and Bishop, S. (2012). Autism Diagnostic Observation Schedule Second Edition (ADOS-2) Manual (Part 1): Modules 1–4. Torrrance, CA: Western Psychological Services.
Lord, C., Rutter, M., and Le Couteur, A. (1994). Autism diagnostic interview-revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders. J. Autism Dev. Disord. 24, 659–685. doi: 10.1007/bf02172145
MacNeil, L. K., and Mostofsky, S. H. (2012). Specificity of dyspraxia in children with autism. Neuropsychology 26, 165–171. doi: 10.1037/a0026955
Markram, K., and Markram, H. (2010). The intense world theory–a unifying theory of the neurobiology of autism. Front. Hum. Neurosci. 4:224. doi: 10.3389/fnhum.2010.00224
Moll, J., de Oliveira-Souza, R., Passman, L. J., Cunha, F. C., Souza-Lima, F., and Andreiuolo, P. A. (2000). Functional MRI correlates of real and imagined tool-use pantomimes. Neurology 54, 1331–1336. doi: 10.1212/wnl.54.6.1331
Mostofsky, S. H., Dubey, P., Jerath, V. K., Jansiewicz, E. M., Goldberg, M. C., and Denckla, M. B. (2006). Developmental dyspraxia is not limited to imitation in children with autism spectrum disorders. J. Int. Neuropsychol. Soc. 12, 314–326. doi: 10.1017/s1355617706060437
Mostofsky, S. H., and Ewen, J. B. (2011). Altered connectivity and action model formation in autism is autism. Neuroscientist 17, 437–448. doi: 10.1177/1073858410392381
Murias, M., Webb, S. J., Greenson, J., and Dawson, G. (2007). Resting state cortical connectivity reflected in EEG coherence in individuals with autism. Biol. Psychiatry 62, 270–273. doi: 10.1016/j.biopsych.2006.11.012
Murphy, J. W., Foxe, J. J., Peters, J. B., and Molholm, S. (2014). Susceptibility to distraction in autism spectrum disorder: probing the integrity of oscillatory alpha-band suppression mechanisms. Autism Res. 7, 442–458. doi: 10.1002/aur.1374
Neuper, C., Wörtz, M., and Pfurtscheller, G. (2006). ERD/ERS patterns reflecting sensorimotor activation and deactivation. Prog. Brain Res. 159, 211–222. doi: 10.1016/s0079-6123(06)59014-4
Oldfield, R. C. (1971). The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia 9, 97–113. doi: 10.1016/0028-3932(71)90067-4
Orekhova, E. V., Stroganova, T. A., Nygren, G., Tsetlin, M. M., Posikera, I. N., Gillberg, C., et al. (2007). Excess of high frequency electroencephalogram oscillations in boys with autism. Biol. Psychiatry 62, 1022–1029. doi: 10.1016/j.biopsych.2006.12.029
Pfurtscheller, G., Stancák, A., Jr., and Neuper, C. (1996). Event-related synchronization (ERS) in the alpha band–an electrophysiological correlate of cortical idling: a review. Int. J. Psychophysiol. 24, 39–46. doi: 10.1016/s0167-8760(96)00066-9
Reich, W., Welner, Z., and Herjanic, B. (1997). The Diagnostic Interview for Children and Adolescents-IV. North Tonawanda: Multi-Health Systems.
Rouder, J. N., Speckman, P. L., Sun, D., Morey, R. D., and Iverson, G. (2009). Bayesian t-tests for accepting and rejecting the null hypothesis. Psychon. Bull. Rev. 16, 225–237. doi: 10.3758/PBR.16.2.225
Rubenstein, J. L. (2010). Three hypotheses for developmental defects that may underlie some forms of autism spectrum disorder. Curr. Opin. Neurol. 23, 118–123. doi: 10.1097/WCO.0b013e328336eb13
Rubenstein, J. L., and Merzenich, M. M. (2003). Model of autism: increased ratio of excitation/inhibition in key neural systems. Genes Brain Behav. 2, 255–267. doi: 10.1034/j.1601-183x.2003.00037.x
Tuladhar, A. M., ter Huurne, N., Schoffelen, J. M., Maris, E., Oostenveld, R., and Jensen, O. (2007). Parieto-occipital sources account for the increase in alpha activity with working memory load. Hum. Brain Mapp. 28, 785–792. doi: 10.1002/hbm.20306
Uhlhaas, P. J., Pipa, G., Lima, B., Melloni, L., Neuenschwander, S., Nikolic, D., et al. (2009). Neural synchrony in cortical networks: history, concept and current status. Front. Integr. Neurosci. 3:17. doi: 10.3389/neuro.07.017.2009
Uhlhaas, P. J., Roux, F., Rodriguez, E., Rotarska-Jagiela, A., and Singer, W. (2010). Neural synchrony and the development of cortical networks. Trends Cogn. Sci. 14, 72–80. doi: 10.1016/j.tics.2009.12.002
Uhlhaas, P. J., and Singer, W. (2006). Neural synchrony in brain disorders: relevance for cognitive dysfunctions and pathophysiology. Neuron 52, 155–168. doi: 10.1016/j.neuron.2006.09.020
Uhlhaas, P. J., and Singer, W. (2010). Abnormal neural oscillations and synchrony in schizophrenia. Nat. Rev. Neurosci. 11, 100–113. doi: 10.1038/nrn2774
Uttal, W. (2012). Reliability in Cognitive Neuroscience: A Meta-Meta-Analysis. Cambridge, MA: The MIT Press.
Viscidi, E. W., Triche, E. W., Pescosolido, M. F., McLean, R. L., Joseph, R. M., Spence, S. J., et al. (2013). Clinical characteristics of children with autism spectrum disorder and co-occurring epilepsy. PLoS One 8:e67797. doi: 10.1371/journal.pone.0067797
Webb, S. J., Bernier, R., Henderson, H. A., Johnson, M. H., Jones, E. J., Lerner, M. D., et al. (2015). Guidelines and best practices for electrophysiological data collection, analysis and reporting in autism. J. Autism Dev. Disord. 45, 425–443. doi: 10.1007/s10803-013-1916-6
Wechesler, D. (2003). Wechsler Intelligence Scale for Children, Fourth Edition. San Antonio, TX: The Psychological Corporation.
Wheaton, L., Fridman, E., Bohlhalter, S., Vorbach, S., and Hallett, M. (2009). Left parietal activation related to planning, executing and suppressing praxis hand movements. Clin. Neurophysiol. 120, 980–986. doi: 10.1016/j.clinph.2009.02.161
Wheaton, L. A., and Hallett, M. (2007). Ideomotor apraxia: a review. J. Neurol. Sci. 260, 1–10. doi: 10.1016/j.jns.2007.04.014
Wheaton, L. A., Nolte, G., Bohlhalter, S., Fridman, E., and Hallett, M. (2005a). Synchronization of parietal and premotor areas during preparation and execution of praxis hand movements. Clin. Neurophysiol. 116, 1382–1390. doi: 10.1016/j.clinph.2005.01.008
Wheaton, L. A., Shibasaki, H., and Hallett, M. (2005b). Temporal activation pattern of parietal and premotor areas related to praxis movements. Clin. Neurophysiol. 116, 1201–1212. doi: 10.1016/j.clinph.2005.01.001
Williams, J. H., Whiten, A., and Singh, T. (2004). A systematic review of action imitation in autistic spectrum disorder. J. Autism Dev. Disord. 34, 285–299. doi: 10.1023/b:jadd.0000029551.56735.3a
Wilson, T. W., Rojas, D. C., Reite, M. L., Teale, P. D., and Rogers, S. J. (2007). Children and adolescents with autism exhibit reduced MEG steady-state gamma responses. Biol. Psychiatry 62, 192–197. doi: 10.1016/j.biopsych.2006.07.002
Keywords: dyspraxia, motor planning, event-related desynchronization, autism, praxis
Citation: Ewen JB, Lakshmanan BM, Pillai AS, McAuliffe D, Nettles C, Hallett M, Crone NE and Mostofsky SH (2016) Decreased Modulation of EEG Oscillations in High-Functioning Autism during a Motor Control Task. Front. Hum. Neurosci. 10:198. doi: 10.3389/fnhum.2016.00198
Received: 25 January 2016; Accepted: 19 April 2016;
Published: 06 May 2016.
Edited by:
John J. Foxe, University of Rochester School of Medicine and Dentistry, USAReviewed by:
Matthew K. Belmonte, The Com DEALL Trust, IndiaWilliam C. Gaetz, The Children’s Hospital of Philadelphia, USA
Copyright © 2016 Ewen, Lakshmanan, Pillai, McAuliffe, Nettles, Hallett, Crone and Mostofsky. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution and reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Joshua B. Ewen, ewen@kennedykrieger.org