Next Article in Journal
Modeling and Performance Optimization of an Irreversible Two-Stage Combined Thermal Brownian Heat Engine
Next Article in Special Issue
Bispectral Analysis of Heart Rate Variability to Characterize and Help Diagnose Pediatric Sleep Apnea
Previous Article in Journal
Research on the Prediction of A-Share “High Stock Dividend” Phenomenon—A Feature Adaptive Improved Multi-Layers Ensemble Model
Previous Article in Special Issue
Short-Term HRV Analysis Using Nonparametric Sample Entropy for Obstructive Sleep Apnea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Complexity of Body Movements during Sleep in Children with Autism Spectrum Disorder

1
Department of Psychiatry and Neurobiology, Graduate School of Medical Science, Kanazawa University, Kanazawa 920-8640, Japan
2
Research Center for Child Mental Development, Kanazawa University, Kanazawa 920-8640, Japan
3
Department of Neuropsychiatry, University of Fukui, Fukui 910-1193, Japan
4
Uozu Shinkei Sanatorium, Uozu 937-0017, Japan
5
Department of Earth and Planetary Sciences, Graduate School of Sciences, Kyushu University, Fukuoka 819-0395, Japan
6
Institute of Human and Social Sciences, Kanazawa University, Kanazawa 921-1192, Japan
*
Authors to whom correspondence should be addressed.
Entropy 2021, 23(4), 418; https://doi.org/10.3390/e23040418
Submission received: 12 March 2021 / Revised: 30 March 2021 / Accepted: 30 March 2021 / Published: 31 March 2021
(This article belongs to the Special Issue Entropy and Sleep Disorders II)

Abstract

:
Recently, measuring the complexity of body movements during sleep has been proven as an objective biomarker of various psychiatric disorders. Although sleep problems are common in children with autism spectrum disorder (ASD) and might exacerbate ASD symptoms, their objectivity as a biomarker remains to be established. Therefore, details of body movement complexity during sleep as estimated by actigraphy were investigated in typically developing (TD) children and in children with ASD. Several complexity analyses were applied to raw and thresholded data of actigraphy from 17 TD children and 17 children with ASD. Determinism, irregularity and unpredictability, and long-range temporal correlation were examined respectively using the false nearest neighbor (FNN) algorithm, information-theoretic analyses, and detrended fluctuation analysis (DFA). Although the FNN algorithm did not reveal determinism in body movements, surrogate analyses identified the influence of nonlinear processes on the irregularity and long-range temporal correlation of body movements. Additionally, the irregularity and unpredictability of body movements measured by expanded sample entropy were significantly lower in ASD than in TD children up to two hours after sleep onset and at approximately six hours after sleep onset. This difference was found especially for the high-irregularity period. Through this study, we characterized details of the complexity of body movements during sleep and demonstrated the group difference of body movement complexity across TD children and children with ASD. Complexity analyses of body movements during sleep have provided valuable insights into sleep profiles. Body movement complexity might be useful as a biomarker for ASD.

1. Introduction

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by difficulties with social interaction and communication [1]. Many results of earlier studies have demonstrated that children with ASD commonly experience sleep problems [2,3,4], and that they are more likely to persist than in typically developing (TD) children [5,6,7,8,9]. Sleep disturbances reportedly correlate with core ASD symptoms [4,10,11,12,13,14,15,16] and other co-occurring symptoms and behaviors [3,10,15,16,17,18], such as anxiety [13,14,19], withdrawal [20,21], attention, hyperactivity, and aggression [12,14,15,22,23].
Polysomnography (PSG), regarded as the gold standard for sleep quality assessment, records many physiological parameters, including electroencephalography, electrooculography, electromyography, and electrocardiography, respiratory measures, and leg muscle activity. However, many ASD children with sensory abnormalities might not tolerate multiple electrodes on their body during sleep [24,25,26]. As an alternative, actigraphy, which measures body movement acceleration, has been used as a non-invasive and objective tool to assess the sleep–wake cycle. Recently, good agreement between actigraphy and PSG has been reported from sleep assessments of children with ASD [25]. Earlier studies of children with ASD have mainly found longer sleep latency, circadian rhythm disruption, lower sleep efficiency, and shorter total sleep time [2,3,4,17,27]. However, across individuals, children with ASD experience different sleep difficulties such as insomnia, circadian rhythm disorder, and parasomnia. Such difficulties might be attributable to various factors such as circadian-relevant gene anomalies, abnormal melatonin rhythms, brain wave organization, hyper-arousal, sensory hyper-reactivity, and various stresses, suggesting great heterogeneity of sleep problems in children with ASD [3,16]. Given that context, assessing the complex dynamics of body movements is expected to be useful in addition to assessing simple sleep parameters. In fact, the importance of nonlinear analyses and feature extraction have been reported from an actigraphy study [28]. Moreover, some earlier studies have shown significant alterations of movement complexity in patients with bipolar disorder [29,30,31,32,33], dementia [34,35], and insomnia [36].
Physiological systems are complex dynamical systems that integrate various mechanisms with various spatial and temporal scales. To date, the complexities of various physiological data such as EEG, ECG, respiration, pulse, and DNA sequences have been investigated. An exceedingly common complexity analysis of a one-dimensional time series is a deterministic approach that evaluates chaotic attractors. According to Takens’ embedding theorem, when a sequence consists of scalar measurement of a deterministic dynamical system [37], the time delay coordinate map with a sufficiently large dimension of the sequence provides a one-to-one image of the overall system behavior. Chaotic attractors are evaluated using two approaches with the delay coordinate: geometrically by fractal dimensions of the attractors and dynamically by divergence of nearby trajectories [38]. The latter one characterizes the sensitivity to initial conditions, which is an important aspect of complex systems. Another approach to characterize this sensitivity is the use of information-theoretic measures [39]. This type of nonlinear analysis, such as approximate entropy (ApEn) and sample entropy (SampEn), is a stochastic approach to quantify the diffusion of nearby trajectories. In other words, these measures evaluate signal variation by quantifying the irregularity and unpredictability of the time series. It is important to measure the entropy on multiple time scales to consider the hierarchical structures of physiological bases [40,41]. Compared to investigation of strange attractors, the information-theoretic measures are better suited for analyzing short, noisy, and non-stationary sequences, which are common in physiological data [42]. Physiological processes are non-stationary. Therefore, it is also meaningful to examine temporal changes in complexity. Our earlier work showed that our proposed expanded sample entropy (expSampEn), which quantifies the unpredictability of each time point, improved neural decoding accuracy [43]. In addition, other nonlinear analyses have specifically examined longer-range temporal relations. Fluctuation analysis quantifies the fractality of a cumulative summed sequence to ascertain the long-range temporal correlation in the original time series [44,45].
For this study, we applied several complexity analyses to data of body movements during sleep. This report is the first of a study characterizing the complexity of body movement in TD children and children with ASD. In summary, our aims were two-fold: (1) describe details of the complex dynamics of body movements, and (2) demonstrate group differences of complexity of body movements across TD children and children with ASD. To this end, we first examined the presence of determinism using the false nearest neighbor (FNN) method. Then, we evaluated sensitivity to initial conditions using ApEn, SampEn, and expSampEn, and evaluated temporal correlation using detrended fluctuation analysis (DFA).

2. Materials and Methods

2.1. Participants, Measurements, and Data Preprocessing

