Altered functional connectivity of the amygdaloid input nuclei in adolescents and young adults with autism spectrum disorder: a resting state fMRI study
© Rausch et al. 2016
Received: 10 April 2015
Accepted: 7 December 2015
Published: 28 January 2016
Amygdala dysfunction is hypothesized to underlie the social deficits observed in autism spectrum disorders (ASD). However, the neurobiological basis of this hypothesis is underspecified because it is unknown whether ASD relates to abnormalities of the amygdaloid input or output nuclei. Here, we investigated the functional connectivity of the amygdaloid social-perceptual input nuclei and emotion-regulation output nuclei in ASD versus controls.
We collected resting state functional magnetic resonance imaging (fMRI) data, tailored to provide optimal sensitivity in the amygdala as well as the neocortex, in 20 adolescents and young adults with ASD and 25 matched controls. We performed a regular correlation analysis between the entire amygdala (EA) and the whole brain and used a partial correlation analysis to investigate whole-brain functional connectivity uniquely related to each of the amygdaloid subregions.
Between-group comparison of regular EA correlations showed significantly reduced connectivity in visuospatial and superior parietal areas in ASD compared to controls. Partial correlation analysis revealed that this effect was driven by the left superficial and right laterobasal input subregions, but not the centromedial output nuclei.
These results indicate reduced connectivity of specifically the amygdaloid sensory input channels in ASD, suggesting that abnormal amygdalo-cortical connectivity can be traced down to the socio-perceptual pathways.
KeywordsAmygdala Autism spectrum disorder Centromedial Connectivity Input-output Laterobasal Nuclei Social perception Superficial
Autism spectrum disorders (ASD) are a group of neurodevelopmental disorders characterized by severe impairments of reciprocal social interaction and verbal and nonverbal communication and by repetitive and stereotyped behaviors [1, 2]. Structural and functional neuroimaging studies have linked a number of brain structures to ASD symptoms [3–6], one of which is the amygdala. The amygdala theory of autism describes this structure as potential key component in the pathogenesis of ASD [7, 8], since it is involved inwfi various aspects of the social brain, such as social cognition, emotion recognition, socio-communicative perception, and the regulation of emotional responses . In line with this, individuals with ASD tend to show abnormal volume enlargements of the amygdala [10, 11] and have overactive amygdalae in response to mildly aversive stimuli  and faces , while symptom severity in ASD has been found to correlate with amygdala size [10, 14, 15]. Although amygdala impairments likely relate to pathophysiological socio-emotional processes, its subregion-specific amygdalo-cortical abnormalities have not been stratified in ASD. The aim of this study is to advance our understanding of the pathway-specific amygdala involvement and discern sensory input and response output channels separately.
With the last decade’s paradigm shift in neuroimaging from activity assessment within brain structures to connectivity within neural networks, increasing evidence supports the notion of atypical large-scale neural connectivity in ASD. Some authors hypothesized that the brain in ASD is characterized by long distance underconnectivity and local overconnectivity . Indeed, a number of studies reported large-scale underconnectivity with decreased structural, functional, and interhemispheric connectivity [4, 17, 18], and a few functional magnetic resonance imaging (fMRI) studies found local overconnectivity patterns [19–23]. To date, large-scale functional connectivity abnormalities in the autistic brain have been described most conclusively with respect to hyper- or hyposensitivity to perceptual stimuli , which may in part occur due to a lack of sensory integration in ASD . Abnormal sensory processing has primarily been found in the auditory [25–27] and visual system [28–30]. Because social functioning requires the selection and integration of many socio-perceptual stimuli simultaneously, perceptual processing abnormalities may in part account for some of the social difficulties in ASD.
Previous fMRI studies on the role of the amygdala in ASD have generally treated it as a single structure, while in fact it is comprised of at least 13 functionally and structurally distinct nuclei, in which three major input and output units can be discerned [31–33]: the centromedial (CM), laterobasal (LB), and superficial (SF) nuclei. Prior work in animals identified the centromedial part as an output area, which regulates cardiovascular control via projections to the brainstem, cerebellum, and hypothalamus . More specifically, the CM subregion generates ascending projections via the forebrain throughout the cortex and descending projections via the hypothalamus to the brainstem . Via these complex pathways, the CM subregion is thought to modulate autonomic, somatic, and endocrine responses to facilitate appropriate behavioral outcome . A recent study in humans mapped projections from the cortical social "aversion network", a network that is situated around the anterior cingulated cortex, onto a functionally defined amygdaloid subdivision that corresponds to the CM subcompartment . This subdivision may therefore be associated with emotion regulation and response preparation in humans as well . The CM subcompartment receives and integrates most of its projections from the LB nuclei, which maintains broad axonal connections to sensory areas. The LB has been linked to multisensory input and emotional learning [9, 38, 39], especially emotional memory . The SF subregion primarily maintains axonal projections to olfactory cortex [38, 40, 41], and it was found to be more sensitive to emotional face recognition than LB and CM , as well as to maintain the most behavioral correlates of the three amygdaloid subregions . The LB and SF comprise structurally and functionally clearly differentiable properties: while the LB mainly receives multisensory environmental input, the SF is thought to receive socially relevant information . Yet, they are often mentioned together as the olfactory/multimodal pallial section of the amygdala as both structures generally process incoming stimuli , so as to facilitate social-perceptual processing .
Our study specifically aimed to investigate global networks in ASD. Since we consider amygdalo-cortical connectivity long-range connectivity, our hypotheses were aimed at underconnectivity. Given the prominent role of the amygdala in the social brain and previous findings of long-range underconnectivity in sensory areas in ASD, we hypothesized to find reduced functional amygdalo-cortical connectivity in ASD, especially among the projections from sensory cortex to the amygdala. We applied dual-echo imaging and a stringent correction for heart rate and respiratory signals to ensure optimal sensitivity in both the amygdala and neocortical structures.
N = 19 (95 %)
N = 22 (88 %)
N = 1 (5 %)
N = 3 (12 %)
Autism Questionnaire (AQ)
Parents about participant
Autism Diagnostic Interview (ADI-R)
ADI-R A (10)
ADI-R B (8)
ADI-R C (3)
ADI-R D (1)
The study (including the informed consent procedure and all information brochures) was approved by both the regional ethics committee (Commissie Mensgebonden Onderzoek Arnhem Nijmegen) and Karakter’s review board. All participants provided verbal and written informed consent. For underage participants, parents/guardians also provided verbal and written informed consent. The signed consent forms are kept at Karakter, Child and Adolescent Psychiatry Centre Nijmegen, The Netherlands. To ensure an adequate consent procedure, (1) potential participants were provided with simple language brochures (parents/guardians were provided with regular brochures), (2) the study was not advertised with financial or other incentives other than travel reimbursement (after the scanning procedure, all participants did receive 20 euros for participating irrespective of completion), (3) participants were reminded at each phase of the study that they were free to withdraw from participating, (4) only participants with a total IQ of 80 or higher could participate, and (5) all participants practiced once with the scanning procedure in a replicate (dummy) scanner so that they could experience the scanning procedure and make an informed decision on whether or not to participate.
Image data acquisition
For each participant, we acquired magnetic resonance imaging (MRI) data at the Donders Institute for Brain, Cognition and Behaviour, Centre for Cognitive Neuroimaging in Nijmegen, The Netherlands, using a 3 Tesla Magnetom Trio (Siemens, Erlangen, Germany) with a 32-channel head coil. The entire scanning session lasted approximately 45 min. For each participant, we collected a T1-weighted whole-brain scan (magnetization-prepared rapid acquisition with gradient echo (MPRAGE), inversion time (TI) = 1100 ms, repetition time (TR) = 2300 ms, echo time (TE) = 3.03 ms, flip angle = 8°, field of view (FOV) = 256 × 256 × 192 mm3, voxel size = 1 × 1 × 1 mm3) and a resting state scan using T2*-weighted dual-echo planar imaging (EPI, TR = 2510 ms, TE1 = 16 ms, TE2 = 36 ms, flip angle = 83°, FOV = 212 × 212 × 119 mm3, voxel size = 2 × 2 × 2.5 mm3, number of volumes = 400, imaging bandwidth = 1814 Hz/px, grappa acceleration factor = 4). Note that the usage of dual-echo imaging provides optimal sensitivity for blood-oxygen-level-dependent (BOLD) imaging in both subcortical structures such as the amygdala and the neocortex . For the 16-min resting state scan, participants were instructed to lie still within the scanner with their eyes open, while staying awake and focusing on a small white cross presented at the centre of a projection screen. The first five volumes (12.55 s) were discarded to reduce magnetization equilibration effects. Gradient echo field mapping data were also acquired with identical geometry to the EPI data for EPI off-resonance distortion correction (TR = 1020 ms, TE1 = 10 ms, TE2 = 12.46 ms, flip angle = 90°, FOV = 224 × 224 × 191 mm3, voxel size = 3.5 × 3.5 × 2 mm3). All participants were able to familiarize themselves with scanner setup and scanning procedure through rehearsal in a replicate (dummy) scanner before actual image acquisition. The scanning session further also included a diffusion tensor imaging (DTI) scan (not reported here).
We recorded participants’ heartbeats using the scanner’s built-in photoplethysmograph, placed on the right index finger. Respiration was measured with a pneumatic belt positioned at the level of the abdomen. In order to reduce the potential bias that the heartbeat and respiration have in resting state BOLD correlation studies [52, 53], we used cardiac and respiratory phase regressors, as well as other nuisance regressors in the fMRI time series analysis.
All image preprocessing and analyses were performed using FSL (FMRIB Software Library, http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/) . The following pre-statistical processes were applied to the fMRI data: nonbrain removal using Brain Extraction Tool (BET); rigid-body motion correction using MCFLIRT; high-pass temporal filtering (Gaussian-weighted least-squares fitting with frequency cutoff point = 100 s); correction of off-resonance geometric distortions in the EPI data using PRELUDE and FUGUE, using B0 field maps derived from the dual-echo gradient echo dataset; artifact removal based on probabilistic ICA (Independent Component Analysis) using MELODIC; spatial normalization to Montreal Neurological Institute (MNI152) 2 mm isotropic atlas space using boundary-based registration (BBR); and FNIRT and Gaussian filtering (full width at half maximum (FWHM) = 6 mm; see the the “Statistical analysis” section). The dual-echo images (TE = 16 and TE = 36) were combined by averaging both echo times. We excluded one participant with ASD due to excessive head movement in terms of frame-wise displacement (max. frame-wise displacement (FD) = 8.7 mm, M fd = 0.89 mm), resulting in 20 datasets from the ASD group and 25 datasets from the control group for further analysis (see Additional file 1 for relative frame-wise displacement). To rule out the possibility that differences in movement between the ASD and control group could contribute to the results, we calculated the mean value of frame-wise movement (i.e., the movement of one TR relative to previous TR) for each participant and compared it between the two groups. No group difference was found (M asd = 0.10, SDasd = 0.10; M ctrl = 0.07, SDctrl = 0.42; t(25) = 1.73, p = 0.1).
Controlling for structured noise
Our preprocessing stream included several steps to limit the influence of structured noise, such as motion artifacts , heartbeat , and respiration . First, we conducted manual ICA-based artifact removal. The first author visually inspected all the independent component maps for each participant to identify noise components based on the spatial layout of the component maps and the power spectra of the associated time series . We applied nonaggressive denoising with FSL’s fsl_regfilt, i.e., only variance that was uniquely related to the components labeled as noise component (approx. 70 %) was removed.
After ICA-based noise removal and further preprocessing, we conducted nuisance regression modeling the potential effect from motion and physiological noise on the resting state fMRI data. Specifically, we included six rigid body parameters and the eigenvariate of signals over the entire white matter and the cerebrospinal fluid (CSF) in our GLM. Moreover, we calculated 10 cardiac phase regressors, 10 respiratory phase regressors, and 6 other nuisance regressors including heart rate fluctuation (HRF), heart rate variability (HRV), respiration raw data averaged per TR, respiratory amplitude in 9-s window, respiratory frequency in 9-s window, and RVT (frequency times amplitude of respiration, averaged per TR) that are derived from retrospective image correction (RETROICOR) method .
Region of interest selection
First level analyses were carried out using FSL’s seed-based correlation analysis (SBCA)  to calculate the partial correlation between the average time series of the voxels in one ROI (i.e., one of the three amygdala subregions) and the time series of every voxel of the whole brain, corrected for the average time series of the other two amygdala seeds. This tool incorporates the option to add data with different smoothing kernels for the seed and target areas. For the amygdala, cortical correlation analysis we entered spatially unfiltered data for the small amygdaloid areas and FWHM = 6 mm Gaussian filtered functional images for the whole brain. Thus, one single-subject partial correlation map of the brain for each subregion (left and right CM, LB, and SF) was obtained, yielding each subregion’s unique connectivity with the rest of the brain. In addition to the partial correlation analyses, we also performed a regular correlation analysis of the entire amygdala (EA) (left and right) with every voxel in the brain to serve as a reference for the partial correlation results, increasing the interpretability of the partial correlation results.
To test for between-group differences, we performed a nonparametric test with Randomise . As this analysis does not require the data to be normally distributed, an r to z transform is not necessary. Thus, 5000 random permutations of a threshold-free cluster enhancement statistic (TFCE)  against the null hypothesis were conducted for each ROI separately. The p values were extracted with FSL’s cluster command, where cluster peaks and local maxima with p < 0.05 were acquired from the threshold-free cluster enhancement FWE corrected 1-p statistical images of the Randomise output. The permutation method strongly controls for the family-wise error (FWE) rate when a large amount of voxels is tested. Six contrasts (positive main effects in ASD; positive main effects in controls, negative main effects in ASD; negative main effects in controls; ASD > controls; controls > ASD) were tested with an unpaired samples t test. The same approach was used for entire amygdala analysis.
Entire amygdala connectivity
Figure 2a presents the positive entire amygdala correlations in the ASD and control sample. This reveals a number of symmetric left and right frontal, occipital, temporal, and sensorimotor networks including the cingulate gyrus, which is consistent with the amygdala’s role in socio-emotive circuits . Left amygdala correlations yielded one large cluster with its peak in the left anterior parahippocampal gyrus (−28,−8,−36; p < 0.001). Local maxima extended dorsally into the left precentral gyrus, superior parietal lobe, and both sides of the postcentral gyrus and ventrally into the ipsilateral superior temporal gyrus, planum temporale, and temporal pole. Right amygdala correlations were found in the anterior and posterior cingulate gyrus (0,−2,44; p = 0.012), including the supplementary motor cortex and clusters in the right temporal fusiform cortex (32,−6,−36; p < 0.001), middle temporal gyrus, central opercular cortex, supramarginal gyrus, parietal operculum cortex, a cluster in the left superior temporal gyrus (−62,−8,−4; p = 0.002), and a number of smaller clusters in the left inferior temporal gyrus (−46,−58,−16; p = 0.036).
The results of the ASD group are comparable with the control group’s results with symmetric left and right frontal, occipital, temporal, and sensorimotor networks underlying emotion regulation circuits (see EA patients in Fig. 2a). The left amygdala revealed a large cluster in the left temporal pole (−26,4,−46; p < 0.001) with local maxima in the bilateral pre- and postcentral gyrus. The right amygdala yielded a large cluster in the right temporal pole (28,6,−42; p < 0.001) including the local maxima in the temporal areas, the brainstem, and the right hippocampus.
Negative correlations were largely absent in our sample. This is probably a consequence of (1) the participants’ age because during adolescence, negative (inhibitory) mechanisms are not yet as exuberant as after puberty’s transition phase [63–65] and (2) no global signal regression was performed, which could have introduced spurious negative correlations [66, 67] in previous work that did subtract the global mean signal.
Next, to investigate the individual contributions of the nuclei group-specific seeds, we mapped the intrinsic whole brain connectivity of each ROI (Fig. 2b–d). As before, only significant clusters surviving family-wise error correction with alpha < 0.05 are reported (Additional files 3, 4, 5, and 6).
Superficial amygdala connectivity
Left SF connectivity in controls (SF controls, yellow) showed extensive unique bilateral positive correlations throughout frontal (right frontal medial cortex (4,44,−20; p = 0.029), left frontal pole (−2,62,−6; p = 0.031), temporal (left parahippocampal gyrus (−14,−6,−26; p < 0.001), occipital, and parietal lobe and limbic areas including anterior and posterior cingulate gyrus, parahippocampal gyrus, hippocampus, amygdala, thalamus, putamen as well as in the brainstem and cerebellum (Additional files 3 and 6). The right SF showed connectivity to smaller, more selective sensory and limbic areas. In ASD (Fig. 2b, SF patients), SF connectivity yielded a similar pattern, except for a larger cluster in the frontal lobe area, with sparse overall occipital lobe connectivity.
Note that while partial correlation analysis reported for the SF yielded a relatively small frontal lobe cluster in our healthy subjects, a supplementary full correlation analysis revealed a large frontal lobe involvement, indicating that the overlapping signal between SF and CM/LB time series leads to partialling out some of the SF-frontal connectivity (Fig. 1b). Thus, the results indicate that (1) SF functional connectivity in our adolescent sample is consistent with known limbic-striatal-frontal reward related and auditory-parietal-visual valence evaluating circuits  and (2) the SF subregion maintains unique whole brain connectivity.
Laterobasal amygdala connectivity
The laterobasal subregion exhibited mainly unique connectivity with temporal regions and regions along the lateral superior cortical axis. Although left and right LB connectivity maps overlapped in the parietal operculum cortex, frontal orbital cortex, and temporal pole, strong lateralization effects were also observed (Fig. 2c, LB controls). Left LB connectivity peaks were bilaterally present in the parahippocampal gyrus (left (−26,−8,−34; p < 0.001); right (32,−28,−30; p = 0.047)), precentral gyrus (left (−44,−14,56; p = 0.014)); right (42,−14,52; p = 0.017)), and lateral occipital cortex (left (−20,−80,42; p = 0.020); right (30,−84,28; p = 0.021)) (Additional file 4). Except from one large cluster of connectivity in the right temporal fusiform cortex (32,−6,−36; p < 0.001), right LB showed many localized clusters of connectivity, such as in the right supplementary motor area (2,−8,62; p = 0.033), right frontal orbital cortex (44,28,−18; p = 0.048) and bilateral anterior cingulate gyrus (left (6,−12,44; p = 0.031); right (0,−2,44; p = 0.033)). In the ASD group, similar patterns of connectivity were observed (Fig. 2c, LB patients). Thus, the LB connectivity maps are in line with its putative involvement in emotional learning through its connection with the (para)hippocampal area and its involvement in multisensory processing via its projections to the sensory systems along the superior temporal gyrus, sensorimotor areas, and visual areas, in combination with cingulate-frontal connectivity .
Centromedial amygdala connectivity
In the control group, only a small cluster in the right occipital fusiform gyrus (28,−64,−10; p = 0.040) from the right centromedial subregion showed a negative correlation, while positive unique functional connectivity did not reach the significance threshold for the left or right CM in the control group using partial correlation analysis (Fig. 2c, CM controls). We therefore performed an additional, regular correlation analysis for the CM to test whether some of its unique signal was “partialled out” due to a functional overlap between CM and LB/SF nuclei (Fig. 1b). As expected, and in line with previous reports , positive correlations in anterior cingulate cortex were significantly higher for the CM nucleus compared to SF and LB in healthy controls. This suggests that the partial correlation approach had corrected for the strong overlap with time series from LB and SF to such an extent that the CM’s unique contribution did not reach the FWE-corrected threshold (Additional file 7). The relative lack of unique CM correlations in controls may well reflect the lower degree of functional specialization of the amygdaloid nuclei in adolescents . In the ASD group, however, left CM demonstrated unique partial correlations with bilateral primary sensorimotor areas (postcentral gyrus (left, −36,−22,40; p = 0.004), precentral gyrus (right, 40,−14,60; p = 0.007)), and left insular cortex in the partial correlation analysis, while right CM correlations with the striatum (including bilateral thalamus, left putamen (−24,0,−10; p = 0.046) and right pallidum) and right hemispheric speech processing areas (including central and frontal opercular cortex, supramarginal gyrus, and Heschl’s gyrus) were present (Fig. 2c, CM patients; Additional file 5).
Difference between ASD and controls
Entire amygdala connectivity
To test our hypothesis of reduced amygdala connectivity along sensory input channels in ASD, we directly compared the SF, LB, and CM connectivity between groups (Fig. 3b; Additional file 8). In line with our hypothesis, the ASD group showed reduced left SF connectivity with bilateral precuneus cortex (left, 16,−58,8; p = 0.003), cuneal cortex, intracalcarine cortex, lateral occipital cortex (right, 44,−60,12; p = 0.023), and supracalcarine cortex in the right hemisphere, the occipital pole, the superior parietal lobe (28,−54,62; p = 0.004), and the pre- and postcentral and angular gyri (44,−16,64; p = 0.047) when compared with controls. Furthermore, functional connectivity between the right LB and right superior parietal lobe (14, −54, 62; p = 0.023) were reduced in ASD. There were no significant between-group differences in the CM correlation maps. Statistical tests in both directions (ASD > controls; controls > ASD) were included in our analysis, but yielded nonsignificant results for ASD overconnectivity patterns. When the analysis was repeated including age as covariate, the group effects were very similar to the previous between-group results and no age-by-group interactions were found. Age showed a negative main effect with connectivity strength in the left SF, indicating that connectivity strength decreases during development in the temporal pole and lateral occipital cortex (Additional file 9).
Relation between reduced connectivity strength and AQ scores
To investigate whether the subregion-specific group differences in functional connectivity can be related to specific behavioral/symptomatic characteristics, we quantified the relationship between the participants’ scores within each group on each AQ subdomain and the (z-transformed) correlation between the mean fMRI time series from the right LB and left SF nuclei and the cortical regions that exhibited a significant group difference in functional connectivity. For either group, however, none of the specific AQ sub-domains were significantly related to the subregion-specific functional connectivity strengths. This absence of significance was also observed for the relationship between total AQ score and subregion-specific connectivity strength.
Subregion signal-to-noise ratios
Analysis of subregion signal-to-noise ratios
In the present study, we investigated amygdalo-cortical connectivity in adolescents and young adults with ASD and controls during resting state fMRI. We found reduced cortical connectivity of both amygdaloid input subregions (SF and LB) with the prefrontal, parietal, and occipital cortices in participants with ASD, while CM output connectivity was spared.
The positive correlation maps from both the entire amygdala and the three subregions in our control group of healthy adolescents revealed large overall overlap with the spontaneous activation maps previously reported in healthy adults . The present results indicate that most of the entire amygdala connectivity main effects can be disentangled into subunit functionality in adolescents and young adults with and without ASD, while some of the global effects (e.g., frontal lobe activation) could not be traced down to one particular subunit. As positive partial CM correlations did not reach FWE-corrected significance in the control group, one interpretation of these findings could entail a lower degree of functional specialization in healthy adolescents. While this may reflect the not yet fully differentiated amygdala in adolescence, the conservative partial correlation approach may also have contributed to the negative finding because normal CM anterior cingulated gyrus connections were found with the full correlation approach. The absence of significant partial correlations for the CM compartment in the controls is most likely due to commonalities (i.e., shared variance) between the CM signals and those from the other subregions, leading them to be partialled out. Indeed, a direct comparison revealed no significant differences between patients and controls. Positive laterobasal-cortical correlations in the prefrontal, parietal, and temporal cortex further confirmed previous findings that connect the LB with associative learning processes across sensory modalities [33, 71], while SF-cortical correlations were in line with known limbic lobe connectivity .
Importantly, we found reduced EA connectivity in the ASD group when compared to controls. In line with our hypothesis of amygdala underconnectivity in ASD, we found reduced left EA connectivity with portions of the occipital pole, cuneal cortex, intra- and supracalcarine cortex and lateral occipital cortex in the ASD group. Partial correlation analysis revealed that these differences were largely driven by the left SF subregion, with reduced left SF connectivity in the precuneus, cuneus, angular gyrus, precentral, and postcentral gyrus and superior parietal cortex. The ASD group also exhibited reduced right entire amygdala connectivity with the right superior parietal lobe. Partial correlation analysis revealed that this difference was mainly driven by the right LB.
Prior research showed that functional specialization of amygdaloid subregions continues throughout adolescence [70, 72] and that volumetric abnormalities in ASD are age specific [10, 11]. The possibility that our between-group effects were driven by developmental differences was however not supported by our data. There were no group-by-age interactions for EA and subregion-specific connectivity, and our reported between-group effects were similar to the between-group effects with age as covariate. The negative main effects of age showed no overlap with the abnormalities found for the left SF in our ASD group. Thus, while the age effect in the left superficial area might reflect the developmental changes of the amygdaloid subcompartments throughout adolescence as reported by Gabard-Durnam and colleagues  and Qin et al. , the amygdala-cortical abnormalities in ASD are not significantly age dependent, at least within the age range examined in the present study.
Although our amygdalo-cortical analysis was not specifically designed for investigating local overconnectivity, increased activation patterns in ASD were also tested and yielded nonsignificant results. Investigating the local overconnectivity account remains challenging in fMRI research , and only few studies particularly investigated local overconnectivity in fMRI resting state [19, 21–23] and generated inconsistent results.
Previous task-based fMRI studies have reported some support for pathway-specific deficits of the parietal visuospatial domain in ASD [28, 29, 74] and its connections to the amygdala . Since reduced functional amygdalo-cortical connectivity in our ASD sample was mainly present in the dorso-dorsal and ventro-dorsal pathway, our results suggest that abnormal amygdaloid connectivity in ASD is pathway specific. That is, the fact that we only found abnormal connectivity with the left SF and the right LB, i.e., the amygdaloid input areas, but not with the CM output areas, supports the notion that social-emotional deficits in ASD may be reflected in reduced connectivity along amygdalo-sensory input pathways. Thus, a deficiency specific to the amygdalo-cortical input pathway may account for the social perceptual deficits in ASD.
Our results also showed lateralized subregion-specific amygdaloid connectivity, which contradicts a previous finding of bilateral homogenous amygdala connectivity . However, a number of studies associated the left amygdala with slower explicit emotion appraisal processes, while the right amygdala is more involved with faster implicit threat detection . Another hemispheric lateralization account distinguishes the left hemispheric abstract category subsystem and the right hemispheric whole-based subsystem: the left subsystem uses a parts-based processing strategy to represent smaller features of larger whole objects, while the right hemispheric whole-based subsystem serves visual discrimination of similar objects (see ). As the amygdalo-cortical deficits in our ASD sample accumulated along the left SF and right LB sensory input pathways, while other amygdalo-cortical pathways were spared, the results might indicate that the cortical processing of visual object features may be affected in ASD. That is, deficits in the left SF may account for abnormal parts-based perceptual processes along the ventro-dorsal and dorso-dorsal perceptual pathway, while abnormal functional connectivity in the right LB might for instance reflect whole-object face processing difficulties in ASD caused by abnormal amygdalo-cortical connectivity with the right superior parietal lobe.
Within groups, we did not observe a significant relationship between the subregion-specific reductions in functional connectivity strength and the participants’ AQ scores. Given the observed differences in subregion-specific functional connectivity strength across groups, and the fact that ASD status and AQ score are clearly related (see Table 1), it is likely that the absence of statistical significance for these within-group comparisons is related to insufficient statistical power. Future work based on larger sample sizes may therefore be able to tease apart the effects of subregion-specific reductions in functional connectivity on behavior.
Our investigative approach draws on Roy and colleagues’ paper that describes resting state analysis of probabilistic cytoarchitectonically defined amygdala nuclei . In our study, however, partial correlation analysis was used instead of regression analysis with statistical mean testing to compute the unique contributions of each of the three amygdaloid subdivisions, and we studied adolescent brains rather than healthy adult brains. As such, small differences between the two studies may be expected. For instance, in our study, the CM subcompartment showed connectivity with the striatal circuitry in the ASD group, while positive partial correlations did not reach FWE-corrected significance in the control group. Because we did find normal positive correlations in a full correlation analysis for the CM (Fig. 1b), we interpret the lack of partial CM-whole-brain correlations in controls as a reflection of the lower degree of functional specialization of the amygdaloid nuclei in the adolescent brain . Further support for the notion of functional overlap between our three subcompartments was found during supplementary analysis of colinearity (Additional file 7). A putative faster maturation of the CM region in ASD as compared to controls was not supported by the data; we directly tested for increased connectivity patterns in ASD compared to controls but found no difference. Overall, we regard our method as a valid approach for detecting true functional connectivity differences between healthy and clinical populations , since (1) the entire amygdala results showed strong overlap with the subregion-specific outcomes from the partial correlation analysis (Additional file 11), (2) the current results clearly demonstrate the expected functional connectivity with known amygdala circuits, and (3) sensitivity to between-group effects increases with partial correlation analysis (Fig. 3, Additional file 8).
To conclude, the current findings provide further evidence for underconnectivity in socio-emotive circuits in adolescents and young adults with ASD. As we found abnormal connectivity in the amygdala’s input areas but not in the output areas, the findings support the notion that deficient/impaired amygdaloid sensory input mechanisms may underlie ASD. This might indicate that therapeutic interventions should target sensory input channels.
Adult Behavior Check List/age 18-59
Autism Diagnostic Interview-Revised
attention deficit/hyperactivity disorder
analysis of variance
autism spectrum disorder
autism spectrum quotient
Brain Extraction Tool
symbol for the constant magnetic field
Child Behavior Check List/age 6-18
centromedial amygdaloid subcompartment
Diagnostic and Statistical Manual of Mental Disorders (4th edition)
diffusion tensor imaging
functional magnetic resonance imaging
(Oxford centre for) functional magnetic resonance imaging of the brain
FMRIB’s nonlinear image registration tool
field of view, spatial encoding area of the image
FMRIB Software Library
FMRIB's Utility for Geometrically Unwarping EPIs
full width at half maximum
heart rate fluctuation
heart rate variability
independent component analysis
laterobasal amygdaloid subcompartment
FMRIB's Linear registration and motion correction tool
Multivariate Exploratory Linear Optimized Decomposition into Independent Components
magnetization-prepared rapid acquisition with gradient echo
magnetic resonance imaging
Phase Region Expanding Labeller for Unwrapping Discrete Estimates
retrospective image correction
region of interest
respiration volume per time
seed-based correlation analysis
superficial amygdaloid subcompartment
- SPSS 20:
Statistical Package for the Social Sciences (20th edition)
an image whose contrast and brightness are predominately determined by T1 signals (T1: time constant at which spins realign themselves with the external magnetic field after excitation)
T-two-star, time constant for the loss of phase coherence among spins oriented at an angle to the static magnetic field
Wechsler Adult Intelligence Scale (3rd edition)
Wechsler Intelligence Scale for Children (3rd edition)
This work was supported by the Hersenstichting Fellowship (F2010(1)-13). The authors would like to thank Marieke W.M. Kuiper for her assistance with the data collection.
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.
- American Psychiatric Association. Diagnostic and statistical manual of mental disorders: DSM-IV-TR®. Washington (DC): American Psychiatric Publisher; 2000.Google Scholar
- American Psychiatric Association. Diagnostic and statistical manual of mental disorders. 5th ed. Washington (DC): American Psychiatric Publisher; 2013.Google Scholar
- Verhoeven JS, De Cock P, Lagae L, Sunaert S. Neuroimaging of autism. Neuroradiology. 2010;52(1):3–14.PubMedView ArticleGoogle Scholar
- Just MA, Cherkassky VL, Keller TA, Kana RK, Minshew NJ. Functional and anatomical cortical underconnectivity in autism: evidence from an FMRI study of an executive function task and corpus callosum morphometry. Cereb Cortex. 2007;17(4):951–61.PubMedPubMed CentralView ArticleGoogle Scholar
- Anagnostou E, Taylor MJ. Review of neuroimaging in autism spectrum disorders: what have we learned and where we go from here. Mol Autism. 2011;2:4.PubMedPubMed CentralView ArticleGoogle Scholar
- Amaral DG, Schumann CM, Nordahl CW. Neuroanatomy of autism. Trends Neurosci. 2008;31(3):137–45.PubMedView ArticleGoogle Scholar
- Ashwin C, Chapman E, Colle L, Baron-Cohen S. Impaired recognition of negative basic emotions in autism: a test of the amygdala theory. Soc Neurosci. 2006;1(3-4):349–63.PubMedView ArticleGoogle Scholar
- Baron-Cohen S, Ring HA, Bullmore ET, Wheelwright S, Ashwin C, Williams SC. The amygdala theory of autism. Neurosci Biobehav Rev. 2000;24(3):355–64.PubMedView ArticleGoogle Scholar
- Phelps EA, LeDoux JE. Contributions of the amygdala to emotion processing: from animal models to human behavior. Neuron. 2005;48(2):175.PubMedView ArticleGoogle Scholar
- Bellani M, Calderoni S, Muratori F, Brambilla P. Brain anatomy of autism spectrum disorders II. Focus on amygdala. Epidemiol Psychiatr Sci. 2013;22:309–12.PubMedView ArticleGoogle Scholar
- Groen W, Teluij M, Buitelaar J, Tendolkar I. Amygdala and hippocampus enlargement during adolescence in autism. J Am Acad Child Adolesc Psychiatry. 2010;49(6):552–60.PubMedGoogle Scholar
- Green SA, Rudie JD, Colich NL, Wood JJ, Shirinyan D, Hernandez L, et al. Overreactive brain responses to sensory stimuli in youth with autism spectrum disorders. J Am Acad Child Adolesc Psychiatry. 2013;52(11):1158–72.PubMedView ArticleGoogle Scholar
- Tottenham N, Hertzig ME, Gillespie-Lynch K, Gilhooly T, Millner AJ, Casey B. Elevated amygdala response to faces and gaze aversion in autism spectrum disorder. Soc Cogn Affect Neurosci. 2014;9(1):106–17.PubMedPubMed CentralView ArticleGoogle Scholar
- Nacewicz BM, Dalton KM, Johnstone T, Long MT, McAuliff EM, Oakes TR, et al. Amygdala volume and nonverbal social impairment in adolescent and adult males with autism. Arch Gen Psychiatry. 2006;63(12):1417–28.PubMedView ArticleGoogle Scholar
- Schumann CM, Barnes CC, Lord C, Courchesne E. Amygdala enlargement in toddlers with autism related to severity of social and communication impairments. Biol Psychiatry. 2009;66(10):942–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Courchesne E, Pierce K. Why the frontal cortex in autism might be talking only to itself: local over-connectivity but long-distance disconnection. Curr Opin Neurobiol. 2005;15(2):225–30.PubMedView ArticleGoogle Scholar
- Anderson JS, Druzgal TJ, Froehlich A, DuBray MB, Lange N, Alexander AL, et al. Decreased interhemispheric functional connectivity in autism. Cereb Cortex. 2011;21(5):1134–46.PubMedPubMed CentralView ArticleGoogle Scholar
- Travers BG, Adluru N, Ennis C, Tromp DP, Destiche D, Doran S, et al. Diffusion tensor imaging in autism spectrum disorder: a review. Autism Res. 2012;5(5):289–313.PubMedPubMed CentralView ArticleGoogle Scholar
- Di Martino A, Zuo X-N, Kelly C, Grzadzinski R, Mennes M, Schvarcz A, et al. Shared and distinct intrinsic functional network centrality in autism and attention-deficit/hyperactivity disorder. Biol Psychiatry. 2013;74(8):623–32.PubMedPubMed CentralView ArticleGoogle Scholar
- Khan S, Gramfort A, Shetty NR, Kitzbichler MG, Ganesan S, Moran JM, et al. Local and long-range functional connectivity is reduced in concert in autism spectrum disorders. Proc Natl Acad Sci U S A. 2013;110(8):3107–12.PubMedPubMed CentralView ArticleGoogle Scholar
- Maximo JO, Keown CL, Nair A, Müller R-A. Approaches to local connectivity in autism using resting state functional connectivity MRI. Front Hum Neurosci. 2013;7:605.PubMedPubMed CentralView ArticleGoogle Scholar
- Paakki J-J, Rahko J, Long X, Moilanen I, Tervonen O, Nikkinen J, et al. Alterations in regional homogeneity of resting-state brain activity in autism spectrum disorders. Brain Res. 2010;1321:169–79.PubMedView ArticleGoogle Scholar
- Shukla DK, Keehn B, Müller RA. Regional homogeneity of fMRI time series in autism spectrum disorders. Neurosci Lett. 2010;476(1):46–51.PubMedPubMed CentralView ArticleGoogle Scholar
- Mottron L, Dawson M, Soulieres I, Hubert B, Burack J. Enhanced perceptual functioning in autism: an update, and eight principles of autistic perception. J Autism Dev Disord. 2006;36(1):27–43.PubMedView ArticleGoogle Scholar
- Samson F, Mottron L, Jemel B, Belin P, Ciocca V. Can spectro-temporal complexity explain the autistic pattern of performance on auditory tasks? J Autism Dev Disord. 2006;36(1):65–76.PubMedView ArticleGoogle Scholar
- Groen WB, van Orsouw L, Ter Huurne N, Swinkels S, van der Gaag R-J, Buitelaar JK, et al. Intact spectral but abnormal temporal processing of auditory stimuli in autism. J Autism Dev Disord. 2009;39(5):742–50.PubMedView ArticleGoogle Scholar
- Foxton JM, Stewart ME, Barnard L, Rodgers J, Young AH, O’Brien G, et al. Absence of auditory ‘global interference’in autism. Brain. 2003;126(12):2703–9.PubMedView ArticleGoogle Scholar
- Pellicano E, Gibson L, Maybery M, Durkin K, Badcock DR. Abnormal global processing along the dorsal visual pathway in autism: a possible mechanism for weak visuospatial coherence? Neuropsychologia. 2005;43(7):1044–53.PubMedView ArticleGoogle Scholar
- Braddick O, Atkinson J, Wattam-Bell J. Normal and anomalous development of visual motion processing: motion coherence and ‘dorsal-stream vulnerability’. Neuropsychologia. 2003;41(13):1769–84.PubMedView ArticleGoogle Scholar
- Bertone A, Mottron L, Jelenic P, Faubert J. Enhanced and diminished visuo-spatial information processing in autism depends on stimulus complexity. Brain. 2005;128(10):2430–41.PubMedView ArticleGoogle Scholar
- Amunts K, Kedo O, Kindler M, Pieperhoff P, Mohlberg H, Shah N, et al. Cytoarchitectonic mapping of the human amygdala, hippocampal region and entorhinal cortex: intersubject variability and probability maps. Anat Embryol (Berl). 2005;210(5-6):343–52.View ArticleGoogle Scholar
- Ball T, Rahm B, Eickhoff SB, Schulze-Bonhage A, Speck O, Mutschler I. Response properties of human amygdala subregions: evidence based on functional MRI combined with probabilistic anatomical maps. PLoS One. 2007;2(3), e307.PubMedPubMed CentralView ArticleGoogle Scholar
- Roy AK, Shehzad Z, Margulies DS, Kelly AC, Uddin LQ, Gotimer K, et al. Functional connectivity of the human amygdala using resting state fMRI. Neuroimage. 2009;45(2):614.PubMedPubMed CentralView ArticleGoogle Scholar
- Pitkänen A, Savander V, LeDoux JE. Organization of intra-amygdaloid circuitries in the rat: an emerging framework for understanding functions of the amygdala. Trends Neurosci. 1997;20(11):517–23.PubMedView ArticleGoogle Scholar
- Pessoa L. Emotion and cognition and the amygdala: from “what is it?” to “what’s to be done?”. Neuropsychologia. 2010;48(12):3416–29.PubMedPubMed CentralView ArticleGoogle Scholar
- Moreno N, González A. Evolution of the amygdaloid complex in vertebrates, with special reference to the anamnio‐amniotic transition. J Anat. 2007;211(2):151–63.PubMedPubMed CentralView ArticleGoogle Scholar
- Bickart KC, Hollenbeck MC, Barrett LF, Dickerson BC. Intrinsic amygdala–cortical functional connectivity predicts social network size in humans. J Neurosci. 2012;32(42):14729–41.PubMedPubMed CentralView ArticleGoogle Scholar
- Bzdok D, Laird AR, Zilles K, Fox PT, Eickhoff SB. An investigation of the structural, connectional, and functional subspecialization in the human amygdala. Hum Brain Mapp. 2013;34(12):3247–66.PubMedView ArticleGoogle Scholar
- LeDoux J. The emotional brain, fear, and the amygdala. Cell Mol Neurobiol. 2003;23(4-5):727–38.PubMedView ArticleGoogle Scholar
- Heimer L, Van Hoesen GW. The limbic lobe and its output channels: implications for emotional functions and adaptive behavior. Neurosci Biobehav Rev. 2006;30(2):126–47.PubMedView ArticleGoogle Scholar
- Price JL. Comparative aspects of amygdala connectivity. Ann N Y Acad Sci. 2003;985(1):50–8.PubMedView ArticleGoogle Scholar
- Goossens L, Kukolja J, Onur OA, Fink GR, Maier W, Griez E, et al. Selective processing of social stimuli in the superficial amygdala. Hum Brain Mapp. 2009;30(10):3332–8.PubMedView ArticleGoogle Scholar
- Pessoa L, Adolphs R. Emotion processing and the amygdala: from a ‘low road’ to ‘many roads’ of evaluating biological significance. Nat Rev Neurosci. 2010;11(11):773–83.PubMedPubMed CentralView ArticleGoogle Scholar
- Adolphs R. The neurobiology of social cognition. Curr Opin Neurobiol. 2001;11(2):231–9.PubMedView ArticleGoogle 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(5):659–85.PubMedView ArticleGoogle Scholar
- Kort W, Schittekatte M, Dekker P, Verhaeghe P, Compaan E, Bosmans M. Wechsler intelligence scale for children-dutch version. London: Psychological Corporation; 2005.Google Scholar
- Wechsler D. WAIS-III Nederlandstalige bewerking: Afname-en scoringshandleiding [WAIS-III Dutch version: User manual]. Lisse, The Netherlands: Swets and Zeitlinger; 2000.Google Scholar
- Oldfield RC. The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia. 1971;9(1):97–113.PubMedView ArticleGoogle Scholar
- Baron-Cohen S, Wheelwright S, Skinner R, Martin J, Clubley E. The Autism-Spectrum Quotient (AQ): evidence from Asperger Syndrome/high-functioning autism, males and females, scientists and mathematicians. J Autism Dev Disord. 2001;31(1):5–17.PubMedView ArticleGoogle Scholar
- Woodbury-Smith MR, Robinson J, Wheelwright S, Baron-Cohen S. Screening adults for Asperger Syndrome using the AQ: a preliminary study of its diagnostic validity in clinical practice. J Autism Dev Disord. 2005;35(3):331–5.PubMedView ArticleGoogle Scholar
- Poser BA, Versluis MJ, Hoogduin JM, Norris DG. BOLD contrast sensitivity enhancement and artifact reduction with multiecho EPI: parallel‐acquired inhomogeneity‐desensitized fMRI. Magn Reson Med. 2006;55(6):1227–35.PubMedView ArticleGoogle Scholar
- Chang C, Metzger CD, Glover GH, Duyn JH, Heinze H-J, Walter M. Association between heart rate variability and fluctuations in resting-state functional connectivity. Neuroimage. 2013;68:93–104.PubMedPubMed CentralView ArticleGoogle Scholar
- Birn RM, Smith MA, Jones TB, Bandettini PA. The respiration response function: the temporal dynamics of fMRI signal fluctuations related to changes in respiration. Neuroimage. 2008;40(2):644–54.PubMedPubMed CentralView ArticleGoogle Scholar
- Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, Johansen-Berg H, et al. Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage. 2004;23:S208–S19.PubMedView ArticleGoogle Scholar
- Power JD, Barnes KA, Snyder AZ, Schlaggar BL, Petersen SE. Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage. 2012;59(3):2142–54.PubMedPubMed CentralView ArticleGoogle Scholar
- Kelly Jr RE, Alexopoulos GS, Wang Z, Gunning FM, Murphy CF, Morimoto SS, et al. Visual inspection of independent components: defining a procedure for artifact removal from fMRI data. J Neurosci Methods. 2010;189(2):233–45.PubMedPubMed CentralView ArticleGoogle Scholar
- Glover GH, Li TQ, Ress D. Image‐based method for retrospective correction of physiological motion effects in fMRI: RETROICOR. Magn Reson Med. 2000;44(1):162–7.PubMedView ArticleGoogle Scholar
- O’Reilly JX, Beckmann CF, Tomassini V, Ramnani N, Johansen-Berg H. Distinct and overlapping functional zones in the cerebellum defined by resting state functional connectivity. Cereb Cortex. 2010;20(4):953–65.PubMedPubMed CentralView ArticleGoogle Scholar
- Nichols TE, Holmes AP. Nonparametric permutation tests for functional neuroimaging: a primer with examples. Hum Brain Mapp. 2002;15(1):1–25.PubMedView ArticleGoogle Scholar
- Smith SM, Nichols TE. Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. Neuroimage. 2009;44(1):83–98.PubMedView ArticleGoogle Scholar
- IBM Corp. IBM SPSS Statistics for Windows, Version 20.0. Armonk: IBM Corp; 2011.Google Scholar
- Phillips ML, Drevets WC, Rauch SL, Lane R. Neurobiology of emotion perception I: the neural basis of normal emotion perception. Biol Psychiatry. 2003;54(5):504–14.PubMedView ArticleGoogle Scholar
- Durston S, Thomas KM, Yang Y, Uluğ AM, Zimmerman RD, Casey B. A neural basis for the development of inhibitory control. Dev Sci. 2002;5(4):F9–F16.View ArticleGoogle Scholar
- Paus T. Mapping brain maturation and cognitive development during adolescence. Trends Cogn Sci. 2005;9(2):60–8.PubMedView ArticleGoogle Scholar
- Williams BR, Ponesse JS, Schachar RJ, Logan GD, Tannock R. Development of inhibitory control across the life span. Dev Psychol. 1999;35(1):205.PubMedView ArticleGoogle Scholar
- Murphy K, Birn RM, Handwerker DA, Jones TB, Bandettini PA. The impact of global signal regression on resting state correlations: are anti-correlated networks introduced? Neuroimage. 2009;44(3):893–905.PubMedPubMed CentralView ArticleGoogle Scholar
- Cole DM, Smith SM, Beckmann CF. Advances and pitfalls in the analysis and interpretation of resting-state FMRI data. Front Syst Neurosci. 2010;4:8.PubMedPubMed CentralGoogle Scholar
- Koelsch S, Skouras S, Fritz T, Herrera P, Bonhage C, Küssner MB, et al. The roles of superficial amygdala and auditory cortex in music-evoked fear and joy. Neuroimage. 2013;81:49–60.PubMedView ArticleGoogle Scholar
- Sah P, Faber E, De Armentia ML, Power J. The amygdaloid complex: anatomy and physiology. Physiol Rev. 2003;83(3):803–34.PubMedView ArticleGoogle Scholar
- Qin S, Young CB, Supekar K, Uddin LQ, Menon V. Immature integration and segregation of emotion-related brain circuitry in young children. Proc Natl Acad Sci U S A. 2012;109(20):7941–6.PubMedPubMed CentralView ArticleGoogle Scholar
- Schoenbaum G, Chiba AA, Gallagher M. Changes in functional connectivity in orbitofrontal cortex and basolateral amygdala during learning and reversal training. J Neurosci. 2000;20(13):5179–89.PubMedGoogle Scholar
- Gabard-Durnam LJ, Flannery J, Goff B, Gee DG, Humphreys KL, Telzer E, et al. The development of human amygdala functional connectivity at rest from 4 to 23 years: a cross-sectional study. Neuroimage. 2014;95:193–207.PubMedPubMed CentralView ArticleGoogle Scholar
- Vissers ME, Cohen MX, Geurts HM. Brain connectivity and high functioning autism: a promising path of research that needs refined models, methodological convergence, and stronger behavioral links. Neurosci Biobehav Rev. 2012;36(1):604–25.PubMedView ArticleGoogle Scholar
- Bertone A, Mottron L, Jelenic P, Faubert J. Motion perception in autism: a “complex” issue. J Cogn Neurosci. 2003;15(2):218–25.PubMedView ArticleGoogle Scholar
- Bonda E, Petrides M, Ostry D, Evans A. Specific involvement of human parietal systems and the amygdala in the perception of biological motion. J Neurosci. 1996;16(11):3737–44.PubMedGoogle Scholar
- Gläscher J, Adolphs R. Processing of the arousal of subliminal and supraliminal emotional stimuli by the human amygdala. J Neurosci. 2003;23(32):10274–82.PubMedGoogle Scholar
- McMenamin BW, Marsolek CJ. Can theories of visual representation help to explain asymmetries in amygdala function? Cogn Affect Behav Neurosci. 2013;13(2):211–24.PubMedPubMed CentralView ArticleGoogle Scholar
- Smith SM, Miller KL, Salimi-Khorshidi G, Webster M, Beckmann CF, Nichols TE, et al. Network modelling methods for FMRI. Neuroimage. 2011;54(2):875–91.PubMedView ArticleGoogle Scholar