Using urine FTIR spectra to screen autism spectrum disorder

Autism spectrum disorder (ASD) is a heterogeneous neurodevelopmental disorder caused by multiple factors, lacking clear biomarkers. Diagnosing ASD still relies on behavioural and developmental signs and usually requires lengthy observation periods, all of which are demanding for both clinicians and parents. Although many studies have revealed valuable knowledge in this field, no clearly defined, practical, and widely acceptable diagnostic tool exists. In this study, 26 children with ASD (ASD+), aged 3–5 years, and 26 sex and age-matched controls are studied to investigate the diagnostic potential of the Attenuated Total Reflectance-Fourier Transform Infrared (ATR-FTIR) spectroscopy. The urine FTIR spectrum results show a downward trend in the 3000–2600/cm region for ASD+ children when compared to the typically developing (TD) children of the same age. The average area of this region is 25% less in ASD+ level 3 children, 29% less in ASD+ level 2 children, and 16% less in ASD+ level 1 children compared to that of the TD children. Principal component analysis was applied to the two groups using the entire spectrum window and five peaks were identified for further analysis. The correlation between the peaks and natural urine components is validated by artificial urine solutions. Less-than-normal levels of uric acid, phosphate groups, and ammonium (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${NH}_{4}^{+}$$\end{document}NH4+) can be listed as probable causes. This study shows that ATR-FTIR can serve as a practical and non-invasive method to screen ASD using the high-frequency region of the urine spectrum.


Sample collection
This study was approved by the Atilim University Human Studies Ethics Committee (Ref.Number: 59394181-604.01.01-8509).All experiments were performed in accordance with relevant guidelines and regulations.After the participants' parents signed an informed consent form, first-morning midstream urine samples were collected from 26 Turkish children aged 3-5 years newly diagnosed with ASD, as well as from 26 sex-, BMI-, and age-matched controls showing typical development after at least an 8 h fasting period.The parents answered a list of questions about their children regarding their general eating habits, food their children are obsessed with, and the use of any prescribed medication.Children who did not use any medications, have a healthy diet with no particular food obsessions, and who had toilet training were included in the study.Firstly, the urine samples were collected in clean, dry specimen tubes according to the urine collection guidelines.Then, a dipstick test (Mission Acon, San Diego, California) was used to screen the urinary parameters.Each urine sample was aliquoted into 2 mL tubes and stored at -20 °C until measured with ATR-FTIR.The BMI was calculated by weight/height 2 for all participants.
The Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) criteria were used to diagnose children with ASD.The same scale was also used to clarify the severity level of ASD, and a paediatric psychiatrist performed the clinical evaluation.In addition, the Autism Behaviour Checklist (ABC) as a tool for ASD assessment was also applied.ABC is a list of 57 questions assessing the behaviours and symptoms of ASD and rating their scale for children above 3 years old.The ABC consists of questions in five categories: Sensory, Relating, Body and Object Use, Language and Social, and Self-Help skills.
Typically developing sex-, BMI-, and age-matched children with no known health problems were enrolled to form the control group.All control participants were determined healthy based on a questionnaire, standard clinical evaluation, and urinary dipstick test.
Having another neuropsychiatric disorder, chronic diseases such as diabetes, hypertension, etc., abnormal urine test results, complaints related to the urinary system, and using any medication were determined as exclusion criteria.In addition, the presence of leucocytes, protein, blood, glucose, bilirubin, nitrite, and urine pH being less than 5 or more than 7, were considered abnormal after the urine test result.

