The Metabolic Signature of Autism Spectrum Disorders Using Dried Blood Spots at Birth

Article information

Psychiatry Investig. 2025;22(6):678-686
Publication date (electronic) : 2025 June 16
doi : https://doi.org/10.30773/pi.2024.0293
Department of Child Healthcare, Changzhou Maternal and Child Health Care Hospital, Changzhou, China
Correspondence: Yuping Zhang, MM Department of Child Healthcare, Changzhou Maternal and Child Health Care Hospital, Changzhou 213000, China Tel: +86-138-1355-6560, E-mail: yuping6733123@163.com
Correspondence: Bin Yu, MM Department of Child Healthcare, Changzhou Maternal and Child Health Care Hospital, Changzhou 213000, China Tel: +86-138-6125-1515, E-mail: binyu@njmu.edu.cn
Received 2024 September 23; Revised 2025 March 9; Accepted 2025 April 6.

Abstract

Objective

This study aims to evaluate the metabolic impact of the autism spectrum disorder (ASD) at birth, while such insight may lead to increased understanding of the etiology.

Methods

We performed tandem mass spectrometry (TMS) in a large sample of dried blood spots (DBS) derived at birth from 106 autistic patients and 401 controls, to identify a metabolic signature for ASD. We also examined trait-specific metabolomic patterns within ASDs, focusing on the age, sex, and the comorbidities including the language delay (LD) and global developmental delays.

Results

The results showed that there were no significant differences in metabolism data between ASD patients and controls. The forest plot analysis revealed distinct associations between genetic metabolic substances and autism in male and female populations. Among males, the variable C0 demonstrated a statistically significant association (odds ratio [OR]=1.05, 95% confidence interval [CI]: 1.006–1.096, p=0.024). For females, a significant association was observed with C3 (OR=2.541, 95% CI: 1.089–6.140, p=0.032). Based on their chronological ages of the first diagnosis, autistic individuals were divided in two groups: younger (n=59) or older than 3 years (n=47). The metabolism of succinic acid is increased (p<0.05), as well as carnitines C5:1.

Conclusion

Most analytes included in the TMS screen had no significant differences between the autism group and the control group at birth; however, sex, the age of first diagnose for ASD, and comorbidities may be the important factors affecting metabolic characteristics, as well as the genetic metabolic changes arise after birth.

INTRODUCTION

Autism spectrum disorder (ASD) is a neurodevelopmental disorder (NDD) that is a collection of early social intercourse dysfunction and restrictive and repetitive behavior patterns, interests, activities and sensory processing anomalies [1]. The prevalence of ASD has been steadily increasing in recent decades which is likely associated with changes in diagnostic criteria, improved performance of screening and diagnostic tools, and increased public awareness [2-4]. The recent overall ASD prevalence was 27.6 per 1,000 (one in 36) children aged 8 years and was 3.8 times as prevalent among boys as among girls based on estimates from the Centers for Disease Control and Prevention’s Autism and Developmental Disabilities Monitoring Network [5]. More than 70% of individuals with autism have concurrent developmental conditions, including intellectual disability (ID), language disorders (LD), and attention-deficit/hyperactivity disorder, as well as medical or psychiatric conditions [6,7].

The etiology of ASD is not fully understood, both genetics and environmental factors early in development play a vital role, among it the genetic basis of the disease accounts for about 50% of cases with ASDs, twin studies have provided strong proof of the heritability of ASD, which is estimated to be even as high as 85% [8-10]. Although its exact cause is not known, several factors have been implicated in its etiology, including inborn errors of metabolism (IEMs), defined by a vast array of disorders that are caused by specific enzyme deficiencies or transport protein defects, is as frequent as in 1 in 800 births [11]. IEMs are particularly frequent as diseases of the nervous system, a review reported that clinical symptoms deriving from central nervous system occur in more than 50% of patients with IEMs [12]. The abnormal functioning of molecules involved in IEM can disrupt neurodevelopment at different stages producing a wide repertoire of clinical manifestations that oscillate from severe brain malformations to mild neurological signs such as ASD relating to amino acids (AAs) [13,14] and carnitine [15-17]. Meanwhile, it has been recently proposed that a significant number of ASD cases could be associated with various metabolic abnormalities, some of them identifiable through untargeted metabolomic profiling [18], simultaneously opening additional space for therapeutic attempts [19].

Recent reports emphasized the causal role of IEMs in ASD, underlying autistic symptoms in less than 5% of cases. The literature on the association between ASD and respiratory chain abnormalities is growing, including complex III/IV deficiency and mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes syndrome, as well as glucose-6-phosphate dehydrogenase deficiency [20]. The discovery of this association and the potential opportunity for therapeutic intervention have prompted the renewed search for characteristic neurometabolic biomarkers to aid in early screening and diagnosis of patients. However, previous studies identified metabolic abnormalities (e.g., AA imbalances, carnitine deficiencies) in ASD children at diagnosis [21], which were confounded by postnatal factors like diet or medication. Neonatal metabolic profiling to distinguish inborn vs. acquired metabolic changes was needed. Our core hypothesis proposes that comprehensive metabolic profiling of a long-term stored dried blood spot (DBS) collected at birth was different in ASDs, enabling us to explore early metabolic signatures while minimizing confounding postnatal influences. Such insight may lead to increased understanding of the etiology, as well as provide a theoretical basis for the early detection and early intervention of ASD.

METHODS

Study design

