Sample
Participants were part of the QUEST follow-up study (25), a longitudinal population-based sample recruited at age 4–8 years (Wave 1; N = 277) and followed up at ages 11–15 years (Wave 2; N = 211) and 13–17 years (Wave 3; N = 214), as part of the wider IAMHealth project. The original target population for the study was all children born in a four-year period, living in two London boroughs, who had a clinical diagnosis of ASD. 277 children were recruited into the study upon entry and selectively stratified into an ‘intensively studied’ (hereafter intensive; n = 101) and ‘extensively studied’ group (hereafter extensive; n = 176) and this sampling structure was maintained at subsequent waves of data collection. The current study focuses on the Wave 3 intensive group only. Although all participants had a clinical diagnosis of autism spectrum disorder, the intensive group had their diagnosis confirmed at Wave 2 with the Autism Diagnostic Observation Schedule-2 (ADOS-2; 26), and a subset also with the Autism Diagnostic Interview-Revised (ADI-R; 27). All participating families gave their written informed consent (from young people themselves if ≥ 16 years in age and were deemed to have capacity, otherwise from parents or caregivers) and the study was approved by Camden and King’s Cross Ethics Sub-Committee (17/LO/2098 for Wave 2, 17/LO/0397 for Wave 3). Table 1 gives a comparison of key measures between the full Wave 3 intensive sample (n = 77) versus the Wave 3 intensive subsample who completed the experimental frustration task (n = 52).
Table 1
Wave 3 Full Sample and Subsample Demographic Information.
Mean (standard deviation; range)
|
Full Intensive Sample (n = 77)
|
Sample Who Completed Frustration Task (n = 52)
|
t-test of group differences
|
Age
|
15.38
(1.16; 13.2–17.8)
|
15.40
(1.10; 13.2–17.3)
|
p = .61
|
% male (n)
|
60% (46)
|
63% (33)
|
p = .61
|
IQ^
|
69.88
(31.36; 19–129)
|
84.54
(21.74; 33–129)
|
p < .001
|
Autism Severity^ (ADOS-CSS)
|
6.72
(2.66; 1–10)
|
6.31
(2.83; 1–10)
|
p = .32
|
ARI Total
|
3.74
(3.26; 0–12)
|
4.02
(3.37; 0–12)
|
p = .56
|
^measured at Wave 2, approximately two years previously
ADOS-CSS indicates Autism Diagnostic Observation Schedule-2 calibrated severity score; ARI Affective Reactivity Index
|
Measures
Psychiatric Symptoms
Clinical Interview
The Child and Adolescent Psychiatric Assessment-parent version (CAPA; 28, 29) is an interviewer-based structured diagnostic interview for use with children aged 9 -17 years. This was used to identify symptoms of psychiatric disorders that had been present in the past three months. The current study uses the total count of oppositional defiant disorder (ODD) symptoms (aside from the ‘spiteful/vindictive behaviour’ and ‘blames others’ items as these had <5 endorsements across the whole sample).
Parent-Rated Questionnaires
Affective Reactivity Index (ARI)
The ARI (30) was used to assess participants’ level of irritability and includes six items relating to feelings/behaviours specific for irritability and one question assessing impairment due to irritability, with a higher score indicative of a higher level of irritability. The internal consistency was examined in the full QUEST sample (Wave 2 intensive + extensive; n=201) and found to be excellent (α = 0.90), and comparable to that reported previously in samples of autistic young people (α = 0.82) (19).
Aberrant Behavior Checklist (ABC) – Irritability subscale
The ABC (31) is a measure developed to assess behaviour problems in children with developmental and intellectual disabilities. The 15-item Irritability subscale used currently is often used as an outcome measure in clinical trials (32). The internal consistency of the subscale was examined in the current intensive sample and found to be excellent (α = 0.93).
Direct Assessments
Baseline HR
Prior to beginning the task battery, participants watched relaxing videos for five minutes to obtain an estimate of their baseline HR. Average HR was calculated across the four consecutive 30 second segments of data collected during this period, the first 60 seconds and the last 120 seconds were excluded to ensure data quality.
Frustration Task
A novel task was designed and programmed in E-Prime 2.0 (33), based on a previously described delay frustration task (34, 35). The task was simplified to allow maximum participation in our sample. Participants were asked to select the smallest square from a choice of three. To motivate participation, participants were informed that most people their age completed around 60 trials, and a pie was shown for each trial to indicate how much time they left (see Figure 1 for a schematic of the task). However, during the task participants experienced several unexpected delays (14 delay trials out of a total of 50 trials), where the computer became unresponsive to their button presses for six seconds. These were pseudo-randomly presented, in that the first six trials were always non-delay trials, and the order of presentation was the same for each participant. The number of button presses was recorded during each six-second delay trial. The task lasted approximately five minutes and was part of a wider task battery.
Physiological Data Extraction and Processing
Electrocardiogram (ECG) data were recorded at 2000Hz using BIOPAC systems MP160 with BioNomadix wireless transmitters. Data was collected and processed using AcqKnowledge 5.0.1 (36). ECG measurements were collected using electrodes placed in the lead-II position on the back. The ECG signal was filtered using a Comb Band 50Hz filter to remove electrical noise and a 1Hz High Pass filter to remove baseline drift and movement artefact. R wave peaks, each representing a heartbeat, were automatically identified and labelled using the AcqKnowledge find cycle protocol. The signal was visually inspected to ensure that R wave peaks had been correctly identified and any movement artefact removed. HR was extracted for each inter-beat interval during the six-second delay periods. Digital markers indicating the beginning and end of each delay trial were sent via E-Prime, and these were used to demarcate the segments of data extraction.
For both the baseline and experimental task recording, four participants from the 52 who completed the direct assessments had no usable ECG data due to electrode refusal (n=2) and corrupted data files (n=2). For the baseline recording, ECG segments with more than three consecutive missing peaks or 10% of data missing were excluded (as in previous studies of autistic populations) (37), and participants with ≥ 50% missing task data overall were excluded (n=4), leaving a final sample of n=44 for the baseline HR data. For the experimental task, ECG task segments with more than one peak missing were excluded, and as before, participants with ≥ 50% missing task data were excluded (n=1), leaving a final sample of n=47 for the frustration task HR data.
Statistical Analysis
All analyses were conducted in Stata 16. First, bivariate correlations were run between the ARI and demographic characteristics (age, sex, IQ and autism severity), and other measures of irritability (parent-rated ABC irritability subscale and the number of ODD symptoms from the CAPA), to confirm the construct validity of the ARI. Next, multilevel mixed-effect regression was used to test associations between ARI and trajectories of behavioural and physiological responses during the frustration task. The key term of interest was the time-by-irritability interaction, but we also tested for main effects of irritability, equating to an association with the overall number of presses/HR (rather than the slope of change). Any significant interactions were explored graphically using Empirical Bayes’ estimates of behavioural and physiological response and low vs. high irritability, defined using a median split on the ARI. Age was included as a covariate in all task analyses, along with baseline physiology in HR analyses. As the behavioural data was the count of presses during each delay trial, a negative binomial model was specified. Likelihood ratio (LR) tests suggested a model with random intercept and slope was adequate for HR data, but the addition of a quadratic term of time (time2) was necessary for the press data (LR χ2(1) = 18.21, p<.01). After primary analyses, the following variables were added as covariates; IQ, autism severity, and medication status (coded as a binary variable of currently taking medication yes/no; split 45/55% (n=23/28); made up of 13% (n=3) minor tranquilizers/sedatives, 22% (n=5) stimulants, 4% (n=1) non-stimulants (e.g, atomoxetine, guanfacine, clonidine), 17% (n=4) anti-depressant, 9% (n=2) anti-convulsant, 22% (n=5) asthma medication, 52% (n=12) other medication) to assess the evidence for potential confounders (especially those which may index difficulties in understanding the task) on motor and physiological response. We report unstandardized coefficients throughout (b).