2.1 Participants
Participants were recruited as part of a larger study examining OXT and cortisol levels in autistic and non-autistic adults [35]. Recruitment sources included the “Outpatient and Day Clinic for Disorders of Social Interaction” at the Max Planck Institute of Psychiatry (MPIP) for autistic individuals and an online study application system on the Institute’s website, as well as public advertisements for participants in the comparison group. Of the sixty-four participants included in the larger project, fifty-nine underwent structural MRI. After exclusion of three scans following a quality check protocol (2 due to poor image quality, 1 due to structural abnormalities, see Appendix), the final dataset included brain scans from fifty-six adults aged 18–60 years: twenty-nine autistic adults in the ASD group (17 men; mean age = 36.03 ± 11.0 years) and twenty-seven non-autistic adults in the CG (comparison group) (9 men; mean age = 30.96 ± 11.2 years). Autistic individuals met DSM-V criteria for ASD and had been previously diagnosed in accordance with current guidelines [36]. Autistic participants had no intellectual impairment (IQ scores > 70) and were thus regarded as individuals with high-functioning autism (HFA). A detailed description of the diagnostic assessment is reported in [35]. Test for normal distribution of demographic data using a Shapiro-Wilk test revealed a non-normal distribution for age, handedness scores and verbal IQ (WST [37]). Group comparison for these variables was, therefore, performed with a Man-Whitney U test, whereas Chi2 tests were used for categorical variables and t-tests for continuous variables. There were no significant differences (all p > 0.05) between groups in age, sex distribution, handedness score, BMI, verbal IQ, or lifestyle factors (smoking, alcohol, exercise) (Tbl. 1). Thirteen subjects in the ASD group took psychiatric medication regularly. These patients were asked not to take the medication in the morning before the OXT measurements, but after the experiments. General exclusion criteria were severe somatic illness, a current or previous schizophrenia diagnosis, breastfeeding, pregnancy, hormonal contraception, and a contraindication to MRI. The study protocol followed the guidelines of the Declaration of Helsinki and was approved by the ethics committee of the Ludwig-Maximilians-University of Munich. All participants gave written informed consent before participating in the study and received fixed monetary compensation at the end of the experiment.
Table 1 Demographic data.
|
CG (n = 27)
|
ASD (n = 29)
|
group comparison
|
|
mean (± SD) / N (%)
|
mean (± SD) / N (%)
|
t/χ2/U-test
|
p
|
Age
|
30.96 (11.2)
|
36.03 (11.0)
|
275.5
|
0.06
|
Sex (male:female)
|
9:18
|
17:12
|
3.60
|
0.06
|
BMI b
|
22.61 (4.1)
|
24.2 (4.3)
|
-1.39
|
0.17
|
Handedness score b
|
89.62
|
66.10
|
281.50
|
0.13
|
Smoker
|
2 (7.40%)
|
5 (17.24%)
|
1.24
|
0.27
|
Alcohol
|
10 (37.03%)
|
12 (41.37%)
|
0.11
|
0.74
|
Exercise
|
20 (74.07%)
|
17 (58.62%)
|
1.49
|
0.22
|
Psychiatric medication
|
0 (0%)
|
13 (44.82%)
|
15.76
|
< 0.001
|
AQ-scores a
|
13.89 (5.47)
|
35.07 (10.27)
|
-9.6
|
< 10− 11
|
a Data available for 55 participants .
b Data available for 54 participants.
2.2 Experimental procedures and analyses
2.2.1 Questionnaires
For the quantification of autistic traits, the Autism Spectrum Quotient AQ [38] was used. The AQ is a well-established, self-report questionnaire that provides a scaled measure of the characteristics associated with ASD on a scale of 0 to 50. The autistic traits themselves are regarded as a dimensional construct, which reflects both the autistic and the non-autistic population [39, 40]. Groups were compared using an independent samples t-test in SPSS (IBM Corp. Released 2020. IBM SPSS Statistics for Windows, Version 27.0. Armonk, NY: IBM Corp). AQ scores were available for n = 55 participants (ASD: n = 28; CG: n = 27).
2.2.2 Oxytocin measurements
For a detailed description of sample acquisition and OXT extraction please refer to the relevant publication [35]. In brief, the participants were asked to abstain from food (> 12h), water (> 1h) and sports the day before the study. After arriving at the outpatient unit of the MPIP at 8:30 am, blood samples were obtained at rest. OXT concentrations were quantified in an external laboratory (RIAgnosis, Sinzing, Germany) using radioimmunoassay (RIA) as previously described [41]. For VBM correlation analyses, which included peripheral OXT baseline levels, the concentrations in pg/ml were used. To confirm that the present subsample of participants showed OXT characteristics similar to the sample included in Albantakis et al. 2021 [35], univariate analyses were performed to test for effects of diagnostic group on baseline OXT concentrations, including age and sex as covariates. Data on baseline OXT levels were available for n = 53 participants (ASD: n = 26; CG: n = 27).
2.2.3 Image Acquisition
Structural brain scans were obtained using a 3 Tesla MRI scanner by GE, model ‘Discovery MR750’. In the scanner room participants received earplugs and an emergency bell in their right hand, the head was then secured firmly in the head coil with inflatable cushions. 3D anatomical T1-weighted MPRAGE (magnetization-prepared rapid gradient echo) magnetic resonance images were acquired with the following parameters: Flip angle 12°, Prep Time 450, echo time (TE) 2.3 ms, frequency 256, voxel resolution 1×1×1 mm.
2.2.4 Image processing
First, T1 images were converted from DICOM into NIfTI format. After checking for correct alignment to the AC/PC axis, images were pre-processed using the publicly available toolbox Cat12 Version 12.7 (1600) (http://www.neuro.uni-jena.de/cat/) implemented in SPM12 (Statistical Parametric Mapping software, http://www.fil.ion.ucl.ac.uk/spm/) using Matlab version R2019a (The MathWorks, Inc., Natick, Massachusetts, USA). Pre-processing was carried out using the standard pipeline and pre-set parameters as suggested in the Cat12 manual (http://www.neuro.uni-jena.de/cat12/CAT12-Manual.pdf). Pre-processing with Cat12 involved bias field inhomogeneity correction and denoising, using the Spatially Adaptive Non-Local Means (SANLM) Filter [42], segmentation into grey matter (GM), white matter (WM) and cerebrospinal fluid (CSF) in accordance with the unified registration approach [43] and spatial normalization and affine registration to MNI space using a template for high-dimensional DARTEL registration derived from 555 healthy subjects of the IXI- database (http://brain-development.org/) with a final voxel size of 1.5 × 1.5 × 1.5 mm. In version 12.7 of CAT12 used here, this process is extended by refined voxel-based processing using adaptive maximum a posteriori (AMAP) estimation, a Markov Random Field approach (MRF) [44] and accounting for partial volume effects [45]. Finally, segmentations were modulated by multiplication with the Jacobian determinant derived from spatial registrations. This step preserves the original volumes within a voxel, which are altered during registration [46] and is recommended by default. For detailed description of the individual steps in Cat12 we refer to the publishers website. Prior to smoothing, images were checked for correct pre-processing in accordance with the quality check protocol suggested in the Cat12 manual (see Appendix). Following suggestions of applying comparably small kernels for analyses in the HTH due to its small size and size of expected effects [30, 47] images were smoothed with a Gaussian kernel of 4 mm (FWHM). An absolute grey matter threshold masking of 0.1 was applied to account for a possible misclassification of tissues.
2.2.5 ROI-based VBM analysis
For statistical analysis, the general linear model (GLM) approach as implemented in SPM12 was used. Smoothed and modulated GM images were entered in a region of interest (ROI)- based VBM analysis using a HTH mask derived from the subcortical brain nuclei atlas by Pauli et al. (2018) (https://neurovault.org/collections/3145/). The mask was resliced to fit the template space in SPM12 and encompassed 1085 voxels. Clusters were regarded as significant when falling below an initial uncorrected voxel threshold of .001 and an FWE-corrected cluster threshold of .05 inside the ROI. To account for variance due to differences in total intracranial volume (TIV), age and sex, these variables were included as nuisance parameters. Due to co-linearity between TIV and OXT (cos (θ) = r =- 0.35), the correction for TIV in the analysis including OXT levels was applied by global scaling with TIV, as is recommended [49]. Voxel wise group comparison of hypothalamic GMV was conducted using a two-sample t-test for the contrasts ASD > CG and ASD < CG. To test for group differences in associations of GMV and OXT, (i.e. interaction effects), OXT concentrations were included in a full-factorial model and tested for significant positive and negative contrasts inside the HTH mask, i.e. (GMV[ASD] × OXT > GMV[CG] × OXT) and the other way round. To test if hypothalamic GMV is associated with autistic traits, a multiple regression analysis on all subjects was performed with AQ scores as covariate of interest and tested for significant positive and negative associations. Since we did not find a significant association here, we subsequently tested associations of GMV and AQ scores in both groups separately. To investigate whether significant results at the voxel level were also evident at the level of the overall mean hypothalamic GMV, significant statistical models were repeated using the mean hypothalamic GMV. To this end, the unadjusted eigenvariate within the HTH mask was extracted in SPM and implemented in SPSS for further analysis.