A total of 507 subjects whose DBSs were collected for newborn screening (NBS) at 3 days of age in Changzhou Maternal and Child Health Care Hospital, China, were included in this study. Participants comprised 106 autistic individuals (85 males, 21 females; diagnostic age 2.85±0.79 years) and 401 control participants selected from healthy newborns who received NBS during the same period (197 males, 204 females) (Supplementary Table 1). Among ASD individuals, exclusion criteria comprised the following conditions verified through a three-tiered screening protocol: 1) monogenic disorders including Fragile-X syndrome and tuberous sclerosis were systematically excluded through clinical whole-exome sequencing; 2) comorbid medical conditions including autoimmune diseases were ascertained through electronic health record review. Of 112 initially recruited participants, 6 (5.36%) were excluded through this protocol. We also evaluated the co-occurring conditions of ASDs, and excluded NDDs except GDD and LD. The diagnosis of ASD, GDD, and LD were diagnosed through comprehensive clinical evaluations conducted by a multidisciplinary team consisting of a psychiatrist, three developmental pediatricians, and other professionally trained assessors with ≥5 years of ASD-specific practice. Diagnostic confirmation strictly adhered to Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) criteria [1].

The study design and protocol were reviewed and approved by the ethics committee of the Changzhou Maternity and Child Health Care Hospital affiliated with Nanjing Medical University (Ethics Committee of Changzhou Maternal and Child Health Care Hospital No. 2024[41]). Written informed consent was obtained from all participants’ parent or legal guardian in order to enter clinical and laboratory data from the clinical files into the present study.

NBS by tandem mass spectrometry

DBS were collected from all infants on 903 filter paper (Wallac OY) in 72 hours after birth. In our laboratory, the levels of AAs, carnitines, and ketones in DBS were detected by tandem mass spectrometry (TMS) [21] using the NeoBase Nonderivatized Amino Acid kit (PerkinElmer). Initially, 3.2-mm spots were punched into the V-bottom plate, and 100 μL of the internal standard working solution was added. After incubating for 45 minutes, we transferred the 75 μL solution to the new V-bottomed microporous plate. The AAs, carnitines, and ketones were analyzed using a Waters AcquityTM TQD tandem mass spectrometer (Waters) after 2 hours.

Statistical analysis

Summary statistics for normally distributed quantitative variables are presented as means and standard deviations. For non-normally distributed variables, we used the median and interquartile range. Categorical data were summarized using numbers and percentages. Comparisons between groups were conducted using either a t-test, one-way analysis of variance, or Kruskal-Wallis rank sum test, depending on the data distribution for measurement data, and a χ2 test or Fisher’s exact test for enumeration data. A two-sided p-value of 0.05 was considered statistically significant, subsequently adjusted using the Benjamini-Hochberg method to control the false discovery rate at q<0.05.

The associations between genetic metabolic substances and autism were evaluated using univariate logistic regression analyses, stratified by sex. Odds ratios (ORs) with corresponding 95% confidence intervals (CIs) were calculated to quantify the magnitude and direction of associations. Statistical significance was assessed using two-sided Wald tests, with a threshold of p<0.05. Analyses were performed separately for males and females to explore sex-specific effects. Forest plots were generated to visually summarize the ORs, CIs, and significance levels for each variable, with male and female results displayed in distinct panels.

All statistical analyses were conducted using R software (version 4.4.2, R Foundation for Statistical Computing) with the ggplot2 package for visualization.

RESULTS

Demographic and metabolism characteristics of ASD and controls divided by sex

As it was shown in Supplementary Table 1, ASD patients (80.2% male vs. 19.8% female) had disproportion of male and female compared to controls (49.1% male vs. 50.9% female) (p<0.001). There were no significant differences in other demographic and metabolism data between ASD patients and controls. Based on their sex, participants were further divided in two groups: male (ASD, n=82; control, n=197) or female (ASD, n=21; control, n=204). Demographic data and statistically significant metabolites in ASD participants with respect to controls divided in sex are shown in Table 1. After adjusted by Benjamini-Hochberg, two carnitines, including C16:1 OH and C18:2, were significantly different in ASD patients compared to sex-matched controls. Univariate logistic regression analyses found no differences between ASD and controls except the sex (Table 2). The forest plot analysis revealed distinct associations between genetic metabolic substances and autism in male and female populations. Among males, the variable C0 demonstrated a statistically significant association (OR=1.05, 95% CI: 1.006–1.096, p=0.024), suggesting a marginal increase in odds per unit change. For females, a significant association was observed with C3 (OR=2.541, 95% CI: 1.089–6.140, p=0.032), indicating a pronounced elevation in odds (Figure 1).

Demographic and metabolism characteristics of ASD patients and controls divided by sex (N=507)

Multivariate regression analysis of the associations between genetic metabolic substances and autism

Figure 1.

The vertical line (OR=1) represents the null effect; confidence intervals crossing this line are non-significant. *p<0.05. ASD, autism spectrum disorder; OR, odds ratio; CI, confidence interval.

Metabolism characteristics of ASD divided by co-occurring conditions

Based on their diagnosis, ASD patients were divided in three groups: 46.22% (49/106) of autistic individuals without any other NDDs were classified into group ASD simply, 30.19% (32/106) of that co-occur with GDD, as well as 23.58% (25/106) with LD. The results showed that autistic individuals co-occurring LD were diagnosed earliest, and the diagnostic age for ASD without NDDs was the largest (p=0.005). No significant differences were found in other demographic features and metabolites (p>0.05) (Table 3). We further divided the three co-occurring groups by sex, showing that in the group of ASD with GDD there were four carnitines, including C3, C4, C5, and C6DC were significantly different between male and female (Supplementary Table 2).

Demographic and metabolism characteristics of autistic individuals divided by co-occurring conditions

Metabolic profile for participants divided by age

