Participants
Fifty-four autistic participants (16 females) and 31 neurotypical (NT) controls (18 females) matched on age (range 10–25 years) completed tests of precision gripping with their dominant hand (Table 1). Autistic participants were recruited through our research registries comprised of individuals evaluated through the University of Kansas Health System who have consented to be contacted for research purposes, and though community advertisements. NT controls were recruited through community advertisements. ASD diagnoses were confirmed based on Diagnostic and Statistical Manual of Mental Disorders, Edition 5 (DSM-5) (40) criteria and classification criteria from the Autism Diagnostic Observation Schedule, Second Edition (ADOS-2) (41) and Autism Diagnostic Interview – Revised (ADI-R) (42). Autistic participants were excluded if they had a known genetic or metabolic disorder associated with ASD (e.g., Fragile X syndrome) or a full scale IQ (FSIQ) below 60 as measured using the Wechsler Abbreviated Scales of Intelligence, Second Edition (WASI-II) (43). NT participants were excluded if they scored ≥ 8 on the Social Communication Questionnaire (44), reported a history of psychiatric or neurologic disorders, had a family history of ASD in first- or second-degree relatives, had a family history of a developmental or learning disorder, psychosis, or obsessive compulsive disorder in first-degree relatives, or had a FSIQ below 85 as measured using the WASI-II. Participants also were excluded if they had a history of head injury with neurological sequelae, birth injury, or seizure disorder. No participants were taking medications known to affect sensorimotor behavior, including antipsychotics, stimulants, or anticonvulsants at the time of testing (45). All participants had corrected or uncorrected visual acuity of at least 20/40. Adult participants provided written informed consent after a complete description of the study, in accordance with the Declaration of Helsinki and the approved Institutional Review Board study protocol (IRB#: STUDY00140269). For participants under the age of 18 and adults who were under legal guardianship, a parent or legal guardian provided written informed consent, and the participant provided written assent. All study procedures were approved by the local Institutional Review Board.
Table 1
Demographic and clinical characteristics of autistic individuals (ASD) and neurotypical controls (NT)
| | ASD | NT | |
| N | Ratio | | | N | Ratio | | OR |
Sex1 | 54 | 38M:16F | - | | 31 | 13M:18F | - | .309* |
Handedness2 | 54 | 6L:48R | - | | 31 | 4L:27R | - | .845* |
| N | Mean | SD | | N | Mean | SD | t |
Age | 54 | 14.87 | 3.68 | | 31 | 15.90 | 4.11 | -1.16 |
ADOS-CSS | 54 | 6.26 | 2.13 | | - | - | - | - |
RBS-R | 52 | 32.15 | 21.75 | | - | - | - | - |
FSIQ3 | 54 | 97.39 | 16.47 | | 31 | 110.51 | 10.27 | -4.52* |
VIQ3 | 52 | 95.40 | 17.82 | | 31 | 108.45 | 10.63 | -4.18* |
PIQ3 | 53 | 99.32 | 17.09 | | 31 | 110.45 | 12.47 | -3.43* |
BOT-2: Fine Manual Control | 48 | 41.31 | 9.35 | | 23 | 48.43 | 10.21 | -2.83* |
MVC | 54 | 40.17 | 17.41 | | 31 | 49.59 | 17.65 | -2.38* |
ASD: Autism spectrum disorder; NT: Neurotypical; OR: Odds ratio from Fisher’s Exact test; M: Male; F: Female; L: Left-handed; R: Right-handed; SD: Standard deviation; ADOS-CSS: Autism Diagnostic Observation Scale – Composite Severity Score; RBS-R: Repetitive Behavior Scale – Revised; FSIQ: Full-scale intelligence quotient; VIQ: Verbal intelligence quotient; PIQ: Perceptual intelligence quotient; BOT-2: Bruininks-Oseretsky Test of Motor Proficiency, 2nd Edition; MVC: Maximum voluntary contraction.
1 Biological sex is used here. Three autistic participants did not identify as the sex they were assigned at birth (all assigned female). One identified as a transgender male; one as non-binary, and one as gender fluid.
2 Handedness here refers to which hand the participant used for precision grip testing. Five participants either did not complete the Annett (two autistic participants) or had scores on the Annett that did not match the hand they used for writing (three controls).
3 Two autistic participants did not complete all four subscales of the WASI-II, so the two-subscale FSIQ was used for those participants, and PIQ and VIQ subscales were omitted for participants who did not complete both subscales required for those scores.
Clinical Assessments
Participants completed the Wechsler Abbreviated Scales of Intelligence, Second Edition (WASI-II) to assess verbal IQ (VIQ), perceptual IQ (PIQ), and full-scale IQ. The WASI-II is validated for individuals aged 6–89 years. For this study we report scores for the full-scale IQ value that is calculated from all four of the administered subscales, with the exception of two autistic participants who did not complete all four subtests. For these two participants, the two subtest scores are used.
Autistic participants completed the Autism Diagnostic Observation Schedule, Second Edition (ADOS-2) (41) to confirm diagnosis and to quantify severity of autism for analysis. The ADOS-2 is a semi-structured play-based assessment of autistic traits that is the gold-standard diagnostic assessment for autism. It was conducted by a trained research reliable study clinician. Participants in our study were administered module 2 (phrase speech), 3 (verbally fluent children), or 4 (verbally fluent adolescents and adults) according to age and language abilities. The composite severity score (ADOS-CSS) is a standardized score indicating severity of autism that can be compared across modules (higher scores indicate greater severity). The ADOS-CSS are reported and used for analyses.
Autistic participants completed the Repetitive Behavior Scale – Revised (RBS-R) (46, 47), to assess restricted and repetitive behaviors associated with autism. The RBS-R is a questionnaire that asks individuals to rate items from five categories of repetitive behavior (motor stereotypy, self-injurious behavior, compulsions, routines/sameness, and restricted interests). Parents or caregivers completed the RBS-R for participants under 18 years of age, and adult participants completed it as a self-report questionnaire. Higher scores indicate more severe repetitive behavior.
To determine handedness, participants completed the Annett Hand Preference Questionnaire (Annett)(48). The Annett is a 12-item questionnaire that asks the participant which hand they prefer to use for various daily activities (e.g., writing, throwing, using a hammer, etc.), with endorsement of left hand, right hand, or either hand. If left hand is endorsed more than right hand, the participant is considered left-handed, and vice-versa for a classification of right-handed. If left and right hand are equally endorsed, the participant is classified as mixed-handed. Two neurotypical controls scored as mixed-handed and two autistic participants did not complete the Annett, so their dominant hand for precision grip testing was determined based on which hand they used for writing. One control scored as left-handed on the Annett but self-reported as right-handed and used right hand for writing, so they completed precision grip testing with right hand. Handedness counts in Table 1 are based on the hand used to complete precision grip testing.
Participants completed the Bruininks-Oseretsky Test of Motor Proficiency, Second Edition (BOT-2) to assess motor abilities. The BOT-2 is a structured skill-based motor assessment. Participants completed a series of structured motor tasks from three areas: Fine Manual Control, Manual Coordination, and Body Coordination. Composite scores from the Fine Manual Control tests are reported and analyzed for the present study. Higher scores on the BOT-2 reflect better motor performance.
Precision grip testing
Participants completed tests of precision gripping while seated 52cm from a 67cm (27in) Samsung liquid crystal display monitor with a resolution of 1920x1080 and a 120 Hz refresh rate (Fig. 1). Participants sat with the elbow of their dominant hand comfortably positioned at 90° and their forearm resting in a custom arm brace fixed to the table to provide stability during testing. The participants used their thumb and index finger of their dominant hand to press against two opposing precision load cells that were secured to a custom grip device attached to the arm brace. A Coulbourn (V72-25) resistive bridge strain amplifier received analog signals from the load cells, which were converted to digital signals sampled at 100 Hz with a 16-bit analog-to-digital converter (NI USB-6341; National Instruments Corporation). During the first part of the study, ELFF load cells (ELFF-B4-100N; Entran) 1.27cm in diameter were used. Due to normal wear, the ELFF loadcells were replaced with Honeywell loadcells (Model 53, Honeywell International, Inc.) 1.5cm in diameter during the course of the study. ELFF load cells were used for 35.3% of participants (35.5% of controls, 35.2% of autistic participants), and Honeywell load cells were used for 65.5% of participants (64.5% of controls, 64.8% of autistic participants. The voltage-to-Newton calibration was different for each type of loadcell, so a correction calculated from known weights was applied to the force trace after data collection to correct for calibration errors. Additionally, loadcell type was included as a covariate in our analyses to account for differences in loadcell design and calibration effects, including the visual angles of feedback during the task (described below).
[Insert Fig. 1 here]
Prior to precision grip testing, participants completed an assessment of their maximum grip force, or maximum voluntary contraction (MVC) using their dominant hand. Participants completed three trials in which they were asked to press as hard as they could for three seconds. The average of the participant’s maximum force output across these trials comprised their MVC. For the precision gripping tasks, the target force was set at 45% of the participants’ MVC to account for differences in strength across participants.
During the precision gripping task, participants viewed two horizontal bars on the screen (Fig. 1B). A horizontal white force bar moved upward with increased force and downward with decreased force, and a static bar representing the target force was red during periods of rest. The target bar turned yellow to cue the participant to get ready for the start of the trial, and it turned green to cue the participant to begin pressing at the beginning of each trial. Participants were instructed to press the load cells as quickly as possible when the yellow target bar turned green and to keep pressing so that the white force bar stayed as steady as possible at the level of the green target bar until the target bar turned red, marking the end of the trial.
To test the impact of visual feedback and motor memory processes on grip force behavior, participants completed precision grip testing with and without visual feedback. During visual feedback trials, visual feedback was presented continuously throughout the 15s trial. Due to calibration differences for the two types of loadcells that were used during the study and variance in the distance between the screen and the participants’ eyes during naturalistic viewing, the visual angles ranged from .74 to 1.15 degrees per 1N increase in force output. Visual angles between .623-2.023 degrees result in small changes in the spatial amplitude of visual feedback and small changes in force error (Coombes et al 2010). This range is also associated with stable and optimal variability and regularity of grip force in autistic individuals and NT controls (10).
For the trials without visual feedback (“memory guided” trials), the initial part of the trial was the same as for the visually guided trials – the target bar turned from yellow to green, and the participant pressed on the force transducers to match the white bar to the level of the green target bar. After three seconds, the white force output bar disappeared, and the participant was instructed to continue pressing at the same level until the target bar turned red (12s after the visual feedback was removed). Participants completed blocks of five trials of each condition using their dominant hand (5 trials x 2 conditions = 10 trials). Trials were 15s in duration and alternated with 15s rest periods. Each block was separated by 30s of rest. The target force was set to 45% of the participant’s MVC for all trials. The order of the blocks was pseudo-randomized and counterbalanced across participants.
Data processing
Grip force data were processed using custom applications developed by our lab in R and MATLAB (MathWorks, Inc., Natick, Massachusetts). Trials were excluded if the load cells were not properly re-zeroed between trials or if there were indications that the participant was not following instructions (e.g., the mean force exceeded twice the target force, there was evidence that the participants used fingers other than dominant hand index finger and thumb to press, participants stopped pressing during the trial). For the memory guided condition, trials were excluded if the participant did not reach a stable level of force output within ± 2 Newtons of the target force before visual feedback was removed. Based on these criteria, 19.8% of trials were excluded for the ASD group (12.3% of visually guided and 26.5% of memory guided feedback trials) and 1.3% of trials were excluded for the control group (1.9% of visually guided and .7% of memory guided trials). Trial-level data were averaged for each participant within each condition. Participants needed to have at least two useable trials of a condition for their data to be included. Four autistic participants were excluded from analyses due to insufficient data. An additional five autistic participants had insufficient data for only the memory guided condition, and one had insufficient data for only the visual feedback condition. Final analyses included 50 autistic participants and 31 NT controls with valid data for at least one condition.
To compare force output across conditions, the sustained force output for each trial was analyzed. To account for the differences in trial structure between conditions, only the last 12s of each trial were used for analysis of sustained force output. This 12s phase corresponds to the segment of the memory guided trials where visual feedback was not available and the analogous segment of the trials with visual feedback.
The force traces for each trial were low-pass filtered via a double-pass fourth-order Butterworth filter at a low-pass cutoff of 15 Hz following previous studies from our lab (18, 22). Sustained force data were linearly detrended to account for drift in participants’ force output over the duration of the trial. The mean force of the sustained force data divided by the target force was used as a measure of force accuracy, such that values close to 1 reflect greater accuracy. To assess force variability, the standard deviation (SD) of the force time series was examined. To test the time dependent regularity of the force time series, sample entropy (SampEn) was calculated for each trial (49, 50). SampEn is defined as the natural logarithm of the conditional probability that two similar sequences of m data points in a timeseries of a given length (N) remain similar within a tolerance level (r) at the next data point in the series. SampEn returns a value between 0 and 2. Lower values of SampEn indicate greater regularity of the timeseries (e.g., a sine wave, with its predictable oscillating pattern, would have a SampEn value near 0). Parameter settings for SampEn calculations were m = 2 and r = .2 x SD of the timeseries. The timeseries length was 1200 data points (12s sampled at 100 Hz). The sampenc.m function (for MATLAB) from the PhysioNet Toolbox (51, 52) was used to calculate SampEn values for each trial.
To characterize the trajectory of force output after visual feedback was removed during the memory guided trials, models were fit to the force traces for each trial. The model consisted of three segments: 1) a horizontal line fit to the stable force output at the beginning of the trial, starting before visual feedback was removed, 2) a logarithmic function fit to the decay in force output after visual feedback was removed, 3) for trials where participants reached a stable force output after their force decayed and before the end of the trial, a horizontal line was fit to the data to model this secondary stable force output. Latency, slope, and magnitude of the force decay were analyzed. Decay latency was measured as the difference between the removal of visual feedback and the beginning of the logarithmic decay model segment. The decay slope was calculated as the log slope of the logarithmic decay model segment. For 2.4% of trials (ASD: 2.5% and NT 2.0%), the force showed little to no decay and was best modeled using a linear fit rather than a logarithmic fit. Descriptive statistics were calculated and reported for the trials with linear decay slopes, but these trials were not factored into the participants’ trial averages or the linear regression models, as they are not directly comparable to the trials with logarithmic slopes. The magnitude of decay was calculated by taking the difference between the force output at the onset of the logarithmic (or linear) decay model segment and the end of the decay model segment or the end of the trial, if force did not stabilize before the end of the trial. This value was then converted to percentage of target force by dividing the raw difference by the participant’s target force to account for the differences in initial force.
Statistical Analysis
Force accuracy, SD, and SampEn were analyzed using separate linear multilevel mixed effects models (MLM) (53, 54) with the lme4 package in R version 4.0.0 (53). MLM allows for the analysis of within- and between-subjects fixed effects while allowing within-subjects effects to vary randomly and is robust to missing data. Task condition (visually guided, memory guided) was included as a level 1 predictor. Group (ASD, NT) and age were included as level 2 predictors. For all dependent variables, the models also included a two-factor covariate to account for different load cells used during testing (“load cell type”). Random intercepts of participant also were included in the models.
Initial models for force accuracy, SD, and SampEn included the three-way interaction of group x task condition (visually guided, memory guided) x age, all relevant two-way interactions and main effects terms, as well as the covariate for load cell type. To maintain the most parsimonious models possible, other 3- way interactions were not included. Models were fit using the maximum likelihood approach to allow for model comparisons. Terms were removed systematically, and model fit was compared between the previous model and the model with the removed term using likelihood ratio tests. Terms that did not significantly improve model fit (p < 0.05), based on the model comparisons, were not included in the final models. Satterthwaite’s method was used to calculate degrees of freedom for the final model and post hoc comparisons (55). Due to the inherent challenge in determining denominator degrees-of-freedom and calculating p-values for MLMs, we treated the t-value as a z-value and used a z > 1.96 threshold as an additional guideline for determining whether terms explained significant variance in the model (55).
The latency, slope, and magnitude of force decay following the removal of visual feedback were analyzed using separate linear regression models with the lm (linear model) function in R. Group (ASD, NT), age, and the group x age interaction were included as predictors. For all dependent variables, the models also included a covariate for load cell type.
Simple coding was used for group (NT = -0.5, ASD = 0.5), task condition (memory guided = -0.5, visually guided = 0.5), and sex (male = -0.5, female = 0.5). Age was log10 transformed. SD, SampEn, and decay magnitude were log10 transformed and decay slope and decay latency were square root transformed to correct for skewed distributions. Based on this coding system, the intercept for each model represented the grand mean of the sample.
Pearson correlations were used to assess the relation between experimental variables and ASD symptom severity measured using the ADOS Composite Severity Score (ADOS-CSS) as well as repetitive behaviors measured using the RBS-R. Pearson correlations also were used to assess the relation between visuomotor and motor memory behaviors and IQ for each group.
To determine whether visual feedback guided motor control and motor memory during precision gripping relate to clinically relevant fine motor skills, Pearson correlations were run between the motor variables and BOT-2 Fine Manual Control subscale scores. For each set of correlations, p-values were adjusted using false discovery rate (FDR) to limit Type I error.