Spectral data collection
All the spectra were collected using a Thermo Scientific Nicolet 6700 FTIR spectrometer with a diamond attenuated total reflection (ATR, ConcentratIR2, Harrick) accessory with ten internal reflections.A deuterated triglycine sulphate (DTGS) detector was used for measurements.First, the diamond surface was cleaned with ethanol and rinsed with water.It was then air-dried for 2 min before recording the background interferogram.A frozen urine sample was thawed at room temperature and stirred using a vortex.Five microliters of the urine sample were pipetted onto the diamond surface and dried for 15 min by a gentle stream of N 2 gas prior to data collection to remove the excess water.All samples were recorded in the 4000-600/cm wavenumber range, and each spectrum was obtained by averaging 128 interferograms for a final resolution of 4/cm.A total of 128 scans with the DTGS detector take about 2 min to complete.Each sample was deposited and measured three times to eliminate pipetting errors and experimental variations.The average of the triplet measurement from each sample was used for further analysis.The Thermo Scientific OMNIC Spectrometer software, version 8.2.388, was used for spectra collection.The spectrometer was constantly purged with N 2 during measurements to eliminate atmospheric variations in the water vapor.
For the peak analysis results shown in Figs. 4 and 6, integration limits and baseline points are listed in Table 1.Spectra were not pre-processed before integration; however, limits of integration were chosen to exclude the contribution of other urine components, and a suitable baseline was chosen to eliminate the baseline differences among spectra.In order to see the differences in the analysis results due to different pre-processing methods, the spectra were also area normalized before integrating the 3000-2600/cm region.Normalization to the urea band at 1605/cm was also tested in that regard.

Statistical analysis
All statistical analyses and visualizations were performed using OriginPro 2017 SR2 and SPPS statistics for Windows, version 22.0 (IBM Corp., NY, USA).Shapiro-Wilk test and Kolmogrov-Smirnov with histogram curves were used to investigate the normality of the distribution.The numerical data were presented as mean

Results and discussion
In this study, the infrared spectra of urine from 26 children diagnosed with autism spectrum disorder (ASD+) are compared to those of 26 healthy, typically developing children.Both ASD+ (n = 26) and control (n = 26) groups were the same in terms of gender.The age (42.2 ± 9.1 months in ASD + , 44.7 ± 8.3 months in controls; p = 0.585) and BMI (16.3 ± 2.9 kg/m 2 in ASD + and 15.7 ± 1.9 kg/m 2 , in controls; p = 0.285) were also statistically similar in these groups.Although it was not statistically significant, the parameter 'sign of mental illness in a relative' is slightly more in the ASD+ group (6/26, 23.1%) compared to that in controls (1/26, 3.8%) (p = 0.086).Family income is higher in controls than in ASD+, and it is different with statistical significance (p < 0.0001).The majority (14/26, 55.8%) of ASD+ patients' family income is in the range of 2500-5000 Turkish Lira (TL) ($300-$600), while 15 of 26 (57.7%) in controls are in 10,000-20,000 TL ($900-$1800) range.One-way ANOVA test revealed that only the parameter 'body and object use' was different statistically (p = 0.037).Post hoc analysis showed that this difference was due to the ASD Level 2 group.Other ABC scores were similar statistically.The data are presented briefly in Table 2.
We found that the family income was higher and significantly different in TD controls.The studies relating ASD to family income indicate that its prevalence is higher in families with high household incomes 41 .It is emphasized that the difference may also be due to the delay or lack of diagnosis in low-income households.In this study, firstly, the ASD+ group was formed, and then the TD children who were admitted paediatric psychiatry clinic of the tertiary university hospital were assigned to the control group.Thus, families who take their children to regular paediatric evaluations can be expected to be in a better socio-economic status since they are mostly the families with high awareness and education levels.Irrespective of the reason, family income levels may affect the diet habits, its content, educational status, and the environmental factors of children 42,43 .In addition, some of the possible mechanisms discussed in the literature about ASD are related to diet and microbiota 44  www.nature.com/scientificreports/see the effect of family income difference on the urine spectra, the spectra of TD children from income groups 1 and 5 are compared (Supporting Information Fig. S1).The spectral difference between these two groups is minute compared to that observed between ASD+ and TD group children.Possible reasons for this observation are discussed in Supporting Information.Therefore, presented spectral data indicate disease-related differences.