Information related to participants, measurements, and preprocessing of actigraphy has been reported by Naito et al. [46]. Participants in this cross-sectional study were 17 TD children with no reported behavioral or language problem and 17 children with ASD recruited from Kanazawa University and prefectural hospitals in Kanazawa or Toyama. Children in the TD and ASD groups were excluded from the study if a review of their medical history showed a history of epilepsy or intellectual disability or if they were taking psychotropic medications. Children with ASD were diagnosed using the Diagnostic and Statistical Manual of Mental Disorders (5th edition) (DSM-5) [1], the Japanese translation of the Diagnostic Interview for Social and Communication Disorders (DISCO) [47,48], or the Japanese translation of the Autism Diagnostic Observation Schedule (ADOS) [49,50,51]. All children fulfilled the diagnosis of childhood autism or atypical autism with DISCO and autism spectrum disorder with DSM-5. All except one also met the ADOS criteria for autism or autism spectrum disorder (Table 1). Children with psychiatric disorders other than ASD and ADHD were excluded using the Mini International Neuropsychiatric Interview for Children and Adolescents [52]. As presented in Table 2, no significant difference was found between TD and ASD groups in terms of gender, age, usual sleep duration, or sleep quality (rating scale: 1 = very bad, 2 = bad, 3 = somewhat bad, 4 = somewhat good, 5 = good, and 6 = very good) as reported by caregivers, or in cognitive skills as assessed by the Japanese adaptation of the Kaufman Assessment Battery for Children (K-ABC) [53,54]. Yuko Yoshimura, a speech language hearing therapist, and Chiaki Hasegawa, a psychologist, administered the psychological tests. The Ethics Committee of Kanazawa University Hospital approved the methods and procedures, which were performed in accordance with the Declaration of Helsinki. Parents agreed to their child’s participation in the study with full knowledge of the experimental characteristics of the research. Written informed consent was obtained after a complete explanation of this study, but before participation in the study.
Body movement was measured using a wristwatch-like accelerometer (Gen-2 GSR Wristband; Interuniversity Microelectronics Centre, Leuven, Kingdom of Belgium) attached to the waist [46]. The acceleration was measured and sampled at 32 Hz from immediately before the child entered the bed to immediately after the child exited the bed. Body movements were measured on at least three nights, excluding days when the child had a sickness such as a common cold.
The acceleration time series were preprocessed using software (Brain Vision Analyzer; Brain Products GmbH, Gilching, Germany) and Python. The three-dimensional (3D) time series recorded at a sampling rate of 32 Hz was resampled to 1 Hz by block averaging over 1 s, high-pass filtering at 0.0028 Hz in each dimension to exclude the effect of sustained gravitational acceleration, and conversion to a one-dimensional (1D) time series by calculating the root mean square (RMS) for each time point. We analyzed data of 8 h after the first time the acceleration became less than 0.1 G/s for 10 min.

2.2. Surrogate Analysis

To confirm the body movement nonlinearity, we also performed surrogate analysis using the iterative amplitude-adjusted Fourier transform (IAAFT) algorithm [55,56]. This algorithm is proposed for testing of the null hypothesis that the time series has been generated by a stationary linear stochastic process with Gaussian inputs that has possibly gone through a static monotonic transformation during measurement. This surrogate series has the same linear correlation and probability distribution as the original time series. Rejection of the null hypothesis indicates the existence of nonlinearity in the time series [38]. We generated 20 surrogates for each time series and examined their power spectral density (PSD) using Welch’s method (Figure 1). No significant difference was found for PSD among TD, ASD, and the surrogates.

2.3. Raw and Thresholded Data

Additionally, time series of two types were applied to the complexity analyses. In conventional actigraphs, the activity is usually recorded in each epoch by algorithms of three types: summing all acceleration values (digital integration, DI), computing the number of values above a threshold (time above threshold, TAT), or computing the number of zero crossings (ZC). However, raw data can also be used [57]. The DI and raw data include information about the magnitude of acceleration. For TAT and ZC, the information is converted to the iterations and durations of the activity, not the magnitude of acceleration. Therefore, we also used raw and thresholded data to analyze information of two types: magnitude of activity and appearance of activity. To exclude the influence of inter-subject and inter-day variation in the magnitude of body movements, we binarized body movements using the median value for each subject and each night. No significant difference of the threshold was found between TD children and children with ASD (Figure 2).

2.4. Complexity Analyses

2.4.1. Determinism Detection

