Bioenergetic variation is related to autism symptomatology

Autism spectrum disorder (ASD) has been associated with mitochondrial dysfunction but few studies have examined the relationship between mitochondrial function and ASD symptoms. We measured Complex I and IV and citrate synthase activities in 76 children with ASD who were not receiving vitamin supplementation or medication. We also measured language using the Preschool Language Scales or Clinical Evaluation of Language Fundamentals, adaptive behavior using the Vineland Adaptive Behavioral Scale, social function using the Social Responsiveness Scale and behavior using Aberrant Behavior Checklist, Childhood Behavior Checklist and the Ohio Autism Clinical Impression Scale. Children with ASD demonstrated significantly greater variation in mitochondrial activity compared to controls with more than expected ASD children having enzyme activity outside of the normal range for Citrate Synthase (24%), Complex I (39%) and Complex IV (11%). Poorer adaptive skills were associated with Complex IV activity lower or higher than average and lower Complex I activity. Poorer social function and behavior was associated with relatively higher Citrate Synthase activity. Similar to previous studies we find both mitochondrial underactivity and overactivity in ASD. This study confirms an expanded variation in mitochondrial activity in ASD and demonstrates, for the first time, that such variations are related to ASD symptoms.


Background
Autism spectrum disorder (ASD) is a behaviorally defined disorder which has continued to rise in prevalence over the last two decades and now affects~2% of children in the United States (Zablotsky et al. 2015). Although its etiology remains poorly understood, recent studies suggest a link to mitochondrial dysfunction (Frye and Rossignol 2011;Rossignol and Frye 2012;Legido et al. 2013). The exact nature of mitochondrial dysfunction and its link to symptoms, however, remains an area of active investigation. A recent meta-analysis demonstrated that classic mitochondrial disease is found in approximately 5% (Rossignol and Frye 2012) of children with ASD, yet 30-50% have biomarkers of mitochondrial dysfunction (Frye 2012a, b;Rossignol and Frye 2012) and a higher rate of abnormal electron transport chain (ETC) activity has been documented in lymphocytes and granulocytes (Giulivi et al. 2010;Napoli et al. 2014) and post-mortem brain tissue (Palmieri et al. 2010). This suggests that ASD may have a unique relationship to abnormal mitochondrial function.
Major criteria that define classic mitochondrial disease include unequivocal genetic mutations, severe depressions (i.e., <30%) in ETC function or syndromic presentation (Frye and Rossignol 2011). However, our meta-analysis demonstrated that the minority (~25%) of children with ASD and mitochondrial disease have identifiable genetic mutations that explain the mitochondrial disease (Rossignol and Frye 2012) and many case reports and series have described moderate, rather than severe, deficiencies in ETC activity (Frye 2012a(Frye , b, 2013b. Perhaps more unique is the fact that ETC activity in muscle (Frye and Naviaux 2011;Frye 2012a, b), skin (Frye et al. 2013b), buccal epithelium (Legido et al. 2013;Goldenthal et al. 2015;Frye et al. 2016a) and brain (Palmieri et al. 2010) has been documented to be significantly increased, rather than decreased, in individuals with ASD. This is consistent with our in vitro data showing elevated mitochondrial respiration in lymphoblastoid cell lines (LCLs) derived from children with ASD (Rose et al. 2014a, b;Frye et al. 2016c;Rose et al. 2017).
In order to investigate variation in ETC function in children with ASD without using invasive procedures like a muscle biopsy, Goldenthal et al. (Goldenthal et al. 2015) used a validated buccal swab method  to collect data on the largest sample to date. Goldenthal et al. (Goldenthal et al. 2015) verified the expanded variation in mitochondrial function in children with ASD and demonstrated that children with ASD demonstrated both depressions and elevations in ETC function. Still the relationship between ETC complex activity and cognitive and behavioral symptoms has not been investigated in a large sample. Thus this study uses the buccal swab technique to examine ETC function in a large sample of ASD children and determine the correspondence between variations in mitochondrial function and cognitive development and behavior. Most importantly, participants in this study were not taking any supplements or medications that could interfere with mitochondrial function at the time of testing.

Methods
The study was approved by the Institutional Review Board at the University of Arkansas for Medical Sciences (Little Rock, AR). Parents of participants provided written informed consent. All experiments were performed in accordance with relevant guidelines and regulations.

ASD participants
76 individuals with ASD who met inclusion and exclusion criteria had ETC complex function measured with 61 undergoing cognitive and behavioral evaluations. Clinical characteristics of the participants are provided in Table 1.
Comorbid conditions were both derived from a parent reported medical questionnaire and from review of conditions diagnosed in the medical records. Regression (defined as loss of already obtained skills) was defined in detail in our questionnaire. Questions regarding regression included the timing, specific skills lost, duration of the regression, trigger and whether there were multiple regressions. Since this paper is not specifically on regression, we have not summarized these details in the table. This method for assessing medical comorbidies has been used in several of our previous studies (Frye et al. 2016b, d;Delhey et al. 2017;Frye et al. 2017b).
Since there are numerous medical conditions, we summarized them into categories. Gastrointestinal disorders included constipation, diarrhea, stomach/abdominal pain, gastroesophageal reflux disease, vomiting, and feeding problems. Allergic disorders included allergies, asthma, breathing problems, other lung problems, and allergic skin conditions. Psychiatric disorders included depression, bipolar, anxiety disorder, obsessive compulsive disorder, attention deficit with or without hyperactivity. Immune disorder included ear infections, sinusitis, throat infections and other immune problems. Neurological disorders included hearing problems, headaches, vision problems, microcephaly, macrocephaly, seizures, opthalmoplegia, strabismus, cerebral palsy, tics, and muscle weakness. General health disorders included dental problems, fatigue, lack of sweating, exercise intolerance, heat intolerance, bone and joint problems, and blood problems and anemia. Cardiovascular disorders included heart conditions, syncope, congenital heart disease, and heart failure. Growth disorders included growth hormone deficiency, failure to thrive, short stature, tall stature, overweight and obesity. Endocrine disorders includes thyroid or other endocrine problems. Genitourinary disorders includes kidney, bladder, and genital problems. Inclusion criteria were (i) age 3 to 14 years of age and (ii) ASD diagnosis. Exclusion criteria were (i) chronic treatment with medications that would detrimentally affect mitochondrial function such as antipsychotic medications; (ii) vitamin or mineral supplementation exceeding the recommended daily allowance, and (iii) prematurity.
The ASD diagnosis was defined by one of the following: (i) a gold-standard diagnostic instrument such as the Autism Diagnostic Observation Schedule and/or Autism Diagnostic Interview-Revised; (ii) the state of Arkansas diagnostic standard, defined as agreement of a physician, psychologist and speech therapist; and/or (iii) Diagnostic Statistical Manual (DSM) diagnosis by a physician along with standardized validated questionnaires and diagnosis confirmation by the Principal Investigator. We have validated that this criteria captures an accurate diagnosis of ASD in our previous studies by re-evaluating a portion of the participants with the Autism Diagnostic Interview-Revised and determining that they all were well within the diagnostic criteria for ASD (Frye et al. 2016b, d;Delhey et al. 2017;Frye et al. 2017b).

Historical healthy controls
Historical controls of similar age and gender included 68 healthy individuals without neurological disease as described in previous studies (Goldenthal et al. 2015). Controls ranged in age from 3 to 21 years of age [mean (SD) 10.1 years (4.6 years)] with 33 (49%) being female. These controls were evaluated with the same methodology used in this study so that the data from the controls were truly comparable to the data collected in the current study. In a previous report, it was found that there was no correlation between enzyme activities and age and no difference in protein activities across ethnicity or race in both controls and mitochondrial disease patients .

Contemporaneous healthy controls
Contemporaneous controls of similar and age and gender included 14 generally healthy children without significant chronic health conditions. These controls ranged in age from 1 to 15 years of age [Mean (SD) 9.1 years (4.8 years)] with 6 (42%) being female. These individuals underwent the Vineland Adaptive Behavior Scale 2nd Edition (VABS) evaluation as well as an evaluation for mitochondrial function so their data could be included in the analysis examining the relationship between the VABS and variation in mitochondrial function in order to investigate whether the patterns found for children with ASD were in line with those found in typically developing individuals.

Cognitive and behavioral outcome measures
As mentioned in our previous study (Frye et al. 2016d), our research staff was trained by a multispecialty team consisting of two licensed psychologists and a speech therapist prior to performing assessments. During the study a research psychologist supervised research staff and provided feedback and retraining if necessary.
Observer rated measures included language and the Ohio Autism Clinical Impression Scale (OACIS). Parents completed the Aberrant Behavior Checklist (ABC), Social Responsiveness Scale (SRS) the VABS Survey Interview Form and Childhood Behavioral Checklist (CBCL). Language was assessed by an ability-appropriate instrument, the CELF-preschool-2, the CELF-4 or the Preschool Language Scale-5 (PLS-5). The CELF is a standardized, well-validated instrument that has been used in studies on ASD (Verly et al. 2014;Edgar et al. 2015). The PLS-5 is a standardized, well-validated instrument that measures subtle changes in verbal communication, particularly in pre-verbal children (Volden et al. 2011). The standardized summary score of each instrument (mean 100, standard deviation 15) was used as the outcome measure.
The OACIS is an observer-rated scale sensitive to clinically meaningful changes in ASD symptoms (Butter and Mulick 2006) which has been shown to have good inter-rater and cross-cultural reliability (Choque Olsson and Bolte 2014) and has been used in several ASD clinical trials . The VABS is a reliable and valid measure of the ability to perform age-appropriate everyday skills, including communication, daily living, social and motor skills, through a 20-30 min structured interview with a caretaker (Frye et al. 2013a). The ABC is a 58-item questionnaire (Frye et al. 2013a) that measures disruptive behaviors and has convergent and divergent validity (Kaat et al. 2014). The SRS is a 65-item questionnaire that measures the severity of social skill deficits across five domains (Constantino 2002) which has been shown to have good correspondence to the gold-standard instrument (Murray et al. 2011). The CBCL measures general behavioral problems in children and has been described as robust for use in measuring behavior in children with ASD (Hanratty et al. 2015).

Measures of mitochondrial function
The buccal cells were collected using Catch-All Buccal Collection Swabs (Epicentre Biotechnologies, Madison, WI). Four swabs were collected by firmly pressing a swab against the inner cheek while twirling for 30s. Swabs were clipped and placed in 1.5 ml microcentrifuge tubes that were labeled and placed on dry ice for overnight transportation to the Goldenthal laboratory.
Buccal extracts were prepared using an ice-cold buffered solution (Buffer A, ABCAM) containing protease inhibitor cocktail and membrane solubilizing non-ionic detergent and cleared of insoluble cellular material by high speed centrifugation at 4°C. Duplicate aliquots of the protein extract were analyzed for protein concentration using the bicinchoninic acid method (Pierce Biotechnology, Rockford, IL). Samples were typically stored at −80°C for up to 1 week prior to enzymatic analysis.
Dipstick immunocapture assays measured ETC Complex I activity using 50 μg of extracted protein (Ezugha et al. 2010;Goldenthal et al. 2012;Yorns et al. 2012). Signals were quantified using a Hamamatsu immunochromato reader (MS 1000 Dipstick reader). Raw milliAbsorbance results were corrected for protein concentration and data were expressed as percentages of the values obtained with control extracts run on the same assay. ETC Complex IV and CS activity were assessed using standard spectrophotometric procedures in 0.5 ml reaction volume. Specific activities of respiratory complexes were initially expressed as nanomoles/min/mg protein.
For comparison to controls, Complex I and IV activities were normalized to CS activity.

Statistical analysis
Analyses were performed using SAS 9.4 (SAS Institute Inc., Cary, NC). Graphs were produced using Excel version 14.0 (Microsoft Corp, Redmond, WA). This analysis is modeled on our recently analysis of bioenergetic variation in a genetic syndrome (Frye et al. 2016a).
Normal control values were based upon the established historical controls from the Goldenthal laboratory published previously. The two-tailed 95% normative range was considered within 1.96 +/− standard deviation (SD) of the mean, meaning that abnormal values were above 97.5% or below 2.5% of the normal distribution as defined by the historical control mean and standard deviation for each ETC complex activity measured. For comparison to normal historical controls, ETC complex activity was normalized to CS activity.
Abnormalities in ETC complex activity were analyzed in several ways. First, the number of values outside of the normal range was calculated for each ETC complex measured. The binomial probability was calculated to determine if the number of values outside the normal range were significantly different than expected by chance. Second, the ASD group mean complex activity values were compared to the entire control population using t-tests. Third, a difference in the group variance was compared using F-tests.
To better understand the relationship between cognitive and behavioral indexes and mitochondrial enzyme activity we used a general linear model with the cognitive and behavioral indexes as the dependent measure and mitochondrial enzyme activity as predictor variables. ETC complex activity was not normalized in this analysis so the effect of CS could be differentiated from the effect of individual complex activity. Since there is data suggesting that mitochondrial enzyme activity both abnormal and below normal is related to ASD, a second order polynomial is used to model the predictor variables in order to account for a potential curvilinear relationship. The generalized linear model module 'genmod' in SAS was used with an alpha of 0.05 with a backward stepwise elimination of non-significant variables. Since adaptive behavior is known to affect social (i.e., SRS (Hus et al. 2013)) and behavior (e.g. ABC (Rattaz et al. 2015)) symptoms in ASD, the composite score of the VABS was used as a covariate for the SRS, ABC, CBCL and OASIS outcome measures.

Results
Mitochondrial function in autism spectrum disorder Table 2 outlines the results of the statistical analyses. Overall 62% of ASD individuals showed mitochondrial enzyme activity outside the control range (Fig. 1). Complex I activity was outside of the control range for 39% of ASD participants, which is significantly more than expected by chance alone. Variation in Complex I activity was significantly greater in ASD participants as compared to the control group and mean Complex I activity for ASD participants was significantly lower than the control group. Complex IV activity was outside the normal range for 11% of the participants, which is significantly more than expected by chance alone. Mean Complex IVactivity was significantly lower in ASD participants as compared to the control group, and variation was significantly greater in the ASD group as compared to the control group. The ratio of Complex I to Complex IV was abnormal in 24% of the ASD participants. The average ratio was not significantly different than controls but the variation was greater in the ASD group as compared to the control group. Table 3 outlines the overlap of complex function abnormalities within the 47 ASD participants with atypical mitochondrial function. Overall, most, 57% (27/47), of the ASD participants with mitochondrial enzyme abnormalities had only one enzyme affected, while 32% (15/47) had two enzymes affects and 7% (3/47) had all three enzymes affected. Complex I and Citrate Synthase (CS) were isolated a majority of the time while Complex IV abnormalities were commonly associated with either a Complex I or CS abnormality. Table 4 provides the details of the overlap of the underactivity and overactivity of specific mitochondrial enzymes. Although some interesting patterns emerge, the number of samples limit conclusions that can be made from these patterns.
Mitochondrial enzymes demonstrated moderate correlations with one another (Fig. 2). CS activity was positively correlated with activity of Complex I (r = 0.47, p = 0.0001) and Complex IV (r = 0.65, p < 0.0001) and Complex I activity was positively correlated with activity of Complex IV (r = 0.50, p < 0.0001).

, better Communication [Linear
Mitochondrial enzyme activity was not related to the other behavioral or language assessments.

Discussion
Previous studies have suggested that abnormal mitochondrial function is associated with ASD but the exact nature of this mitochondrial dysfunction and how it is associated with ASD symptoms is not clear. In this study we have measured mitochondrial enzyme activity on the one of the largest cohorts of individuals with ASD studied to date with concurrent measurement of symptoms in a subset. This study demonstrates that there is an expanded variation in mitochondrial enzyme activity in children with ASD, with enzyme activity that is both higher and lower than found in typically developing individuals. These variations are associated with ASD symptoms as well as measures of general development with some of these relationships non-linear such that both relatively higher and lower enzyme activity values are similarly associated with unfavorable alternations in behavior and development. A strength of this study is that the children in this study were not taking any medications or supplements that could potentially influence mitochondrial activity. Indeed, many medications commonly used to treat children with ASD, such as antipsychotics, can inhibit mitochondrial function (Wallace 2014) and many supplements commonly used in ASD can alter mitochondrial enzyme activity (Delhey et al. 2017).
In classic mitochondrial disease, the mitochondrion is thought to be severely dysfunctional with significant depression in mitochondrial respiration. Criteria such as the modified Walker's criterion reflects this notion and are commonly used to diagnose mitochondrial disease, although other criteria like the Morava criteria, which are more clinically based, are also sometimes used (Frye and Rossignol 2011). Classic mitochondrial disease is often found to have a genetic etiology, Fig. 1 Activity of Citrate Synthase and activities of Complex I and IV normalized to Citrate Synthase for children with Autism Spectrum Disorder (ASD). For each graph the control mean and minimal and maximum limits of normal controls are provided on the right.
Individual participant data is provided in the center of each graph with individuals having values outside of the normal range highlighted in red. The mean of the ASD group is depicted on the left of the individual data points as diamond although this seems to be more relevant to adults with classic forms of mitochondrial disease than children. It is also being recognized that many genetic and non-genetic disorders can cause mitochondrial dysfunction that can may require treatments targeting the mitochondria (Niyazov et al. 2016). Additionally, the great majority of children with ASD who are diagnosed with mitochondrial disease do not have a simple genetic alteration to explain their problem in mitochondrial function (Rossignol and Frye 2012). Many of the children with ASD in this study were found to have a decrease in mitochondrial activity, particularly below normal activity in Complex I. This is consistent with muscle biopsies of children with ASD diagnosed with mitochondrial disease, where the most common ETC deficit was Complex I (Rossignol and Frye 2012). Possibly unique to ASD is mitochondrial enzyme overactivity. Increases in ETC complex activity, rather than deficiencies in ETC complex activity, has been reported in several tissues derived from individuals with ASD, including muscle (Frye and Naviaux 2011;Frye 2012a, b), skin (Frye et al. 2013b), buccal (Legido et al. 2013;Goldenthal et al. 2015;Frye et al. 2016a) and brain (Palmieri et al. 2010) tissue. This increase in respiratory function is consistent with our cell line model of mitochondrial function in ASD. Indeed, we have repeatedly demonstrated in the laboratory that a subset of lymphoblastoid cell lines derived from children with ASD have respiratory rates approximately twice that of control cell lines (Rose et al. 2014a(Rose et al. , b, 2015aFrye et al. 2016cFrye et al. , 2017aRose et al. 2017). We recently demonstrated that this alteration in respiration is associated with worse repetitive behaviors . We believe that this increase in respiratory rate may be an adaptation designed to resistant toxicants associated with ASD, perhaps because of previous exposure to environmental toxicants (Rose et al. 2014a(Rose et al. , b, 2017Frye et al. 2017a). We have also demonstrated that this increased respiratory rate results in an increased susceptibility of the mitochondrial to in vitro increases in ROS (Rose et al. 2014a, b). This study confirms the association with increased ETC complex activity and ASD. Further research will be needed to better understand this unique alteration in mitochondrial physiology.
The above speculations must be tempered with the limitations of this study. We have only demonstrated an association between behavior and biomarkers of mitochondrial activity. This association does not imply causation. Indeed, there may be many other factors that influence both mitochondrial activity and behavior. In addition, both behavior and mitochondrial function demonstrate significant variability. This indeed suggests that further research will be needed to determine more details of this relationship. For example, there may be several subgroups of children with specific medical conditions and ASD which have abnormal mitochondrial function such as those with epilepsy or Downs syndrome. Despite these limitations, this report helps defined a connection between behavior and biology which is very intriguing. Further studies, using larger samples sizes will be needed.
Author's contribution Leanna Delhey, John Slattery, Stephen Kahler and Richard E. Frye designed the study and wrote the protocol, collected the data, analyzed the data, drafted the manuscript. Ekim Nur Kilinc, Li Yin and Marie Tippett analyzed and collected the data. Rebecca Wynne, Shannon Rose, Shirish Damle and Agustin Legido performed laboratory analysis. Michael J Goldenthal designed the study and supervised the laboratory analysis of mitochondrial function. All authors contributed to and have approved the final manuscript.

Compliance with ethical standards
Conflict of interest All authors' have no conflicts of interest to declare.
Open Access This article is distributed under the terms of the Creative Comm ons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.