FTIR spectra of urine
Figure 1 shows all the collected spectra overlaid with the average spectrum of all ASD + children and the average of controls.As can be seen, the two average spectra are comparable in terms of the general spectral profile and prominent peak positions.Nevertheless, there are marked differences that are common among all pair-wise comparisons.The two average spectra are significantly different at the 0.05 level based on the Mann-Whitney U test (Z = -8.6,U = 5.48 × 10 6 ).
When individual spectra from the ASD+ group are pair-wise compared to their controls, it is observed that almost all the spectra from ASD+ children have lower absorbance at all points than those from their control children (data not shown).However, the deviation of each spectrum with respect to its control is not uniform, being high at some points as 1071 and 930/cm compared to others.To demonstrate the absorbance differences between the spectra of ASD+ and control children, the difference spectrum of each ASD+ participant from the average spectrum of controls is shown in Fig. 2. The positive peaks indicate higher absorbance with respect to the average control, whereas the negative peaks indicate lesser absorbance.In a previous study, it has been shown that the spectra of the healthy children population show the least intensity variations compared to other age groups 32 .In the same study, the urine spectra collected from 41 children aged 3-10 were shown to be similar despite the wide age range and the difference in gender, dietary habits, and their parents' socio-economic status.In this study, the spectra of the control children are also very similar in intensity (Fig. 1-gray lines).The intensities of the difference spectra are particularly low in the high wavenumber region for these children (Fig. 2-upper panel).On the other hand, the absorbance difference from the average spectrum of controls is substantial for ASD+ children on the same scale (Fig. 2-lower panel).It can be observed that most ASD+ urine spectra have lower absorbance than the average spectrum of controls.It can be due to the lack of some urine components excreted consistently by ASD+ children.Although the low-frequency region (1800-600/cm) shows this difference more www.nature.com/scientificreports/dramatically, in essence, the high-frequency region (4000-2400/cm) is more interesting since the control group has the least variation in the same region.Thus, the high-wavenumber region is further analysed in an attempt to quantify and explain the observed deviation.
The PCA is performed to determine any relationships among the two groups using the entire spectral region.Figure 3a and b present a two-dimensional PCA scatter plot and the loadings plot of PC1 and PC2 as a function of wavenumber for the 3800-600/cm region, respectively.The score plots present the distribution of ASD+ (red stars) and controls (black stars).Most of the controls are distinguished from ASD+ along the PC1 axis.The PC1 accounts for 87.3% while PC2 accounts for 5.8% of the total variance, both within 95% confidence ellipses.The scores of controls lie on the positive side of the PC1 scores, but most ASD+ are on the negative side.The ASD+ samples show a wider distribution than the controls.
The positive and negative peaks can be observed in the loading plots.The positive peaks indicate the positive score values due to increased absorbance in the spectra.In contrast, we can relate the negative loading values to the increased absorbance in the spectra of the negative score values 45 .A loading plot is used to visualize the relationship between variables and PCs extracted from a dataset.Although PCA does not show a net separation of ASD+ data from that of TD children, it identifies the bands and regions that contributed to the partial dissociation observed in the score plot.The absorbance differences at 781, 930, 1065, 2810, and 3019/cm contributed to clustering the two groups.While the peak at 781/cm can be mainly attributed to uric acid, others are affected

A closer look at the 3000-2600/cm region
In the high-frequency region, ASD+ children have lower absorbance with respect to their controls in the 3000-2600/cm region (Figs. 1 and 2).The absorbance difference is more dramatic in children with ASD severity levels 2 and 3.In order to quantify the absorbance difference, this region is integrated using raw spectra to find the area under the curve (Table 1).The results are tabulated in Table S1 (Supporting Information).The spectra of the ASD+ level 3 children show 24% less absorbance with respect to controls, while level 2 children show 29% less absorbance in the same region.Although the ASD+ level 1 children also show the same trend (16% less), the difference is less in comparison to levels 2 and 3. Figure 4 compares the mean of the peak area for all ASD levels using ANOVA with Fisher's LSD and Levene's test.The results show that the population variances are similar for all ASD levels, but the population means are significantly different (Fig. 4a).The difference between ASD+ and TD groups is significant for levels 2 and 3, except level 1 (Fig. 4b-d).
Interestingly, the average area of the 3000-2600/cm region in the TD children, compared to that of the ASD+ group at each level, is nearly the same (Table S1).This integrated region shows slight variance among healthy children, as previously shown in a study with 41 children aged 3-10 32 .In the same study, the two genders showed no difference within the same age group in the high-frequency region.
Two other common normalization methods in the literature, namely, the area normalization and the min-max normalization (1605/cm-urea normalization), are also tested for the same area calculation.The results of these   3.The area of the 3000-2600/cm region is consistently less in all ASD+ children with respect to their controls.Different normalization methods reduced, but did not eliminate in full, the absorbance difference between the two groups.Therefore, it can be suggested that possible urine concentration differences among the children are not a factor in the observed results.