Based on their chronological ages of the first diagnosis, autistic individuals were divided in two groups: younger (n=59) or older than 3 years (n=47). Demographic data and metabolic profile of all participants in the two age groups are presented in Table 4. After adjusted by Benjamini-Hochberg, the metabolism of succinic acid (SA) is increased (p<0.05), as well as carnitines C5:1. In younger group, the metabolism of C14OH is significantly different in three groups divided by co-occurring conditions (p=0.049) (Supplementary Table 3).

Metabolic profiles for participants divided by the age of first diagnosis for ASD

DISCUSSION

Features of autism may be detected in early childhood, but the diagnosis of autism is usually not made until much later. Current research has not identified clear biomarkers of autism that can provide a basis for diagnostic criteria, and at this stage, metabolomics has been used to characterize specific perturbations of the metabolome in children with autism compared to neurotypical controls, including AAs, acylcarnitines, and various aromatic and phenolic compounds [22,23]. Given that the metabolome of different compounds, subtle changes identified in a clinical condition can also lead to clues to its mechanism or even treatment, but it’s hard to distinguish IEMs as a primary or secondary pathology [24] because of most of existing studies reflecting cross-sectional differences between autism cases and controls. In order to identify early pathophysiological mechanisms of autism, we analyzed the DBS of future ASDs and controls on birth.

Recent reports find that disorders of IEMs caused by mutations in genes coding for enzymes and other important components in genes coding for enzymes and other important components of various metabolic pathways which are known to be associated with ASD, emphasizing the causal role of IEMs in autism [25,26]. However, in our research, few differences in the neonatal metabolomic profiles were found in ASDs at birth, which was opposite to the metabolic abnormities in autistic children at the diagnostic ages in previous researches. It is well-known that metabolic profiles are affected by several factors, such as diet and gut bacteria, which differ between ASD and typically developing populations and can influence each other. Diet is a major source of circulating metabolites, impacting the metabolome directly or indirectly through chemical transformation by the trillions of gut microbes, the microbiome, which has been proposed to modulate complex behaviors [23]. A recent study supported a model whereby ASD-related restricted interests are associated with less-diverse diet, and in turn reduced microbial taxonomic diversity and looser stool consistency [27]. Our results have implications for the interpretation of cause and effect in relation to the analyses of metabolic profiling in ASDs. There is growing interest in the signature of metabolism to ASDs, but our results emphasize the need to consider the (arguably more intuitive) impact of environmental factors after birth.

As ASD is a NDD with a male/female ratio of 4:1, we considered possible interactions of measured metabolites with sex. We split the sample in two age categories (male and female) to understand possible predictive metabolic signatures capable of distinguishing ASD and controls. The present study identified sex-specific associations between metabolic markers and ASD risk, with elevated C0 (free carnitine) levels linked to ASD in males and higher C3 (propionylcarnitine) levels associated with ASD in females. These findings align with emerging evidence highlighting sexually dimorphic metabolic pathways in NDDs. In males, the role of C0 as a risk factor may reflect disruptions in mitochondrial fatty acid oxidation (FAO), a pathway critical for neurodevelopment. Carnitine deficiency has previously been implicated in ASD pathophysiology, as carnitine facilitates the transport of long-chain fatty acids into mitochondria for energy production [28]. Another study identified mutations in carnitine biosynthesis genes (e.g., SLC6A14) in ASD patients, suggesting impaired carnitine metabolism could contribute to neuronal dysfunction [29]. The male-specific association observed here might be amplified by testosterone’s regulatory effects on carnitine palmitoyltransferase activity, potentially exacerbating metabolic vulnerabilities in males [30].

Conversely, the female-specific link between C3 and ASD may relate to propionate metabolism. Elevated C3 levels indicate propionyl-CoA accumulation, which disrupts mitochondrial function and neurotransmitter synthesis. Propionate overload has been shown to induce ASD-like behaviors in animal models [31], and sex differences in propionate detoxification pathways (e.g., estrogen-enhanced urea cycle activity) might explain the female-specific risk [32]. Additionally, epigenetic modifications in ACADS (acyl-CoA dehydrogenase gene), which are influenced by X-chromosome inactivation patterns, could further drive sex-biased propionylcarnitine accumulation [33].

These results underscore the importance of sex-stratified analyses in ASD research. Future studies should explore hormonal interactions with metabolic pathways and validate these biomarkers in longitudinal cohorts [27].

Autism has historically been considered an impairing condition. Based on the current DSM-5 criteria, ASD includes three levels of severity ranging from “requiring support” to “requiring very substantial support.” Judgments of severity are based solely on the characteristics of the two core domains that make up the diagnostic criteria. However, DSM-5 levels are not sufficient to closely map on to an individual’s real-life functioning, including their experiences and challenges, Waizbard-Bartov et al. [34] pointed that common, co-occurring conditions such as ID, and language delays are as impairing to functioning and wellbeing for many ASD individuals as are the core symptoms themselves, when we analyze the metabolism characteristics of ASDs, we split them in three categories. Although no significant differences were found between the three categories, the American Academy of Pediatrics recommended in favor of an initial metabolic workup in patients with GDD, as these conditions are thought to be more likely of a metabolic nature. This dichotomous approach has raised some concerns including difficulties differentiating the early signs of an isolated ASD from GDD and the paucity of medical providers comfortable with recognizing the signs of an underlying IEM [19].

