- 1Graduate School of Engineering, School of Science and Engineering, Kokushikan University, Setagaya, Japan
- 2Graduate School of Biomedical Sciences, Nagasaki University, Nagasaki, Japan
- 3Graduate School of Engineering, Chiba University, Chiba, Japan
- 4Medical Institute of Developmental Disabilities Research, Showa University, Tokyo, Japan
- 5Faculty of Humanities, Wayo Women's University, Chiba, Japan
People with autism spectrum disorder (ASD) exhibit atypicality in various domains of behavior. Previous psychophysiological studies have revealed an atypical pattern of autonomic nervous system (ANS) activation induced by psychosocial stimulation. Thus, it might be feasible to develop a novel assessment tool to evaluate the risk of ASD by measuring ANS activation in response to emotional stimulation. The present study investigated whether people with ASD could be automatically classified from neurotypical adults based solely on physiological data obtained by the recently introduced non-contact measurement of pulse wave. We video-recorded faces of adult males with and without ASD while watching emotion-inducing video clips. Features reflective of ANS activation were extracted from the temporal fluctuation of facial skin coloration and entered into a machine-learning algorithm. Though the performance was modest, the gradient boosting classifier succeeded in classifying people with and without ASD, which indicates that facial skin color fluctuation contains information useful for detecting people with ASD. Taking into consideration the fact that the current study recruited only high-functioning adults who have relatively mild symptoms and probably developed some compensatory strategies, ASD screening by non-contact measurement of pulse wave could be a promising assessment tool to evaluate ASD risk.
Introduction
Autism spectrum disorder (ASD) is a group of closely related developmental disabilities characterized by a collection of symptoms including persistent deficits in social communication and reciprocal interaction, repetitive behavior, and restricted interest (1). Although genetics is considered to play major role in the development of ASD, despite years of effort by researchers, the exact cause of ASD still remains elusive (2). As such, there exists no definitive cure for ASD.
Public awareness of ASD has risen significantly in recent years (3), but several surveys suggest the possibility that a significant portion of the ASD population remains undiagnosed (4). Importantly, some of them remain undiagnosed well into adulthood and experience great distress due to failure in adaptation to occupational and social environments (5–7). The difficulty in detecting ASD risk in adults is attributable to several causes. First, adults with milder cases of ASD can sometimes compensate for their atypicality to some degree by adjusting their behavior with deliberate effort. Second, misdiagnosis is prone to be made because of the comorbidity of other psychiatric conditions.
One promising way to reduce the number of undiagnosed adults with ASD is to develop a supplementary assessment tool for ASD. The currently accepted protocol for ASD diagnosis requires hours of interviews and assessment by clinicians and is quite labor-intensive. Such strict protocols are relevant from the perspective of suppressing overdiagnosis of ASD. At the same time, it is also necessary to establish a convenient and handy assessment tool for primary screening to broadly catch those who “might” have ASD.
Recent studies on digital phenotyping (8) have pointed to the potential of semiautomated quantification of behavioral phenotypes in the primary screening of ASD. In these studies, behavioral features of people with ASD were captured by information technology, such as video analysis of facial movement by computer vision (9–11) and bodily movement analysis using cost-effective sensors [Kinect sensors and touch sensors implemented in tablet computers; (12, 13)]. Some studies have entered the extracted features into machine learning and succeeded in classifying ASD and neurotypical counterparts (11–13).
Classical theories of ASD have claimed that atypicality in emotional function causes ASD-like behaviors (14, 15). Although not in its original form, some of these theories have gained support from studies showing poor performance in emotion recognition (16, 17) and atypical patterns of neural activation and morphology in emotion-associated regions such as the amygdala (14, 18). In addition to atypical development of emotion-related neural regions, a number of studies indicate that people with ASD, especially pediatric patients, show poor ability of emotional regulation (19) as can be seen in occurrences of sudden burst of emotion, so-called “meltdown.” Taken together, people with ASD have atypicality in both the induction of emotion and emotional regulation, which leads to an atypical pattern of emotional responses.
Taking such atypical emotional responses in ASD people into consideration, measurement of emotional response might be useful in screening adults with ASD. Emotional response can be objectively quantified by measuring the activation of the autonomic nervous system (ANS), which comprises the sympathetic and parasympathetic nervous system. Psychophysiological studies have linked the relative increase of sympathetic over the parasympathetic nervous system to a high arousal state (20). In addition, several studies have found atypical patterns of ANS activation in people with ASD, especially in response to social situations. For example, studies have revealed reduced sympathetic activation in adults with ASD in response to social interaction (21) and psychosocial stress (22). Similar atypicality has also been observed in pediatric samples (23), although the exact pattern of atypicality differs among studies (24, 25). Previous studies on ANS activation in people with ASD also revealed increased heart rate and reduced respiratory sinus arrhythmia during rest (26, 27) and sleep (28) in individuals with ASD.
Currently, measurement using a polygraph system is deemed the gold standard for ANS activity measurement. However, polygraph systems are relatively expensive and thus difficult to introduce in many facilities. In addition, attachment of sensors onto the skin surface makes people with ASD, especially those with hypersensitivity to tactile stimulation (29), uncomfortable. A recently introduced method of non-contact measurement of pulse waves (30, 31), an index of ANS activation, could be an alternative technique to a polygraph. In this method, the timing of the pulse wave is detected by analyzing a slight change in facial skin coloration. Immediately after cardiac pulsation, arterial blood rich in fresh-red oxygenated hemoglobin circulates through the blood vessels. Hemoglobin pigmentation on the skin surface is slightly intensified synchronously with this event. The temporal interval between neighboring peaks of pigmentation intensity reveals the balance between sympathetic and parasympathetic nervous system activations (32, 33) and, thus, certain aspects of emotional response.
The present study examined whether this non-contact measurement of ANS activation is assistive in distinguishing between neurotypical adults and adults with ASD. We extracted features reflective of ANS activation from video-recorded images of the participant's face while experiencing film-induced emotions (34, 35). These features were fed to machine-learning algorithms to ascertain whether pulse waves extracted from video-recorded images contain enough information to distinguish ASD adults from neurotypical ones. We focused on the phasic change of ANS activation during emotion-inducing film viewing from baseline.
Method
Participants
A total of 21 typically developed (TD) adult males (M = 31.7 years old; SD = 6.03) and 17 adult males with ASD (M = 31.5 years old; SD = 5.2) participated in the present study after giving written informed consent. We excluded female participants so as not to increase variance in responses in the ASD group because previous studies have found phenotypic differences between males and females with ASD (36, 37). The TD participants were recruited by word of mouth and recruitment agency. They had no ties or relations with the ASD participants. The average age did not differ significantly between the ASD and TD groups, t(36) = 0.07, p = 0.94. All participants with ASD were referred to Showa University Hospital by physicians from other clinics. The inclusion criteria of the present study were as follows: males, Wechsler Adult Intelligence Scale-III, full intelligence ≥70 (high functioning evaluated by WAIS), and a formal diagnosis of ASD based on DSM-5 (1). The exclusion criteria were comorbid personality disorders. Three of them had a comorbidity of mood disorders. Out of 17 ASD participants, nine were on medication (antidepressant, n = 5; anxiolytic, n = 4; soporific, n = 3; antipsychotic, n = 4; anticonvulsant, n = 3). None was taking behavioral therapy. A psychiatrist made a diagnosis according to the diagnostic criteria of the DSM-5 for ASD based on a consensus among a team comprised of experienced psychiatrists and a clinical psychologist who carried out clinical interviews. All participants satisfied the Autism Diagnostic Observation Schedule diagnostic criteria at the time of the interview (M = 11.4, SD = 2.8). The diagnosis was reconfirmed after at least a 2-month follow-up period. The protocol of this research was approved by the ethical committees of the Graduate School of Biomedical Sciences in Nagasaki University and the Graduate School of Engineering in Chiba University, according to the Declaration of Helsinki.
Apparatus and Stimuli
The stimuli included six video clips with lengths of about 1.5–4.5 min intended to induce happy, sad, fearful, and neutral emotional states. These video clips were taken from a TV show or films to induce these emotional states, the details of which are described in Table 1. Four of the video clips depicted social interactions among several people. Two more video clips, that is, popular line drawing animation and a video clip of factory machines without human figures, were included in the stimulus set because previous studies have indicated that some people with ASD are prone to show interest in these contents (38, 39). Due to the scarcity of studies on emotion induction by films in Japanese samples, a new set of emotion-inducing films was created for the present study. We first screened Asian films (most of them, Japanese ones) that seemed suitable for emotion induction. Scenes including graphic content and violence were avoided so as not to elicit shock from the participants. Then, a small pilot study was carried out in which our lab members (one female and five males) evaluated 13 Asian films1 intended to induce sad (two films), fear (three films), happy (two films), and neutral (two films) emotions and two machine factory scenes (two films). Based on the selectivity and intensity of emotional state evaluation, we chose the six films summarized in Table 1.
The participants viewed a 19-inch monitor with their head stabilized on a chin rest. Their faces were lit up using two fluorescent lights. The viewing distance was 75 cm. An LED light located near the participant's chin lit up at the start and the end of each video clip. The participants' faces and the LED lights were recorded by an RGB camera (situated between the monitor and the chin rest) for industrial use (ARGO Corporation, Suita, Osaka, Japan). A picture of the experimental setup is shown in Figure 1A.
Figure 1. (A) A picture of experimental setting taken during preparation. The man in the picture is one of the authors, not a participant. During the actual experiment, participant held a keypad with their head stabilized on a chin rest. (B) Typical pattern of pulsatile fluctuation in hemoglobin intensity recorded from one ASD participant.
Procedure
Face Movie Recording
At the start of each trial, a black background was shown for 40 s. Then, a video clip was presented. The video clip included 10 s presentation of a brief explanation about the context of the scene. The participants were instructed to view the video clip without overt movement quietly. Immediately after the end of the video clip, seven track bars with 11 tics were presented on the screen. The participants' task was to answer their emotional states, that is, the levels of happiness, anger, sadness, fear, disgust, unpleasantness, and arousal, using the track bar. After the completion of the subjective evaluation of emotional state, the next trial started. Six video clips were separated into three blocks with interblock rest intervals. The order of the video clips was determined randomly.
Self-Administered Questionnaires
After the completion of the facial movie recording, the participants answered a battery of self-administered paper-and-pencil questionnaires in front of the experimenter. The questionnaires included Japanese translations of the IRI [Interpersonal Reactivity Index; (40)] and TAS-20 [The 20-item Toronto Alexithymia Scale; (41)]. The IRI measures empathic tendencies, and it comprises four subscales: perspective-taking, fantasy, empathic concern, and personal distress (42), while TAS measures the level of difficulty in identifying feelings, difficulty in describing feelings, and externally oriented thinking (41, 43).
Image Processing
Recorded movies of the participants' faces were analyzed offline according to the algorithm described by Fukunishi et al. (31, 44). The skin color of the small facial region near the cheekbone was analyzed. First, skin coloration in each frame was separated into hemoglobin pigmentation, melanin pigmentation, and shade information by independent component analysis (45). Because shade information is separated from hemoglobin intensity, the intensity of hemoglobin pigmentation can be assessed independently of lighting condition. Then, the temporal fluctuation of hemoglobin intensity across video frames was analyzed.
The temporal course of hemoglobin pigmentation intensity shows pulsatile fluctuation as shown in Figure 1B. From the temporal course of hemoglobin pigmentation fluctuation, 20 indices of ANS activation were computed for each video clip for each participant. The details of the extracted indices are summarized in Table 2. These indices were measured during 30 s baseline just before video clip presentation and during video clip viewing. We used Kubios-HRV software (46) for computing some of the indices. The phasic change in these indices from baseline was calculated by subtracting index values during baseline from those during video clip viewing. Thus, six change-from-baseline values were obtained for each of the 20 indices, yielding a total of 120 features for machine learning.
Group Classification by Machine Learning
To ascertain whether participants with and without ASD can be automatically classified by 120 features reflective of ANS activation obtained by non-contact measurement, data were analyzed using a supervised machine-learning algorithm. First, 120 features were standardized and then compressed to principal components (PCs) by principal component analysis (PCA). Then, the scores of PCs were fed into a gradient boosting (GB) classifier. GB is one of the ensemble learning algorithms, which is known to show good prediction performance and generalizability, i.e., the ability of the classifier to predict the class of unknown data, over the other types of machine learning algorithms in many classification problems. In ensemble learning, a set of weak predictors with limited prediction ability is created, instead of a single strong predictor with high prediction ability. Then, the final prediction is made by aggregating the predictions of these weak predictors. Intuitively, in GB, weak predictors are added sequentially to the decision trees so that the addition of a new predictor gradually minimizes the discrepancy between observed data and prediction. There are two popular ensemble learning algorithms, GB and random forest. Though random forest excels at generalizability in some cases, GB can show superior performance when its hyperparameters are tuned appropriately.
The GB classifier was cross-validated by a leave-one-out procedure. Specifically, data from one participant were used as the unknown test data, and the data of the remaining 37 participants were used to train the classifier. After training, the classifier predicted which group, i.e., TD or ASD, the unknown test data belonged to. The performance indicators, i.e., accuracy, sensitivity, specificity, precision, and F1 score, were calculated according to the formula below after repeating this procedure 38 times by using the data of every participant as test data.
where TP (true positive), FP (false positive), TN (true negative), and FN (false negative) represent the frequencies with which the ASD participant was correctly identified as ASD, the TD participant was misclassified as ASD, the TD participant was correctly classified as TD, and the ASD participant was misclassified as TD, respectively.
To obtain a robust estimate of performance, we repeated this process 20 times by changing the random number seed.
Basic Emotion Classification by Machine Learning
Among the six video clips, happy, fear, sad, and neutral clips are intended to induce a particular emotional state. In order to examine whether these four clips induced differential physiological states in ASD and neurotypical adults, a GB classifier was trained to discriminate these four clips based on the 20 ANS indices. The features were first compressed by PCA, and then the GB classifier was trained for the TD and ASD groups separately.
In the ASD group, there were emotional category (4) × participant (17) = 68 data. In the leave-one-out procedure, one of the 68 data was treated as unknown test data and the GB classifier was trained on the basis of the remaining 67 data. Accuracy rate and confusion matrix were computed after repeating this procedure 68 times using every one of the 68 data as the test data. Accuracy and confusion matrix were computed in the TD group according to essentially the same procedure using emotional category (4) × participant (21) = 84 data. Computation of accuracy and confusion matrix was repeated 20 times by changing the random seed.
In the next step of analysis, we tested whether emotional states are less clearly differentiated in ASD than TD participants (47). If this hypothesis is correct, the pairwise similarities among physiological responses induced by the four emotion-inducing clips should be higher in the ASD group than in the TD group. To examine this possibility, the cosine similarity of physiological states induced by the four video clips was calculated based on the confusion matrices. The cosine similarity between emotion j and emotion k, CSjk, is computed according to the following formula where Nij represents the frequency with which emotion j was classified as emotion i (happy, fear, sad, neutral).
Statistical Analysis
Subjective ratings of emotional states were entered into an analysis of variance (ANOVA) with the between-participant factor of group (2; TD, ASD) and the within-participant factors of clip (6; happy, fear, sad, neutral, line drawing, machine) and emotional state (7; happiness, anger, fear, sadness, disgust, surprise, unpleasantness, arousal). When an interaction reached significance, simple main effect analysis was conducted to clarify its source.
In leave-one-out, the GB classifier gives the probability estimate that the test data belong to the ASD class. In the present study, leave-one-out cross-validation was repeated 20 times, yielding 20 probability values for each participant. We tested whether the average probabilities assigned to ASD participants were higher than those assigned to TD participants as expected by the Mann–Whitney U-test. The chi-squared test was used to examine if the proportions of participants classified into the ASD group differed between the ASD and TD groups. It was tested whether the mean accuracy of classification was significantly above chance (50%) by two-tailed t-test.
Group difference in subscores of IRI and TAS-20 was tested by two-tailed t-tests. To see if these traits are associated with classification results by the GB classifier, Spearman's rank-order correlation was tested between averaged probabilities and the scores of IRI and TAS-20.
Results
Group Difference in Subjective Emotional States
The mean and standard deviation of the subjective ratings in each condition are shown in Figure 2.
Figure 2. Subjective rating of each emotional state given to each of the six video clips in the ASD (left) and TD (right) groups. The center of the radar chart is zero, and the outermost point of each axis is eight. The thick line and the width of the shaded area along each axis represent mean and standard deviation, respectively.
The ANOVA revealed significant main effects of clip, F(5, 180) = 14.53, p < 0.0001, = 0.29, and emotional state, F(6, 216) = 30.53, p < 0.0001, = 0.46. There were also significant interactions between group and emotional state, F(6, 216) = 4.08, p < 0.001, = 0.10, and between clip and emotional state, F(30, 1, 080) = 35.04, p < 0.0001, = 0.49. All these main effects and interactions were further qualified by a significant three-way interaction between group, clip, and emotional state, F(30, 1, 080) = 1.56, p = 0.028, = 0.042. No other effects reached significance, Fs < 1.0, ps > 0.4. In order to clarify the source of the significant three-way interaction, two-way ANOVA with the between-participant factor of group (2) and within-participant factor of emotional state (7) was carried out for each clip.
Happy Clip
The ANOVA revealed a significant main effect of emotional state, F(6, 216) = 59.67, p < 0.0001, = 0.62, and a significant interaction between group and emotional state, F(6, 216) = 4.48, p = 0.0003, = 0.11. The main effect of group failed to reach significance, F(1, 36) = 0.03, p = 0.85, < 0.001. Simple main effect analysis revealed significantly higher “happy” rating in TD than in the ASD group, F(1, 36) = 5.72, p = 0.02, = 0.14. People with ASD also gave higher “unpleasantness” ratings than TD people, F(1, 36) = 4.87, p = 0.034, = 0.12. A simple main effect of group failed to reach significance in the other emotional states, Fs < 3.4, ps > 0.07.
Fear Clip
The main effect of emotional state reached significance, F(6, 216) = 36.34, p < 0.0001, = 0.50. The main effect of group failed to reach significance, F(1, 36) = 0.39, p = 0.53, = 0.01. The interaction between group and emotional state failed to reach significance either, F(6, 216) = 1.88, p = 0.086, = 0.05.
Neutral Clip
The main effect of emotional state reached significance, F(6, 216) = 13.43, p < 0.0001, = 0.27. Neither the main effect of group, F(1, 36) = 0.20, p = 0.66, = 0.006, nor the interaction between group and emotional state, F(6, 216) = 0.29, p = 0.94, = 0.01, reached significance.
Sad Clip
The main effect of emotional state, F(6, 216) = 36.45, p < 0.0001, = 0.50, and the interaction between group and emotional state reached significance, F(6, 216) = 3.36, p = 0.0035, = 0.085. The main effect of group failed to reach significance, F(1, 36) = 0.93, p = 0.34, = 0.025. Simple main effect analysis revealed a significantly higher “fear” rating in the ASD than in the TD group, F(1, 36) = 9.12, p = 0.005, = 0.20. A simple main effect of group reached significance for no other emotional state, Fs < 1.9, ps > 0.17.
Line Drawing Animation
The main effect of emotional state reached significance, F(6, 216) = 34.37, p < 0.0001, = 0.49. This main effect was qualified by a significant interaction between group and emotional state, F(6, 216) = 2.62, p = 0.018, = 0.07. The main effect of group failed to reach significance, F(1, 36) = 0.06, p = 0.80, = 0.002. Simple main effects of group failed to reach significance, Fs < 3.5, ps > 0.07.
Machine Factory Clip
Only the main effect of emotional state was significant, F(6, 216) = 22.34, p < 0.0001, = 0.38. Neither the main effect of group, F(1, 36) = 0.61, p = 0.45, = 0.017, nor the interaction between group and emotional state reached significance, F(6, 216) = 0.095, p = 0.99, = 0.003.
Automatic Group Classification
We first tested whether the performance of the GB classifier was dependent on the number of PCs used as features. Thus, GB classifiers were trained and tested by using scores of 1–10 PCs as features. The PCs were entered into the classifier in descending order of the amount of variance explained. For example, when the number of PCs was two, the PCs with the first and second largest variance explained were used as the features. The performance of GB classifiers was compared by the accuracy and F1 score as performance indicator. The classifier showed its best performance when the number of PCs was eight.
The Mann–Whitney U-test revealed a significantly higher probability of being classified into ASD by the classifier in the ASD than in the TD group, U = 117, p = 0.037. When the threshold was set to 0.5, 13 out of 17 (76.4%) ASD participants and 8 out of 21 (38.1%) TD participants were classified into the ASD group. The proportions of participants classified into the ASD group differed significantly between the ASD and TD groups, = 5.59, p = 0.018.
Mean accuracy across 20 repetitions of cross-validation was 0.68 (SD = 0.015). Accuracy was significantly higher than chance (0.5), t(19) = 54.46, p < 0.001. Averages of sensitivity, specificity, precision, and F1 score were 0.76 (SD = 0.021), 0.63 (SD = 0.021), 0.63 (SD = 0.016), and 0.68 (SD = 0.015).
Feature importance of each PC is shown in Figure 3. As can be seen in this figure, the fourth and sixth PCs contributed most to the classification followed by the first and second PCs. For completeness, we compared scores of the eight PCs between the ASD and TD groups. No group difference was found, ts < 1.7, ps > 0.10.
Figure 3. Mean feature importance of the right PCs used as features in group classification by GB classifier. PCn is the PC that explained the n-th largest amount of variance.
Relationship Between Classification Results and Scores of IRI and TAS-20
The average and standard deviations of IRI and TAS-20 are summarized in Table 3. Two ASD participants failed to complete TAS-20. In IRI, adults with ASD showed lower perspective-taking and higher personal distress scores than their neurotypical counterparts. The results also revealed a higher score in every factor of TAS-20 in the ASD than in the TD group.
Table 3. Mean and standard deviation of IRI and TAS-20 scores in each group, and the results of group comparison of each score by two-tailed t-test.
The results of correlation coefficients between questionnaire results and the probabilities of being classified as ASD assigned by the GB classifier are summarized in Table 4. After significance threshold correction by Bonferroni's procedure, there was no significant correlation.
Table 4. Spearman's rank-order correlation coefficients between averaged probability and scores of IRI and TAS-20.
Automatic Classification of Emotional States Within Each Group
The performance of the GB classifier was dependent on the number of PCs. Thus, we evaluated the performance of the GB classifier by changing the number of PCs in the same manner as described above. In the TD groups, the best performance was achieved when the number of PCs was six, while one in the ASD group. Thus, the number of PCs was set to these values.
The averaged accuracy was 0.34 in the TD group and 0.31 in the ASD group. The confusion matrix for each group is shown in Table 5 together with the hit rate, i.e., the probability that the classifier made a correct classification of each emotion.
The cosine similarities between each pair of emotional states are listed in Table 6. As can be seen in this table, cosine similarities between happy and neutral, and between happy and sad, were relatively low compared with the other pairwise combinations in both ASD and TD groups.
Discussion
Recent advances in image processing technology have enabled the quantification of ANS activation by non-contact measurement of pulse waves from images of the skin surface (30, 31, 44). Here, we used this technology to measure ANS activity in adults with ASD in response to six types of video clips.
At the behavioral level, ASD people gave lower “happy” rating (for happy clip) and higher “fear” (for sad clip) and “unpleasant” (for happy clip) ratings than TD people. The overall pattern indicates a stronger negativity bias in emotional reactions to video clips in ASD individuals. Anxiety disorder is a common comorbidity of ASD (24, 48) and sometimes leads to misdiagnosis of patients with mood disorders as ASD (7). The present finding of negativity bias in individuals with ASD is in line with this clinical observation. At the same time, the present results tell nothing about the cause of this bias. This could be due to hyperactivation in neural centers of emotion induction, such as the amygdala and anterior cingulate cortex (49), or the atypical interpretation of the physiological state induced by video clips.
We tried a semiautomatic classification of TD and ASD groups by submitting the features reflective of ANS activation into a gradient boosting classifier. The classifier succeeded in classifying unknown participants into either the TD or ASD group above the chance level. Furthermore, the probability of belonging to the ASD group estimated by the best performance classifier was significantly higher in the ASD group than in the TD group, which further supports the proposition that non-contact measurement of pulse waves captures atypical responses useful for classifying ASD from TD adults.
At the same time, we have to point out that our success was only modest. Although the accuracy rate was significantly above chance, the values of performance indicators were far from being acceptable for practical use. The relatively poor performance of the classifier is not surprising considering the fact that we focused on a single facet of peripheral nervous system activation, i.e., pulse wave. Needless to say, faces covey abundant information, and some of the previous studies on digital phenotyping of ASD made efficient use of facial information in automatic discrimination of ASD and TD samples (10, 11, 50). Thus, a combination of the other sources of facial information together with non-contact measurement of pulse waves might lead to the establishment of a classifier with better prediction performance.
The relatively low performance of our classifier is also attributable to heterogeneity in ASD participants. It is now well-acknowledged that people with ASD are a heterogeneous group and the symptomatic profiles differ greatly among patients (51). This observation gained further support from the emerging view that each symptom of ASD derives from a different set of etiological factors (52, 53). Given this, it is highly possible that participants in our ASD group had heterogeneous characteristics, which might have yielded a variance in ANS responses to emotional clips. Of particular relevance to the present study, it is possible that participants with hypersensitivity to visual and auditory stimulation might have shown different responses to the video clips from those without hypersensitivity in the present study. Unfortunately, the sample size of the present study is not large enough to examine the effects of symptomatic heterogeneity on physiological responses to emotional clips. This point is surely an important topic for future study.
The results of self-administered questionnaires revealed lower perspective-taking in ASD. Considering that the vicarious experience of the characters' emotion is at the core of emotional experience in film viewing, one might conceive that a lower tendency of perspective-taking might explain a different pattern of ANS activation in ASD compared with TD participants. However, the results of correlation analysis do not support this view; correlation analyses did not show a reliable association between IRI scores and averaged probability belonging to the ASD group generated by the classifier. All of the three scores measured by TAS-20 were higher in the ASD group, which is consistent with the reported strong alexithymic tendency in ASD (54).
The classification results of basic emotional states revealed comparable accuracy in the classification of four emotional states by gradient boosting classifier in the TD and ASD groups. The averaged accuracy rate of classification across all the emotional categories was numerically low in the ASD than in the TD group, but this is probably because the sample size of adults with ASD was smaller than that of TD adults; the classifier for ASD was disadvantaged due to the smaller size of the training dataset. In both ASD and TD groups, cosine similarity was high for the combination between happy and fear, but was relatively low for the combinations of happy and two low arousal emotions (neutral, sad). This indicates that indices extracted from pulse wave are sensitive to arousal level, consistent with previous studies that show the association between arousal state and ANS activation (20). Interestingly, hit rate for happy emotion was drastically lower in the ASD than in the TD group. Actually, the pattern of ANS activation induced by the happy clip was more likely to be confused with that induced by the fearful clip in the ASD than in the TD group. This might indicate a less clear differentiation of emotional state in ASD (47). This finding also dovetails with the proposal that people with ASD show atypical pattern of activation in fear-related neural circuits (18, 55) by emotional stimulation. However, definitive conclusion awaits further empirical studies. An alternative explanation for this result is that the film set used in the present study was not optimized to induce discrete emotional states in ASD participants. Given the prevalent use of films for emotion induction (34, 56) in TD population, it is desirable to establish a standardized set of emotion-inducing films suitable for people with ASD.
The present study showed the potential of non-contact ANS measurement in the automatic classification of adult males with and without ASD. The number of studies reporting the successful application of affordable non-contact measurement to digital phenotyping of atypical behavior in ASD has increased recently (8, 11–13). Combined with these previously established methodologies, the current technique could serve as an assessment tool to screen hitherto undiagnosed or misdiagnosed adults with ASD (5–7). However, there are several limitations that qualify our conclusion. First, many of the relevant studies, including the current one, fail to include participants with other psychiatric conditions, and as such, it remains uncertain whether non-contact measurement of pulse wave can discriminate ASD from patients with other psychiatric conditions. Second, the present study recruited only adult males with and without ASD. Thus, it is unclear whether the current technology could contribute to the classification of females or pediatric samples with ASD from their neurotypical counterparts. Third, the distribution of ASD and TD people in the present study does not mirror the prevalence of ASD in the population. The prevalence of ASD among males is estimated to be around 1 in 50–60 (57), though the reported number varies widely among studies. Given the relatively low prevalence rate, the sensitivity/specificity of our classifier is unclear when it is applied to actual clinical settings. Fourth, we used a video camera with a 30-Hz frame rate for cost-effectiveness. The use of a low frame rate might have made the estimation of some of the ANS activation indicators somewhat unreliable. In other words, the performance of the classifier might be improved by integrating video cameras with a higher frame rate. Considering these, it would be desirable to test the reliability of the current technique by recruiting a more diverse population and using more suitable equipment in future studies.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by The ethical committee of Graduate School of Biomeidcal Sciences in Nagasaki University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
HD, NT, and CK conceived of the research. HD, CK, KM, RM, and TN collected data. NT, KM, RM, and TN carried out the analysis of video images. HD analyzed the data and wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was financed by Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for Scientific Research(C) Grant No. 17K01904 to HD and was supported by a grant from the Joint Usage/Research Program of Medical Institute of Developmental Disabilities Research, Showa University to HD.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Kazuyuki Shinohara for his generosity allowing us to use the facilities in Nagasaki University.
Footnotes
1. ^The 13 films tested included “Lullaby of Narayama (1983)” (sad), “Dark Water (Japanese Version, 2002)” (fear), “A Tale of Two Sisters (2003)” (fear), “Gottu-Ee-Kanji (Japanese TV comedy show)” (an episode other than the one used in the experiment; happy), and “Maboroshi (1995)” (neutral) and a video of machine factory (other than the one used in the experiment) in addition to six clips selected.
References
1. American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders 5th ed. Arlington, VA: American Psychiatric Association (2013). doi: 10.1176/appi.books.9780890425596
2. Anderson GM. Autism biomarkers: challenges, pitfalls and possibilities. J Autism Dev Disord. (2015) 45:1103–13. doi: 10.1007/s10803-014-2225-4
3. Matson J, Kozlowski AM. The increasing prevalence of autism spectrum disorders. Res Autism Spectr Disord. (2011) 5:418–25. doi: 10.1016/j.rasd.2010.06.004
4. Baron-Cohen S, Scott FJ, Allison C, Williams J, Bolton P, Matthews FE, et al. Prevalence of autism-spectrum conditions: UK school-based population study. Br J Psychiatry. (2009) 194:500–9. doi: 10.1192/bjp.bp.108.059345
5. Foran LL. Identifying autism spectrum disorder in undiagnosed adults. The Nurse Pract. (2018) 43:14–8. doi: 10.1097/01.NPR.0000544285.02331.2c
6. van Niekerk ME, Groen W, Vissers CT, van Driel-de Jong D, Kan CC, Oude Voshaar RC. Diagnosing autism spectrum disorders in elderly people. Int Psychogeriatr. (2011) 23:700–10. doi: 10.1017/S1041610210002152
7. Stagg SD, Belcher H. Living with autism without knowing: receiving a diagnosis in later life. Health Psychol Behav Med. (2019) 7:348–61. doi: 10.1080/21642850.2019.1684920
8. Washington P, Park N, Srivastava P, Voss C, Kline A, Varma M, et al. Data-driven diagnostics and the potential of mobile artificial intelligence for digital therapeutic phenotyping in computational psychiatry. Biol Psychiatry Cogn Neurosci Neuroimaging. (2019) 5:759–69. doi: 10.1016/j.bpsc.2019.11.015
9. Dawson G, Campbell K, Hashemi J, Lippmann SJ, Smith V, Carpenter K, et al. Atypical postural control can be detected via computer vision analysis in toddlers with autism spectrum disorder. Sci Rep. (2018) 8:17008. doi: 10.1038/s41598-018-35215-8
10. Grossard C, Dapogny A, Cohen D, Bernheim S, Juillet E, Hamel F, et al. Children with autism spectrum disorder produce more ambiguous and less socially meaningful facial expressions: an experimental study using random forest classifiers. Mol Autism. (2010) 11:5. doi: 10.1186/s13229-020-0312-2
11. Drimalla H, Scheffer T, Landwehr N, Baskow I, Roepke S, Behnia B, et al. Towards the automatic detection of social biomarkers in autism spectrum disorder: introducing the simulated interaction task (SIT). NPJ Digit Med. (2020) 3:25. doi: 10.1038/s41746-020-0227-5
12. Anzulewicz A, Sobota K, Delafield-Butt JT. Toward the autism motor signature: gesture patterns during smart tablet gameplay identify children with autism. Sci Rep. (2016) 6:31107. doi: 10.1038/srep31107
13. Ardalan A, Assadi AH, Surgent OJ, Travers BG. Whole-body movement during videogame play distinguishes youth with autism from youth with typical development. Sci Rep. (2019) 9:20094. doi: 10.1038/s41598-019-56362-6
14. Bachevalier J. Medial temporal lobe structures and autism: a review of clinical and experimental findings. Neuropsychologia. (1994) 32:627–48. doi: 10.1016/0028-3932(94)90025-6
15. Bachevalier J. Brief report: medial temporal lobe and autism: a putative animal model in primates. J Autism Dev Disord. (1996) 26:217–20. doi: 10.1007/BF02172015
16. Doi H, Fujisawa TX, Kanai C, Ohta H, Yokoi H, Iwanami A, et al. Recognition of facial expressions and prosodic cues with graded emotional intensities in adults with asperger syndrome. J Autism Dev Disord. (2013) 43:2099–113. doi: 10.1007/s10803-013-1760-8
17. Loth L, Garrido L, Ahmad J, Watson E, Duff A, Duchaine B. Facial expression recognition as a candidate marker for autism spectrum disorder: how frequent and severe are deficits? Mol Autism. (2018) 30:9:7. doi: 10.1186/s13229-018-0187-7
18. Baron-Cohen S, Ring HA, Bullmore ET, Wheelwright S, Ashwin C, Williams SCR. The amygdala theory of autism. Neurosci Biobehav Rev. (2000) 24:355–64. doi: 10.1016/S0149-7634(00)00011-7
19. Mazefsky CA, Herrington J, Siegel M, Scarpa A, Maddox BB, Scahill L, et al. The role of emotion regulation in autism spectrum disorder. J Am Acad Child Adolesc Psychiatry. (2013) 52:679–88. doi: 10.1016/j.jaac.2013.05.006
20. Kreibig SD. Autonomic nervous system activity in emotion: a review. Biol Psychol. (2010) 84:394–421. doi: 10.1016/j.biopsycho.2010.03.010
21. Smeekens I, Didden R, Verhoeven EW. Exploring the relationship of autonomic and endocrine activity with social functioning in adults with autism spectrum disorders. J Autism Dev Disord. (2015) 45:495–505. doi: 10.1007/s10803-013-1947-z
22. Jansen LM, Gispen-de Wied CC, Wiegant VM, Westenberg HG, Lahuis BE, van Engeland H. Autonomic and neuroendocrine responses to a psychosocial stressor in adults with autistic spectrum disorder. J Autism Dev Disord. (2006) 36:891–9. doi: 10.1007/s10803-006-0124-z
23. Watson LR, Roberts JE, Baranek GT, Mandulak KC, Dalton JC. Behavioral and physiological responses to child-directed speech of children with autism spectrum disorders or typical development. J Autism Dev Disord. (2012) 42:1616–29. doi: 10.1007/s10803-011-1401-z
24. Kushki A, Drumm E, Pla Mobarak M, Tanel N, Dupuis A, Chau T, et al. Investigating the autonomic nervous system response to anxiety in children with autism spectrum disorders. PLoS ONE. (2013) 8:e59730. doi: 10.1371/journal.pone.0059730
25. Louwerse A, Tulen JH, van der Geest JN, van der Ende J, Verhulst FC, Greaves-Lord K. Autonomic responses to social and nonsocial pictures in adolescents with autism spectrum disorder. Autism Res. (2014) 7:17–27. doi: 10.1002/aur.1327
26. Guy L, Souders M, Bradstreet L, DeLussey C, Herringto JD. Brief report: Emotion regulation and respiratory sinus arrhythmia in autism spectrum disorder. J Autism Dev Disord. (2014) 44:2614–20. doi: 10.1007/s10803-014-2124-8
27. Matsushima K, Matsubayashi J, Toichi M, Funabiki Y, Kato K, Awaya T, et al. Unusual sensory features are related to resting-state cardiac vagus nerve activity in autism spectrum disorders. Res Autism Spectr Disord. (2016) 25:37–46. doi: 10.1016/j.rasd.2015.12.006
28. Harder R, Malow BA, Goodpaster RL, Iqbal F, Halbower A, Goldman SE, et al. Heart rate variability during sleep in children with autism spectrum disorder. Clin Auton Res. (2016) 26:423–32. doi: 10.1007/s10286-016-0375-5
29. Schaffler MD, Middleton LJ, Abdus-Saboor I. Mechanisms of tactile sensory phenotypes in autism: current understanding and future directions for research. Curr Psychiatry Rep. (2019) 21:134. doi: 10.1007/s11920-019-1122-0
30. Poh MZ, McDuff DJ, Picard RW. Non-contact, automated cardiac pulse measurements using video imaging and blind source separation. Opt Express. (2010) 18:10762–74. doi: 10.1364/OE.18.010762
31. Fukunishi M, Kurita K, Yamamoto S, Tsumura N. Non-contact video based estimation of heart rate variability spectrogram from hemoglobin composition. Artif Life Robot. (2017) 22:457–63. doi: 10.1007/s10015-017-0382-1
32. Porges SW. The polyvagal perspective. Biol Psychol. (2007) 74:116–43. doi: 10.1016/j.biopsycho.2006.06.009
33. Shaffer F, Ginsberg JP. An overview of heart rate variability metrics and norms. Front Public Health. (2017) 5:258. doi: 10.3389/fpubh.2017.00258
34. Gross JJ, Levenson RW. Emotion elicitation using films. Cogn Emot. (1995) 9:87–108. doi: 10.1080/02699939508408966
35. Samson AC, Kreibig SD, Soderstrom B, Wade AA, Gross JJ. Eliciting positive, negative and mixed emotional states: a film library for affective scientists. Cognit Emot. (2016) 30:827–56. doi: 10.1080/02699931.2015.1031089
36. Werling DM, Geschwind DH. Sex differences in autism spectrum disorders. Curr Opin Neurol. (2013) 26:146–53. doi: 10.1097/WCO.0b013e32835ee548
37. Van Wijngaarden-Cremers PJM, van Eeten E, Groen WB, Van Deurzen PA, Oosterling IJ, Van der Gaag RJ. Gender and age differences in the core triad of impairments in autism spectrum disorders: a systematic review and meta-analysis. J Autism Dev Disord. (2014) 44:627–35. doi: 10.1007/s10803-013-1913-9
38. Parsons OE, Bayliss AP, Remington A. A few of my favorite things: circumscribed interests in autism are not accompanied by increased attentional salience on a personalized selective attention task. Mol Autism. (2017) 8:20. doi: 10.1186/s13229-017-0132-1
39. Brosnan M, Johnson H, Grawmeyer B, Chapman E, Benton L. Emotion recognition in animated compared to human stimuli in adolescents with autism spectrum disorder. J Autism Dev Disord. (2015) 45:1785–96. doi: 10.1007/s10803-014-2338-9
40. Himichi T, Osanai H, Goto T, Fujita H, Kawamura Y, Davis MH, et al. Development of a Japanese version of the interpersonal reactivity index. Japan J Psychol. (2017) 88:61–71. doi: 10.4992/jjpsy.88.15218
41. Komaki G, Maeda M, Arimura T, Nakata A, Shinoda H, Ogata I, et al. The reliability and factorial validity of the Japanese version of the 20-item toronto alexithymia scale (TAS-20). Jpn J Psychosom Med. (2003) 43:839–46. doi: 10.1016/S0022-3999(03)00360-X
42. Davis MH. Measuring individual differences in empathy: Evidence for a multidimensional approach. J. Pers. Soc. Psychol. (1983) 44:113–26. doi: 10.1037/0022-3514.44.1.113
43. Bagby RM, Taylor GJ, Parker JDA. The twenty-item toronto alexithymia scale: II. Convergent, discriminant, concurrent validity. J Psychosom Res. (1994) 38:33–40. doi: 10.1016/0022-3999(94)90006-X
44. Fukunishi M, Kurita K, Yamamoto S, Tsumura N. Video Based measurement of heart rate and heart rate variability spectrogram from estimated hemoglobin information. In: 2018 IEEE/CVF Conference on Computer Vision and Pattern Recognition Workshops (CVPRW). Salt Lake City, UT. (2018). doi: 10.1109/CVPRW.2018.00180
45. Tsumura N, Haneishi H, Miyake Y. Independent-component analysis of skin color image. J Opt Soc Am A Opt Image Sci Vis. (1999) 16:2169–76. doi: 10.1364/JOSAA.16.002169
46. Tarvainen MP, Niskanen JP, Lipponen JA, Ranta-Aho PO, Karjalainen PA. Kubios HRV–heart rate variability analysis software. Comput Methods Programs Biomed. (2014) 113:210–20. doi: 10.1016/j.cmpb.2013.07.024
47. Erbas Y, Ceulemans E, Boonen J, Noens I, Kuppens P. Emotion differentiation in autism spectrum disorder. Res Autism Spectr Disord. (2013) 7:1221–7. doi: 10.1016/j.rasd.2013.07.007
48. Cath DC, Ran N, Smit JH, van Balkom AJLM, Comijs HC. Symptom overlap between autism spectrum disorder, generalized social anxiety disorder and obsessive-compulsive disorder in adults: A preliminary case-controlled study. Psychopathology. (2008) 41:101–10. doi: 10.1159/000111555
49. Hartley CA, Phelps EA. Changing fear: the neurocircuitry of emotion regulation. Neuropsychopharmacology. (2010) 35:136–46. doi: 10.1038/npp.2009.121
50. Metallinou A, Grossman RB, Narayanan S. Quantifying atypicality in affective facial expression of children with autism spectrum disorders. Proc. (2013) 2013:1–6. doi: 10.1109/ICME.2013.6607640
51. Lenroot RK, Yeung PK. Heterogeneity within autism spectrum disorders: what have we learned from neuroimaging studies? Front Hum Neurosci. (2013) 7:733. doi: 10.3389/fnhum.2013.00733
52. Doi H, Fujisawa TX, Iwanaga R, Matsuzaki J, Kawasaki C, Tochigi M, et al. Association between single nucleotide polymorphisms in estrogen receptor 1/2 genes and symptomatic severity of autism spectrum disorder. Res Dev Disabil. (2018) 82:20–6. doi: 10.1016/j.ridd.2018.02.014
53. Asif M, Martiniano HFMC, Marques AR, Santos JX, Vilela J, Rasga C, et al. Identification of biological mechanisms underlying a multidimensional ASD phenotype using machine learning. Transl Psychiatry. (2020) 10:43. doi: 10.1038/s41398-020-0721-1
54. Kinnaird E, Stewart C, Tchanturia K. Investigating alexithymia in autism: a systematic review and meta-analysis. Eur Psychiatry. (2019) 55:80–9. doi: 10.1016/j.eurpsy.2018.09.004
55. Top DN, Stephenson KG, Doxey CR, Crowley MJ, Kirwan CB, South M. atypical amygdala response to fear conditioning in autism spectrum disorder. Biol Psychiatry Cogn Neurosci Neuroimaging. (2016) 1:308–15. doi: 10.1016/j.bpsc.2016.01.008
56. Hewig J, Hagemann D, Seifert J, Gollwitzer M, Naumann E, Bartussek D. A revised film set for the induction of basic emotions. Cogn Emot. (2005) 19:1095–109. doi: 10.1080/02699930541000084
57. Baio J, Wiggins L, Christensen DL, Maenner MJ, Daniels J, Warren Z, et al. Prevalence of autism spectrum disorder among children aged 8 years — Autism and developmental disabilities monitoring network, 11 sites, United States, 2014. MMWR Surveil Summaries. (2018) 67:1–23. doi: 10.15585/mmwr.ss6706a1
Keywords: ASD, autonomic nervous system, emotion, non-contact measurement, digital phenotyping, pulse wave, color
Citation: Doi H, Tsumura N, Kanai C, Masui K, Mitsuhashi R and Nagasawa T (2021) Automatic Classification of Adult Males With and Without Autism Spectrum Disorder by Non-contact Measurement of Autonomic Nervous System Activation. Front. Psychiatry 12:625978. doi: 10.3389/fpsyt.2021.625978
Received: 04 November 2020; Accepted: 01 April 2021;
Published: 17 May 2021.
Edited by:
Asma Bouden, Tunis El Manar University, TunisiaReviewed by:
Claudia Civai, London South Bank University, United KingdomPrabitha Urwyler, University of Bern, Switzerland
Copyright © 2021 Doi, Tsumura, Kanai, Masui, Mitsuhashi and Nagasawa. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hirokazu Doi, hdoi_0407@yahoo.co.jp