Origin of the absorbance difference in the 3000-2600/cm region
The source(s) of this difference in absorbance between the ASD+ and TD children is not certain; however, based on the individual urine component measurements in previous studies 31 , it is possible to comment on the probable sources.Less-than-normal levels of uric acid, phosphate groups, ammonium ( NH + 4 ), or more-than-normal levels of urea and sulphate groups can be listed as probable causes.To test the contribution of these components, an artificial urine solution is prepared 31 .The spectrum of this artificial urine is compared to that of another solution, including twice the uric acid than normal, so that the spectral contribution of uric acid and its marker bands can be observed.The difference between these two spectra for the high-frequency region is given in Fig. 5.It can be seen that the increased amount of uric acid has a positive contribution to the 3300-2500/cm region of the spectrum.The same procedure is repeated for ammonium and phosphate separately, each of which component positively contributes to the absorbance in the same region.On the other hand, increased urea or sulphate decreases the absorbance of the same region and, thus, has a negative effect.

Examination of spectra in terms of ammonium, phosphate, uric acid, sulphate, and urea
These components, while affecting the absorbance in the high-frequency region, also have other peaks in the low-frequency region.In the subsequent sections, each of these components is discussed separately in terms of the levels detected in our cohort, their importance, and their possible relation with ASD using previous studies.For each urine component discussed below, the effect of dietary habits is not discussed since FTIR spectroscopy is less sensitive to slight changes in the urinary level of these components compared to biochemical analyses.Additionally, dietary choices are relatively healthier and less variant among children of age 3 to 5 compared to adults.In fact, the distribution of variation seen within the TD group can be attributed to different dietary habits.The degree of differences among the spectra of TD children is far less than the observed changes between the two groups in this study, which mostly rules out the effect of diet.

Phosphate group
The urinary phosphate levels can be followed by peaks at ~ 1050 and 930/cm.Upon second derivative calculations, the peak at 1050/cm is seen to overlap highly with those of the sulphate groups and the urea (data not shown).Thus, the peak at 930/cm, attributed to the P-OH bending mode, is used to track changes in the amount Table 3.The percentage difference of the area of the 3000-2600/cm region between ASD+ and TD children using raw spectra (no pre-processing applied), area normalized spectra, and min-max normalized spectra (normalized to 1605/cm-urea band).Area values are higher in the control group spectra than in the spectra of ASD + children.  of phosphate excreted in the urine 31 .The calculated area of this peak is compared between the two groups in Fig. 6b.Although the diversity of peak areas is similar in both groups, the mean value is observed to be less in ASD+ children compared to that of the control group.The difference in the means is significant at the 0.05 level.Therefore, it can be concluded that low urinary phosphate levels have contributed to the lesser absorbance in the 3000-2600/cm region of the spectra of ASD+ children.