We also considered possible interactions of measured metabolites with participant ages because of the developmental nature for ASD. Profiles of carnitine and acyl-carnitines change significantly during the first year of life, but kept at the same level between 2 and 15 years [35]. We split the sample in two age categories (<3 and ≥3 years) to understand possible predictive metabolic signatures capable of distinguishing ASD and individuals at early stages of ASD development. This threshold is consistent with reliable ASD diagnosis and effectiveness of early intervention. Indeed, the definite diagnosis of ASD is generally made between 3 and 5 years [36]. Moreover, increasing evidences support the effectiveness of early interventions (behavioral, developmental, and educational approaches) in pre-schoolers (aged 24–71 months) with ASD [36]. In our study, the results showed that compared with autistic children younger than 3 years of age, the concentrations of succinylacetone were higher in ASDs older than 3 years of age at diagnosis (p=0.013), and C5:1 also differed between the two ages of diagnosis (p=0.006). The observed age-related differences in SA and carnitine C5:1 metabolism between children diagnosed with ASD before and after age 3 years may reflect distinct developmental trajectories in mitochondrial and lipid metabolic pathways [37,38]. The observation that SA and carnitine C5:1 levels measured at birth differ between children later diagnosed with ASD before and after age 3 years suggests that early metabolic dysregulation may contribute to heterogeneous ASD trajectories. Elevated neonatal SA levels in children diagnosed after age 3 years could reflect congenital mitochondrial dysfunction or disrupted TCA cycle flux, predisposing to progressive neurodevelopmental challenges. Rossignol and Frye [37] hypothesized that mitochondrial defects in ASD may originate prenatally, with SA accumulation acting as a compensatory response to impaired oxidative phosphorylation. This aligns with studies linking neonatal metabolic markers to later ASD risk [38]. Lower neonatal C5:1 levels in children diagnosed <3 years may indicate inborn errors of FAO or carnitine shuttle defects, disrupting energy supply during critical neurodevelopmental windows. However, postnatal environmental factors (e.g., gut microbiota maturation) might modulate these neonatal metabolic risks, explaining later diagnosis in some cases [39]. These findings underscore the potential of neonatal metabolic profiling to stratify ASD subtypes and guide early interventions. Future studies should validate whether SA/C5:1 at birth predict ASD severity or developmental trajectories, while controlling for confounders like maternal metabolic health or perinatal stressors.

Our research, taking into account factors such as sex, comorbidities, and the age of diagnosis, found a few significant differences in the metabolic spectra of those born with ASD. However, our study had several limitations including a small sample size, blood samples collected at the time of diagnosis for new screening data analysis and less stringent inclusion criteria for controls (e.g., healthy controls matched by sex and birth date without NDDs or inherited metabolic diseases). Besides, the metabolic signature of ASDs at diagnostic ages is not considered in this study. Finally, correlation analyses between metabolite levels and clinical scores were not conducted due to the absence of granular clinical data (e.g., Autism Diagnostic Observation Schedule scores) in our cohort. In the future, more samples and more detailed study designs are needed in our study.

In conclusion, relatively little is known about the metabolic impact of the ASD, while such insight may lead to increased understanding of the etiology. We analyzed the metabolic nature at birth, finding that most analytes included in the TMS screen had no significant differences between the autism group and the control group at birth; however, sex, the age of first diagnose for ASD, and comorbidities may be the important factors affecting metabolic characteristics, as well as the genetic metabolic changes arise after birth.

Supplementary Materials

The Supplement is available with this article at https://doi.org/10.30773/pi.2024.0293.

Supplementary Table 1.

Demographic and clinical characteristics of ASD patients and controls

pi-2024-0293-Supplementary-Table-1.pdf
Supplementary Table 2.

Demographic and clinical characteristics of autistic individuals divided by co-occurring conditions and sex

pi-2024-0293-Supplementary-Table-2.pdf
Supplementary Table 3.

Demographic and clinical characteristics of autistic individuals divided by diagnostic age and co-occurring conditions

pi-2024-0293-Supplementary-Table-3.pdf

Notes

Availability of Data and Material

The datasets generated or analyzed during the study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors have no potential conflicts of interest to disclose.

Author Contributions

Conceptualization: Yuping Zhang, Bin Yu. Data curation: Haixin Li, Yuqi Yang, Yuping Zhang. Formal analysis: Yuping Zhang. Funding acquisition: Haixin Li. Investigation: Yuqi Yang. Methodology: Bin Yu, Yuqi Yang. Project administration: Yuping Zhang, Bin Yu. Resources: Bin Yu. Software: Haixin Li. Supervision: all authors. Validation: all authors. Visualization: Haixin Li. Writing—original draft: Haixin Li. Writing—review & editing: all authors.

Funding Statement

The study was supported by the science and technology development fund of Nanjing Medical University-General project (NMUB20210054).

Acknowledgments

We thank all the children and their caregivers for their participation, and we are grateful for the valuable contributions of the colleagues from department of child psychology in Changzhou maternal and child health care hospital, including Yan Sun, Xiaojie Yuan, Jiuling Li, Xueting Wu, and Jinyu Wu, for diagnosing and reporting the autistic children.

References

