1. Introduction
Typical neurodevelopment is highly dynamic and variable, with transient changes that tend to occur asynchronously across a given age group. Starting at early infancy, some babies may grow very fast and develop mature neuromotor control before others of the same age group. Other babies may in contrast, be stunted in their growth and/or neuromotor control development [
1]. These disparate processes span non-uniform trajectories that vary highly from person to person, across the general population. Despite these known features, developmental research tends to assume an overall consistency of the developmental paths characterizing aging trajectories, whereby linear, parametric models use normal distributions and enforce precisely expected developmental milestones. When faced with a neurological condition like Autism Spectrum Disorders (ASD), the high heterogeneity of any given random draw of this population, then poses a challenge to science, since the assumptions of normality and linearity fail in ways that we do not fully understand.
Open access data can provide information to address these questions on neurodevelopment, so elusive today to basic science. In autism, we know very little about the structure of the probability distribution landscape spanned by the parameters that we can derive from the nervous systems’ self-generated biorhythms. There is a paucity of work with limited sample size addressing the stochastic shifts and maturation trends across the population that is aging with autism [
2,
3,
4]. Consequently, we have virtually no knowledge about typical trends in the maturation of the nervous systems in large cross-sections of the aging population; let alone on the population aging with autism spectrum disorders (ASD.) Analytical tools in ASD research are largely dominated by a “one size fits all” static model that does not yet consider age-appropriate shifts in probability spaces derived from neuromotor variability.
There is an increasing need to stratify these heterogeneous data generated by clinical diagnoses rooted in observation of behavior [
3]. Some efforts have been recently appreciated in the literature [
5,
6,
7,
8,
9]. Nevertheless, the recent efforts have yet to consider age-appropriate shifts in statistical landscape and provide a framework that encourages doing so in relation to normative data, also absent from all clinical diagnostics tools that inform and largely steer today basic research in autism [
10]. The exclusion of the motor axis from the current inventories by the Diagnostics Statistical Manual DSM-5 from the American Psychiatric Association APA [
11] and by other tools from psychology (e.g., ADOS [
10,
12]) makes it difficult to use movements and their kinesthetic sensations as aids to objectively categorize different forms of autism as they evolve with aging. Large open access repositories offer a new avenue to complement clinical criteria with physical biometrics derived from variations in human biorhythms from the nervous systems.
In recent years, the micro-movements spikes (MMS) fluctuations inherently present in and underlying all natural behaviors that clinicians observe and describe using pencil and paper means, have been assessed across a proposed taxonomy of neuromotor control [
13] (
Figure 1). This phylogenetically orderly taxonomy proposes different levels of maturation in neuromotor control under typical neurodevelopment. Because e.g., autonomic function matures earlier than voluntary control, and likewise, reflexes evolve in early infancy, etc. the model also proposes that studying autism across these levels may help us stratify the heterogeneity of autism across aging. Over years of peer reviewed work, using this model has served to reveal several stochastic features of voluntary [
2,
14], spontaneous [
15], involuntary [
15,
16,
17,
18] and autonomic [
19,
20] motions, offering classification power in autism and other neuropsychiatric and neurological conditions (e.g., schizophrenia [
21], Parkinson’s disease [
22,
23,
24,
25], neuropathies [
26] and impairments of the nervous systems due to traumatic brain injury inducing coma [
27] or stroke [
28]).
In autism, voluntary MMS have provided evidence for a maturational process absent in the autistic phenotype and served as aid to (cross-sectionally) classify different developmental stages according to levels of motor noise evolution and probability distribution shape well-correlated with measures of IQ [
2,
4]. Since we always sample neurotypical controls and characterize their signatures first, we can always go back to this taxonomy and measure the departure of autistics from the neurotypical signatures. In this paper, we focus on the involuntary levels of neuromotor control, because involuntary motions contribute to the peripheral feedback that the central nervous system continuously senses as the person behaves. When in excess, or when random and noisy, this source of self-generated neurofeedback may interfere with overall neuromotor control and autonomy of the brain over the body in motion. Yet, because these micro-motions are so subtle, they may scape the naked eye of the observant clinician. We have previously found excess motor noise in autism, but have yet to characterize it as a function of aging [
16].
Indeed, this recent body of work points at the use of MMS criteria as a potential tool to track the evolution of some of the biorhythms of the nervous systems as the neurotypical person ages. We could then use the normative data to build a proper notion of stochastic change, in a probabilistic sense. This would enable us to better understand the aging process of the autistic nervous systems in cross-sections of the aging population, beyond developmental periods and define age-dependent biometrics capturing different rates of change in young highly plastic systems vs. in aging systems whereby neurodegenerative processes manifest. A pressing question in any lifelong neurological condition like ASD is, when does the plasticity of a young system phase-transitions into a neurodegenerative process? Forecasting this transition would allow us to provide neuroprotective therapies in the rapidly emerging ASD aging population.
In this paper, we present new analytical methods to automatically, in a heuristics-free manner, stratify patterns of variability inherently present in the involuntary head motions of neurotypical controls and of participants with ASD. Their involuntary micro-movements spikes data are reconstructed from brain images, registered as they try to volitionally control their body to prevent it from moving. As they deliberately attempt to dampen movements’ fluctuations upon being instructed to remain still in the fMRI setting, their involuntary motion data reveal the probability landscape of undesirable motor noise across sections of the ASD developing children and the ASD aging population. We explore data drawn at random from over 2144 records of the ABIDE repository and report on fundamental departures of the person aging with autism from the typically aging person.
2. Materials and Methods
2.1. Demographics and Boot Strapping Technique to Form Age Groups of Equal Size
Across all sites of ABIDE, informed consent was obtained from all participants and/or their legal guardians. The following method have been published previously in various papers [
16,
17,
18], but we provide them here as well to facilitate the flow of the paper. All datasets included in this study are from the Autism Brain Imaging Data Exchange (ABIDE) databases:
ABIDE obeys the following guideline on the use of human subject’s data: “In accordance with HIPAA guidelines and 1000 Functional Connectomes Project/INDI protocols, all datasets have been anonymized, with no protected health information included.”
The study includes two main comparisons:
Autism Spectrum Disorder (ASD), and Typical Development (TD), using estimation of stochastic signatures of involuntary head micro-movements of individuals with a formal DSM-ASD [
11] diagnosis of ASD and TD controls.
Ranges of age. Each group (ASD and TD) was split in seven different groups according to their age to assess how the stochastic signature of involuntary head micro-movements evolves with aging. The ranges of age used to that end were the following: from 5 to 10 years old, from 11 to 15 years old, from 16 to 20 years old, from 21 to 25 years old, from 26 to 30 years old, from 31 to 40 years old and from 41 to 65 years old.
Inclusion/Exclusion Criteria
This study includes all sites publicly available through ABIDE I and ABIDE II. They were comprised of 1127 TD and 1017 ASD. As we explained above, those groups were divided by age.
Table 1 provides the number of participants with ASD or TD are in each range of age in ABIDE dataset.
Bootstrapping Method
ABIDE has data deposited from various labs whereby two types of data are available: cleaned up data sets (whereby noisy images have been removed according to some threshold of head motion amplitude) and raw data sets (whereby no images have been removed). In this work we use the raw data sets with no removal, because we are precisely interested in the head motion data. To that end, we aim at examining all frames continuously acquired in each scan session. Further, we note that different scanner sampling rates influence the nature of the noise in the images, which in turn is reflected in the noise derived from the variability of the raw extracted from the images reflecting the patterns of head motions. Our lab previously reported these effects on the linear (mm/s) and angular (rad/s) speed data [
18]. Here we focus on the standardized MMS waveform (see explanation below). We analyze the patters of variability in standardized deviations of the raw signal relative to the empirically estimated mean amplitude, considering randomization across all possible samples. To that end, we use bootstrapping and randomize across our original sample sets.
Furthermore, given the inconsistent group sizes per age groups extracted from the ABIDE datasets (see
Table 1), we used a bootstrapping method to ensure uniform group numbers for pairwise statistical comparisons across ages. To that end, we used random sampling with replacement and created subgroups drawn from the original size group while considering the minimum number n = 25 at a time. The time series waveforms from these 25 randomly chosen participants were pooled into a representative point of a given age group. We chose 25 because the smallest age’s sub-groups size was n = 30. Thus, after dividing the groups by age, we did extract 100 and 500 random representative points (as described) to form sub-groups of 100 and 500 group sizes from all the age’s sub-groups. The five-year increments were motivated by prior work highlighting medication use and clinical characterizations of this ABIDE ASD cohort. The standardized amplitude using the MMS waveform ensured shuffling of data across all types of sampling resolution in ABIDE and accounted for anatomical disparities across the different age groups [
30]. In this paper we focus on the 500 size groups, but consistent trends were found with 100 size groups.
2.2. Data Processing
2.2.1. Motion Extraction
Head motion patterns were extracted from imaging data during (rs) fMRI experiments. Motion extraction was performed using the Analysis of Functional NeuroImages (AFNI) software packages [
31]. Single-subject processing scripts were generated using the afni_proc.py interface [
31]. Skull stripping was performed on anatomical data and functional EPI data were co-registered to anatomical images. The median was used as the EPI base in alignment. Motion parameters: three translational (
x,
y, and
z) and three rotational (pitch-about the x axis, roll-about the y axis, and yaw- about the
z axis), from EPI time-series registration was saved.
2.2.2. Statistical Analyses
We assess the scan-by-scan speed-dependent variations in the linear displacement of the head during resting-state functional magnetic resonance imaging (rs-fMRI) sessions. The analyses specifically refer to the stochastic signatures of micro-movement spikes (the MMS, defined in prior peer reviewed work [
2,
16,
17,
18]) and their treatment as a continuous random process, under the general rubric of Poison random processes. Specifically, we use a Gamma process and estimate the continuous family of probability distributions using maximum likelihood estimation (MLE) with 95% confidence.
2.2.3. Micro-Movement Spikes Data Type
The maximum amplitude of the speed (mm/s) was obtained (
Figure 2A shows the raw waveforms and
Figure 2B the MMS, which is the normalized version of the raw peaks). The empirically estimated mean speed of each trial was also obtained and used as reference to determine the maximal absolute amplitude deviations from it.
The time-series of these fluctuations in maximal amplitude deviations from the empirically estimated mean provides the waveform of interest for our analyses. To remove allometric effects of body-size across ages in each trial we computed the normalized peak amplitude (the peak speed amplitude is divided by the sum of the peak speed amplitude and the averaged speed amplitude value comprising points between the two speed minima surrounding the local speed peak amplitude) [
32,
33]. These normalized peaks spanning real values in the [0,1] interval are the spike trains of random fluctuations in signal amplitude (speed in this case) (
Figure 2C). They are assumed to characterize a continuous random process, where events in the past may (or may not) accumulate evidence towards the prediction of future events. The normalized fluctuations define the MMS of the waveform [
2] (
Figure 2C).
2.3. Distance in Probability Space: The Earth Mover’s Distance
The Earth Mover’s Distance (EMD) [
34,
35], also known as the Kantarovich-Wasserstein distance [
36], measures the distance between two discrete probability distributions. Given two discrete distributions
Pr and
Pθ, each with m possible states x or y, respectively; the EMD computes how much mass needs to be moved how far, to turn one distribution into the other.
Figure 2D shows two sample distributions derived from the linear speed that we obtained from head involuntary motions of two participants, one TD and one ASD, while
Figure 3A shows the matrix of pairwise EMD quantities across all 500 TD participants in the 5–10-year-old group.
The problem of finding the amount of work that it takes to turn one distribution into another has been posed as an optimization problem that can be resolved using linear programming [
35,
37] in O(n
2log(n)) vs. a more recent version in O(n) [
38]. Briefly, if one considers these histograms as two piles of earth, there are infinite ways to move the earth from one pile to the other. The goal is to find the optimal one (a unique value despite non-unique ways to do it). These are known as transport problems, with transport plan
γ(
x,y), to distribute the amount of earth from one place x over the domain of y (or vice versa.)
Here the constraints
and
apply to ensure that the plan yields the correct distributions (sum over rows and columns respectively). The
γ is the joint probability distribution such that
and
is the set of all distributions whose marginals are
Pr and
Pθ respectively. To obtain the EMD, every value of the matrix
is multiplied with the Euclidean distance between x and y:
where
is the Frobenius inner product (the sum of all the element-wise products) and inf is the infimum (minimizing the expression on the right). This represents an expectation with respect to the joint distribution γ where the outcome is unique but the flow from the minimization process is not unique. The product can be physically interpreted as mass times distance, (work), so Equation (1) is thought of as minimizing the work to transform one distribution into another [
35].
We use the EMD to compute pairwise for each age group of each cohort, the cost of turning one probability distribution into the other.
2.4. Stochastic Analyses
The foundation of these methods has been explained elsewhere (e.g., [
2,
3,
16,
17,
18,
39] among others). The empirically estimated mean was obtained using the continuous Gamma family of probability distributions for every group because it gave the best fit according to maximum likelihood estimation, MLE). Then, the MMS served as input to the Gamma process and a stochastic characterization of their fluctuations in amplitude was used to characterize the signature of involuntary head motions in the ASD vs. typically developing TD groups.
We examined the frequency histograms of the MMS waveform and use MLE to approximate the best fitting distribution encompassing all cases. To that end, we compared different families of probability distributions (e.g., the Gaussian, Normal, Lognormal, Exponential and Gamma) and chose the best fit in an MLE sense. Owing to our prior works using the ABIDE sets [
16,
17] we could determine that the Gamma had the best fit in an MLE sense. As such, we settled on the continuous Gamma family of probability distributions [
40]. The estimated parameters were plotted on a Gamma parameter plane, where the
x-axis represents the shape parameter value and the
y-axis represents the scale parameter value.
Figure 2E shows the representation of the estimated Gamma parameters (shape and scale) on the Gamma parameter plane with 95% confidence intervals.
Figure 2F shows the empirically estimated probability density functions (PDFs), while
Figure 2G shows their representation as points on a Gamma moments multidimensional parameter space. There, the mean, variance and skewness are represented as three dimensional points along the x, y, z axes, respectively. The fourth moment, kurtosis, is proportional to the size of the marker (higher kurtosis is represented by larger marker size) and the color conveys other information. For example, here we use the color to represent the participant type (ASD vs. TD) and the arrow marks the shift in ASD signature relative to the normative value. The color can also provide information on pairwise distance in some other space. We examined the cross-sectional trajectories representing the evolution of these parameters using canonical Cartesian coordinates and Euclidean distance to quantify change and their rates, along with the cumulative distances traveled by each cohort, and other kinematic parameters that such a representation of the probability data affords us. We also represented the EMD as explained above.
2.5. Network Connectivity Analyses
The MMS derived from the linear and angular speed were used to group the normalized peak amplitudes into a frequency histogram for each participant. Then, the EMD [
34,
35,
37] was obtained pairwise across all 500 records to estimate the cost of transforming one probability distribution into another. An adjacency matrix derived from this pairwise distance quantities (
Figure 3A) was used to obtain an undirected weighted network representation of each group, where the EMD was the weight of the links at entry (
i,j) connecting two records from participants
i and
j. Common network connectivity metrics (e.g., the clustering coefficient, the network’s node-to-node distances and the characteristic pathlength) were used to characterize the neurotypical cohort spanning the 7 age groups with similar number of participants (500 selected at random through bootstrapping, see above) and to provide the normative states of the network. Then, network states generated from the data of the ASD cohort were compared against the normative states.
The distributions of EMD values was also examined (
Figure 3B) to ascertain if the values spanned a multimodal distribution or a unimodal one. The Hartigan’s dip test of unimodality was used to that end [
41]. Then, the count per bin and the edges information was used to plot the count curve (as in
Figure 3C). The curve thus obtained was smoothed using loess algorithm set to 0.1 value using the smooth function in MATLAB version R2018b. The loess algorithm was used to preserve the peaks and valleys. Other smoothing procedures failed to preserve the original local minima and local maxima locations of the original histogram. To automatically separate self-emerging clusters of values using the local minima and local maxima spanned by the smooth count curve, the
findpeaks function in MATLAB was used.
Figure 3D, E show one example for the group of neurotypical controls between 31–40 years of age where the procedure was applied to select the peaks and valleys. The peak amplitudes (local maxima) were located and the values of the EMD that they comprised obtained (see arrows delimiting the clusters bounded by local minima values surrounding the local maximum value.) Then all local maxima and minima were pooled, to determine the full range of values across all age groups and both cohorts. This set of EMD values from the entire data set (7000 points from 500 participants in each of the 7 groups of the neurotypical controls and of the ASD) was sorted and used to build a color code for each participant. To that end, each
i-participant’s EMD value from the
i,j matrix entry was mapped (as in
Figure 3D–E) to the corresponding distribution bump and the corresponding local maxima of the bump was used to color the point representing the participant in the scatter of probability points plotted on the Gamma moments parameter space.
Figure 3F locates on the normalized color scale the value corresponding to the local EMD maximum for participant
i (matrix entry
i,j) taken across all other participants. Then, the point located on the Gamma moments parameter space is colored using this scale.
The idea was to uncover self-emerging clusters of EMD values in each age-group. Low values represent distributions that are similar, while high values represent distributions that are far apart in probability space. Possible scenarios are to have within each age group, highly heterogeneous values, or highly homogeneous values. Further, across the full scatter representing all age groups, we may have highly heterogeneous values, or homogeneous values. We set to interrogate this probability landscape about the evolution of the age groups, the heterogeneity/homogeneity of the groups and the nature of the cross-sectional trajectories of the groups.
To interrogate about trajectories, we used the mean values of the three first Gamma moments in each cohort and traced a line connecting each representative point of each age group for each cohort. We also did this for the mode of the age-group scatter, representing the most frequent value and for the median (not shown). The coordinates of the mean values along each axis, were set as the representative point of the scatter and tracked over time. The departure of the coordinates of the point representing the ASD age-group from the neurotypical control age group were obtained using the norm of the vector difference and plotted to ascertain the trends of the trajectories and to identify inflection points of the evolving curve representing the point by point distance between the two cohorts. Lastly, the trajectories were plotted in the scatter of points representing the probability landscape to ascertain their locations within the scatters.
A measure of the scatter dispersion was built using the Delaunay Triangulation (DT) of each age-group scatter and the area of each of the triangles obtained. The fluctuations in the values of the area were treated as a point process and their MMS obtained as well. The values of the DT-area MMS were used to empirically fit the Gamma distribution (given that this was the best distribution using MLE [
21]) and the estimated shape and scale parameters were plotted on the Gamma parameter plane with 95% confidence intervals for the empirically estimated shape and scale parameters. We used the estimated PDFs to ascertain differences of statistical significance across all age groups of the cohorts and within each age group for the ASD and neurotypical cases.
4. Discussion
This work aimed at using a distribution-free approach to characterize the cross-sectional aging trajectories of the stochastic signatures of involuntary micro-movements in neurotypicals from 5–65 years of age. The work further aimed at quantifying the departure from the normative data of age-matched participants with ASD. To that end, instead of applying a grand average method which would eliminate as noise the important motor fluctuations from the involuntary motions, we examined the MMS in two different spaces. One space was the canonical representation of points located on a Gamma moments parameter space, where differences are measured using Euclidean distance. The other was the probability space, where we use the EMD to measure similarities/differences between points. Each point represents the probability distribution of a person’s MMS derived from involuntary linear displacements of the head, registered as they try to remain still. We used a Gamma process to characterize the spike data. Together, at each age group, these points represent a family of continuous probability distribution functions. This family was well fit by the continuous Gamma family using Maximum Likelihood Estimation (MLE) with 95% confidence [
2].
We found that each age group could be well characterized by a different family of probability distributions, with different scatter dispersion. The methods provide a new way to automatically stratify the autistic cohort using the EMD on the probability distributions generated by the MMS, and to characterize the differences from the neurotypical cohort by representing the data as weighted connected graphs. The network analyses provided several connectivity metrics amenable to automatically separate atypical autistic development from neurotypical development. More importantly, we were able to automatically identify critical differences in the aging trajectories of the autistic cohort that were absent in the neurotypical cohort for the eldest group (41–65 years of age).
These results are relevant to autism research in more than one way. They emphasize the importance of not enforcing a priori any statistical assumptions, but rather allowing the variability inherently present in the data to automatically reveal self-emerging patterns. Further, the results show that a “one size fits all” approach is inefficient to cope with the heterogeneity of autism. Each age group in the autistic cohort had a different stochastic signature unique to the group, spanning a family of PDFs. These were distinguishable from other age groups within the ASD cohort. This was also the case for the normative data, which was nonetheless, overall, visibly more homogeneous in nature than the ASD set. At an individual level, each person in each age-group of the normative cohort was unambiguously statistically different from every other person. This was captured by the large values of the EMD, signaling high effort to transform one frequency histogram into any other frequency histogram of the cohort.
In stark contrast to the normative trends, the ASD group had non-uniform changes in the first 20 years of life. On the Gamma moments parameter space (
Figure 8A) a representative 11–15-year-old ASD child has stochastic signatures closer to a representative 5–10-year-old neurotypical child, and then accelerates to the 16–20 range, which remains behind normative values. The irregular trends in the acceleration and deceleration period separating the cohorts are visible in
Figure 8B,C across ages. As the two cross-sectional cohorts age, these differences increase, suggesting fundamental neurological disparity in the involuntary motions. These results are congruent with prior cross-sectional data from goal-directed voluntary motions, where the dramatically different trajectory in the maturation of the somatic-sensory-motor systems had been revealed and are well characterized by a power law [
2,
3] and by a parameter space defined by the R-metric [
4]. Here, using new methods, we extend the results from voluntary to involuntary motions across a much larger range of ages and larger cohorts that open-access data can afford us.
Beyond the MMS related work in ASD, the broader literature of motor control problems studying movements and their sensations in ASD has by now grown considerably [
42,
43], including problems with gait disturbances and vestibular dysfunction from an early age [
44,
45]. The current clinical diagnostics criteria do not include these problems as core symptoms. Yet, for well over two decades now, basic science has provided objective, physical quantification of neurological problems in ASD that supports an earlier neurological model [
46,
47] amenable to explain many of the social and communication issues defining the condition. Patterns of cerebellar dysfunction [
48,
49], motor and vestibular coordination [
44] and sensory-motor integration issues [
2] have been reproduced across tractable biorhythms using non-invasive means [
50,
51,
52,
53,
54,
55,
56]. Perhaps these data can help us derive further criteria to stratify ASD and help design treatments tailored to address the phenotypic features that each subgroup may reveal. More importantly, using the neurological model in combination with the current behavioral descriptors may be effective to tackle these somatic sensory motor issues from an early age.
Under the current behavioral criteria for diagnosis and treatments of ASD, these neurological issues are treated as comorbid and secondary to the observed behaviors defining the condition. Under such criteria, the adult autistic population, which is rather large today, never received neuroprotective therapies aimed to slow down the progression of neurological issues. Yet, Parkinsonism is more prevalent and has earlier onset in autistic adults than in neurotypicals of comparable age [
57]. Other medical conditions are also more prevalent in autistic adults than in controls of similar age [
58]. According to a recent report by US-Kaiser, medical issues include all major psychiatric disorders such as depression, anxiety, bipolar disorder, obsessive-compulsive disorder, schizophrenia, and suicide attempts. Further the report cites immune conditions, gastrointestinal and sleep disorders, seizure, obesity, dyslipidemia, hypertension, and diabetes with higher prevalence in adults with autism, along with stroke and Parkinson′s disease.
The results presented here strongly suggest that the involuntary micro-movements can not only provide a way to characterize the lack of volitional control of the brain over the body (i.e., control the body at will to remain still). More importantly, the results offer for the first time the means to differentiate aging with autism from typically aging in this cross section of the population. The methods allowed us to automatically stratify the broad spectrum of autism as a lifelong condition. The cross-sectional aging trajectories of these stochastic signatures pinpoint critical points of change in the statistical landscape of development and aging. These may perhaps forecast important milestones signaling physiological changes in the autistic system.
Along those lines, a striking feature of the data concerns the separation of the 41–65-year-old cluster from the rest of the scatter. It may be possible that the high penetrance of Fragile X Syndrome in autism accounts for such differentiation, as motor symptoms are common among this part of the population and often result in Parkinsonism after 40 years of age [
59,
60,
61]. It will be important to address this question, as the dramatic change in the involuntary motions for this ASD age group raises concerns about the large number of aging adults with ASD who may face motor issues impeding independent living.
Limitations
We note that among the autistic participants of ABIDE, there are those who have taken psychotropic medications and those who have not. We have previously shown that the involuntary MMS are present in both types of participants, but that as the number of meds increases, so does the noise levels and the randomness of these involuntary motions [
4,
16]. This is relevant to these types of analyses, as they are very sensitive to medication effects. Further, we have shown that these MMS can differentiate those individuals with ASD who never took meds from those do did and identify patterns of medication intake across ages in the ABIDE repository. Likewise, the MMS have differentiated ASD with and without comorbidities (e.g., ADHD) and identified levels of severity according to clinical scales commonly used for research (e.g., ADOS and IQ test scores reported on ABIDE.) Although these may be possible advantages of the methods and MMS data type, we note that one limitation is the novelty of the approach, as very few labs have used these methods. Nevertheless, our lab is making them available through various means to invite reproducibility of results across open access databases such as ABIDE, Kaggle.com, the connectome project and others.
Along those lines, the sensory-motor axes [
62] were only recently added on January of 2019, to the Research Domain Criteria (RDoC) matrix created by the National Institute of Mental Health (NIMH) [
63]. As such, the sensory motor issues discussed here have not been broadly explored in autism, although there is a nascent field researching the neurological issues underlying behaviors across the spectrum of autism (e.g., [
42,
46,
49,
50,
51,
54,
56] among others) is beginning to take advantage of the wearable sensors revolution and measure the biorhythms of the peripheral nervous systems signals [
2,
15]. This is important as this form of feedback, which is disrupted in autism, is critical to continuously guide the central controllers of the brain generating observable behaviors that define autism today. In this sense, the present results warrant further research in autism across the lifespan. This type of work may mark the start of a new medical era for ASD, a prospect that to some may appear as a limitation, while to others, may seem advantageous.