ASD severity level
It is well known that phosphate (Pi) has essential roles in the intracellular compartment and the skeleton.Many factors may affect its absorption and reabsorption in the kidneys, such as the amount ingested with the diet, parathyroid hormone, other hormones (insulin, growth hormone/insulin-like growth factor 1, calcitonin, glucocorticoids, atrial natriuretic peptide, phosphatonin, prostaglandins, stanniocalcin, etc.), and non-hormonal mechanisms (fasting, plasma calcium, acid-base status, volume status, and medication) 46 .The urinary phosphate excretion is also linked to vitamin D levels, which depend on calcium and phosphate for bone formation.Some ASD+ children were also reported to have vitamin D deficiency, and the amount of serum 25 (OH) D 3 was significantly low 47,48 .In another study, urinary hippurate and phosphate levels were shown to be low in ASD+ children compared to their controls (n = 26 ASD+, 24 TD) 20 .Another previous study also showed lower amounts of phosphoric acid, another phosphorus source in urine, in ASD+ children compared to controls (n = 156 ASD+, 64 TD) 13 .The results of this study are in line with the previous findings regarding the phosphate.

Uric acid
FTIR has several signatures for uric acid; however, many are overlapped by stronger urea absorptions.A small peak at 1345/cm or at 780/cm can be used to monitor uric acid changes.The 780/cm peak is more isolated from the contribution of other urine components, and thus is preferred over the 1345/cm peak, although the analysis of both gives similar results (data not shown).The population variances are significantly different when the 780-peak area is calculated and compared (Fig. 6a).The control group displays close results, but the ASD+ group shows a diversity.The difference between the mean of the peak areas is also significant at the 0.05 level.The 780-peak is less intense in the ASD+ group, meaning lower uric acid levels in the same group than in the controls.As a result, low uric acid levels could have contributed to the absorbance difference measured at the 3000-2600/cm region.
In ASD, the level of antioxidants excreted in urine was shown to be significantly lower than TD children 49 .The same study also showed that the amount of urinary uric acid excreted by children with low-functioning autism (1.23 g/day), medium-functioning autism (0.99 g/day), and high-functioning autism (0.69 g/day) was higher than that by TD children (0.62 g/day).However, two other studies showed that the urinary uric acid level in ASD + children is low compared to the TD children 4,22 .In the normal state, the sulphate excretion rates are lower in children than in adults, and very low in infants, implying the conservation of sulphate during growth 50,51.Sulphate has an inverse effect in the high-frequency region of the spectrum; namely, as the urinary sulphate level increases, the area under the curve decreases in the high-frequency region.Urinary sulphate levels can be followed by the intensity changes at 992 and 1090/cm that are tentatively assigned to R-SO − 2 and S=O, respectively.The relatively more isolated peak at 992/cm is preferred for comparing two groups in terms of the sulphate levels.Figure 6c shows the 992/cm peak area changes using One-way ANOVA for ASD+ and TD.The results show that the population variances and the difference between the means are not significantly different.Therefore, urinary sulphate is not considered as one of the contributors to the observed difference in the high-frequency region.Contrary to the findings in this cohort, high or low sulphate levels have been associated with ASD in previous studies.
Renal sulphate wasting due to mutations in the sulphate reabsorbing carrier protein is found in some cases of autism 52,53 .The low urinary levels of 4-cresol sulphate, hippurate, and creatine in ASD+ children were documented previously 28 .On the other hand, in other studies, the urinary indoxyl sulphate levels were significantly higher in ASD+ children compared to controls 6,17,27 .Another research showed that ASD+ children (average age 7.6) had high urinary sulphate, thiosulphate, and sulphite 54 .A study conducted in the U.S. documented lower sulphate levels in water samples taken from regions where autism is diagnosed more frequently 55 .