1. American Psychiatric Association. Diagnostic and statistical manual of mental disorders (DSM-5) (5th ed) Arlington: American Psychiatric Association; 2013.
2. Hirota T, King BH. Autism spectrum disorder: a review. JAMA 2023;329:157–168.
3. Blumberg SJ, Bramlett MD, Kogan MD, Schieve LA, Jones JR, Lu MC. Changes in prevalence of parent-reported autism spectrum disorder in school-aged U.S. children: 2007 to 2011-2012. Natl Health Stat Report 2013;65:1–11.
4. Bachmann CJ, Gerste B, Hoffmann F. Diagnoses of autism spectrum disorders in Germany: time trends in administrative prevalence and diagnostic stability. Autism 2018;22:283–290.
5. Maenner MJ, Warren Z, Williams AR, Amoakohene E, Bakian AV, Bilder DA, et al. Prevalence and characteristics of autism spectrum disorder among children aged 8 years—autism and developmental disabilities monitoring network, 11 sites, United States, 2020. MMWR Surveill Summ 2023;72:1–14.
6. Lai MC, Lombardo MV, Baron-Cohen S. Autism. Lancet 2014;383:896–910.
7. Lai MC, Kassee C, Besney R, Bonato S, Hull L, Mandy W, et al. Prevalence of co-occurring mental health diagnoses in the autism population: a systematic review and meta-analysis. Lancet Psychiatry 2019;6:819–829.
8. Wang L, Wang B, Wu C, Wang J, Sun M. Autism spectrum disorder: neurodevelopmental risk factors, biological mechanism, and precision therapy. Int J Mol Sci 2023;24:1819.
9. Höfer J, Hoffmann F, Kamp-Becker I, Poustka L, Roessner V, Stroth S, et al. Pathways to a diagnosis of autism spectrum disorder in Germany: a survey of parents. Child Adolesc Psychiatry Ment Health 2019;13:16.
10. Havdahl A, Niarchou M, Starnawska A, Uddin M, van der Merwe C, Warrier V. Genetic contributions to autism spectrum disorder. Psychol Med 2021;51:2260–2273.
11. Pampols T. Inherited metabolic rare disease. Adv Exp Med Biol 2010;686:397–431.
12. Saudubray JM, Garcia-Cazorla A. An overview of inborn errors of metabolism affecting the brain: from neurodevelopment to neurodegenerative disorders. Dialogues Clin Neurosci 2018;20:301–325.
13. Tărlungeanu DC, Deliu E, Dotter CP, Kara M, Janiesch PC, Scalise M, et al. Impaired amino acid transport at the blood brain barrier is a cause of autism spectrum disorder. Cell 2016;167:1481–1494.e18.
14. Ghanizadeh A. Increased glutamate and homocysteine and decreased glutamine levels in autism: a review and strategies for future studies of amino acids in autism. Dis Markers 2013;35:281–286.
15. Demarquoy C, Demarquoy J. Autism and carnitine: a possible link. World J Biol Chem 2019;10:7–16.
16. Celestino-Soper PB, Violante S, Crawford EL, Luo R, Lionel AC, Delaby E, et al. A common X-linked inborn error of carnitine biosynthesis may be a risk factor for nondysmorphic autism. Proc Natl Acad Sci U S A 2012;109:7974–7981.
17. Ziats MN, Comeaux MS, Yang Y, Scaglia F, Elsea SH, Sun Q, et al. Improvement of regressive autism symptoms in a child with TMLHE deficiency following carnitine supplementation. Am J Med Genet A 2015;167A:2162–2167.
18. Liu N, Xiao J, Gijavanekar C, Pappan KL, Glinton KE, Shayota BJ, et al. Comparison of untargeted metabolomic profiling vs traditional metabolic screening to identify inborn errors of metabolism. JAMA Netw Open 2021;4e2114155.
19. Glinton KE, Elsea SH. Untargeted metabolomics for autism spectrum disorders: current status and future directions. Front Psychiatry 2019;10:647.
20. Ahmadabadi F, Nemati H, Abdolmohammadzadeh A, Ahadi A. Autistic feature as a presentation of inborn errors of metabolism. Iran J Child Neurol 2020;14:17–28.
21. Barone R, Alaimo S, Messina M, Pulvirenti A, Bastin J, Ferro A, et al. A subset of patients with autism spectrum disorders show a distinctive metabolic profile by dried blood spot analyses. Front Psychiatry 2018;9:636.
22. West PR, Amaral DG, Bais P, Smith AM, Egnash LA, Ross ME, et al. Metabolomics as a tool for discovery of biomarkers of autism spectrum disorder in the blood plasma of children. PLoS One 2014;9e112445.
23. Needham BD, Adame MD, Serena G, Rose DR, Preston GM, Conrad MC, et al. Plasma and fecal metabolite profiles in autism spectrum disorder. Biol Psychiatry 2021;89:451–462.
24. Spilioti M, Evangeliou AE, Tramma D, Theodoridou Z, Metaxas S, Michailidi E, et al. Evidence for treatable inborn errors of metabolism in a cohort of 187 Greek patients with autism spectrum disorder (ASD). Front Hum Neurosci 2013;7:858.
25. Ranganath P, Dalal A. Does every child with autism need investigations for inborn errors of metabolism? Indian Pediatr 2023;60:177–178.
26. Yu TW, Chahrour MH, Coulter ME, Jiralerspong S, Okamura-Ikeda K, Ataman B, et al. Using whole-exome sequencing to identify inherited causes of autism. Neuron 2013;77:259–273.
27. Yap CX, Henders AK, Alvares GA, Wood DLA, Krause L, Tyson GW, et al. Autism-related dietary preferences mediate autism-gut microbiome associations. Cell 2021;184:5916–5931.e17.
28. Voeller KK. Dyslexia. J Child Neurol 2004;19:740–744.
29. Carayol J, Schellenberg GD, Dombroski B, Genin E, Rousseau F, Dawson G. Autism risk assessment in siblings of affected children using sex-specific genetic scores. Mol Autism 2011;2:17.
30. Wu Q, Jiang D, Schaefer NR, Harmacek L, O’Connor BP, Eling TE, et al. Overproduction of growth differentiation factor 15 promotes human rhinovirus infection and virus-induced inflammation in the lung. Am J Physiol Lung Cell Mol Physiol 2018;314:L514–L527.
31. MacFabe DF, Cain DP, Rodriguez-Capote K, Franklin AE, Hoffman JE, Boon F, et al. Neurobiological effects of intraventricular propionic acid in rats: possible role of short chain fatty acids on the pathogenesis and characteristics of autism spectrum disorders. Behav Brain Res 2007;176:149–169.
32. Leinroth AP, Mirando AJ, Rouse D, Kobayahsi Y, Tata PR, Rueckert HE, et al. Identification of distinct non-myogenic skeletal-muscle-resident mesenchymal cell populations. Cell Rep 2022;39:110785.
33. Ebrahimi-Fakhari D, Wahlster L, Bartz F, Werenbeck-Ueding J, Praggastis M, Zhang J, et al. Reduction of TMEM97 increases NPC1 protein levels and restores cholesterol trafficking in Niemann-pick type C1 disease cells. Hum Mol Genet 2016;25:3588–3599.
34. Waizbard-Bartov E, Fein D, Lord C, Amaral DG. Autism severity and its relationship to disability. Autism Res 2023;16:685–696.
35. Cavedon CT, Bourdoux P, Mertens K, Van Thi HV, Herremans N, de Laet C, et al. Age-related variations in acylcarnitine and free carnitine concentrations measured by tandem mass spectrometry. Clin Chem 2005;51:745–752.
36. Narzisi A, Costanza C, Umberto B, Filippo M. Non-pharmacological treatments in autism spectrum disorders: an overview on early interventions for pre-schoolers. Curr Clin Pharmacol 2014;9:17–26.
37. Rossignol DA, Frye RE. A review of research trends in physiological abnormalities in autism spectrum disorders: immune dysregulation, inflammation, oxidative stress, mitochondrial dysfunction and environmental toxicant exposures. Mol Psychiatry 2012;17:389–401.
38. Ruden JB, Dugan LL, Konradi C. Parvalbumin interneuron vulnerability and brain disorders. Neuropsychopharmacology 2021;46:279–287.
39. Sharon G, Cruz NJ, Kang DW, Gandal MJ, Wang B, Kim YM, et al. Human gut microbiota from autism spectrum disorder promote behavioral symptoms in mice. Cell 2019;177:1600–1618.e17.

Article information Continued

Figure 1.

The vertical line (OR=1) represents the null effect; confidence intervals crossing this line are non-significant. *p<0.05. ASD, autism spectrum disorder; OR, odds ratio; CI, confidence interval.

Table 1.

Demographic and metabolism characteristics of ASD patients and controls divided by sex (N=507)

Control
ASD
Statistics p p-adjusted
Total (N=401) Male (N=197) Female (N=204) Total (N=106) Male (N=85) Female (N=21)
Birth weight (g) 3,358±461 3,376±443 3,341±479 3,320±474 3,325±503 3,300±344 F=0.388 0.762
Gestational age (wk) 39.10±1.34 39.04±1.32 39.22±1.35 38.90±1.67 38.75±1.78 39.32±1.03 F=2.466 0.061
Delivery mode χ2=3.043 0.385
 Vaginal delivery 204 (50.9) 103 (52.3) 101 (49.5) 50 (47.2) 37 (43.5) 13 (61.9) 0.250
 Caesarean section 197 (49.1) 94 (47.7) 103 (50.5) 56 (52.8) 48 (56.5) 8 (38.1)
C3 (μmol/L) 1.56 [1.23–1.98] 1.62 [1.31–2.02] 1.45 [1.18–1.94] 1.56 [1.23–1.97] 1.48 [1.18–1.88] 1.96 [1.69–2.24] H=12.32 0.006*§ 0.065
C4DC+C5OH (μmol/L) 0.20 [0.17–0.24] 0.21 [0.17–0.25] 0.19 [0.16–0.23] 0.19 [0.16–0.23] 0.19 [0.16–0.23] 0.21 [0.18–0.27] H=8.72 0.033*§ 0.177
C5 (μmol/L) 0.10 [0.08–0.12] 0.10 [0.08–0.12] 0.10 [0.08–0.13] 0.10 [0.08–0.13] 0.09 [0.08–0.12] 0.11 [0.10–0.13] H=9.62 0.022*§ 0.135
C10:1 (μmol/L) 0.07 [0.06–0.09] 0.08 [0.06–0.09] 0.07 [0.05–0.08] 0.07 [0.05–0.09] 0.07 [0.06–0.09] 0.06 [0.05–0.08] H=9.91 0.019*§ 0.135
C14 (μmol/L) 0.16 [0.13–0.20] 0.17 [0.14–0.22] 0.15 [0.12–0.19] 0.16 [0.13–0.20] 0.16 [0.13–0.20] 0.16 [0.13–0.18] H=11.18 0.011*§ 0.090
C16:1 (μmol/L) 0.16 [0.11–0.20] 0.16 [0.11–0.21] 0.15 [0.10–0.19] 0.15 [0.10–0.21] 0.14 [0.10–0.22] 0.15 [0.10–0.17] H=8.06 0.045*§ 0.215
C16:1OH (μmol/L) 0.03 [0.03–0.04] 0.04 [0.03–0.05] 0.03 [0.03–0.04] 0.03 [0.03–0.04] 0.04 [0.03–0.04] 0.03 [0.03–0.04] H=13.79 0.003*§ 0.043*
C18:2 (μmol/L) 0.27 [0.19–0.35] 0.28 [0.20–0.38] 0.24 [0.18–0.33] 0.27 [0.20–0.35] 0.25 [0.19–0.35] 0.30 [0.25–0.34] H=13.63 0.003*§ 0.040*
*

p<0.05;

mean±standard deviation, one-way analysis of variance;

N (%), Pearson’s chi-squared test;

§

median [interquartile range], Kruskal-Wallis rank sum test; p-adjusted (false discovery rate-adjusted p-values [Benjamini-Hochberg] are reported, with significance set at q<0.05).

ASD, autism spectrum disorder

Table 2.

Multivariate regression analysis of the associations between genetic metabolic substances and autism

Variables Control (N=401) ASD (N=106) OR 95% CI p
Sex
 Male 197 (49.1) 85 (80.2) 1 (ref)
 Female 204 (50.9) 21 (19.8) 0.232 0.132–0.392 0.001*
Birth weight (g) 3,358±461 3,320±474 0.999 0.999–1.001 0.727
Gestational age (wk) 39.10±1.34 38.90±1.67 0.912 0.743–1.116 0.373
Delivery mode
 Vaginal delivery 204 (50.9) 50 (47.2) 1 (ref)
 Caesarean section 197 (49.1) 56 (52.8) 1.039 0.655–1.648 0.871
C0 (μmol/L) 21.8 [16.7–26.3] 22.0 [18.3–26.4] 1.029 0.992–1.068 0.123
C3 (μmol/L) 1.56 [1.23–1.98] 1.56 [1.23–1.97] 1.077 0.699–1.640 0.731

Analyses were adjusted for C4DC+C5OH, C5, C10:1, C14, C16:1, C16:1OH, and C18:2. Values are presented as mean±standard deviation, number (%), or median [interquartile range].

*

p<0.05.

ASD, autism spectrum disorder; OR, odds ratio; CI, confidence interval

Table 3.

Demographic and metabolism characteristics of autistic individuals divided by co-occurring conditions

Total (N=106) ASD simply (N=49) ASD with GDD (N=32) ASD with LD (N=25) Statistics p
Sex 0.633
 Male 85 (80.2) 41 (83.7) 24 (75.0) 20 (80.0)
 Female 21 (19.8) 8 (16.3) 8 (25.0) 5 (20.0)
Diagnostic age (yr) 2.85±0.79 3.08±0.83 2.79±0.74 2.46±0.62 F=5.666 0.005*
Birth weight (g) 3,320±474 3,320±509 3,323±432 3,316±475 F=0.002 0.998
Gestational age (wk) 38.90±1.67 39.00±1.69 39.10±1.04 38.40±2.19 F=1.32 0.272
Delivery mode 0.494
 Vaginal delivery 44 (47.2) 26 (53.1) 14 (43.8) 10 (40.0)
 Caesarean section 56 (52.8) 23 (46.9) 18 (56.2) 15 (60.0)
ALA (μmol/L) 308 [254–357] 312 [258–367] 280 [251–349] 310 [244–333] H=1.720 0.423§
ARG (μmol/L) 9.93 [6.18–17.5] 9.67 [5.41–16.2] 10.5 [7.62–18.6] 9.95 [3.62–13.8] H=2.845 0.241§
CIT (μmol/L) 12.5 [10.1–15.6] 13.4 [10.3–15.8] 12.2 [10.0–16.6] 11.5 [10.1–14.6] H=0.976 0.614§
GLY (μmol/L) 447 [372–517] 446 [373–494] 440 [375–486] 481 [350–563] H=0.495 0.781§
LEU.ILE.PRO.OH (μmol/L) 148 [130–184] 148 [131–186] 143 [127–170] 151 [140–183] H=1.264 0.532§
MET (μmol/L) 22.8 [17.4–26.5] 24.0 [17.3–26.8] 22.1 [17.9–26.5] 22.0 [18.1–26.0] H=0.208 0.901§
ORN (μmol/L) 118 [103–146] 119 [103–146] 115 [105–149] 118 [91.5–140] H=0.433 0.805§
PHE (μmol/L) 54.4 [46.3–59.6] 53.0 [45.7–57.8] 52.4 [46.4–56.7] 58.3 [48.9–63.5] H=4.350 0.114§
PRO (μmol/L) 177 [155–207] 188 [153–211] 176 [158–196] 174 [157–199] H=0.745 0.689§
SA (μmol/L) 0.57 [0.49–0.66] 0.58 [0.49–0.69] 0.54 [0.48–0.62] 0.56 [0.51–0.63] H=0.339 0.844§
TYR (μmol/L) 106±31.6 109±31.4 105±29.5 104±35.1 F=0.327 0.722
VAL (μmol/L) 134 [115–155] 134 [113–160] 136 [110–148] 130 [123–147] H=0.045 0.978§
C0 (μmol/L) 22.0 [18.3–26.4] 22.2 [19.0–28.3] 21.1 [17.9–24.6] 23.1 [16.4–28.6] H=2.243 0.326§
C2 (μmol/L) 18.1 [14.3–22.8] 18.3 [14.1–22.5] 18.5 [15.0–24.2] 17.0 [14.0–22.6] H=0.667 0.716§
C3 (μmol/L) 1.56 [1.23–1.97] 1.52 [1.18–1.93] 1.71 [1.39–1.96] 1.58 [1.14–2.40] H=1.086 0.581§
C3DC+C4OH (μmol/L) 0.09 [0.07–0.15] 0.09 [0.06–0.12] 0.11 [0.07–0.16] 0.09 [0.08–0.12] H=0.735 0.693§
C4 (μmol/L) 0.19 [0.16–0.24] 0.19 [0.16–0.22] 0.21 [0.16–0.26] 0.20 [0.17–0.25] H=1.087 0.581§
C4DC+C5OH (μmol/L) 0.19 [0.16–0.23] 0.19 [0.17–0.23] 0.19 [0.16–0.23] 0.19 [0.18–0.23] H=0.239 0.887§
C5 (μmol/L) 0.10 [0.08–0.13] 0.10 [0.08–0.12] 0.11 [0.08–0.12] 0.09 [0.08–0.13] H=0.158 0.924§
C5:1 (μmol/L) 0.01 [0.01–0.01] 0.01 [0.01–0.01] 0.01 [0.01–0.01] 0.01 [0.01–0.01] H=0.001 0.999§
C5DC+C6OH (μmol/L) 0.08 [0.07–0.10] 0.08 [0.08–0.10] 0.09 [0.07–0.10] 0.08 [0.07–0.10] H=0.101 0.951§
C6 (μmol/L) 0.04 [0.03–0.05] 0.04 [0.03–0.05] 0.04 [0.03–0.05] 0.04 [0.03–0.04] H=0.245 0.885§
C6DC (μmol/L) 0.08 [0.06–0.10] 0.08 [0.07–0.10] 0.08 [0.07–0.09] 0.08 [0.06–0.11] H=0.561 0.755§
C8 (μmol/L) 0.06 [0.04–0.08] 0.05 [0.04–0.07] 0.06 [0.04–0.07] 0.06 [0.04–0.08] H=0.156 0.925§
C8:1 (μmol/L) 0.12 [0.10–0.15] 0.12 [0.10–0.15] 0.12 [0.09–0.16] 0.11 [0.09–0.14] H=1.486 0.476§
C10 (μmol/L) 0.07 [0.05–0.10] 0.07 [0.06–0.08] 0.08 [0.05–0.09] 0.07 [0.05–0.10] H=0.174 0.917§
C10:1 (μmol/L) 0.07 [0.05–0.09] 0.07 [0.06–0.09] 0.07 [0.06–0.09] 0.06 [0.05–0.08] H=5.869 0.053§
C10:2 (μmol/L) 0.07 [0.05–0.09] 0.01 [0.01–0.01] 0.01 [0.01–0.01] 0.01 [0.01–0.02] H=0.071 0.965§
C12 (μmol/L) 0.05 [0.03–0.07] 0.07 [0.05–0.10] 0.06 [0.05–0.09] 0.06 [0.05–0.09] H=0.344 0.842§
C12:1 (μmol/L) 0.05 [0.03–0.07] 0.04 [0.03–0.07] 0.05 [0.03–0.07] 0.05 [0.03–0.06] H=0.614 0.736§
C14 (μmol/L) 0.17±0.06 0.17±0.06 0.16±0.06 0.17±0.07 F=0.179 0.673
C14:1 (μmol/L) 0.07 [0.05–0.10] 0.07 [0.05–0.10] 0.08 [0.05–0.11] 0.06 [0.05–0.08] H=0.913 0.634§
C14:2 0.02 [0.01–0.02] 0.02 [0.01–0.02] 0.02 [0.02–0.02] 0.02 [0.01–0.02] H=5.072 0.079§
C14OH 0.01 [0.01–0.01] 0.01 [0.01–0.01] 0.01 [0.01–0.01] 0.01 [0.01–0.01] H=0.110 0.947§
C16 (μmol/L) 2.88±1.20 2.86±1.20 2.96±1.18 2.84±1.28 F=0.196 0.659
C16:1 (μmol/L) 0.15 [0.10–0.21] 0.14 [0.10–0.20] 0.16 [0.10–0.23] 0.14 [0.11–0.19] H=1.041 0.594§
C16:1OH (μmol/L) 0.03 [0.03–0.04] 0.04 [0.03–0.04] 0.04 [0.03–0.04] 0.03 [0.03–0.04] H=0.063 0.969§
C16OH (μmol/L) 0.02 [0.01–0.02] 0.02 [0.01–0.02] 0.02 [0.01–0.02] 0.01 [0.01–0.02] H=0.862 0.650§
C18 (μmol/L) 0.74 [0.61–1.00] 0.71 [0.61–1.02] 0.78 [0.65–0.95] 0.76 [0.62–0.91] H=0.043 0.979§
C18:1 (μmol/L) 1.52 [1.25–1.71] 1.51 [1.31–1.71] 1.55 [1.21–1.61] 1.45 [1.25–1.73] H=0.192 0.908§
C18:1OH (μmol/L) 0.02 [0.02–0.03] 0.02 [0.02–0.03] 0.02 [0.02–0.02] 0.02 [0.02–0.03] H=1.102 0.576§
C18:2 (μmol/L) 0.27 [0.20–0.35] 0.28 [0.23–0.36] 0.26 [0.21–0.31] 0.29 [0.18–0.34] H=1.182 0.554§
C18OH (μmol/L) 0.01 [0.01–0.01] 0.01 [0.01–0.01] 0.01 [0.01–0.01] 0.01 [0.01–0.01] H=0.292 0.864§
*

p<0.05;

N (%), Fisher’s exact test;

mean±standard deviation, one-way analysis of variance;

§

median [interquartile range], Kruskal-Wallis rank sum test.

ASD, autism spectrum disorder; GDD, global developmental delay; LD, language delay

Table 4.

Metabolic profiles for participants divided by the age of first diagnosis for ASD

<3 years (N=59) ≥3 years (N=47) Statistics p p-adjusted
Sex 0.048* 0.096
 Male 43 (72.9) 42 (89.4)
 Female 16 (27.1) 5 (10.6)
Birth weight (g) 3,341±463 3,294±491 t=0.499 0.618 0.694
Gestational age (wk) 38.90±1.73 38.80±1.61 t=0.394 0.694 0.694
Delivery mode 0.438 0.657
 Vaginal delivery 30 (50.8) 20 (42.6)
 Caesarean section 29 (49.2) 27 (57.4)
Succinic acid (μmol/L) 0.54 [0.48–0.61] 0.59 [0.52–0.78] w=997 0.013*§ 0.039*
C5:1 0.006* 0.036*
 0 4 (6.8) 12 (25.5)
 0.01 52 (88.1) 35 (74.5)
 0.02 3 (5.1) 0 (0.0)
*

p<0.05;

N (%), Fisher’s exact test;

mean±standard deviation, t-test;

§

median [interquartile range], Wilcoxon rank sum test; p-adjusted (false discovery rate-adjusted p-values [Benjamini-Hochberg] are reported, with significance set at q<0.05).

ASD, autism spectrum disorder