Exploring the heterogeneity of neural social indices for genetically distinct etiologies of autism
© The Author(s). 2017
Received: 13 October 2016
Accepted: 10 May 2017
Published: 26 May 2017
Autism spectrum disorder (ASD) is a genetically and phenotypically heterogeneous disorder. Promising initiatives utilizing interdisciplinary characterization of ASD suggest phenotypic subtypes related to specific likely gene-disrupting mutations (LGDMs). However, the role of functionally associated LGDMs in the neural social phenotype is unknown.
In this study of 26 children with ASD (n = 13 with an LGDM) and 13 control children, we characterized patterns of mu attenuation and habituation as children watched videos containing social and nonsocial motions during electroencephalography acquisition.
Diagnostic comparisons were consistent with prior work suggesting aberrant mu attenuation in ASD within the upper mu band (10–12 Hz), but typical patterns within the lower mu band (8–10 Hz). Preliminary exploration indicated distinct social sensitization patterns (i.e., increasing mu attenuation for social motion) for children with an LGDM that is primarily expressed during embryonic development. In contrast, children with an LGDM primarily expressed post-embryonic development exhibited stable typical patterns of lower mu attenuation. Neural social indices were associated with social responsiveness, but not cognition.
These findings suggest unique neurophysiological profiles for certain genetic etiologies of ASD, further clarifying possible genetic functional subtypes of ASD and providing insight into mechanisms for targeted treatment approaches.
KeywordsAutism spectrum disorders (ASD) Likely gene-disrupting mutations Electroencephalography (EEG) Social cognition Mu rhythm attenuation Social perception Molecular subtyping
The significant etiologic and phenotypic heterogeneity of autism spectrum disorder (ASD)  has made it challenging to target underlying mechanisms of ASD pathology. Considering that more than 1000 genes have been implicated in ASD [1, 2], recent initiatives have targeted genetic pathways [3, 4] and rare de novo likely gene-disrupting mutations (LGDMs) . As such, a burgeoning “genetics-first” approach has been proposed to improve identification and characterization of genetic subtypes of individuals with ASD . For instance, genetics-first studies have identified phenotypically distinct subtypes of autism for CHD8  and DYRK1A [8, 9] based upon behavioral and physical features within both children and animal models. However, the relevant contribution of genetic risk to aspects of the ASD phenotype (i.e., social communicative impairments) is poorly understood, especially for low-functioning individuals with ASD.
Recent work supports social perception as a possible neural index related to the hallmark social deficits in ASD [10–15]. Although the neural indices have been targeted in relation to copy number variations, such as the 16p11.2 locus [16, 17], little is known about neural patterns associated with LGDMs, likely due to the wide range of variability of specific gene expression across LGDMs. Our objective was to examine patterns of neural heterogeneity associated with LGDMs by completing a series of diagnostic and genetics-guided analyses of social perception. We hypothesized a diagnostic approach would indicate atypical social perception in ASD, consistent with theories of social brain dysfunction in ASD . Then, as a preliminary exploration, we predicted that the heterogeneity associated with LGDMs would indicate potentially divergent patterns of social perception based upon LGDM function. Following work suggesting distinct functional roles for genes strongly expressed during embryonic development [5, 18], we tested children with and without a LGDM associated with embryonic development as a possible functional neurodevelopment pathway that contributes to a shared phenotype. There is a growing body of evidence suggesting early embryonic disruptions may be related to impairments in social behavior (e.g., lack of interest in conspecific proximity) [19, 20] and/or dysfunctional information encoding (i.e., intellectual or developmental delays) [21, 22]. The current study sought to add to this literature by addressing whether individuals with an embryonically expressed LGDM exhibit dysfunctional information habituation within the social domain.
We opted to target mu attenuation, which is specifically sensitive to detecting the movements associated with biological motion and is known as a reliable index of social perception in typical populations [23, 24]. Mu rhythm is typically defined as neural activity oscillations within the 8–12-Hz frequency range of electroencephalography (EEG) across electrodes above the sensorimotor cortex. During the observation and execution of biological motion, the underlying neural assemblies of the mu rhythm desynchronize [25, 26]. This desynchronization results in the reduction of oscillatory power (i.e., attenuation of the signal), with a greater reduction for conditions with social significance (i.e., biological motion relative to nonbiological motion) in children and adults; for a review, see .
In ASD, several studies suggest atypical mu attenuation in ASD (e.g., no discrimination for social relative to nonsocial observed motion) [12, 28–30] while other studies indicate no difference in ASD compared to typical controls [31–33]. Recent work by Dumas and colleagues  suggests that this discrepancy may be driven by the functional significance of the lower and upper mu rhythm bands. Notably, there is evidence that the lower mu rhythm (8–10 Hz) is more responsive to observed motion than the upper mu rhythm (10–12 Hz) , which may be indicative of bottom-up sensory processing [36, 37]. Previous studies also implicate that upper mu (or alpha) is more sensitive to top-down cognitive processing, such as self-monitoring within social contexts  or increasing cognitive demands [39, 40].
It is also possible that conflicting mu attenuation results reflect the underlying heterogeneity in ASD, potentially driven by genetic etiology. For instance, both disrupted social cognition and information habituation are associated with embryonically expressed LGDMs (e.g., ADNP , POGZ ). Yet, it is unclear whether social information habituation is also disrupted and the extent to which this profile is unique to children with an embryonically expressed LGDM. To date, only one study has tested the rate at which mu attenuation is modulated (i.e., habituates, sensitizes) in ASD . In that study, children with ASD and an ASD-associated deletion or duplication at the 16p11.2 locus demonstrated divergent dynamic patterns of mu attenuation providing additional insight into the relationship between ASD-associated copy number variations (CNVs) and social neural phenotypes.
This study sought to characterize social motion discrimination in ASD within the upper and lower mu bands continuously over time to capture dynamic neural social indices that may be associated with LGDMs expressed in embryonic development. We tested mu attenuation and habituation first via diagnostic comparisons between typically developing (TYP) and ASD children and, second, via genetics-guided comparisons between children with and without LGDMs expressed preferentially in embryonic development (LGDM E+ vs. LGDM E−). Based upon prior work [12, 28–30], we predicted a lack of social motion discrimination in ASD relative to TYP and anticipated no habituation to either condition in ASD, consistent with . We predicted that children with an embryonically expressed ASD-associated LGDM might have a more severely impacted social profile relative to LGDM E−, in part due to embryonic development as a (more) critical period for regulation of gene expression in support of brain development . Lastly, we evaluated relationships between the neural social indices and individual predictors of social and cognitive behavioral features to better assess the specificity of mu attenuation to capture social processing.
Participants and clinical procedures
N (n female)
See Additional file 1: Table S1 for full clinical characterization of ASD participants. Diagnoses of autism were confirmed using the Autism Diagnostic Interview-Revised (ADI-R) [50, 51] and the Autism Diagnostic Observation Schedule-2 (ADOS-2) [49, 52]. Verbal and nonverbal IQ (VIQ, NVIQ) was assessed using the Wechsler Abbreviated Scale of Intelligence  or the Differential Ability Scales-Second Edition , depending on age. Due to the severe intellectual disabilities within the ASD-LGDM cases, IQ ratio scores (n = 3) were substituted when standard IQ deviation scores were not available. The ASD-LGDM cases were more cognitively impaired than both comparison groups in VIQ and NVIQ, F(1,24)’s > 23.21, p’s < .001. Thus, VIQ and NVIQ were included in the statistical models to account for known variation in IQ and explicitly tested as part of our third objective. We recognize that cognitive differences between LGDM and comparison groups (TYP, ASD-NON) is a limitation of our study; however, evidence suggests that mu attenuation during passive social perception is not linked to cognitive abilities . In addition, the genetics-guided comparisons between LGDM E+ and LGDM E− were matched on age, IQ, autism severity (i.e., via the ADOS-2 score), and adaptive behavior (i.e., via the Vineland Adaptive Behavior Scales-2 ), F(1,11)’s < 2.15, p > .17.
Identification of genetic variants
Small-molecule molecular inversion probes (smMIPs)  were designed to the coding portions of CHD8, DSCAM, DYRK1A, GRIN2B, SCN2A, SETBP1, SETD2, ADNP, and MED13L with a 5-bp single-molecule tag using a scoring algorithm described previously  in order to identify single-nucleotide variants (SNVs) and insertions/deletions (INDELs). Oligonucleotides (IDT, Coralville, IA) were ordered, and probes were pooled at an equal molar ratio and phosphorylated (1X pool). After initial testing, poor-performing smMIPs were repooled and phosphorylated in either a 10X or a 50X pool. A final working probe pool was created by combining the three pools so that the final concentration of each smMIP in the 10X and 50X initial pools was a 10- or 50-fold excess relative to the 1X pool. Genomic DNA capture, exonuclease treatment, and PCR amplification of each library were performed as previously described  with 120 ng of genomic DNA input. smMIP concentration was based on a ratio of 800 copies of each MIP to each haploid genome copy, based on the 1X pool concentration. We pooled barcoded libraries together and purified the pools with 0.8X AMPure XP beads (Beckman Coulter, Brea, CA) according to the manufacturer’s protocol. Pools were quantified in duplicate using the Qubit dsDNA HS Assay (Life Technologies, Grand Island, NY). All samples were sequenced on an Illumina MiSeq (Reagent Kit 300V2) or HiSeq 2000 according to the manufacturer’s instructions. Sequencing reads were analyzed using the mipgen analysis pipeline as described previously . SNVs and INDELs were called using freebayes/0.9.14 and required a minimum of 8X coverage with a variant quality (QUAL) score greater than 20. Severe events (nonsense, frameshift, splice, and INDELs) were validated by Sanger sequencing.
In order to detect large CNVs, all samples (ASD-LGDM and ASD-NON) were run on custom genome-wide arrays (Agilent Technologies, Santa Clara, CA). Events for all ASD-LGDM samples have been previously published . See Additional file 2: Table S2 for full genetic characterization of ASD participants.
Social motion task
The objective of the social motion task was to examine the neural response to moving stimuli as it pertains specifically to social, biological agents more so than nonsocial, nonbiological objects. In the same paradigm and procedures as Hudac et al. , each child watched 12 total minutes of silent motion, alternating between conditions of social motion (i.e., hands clapping, animated person dancing), nonsocial motion (i.e., tubes swinging, animated ball bouncing), and no motion (i.e., the two empty backgrounds of motion videos). In this way, we can distinguish between the neural response for social and nonsocial motions, both relative to a baseline without motion. Each of the six 60-s videos was observed twice in one of two possible stimuli presentation orders. Between videos, participants were directed to take a break and the experimenter initiated the next video after confirmation that the participant was ready. Children were seated approximately 75 cm from a video monitor and were instructed to sit still and attend to the videos. Video stimuli were displayed using E-Prime 2.0 software (Psychology Software Tools, Inc., Pittsburgh, PA) at a size of 27 cm by 36.8 cm and subtended a visual angle of 20.4° by 27.6°.
Continuous electroencephalogram (EEG) was recorded from a high-density 128-channel geodesic net using Net Station 4.3.1 software integrated with a 200-series high-impedance amplifier (Electrical Geodesics, Inc,, Eugene, OR). Electrode impedances were below 50 kΩ to maximize the signal-to-noise ratio, within the standard range for high-impedance amplifiers. During collection, EEG signals were referenced to the vertex electrode, analog filtered (0.1 Hz high-pass, 100 Hz elliptical low-pass), amplified, and digitized with a sampling rate of 500 Hz. A photocell recorded and marked the precise onset time of each video. During acquisition, researchers observed and marked periods containing movement and/or improper attention (e.g., participant looking away).
Methodological decisions were consistent with our previous study  and standard practices for processing EEG data . Following data collection, continuous EEG was segmented into 2-s epochs starting with the onset of each 1-min video in order to generate 30 epochs per video. Epochs marked during acquisition as contaminated by movement or improper attention were removed. Automatic artifact detection rejected channels containing voltage shifts greater than 100 μV for each trial. In addition, trained research assistants reviewed and verified the automated artifact detection to ensure the data were sufficiently clean. If the channel was rejected for more than 50% of epochs (i.e., signifying poor data recording for that channel specifically), the channel was marked as a bad channel. After artifact detection, bad channels were corrected via interpolation from neighboring channels. Epochs were re-referenced to the average reference, excluding the rim channels due to the increased amount of noise from channels consistent with prior work (e.g., [58, 59]) in order to reduce the contribution of noise to the average reference.
Artifact-free EEG data by group and condition
TYP vs. ASD-LGDM
t(24) = −1.59, p = .13
t (24) = −2.55, p = .018
TYP vs. ASD-NON
t(24) = 1.03, p = .32
t (24) = 2.27, p = .032
ASD-NON vs. ASD-LGDM
t(24) = −.78, p = .44
t(24) = −.72, p = .48
To create neural indices of social perception, we computed power attenuation relative to the average baseline (no motion) for the first 30 artifact-free epochs in each condition (social, nonsocial) across central electrodes. Spectral power was calculated using fast Fourier transforms (FFTs) in MATLAB (version 7.12.0, R2011a; Natick, MA) on each 2-s epoch. Each power spectra was averaged across standard [12, 24] central electrodes clustered around the C3 (31, 32, 37, 38, 42, 43, 53, and 54) and C4 (80, 81, 87, 88, 94, 104, 105, and 106) positions. These electrodes are across to the sensorimotor region of the brain, thought to correspond to the mirror neuron system. Power attenuation was computed as the natural log of the ratio between the power of social motion or nonsocial motion epoch and the power of the individual’s average response during the no-motion condition. Subsequently, a value of zero represents no power attenuation relative to baseline [e.g., ln(condition/rest) = 0], and larger negative values represent greater power attenuation (i.e., social motion < rest). Only the first 30 artifact-free epochs in each condition (social, nonsocial) were included to ensure that each individual contributed an equal number of epochs. Power spectra included lower mu (8–10 Hz) and upper mu (10–12 Hz).
Data analysis strategies
All analyses were conducted via SAS 9.3 (SAS Institute) using restricted maximum likelihood (REML) and Satterthwaite denominator degrees of freedom. A series of multilevel models were generated using PROC MIXED to describe the variances and covariances of attenuation, separately for each analysis and the two mu bands. All models included a random intercept for each individual.
As we were interested in dynamic changes in power attenuation varying by condition across exposure to the video (i.e., trial order), the final models included fixed effects for condition (0 = nonsocial), ASD diagnosis (diagnostic comparison model) or genetic group (genetics-guided comparison model), time (0 = trial 30, the last epoch), and subsequent interactions. Subject predictors were included as fixed effects to test additional contribution by each child’s age (0 = 12 years), FSIQ (0 = 100), VIQ (0 = 100), NVIQ (0 = 100), and gender (0 = male). None of the subject predictors significantly contributed to the models, F’s < 3.14, p’s > .086. However, all subject predictor fixed effects remained in the model to account for potential contributions to group variances.
By modeling dynamic changes over time (i.e., trial order), we were able to measure the rate by which mu attenuation habituates (becomes more positive over time) or sensitizes (becomes more negative over time). Specifically, neural indices were measured as the power attenuation difference between social and nonsocial conditions with the Tukey correction.
Diagnostic comparisons (TYP and ASD)
Diagnostic comparison MLM results
F (1,4635) = 27.88, p < .0001
F (1,4635) = 9.43, p = .002
F(1,35) = 1.35, p = .254
F(1,34.9) = 0.66, p = .423
Condition by group
F (1,4635) = 4.75, p = .029
F (1,4635) = 6.32, p = .012
F (1,4635) = 19.64, p < .0001
F (1,4635) = 11, p < .001
Slope by condition
F(1,4635) = 3.21, p = .073
F (1,4635) = 3.86, p = .049
Slope by group
F (1,4635) = 4.85, p = .028
F(1,4635) = 2.43, p = .119
Slope by condition by group
F(1,4635) = 0.47, p = .495
F(1,4635) = 0.13, p = .717
F(1,33) = 0.24, p = .631
F(1,33) = 1.83, p = .186
F(1,33) = 1.66, p = .207
F(1,33) = 1.92, p = .175
F(1,33) = 1.16, p = .290
F(1,33) = 1.71, p = .200
F(1,33) = 2.07, p = .160
F(1,33) = 2.9, p = .098
Genetics-guided comparisons of the LGDM group (E− vs. E+)
Next, we sought to explore whether these patterns were consistent within a subsample of the ASD group with a known genetic etiology of ASD. In particular, we were interested in post hoc clustering comparisons based upon gene function given prior work targeting functional classes of LGDMs [5, 18]. To this extent, we evaluated mu attenuation related to social discrimination within the LGDM group to compare LGDMs associated with and without strong gene expression during embryonic development (LGDM E+ vs. LGDM E−).
LGDM comparison MLM results
F (1,1541) = 12.18, p < .001
F(1,1541) = 0.69, p = .406
F(1,9.33) = 0.24, p = .634
F(1,7.79) = 0.03, p = .865
Condition by group
F(1,1541) = 0.01, p = .914
F(1,1541) = 0.06, p = .809
F(1,1541) = 0.01, p = .904
F(1,1541) = 0.52, p = .469
Slope by condition
F (1,1541) = 4.14, p = .042
F(1,1541) = 0.04, p = .849
Slope by group
F (1,1541) = 4.11, p = .043
F(1,1541) = 2.98, p = .084
Slope by condition by group
F(1,1541) = 1.55, p = .213
F(1,1541) = 0.16, p = .687
F(1,7) = 0.71, p = .426
F(1,7) = 2.71, p = .144
F(1,7) = 1.6, p = .246
F(1,7) = 0.86, p = .385
F(1,7) = 0.09, p = .776
F(1,7) = 0.08, p = .785
F(1,7) = 1.56, p = .252
F(1,7) = 4.41, p = .074
Omnibus tests (see Table 4) indicated slopes differed by condition and by group. Both groups exhibited different social and nonsocial slopes within lower mu [LGDM E−: F(1,1541) = 4.89, p = .027; LGDM E+: F(1,1541) = 7.52, p = .0062]. Neither group habituated or sensitized to nonsocial motion (slopes = 0), and the LGDM E− group also did not habituate or sensitize to social motion (slope = 0). However, the LGDM E+ group sensitized to social motion (slope = −.007), such that condition discrimination initially indicated more mu attenuation for nonsocial motion and ended with more mu attenuation for social motion at the end of the experiment. In other words, similar to typical development, both LGDM groups exhibited more social mu attenuation by the end of the experiment, but this effect was only shown after a delayed amount of time in the LGDM E+ group.
Relationships between neural indices and individual differences
Social impairments are a hallmark of ASD, yet phenotypic and genetic heterogeneity is thought to contribute to discrepant evidence in the literature. We explore a neural mechanism associated with ASD by considering patterns of social discrimination as measured by mu attenuation over time for children with different genetic etiologies. From the diagnostic comparisons, we show aberrant patterns of mu attenuation in ASD are specific to the upper mu band, while the lower mu band reflects less atypical patterns, consistent with prior work . The dynamic patterns indicate that children with ASD show an increasing lower mu difference between social and nonsocial motions, which may help resolve diagnostic inconsistencies within the literature. For instance, prior evidence of atypical mu attenuation in ASD between observed motion conditions (i.e., social relative to nonsocial motion, as in the current study) [12, 28–30] has relied on individual averages. Our findings suggest that the discrimination pattern may not be evident if there are too few trials (i.e., before children with ASD habituated to nonsocial motion observations). Although it is a concern that children with ASD did contribute fewer trials than the typical controls, the eventual condition differentiation in ASD (i.e., noted by approximately trial 15 in Fig. 2) indicates that our study had a sufficient number of trials.
Mu attenuation has been proposed to reflect a human corollary to the mirror neuron system [24, 61], which describes activation recorded over the sensorimotor cortex during both action execution and observation of human actions. Although it is possible that mu attenuation reflects the conductance of occipital or posterior alpha rhythm more broadly (i.e., responsivity to general motion information) , our results indicate differentiation of social and nonsocial motions. One working hypothesis of ASD suggests mirror neuron system deficits that disrupt neural correlates supporting the action/observation system, subsequently eliciting atypical mu attenuation . Aligned with this theory, other evidence suggest that atypical functioning of the mirror neuron system may lead to a downstream effect of poor imitative abilities  or disrupted higher order social cognitive abilities (i.e., theory of mind) . However, it is important to note that similar to prior work by Dumas and colleagues, we found mu attenuation diagnostic differences within the upper mu band (10–12 Hz), but no group difference within the lower mu band (8–10 Hz). This is consistent with prior work suggesting that this lower frequency may reflect primary sensory processing [36, 37] that habituates over the course of the exposure. Yet, sensory processing of biological motion occurred more rapidly in the TYP group compared to longer processing in the ASD group, perhaps indicative of functional connectivity reductions related to social cognition . Our results offer further evidence of atypical mu attenuation patterns in ASD, although unique neural mechanisms underlying atypical social discrimination may be derived from specific genetic etiologies. In other words, a mirror neuron hypothesis may indeed describe a subset of children with ASD, while a more general, distributed network of neural correlates may be impacted in other ASD subgroups.
As part of a preliminary analysis, we examined the neural social indices associated with different functional genetic roles of LGDMs as a first step to explore a possible shared neural social phenotype. We implemented a post hoc clustering strategy in order to examine potential convergent pathways between LGDMs that are and are not functionally expressed during embryonic development [18, 66]. The choice to cluster LGDMs around functional expression during embryonic development is based on early genetic regulatory control supporting regional differentiation within the embryonic brain [67, 68], including key social neural structures (e.g., amygdala). We had predicted that the LGDM within the embryonic development group would have a more severely disrupted neural social index due to evidence from animal and human models indicating significant impairments related to social behavior [19, 20] and/or information encoding [21, 22]. The results indicated that children with an LGDM primarily expressed during embryonic development exhibit sensitization of lower mu attenuation to social motion. In other words, these children initially exhibited more mu attenuation for nonsocial motion, but eventually demonstrate more for social motion. This pattern was distinct from children with an LGDM not primarily expressed during embryonic development that exhibited greater lower mu attenuation discrimination throughout the entire experiment (i.e., greater mu attenuation to social than nonsocial motion beginning at the first few trials).
Our results suggest that social motion perception may be conserved despite early genetic disruption, though the delayed processing supports the notion of potentially delayed information processing. It is important to note that this delay was specific to the social motion condition (increasing neural response over time) but not the nonsocial motion condition (i.e., no change over time), which may help clarify the mechanism by which prior models [19, 20] derive impaired social behavior. An interpretation of the results may be that children with an LGDM primarily expressed during embryonic development are increasing their attention to, or interest in, social stimuli after an initial period, which may reflect a delayed social engagement (e.g., motivation or salience). One explanation may be that the impact of embryonic genes on social perception is greater , suggesting that functional timing of genetic expression may differentially affect the neural social phenotype. Importantly, these findings align with genetics research indicating that ASD genes converge on several select pathways [70, 71], which may help to further explain the underlying neural social heterogeneity.
An important limitation of the current study is the continued genetic heterogeneity despite functionally classifying the expression of LGDM within early development. Within our LGDM groups, there are only several children with a shared LGDM (i.e., SETD2, n = 2; DYRK1A, n = 2; CHD8, n = 2). Thus, the discoveries of this work are not to be taken as firm conclusions, but rather considered in order to motivate and guide continued use of a genetics-first approach to elucidate potential etiological mechanisms of ASD. For instance, most of the children within the early embryonic LGDM group exhibit the social sensitization pattern described here (six out of eight cases; see Additional file 3: Figure S1 for individual patterns), except for one child with MED13L and one child with DYRK1A. In part, this qualitative finding is consistent with the overall group clustering approach indicating delayed social processing, suggesting a potential neural index associated with this particular genetic etiology. However, the specificity for specific LGDMs may be poor, considering that only two out of three children with a DYRK1A LGDM exhibited this pattern. Similar to prior work linking core social symptoms to biomarkers of ASD [11, 72–76], we encourage the use of this data as a way to bridge the gap between genetic and phenotypic characterization as a means to facilitate the discovery of ASD etiological mechanisms and accelerate progress for ASD therapeutic interventions.
It may be surprising that our task elicited mu attenuation during nonsocial motion observation (i.e., ball bouncing, tubes swinging) that is not biological and subsequently should not be simulated within the action/observation system. However, to a large extent, the majority of studies implementing mu attenuation as an outcome utilized comparisons between self-executed, social observed, and nonsocial observed motion. It may be the case that by engaging the motor execution system during these tasks, the threshold for the action/observation system is elevated, reducing the amount of mu attenuation for nonsocial comparisons. In fact, neural regions implicated in mu suppression during execution vs. observation  involve regions that also play a role in general motion perception, including the occipital, premotor, and somatosensory cortices. Moreover, this study replicated prior work with this same task that indicated a modest degree of mu attenuation to nonsocial motion, in addition to social motion . We posit that our task measured more globally distributed neural differences between social and nonsocial motions compared to other tasks that have used self-initiated actions to target the premotor cortex. Of note, this passive viewing task is more conducive for children with reduced capacity for following behavioral instructions (i.e., to make self-initiated motions), while still providing a robust neural index, which specifies individual patterns.
The neural social indices were correlated with features of social cognition (i.e., social responsiveness), particularly with the lower mu band. This finding is compelling evidence that these indices accurately capture subtle levels of social impairments in vivo, as opposed to relying on parental reports (e.g., SRS-2). Additionally, average patterns of mu attenuation were unaffected by general cognition, despite drastic cognitive differences for children with a LGDM. Although this may not negate a contributory role of cognitive ability for higher-order operations related to social motion (e.g., action prediction), this evidence from this study suggests that motion perception is intact for children with lower cognitive abilities (i.e., cognitive scores under 50). Much of the existing research investigating neural social indices is restricted to children and adults with moderate to average cognitive capabilities. The majority of ASD-LGDM cases with low verbal IQ show typical mu attenuation patterns (i.e., greater for social motion in five out of eight cases with verbal IQ < 50). Taken together, these neural social indices can provide a robust characterization of the underlying neural mechanisms supporting social cognition, regardless of level of cognitive function, thereby improving our understanding of the social phenotype.
This study is the first to use a genetics-first approach to explore the genetic etiologies of autism associated with severe LGDMs in the context of neural social indices. Our use of a unique statistical method to measure ongoing dynamic changes associated with social motion perception demonstrates the utility of this method to better understand underlying processes relevant to ASD and LGDMs. Although this study is limited by a small sample size and thus should be considered exploratory, the analysis of neural social phenotypes based on functional clustering offers a promising approach for narrowing in on convergent pathways that may reflect shared phenotypes and provide insight for targeted treatment [5, 18]. Future research will need to take into account the variety and combination of genetic functional roles. Ongoing efforts to recruit a larger, more genetically homogenous group will help target specific functional outcomes during early childhood and adolescence. However, due to the rarity of this population, these preliminary results are informative and can help guide future research by better describing the functional processes during social motion perception and similar processes that are impaired in ASD.
In this study, we demonstrated distinct neural social indices for genetic etiologies of ASD, providing critical insight into the underlying mechanisms of ASD pathology. A unique mechanism was identified for children with ASD and genetic etiology associated with early embryonic development, based upon level of mu attenuation related to social discrimination and patterns over time (i.e., habituation). Our findings implicate genetic heterogeneity as a possible reason for divergent findings in the literature and distinguish the manner by which neural social indices differ between groups and over time. We emphasize the need to continue to discover how phenotypic profiles align within children in specific genotypic subgroups of ASD. Taken together, we predict that future work pursuing phenotypic characterization via the integration of genetic, neural, and behavioral information will continue to inform our understanding of ASD subtypes and will have broad implications for our ability to adopt precision medicine strategies.
We would like to thank the children and families for their participation in this study. We are grateful to all of the families at the participating Simons Simplex Collection (SSC) sites, as well as the principal investigators (A. Beaudet, R. Bernier, J. Constantino, E. Cook, E. Fombonne, D. Geschwind, R. Goin-Kochel, E. Hanson, D. Grice, A. Klin, D. Ledbetter, C. Lord, C. Martin, D. Martin, R. Maxim, J. Miles, O. Ousley, K. Pelphrey, B. Peterson, J. Piggot, C. Saulnier, M. State, W. Stone, J. Sutcliffe, C. Walsh, Z. Warren, E. Wijsman). We appreciate obtaining access to phenotypic data on SFARI Base. E.E.E. is an investigator of the Howard Hughes Medical Institute. The computerized dancer stimuli were provided courtesy of Nick Neave and Kristofor McCarty at Northumbria University.
Research reported in this publication was supported by the National Institute of Mental Health (MH100047, Bernier; MH10028, Webb/Bernier; and MH101221, Eichler), by the National Institute of Child Health and Human Development to the University of Washington’s Center on Human Development and Disability (U54 HD083091), and in part by a grant from the Simons Foundation (SFARI 303241, Eichler). H.A.F.S. was supported in part by the NHGRI Interdisciplinary Training in Genome Science Grant (T32HG00035).
Availability of data and materials
Approved researchers can obtain the SSC population dataset described in this study (http://sfari.org/resources/simons-simplex-collection) by applying at https://base.sfari.org.
RAB contributed to the conceptualization and methodology. CMH and HAFS contributed to the formal analysis. TDD, AK, SF, and EN contributed to the investigation. RAB and EEE contributed to the resources. CMH and RAB wrote the original draft. All authors performed the review and editing. CMH contributed to the visualization. SJW, EEE, and RAB supervised the study. EEE and RAB contributed to the project administration and funding acquisition. All authors read and approved the final manuscript.
E.E.E. is on the scientific advisory board (SAB) of DNAnexus, Inc. and is a consultant for Kunming University of Science and Technology (KUST) as part of the 1000 China Talent Program. All other authors do not have competing interests.
Consent for publication
Ethics approval and consent to participate
All research procedures conformed to regulations in accordance with the local ethical review board at the University of Washington, which approved this project. Written informed consent was obtained from each adult participant or parental representative(s). All children verbally assented to participate in the procedures, and written assent was obtained from children with a mental age of 7 or greater.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Geschwind DH. Advances in autism. Annu Rev Med. 2009;60:367–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Betancur C. Etiological heterogeneity in autism spectrum disorders: more than 100 genetic and genomic disorders and still counting. Brain Res. 2011;1380:42–77.View ArticlePubMedGoogle Scholar
- Bill BR, Geschwind DH. Genetic advances in autism: heterogeneity and convergence on shared pathways. Curr Opin Genet Dev. 2009;19:271–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Geschwind DH, State MW. Gene hunting in autism spectrum disorder: on the path to precision medicine. Lancet Glob Health. 2015;14(11):1109–20.Google Scholar
- Iossifov I, O’Roak BJ, Sanders SJ, Ronemus M, Krumm N, Levy D, et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature. 2014;515:216–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Higdon R, Earl RK, Stanberry L, Hudac CM, Montague E, Stewart E, et al. The promise of multi-omics and clinical data integration to identify and target personalized healthcare approaches in autism spectrum disorders. OMICS. 2015;19:197–208. Available from: http://online.liebertpub.com/doi/https://doi.org/10.1089/omi.2015.0020.View ArticlePubMedPubMed CentralGoogle Scholar
- Bernier R, Golzio C, Xiong B, Stessman HA, Coe BP, Penn O, et al. Disruptive CHD8 mutations define a subtype of autism early in development. Cell. 2014;158:1–14. Elsevier Inc.View ArticleGoogle Scholar
- van Bon B, Hoischen A, Hehir-Kwa J, de Brouwer A, Ruivenkamp C, Gijsbers A, et al. Intragenic deletion in DYRK1A leads to mental retardation and primary microcephaly. Clin Genet. 2011;79:296–9.View ArticlePubMedGoogle Scholar
- van Bon BWM, Coe BP, Bernier R, Green C, Gerdts J, Witherspoon K, et al. Disruptive de novo mutations of DYRK1A lead to a syndromic form of autism and ID. Molecular Psychiatry. 2016;21:126–32. Nature Publishing Group.View ArticlePubMedGoogle Scholar
- Wechsler D. Wechsler Abbreviated Scale of Intelligence (WASI). San Antonio: The Psychological Corporation; 1999.Google Scholar
- Kaiser MD, Hudac CM, Shultz S, Lee SM, Cheung C, Berken AM, et al. Neural signatures of autism. Proc Natl Acad Sci. 2010;107:21223–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Bernier R, Dawson G, Webb S, Murias M. EEG mu rhythm and imitation impairments in individuals with autism spectrum disorder. Brain Cogn. 2007;64:228–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Webb SJ, Merkle K, Murias M, Richards T, Aylward E, Dawson G. ERP responses differentiate inverted but not upright face processing in adults with ASD. Soc Cogn Affect Neurosci. 2012;7:578–87. Oxford University Press.View ArticlePubMedGoogle Scholar
- Gordon I, Vander Wyk BC, Bennett RH, Cordeaux C, Lucas MV, Eilbott JA, et al. Oxytocin enhances brain function in children with autism. Proc Natl Acad Sci. 2013;110:20953–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Pelphrey KA, Shultz S, Hudac CM, Vander Wyk BC. Research review: Constraining heterogeneity: the social brain and its development in autism spectrum disorder. J Child Psychol Psychiatry. 2011;52:631–44. Available from: http://doi.wiley.com/https://doi.org/10.1111/j.1469-7610.2010.02349.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Berman JI, Chudnovskaya D, Blaskey L, Kuschner E, Mukherjee P, Buckner R, et al. Abnormal auditory and language pathways in children with 16p11.2 deletion. Neuroimage Clin. 2015;9:50–7. Elsevier B.V.View ArticlePubMedPubMed CentralGoogle Scholar
- Hudac CM, Kresse A, Aaronson B, DesChamps TD, Webb SJ, Bernier RA. Modulation of mu attenuation to social stimuli in children and adults with 16p11. 2 deletions and duplications. J Neurodevelopmental Dis. 2015;7(1):25.Google Scholar
- Iossifov I, Levy D, Allen J, Ye K, Ronemus M, Lee Y-H, et al. Low load for disruptive mutations in autism genes and their biased transmission. Proc Natl Acad Sci. 2015;112(41):E5600-7.Google Scholar
- Rampersad M, Gerlai R. Impairment of social behaviour persists two years after embryonic alcohol exposure in zebrafish: a model of fetal alcohol spectrum disorders. Behav Brain Res. 2015;292:102–8. Elsevier B.V.View ArticlePubMedPubMed CentralGoogle Scholar
- Selemon LD, Zecevic N. Schizophrenia: a tale of two critical periods for prefrontal cortical development. Translational Psychiatry. 2015;5:e623–11. Nature Publishing Group.View ArticlePubMedPubMed CentralGoogle Scholar
- Malishkevich A, Amram N, Hacohen-Kleiman G, Magen I, Giladi E, Gozes I. Activity-dependent neuroprotective protein (ADNP) exhibits striking sexual dichotomy impacting on autistic and Alzheimer’s pathologies. Nat Publ Group. 2015;5:e501–9.Google Scholar
- Stessman HAF, Willemsen MH, Fenckova M, Penn O, Hoischen A, Xiong B, et al. Disruption of POGZ is associated with intellectual disability and autism spectrum disorders. Am J Hum Genet. 2016;98:541–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Muthukumaraswamy SD, Johnson BW. Changes in rolandic mu rhythm during observation of a precision grip. Psychophysiol. 2004;41:152–6. Available from: http://doi.wiley.com/https://doi.org/10.1046/j.1469-8986.2003.00129.x.View ArticleGoogle Scholar
- Muthukumaraswamy SD, Johnson BW, McNair NA. Mu rhythm modulation during observation of an object-directed grasp. Cogn Brain Res. 2004;19:195–201.View ArticleGoogle Scholar
- Arroyo S, Lesser RP, Gordon B, Uematsu S, Jackson D, Webber R. Functional significance of the mu rhythm of human cortex: an electrophysiologic study with subdural electrodes. Electroencephalogr Clin Neurophysiol. 1993;87:76–87. Available from: http://linkinghub.elsevier.com/retrieve/pii/001346949390114B.View ArticlePubMedGoogle Scholar
- Pfurtscheller G, Neuper C. Motor imagery activates primary sensorimotor area in humans. Neurosci Lett. 1997;239:65–8.View ArticlePubMedGoogle Scholar
- Fox NA, Bakermans-Kranenburg MJ, Yoo KH, Bowman LC, Cannon EN, Vanderwert RE, et al. Assessing human mirror activity with EEG mu rhythm: a meta-analysis. Psychol Bull. 2016;142:291–313. American Psychological Association.View ArticlePubMedGoogle Scholar
- Pineda JA, Brang D, Hecht E, Edwards L, Carey S, Bacon M, Futagaki C, Suk D, Tom J, Birnbaum C, Rork A. Positive behavioral and electrophysiological changes following neurofeedback training in children with autism. Res Autism Spectr Disord. 2008;2(3):557–81.View ArticleGoogle Scholar
- Oberman LM, McCleery JP, Hubbard EM, Bernier R, Wiersema JR, Raymaekers R, et al. Developmental changes in mu suppression to observed and executed actions in autism spectrum disorders. Soc Cogn Affect Neurosci. 2013;8:300–4. Oxford University Press.View ArticlePubMedGoogle Scholar
- Oberman LM, Hubbard EM, McCleery JP, Altschuler EL, Ramachandran VS, Pineda JA. EEG evidence for mirror neuron dysfunction in autism spectrum disorders. Brain Res Cogn Brain Res. 2005;24:190–8.View ArticlePubMedGoogle Scholar
- Raymaekers R, Wiersema JR, Roeyers H. EEG study of the mirror neuron system in children with high functioning autism. Brain Res. 2009;1304:113–21.View ArticlePubMedGoogle Scholar
- Fan Y-T, Decety J, Yang C-Y, Liu J-L, Cheng Y. Unbroken mirror neurons in autism spectrum disorders. J Child Psychol Psychiatry. 2010;51:981–8.View ArticlePubMedGoogle Scholar
- Bernier R, Aaronson B, McPartland J. The role of imitation in the observed heterogeneity in EEG mu rhythm in autism and typical development. Brain Cogn. 2013;82:69–75.View ArticlePubMedGoogle Scholar
- Dumas G, Soussignan R, Hugueville L, Martinerie J, Nadel J. Revisiting mu suppression in autism spectrum disorder. Brain Res. 2014;1585:108–19.View ArticlePubMedGoogle Scholar
- Frenkel-Toledo S, Bentin S, Perry A, Liebermann DG, Soroker N. Dynamics of the EEG power in the frequency and spatial domains during observation and execution of manual movements. Brain Res. 2013;1509:43–57.View ArticlePubMedGoogle Scholar
- Sacchet MD, LaPlante RA, Wan Q, Pritchett DL, Lee AKC, Hamalainen M, et al. Attention drives synchronization of alpha and beta rhythms between right inferior frontal and primary sensory neocortex. J Neurosci. 2015;35:2074–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Başar E. A review of alpha activity in integrative brain function: fundamental physiology, sensory coding, cognition and pathology. Int J Psychophysiol. 2012;86:1–24. Elsevier B.V.View ArticlePubMedGoogle Scholar
- Naeem M, Prasad G, Watson DR, Kelso JAS. Electrophysiological signatures of intentional social coordination in the 10-12 Hz range. Neuroimage. 2012;59:1795–803.View ArticlePubMedGoogle Scholar
- Fink A, Schwab D, Papousek I. Sensitivity of EEG upper alpha activity to cognitive and affective creativity interventions. Int J Psychophysiol. 2011;82:233–9.View ArticlePubMedGoogle Scholar
- Fink A, Grabner RH, Neuper C, Neubauer AC. EEG alpha band dissociation with increasing task demands. Cogn Brain Res. 2005;24:252–9.View ArticleGoogle Scholar
- Patapoutian A, Reichardt LF. Roles of Wnt proteins in neural development and maintenance. Curr Opin Neurol. 2000;10:392–9. NIH Public Access.View ArticleGoogle Scholar
- O'Roak BJ, Stessman HA, Boyle EA, Witherspoon KT, Martin B, Lee C, et al. Recurrent de novo mutations implicate novel genes underlying simplex autism risk. Nat Commun. 2014;5:5595.View ArticlePubMedPubMed CentralGoogle Scholar
- O’Roak BJ, Deriziotis P, Lee C, Vives L, Schwartz JJ, Girirajan S, et al. Exome sequencing in sporadic autism spectrum disorders identifies severe de novo mutations. Nat Genet. 2011;43:585–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Vandeweyer G, Helsmoortel C, Van Dijck A, Vulto-van Silfhout AT, Coe BP, Bernier R, et al. The transcriptional regulator ADNP links the BAF (SWI/SNF) complexes with autism. Kosho T, Miyake N, editors. Am J Med Genet C Semin Med Genet. 2014;166:315–26Google Scholar
- Coe BP, Witherspoon K, Rosenfeld JA, van Bon BWM, Vulto-van Silfhout AT, Bosco P, et al. Refining analyses of copy number variation identifies specific genes associated with developmental delay. Nat Genet. 2014;46:1063–71. Available from: http://dx.doi.org/https://doi.org/10.1038/ng.3092.View ArticlePubMedPubMed CentralGoogle Scholar
- Lumish HS, Wynn J, Devinsky O, Chung WK. Brief report: SETD2 mutation in a child with autism, intellectual disabilities and epilepsy. J Autism Dev Disord. 2015;45(11):3764–70.View ArticlePubMedGoogle Scholar
- O’Roak BJ, Vives L, Girirajan S, Karakoc E, Krumm N, Coe BP, et al. Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations. Nature. 2013;485:246–50. Nature Publishing Group.View ArticleGoogle Scholar
- Weiss LA, Escayg A, Kearney JA, Trudeau M, MacDonald BT, Mori M, et al. Sodium channels SCN1A, SCN2A and SCN3A in familial autism. Mol Psychiatry. 2003;8:186–94.View ArticlePubMedGoogle Scholar
- Constantino JN, Przybeck T, Friesen D, Todd RD. Reciprocal social behavior in children with and without pervasive developmental disorders. LWW; 2000;21:2–11Google Scholar
- Elliott CD. Differential Ability Scales-ll. Pearson: San Antonio; 2007.Google Scholar
- Lord C, Rutter M, Le Couteur A. Autism Diagnostic Interview-Revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders. J Autism Dev Disord. 1994;24:659–85.View ArticlePubMedGoogle Scholar
- Lord C, Rutter M, Goode S, Heemsbergen J, Jordan H, Mawhood L, et al. Autism diagnostic observation schedule: a standardized observation of communicative and social behavior. J Autism Dev Disord. 1989;19:185–212.View ArticlePubMedGoogle Scholar
- Sparrow SS, Cicchetti DV, Balla DA, Doll EA. Vineland-II, Vineland Adaptive Behavior Scales: Teacher Rating Form. 2006.Google Scholar
- Hiatt JB, Pritchard CC, Salipante SJ, O’Roak BJ, Shendure J. Single molecule molecular inversion probes for targeted, high-accuracy detection of low-frequency variation. Genome Res. 2013;23:843–54. Cold Spring Harbor Lab.View ArticlePubMedPubMed CentralGoogle Scholar
- Boyle EA, O’Roak BJ, Martin BK, Kumar A, Shendure J. MIPgen: optimized modeling and design of molecular inversion probes for targeted resequencing. Bioinformatics. 2014;30:2670–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Girirajan S, Dennis MY, Baker C, Malig M, Coe BP, Campbell CD, et al. Refinement and discovery of new hotspots of copy-number variation associated with autism spectrum disorder. Am J Hum Genet. 2013;92:221–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Webb SJ, Bernier R, Henderson HA, Johnson MH, Jones EJH, Lerner MD, et al. Guidelines and best practices for electrophysiological data collection, analysis and reporting in autism. J Autism Dev Disord. 2015;45:425–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Neuhaus E, Jones EJH, Barnes K, Sterling L, Estes A, Munson J, et al. The relationship between early neural responses to emotional faces at age 3 and later autism and anxiety symptoms in adolescents with autism. J Autism Dev Disord. 2016;46:1–14. Springer US.View ArticleGoogle Scholar
- Ding J, Sperling G, Srinivasan R. Attentional modulation of SSVEP power depends on the network tagged by the flicker frequency. Cerebral Cortex. 2006;16:1016–29.View ArticlePubMedGoogle Scholar
- Hobson HM, Bishop DVM. Mu suppression—a good measure of the human mirror neuron system? Cortex. 2016;82:290–310.View ArticlePubMedPubMed CentralGoogle Scholar
- Pineda JA. The functional significance of mu rhythms: translating “seeing” and “hearing” into “doing”. Brain Res Brain Res Rev. 2005;50:57–68.View ArticlePubMedGoogle Scholar
- Varela FJ, Toro A, John ER, Schwartz EL. Perceptual framing and cortical alpha-rhythm. Neuropsychologia. 1981;19:675–86.View ArticlePubMedGoogle Scholar
- Pineda JA, Hecht E. Mirroring and mu rhythm involvement in social cognition: are there dissociable subcomponents of theory of mind? Biol Psychol. 2009;80:306–14. Available from: http://linkinghub.elsevier.com/retrieve/pii/S0301051108002287.View ArticlePubMedGoogle Scholar
- Hamilton AF De C, Brindley RM, Frith U. Imitation and action understanding in autistic spectrum disorders: how valid is the hypothesis of a deficit in the mirror neuron system? Neuropsychologia. 2007;45:1859–68.View ArticleGoogle Scholar
- von dem Hagen EA, Stoyanova RS, Baron-Cohen S, Calder AJ. Reduced functional connectivity within and between “social” resting state networks in autism spectrum conditions. Soc Cogn Affect Neurosci. 2012;8:694–701.View ArticleGoogle Scholar
- Iossifov I, O’Roak BJ, Sanders SJ, Ronemus M, Krumm N, Levy D, et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature. Nature Publishing Group; 2014;515:216–21.Google Scholar
- Puelles L, Kuwana E, Puelles E, Bulfone A, Shimamura K, Keleher J, et al. Pallial and subpallial derivatives in the embryonic chick and mouse telencephalon, traced by the expression of the genes Dlx-2, Emx-1, Nkx-2.1, Pax-6, and Tbr-1. J Comp Neurol. 2000;424:409–38.View ArticlePubMedGoogle Scholar
- Puelles L, Rubenstein JLR. Forebrain gene expression domains and the evolving prosomeric model. Trends Neurosci. 2003;26:469–76.View ArticlePubMedGoogle Scholar
- Chang J, Gilman SR, Chiang AH, Sanders SJ, Vitkup D. Genotype to phenotype relationships in autism spectrum disorders. Nat Neurosci. 2014;18:191–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Robinson EB, St Pourcain B, Anttila V, Kosmicki JA, Bulik-Sullivan B, Grove J, et al. Genetic risk for autism spectrum disorders and neuropsychiatric variation in the general population. Nat Genet. 2016;48:552–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Krishnan A, Zhang R, Yao V, Theesfeld CL, Wong AK, Tadych A, et al. Genome-wide prediction and functional characterization of the genetic basis of autism spectrum disorder. Nat Neurosci. 2016;19:1454–62.View ArticlePubMedGoogle Scholar
- Veenstra-VanderWeele J, Blakely RD. Networking in autism: leveraging genetic, biomarker and model system findings in the search for new treatments. Neuropsychopharmacology. 2012;37:196–212. Nature Publishing Group.View ArticlePubMedGoogle Scholar
- Key AP, Ibanez LV, Henderson HA, Warren Z, Messinger DS, Stone WL. Positive affect processing and joint attention in infants at high risk for autism: an exploratory study. J Autism Dev Disord. 2015;45:4051–62. Springer US.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang F, Zhu AJ, Lajiness-O'Neill R, Bowyer S. Functional network connectivity: possible biomarker for autism spectrum disorders (ASD). 2015.Google Scholar
- Neuhaus E, Bernier R, Beauchaine TP. Brief report: social skills, internalizing and externalizing symptoms, and respiratory sinus arrhythmia in autism. J Autism Dev Disord. 2014;44:730–7.View ArticlePubMedGoogle Scholar
- Spencer MD, Holt RJ, Chura LR, Suckling J, Calder AJ, Bullmore ET, et al. A novel functional brain imaging endophenotype of autism: the neural response to facial expression of emotion. Transl Psychiatry. 2011;1:e19–7. Nature Publishing Group.View ArticlePubMedPubMed CentralGoogle Scholar
- Braadbaart L, Williams JHG, Waiter GD. Do mirror neuron areas mediate mu rhythm suppression during imitation and action observation? Int J Psychophysiol. 2013;89:99–105. Elsevier B.V.View ArticlePubMedGoogle Scholar