Ammonium
The urinary ammonium ion (NH 4 + ) indicates eliminated free hydrogen ions (protons) to balance kidney acidosis.Most of its overall concentration comes from tubular cells, especially the proximal tubular cells-and not due to glomerular filtration rate-as opposed to sulphate.Various factors may affect the activity of the NH 3 -producing enzyme glutaminase, and these factors can be summarized under two categories: systemic and intra-renal.The systemic factors are low levels of plasma bicarbonate, low plasma pH, and the action of adrenal steroids.The intra-renal factors are intracellular potassium deficiency and low intracellular pH.Under the influence of one or more of these factors, glutamine is cleaved by glutaminase to glutamate and NH 3 in tubular cells.NH 3 reacts with the H + ions to form NH 4 in the intra-tubular fluid, thereby neutralizing the H + ions 56 .
The N-H vibrational modes absorb strongly in the high-frequency region, whereas only small changes can be observed in the low-frequency region, particularly at 1660, 1430, and 930/cm.Since the first two of these peaks overlap with strong urea absorptions and the latter is dominated by the phosphate group absorptions, another peak at 3200/cm in the high-frequency region had to be selected.However, urea also has a contribution to the same position.Thus, changes in urea absorption, measured at 3335/cm, are used to eliminate its contribution to the 3200/cm peak.The ratio of the calculated peak areas, i.e., 3200/3335/cm, shows that the ASD+ children have lower values compared to the controls (Fig. 6e).The difference between the means and the population variances are significantly different at the 0.05 level.Low ammonium levels due to the neutral-to-alkaline pH of urine can contribute to the observed smaller area in the 3000-2600/cm region.The average urine pH of our ASD+ children (6.3) is slightly more than that of the TD children (5.6).Although this slight difference in urinary pH could be the cause of the presented spectral findings, to the best of our knowledge, there is no other study showing the same urinary pH trend to back up the current hypothesis.The urea/ornithine cycle is a sequential multiple biochemical reaction that ends with the production of urea (NH 2 ) 2 CO from ammonia (NH 3 ) in the mitochondria of liver cells.The literature shows that there is a relation between urea cycle disturbances and mental retardation 57 .According to a previous study, accumulated ammonia produced by H. pylori causes the observation of elevated serum ammonia among ASD+ children 58 .Several studies have shown that the serum ammonium level is high in ASD+ children with respect to controls [59][60][61] .

Urea
An increase in urea can be directly followed in the spectrum by the growing intensity in the 3420, 3335, and 1609/cm peaks.However, urea has the opposite effect in the 3000-2600/cm region; as the amount increases, the region intensity decreases.Urea (nitrogen compounds) is the end-product of amino acid metabolism; its level in the urine is directly proportional to protein intake.A protein-rich diet tends to decrease the absorbance in the 3000-2600/cm region; however, no such dietary habit was reported in the questionnaires by the families of ASD+ children.In fact, when 3335-peak areas are compared, the difference between the means is not significant and the population variances are similar (Fig. 6d).Therefore, it is suggested that urea does not contribute to the lower intensities observed in the 3000-2600/cm region in the ASD+ group.

Other possible urine components
Previous studies with autistic children documented low levels of phenylactic acid, aconitic acid, 3-oxoglutaric acid, vitamin B6, and kynurenine 6,13 .These components also have large-base and low-amplitude absorbances in the 2600-3000/cm region of the spectrum 62 .Each component listed in the previous sections has typical concentrations ranging from µmol to mmol in the average urine 26 .Therefore, the spectroscopic results provided in Figs. 1 and 2 cannot be the results of one component alone; rather, they are the consequence of a probable common habit or metabolism among autistic children, eventually altering several components within the urine.