To detect the determinism of body movement dynamics, we used the FNN algorithm with an appropriate time delay [38,58]. One can consider a time series = {xn, n = 1, …, N}, for which n represents each time point from 1 to N time points. As the first step in determinism detection, we computed 1/e decay time of the autocorrelation function (autocorrelation time), which is a common choice of appropriate time delay [38]. We used the median value of autocorrelation time as the time delay for the FNN method (Figure 3). Then, we consider d-dimensional delay vectors as:
y i d   =   x i ,   x i + τ ,   x i + 2 τ ,   ,   x i + d 1 τ
where d represents the embedding dimension and τ stands for the time delay. We computed the fraction of FNN for each d that satisfies the following two conditions:
x i + d τ x n i , d + d τ y i d y n i , d d   > A  
y i d + 1 y n i , d d + 1 σ   > B
Therein, n(i,d) stands for the time point of the nearest neighbor of yi(d), σ signifies the standard deviation (SD) of the scalar time series {xi}, ||∙|| represents the Euclidean norm, and A and B denote suitable thresholds. The fraction of i that satisfies (2) or (3) among all nearest neighbor combinations was defined as the fraction of FNN. If the time series is generated using a deterministic process, then it is expected that the chaotic attractor can be unfolded in the dimension where the fraction of FNN converges to zero. However, for a stochastic time series, it is expected that the fraction of FNN will never become zero. Herein, A and B were set as 10.0 and 2.0 [58]. The Theiler window was set to the same value as the delay, τ [38]. In addition, the inputs were moving-averaged over delay, τ, before counting the FNNs. The analyses in this section were performed using NoLiTSA, a Python module (https://github.com/manu-mannattil/nolitsa) assessed on 12 March 2021.

2.4.2. Information-Theoretic Analyses

To characterize the irregularity and unpredictability of body movement dynamics, we conducted entropy-based analyses: ApEn [59], SampEn [42], and expSampEn [43]. The following is a brief summary of these algorithms [43,60]. Consider a time series = x n ,   n = 1 , , N , in which n represents each time point and N denotes the total length. Defining x n = x n m x n 1 as the m-dimensional past vector, then x n = x n m x n is the (m + 1)-dimensional vector obtained by concatenating current value xn to x n . Also, p x n | x n is the conditional probability that the current value is close to x n when the m-dimensional past vector is close to x n . Also, p x n and p x n are the joint probabilities that m + 1 and m consecutive points in the time series are respectively close to x n and x n . Generally, points or vectors within r × SD of Chebyshev distance are regarded as “close”. Then, ApEn, SampEn, and expSampEn are represented as:
ApEn = log p x n | x n = log p x n p x n
SampEn = log p x n p x n
expSampEn   n   = log p x n | x n
where <∙> represents the average over time. Here, ApEn and SampEn are single values because they represent the irregularity of the entire time series, whereas expSampEn is a time series because it represents the irregularity of each time point with respect to the entire time series. In other words, the idea of expSampEn is based on the local or pointwise use of information theory [61,62]. In fact, the average of expSampEn over the entire time is equal to ApEn. As in our earlier work [43,63], we used m = 2 and r = 0.2 to calculate the entropy.
Additionally, the complex dynamics are expected to be described in greater detail by examining the overall picture of these information-theoretic measures over multiple time scales [40,41]. Therefore, we calculated ApEn, SampEn, and expSampEn for three time scales (30, 100, and 300 s). The inputs were moving-averaged over the time scales before calculating ApEn, SampEn, and expSampEn. The threshold of nearest neighbor, r × SD, was defined for the moving-averaged time series of each time scale. The analyses explained in this section were performed using EntroPy, a Python module (https://github.com/raphaelvallat/entropy) assessed on 12 March 2021.

2.4.3. Fluctuation Analysis

To examine the long-range temporal correlation of body movement dynamics, we performed fluctuation analysis that is able to quantify the fractality of a cumulative summed sequence. We applied DFA (detrended fluctuation analysis): a very common fluctuation analyses [44,45]. In this algorithm, the entire series is first split by size, n. Then, the root mean square of the deviation from the local trend is computed to show typical fluctuations of the series, F(n):
F n = 1 N k = 1 N y k y n k 2
In that equation, y is the cumulative summed sequence of x, and yn(k) is the local trend. F(n) usually exhibits a power law scaling F n   n α . Temporal correlation is evaluated by its exponent α. When α ≈ 1/2, x is temporally uncorrelated, and when α > 1/2, x is temporally correlated. When α ≈ 1, F(n) exhibits 1/f scaling, and x is temporally long-range correlated [36]. We examined the scaling properties of y between 5 min and 2 h.

2.5. Statistical Analyses

Student’s t-tests were applied to all two-group comparisons. ExpSampEn is a time series. Therefore, we did not correct for multiple comparisons in expSampEn. To ensure the same conditions as those used for expSampEn, we also did not correct for multiple comparisons in any other two-group comparison. Additionally, because of the large variation of expSampEn, a moving average for the 60 min preceding the t-tests was taken.
Spearman’s correlation coefficients were calculated to examine the correlation between the original time series (raw data and thresholded data) and expSampEn (Table 3). All time series were moving-averaged for the 60 min before calculating the correlation coefficients.
For results of all statistical analyses, p < 0.05 was inferred as indicating statistical significance.

3. Results

3.1. Stationarity of Body Movements Overnight

Initially, to test the stationarity of the body movement time series, we compared the distribution of the first half of the time series with the distribution of the entire time series [38,64]. We conducted a chi-square test using five bins for the raw data and two bins for the thresholded data. Results showed that all 34 subjects × 3 times of the raw data and 92.2% (except eight data) of the thresholded data were rejected for stationarity at the 0.05 level of significance. These results indicate that the distributions of both the magnitude and the presence of body movements are significantly different between the first and latter half of sleep, i.e., the magnitude and the presence of body movements were non-stationary overnight.

3.2. Deterministic Chaos

Next, to detect the deterministic chaos of body movement dynamics, we used the FNN algorithm. According to Takens’ embedding theorem, detecting determinism requires an appropriate time delay [38]. We applied median values of the autocorrelation time as the time delay (raw data, 116 s; thresholded data, 347 s; Figure 3, left panels). As shown in the right panels of Figure 3, the fraction of FNN did not converge to zero for either the raw and thresholded data, indicating that deterministic chaos was not detected using the FNN algorithm. However, this lack of convergence does not signify that the body movement is not a complex system. Rather, the chaos cannot be detected because the deterministic approach is vulnerable to noise and non-stationarity. Actually, stationarity was rejected for much of the time series, as described in Section 3.1. Therefore, we investigated body movement complexity using information theoretic measures and fluctuation analysis, which are robust to noise and non-stationarity, as described in the following sections.

3.3. Information-Theoretic Measures

We then applied information theoretic measures to quantify the irregularity and unpredictability associated with sensitivity to initial conditions [39]. First, we examined the irregularity of the time series of data for the entire night using ApEn and SampEn. No significant difference was found in ApEn or SampEn between TD children and children with ASD for either raw or thresholded data. However, on almost all scales, ApEn and SampEn were significantly lower in the original data than in the surrogates, for which linearity was assumed (Figure 4). These results demonstrated nonlinearity in the irregularity of body movements, but the irregularity of overnight body movements was comparable between TD children and children with ASD.
Next, we examined the irregularity and unpredictability at each time point using expSampEn, which is the local version of ApEn and SampEn [43,61]. We initially checked for temporal changes in irregularity of each subject of each night. Although both the expSampEn of raw data and of thresholded data showed phasic behaviors, the behaviors are apparently uncorrelated (Figure 5A). As described in Table 3, although the raw data and thresholded data were highly correlated (R = 0.88–0.89), expSampEn of the thresholded data showed weaker correlation with raw data (R = 0.14–0.25) and with expSampEn of raw data (R = 0.33–0.51). This result is apparently mainly attributable to the weak correlation between the thresholded data and their expSampEn (R = 0.23–0.46). Subsequently, we investigated group differences in expSampEn between TD and ASD (Figure 5B). Although almost no significant difference was found in the raw data, expSampEn in the thresholded data yielded significantly lower results for ASD than for TD for 300 s from sleep onset to 2 h after sleep onset, and for 30, 100, and 300 s at approximately 6 h after sleep onset. Additionally, for several time scales, expSampEn yielded significantly lower results for the original data than for the surrogates.
To evaluate group differences of irregularity, we compared the representative values of TD and ASD obtained for high- and low-irregularity states. The 90th and 10th percentiles of the moving-averaged expSampEn for thresholded data were defined respectively as the representative values for the high-irregularity and low-irregularity states because the degree of irregularity behaves phasically (Figure 5A), and because the timings of switching between the high-irregularity and low-irregularity states differ among subjects. The high-irregularity values of ASD on the time scales of 100 and 300 s were significantly lower than those of TD, although no significant difference was found for the low-irregularity value (Table 4).
In summary, children with ASD moved less irregularly than TD children in the early (up to 2 h) and late (after 5 h) periods of sleep. From a periodical view of irregularity, differences in irregularity between TD and ASD were observed at the high-irregularity state. Additionally, nonlinearity was involved in the irregularity of body movements.

3.4. DFA

To characterize the temporal correlation of body movement dynamics, we applied DFA to the cumulative summed sequence of raw and thresholded data (Figure 6). We confirmed power-law scaling for at least one order of time duration and evaluated exponent α by fitting it. No significant difference was found between exponent α of TD children and children with ASD, although exponent α of surrogates was significantly lower than that of TD children and children with ASD for thresholded data (p < 0.001). Exponent α of the original data (especially the thresholded data) was approximately 1.0, suggesting that fluctuation of the original data exhibits 1/f scaling. In other words, body movements during sleep are not temporally discrete but are rather temporally correlated. This temporal correlation was similar between TD children and children with ASD.

4. Discussion

Earlier reports of some studies have described high prevalence of sleep difficulties of various types across children with ASD [2,3] and have suggested the involvement of various neurophysiological mechanisms [3,16]. To characterize this diverse ASD-related sleep difficulty, assessment of body movement complexity during sleep might add another dimension to sleep difficulties already identified from conventional sleep assessments, such as sleep latency and total sleep time [24]. Many physiological phenomena can be described in terms of complex dynamics. Body movements are also expected to have complexity. In fact, complexity analyses have been applied to some earlier actigraphy studies [29,30,31,32,33,34,35,36,65,66]. Nevertheless, earlier studies have yielded no consistent findings across various analyses of complexity. This report is the first of a detailed study comprehensively investigating body movement complexity during sleep using different complexity measures and highlighting their relevance to ASD diagnosis. To examine the validity of complexity analyses, we first discuss the complex dynamics of body movements during sleep, irrespective of the psychiatric disorder. Subsequently, we characterized the ASD-related complexity profiles of body movements. Finally, the study limitations are listed.

4.1. Complex Dynamics in Body Movements

Initially, body movement complexity during sleep, irrespective of psychiatric disorders, is discussed. As explained earlier in Section 2.3, raw and thresholded data were used respectively, to examine the complexity profile of the magnitude and presence of body movement. Additionally, body movement nonlinearity was investigated using surrogates. In summary, although we were unable to detect determinism (Figure 3), results showed that nonlinear processes reduced the irregularity and unpredictability (Figure 4 and Figure 5B) and enhanced long-range temporal correlation (Figure 6) in both the magnitude and presence of body movement. Strong non-stationarity, as described in Section 3.1, and the influence of nonlinear processes were observed in body movements during sleep, suggesting that information-theoretic analyses and fluctuation analyses, which are robust to non-stationary data, are effective for investigating complex systems of body movement.
Regarding temporal correlation, earlier studies applying DFA analysis to actigraphy have also revealed the existence of 1/f fluctuations in body movements during daytime [67], night [36], and 24 h [29,34,35]. These findings suggest that body movement during sleep shows long-range temporal correlation and suggest that it can be modeled by self-similar stochastic processes such as fractional Brownian motion [68]. Additionally, although entropy-based measures cannot be compared directly because of their relativity [69], altered irregularity in body movement has been reported for various psychiatric disorders, as described specifically hereinafter. These observations also suggest that body movement complexity can be characterized by irregularity. Surrogate analyses in both information-theoretic measures and DFA suggest the influence of nonlinear processes in body movement. The only report describing a study that compares body movement time series with their surrogates using fluctuation analysis explains the existence of nonlinear processes in body movement fluctuations in healthy subjects [70]. Consideration of these findings together suggests that body movements can be modeled using a nonlinear process with self-similarity that is more complex than fractional Brownian motion.

4.2. Group Differences (TD vs. ASD) in Complexity of Body Movement

Results of group comparison of the body movement complexity demonstrated less irregularity in ASD than in TD, in a time-specific manner (Figure 5B). Earlier reports of some studies have described increased wake after sleep onset based on caregiver’s reports [8], PSG [71,72,73], and actigraphy [74], suggesting increased body movement after sleep onset in ASD. Use of the same dataset as that used for the present study revealed that the duration of body movements was longer in children with ASD over similar time periods to those examined for our study [46]. Considering these findings collectively, the behavior of increased body movement in ASD is characterized with regularity, specifically in the early and late stages of sleep. It is noteworthy that, despite the high correlation across original datasets (raw vs. thresholded data), complexity analysis applied to them indicated different and distinct characteristics (Table 3). This finding might be the reason for the low complexity of thresholded data in spite of the increased body movement in children with ASD. Additionally, when we defined high-irregularity and low-irregularity states to incorporate consideration of the phasic behavior of irregularity (Figure 5A) and compared the respective states for TD and ASD, reduced irregularity of body movement in ASD was found only in the high-irregularity state. This result suggests the importance of phase (e.g., sleep stage or gene expression) for elucidating the sleep characteristics of ASD.
Few reports of the relevant literature have described studies examining the complexity of body movements in psychiatric disorders. For example, irregularities assessed by entropy-based methods such as SampEn are reportedly higher for mania than for bipolar depression [30,33]. Another study of schizophrenia showed a disease-associated increase of body movement [65]. However, most actigraphy studies using these information-theoretic analyses have specifically examined the active periods. No report of the relevant literature has specifically examined the sleep period. Additionally, no report of the relevant literature has described an examination of the irregularity and unpredictability of body movements in children with ASD. Compared to entropy analyses, the DFA algorithm requires long time series data aimed at evaluating long-range temporal correlation. In the only reported study analyzing fluctuations during sleep, exponent α was found to increase, showing long-range temporal correlation for insomnia patients [36]. Other reports of the literature have described that exponent α decreases with age and in Alzheimer’s disease patients [34,67], and that scale-dependent changes in exponent α are found in cases of bipolar disorder [29]. Earlier actigraphy reports have described altered irregularity and unpredictability during active periods in bipolar disorder and schizophrenia, along with altered temporal correlation in insomnia, aging, Alzheimer’s disease, and bipolar disorder. Novel findings of the present study are that the appearance of body movements during sleep in ASD is less irregular than in TD, suggesting the usefulness of entropy-based analyses of body movement during sleep for the diagnosis of ASD. However, it seems contradictory that no disruption was found in the fractality of body movement during sleep in children with ASD, who experience many sleep difficulties, despite altered fractality in the fluctuation of body movements in insomnia patients. This finding might be explained by the heterogeneity of sleep problems in children with ASD [3,16], or by differences in the ages of subjects (18–65 years in the study of Holloway et al. [36], 5–8 years in this study). Therefore, additional studies must be conducted to assess the relevance of ASD itself and sleep problems in children with ASD to the complexity of body movements.

4.3. Limitations

Body movement complexity exhibited phasic behavior (Figure 5A) that resembles a periodic pattern in the sleep architecture. To date, no consistent conclusion has been reported for the sleep architecture of children with ASD [75,76]. For instance, reports of earlier studies have described both increased and decreased rapid eye movement (REM) sleep [20,26,71,77,78] and slow wave sleep (SWS) [20,26,73,78]. Moreover, alteration of the cyclic alternation pattern (CAP) [27] and increased undifferentiated sleep [26,77] have been reported. For the present study, we detected a difference in the certain state of the periodic pattern during sleep between children with ASD and TD children (Table 4). It is particularly interesting, however, that two aspects of complexity (magnitude of body movement (raw data) and presence of body movement (thresholded data)) oscillated differently (Table 3). This result suggests that body movement has at least two oscillatory components. Although accelerometry alone is not regarded as able to classify sleep stages because it is an indirect measure of sleep [79,80], a recent report has described that actigraphy combined with cardiorespiratory cues can classify wake/REM/non-REM with high accuracy [81]. However, because we did not measure PSG, the relation between REM/non-REM cycle and the two types of periodic behavior of irregularity we detected remains unknown. Additional studies using both actigraphy and PSG must be conducted to identify the meanings of the periodic behaviors of complexity.
Additionally, surrogate analyses suggest that body movement is a complex system generated by nonlinear processes rather than by linear stochastic processes. However, because the null hypothesis of IAAFT assumes a stationary linear stochastic process in a precise sense, the non-stationarity of body movements (see Section 3.1) might affect the surrogate behavior. Therefore, caution is necessary when interpreting the results of surrogate analyses.
For this study, the accelerometer was attached to the waist because the comfort of children was prioritized. Another report described that even 1.5-year-old TD children placed the accelerometer on their waists for 7 consecutive days [82]. However, some reports of earlier studies have described differences of actigraphy depending on the accelerometer attachment site [83,84,85]. Additional studies must be conducted to validate the actigraphy data obtained using waist-attached accelerometers.

5. Conclusions

As described herein, we were unable to prove the presence of determinism using the FNN algorithm. However, the irregularity and unpredictability evaluated using entropy-based analyses and the exponent α of fluctuation analysis showed significant differences between the original data and its surrogates, suggesting that body movements become less irregular and become temporally long-range correlated by nonlinear processes. Additionally, the irregularity of the magnitude and the presence of body movement behaved phasically, but differently from one another. Within these complexities, the irregularity and unpredictability of the presence of body movements was found to be significantly lower in children with ASD than in TD children, especially during high-irregularity periods. Further investigation using PSG will be needed to determine whether body movement complexity in children with ASD is reduced only during sleep, and if so, at which sleep stage.

Author Contributions

Conceptualization, Formal analysis, Software, Validation, Visualization, Writing-original draft, N.F.; Data curation and Resources, N.N., Y.Y., C.H., T.H. and M.K.; Funding acquisition, Project acquisition and Supervision, T.T. and M.K.; Investigation, N.F., N.N., Y.Y., C.H., T.H. and M.K.; Methodology, N.F., T.T. and T.M.; Writing-review and editing, N.F., T.T., T.M., Y.Y. and M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Center of Innovation Program and CREST (Grant Number JPMJCR17A4) from the Japan Science and Technology Agency (https://www.coistream.osaka-u.ac.jp/en).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and the protocol was approved by the Medical Ethics Committee of Kanazawa University (protocol code 2014-035 (1641), September 17, 2014).

Informed Consent Statement

Informed consent was obtained from the parents of all subjects involved in the study.

Data Availability Statement

The data underlying this article cannot be shared publicly because of a non-disclosure agreement. The data will be shared on reasonable request to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ADOS-GAutism Diagnostic Observation Schedule–Generic
ApEnapproximate entropy
ASDautism spectrum disorder
CAPcyclic alternation pattern
DFAdetrended fluctuation analysis
DIdigital integration
DISCODiagnostic Interview for Social and Communication Disorders
DSM-5Diagnostic and Statistical Manual of Mental Disorders (5th edition)
expSampEnexpanded sample entropy
FNNfalse nearest neighbor
IAAFTiterative amplitude-adjusted Fourier transform
K-ABCKaufman Assessment Battery for Children
PSDpower spectral density
PSGPolysomnography
REMrapid eye movement
RMSroot mean square
SampEnsample entropy
SDstandard deviation
SWSSlow-wave sleep
TATtime above threshold
TDtypically developing children
ZCnumber of zero crossings

References

  1. Association, A.P. Diagnostic and Statistical Manual of Mental Disorders, 5th ed; American Psychiatric Publishing: Washington, DC, USA, 2013. [Google Scholar]
  2. Elrod, M.G.; Hood, B.S. Sleep differences among children with autism spectrum disorders and typically developing peers: A Meta-analysis. J. Dev. Behav. Pediatr. 2015, 36, 166–177. [Google Scholar] [CrossRef] [PubMed]
  3. Souders, M.C.; Zavodny, S.; Eriksen, W.; Sinko, R.; Connell, J.; Kerns, C.; Schaaf, R.; Pinto-Martin, J. Sleep in Children with Autism Spectrum Disorder. Curr. Psychiatry Rep. 2017, 19, 34. [Google Scholar] [CrossRef] [PubMed]
  4. Carmassi, C.; Palagini, L.; Caruso, D.; Masci, I.; Nobili, L.; Vita, A.; Dell’Osso, L. Systematic review of sleep disturbances and circadian sleep desynchronization in autism spectrum disorder: Toward an integrative model of a self-reinforcing loop. Front. Psychiatry 2019, 10, 366. [Google Scholar] [CrossRef]
  5. Hodge, D.; Carollo, T.M.; Lewin, M.; Hoffman, C.D.; Sweeney, D.P. Sleep patterns in children with and without autism spectrum disorders: Developmental comparisons. Res. Dev. Disabil. 2014, 35, 1631–1638. [Google Scholar] [CrossRef] [PubMed]
  6. Goldman, S.E.; Richdale, A.L.; Clemons, T.; Malow, B.A. Parental sleep concerns in autism spectrum disorders: Variations from childhood to adolescence. J. Autism Dev. Disord. 2012, 42, 531–538. [Google Scholar] [CrossRef]
  7. May, T.; Cornish, K.; Conduit, R.; Rajaratnam, S.M.W.; Rinehart, N.J. Sleep in High-Functioning Children With Autism: Longitudinal Developmental Change and Associations With Behavior Problems. Behav. Sleep Med. 2015, 13, 2–18. [Google Scholar] [CrossRef]
  8. Humphreys, J.S.; Gringras, P.; Blair, P.S.; Scott, N.; Henderson, J.; Fleming, P.J.; Emond, A.M. Sleep patterns in children with autistic spectrum disorders: A prospective cohort study. Arch. Dis. Child. 2014, 99, 114–118. [Google Scholar] [CrossRef] [Green Version]
  9. Verhoeff, M.E.; Blanken, L.M.E.; Kocevska, D.; Mileva-Seitz, V.R.; Jaddoe, V.W.V.; White, T.; Verhulst, F.; Luijk, M.P.C.M.; Tiemeier, H. The bidirectional association between sleep problems and autism spectrum disorder: A population-based cohort study. Mol. Autism 2018, 9. [Google Scholar] [CrossRef] [Green Version]
  10. Veatch, O.J.; Sutcliffe, J.S.; Warren, Z.E.; Keenan, B.T.; Potter, M.H.; Malow, B.A. Shorter sleep duration is associated with social impairment and comorbidities in ASD. Autism Res. 2017, 10, 1221–1238. [Google Scholar] [CrossRef]
  11. Schreck, K.A.; Mulick, J.A.; Smith, A.F. Sleep problems as possible predictors of intensified symptoms of autism. Res. Dev. Disabil. 2004, 25, 57–66. [Google Scholar] [CrossRef] [PubMed]
  12. Mayes, S.D.; Calhoun, S.L. Variables related to sleep problems in children with autism. Res. Autism Spectr. Disord. 2009, 3, 931–941. [Google Scholar] [CrossRef]
  13. Hollway, J.A.; Aman, M.G.; Butter, E. Correlates and risk markers for sleep disturbance in participants of the autism treatment network. J. Autism Dev. Disord. 2013, 43, 2830–2843. [Google Scholar] [CrossRef]
  14. Bangerter, A.; Chatterjee, M.; Manyakov, N.V.; Ness, S.; Lewin, D.; Skalkin, A.; Boice, M.; Goodwin, M.S.; Dawson, G.; Hendren, R.; et al. Relationship Between Sleep and Behavior in Autism Spectrum Disorder: Exploring the Impact of Sleep Variability. Front. Neurosci. 2020, 14, 211. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Goldman, S.E.; McGrew, S.; Johnson, K.P.; Richdale, A.L.; Clemons, T.; Malow, B.A. Sleep is associated with problem behaviors in children and adolescents with Autism Spectrum Disorders. Res. Autism Spectr. Disord. 2011, 5, 1223–1229. [Google Scholar] [CrossRef]
  16. Cohen, S.; Conduit, R.; Lockley, S.W.; Rajaratnam, S.M.; Cornish, K.M. The relationship between sleep and behavior in autism spectrum disorder (ASD): A review. J. Neurodev. Disord. 2014, 6, 44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Goldman, S.E.; Surdyka, K.; Cuevas, R.; Adkins, K.; Wang, L.; Malow, B.A. Defining the sleep phenotype in children with autism. Dev. Neuropsychol. 2009, 34, 560–573. [Google Scholar] [CrossRef] [Green Version]
  18. Lindor, E.; Sivaratnam, C.; May, T.; Stefanac, N.; Howells, K.; Rinehart, N. Problem behavior in autism spectrum disorder: Considering core symptom severity and accompanying sleep disturbance. Front. Psychiatry 2019, 10, 487. [Google Scholar] [CrossRef] [Green Version]
  19. Mazurek, M.O.; Petroski, G.F. Sleep problems in children with autism spectrum disorder: Examining the contributions of sensory over-responsivity and anxiety. Sleep Med. 2015, 16, 270–279. [Google Scholar] [CrossRef]
  20. Aathira, R.; Gulati, S.; Tripathi, M.; Shukla, G.; Chakrabarty, B.; Sapra, S.; Dang, N.; Gupta, A.; Kabra, M.; Pandey, R.M. Prevalence of Sleep Abnormalities in Indian Children With Autism Spectrum Disorder: A Cross-Sectional Study. Pediatr. Neurol. 2017, 74, 62–67. [Google Scholar] [CrossRef]
  21. Yavuz-Kodat, E.; Reynaud, E.; Geoffray, M.-M.; Limousin, N.; Franco, P.; Bonnet-Brilhault, F.; Bourgin, P.; Schroder, C.M. Disturbances of Continuous Sleep and Circadian Rhythms Account for Behavioral Difficulties in Children with Autism Spectrum Disorder. J. Clin. Med. 2020, 9, 1978. [Google Scholar] [CrossRef]
  22. Mazurek, M.O.; Sohl, K. Sleep and Behavioral Problems in Children with Autism Spectrum Disorder. J. Autism Dev. Disord. 2016, 46, 1906–1915. [Google Scholar] [CrossRef] [PubMed]
  23. Allik, H.; Larsson, J.O.; Smedje, H. Insomnia in school-age children with Asperger syndrome or high-functioning autism. BMC Psychiatry 2006, 6, 18. [Google Scholar] [CrossRef] [Green Version]
  24. Moore, M.; Evans, V.; Hanvey, G.; Johnson, C. Assessment of Sleep in Children with Autism Spectrum Disorder. Children 2017, 4, 72. [Google Scholar] [CrossRef] [PubMed]
  25. Yavuz-Kodat, E.; Reynaud, E.; Geoffray, M.-M.; Limousin, N.; Franco, P.; Bourgin, P.; Schroder, C.M. Validity of Actigraphy Compared to Polysomnography for Sleep Assessment in Children With Autism Spectrum Disorder. Front. Psychiatry 2019, 10, 551. [Google Scholar] [CrossRef] [Green Version]
  26. Malow, B.A.; Marzec, M.L.; McGrew, S.G.; Wang, L.; Henderson, L.M.; Stone, W.L. Characterizing sleep in children with autism spectrum disorders: A multidimensional approach. Sleep 2006, 29, 1563–1571. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Miano, S.; Bruni, O.; Elia, M.; Trovato, A.; Smerieri, A.; Verrillo, E.; Roccella, M.; Terzano, M.G.; Ferri, R. Sleep in children with autistic spectrum disorder: A questionnaire and polysomnographic study. Sleep Med. 2007, 9, 64–70. [Google Scholar] [CrossRef]
  28. Fuster-Garcia, E.; Bresó, A.; Miranda, J.M.; García-Gómez, J.M. Actigraphy pattern analysis for outpatient monitoring. Methods Mol. Biol. 2015, 1246, 3–17. [Google Scholar] [CrossRef]
  29. Knapen, S.E.; Li, P.; Der Lek, R.F.R.-V.; Verkooijen, S.; Boks, M.P.M.; Schoevers, R.A.; Scheer, F.A.J.L.; Hu, K. Fractal biomarker of activity in patients with bipolar disorder. Psychol. Med. 2020, 1–8. [Google Scholar] [CrossRef]
  30. Krane-Gartiser, K.; Henriksen, T.E.G.; Morken, G.; Vaaler, A.; Fasmer, O.B. Actigraphic assessment of motor activity in acutely admitted inpatients with bipolar disorder. PLoS ONE 2014, 9, e89574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Krane-Gartiser, K.; Asheim, A.; Fasmer, O.B.; Morken, G.; Vaaler, A.E.; Scott, J. Actigraphy as an objective intra-individual marker of activity patterns in acute-phase bipolar disorder: A case series. Int. J. Bipolar Disord. 2018, 6, 8. [Google Scholar] [CrossRef] [PubMed]
  32. Cuesta-Frau, D.; Schneider, J.; Bakštein, E.; Vostatek, P.; Spaniel, F.; Novák, D. Classification of actigraphy records from bipolar disorder patients using slope entropy: A feasibility study. Entropy 2020, 22, 1243. [Google Scholar] [CrossRef]
  33. Scott, J.; Vaaler, A.E.; Fasmer, O.B.; Morken, G.; Krane-Gartiser, K. A pilot study to determine whether combinations of objectively measured activity parameters can be used to differentiate between mixed states, mania, and bipolar depression. Int. J. Bipolar Disord. 2017, 5, 5. [Google Scholar] [CrossRef] [Green Version]
  34. Hu, K.; Van Someren, E.J.W.; Shea, S.A.; Scheer, F.A.J.L. Reduction of scale invariance of activity fluctuations with aging and Alzheimer’s disease: Involvement of the circadian pacemaker. Proc. Natl. Acad. Sci. USA 2009, 106, 2490–2494. [Google Scholar] [CrossRef] [Green Version]
  35. Huber, S.E.; Sachse, P.; Mauracher, A.; Marksteiner, J.; Pohl, W.; Weiss, E.M.; Canazei, M. Assessment of Fractal Characteristics of Locomotor Activity of Geriatric In-Patients With Alzheimer’s Dementia. Front. Aging Neurosci. 2019, 11, 272. [Google Scholar] [CrossRef]
  36. Holloway, P.M.; Angelova, M.; Lombardo, S.; St. Clair Gibson, A.; Lee, D.; Ellis, J. Complexity analysis of sleep and alterations with insomnia based on non-invasive techniques. J. R. Soc. Interface 2014, 11, 20131112. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Sauer, T.; Yorke, J.A.; Casdagli, M. Embedology. J. Stat. Phys. 1991, 65, 579–616. [Google Scholar] [CrossRef]
  38. Mannattil, M.; Gupta, H.; Chakraborty, S. Revisiting Evidence of Chaos in X-ray Light Curves: The Case of GRS 1915+105. Astrophys. J. 2016, 833, 208. [Google Scholar] [CrossRef] [Green Version]
  39. Burggren, W.W.; Monticino, M.G. Assessing physiological complexity. J. Exp. Biol. 2005, 208, 3221–3232. [Google Scholar] [CrossRef] [Green Version]
  40. Costa, M.; Goldberger, A.L.; Peng, C.-K. Multiscale Entropy Analysis of Complex Physiologic Time Series. Phys. Rev. Lett. 2002, 89, 068102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Costa, M.; Goldberger, A.L.; Peng, C.-K. Multiscale entropy analysis of biological signals. Phys. Rev. E 2005, 71, 021906. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Richman, J.S.; Moorman, J.R. Physiological time-series analysis using approximate and sample entropy. Am. J. Physiol. Hear. Circ. Physiol. 2000, 278, H2039–H2049. [Google Scholar] [CrossRef] [Green Version]
  43. Furutani, N.; Nariya, Y.; Takahashi, T.; Ito, H.; Yoshimura, Y.; Hiraishi, H.; Hasegawa, C.; Ikeda, T.; Kikuchi, M. Neural Decoding of Multi-Modal Imagery Behavior Focusing on Temporal Complexity. Front. Psychiatry 2020, 11, 746. [Google Scholar] [CrossRef] [PubMed]
  44. Peng, C.K.; Buldyrev, S.V.; Goldberger, A.L.; Havlin, S.; Sciortino, F.; Simons, M.; Stanley, H.E. Long-range correlations in nucleotide sequences. Nature 1992, 356, 168–170. [Google Scholar] [CrossRef]
  45. Peng, C.K.; Havlin, S.; Stanley, H.E.; Goldberger, A.L. Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series. Chaos 1995, 5, 82–87. [Google Scholar] [CrossRef]
  46. Naito, N.; Kikuchi, M.; Yoshimura, Y.; Kumazaki, H.; Kitagawa, S.; Ikeda, T.; Hasegawa, C.; Saito, D.N.; Tomiyama, S.; Minabe, Y. Atypical body movements during night in young children with autism spectrum disorder: A pilot study. Sci. Rep. 2019, 9, 6999. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Wing, L.; Leekam, S.R.; Libby, S.J.; Gould, J.; Larcombe, M. The Diagnostic Interview for Social and Communication Disorders: Background, inter-rater reliability and clinical use. J. Child Psychol. Psychiatry Allied Discip. 2002, 43, 307–325. [Google Scholar] [CrossRef] [Green Version]
  48. Wing, L. Diagnostic Interview for Social and Communication Disorders 11th edition. (T. Uchiyama et al., Trans.); Spectrum Publishing: Tokyo, Japan, 2007. [Google Scholar]
  49. Lord, C.; Risi, S.; Lambrecht, L.; Cook, E.H.; Leventhal, B.L.; Dilavore, P.C.; Pickles, A.; Rutter, M. The Autism Diagnostic Observation Schedule-Generic: A standard measure of social and communication deficits associated with the spectrum of autism. J. Autism Dev. Disord. 2000, 30, 205–223. [Google Scholar] [CrossRef] [PubMed]
  50. Lord, C.; Rutter, M.; Dilavore, P.C.; Risi, S.; Gotham, K.; Bishop, S.L. Autism Diagnostic Observation Schedule, 2nd ed.; Western Psychological Services: Torrance, CA, USA, 2012. [Google Scholar]
  51. Lord, C.; Rutter, M.; Dilavore, P.C.; Risi, S.; Gotham, K.; Bishop, S.L.; Luyster, R.J.; Guthrie, W. Autism Diagnostic Observation Schedule Second Edition. (M. Kuroda et al., Trans.); Kanekoshobo: Tokyo, Japan, 2015. [Google Scholar]
  52. Sheehan, D.V.; Sheehan, K.H.; Shytle, R.D.; Janavs, J.; Bannon, Y.; Rogers, J.E.; Milo, K.M.; Stock, S.L.; Wilkinson, B. Reliability and validity of the mini international neuropsychiatric interview for children and adolescents (MINI-KID). J. Clin. Psychiatry 2010, 71, 313–326. [Google Scholar] [CrossRef] [PubMed]
  53. Kaufman, A.S. K-ABC: Kaufman Assessment Battery for Children: Interpretive Manual; American Guidance Service: Circle Pines, MN, USA, 1983. [Google Scholar]
  54. Matsubara, T.; Fujita, K.; Maekawa, H.; Ishikuma, T.; Kaufman, A.S.; Kaufman, N.L. Interpretive Manual for the Japanese K-ABC.; Maruzen Mates: Tokyo, Japan, 1994. [Google Scholar]
  55. Schreiber, T.; Schmitz, A. Improved surrogate data for nonlinearity tests. Phys. Rev. Lett. 1996, 77, 635–638. [Google Scholar] [CrossRef] [Green Version]
  56. Schreiber, T.; Schmitz, A. Surrogate time series. Phys. D Nonlinear Phenom. 2000, 142, 346–382. [Google Scholar] [CrossRef] [Green Version]
  57. Lüdtke, S.; Hermann, W.; Kirste, T.; Beneš, H.; Teipel, S. An algorithm for actigraphy-based sleep/wake scoring: Comparison with polysomnography. Clin. Neurophysiol. 2021, 132, 137–145. [Google Scholar] [CrossRef]
  58. Kennel, M.B.; Brown, R.; Abarbanel, H.D.I. Determining embedding dimension for phase-space reconstruction using a geometrical construction. Phys. Rev. A 1992, 45, 3403–3411. [Google Scholar] [CrossRef] [Green Version]
  59. Pincus, S. Approximate entropy (ApEn) as a complexity measure. Chaos 1995, 5, 110–117. [Google Scholar] [CrossRef]
  60. Porta, A.; Bari, V.; De Maria, B.; Cairo, B.; Vaini, E.; Malacarne, M.; Pagani, M.; Lucini, D. On the Relevance of Computing a Local Version of Sample Entropy in Cardiovascular Control Analysis. IEEE Trans. Biomed. Eng. 2019, 66, 623–631. [Google Scholar] [CrossRef]
  61. Lizier, J.T. Measuring the Dynamics of Information Processing on a Local Scale in Time and Space. In Directed Information Measures in Neuroscience; Springer: Berlin/Heidelberg, Germany, 2014; pp. 161–193. [Google Scholar]
  62. Martinez-Cancino, R.; Heng, J.; Delorme, A.; Kreutz-Delgado, K.; Sotero, R.C.; Makeig, S. Measuring transient phase-amplitude coupling using local mutual information. Neuroimage 2018, 185, 361–378. [Google Scholar] [CrossRef]
  63. Furutani, N.; Nariya, Y.; Takahashi, T.; Noto, S.; Yang, A.C.; Hirosawa, T.; Kameya, M.; Minabe, Y.; Kikuchi, M. Decomposed Temporal Complexity Analysis of Neural Oscillations and Machine Learning Applied to Alzheimer’s Disease Diagnosis. Front. Psychiatry 2020, 11, 531801. [Google Scholar] [CrossRef]
  64. Isliker, H.; Kurths, J. A test for stationarity: Finding parts in time series APT for correlation dimension estimates. Int. J. Bifurc. Chaos 1993, 03, 1573–1579. [Google Scholar] [CrossRef]
  65. Osipov, M.; Behzadi, Y.; Kane, J.M.; Petrides, G.; Clifford, G.D. Objective identification and analysis of physiological and behavioral signs of schizophrenia. J. Ment. Health 2015, 24, 276–282. [Google Scholar] [CrossRef] [Green Version]
  66. Hauge, E.R.; Berle, J.Ø.; Oedegaard, K.J.; Holsten, F.; Fasmer, O.B. Nonlinear analysis of motor activity shows differences between schizophrenia and depression: A study using fourier analysis and sample entropy. PLoS ONE 2011, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Raichlen, D.A.; Klimentidis, Y.C.; Hsu, C.H.; Alexander, G.E.; Newman, A. Fractal Complexity of Daily Physical Activity Patterns Differs with Age over the Life Span and Is Associated with Mortality in Older Adults. J. Gerontol. Ser. A Biol. Sci. Med. Sci. 2019, 74, 1461–1467. [Google Scholar] [CrossRef]
  68. Mandelbrot, B.B.; Van Ness, J.W. Fractional Brownian Motions, Fractional Noises and Applications. SIAM Rev. 1968, 10, 422–437. [Google Scholar] [CrossRef]
  69. Pincus, S.M.; Gladstone, I.M.; Ehrenkranz, R.A. A regularity statistic for medical data analysis. J. Clin. Monit. 1991, 7, 335–345. [Google Scholar] [CrossRef]
  70. Hu, K.; Ivanov, P.C.; Chen, Z.; Hilton, M.F.; Stanley, H.E.; Shea, S.A. Non-random fluctuations and multi-scale dynamics regulation of human activity. Phys. A Stat. Mech. Appl. 2004, 337, 307–318. [Google Scholar] [CrossRef] [Green Version]
  71. Maski, K.; Holbrook, H.; Manoach, D.; Hanson, E.; Kapur, K.; Stickgold, R. Sleep Dependent Memory Consolidation in Children with Autism Spectrum Disorder. Sleep 2015, 38, 1955–1963. [Google Scholar] [CrossRef] [Green Version]
  72. Lázár, A.S.; Lázár, Z.I.; Bíró, A.; Gyori, M.; Tárnok, Z.; Prekop, C.; Keszei, A.; Stefanik, K.; Gádoros, J.; Halász, P.; et al. Reduced fronto-cortical brain connectivity during NREM sleep in Asperger syndrome: An EEG spectral and phase coherence study. Clin. Neurophysiol. 2010, 121, 1844–1854. [Google Scholar] [CrossRef] [PubMed]
  73. Limoges, É.; Mottron, L.; Bolduc, C.; Berthiaume, C.; Godbout, R. Atypical sleep architecture and the autism phenotype. Brain 2005, 128, 1049–1061. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Morgan, B.; Nageye, F.; Masi, G.; Cortese, S. Sleep in adults with Autism Spectrum Disorder: A systematic review and meta-analysis of subjective and objective studies. Sleep Med. 2020, 65, 113–120. [Google Scholar] [CrossRef] [PubMed]
  75. Cortesi, F.; Giannotti, F.; Ivanenko, A.; Johnson, K. Sleep in children with autistic spectrum disorder. Sleep Med. 2010, 11, 659–664. [Google Scholar] [CrossRef]
  76. Devnani, P.A.; Hegde, A.U. Autism and sleep disorders. J. Pediatr. Neurosci. 2015, 10, 304–307. [Google Scholar] [CrossRef] [Green Version]
  77. Elia, M.; Ferri, R.; Musumeci, S.A.; Del Gracco, S.; Bottitta, M.; Scuderi, C.; Miano, G.; Panerai, S.; Bertrand, T.; Grubar, J.C. Sleep in subjects with autistic disorder: A neurophysiological and psychological study. Brain Dev. 2000, 22, 88–92. [Google Scholar] [CrossRef]
  78. Buckley, A.W.; Rodriguez, A.J.; Jennison, K.; Buckley, J.; Thurm, A.; Sato, S.; Swedo, S. Rapid eye movement sleep percentage in children with autism compared with children with developmental delay and typical development. Arch. Pediatr. Adolesc. Med. 2010, 164, 1032–1037. [Google Scholar] [CrossRef] [PubMed]
  79. Chinoy, E.D.; Cuellar, J.A.; Huwa, K.E.; Jameson, J.T.; Watson, C.H.; Bessman, S.C.; Hirsch, D.A.; Cooper, A.D.; Drummond, S.P.A.; Markwald, R.R. Performance of Seven Consumer Sleep-Tracking Devices Compared with Polysomnography. Sleep 2020. [Google Scholar] [CrossRef] [PubMed]
  80. Sundararajan, K.; Georgievska, S.; te Lindert, B.H.W.; Gehrman, P.R.; Ramautar, J.; Mazzotti, D.R.; Sabia, S.; Weedon, M.N.; van Someren, E.J.W.; Ridder, L.; et al. Sleep classification from wrist-worn accelerometer data using random forests. Sci. Rep. 2021, 11. [Google Scholar] [CrossRef]
  81. Kalkbrenner, C.; Brucher, R.; Kesztyüs, T.; Eichenlaub, M.; Rottbauer, W.; Scharnbeck, D. Automated sleep stage classification based on tracheal body sound and actigraphy. GMS Ger. Med. Sci. 2019, 17. [Google Scholar] [CrossRef]
  82. Nakagawa, M.; Ohta, H.; Nagaoki, Y.; Shimabukuro, R.; Asaka, Y.; Takahashi, N.; Nakazawa, T.; Kaneshi, Y.; Morioka, K.; Oishi, Y.; et al. Daytime nap controls toddlers’ nighttime sleep. Sci. Rep. 2016, 6. [Google Scholar] [CrossRef] [Green Version]
  83. Adkins, K.W.; Goldman, S.E.; Fawkes, D.; Surdyka, K.; Wang, L.; Song, Y.; Malow, B.A. A Pilot Study of Shoulder Placement for Actigraphy in Children. Behav. Sleep Med. 2012, 10, 138–147. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Littner, M.; Kushida, C.A.; Anderson, W.M.D.; Bailey, D.; Berry, R.B.; Davila, D.G.; Hirshkowitz, M.; Kapen, S.; Kramer, M.; Loube, D.; et al. Practice parameters for the role of actigraphy in the study of sleep and circadian rhythms: An update for 2002. Sleep 2003, 26, 337–341. [Google Scholar] [CrossRef]
  85. Van Hilten, J.J.; Middelkoop, H.A.M.; Kuiper, S.I.R.; Kramer, C.G.S.; Roos, R.A.C. Where to record motor activity: An evaluation of commonly used sites of placement for activity monitors. Electroencephalogr. Clin. Neurophysiol. Evoked Potentials 1993, 89, 359–362. [Google Scholar] [CrossRef]
Figure 1. Surrogate data generated by the iterative amplitude-adjusted Fourier transform (IAAFT) algorithm: examples of original time series (left) and its surrogate (middle). The mean power spectral density (PSD) of TD, ASD, and surrogates (right). No significant difference of PSD was found among TD (blue), ASD (orange), and surrogates (gray). TD, typically developing children; ASD, children with autism spectrum disorder.
Figure 1. Surrogate data generated by the iterative amplitude-adjusted Fourier transform (IAAFT) algorithm: examples of original time series (left) and its surrogate (middle). The mean power spectral density (PSD) of TD, ASD, and surrogates (right). No significant difference of PSD was found among TD (blue), ASD (orange), and surrogates (gray). TD, typically developing children; ASD, children with autism spectrum disorder.
Entropy 23 00418 g001
Figure 2. Thresholds of TD children and children with ASD. The median value of each time series was used as the threshold. No significant difference was found between TD and ASD. Error bars represent the standard deviation of the thresholds.
Figure 2. Thresholds of TD children and children with ASD. The median value of each time series was used as the threshold. No significant difference was found between TD and ASD. Error bars represent the standard deviation of the thresholds.
Entropy 23 00418 g002
Figure 3. Time delay and false nearest neighbor (FNN). Upper and lower panels respectively show those of raw data and thresholded data. Left panels—Median values of autocorrelation time were set to the time delay for the FNN method (raw data, 116 s; thresholded data, 347 s). Orange lines, boxes, and whiskers respectively, indicate the median, percentiles 25 and 75, and minimum and maximum of autocorrelation time of all subjects. Right panels—Tests I and II respectively signify FNN for Equations (2) and (3) (see Section 2.4.1). The shaded areas represent the standard deviation of the fraction of FNN. For neither raw nor thresholded data did the fraction of FNN converge to zero.
Figure 3. Time delay and false nearest neighbor (FNN). Upper and lower panels respectively show those of raw data and thresholded data. Left panels—Median values of autocorrelation time were set to the time delay for the FNN method (raw data, 116 s; thresholded data, 347 s). Orange lines, boxes, and whiskers respectively, indicate the median, percentiles 25 and 75, and minimum and maximum of autocorrelation time of all subjects. Right panels—Tests I and II respectively signify FNN for Equations (2) and (3) (see Section 2.4.1). The shaded areas represent the standard deviation of the fraction of FNN. For neither raw nor thresholded data did the fraction of FNN converge to zero.
Entropy 23 00418 g003
Figure 4. ApEn and SampEn for raw and thresholded data. Time scales are 30, 100, and 300 s. No significant difference between TD children and children with ASD was found for either raw or thresholded data on any time scale, although significant differences were found between original and surrogate data. ApEn, approximate entropy; SampEn, sample entropy; n.s., not significant (p ≥ 0.05); *, p < 0.001.
Figure 4. ApEn and SampEn for raw and thresholded data. Time scales are 30, 100, and 300 s. No significant difference between TD children and children with ASD was found for either raw or thresholded data on any time scale, although significant differences were found between original and surrogate data. ApEn, approximate entropy; SampEn, sample entropy; n.s., not significant (p ≥ 0.05); *, p < 0.001.
Entropy 23 00418 g004
Figure 5. expSampEn time series are shown (moving-averaged over 60 min). (A) Typical examples of the expSampEn time series of raw and thresholded data. Both are measured on the same day for the same subject. (B) Left and right panels respectively show expSampEn for raw data and thresholded data. Time scales are 30 (upper panels), 100 (middle panels), and 300 s (lower panels). The solid line and the shaded area respectively show the mean and the standard error. Green shading represents that the expSampEn is significantly lower in ASD than in TD (p < 0.05).
Figure 5. expSampEn time series are shown (moving-averaged over 60 min). (A) Typical examples of the expSampEn time series of raw and thresholded data. Both are measured on the same day for the same subject. (B) Left and right panels respectively show expSampEn for raw data and thresholded data. Time scales are 30 (upper panels), 100 (middle panels), and 300 s (lower panels). The solid line and the shaded area respectively show the mean and the standard error. Green shading represents that the expSampEn is significantly lower in ASD than in TD (p < 0.05).
Entropy 23 00418 g005
Figure 6. DFA for raw and thresholded data. The log–log plot shows the mean of F(τ). The bar plot shows the mean ± SD of exponent α. Exponent α of surrogates were significantly lower than that of TD and ASD children for thresholded data. n.s., not significant (p ≥ 0.05); * p < 0.001.
Figure 6. DFA for raw and thresholded data. The log–log plot shows the mean of F(τ). The bar plot shows the mean ± SD of exponent α. Exponent α of surrogates were significantly lower than that of TD and ASD children for thresholded data. n.s., not significant (p ≥ 0.05); * p < 0.001.
Entropy 23 00418 g006
Table 1. Autism Diagnostic Observation Schedule (ADOS) scores of children with ASD (N = 17).
Table 1. Autism Diagnostic Observation Schedule (ADOS) scores of children with ASD (N = 17).
Score (Mean ± SD)ASD Cut-Off
ADOS-G, Module 1 (N = 2)
Communication4.5 ± 0.72
Reciprocal Social Interaction8.0 ± 4.24
Communication + Social Interaction12.5 ± 3.57
Play2.5 ± 2.1
Stereotyped Behaviors and Restricted Interests1.5 ± 0.7
ADOS-G, Module 2 (N = 6)
Communication4.7 ± 2.23
Reciprocal Social Interaction8.8 ± 1.94
Communication + Social Interaction13.5 ± 3.98
Imagination/Creativity1.0 ± 0.6
Stereotyped Behaviors and Restricted Interests1.5 ± 0.8
ADOS-G, Module 3 (N = 1)
Communication42
Reciprocal Social Interaction94
Communication + Social Interaction137
Imagination/Creativity0
Stereotyped Behaviors and Restricted Interests1
ADOS-2, Module 2 (N = 6)
Social Affect9.0 ± 1.4
Restricted and Repetitive Behavior1.5 ± 1.2
Total score10.5 ± 2.28
ADOS-2, Module 3 (N = 2)
Social Affect5.5 ± 2.1
Restricted and Repetitive Behavior1.0 ± 0.0
Total score6.5 ± 2.17
ADOS-G, Autism Diagnostic Observation Schedule–Generic; ADOS-2, Autism Diagnostic Observation Schedule–2.
Table 3. Spearman’s correlation coefficient between original data and expSampEn (mean ± SD).
Table 3. Spearman’s correlation coefficient between original data and expSampEn (mean ± SD).
Data Type–Data TypeTime Scale
30 s100 s300 s
Raw–Thr0.88 ± 0.050.88 ± 0.050.89 ± 0.06
Raw–Raw_expSampEn0.89 ± 0.050.78 ± 0.060.41 ± 0.13
Thr–Thr_expSampEn0.46 ± 0.150.46 ± 0.160.23 ± 0.17
Raw–Thr_expSampEn0.22 ± 0.210.25 ± 0.200.14 ± 0.20
Raw_expSampEn–Thr_expSampEn0.33 ± 0.160.49 ± 0.140.51 ± 0.14
Raw, raw data; Thr, thresholded data. All time series were moving-averaged over 60 min before calculating the correlation coefficients.
Table 4. t-test results for TD and ASD by high and low representative values of expSampEn.
Table 4. t-test results for TD and ASD by high and low representative values of expSampEn.
Time Scale90th Percentile10th Percentile
TpTp
30 s1.520.0690.970.17
100 s2.170.0190.910.18
300 s1.780.0431.460.077
Table 2. Demographic characteristics.
Table 2. Demographic characteristics.
TDASDp-value
Number of participants1717
Gender (male/female)11/613/4n.s.
Age in months, mean (range)71.1 (61−79)77.1 (60−98)n.s.
Usual sleep duration (mean ± SD h)9.45 ± 0.549.51 ± 0.59n.s.
Usual sleep quality 1 (mean ± SD h)5.0 ± 0.715.0 ± 0.87n.s.
K-ABC Mental Processing Scale (mean ± SD)102.8 ± 10.593.9 ± 18.9n.s.
1 Sleep quality was rated as 1 = very bad, 2 = bad, 3 = somewhat bad, 4 = somewhat good, 5 = good, and 6 = very good. K-ABC, Kaufman Assessment Battery for Children; n.s., not significant, TD = typically developing children, ASD = Autism spectrum disorder
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Furutani, N.; Takahashi, T.; Naito, N.; Maruishi, T.; Yoshimura, Y.; Hasegawa, C.; Hirosawa, T.; Kikuchi, M. Complexity of Body Movements during Sleep in Children with Autism Spectrum Disorder. Entropy 2021, 23, 418. https://doi.org/10.3390/e23040418

AMA Style

Furutani N, Takahashi T, Naito N, Maruishi T, Yoshimura Y, Hasegawa C, Hirosawa T, Kikuchi M. Complexity of Body Movements during Sleep in Children with Autism Spectrum Disorder. Entropy. 2021; 23(4):418. https://doi.org/10.3390/e23040418

Chicago/Turabian Style

Furutani, Naoki, Tetsuya Takahashi, Nobushige Naito, Takafumi Maruishi, Yuko Yoshimura, Chiaki Hasegawa, Tetsu Hirosawa, and Mitsuru Kikuchi. 2021. "Complexity of Body Movements during Sleep in Children with Autism Spectrum Disorder" Entropy 23, no. 4: 418. https://doi.org/10.3390/e23040418

APA Style

Furutani, N., Takahashi, T., Naito, N., Maruishi, T., Yoshimura, Y., Hasegawa, C., Hirosawa, T., & Kikuchi, M. (2021). Complexity of Body Movements during Sleep in Children with Autism Spectrum Disorder. Entropy, 23(4), 418. https://doi.org/10.3390/e23040418

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop