Differential Metabolites in Chinese Autistic Children: A Multi-Center Study Based on Urinary 1H-NMR Metabolomics Analysis

Background: Autism spectrum disorder (ASD) is a group of early-onset neurodevelopmental disorders. However, there is no valuable biomarker for the early diagnosis of ASD. Our large-scale and multi-center study aims to identify metabolic variations between ASD and healthy children and to investigate differential metabolites and associated pathogenic mechanisms.

Methods: One hundred and seventeen autistic children and 119 healthy children were recruited from research centers of 7 cities. Urine samples were assayed by ¹H-NMR metabolomics analysis to detect metabolic variations. Multivariate statistical analysis, including principal component analysis (PCA), and orthogonal projection to latent structure discriminant analysis (OPLS-DA), as well as univariate analysis were used to assess differential metabolites between the ASD and control groups. The differential metabolites were further analyzed by receiver operating characteristics (ROC) curve analysis and metabolic pathways analysis.

Results: Compared with the control group, the ASD group showed higher levels of glycine, guanidinoacetic acid, creatine, hydroxyphenylacetylglycine, phenylacetylglycine, and formate and lower levels of 3-aminoisobutanoic acid, alanine, taurine, creatinine, hypoxanthine, and N-methylnicotinamide. ROC curve showed relatively significant diagnostic values for hypoxanthine [area under the curve (AUC) = 0.657, 95% CI 0.588 to 0.726], creatinine (AUC = 0.639, 95% CI 0.569 to 0.709), creatine (AUC = 0.623, 95% CI 0.552 to 0.694), N-methylnicotinamide (AUC = 0.595, 95% CI 0.523 to 0.668), and guanidinoacetic acid (AUC = 0.574, 95% CI 0.501 to 0.647) in the ASD group. Combining the metabolites creatine, creatinine and hypoxanthine, the AUC of the ROC curve reached 0.720 (95% CI 0.659 to 0.777). Significantly altered metabolite pathways associated with differential metabolites were glycine, serine and threonine metabolism, arginine and proline metabolism, and taurine and hypotaurine metabolism.

Conclusions: Urinary amino acid metabolites were significantly altered in children with ASD. Amino acid metabolic pathways might play important roles in the pathogenic mechanisms of ASD.

Introduction

Autism spectrum disorder (ASD) is a group of early-onset neurodevelopmental disorders characterized by social communication difficulties, narrow interests, and repetitive stereotyped behaviors (1). In the United States, the prevalence of ASD has increased in recent years, ranging from 1.57% in 2009 to 1.85% in 2020 (2, 3). The prevalence of ASD among children aged 6–12 years is ~0.70% in China (4). The large cost associated, mainly consisting of special education services and parental productivity loss, has caused a heavy burden to society and families (5). Despite the lack of effective drug treatments, several studies highlight the potential benefits of early diagnosis and parent-mediated interventions, which have to some extent improved children's social and communicative abilities (6, 7). The detection rate of genetic etiology of ASD is about 10–15% (8–10). However, early diagnosis remains a challenge for non-genetic ASD, which is mainly based on combining clinician observation with caregiver reports (11). Currently, there is an urgent need to find valuable biomarkers for the early diagnosis of ASD.

Some metabolomic studies have indicated the presence of elevated biomarkers in blood and urine samples from ASD patients, and these biomarkers include pyruvate, lactate, and mitochondrial-related enzymes (12–14). Major analytical techniques for metabolomics are nuclear magnetic resonance spectroscopy (NMR) and chromatography, including gas chromatography–mass spectrometry (GC-MS) and liquid chromatography–mass spectrometry (LC-MS). The advantages of ¹H-NMR are that the sample preprocessing is simple and non-destructive and that the detection of metabolites is comprehensive (15, 16). In 2010, Yap et al. first used urinary ¹H-NMR analysis to detect potential biomarkers for ASD (17). Over the last few years, with urinary ¹H-NMR analysis, some discriminating metabolites have been identified between ASD and healthy people. Studies indicated that ASD group showed higher levels of hippurate, glycine, tryptophan and D-threitol and lower levels of glutamate, creatine, valine, betaine and 3-methylhistidine. Further analyses indicated possible pathogenic mechanisms involving gut microbial metabolism, oxidative stress conditions and amino acid metabolism (18, 19). Overall, ¹H-NMR analysis shows great potential for the identification of biochemical signatures and investigation of the disease mechanisms of ASD. However, previous ¹H-NMR analyses lack large-scale sample sizes to confirm the significance and connection between metabolites and ASD. We aimed to conduct a large-scale and multi-center study to identify metabolic variations between ASD and control groups through urinary ¹H-NMR metabolomics analysis and to investigate potential biological mechanisms related to differential metabolites.

Methods

Participants

The study was conducted from January 2014 to December 2016. All the participants were recruited from research centers of 7 cities (Shanghai, Guangzhou, Changsha, Chongqing, Chengdu, Wenzhou, and Beijing). Participants were drawn from ASD and control group.

Autistic children from both hospitals and local autism rehabilitation of each research center were enrolled in the ASD group. The inclusion criteria for ASD group were: (a) children aged 2–18 years; (b) no limitation on the gender; (c) the diagnosis of ASD was based on the Diagnostic and Statistical Manual of Mental Disorders 5 (DSM-5) criteria (1) and confirmed with the Autism Diagnostic Observation Schedule (ADOS) and the Autism Diagnostic Interview-Revised (ADI-R) criteria by trained clinical psychiatrists from each research center; (d) urine sample was available. Exclusion criteria for ASD group were: (a) symptomatic autism (such as Rett syndrome and fragile X syndrome); (b) other mental illness (such as attention-deficit hyperactivity disorder); (c) other neurological disorders (such as epilepsy and central nervous system infections); (d) inherited metabolic diseases; (e) history of brain injury; (f) taking non-essential drug or dietary supplement before (72 h) and during sample collection.

Healthy children from the health examination center of each research center were enrolled in the control group. The inclusion criteria for control group were: (a) children with no abnormality in health examination and typical development; (b) age- and sex-matching with ASD group; (c) urine sample was available. Exclusion criteria for control group were: (a) siblings of ASD group; (b) clinical evidence of ASD diagnosis; (c) mental illness (such as attention-deficit hyperactivity disorder); (d) neurological disorders (such as epilepsy and central nervous system infections); (e) inherited metabolic diseases; (f) history of brain injury; (g) taking non-essential drug or dietary supplement before (72 h) and during sample collection.

After all eligible participants and their parents provided informed consent, they were invited to participate in the study. The study was approved by the institutional ethics committee at the Children's Hospital of Fudan University.

Sample Collection

First morning urine specimens were collected from all participants during the research period. Samples were collected in 15 mL urine collection tubes without preservative. Each sample was centrifuged and aliquoted into 1.5 mL EP tubes. Afterwards, samples were numbered (“1” represents the ASD group, and “2” represents the control group) and stored at −80°C immediately until ¹H-NMR analysis.

¹H-NMR Spectroscopy Experiments and Data Processing

Sample Preparation

A 500 μL urine sample was added to a 1.5 mL EP tube prefilled with 14 μL KF (5 M) solution. After vortexing, the sample was allowed to rest for 10 min, followed by centrifugation (12,000 rpm, 4°C) for 10 min. A total of 450 μL liquid supernatant was added to an NMR tube preloaded with 10 μL EDTA-d12 (0.1 M). The NMR tube cap was covered and mixed by hand. Finally, 45 μL Na⁺/K⁺ buffer (1.5 M, pH = 7.40) was added to the NMR tube, which was mixed by hand and then placed in the NMR spectrometer for data collection. The sample preparation process for ¹H-NMR analysis is shown in Supplementary Figure 1.

¹H-NMR Spectroscopy Experiments

All ¹H-NMR spectroscopy experiments were performed at 298 K using a Bruker AVIII 600 MHz NMR spectrometer (Bruker BioSpin, Germany) with a proton resonance frequency of 600.13 MHz.

The NOESYGPPR1D pulse sequence (RD-90°-t1-90°-tm-90°-ACQ) was used to collect the spectra. The parameters were set as follows: spectral width (SW) was 20 ppm, recycle delay (RD) was 2 s, mixing time (t_m) was 80 ms, t₁ was 4 μs, 90° pulse length was 14.8 s, data time was 1.36 s, data points were 32 K, and free induction decay (FID) accumulation was 64 times.

Data Processing

Spectra were processed using MestReNova software (MestReNova 8.1, Spain). The ¹H-NMR FID signals were multiplied by an exponential function equivalent to a line broadening of 1 Hz before performing an automatic Fourier transformation. The phase distortion and baseline of each spectrum were manually adjusted. The internal standard trimethylsilyl propanoic acid (TSP, δ = 0 ppm) was used as the baseline to calibrate the chemical shifts. The concentration of TSP was 0.261 mmol/L, and the spectrum of TSP represented 12 ¹H. The regions of the ¹H-NMR spectra (δ 0.3–9.5 ppm) were divided into consecutive integrated spectral regions of equal width (δ 0.004 ppm). The spectral region of the water (δ 4.71–5.055 ppm) and urea (δ 5.6–6.12 ppm) peaks were removed from each spectrum to minimize variations caused by the presaturation of the residual water and urea resonances. Mnova software was used to correct the spectra with obvious chemical shifts after the integration. Metabolites were assigned by referencing the values for chemical shifts in J-resolved (JRES), COZY, TOCSY, HSQC, and HMBC spectra and literature reports. A series of 2D NMR spectra were acquired for selected samples.

To eliminate the instrument differences of sensitivity and stability and to reduce the analysis errors caused by the concentration differences of the samples, two normalization methods were performed. Creatinine normalization: the ¹H-NMR spectra were normalized by using the creatinine methylene resonance (δ = 4.05 ppm) as a reference. Total area normalization (20): the integrated area in each bucket was normalized by the total sum of peak intensities to eliminate the effects of variable concentration among different samples.

Data and Statistical Analysis

Clinical Characteristics of Participants

Difference in age between the ASD and control groups was evaluated by Student's t-test when the distribution was normal or the Mann-Whitney U test when it was skewed. Difference in sex was investigated using chi-square test. Statistical analyses were performed by using the SPSS statistical package program (version 20, SPSS Inc., Chicago, IL, USA), and P < 0.05 was considered statistically significant.

Multivariate Analysis

The normalized data were imported into the SIMCA-P+ software package (version 13.0, Umetrics, Sweden) for multivariate statistical analysis, including principal component analysis (PCA) and orthogonal projection to latent structure discriminant analysis (OPLS-DA) (21–24). PCA was first used to observe the overall distribution among samples and the stability of the whole analysis process. Abnormal data, which related to sample contamination or improper sampling, were removed according to the overall aggregation trend in all the samples. After which OPLS-DA was used to distinguish the overall difference in the metabolic profile and to find differential metabolites between the groups.

To prevent model overfitting, an internal validation method was used to verify the validity of the model. The OPLS-DA model was validated by a 7-fold cross-validation (CV) (25). R² and Q² were two parameters to assess the quality of the model. OPLS-DA was further validated by variance analysis of the cross-validated residuals (CV-ANOVA) (26), and the model was considered valid at P < 0.05.

After multiplying the loading value of each variable with its standard deviation, backtracking conversion of the data was performed. Then, the data were assessed by multivariate analysis and imported into mapping software based on MATLAB (version 7.1, USA) to plot the loading diagram of the correlation coefficient. The Pearson correlation coefficient represented the linear correlation between the variable and the first principal component score value of the OPLS-DA model and was used to determine whether the variation contributed significantly to the intergroup differentiation. The significance was evaluated by the threshold value of the absolute correlation coefficient, which was determined according to the confidence interval of the sample size. If the absolute value of the Pearson correlation coefficient of the variable was higher than the threshold value (P < 0.05), the content of the variable was considered significantly different between groups.

Univariate Analysis

Differential metabolites between the two groups were selected for univariate analysis. Data from the total area normalized peak area are expressed as the mean ± SD. An independent sample t-test was used for comparing the two groups when the distribution was normal. The non-normal distribution data were evaluated by Mann-Whitney U test. P < 0.05 was considered statistically significant. The false discovery rate (FDR) was used to correct multiple hypothesis testing. FDR P < 0.1 was considered statistically significant. We applied correlation analysis to detect the relationship between differential metabolite levels and age.

Receiver Operating Characteristics (ROC) Curve Analysis

The sensitivity and specificity of metabolites with significant differences between the two groups in the diagnosis of ASD were evaluated using a ROC curve analysis. ROC curves were generated by using MedCalc statistical software (version 19.1.7, Belgium). The area under the curve (AUC) was used to measure the overall degree of identification power. An AUC > 0.7 was considered acceptable. Optimal cut-off points were determined by maximizing the Youden's J index (J = sensitivity + specificity – 1). The sensitivity, specificity, false-negative rate (FNR), false-positive rate (FPR), positive predictive value (PPV), negative predictive value (NPV), likelihood ratio positive (LR+), and likelihood ratio negative (LR–) were calculated to compare the diagnostic accuracy of the metabolites. Further analysis was performed after each group was stratified by sex and age. Logistic regression was used to analyse the combined metabolites.

Metabolic Pathways and Network Analysis