Conclusion
In this study, 26 children aged 3-5 years, newly diagnosed with ASD+, using no medications are studied, and the usability of ATR-FTIR is investigated.The high-frequency region of the urine FTIR spectrum is rarely used in literature due to its complex nature of overlapping bands.This study shows that it can be used as a single window that shows the net effect of multiple urine constituents at once that are important for ASD physiology.It is shown that the infrared spectra of ASD+ children display lower absorbance in the 3000-2600/ cm region of the spectrum compared to TD children.PCA is used to identify the peaks that contributed to the dissociation of two groups.It clusters the healthy children group with 87.3% variance along the PC1 axis.The average spectral area of the high-frequency region is 25% less in ASD+ severity level 3 children, 29% less in ASD+ level 2 children, and 16% less in ASD+ level 1 children compared to that of TD children.Level 1 ASD+ children are more challenging to detect with this method since some of the diagnosed children have the same urine spectrum as their controls.Urinary components as possible sources of this outcome are discussed within the manuscript.In previous studies, most of these components have been detected and linked with autism.Spectroscopic results suggest that lower-than-normal levels of uric acid, ammonium, and phosphate groups may have contributed to the observed results in this study.No evidence is found for the contribution of urea and sulphate groups in this study.Also, the average urinary pH is slightly shifted towards neutral in ASD+ children.The nature of the relation (result vs. cause) between these components and the physiology of autism remains an unanswered question.
Studying urine composition component-by-component is not practical, particularly from a clinical point of view.Instead, infrared spectroscopy offers a general picture of urine composition.It is suggested that a visual comparison of the spectrum with that of a healthy child is enough to know that there may be a problem without a calculation.The physician then at least entertains the possibility of autism during standard followups and can take steps, ask questions, or suggest another evaluation.For this disorder, the etiology of which we still do not know and for which there is no biomarker yet, every trait that will contribute to the referral process to child psychiatry is valuable.Even if FTIR does not make a differential diagnosis, it is enough to raise suspicions.
Based on the experience gained in this study, a reference control spectrum to be used in such comparisons should be constructed (by averaging) using the spectra of children of the same age and gender.Thus, it is proposed that the average spectrum constructed from a pool of 3-year-old boys can be used as the negative of another 3-year-old boy.For the evaluation of a 5-year-old girl, a pool of 5-year-old girls can be used.
In order to prove the usefulness of this spectroscopic method for preventive medicine, larger pools of ASD+ and healthy children should be studied.The lack of blood biochemistry results prevented further clinical interpretation of the urine results obtained in this study.For instance, blood gas analysis could have been performed for a more accurate evaluation of urine pH results.The role of ethnicity should also be studied to discuss the possibility of a universal spectrum pool.It is also unknown whether the observed spectral difference is ASDspecific.Other neurodevelopmental disorders like intellectual disability, global developmental delay, attentiondeficit/hyperactivity disorder, motor and tic disorders, specific learning disorders, and communication disorders should be studied in this regard.

Figure 1 .
Figure 1.Recorded infrared spectra of all ASD+ (pale red) and TD (gray) children.The average spectrum of ASD+ is shown in red, and that of TD is in black (thick lines).The spectra shown are not pre-processed.

Figure 2 .
Figure 2. The difference spectra obtained by subtracting the spectrum of each participant from the calculated average spectrum of the control group.The subtraction is formulated as [S participant -S average control ], where S stands for the spectrum.The results of the control group are compiled in the upper panel (gray) and those of the ASD+ group in the lower panel (red).

Figure 3 .
Figure 3. PCA scores (a) and PCA loading plots (b) of ASD (+) and TD children group for the 3800-600/cm region.

Figure 4 .
Figure 4.The mean values of the peak area of the 3000-2600/cm region for ASD+ and TD children (a).A similar comparison is also shown for each ASD level with respect to their controls in (b), (c), and (d) for ASD level 1, level 2, and level 3, respectively.See "Materials and Methods" for statistical details.

Figure 5 .
Figure 5. Spectral contribution of increased ammonium, phosphate, uric acid, sulphate, and urea in the high wavenumber region.(The spectra are obtained using artificial urine solutions, as explained in the Materials and Methods section).

Figure 6 .
Figure 6.Comparison of selected peak areas between the TD group and the ASD+ group of children.Each peak represents a separate urine component, as discussed in the relevant sections.Statistical analysis was done by one-way ANOVA followed by Fisher's post hoc test.Integration limits and baseline choices are given in Table1. https://doi.org/10.1038/s41598-023-46507-z

Table 1 .
. To List of integrated peaks, integration limits, and baseline points.

Table 2 .
40mographic and clinical data of the ASD + and the control group participants.Significant values in bold.Student T test, *Chi Square test, **Mann-Whithey U test, ***One way ANOVA.TL Turkish Lira, $ US Dollar.The hunger threshold was $340, and the poverty threshold was $1100 during the time data was collected40.

Table 1 .
Many factors influence the sulphate levels in urine, such as age, gender, pregnancy, medication, disorders, etc.