MetaboAnalyst (version 4.0, https://www.metaboanalyst.ca/) was used for pathway and network analyses (27). All differential metabolites were imported into the pathway analysis module to obtain matched pathways according to the P-values from the pathway enrichment analysis and pathway impact values from the pathway topology analysis. Significantly affected metabolic pathways associated with differential metabolites were identified by metabolite set enrichment analysis (MSEA). A metabolite-gene-disease interaction network was established to detect the connections of differential metabolites and associated pathways.

Results

Characteristics of the ASD and Control Groups

A total of 117 children with ASD were enrolled in the ASD group, and 119 healthy children were enrolled in the control group. The mean ages of the ASD and control groups were 10.12 ± 2.60 and 9.90 ± 1.73, respectively, with no significant difference between them (P = 0.468). The male to female ratios were 4.82:1 in the ASD group and 3.58:1 in the control group, with no significant difference between them (P = 0.388).

¹H-NMR Spectrum of Urine Samples

¹H-NMR spectra of urine samples from all participants were collected. A typical ¹H-NMR spectrum is shown in Figure 1. The keys for metabolites in Figure 1 are given in Supplementary Table 1. A total of 39 metabolites were identified in the ¹H NMR spectra of urine samples.

FIGURE 1

Figure 1. A typical 600 MHz ¹H-NMR spectrum of a urine sample. Region of dashed box is vertically magnified 8 times. Metabolites corresponding to each number are listed in Supplementary Table 1.

Multivariate Analysis of ¹H-NMR Spectra of Urine Samples

Creatinine Normalization Analysis

After creatinine normalization, PCA was performed on ¹H-NMR spectra of urine samples. Discriminant variables obtained from the PCA score plot (Figure 2A) were R²X = 0.403 and Q² = 0.135. The OPLS-DA score plot (Figure 2B) showed R²X = 0.257 and Q² = 0.0138. CV-ANOVA of the OPLS-DA model indicated that the model was not valid (P = 0.523).

FIGURE 2

Figure 2. After the creatinine normalization, multivariate analysis results of ¹H-NMR spectra of urine samples from the ASD group (green dots, 1) and control group (blue squares, 2). (A) PCA score plot: R²X = 0 403, Q² = 0.135. (B) OPLS-DA score plot: R²X = 0.257, Q² = 0.0138, CV-ANOVA P = 0.523.

Total Area Normalization Analysis

The ¹H-NMR spectra of the two groups were analyzed by PCA after total area normalization. The parameters of the PCA score plot (Figure 3A) were R²X = 0.366 and Q² = 0.0656. The OPLS-DA score plot (Figure 3B) showed R²X = 0.274, Q² = 0.0565 and CV-ANOVA P = 0.009. These results suggested that the model was valid.

FIGURE 3

Figure 3. After the total area normalization, multivariate analysis results of ¹H-NMR spectra of urine samples from the ASD group (green dots, 1) and control group (blue squares, 2). (A) PCA score plot: R²X = 0.366, Q² = 0.0656. (B) OPLS-DA score plot: R²X = 0.274, Q² = 0.0565, CV-ANOVA P = 0.009.

Differential Metabolites Between the ASD and Control Groups

To distinguish the metabolic differences between the two groups, a polychromatic correlation coefficient loading plot (Figure 4) was drawn. The color of the polychromatic loading plot was encoded by the absolute value of the correlation coefficient. The warmer the color is, the higher the absolute value of the correlation coefficient and the greater the contribution to the intergroup differentiation. The threshold of the absolute value of the Pearson's correlation coefficient was determined to be 0.182. The variables corresponding to the correlation coefficient with an absolute value > 0.182 contributed significantly to the intergroup differentiation (P < 0.05). Differential metabolites identified by the Pearson's correlation coefficient were showed in Table 1. Compared with the control group, the ASD group showed higher levels of glycine, guanidinoacetic acid, creatine, hydroxyphenylacetylglycine, phenylacetylglycine and formate and lower levels of 3-aminoisobutanoic acid, alanine, taurine, creatinine, hypoxanthine and N-methylnicotinamide.

FIGURE 4

Figure 4. OPLS-DA model: a polychromatic correlation coefficient loading plot.

TABLE 1

Table 1. Pearson's correlation coefficient of discriminant metabolites.

The fold changes of differential metabolites are summarized in Figure 5. A total of 6 metabolites, including 3-aminoisobutanoic acid (P = 0.0425), creatine (P = 0.0009), creatinine (P < 0.0001), hypoxanthine (P < 0.0001), formate (P = 0.0267), and N-methylnicotinamide (P = 0.0149), showed significant differences in the normalized peak areas between the two groups. The differences of above 6 metabolites remained significant after correction for multiple hypothesis testing (FDR P < 0.1) (Table 2). Correlation analysis showed that the level of creatinine was positively correlated with age in the ASD group (r = 0.215, 95% CI 0.013 to 0.401, P = 0.0370). In the control group, the level of creatinine was also positively correlated with age (r = 0.215, 95% CI 0.036 to 0.380, P = 0.0190), and creatine was negatively correlated with age (r = −0.277, 95% CI −0.435 to −0.102, P = 0.0023).

FIGURE 5

Figure 5. Metabolic pathway analyses utilizing the MetaboAnalyst functional interpretation tools. (A) Bubble plot of metabolic pathway impact. The metabolic pathways are shown as bubbles. The X coordinate and size of the bubble represent the value of pathway impact in the topology analysis. The Y coordinate and color of the bubble represent the P-value of the enrichment analysis. The darker red color and larger size indicate a more significant metabolite change in the corresponding pathway. (B) Metabolite set enrichment analysis (MSEA) plot. Significantly enriched pathways are represented by bars. The color and length of the bar are based on the P-value and fold enrichment, respectively.

TABLE 2

Table 2. Peak areas of differential metabolites after the total area normalization.

Differential Metabolites and Potential Biological Mechanisms

The Sensitivity and Specificity of Metabolites in the Diagnosis of ASD

The diagnostic accuracies of differential metabolites in the two groups were evaluated by ROC curve analysis (Table 3, Supplementary Figure 2). The ROC curve showed relatively significant diagnostic values of hypoxanthine (AUC = 0.657, 95% CI 0.588 to 0.726), creatinine (AUC = 0.639, 95% CI 0.569 to 0.709), creatine (AUC = 0.623, 95% CI 0.552 to 0.694), N-methylnicotinamide (AUC = 0.595, 95% CI 0.523 to 0.668) and guanidinoacetic acid (AUC = 0.574, 95% CI 0.501 to 0.647) for ASD (Supplementary Figure 3). The AUC of ROC analysis of the creatine/creatinine ratio was 0.6480 (95% CI 0.579 to 0.718). For each metabolite, there was no significant difference in AUCs between males and females (Supplementary Table 2). Compared with age stratification of 7–9 years old, the metabolites guanidinoacetic acid and creatine showed significantly higher diagnostic accuracy for ASD in the age stratification of 13–15 years old (AUC of guanidinoacetic acid = 0.802, 95% CI 0.566 to 0.944, P = 0.0282; AUC of creatine = 0.823, 95% CI 0.589 to 0.955, P = 0.0344. Results shown in Supplementary Table 3). By combining the metabolites of creatine, creatinine and hypoxanthine, the AUC of the ROC curve reached 0.720 (95% CI 0.659 to 0.777), with a sensitivity of 80.34% and specificity of 52.94%.

TABLE 3

Table 3. Diagnostic accuracies of differential metabolites between the ASD and control groups.

Correlated Metabolic Pathways and Networks

The main metabolic pathways associated with differential metabolites are shown in Supplementary Figure 4. According to the bubble plot of the metabolic pathway impact, there were significant metabolite changes in the glycine, serine and threonine metabolism, glyoxylate and dicarboxylate metabolism and taurine and hypotaurine metabolism pathways (Figure 6A). The plot of metabolite set enrichment analysis (MSEA) listed the significant enrichment pathways of differential metabolites (Figure 6B). Glycine, serine and threonine metabolism, primary bile acid biosynthesis and aminoacyl-tRNA biosynthesis were the pathways with the most significant enrichment. The metabolite-gene-disease interaction network provides a global view of the connection of the differential metabolites and the potential functional relationships among metabolites, connected genes, and target diseases (Supplementary Figure 5).

FIGURE 6

Figure 6. Comparation of the fold changes of the differential metabolites between the ASD and control groups. The fold change represents the ratio of average peak area between ASD group and control group.

Discussion

Complex etiologies and atypical symptoms pose very large challenges to the early diagnosis of ASD. Urinary ¹H-NMR analysis provides a fast and comprehensive assessment to detect potential biomarkers of ASD (28). Our study was based on the ¹H-NMR analysis of urine samples from ASD and control groups to identify candidate metabolites and associated pathogenic mechanisms. Compared with other studies of urinary ¹H-NMR metabolomics analysis of ASD, our study sample was large size, multi-center, and representative. Besides, our study used two normalization methods, creatinine normalization and total area normalization. The analysis using creatinine normalization showed no significant differences between the two groups, though the total area normalization did detect a difference. The variation of creatinine was confirmed by the total area normalization. Whiteley et al. also found that excretion of urinary creatinine in the group of pervasive developmental disorders, which included ASD, was significantly lower than controls (29). The abnormal creatinine metabolism might be caused by rigidity in food choice and various exclusion diets associated with ASD (30). Therefore, total area normalization was more suitable for our study.

The results obtained from the ¹H-NMR analysis revealed that the levels of glycine, guanidinoacetic acid, creatine, hydroxyphenylacetylglycine, phenylacetylglycine and formate were higher in the ASD group than those in the control group. Moreover, the levels of 3-aminoisobutanoic acid, alanine, taurine, creatinine, hypoxanthine and N-methylnicotinamide were lower in the ASD group than those in the control group. The normalized peak areas of 3-aminoisobutanoic acid, creatine, creatinine, hypoxanthine, formate and N-methylnicotinamide differed significantly between the two groups. The metabolite levels of glycine, taurine, creatine and creatinine were consistent with those from previous reports (19, 31). Our study is the first to detect variations in the metabolites of hydroxyphenylacetylglycine in urine samples from children with ASD. The results of ROC analysis indicated that creatine, creatinine and hypoxanthine have the potential to be biomarkers for the diagnosis of ASD. The creatine/creatinine ratio slightly improved the diagnostic accuracy. Recent research has reported significant female-related alterations of creatine and creatinine, and the creatinine/creatine ratio might be a good predictor of ASD in female subjects (32). We further found that combining the metabolites creatine, creatinine and hypoxanthine as a potential diagnostic indicator can largely improve the diagnostic accuracy for ASD.

The extension study of the metabolic pathway analysis demonstrated a possible imbalance of amino acid metabolism in ASD children. Differential metabolites between the ASD and control groups involved glycine, serine and threonine metabolism, arginine and proline metabolism, taurine and hypotaurine metabolism, and glutathione metabolism pathways (Supplementary Figure 6). In accordance with previous observations, amino acid metabolism disorder plays an important role in the pathogenesis mechanism of ASD (33–37). Creatine and creatinine, which show significant metabolism alterations in ASD, play an essential role in maintaining a high level of energy supply for the brain (38). Studies have indicated that creatine deficiency occurs in some ASD cases, and creatine may be a therapeutic target for ASD (38–41). Abnormalities in creatinine might closely linked to abnormalities in creatine. Creatine is biosynthesized from glycine and arginine with an intermediate metabolite of guanidinoacetic acid. Glycine acts as an excitatory neurotransmitter in the early developmental stage. As the nervous system matures, it transforms to the major inhibitory neurotransmitter. If the transformation does not occur, an abnormal level of glycine may result in neural disorders, including ASD (42–44). The taurine and hypotaurine metabolism pathway also differed significantly between the two groups. Oxidative stress imbalance is considered to be important components of the pathophysiology of ASD, and antioxidant therapy may improve the prognosis of ASD (45, 46). Park et al. reported that taurine, as an antioxidant and regulator of inflammation, might be a valid biomarker for ASD (47). Combined with vitamin D3, taurine showed benefits in the treatment of ASD (48, 49). All metabolic pathways interact with each other and constitute a complex network with related genes and diseases (Supplementary Figure 5).

Many studies have reported that ASD is associated with abnormal gut microbial metabolism (50–52). Gut microbiota metabolites, including phenylalanine, tyrosine, hippurate and tryptophan, have been reported to be factors in the development of ASD (53–56). In our study, phenylacetylglycine, a gut microbial co-metabolite, had a slight variation between the two groups. Phenylacetylglycine is the end product of the phenylalanine metabolism pathway (57). However, there are no studies on the relationship between phenylacetylglycine and the pathogenesis of ASD.

Limitations

There are some limitations in our study. At the cellular level, biochemical processes such as oxidative phosphorylation, redox reaction, and oxidative stress are regulated by circadian rhythms (58). About 50% of metabolites are thought to be rhythmic, involving metabolic pathways such as nucleotides, energy, oxidation, and carbohydrates (59). Studies indicated that disrupted circadian rhythm is closely related to ASD (60–62). Thus, in our study, inconsistent urine collection times may have affected the quality of metabolic analyses due to the restriction of research centers. A few small-scale sample studies have set healthy sibling group as control to remove the effects of confounding factors, such as heredity and environment (17, 19). However, our study lacks of healthy siblings as controls.

Though ¹H-NMR analysis is the regular NMR method, some studies have used 2D HSQC-NMR to improve urinary screening in ASD. Compared to ¹H-NMR analysis, ¹H-¹³C HSQC-NMR analysis shows the advantage of improving the metabolite detection accuracy and the discrimination ability (18, 31). Moreover, it is necessary to compare metabolite levels that vary with ASD severity to better clarify the pathogenesis of ASD.

Conclusions

In our study, ¹H-NMR metabolomic analysis was used to investigate urinary metabolism patterns in the ASD and control groups. We revealed that urinary amino acid metabolites were significantly altered in children with ASD. A series of variations in amino acid metabolism pathways, including glycine, serine and threonine metabolism, arginine and proline metabolism, and taurine and hypotaurine metabolism might play important roles in the pathogenic mechanisms of ASD. Further studies of differential metabolites are needed to improve the understanding of ASD pathogenesis.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the institutional ethics committee at the Children's Hospital of Fudan University. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.

Author Contributions

YW conceived of the study. YM and HZ contributed to the analysis, synthesis and interpretation of the results, and wrote the manuscript. CL, LZ, and TW contributed to the sample collection. XZ, XL, LW, TL, XC, MM, YH, and EL contributed to the diagnosis of ASD and sample collection from each research center. YA, XX, WY, and YJ provided guidance for the study. All of the authors contributed to the preparation of the manuscript.

Funding

This project was supported by the National Health Commission of the People's Republic of China (201302002, Clinical Trials NCT02200679), Shanghai Hospital Development Center (SHDC12015113), and Guangdong Key Project in Development of new tools for diagnosis and treatment of Autism (2018B030335001).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We sincerely thank all the patients, controls and parents who provided written consent to participate in this study, as well as the clinicians who contributed to the recruitment of the participants.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2021.624767/full#supplementary-material

Supplementary Figure 1. Urine sample preparation process.

Supplementary Figure 2. ROC curve of the metabolites. AUC, area under the ROC curve; SEN, sensitivity; SPE, specificity.

Supplementary Figure 3. ROC curve of the metabolites with relatively significant diagnostic values.

Supplementary Figure 4. Metabolic pathways associated with the differential metabolites.

Supplementary Figure 5. Metabolite-gene-disease interaction network.

Supplementary Figure 6. Differential metabolites associated with amino acid metabolism disorders.

Supplementary Table 1. Main metabolite attributions in the ¹H-NMR spectra of urine samples.

Supplementary Table 2. Sex stratification analysis of the AUC of differential metabolites.

Supplementary Table 3. Age stratification analysis of the AUC of differential metabolites.

Abbreviations

ASD, Autism spectrum disorder; PCA, Principal component analysis; OPLS-DA, Orthogonal projection to latent structure discriminant analysis; NMR, Nuclear magnetic resonance spectroscopy; DSM-5, Statistical manual of mental disorders 5; ADOS, Autism diagnostic observation schedule; ADI-R, Autism diagnostic interview-revised; SW, Spectral width; RD, Recycle delay; FID, Free induction decay; CV, Cross-validation; CV-ANOA, Variance analysis of the cross-validated residuals; ROC, Receiver operating characteristics; AUC, Area under the curve; FNR, False-negative rate; FPR, False-positive rate; PPV, Positive predictive value; NPV, Negative predictive value; LR+, Likelihood ratio positive; LR–, Likelihood ratio negative; MSEA, Metabolite sets enrichment analysis.

References

1. Lai MC, Lombardo MV, Baron-Cohen S. Autism. Lancet. (2014) 383:896–910. doi: 10.1016/S0140-6736(13)61539-1

CrossRef Full Text | Google Scholar

2. Baron-Cohen S, Scott FJ, Allison C, Williams J, Bolton P, Matthews FE, et al. Prevalence of autism-spectrum conditions: UK school-based population study. Br J Psychiatry. (2009) 194:500–9. doi: 10.1192/bjp.bp.108.059345

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Maenner MJ, Shaw KA, Baio J, Washington A, Patrick M, DiRienzo M, et al. Prevalence of autism spectrum disorder among children aged 8 years - autism and developmental disabilities monitoring network, 11 sites, United States, 2016. MMWR Surveill Summ. (2020) 69:1–12. doi: 10.15585/mmwr.ss6904a1

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Zhou H, Xu X, Yan W, Zou X, Wu L, Luo X, et al. Prevalence of autism spectrum disorder in China: a nationwide multi-center population-based study among children aged 6 to 12 years. Neurosci Bull. (2020) 36:961–71. doi: 10.1007/s12264-020-00530-6

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Buescher AV, Cidav Z, Knapp M, Mandell DS. Costs of autism spectrum disorders in the United Kingdom and the United States. JAMA Pediatr. (2014) 168:721–8. doi: 10.1001/jamapediatrics.2014.210

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Sanchack KE, Thomas CA. Autism spectrum disorder: primary care principles. Am Fam Physician. (2016) 94:972–9.

PubMed Abstract | Google Scholar

7. Tachibana Y, Miyazaki C, Ota E, Mori R, Hwang Y, Kobayashi E, et al. A systematic review and meta-analysis of comprehensive interventions for pre-school children with autism spectrum disorder (ASD). PLoS ONE. (2017) 12:e0186502. doi: 10.1371/journal.pone.0186502

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Tsiplova K, Zur RM, Marshall CR, Stavropoulos DJ, Pereira SL, Merico D, et al. A microcosting and cost-consequence analysis of clinical genomic testing strategies in autism spectrum disorder. Genet Med. (2017) 19:1268–75. doi: 10.1038/gim.2017.47

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Jiang YH, Yuen RK, Jin X, Wang M, Chen N, Wu X, et al. Detection of clinically relevant genetic variants in autism spectrum disorder by whole-genome sequencing. Am J Hum Genet. (2013) 93:249–63. doi: 10.1016/j.ajhg.2013.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Tammimies K, Marshall CR, Walker S, Kaur G, Thiruvahindrapuram B, Lionel AC, et al. Molecular diagnostic yield of chromosomal microarray analysis and whole-exome sequencing in children with autism spectrum disorder. JAMA. (2015) 314:895–903. doi: 10.1001/jama.2015.10078

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Lord C, Elsabbagh M, Baird G, Veenstra-Vanderweele J. Autism spectrum disorder. Lancet. (2018) 392:508–20. doi: 10.1016/S0140-6736(18)31129-2

CrossRef Full Text | Google Scholar

12. Khemakhem AM, Frye RE, El-Ansary A, Al-Ayadhi L, Bacha AB. Novel biomarkers of metabolic dysfunction is autism spectrum disorder: potential for biological diagnostic markers. Metab Brain Dis. (2017) 32:1983–97. doi: 10.1007/s11011-017-0085-2

PubMed Abstract | CrossRef Full Text | Google Scholar

13. El-Ansary A, Bjørklund G, Khemakhem AM, Al-Ayadhi L, Chirumbolo S, Ben Bacha A. Metabolism-associated markers and childhood autism rating scales (CARS) as a measure of autism severity. J Mol Neurosci. (2018) 65:265–76. doi: 10.1007/s12031-018-1091-5

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Hollis F, Kanellopoulos AK, Bagni C. Mitochondrial dysfunction in autism spectrum disorder: clinical features and perspectives. Curr Opin Neurobiol. (2017) 45:178–87. doi: 10.1016/j.conb.2017.05.018

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Ryan D, Robards K, Prenzler PD, Kendall M. Recent and potential developments in the analysis of urine: a review. Anal Chim Acta. (2011) 684:8–20. doi: 10.1016/j.aca.2010.10.035

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Kim HK, Choi YH, Verpoorte R. NMR-based plant metabolomics: where do we stand, where do we go? Trends Biotechnol. (2011) 29:267–75. doi: 10.1016/j.tibtech.2011.02.001

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Yap IK, Angley M, Veselkov KA, Holmes E, Lindon JC, Nicholson JK. Urinary metabolic phenotyping differentiates children with autism from their unaffected siblings and age-matched controls. J Proteome Res. (2010) 9:2996–3004. doi: 10.1021/pr901188e

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Mavel S, Nadal-Desbarats L, Blasco H, Bonnet-Brilhault F, Barthélémy C, Montigny F, et al. ¹H-¹³C NMR-based urine metabolic profiling in autism spectrum disorders. Talanta. (2013) 114:95–102. doi: 10.1016/j.talanta.2013.03.064

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Lussu M, Noto A, Masili A, Rinaldi AC, Dessì A, De Angelis M, et al. The urinary ¹H-NMR metabolomics profile of an italian autistic children population and their unaffected siblings. Autism Res. (2017) 10:1058–66. doi: 10.1002/aur.1748

CrossRef Full Text | Google Scholar

20. Craig A, Cloarec O, Holmes E, Nicholson JK, Lindon JC. Scaling and normalization effects in NMR spectroscopic metabonomic data sets. Anal Chem. (2006) 78:2262–7. doi: 10.1021/ac0519312

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Worley B, Powers R. PCA as a practical indicator of OPLS-DA model reliability. Curr Metabolomics. (2016) 4:97–103. doi: 10.2174/2213235X04666160613122429

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Souihi N, Lindegren A, Eriksson L, Trygg J. OPLS in batch monitoring - Opens up new opportunities. Anal Chim Acta. (2015) 857:28–38. doi: 10.1016/j.aca.2014.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Pinto RC. Chemometrics methods and strategies in metabolomics. Adv Exp Med Biol. (2017) 965:163–90. doi: 10.1007/978-3-319-47656-8_7

CrossRef Full Text | Google Scholar

24. Worley B, Powers R. Multivariate analysis in metabolomics. Curr Metabolomics. (2013) 1:92–107. doi: 10.2174/2213235X11301010092

CrossRef Full Text | Google Scholar

25. Triba MN, Le Moyec L, Amathieu R, Goossens C, Bouchemal N, Nahon P, et al. PLS/OPLS models in metabolomics: the impact of permutation of dataset rows on the K-fold cross-validation quality parameters. Mol Biosyst. (2015) 11:13–9. doi: 10.1039/C4MB00414K

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Røraas T, Støve B, Petersen PH, Sandberg S. Biological variation: the effect of different distributions on estimated within-person variation and reference change values. Clin Chem. (2016) 62:725–36. doi: 10.1373/clinchem.2015.252296

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Chong J, Wishart DS, Xia J. Using metaboanalyst 4.0 for comprehensive and integrative metabolomics data analysis. Curr Protoc Bioinformatics. (2019) 68:e86. doi: 10.1002/cpbi.86

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Silva RA, Pereira TCS, Souza AR, Ribeiro PR. ¹H NMR-based metabolite profiling for biomarker identification. Clin Chim Acta. (2020) 502:269–79. doi: 10.1016/j.cca.2019.11.015

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Whiteley P, Waring R, Williams L, Klovrza L, Nolan F, Smith S, et al. Spot urinary creatinine excretion in pervasive developmental disorders. Pediatr Int. (2006) 48:292–7. doi: 10.1111/j.1442-200X.2006.02207.x

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Bandini LG, Curtin C, Phillips S, Anderson SE, Maslin M, Must A. Changes in food selectivity in children with autism spectrum disorder. J Autism Dev Disord. (2017) 47:439–46. doi: 10.1007/s10803-016-2963-6

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Nadal-Desbarats L, Aïdoud N, Emond P, Blasco H, Filipiak I, Sarda P, et al. Combined 1H-NMR and 1H-13C HSQC-NMR to improve urinary screening in autism spectrum disorders. Analyst. (2014) 139:3460–8. doi: 10.1039/C4AN00552J

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Xiong X, Liu D, He W, Sheng X, Zhou W, Xie D, et al. Identification of gender-related metabolic disturbances in autism spectrum disorders using urinary metabolomics. Int J Biochem Cell Biol. (2019) 115:105594. doi: 10.1016/j.biocel.2019.105594

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Smith AM, King JJ, West PR, Ludwig MA, Donley ELR, Burrier RE, et al. Amino acid dysregulation metabotypes: potential biomarkers for diagnosis and individualized treatment for subtypes of autism spectrum disorder. Biol Psychiatry. (2019) 85:345–54. doi: 10.1016/j.biopsych.2018.08.016

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Liu A, Zhou W, Qu L, He F, Wang H, Wang Y, et al. Altered urinary amino acids in children with autism spectrum disorders. Front Cell Neurosci. (2019) 13:7. doi: 10.3389/fncel.2019.00007

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Li C, Shen K, Chu L, Liu P, Song Y, Kang X. Decreased levels of urinary free amino acids in children with autism spectrum disorder. J Clin Neurosci. (2018) 54:45–49. doi: 10.1016/j.jocn.2018.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Delaye JB, Patin F, Lagrue E, Le Tilly O, Bruno C, Vuillaume ML, et al. Post hoc analysis of plasma amino acid profiles: towards a specific pattern in autism spectrum disorder and intellectual disability. Ann Clin Biochem. (2018) 55:543–52. doi: 10.1177/0004563218760351

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Ming X, Stein TP, Barnes V, Rhodes N, Guo L. Metabolic perturbance in autism spectrum disorders: a metabolomics study. J Proteome Res. (2012) 11:5856–62. doi: 10.1021/pr300910n

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Fons C, Campistol J. Creatine defects and central nervous system. Semin Pediatr Neurol. (2016) 23:285–9. doi: 10.1016/j.spen.2016.11.003

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Schulze A, Bauman M, Tsai AC, Reynolds A, Roberts W, Anagnostou E, et al. Prevalence of creatine deficiency syndromes in children with nonsyndromic autism. Pediatrics. (2016) 137:e20152672. doi: 10.1542/peds.2015-2672

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Schulze A. Creatine deficiency syndromes. Handb Clin Neurol. (2013) 113:1837–43. doi: 10.1016/B978-0-444-59565-2.00053-8

CrossRef Full Text | Google Scholar

41. Aydin HI. Creatine transporter deficiency in two brothers with autism spectrum disorder. Indian Pediatr. (2018) 55:67–8. doi: 10.1007/s13312-018-1232-5

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Kaila K, Price TJ, Payne JA, Puskarjov M, Voipio J. Cation-chloride cotransporters in neuronal development, plasticity and disease. Nat Rev Neurosci. (2014) 15:637–54. doi: 10.1038/nrn3819

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Ito S. GABA and glycine in the developing brain. J Physiol Sci. (2016) 66:375–9. doi: 10.1007/s12576-016-0442-7

CrossRef Full Text | Google Scholar

44. Tyzio R, Nardou R, Ferrari DC, Tsintsadze T, Shahrokhi A, Eftekhari S, et al. Oxytocin-mediated GABA inhibition during delivery attenuates autism pathogenesis in rodent offspring. Science. (2014) 343:675–9. doi: 10.1126/science.1247190

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Hu T, Dong Y, He C, Zhao M, He Q. The gut microbiota and oxidative stress in autism spectrum disorders (ASD). Oxid Med Cell Longev. (2020) 2020:8396708. doi: 10.1155/2020/8396708

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Bjørklund G, Meguid NA, El-Bana MA, Tinkov AA, Saad K, Dadar M, et al. Oxidative stress in autism spectrum disorder. Mol Neurobiol. (2020) 57:2314–32. doi: 10.1007/s12035-019-01742-2

CrossRef Full Text | Google Scholar

47. Park E, Cohen I, Gonzalez M, Castellano MR, Flory M, Jenkins EC, et al. Is taurine a biomarker in autistic spectrum disorder? Adv Exp Med Biol. (2017) 975(Pt 1):3–16. doi: 10.1007/978-94-024-1079-2_1

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Omura Y, Lu D, Jones MK, Nihrane A, Duvvi H, Shimotsuura Y, et al. Early detection of autism (ASD) by a non-invasive quick measurement of markedly reduced acetylcholine & DHEA and increased β-amyloid (1-42), asbestos (Chrysotile), titanium dioxide, Al, Hg & often coexisting virus infections (CMV, HPV 16 and 18), bacterial infections etc. in the brain and corresponding safe individualized effective treatment. Acupunct Electrother Res. (2015) 40:157–87. doi: 10.3727/036012915X14473562232941

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Zheng HF, Wang WQ, Li XM, Rauw G, Baker GB. Body fluid levels of neuroactive amino acids in autism spectrum disorders: a review of the literature. Amino Acids. (2017) 49:57–65. doi: 10.1007/s00726-016-2332-y

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Fattorusso A, Di Genova L, Dell'Isola GB, Mencaroni E, Esposito S. Autism spectrum disorders and the gut microbiota. Nutrients. (2019) 11:521. doi: 10.3390/nu11030521

CrossRef Full Text | Google Scholar

51. Ding HT, Taur Y, Walkup JT. Gut microbiota and autism: key concepts and findings. J Autism Dev Disord. (2017) 47:480–9. doi: 10.1007/s10803-016-2960-9

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Bezawada N, Phang TH, Hold GL, Hansen R. Autism spectrum disorder and the gut microbiota in children: a systematic review. Ann Nutr Metab. (2020) 76:16–29. doi: 10.1159/000505363

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Clayton TA. Metabolic differences underlying two distinct rat urinary phenotypes, a suggested role for gut microbial metabolism of phenylalanine and a possible connection to autism. FEBS Lett. (2012) 586:956–61. doi: 10.1016/j.febslet.2012.01.049

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Srikantha P, Mohajeri MH. The possible role of the microbiota-gut-brain-axis in autism spectrum disorder. Int J Mol Sci. (2019) 20:2115. doi: 10.3390/ijms20092115

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Sitkin SI, Tkachenko EI, Vakhitov TY. Metabolic dysbiosis of the gut microbiota and its biomarkers. Eksp Klin Gastroenterol. (2016) 12:6–29. doi: 10.18786/2072-0505-2015-40-12-34

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Agus A, Planchais J, Sokol H. Gut microbiota regulation of tryptophan metabolism in health and disease. Cell Host Microbe. (2018) 23:716–24. doi: 10.1016/j.chom.2018.05.003

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Delaney J, Neville WA, Swain A, Miles A, Leonard MS, Waterfield CJ. Phenylacetylglycine, a putative biomarker of phospholipidosis: its origins and relevance to phospholipid accumulation using amiodarone treated rats as a model. Biomarkers. (2004) 9:271–90. doi: 10.1080/13547500400018570

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Kim Y, Vadodaria KC, Lenkei Z, Kato T, Gage FH, Marchetto MC, et al. Mitochondria, metabolism, and redox mechanisms in psychiatric disorders. Antioxid Redox Signal. (2019) 31:275–317. doi: 10.1089/ars.2018.7606

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Krishnaiah SY, Wu G, Altman BJ, Growe J, Rhoades SD, Coldren F, et al. Clock regulation of metabolites reveals coupling between transcription and metabolism. Cell Metab. (2017) 25:961–74.e4. doi: 10.1016/j.cmet.2017.03.019

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Yenen AS, Çak HT. Melatonin and circadian rhythm in autism spectrum disorders. Turk Psikiyatri Derg. (2020) 31:201–11. doi: 10.5080/u25160

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Logan RW, McClung CA. Rhythms of life: circadian disruption and brain disorders across the lifespan. Nat Rev Neurosci. (2019) 20:49–65. doi: 10.1038/s41583-018-0088-y

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Pinato L, Galina SC, Markus RP, Da SCS. Dysregulation of circadian rhythms in autism spectrum disorders. Curr Pharm Des. (2019) 25:4379–93. doi: 10.2174/1381612825666191102170450

PubMed Abstract | CrossRef Full Text | Google Scholar