Transcript levels of 4 genes in umbilical cord blood are predictive of later autism development: a longitudinal follow-up study =============================================================================================================================== * Qingzheng Jia * Huilin Li * Min Wang * Chongxia Wei * Lichao Xu * Lin Ye * Chunjun Wang * Shaofeng Ke * Ling Li * Paul Yao ## Abstract **Background:** Over recent decades, autism spectrum disorder (ASD) has been of increasing epidemiological importance, given the substantial increase in its prevalence; at present, clinical diagnosis is possible only after 2 years of age. In this study, we sought to develop a potential predictive model for ASD screening. **Methods:** We conducted a longitudinal follow-up study of newborns over 3 years. We measured transcript levels of 4 genes (superoxide dismutase-2 [*SOD2*], retinoic acid–related orphan receptor-α [*RORA*], G protein–coupled estrogen receptor-1 [*GPER*], progesterone receptor [*PGR*]), 2 oxidative stress markers and epigenetic marks at the *RORA* promoter in case–control umbilical cord blood mononuclear cell (UCBMC) samples. **Results:** We followed 2623 newborns; we identified 41 children with ASD, 63 with delayed development and 2519 typically developing children. We matched the 41 children with ASD to 41 typically developing children for UCBMC measurements. Our results showed that children with ASD had significantly higher levels of H3K9me3 histone modifications at the *RORA* promoter and oxidative stress in UCBMC than typically developing children; children with delayed development showed no significant differences. Children with ASD had significantly lower expression of *SOD2*, *RORA* and *GPER*, but higher *PGR* expression than typically developing children. We established a model based on these 4 candidate genes, and achieved an area under the curve of 87.0% (standard deviation 3.9%) with a sensitivity of 1.000 and specificity of 0.854 to predict ASD in UCBMC. **Limitations:** Although the gene combinations produced a good pass/fail cut-off value for ASD evaluation, relatively few children in our study sample had ASD. **Conclusion:** The altered gene expression in UCBMC can predict later autism development, possibly providing a predictive model for ASD screening immediately after birth. ## Introduction Autism spectrum disorder (ASD) has been of increasing epidemiological importance, given the substantial increase in its prevalence over the past 2 decades,1 with heavy economic and social impacts, in addition to a variety of comorbidities.2,3 Unfortunately, the mechanism and causes of ASD remain largely unknown,4–7 and there is no effective cure for ASD, although some interventions can slightly improve the symptoms.8,9 Clinical diagnosis for ASD is usually possible only after 2 years of age, which can be too late for efficient therapeutic treatment; thus, development of a sensitive biomarker or predictive model for early ASD screening is urgently needed.10 We have previously reported that superoxide dismutase-2 (*SOD2*) is suppressed in maternal diabetes–mediated offspring4 and *SOD2* transcripts are reduced in samples of peripheral blood mononuclear cells (PBMC) among children with ASD compared with typically developing children.11 Furthermore, retinoic acid–related orphan receptor-α (*RORA*) has been reported to be an important factor for ASD development,12 and its expression is reduced in samples of PBMC from participants with ASD.13,14 Taken altogether, the evidence suggests that both *SOD2* and *RORA* can be potential biomarkers for ASD diagnosis. We have previously found that estrogen receptor-β (ERβ) expression in the brain is reduced, contributing to ASD development, while the expression in PBMC is too low for detection of significant changes among participants with ASD.5,15 On the other hand, G protein–coupled estrogen receptor-1 (*GPER*), is either directly or indirectly involved in estrogen function through the nongenomic signalling pathway and plays a role in neurologic disorders.16 Levels of GPER in the serum are reduced in participants with ASD, compared with typically developing participants.17 In addition, prenatal hormone exposure is associated with ASD development6 through suppression of ERβ5,18 or oxytocin;15 preliminary data have shown that, in the presence of progestin exposure, the expression of progesterone receptor (*PGR*) increases significantly in human neural progenitor cells and remains high after removal of progestin. We assume that both *GPER* and *PGR* are involved in ASD development and may be potential biomarkers for ASD screening. We sought to develop a sensitive predictive model for ASD screening based on gene expression in mononuclear cells from children with ASD. ## Methods We conducted a longitudinal follow-up study19 involving newborns and their parents, registered in the Hainan Women and Children’s Medical Center, China, from July 2018 to July 2019. Further methodological details are provided in Appendix 1, available at [https://www.jpn.ca/lookup/doi/10.1503/jpn.230046/tab-related-content](https://www.jpn.ca/lookup/doi/10.1503/jpn.230046/tab-related-content). All newborns were required to have regular follow-up behaviour and development assessment every 6 months over the 3 years following enrolment. A participation form was completed to collect the information regarding the newborn child, namely birth date, sex, preterm birth (< 37 wk gestation), Apgar score at 5 minutes to evaluate birth asphyxia (Apgar < 6), family history of mental disorders, father’s age and personality (extroverted, ambivert or introverted), mother’s age and personality, major stressful life events, paternal smoking history and maternal abortion history. We identified typically developing children, children with delayed development (DD) and children with ASD between the ages of 24 and 36 months. We identified children with ASD using the Autism Behaviour Checklist (ABC), in which ABC scores greater than 53 points indicate high probability of ASD and scores greater than 67 points indicate obvious ASD symptoms; children were also assessed by a licensed clinical team.6,20–24 Typically developing children had to have no history of either ASD or DD, with scores greater than 70 on the Mullen Scales of Early Learning and Vineland Adaptive Behavior Scales; children with DD had to have scores less than 70 on these measures.25 We collected 5–20 mL of umbilical cord blood from newborns using EDTA as an anticoagulant. We isolated umbilical cord blood mononuclear cells (UCBMC) and purified RNA to measure transcript levels of 4 genes (*SOD2*, *RORA*, *GPER*, *PGR*) using the primers indicated in Appendix 1, Table S1. We also measured oxidative stress (ratio of glutathione to glutathione disulfide, 8-oxo-2′-deoxyguanosine [8-oxo-dG] staining) and epigenetic marks (histone methylation and DNA methylation) on the *RORA* promoter. ### Population-based case–control study We matched each identified child with ASD with 1 control, randomly selected from the sample of typically developing children. We assessed the intelligence quotient (IQ) for both children with ASD and controls using the Wechsler Intelligence Scale for Children IV, and took care to ensure that the 2 groups had no significant differences in age, sex, full-scale IQ, performance IQ and verbal IQ. We collected 3–5 mL of peripheral blood for transcript analysis from participants at the age of 30 months. We employed transcript levels of the 4 candidate genes (*SOD2*, *RORA*, *GPER* and *PGR*) from both umbilical cord blood and peripheral blood to define the pass/fail cut-off value for ASD screening. ### Pass/fail cut-off value for ASD screening We calculated the pass/fail cut-off value for ASD screening in 5 steps. First, we calculated relative transcript levels for each participant using a quantitative polymerase chain reaction (PCR). Next, we calculated a gene index for *SOD2*, *RORA*, *GPER* and *PGR* (e.g., *SOD2* index = 1 − [mean of *SOD2* transcript level for ASD group/mean of *SOD2* transcript level for typically developing group]). We determined the gene weight for the *SOD2*, *RORA*, *GPER* and *PGR* (e.g., *SOD2* weight = *SOD2* index/[*SOD2* index + *RORA* index + *PGR* index + *GPER* index]). We then determined the ASD index as the total of (*SOD2* weight × *SOD2* transcript level) + (*RORA* weight × *RORA* transcript level) + (*GPER* weight × GPER transcript level) + (*PGR* weight × *PGR* transcript level). Finally, we employed the ASD index to generate the ROC curve, and calculate the pass/fail cut-off value for ASD screening. ### Biomedical analysis Mononuclear cells were isolated from either 3–5 mL of peripheral blood or 5–20 mL of umbilical cord blood using Lymphoprep reagents (#07861, from STEMCELL Technologies). We purified the total RNA from mononuclear cells using either the RNeasy Plus Mini Kit (peripheral blood) or RNeasy Micro Kit (umbilical cord blood) from Qiagen. We evaluated transcript levels by reverse transcription reaction combined with quantitative PCR using primers described in Appendix 1, Table S1. We calculated relative results via the ΔΔCT method using β-actin as a housekeeping gene. We evaluated the redox balance by 8-oxo-dG formation through immunostaining and measurement of the ratio of glutathione to glutathione disulfide. We evaluated the degree of epigenetic modification on the RORA promoter14 by histone methylation through chromatin immunoprecipitation using primers provided in Appendix 1, Table S1. We analyzed DNA methylation on the *RORA* promoter through methylation-specific PCR analysis, combined with bisulfite modification in the region of −1200 to −700 on the *RORA* promoter, as previously described.14 For the analysis of redox balance and epigenetic modifications on the *RORA* promoter, we randomly selected individuals from each group using true random number generators.26 We describe results using means and standard deviations (SDs). The ratio difference was determined using the χ2 test; we used SPSS 22 software to generate the ROC curve and calculate the pass/fail cut-off value.25 ### Ethics approval The protocol was reviewed and approved by Hainan Women and Children’s Medical Center, and all the participants provided written consent forms. ## Results After a 3-year longitudinal follow-up study period, we identified 41 children with ASD, 63 children with DD and 2519 typically developing children (Figure 1). Among the 41 children with ASD, 28 were further identified as having mild ASD and 13 had severe ASD based on ABC scores and opinions from clinical professionals. General demographic and clinical information on participants is summarized in Table 1. Age, paternal age and personality, maternal personality, major stressful life event, paternal smoking and maternal abortion history showed no significant differences by group; sex, family history of mental disorders, preterm birth, birth asphyxia and maternal age showed significant differences among the groups (Table 1).6 In addition, very few women in China are current smokers. Among children with ASD, none of the mothers were smokers; we did not consider maternal smoking as a risk factor for ASD development. Detailed demographic information for the matched case–control subgroup is shown in Appendix 1, Table S2. ![Figure 1](http://jpn.ca/https://www.jpn.ca/content/jpn/48/5/E334/F1.medium.gif) [Figure 1](http://jpn.ca/content/48/5/E334/F1) Figure 1 Research design and participant flow chart for longitudinal follow-up study of autism spectrum disorder (ASD). DD = delayed development; *GPER* = G protein–coupled estrogen receptor 1; MNC = mononuclear cells; PBMC = peripheral blood mononuclear cells; *PGR* = progesterone receptor; ROC = receiver operating characteristic; *RORA* = retinoic acid–related orphan receptor-α; *SOD2* = superoxide dismutase-2; TD = typically developing; UCBMC = umbilical cord blood mononuclear cells. View this table: [Table 1](http://jpn.ca/content/48/5/E334/T1) Table 1 Participant demographic and clinical information ### Epigenetic changes in umbilical cord blood mononuclear cells We investigated potential epigenetic modifications in UCBMC in 9 randomly selected samples from each group. For the ASD group, we randomly selected 6 samples from the mild ASD group and 3 samples from the severe ASD group based on the ratio of these 2 types within the total sample of children with ASD. We found that children in the ASD group had significantly potentiated H3K9me3 histone modifications on the *RORA* promoter at a level 1.87 times higher than the typically developing group, while there was no effect on the DD group. We found no significant difference in H3K27me2, H3K27me3 and H3K9me2, or in DNA methylation, on the *RORA* promoter across those groups (Figure 2A and 2B). The potentiated H3K9me3 modification may contribute to reduced *RORA* expression, while this does not exclude the other possibilities, such as histone H4 modification and histone H3 acetylation. In addition, our data showed that there was no significant difference in epigenetic changes on the *RORA* promoter between children with mild ASD (*n* = 6) and those with severe ASD (*n* = 3). ![Figure 2](http://jpn.ca/https://www.jpn.ca/content/jpn/48/5/E334/F2.medium.gif) [Figure 2](http://jpn.ca/content/48/5/E334/F2) Figure 2 Analysis of epigenetic changes in (A) histone methylation and (B) DNA methylation on the retinoic acid–related orphan receptor-α (*RORA*) promoter region from umbilical cord blood mononuclear cells, isolated from 9 randomly selected participants from each group (typically developing [TD] children, children with delayed development [DD], children with autism spectrum disorder [ASD]). **p* < 0.05 v. TD group. ### Redox balance in umbilical cord blood mononuclear cells We found that the ASD group had a significantly lower ratio of glutathione to glutathione disulfide (62 % of the typically developing group), while the DD group showed no significant differences (Figure 3A). Immunostaining of 8-oxo-dG, and the result showed that 8-oxo-dG formation was 1.96 times higher among children with ASD than the typically developing group; again, the DD group showed no significant differences (Figure 3B and 3C). ![Figure 3](http://jpn.ca/https://www.jpn.ca/content/jpn/48/5/E334/F3.medium.gif) [Figure 3](http://jpn.ca/content/48/5/E334/F3) Figure 3 Analysis of redox balance in umbilical cord blood mononuclear cells, isolated from 9 randomly selected participants from each group (typically developing [TD] children, children with delayed development [DD], children with autism spectrum disorder [ASD]). (A) Intracellular ratio of glutathione (GSH) to glutathione disulfide (GSSG). (B) Quantification of 8-oxo-2′-deoxyguanosine (8-oxo-dG). (C) Representative staining photos for quantification of 8-oxo-dG. DAPI = 4′,6-diamidino-2-phenylindole. **p* < 0.05 v. TD. ### Gene expression in umbilical cord blood mononuclear cells We found that the ASD group showed significantly lower transcript levels of *SOD2* (Figure 4A), *RORA* (Figure 4B) and *GPER* (Figure 4C), at 56 %, 54 % and 74 %, respectively, than the typically developing group. In addition, transcript levels of *PGR* were 233 % higher in the ASD group than the typically developing group (Figure 4D). The DD group showed no significant differences. ![Figure 4](http://jpn.ca/https://www.jpn.ca/content/jpn/48/5/E334/F4.medium.gif) [Figure 4](http://jpn.ca/content/48/5/E334/F4) Figure 4 Analysis of gene expression of (A) superoxide dismutase-2 (*SOD2*), (B) retinoic acid–related orphan receptor-α (*RORA*), (C) G protein–coupled estrogen receptor-1 (*GPER*) and (D) progesterone receptor (*PGR*) in umbilical cord blood mononuclear cells isolated from 2519 typically developing (TD) children, 63 children with delayed development (DD) and 41 children with autism spectrum disorder (ASD). NS = not significant. **p* < 0.05 v. TD group. ### Predictive model for screening of autism spectrum disorder from umbilical cord blood mononuclear cells We analyzed transcript levels of *SOD2*, *RORA*, *GPER* and *PGR* in UCBMC, and calculated the ASD Index for each participant with ASD (*n* = 41) and matched typically developing participants (*n* = 41); we established the related ROC curve for ASD. We achieved an area under the curve (AUC) of 87.0% (SD 3.9%) for identifying ASD, with a 95% confidence interval (CI) of 79.4%–94.7%; the pass/fail cut-off value was 0.3156 when the sensitivity was set to 1.000 and specificity was set to 0.854 (Figure 5 and Appendix 1, Figure S1), which shows high sensitivity and specificity for ASD screening using this method (Table 2). In addition, we calculated the ASD index for individual genes or gene combinations, and our results showed that the related AUCs were less than 85.0 % and that the combination of *SOD2*, *RORA*, *GPER* and *PGR* had the best results (Appendix 1, Table S3). ![Figure 5](http://jpn.ca/https://www.jpn.ca/content/jpn/48/5/E334/F5.medium.gif) [Figure 5](http://jpn.ca/content/48/5/E334/F5) Figure 5 Receiver operating characteristic (ROC) curve of predictive model for screening of autism spectrum disorder based on transcript levels of superoxide dismutase-2, retinoic acid–related orphan receptor-α, G protein–coupled estrogen receptor-1 and progesterone receptor among typically developing children (*n* = 41) and children with autism spectrum disorder (ASD; *n* = 41), showing sensitivity (blue line) and specificity (green line). The ASD index for each participant was employed to calculate the ROC curve. View this table: [Table 2](http://jpn.ca/content/48/5/E334/T2) Table 2 Calculation of pass/fail cut-off value for screening of autism spectrum disorder ### Association between autism spectrum disorder index and symptoms We evaluated the potential association of the ASD index with ASD symptoms in the case–control study. Based on ABC scores and opinions from clinical professionals for each participant, the ASD group was further identified to include 28 children with mild ASD and 13 with severe ASD. The ASD index for each control and ASD participant was calculated based on the transcripts of genes for *SOD2*, *RORA*, *GPER* and *PGR* (Appendix 1, Figure S2). The results showed that the severe ASD group had a significantly lower ASD index than the control and mild ASD groups; the mild ASD group showed nearly statistically significant differences (*p* = 0.056 v. the control group), which may be related to the small sample size. In general, our results indicate that ASD index may reflect, or be associated with, ASD symptoms. ### Gene expression in peripheral blood mononuclear cells We found that the ASD group had significantly lower transcript levels of *SOD2* (Figure 6A), *RORA* (Figure 6B) and *GPER* (Figure 6C) at 51 %, 36 % and 38 %, respectively, compared with typically developing participants. In addition, transcript levels of *PGR* were 202 % higher in the ASD group than the typically developing group (Figure 6D). The DD group showed no significant differences in gene expression. ![Figure 6](http://jpn.ca/https://www.jpn.ca/content/jpn/48/5/E334/F6.medium.gif) [Figure 6](http://jpn.ca/content/48/5/E334/F6) Figure 6 Analysis of gene expression of (A) superoxide dismutase-2 (*SOD2*), (B) retinoic acid–related orphan receptor-α (*RORA*), (C) G protein–coupled estrogen receptor-1 (*GPER*) and (D) progesterone receptor (*PGR*) in peripheral blood mononuclear cells isolated from 41 typically developing (TD) children and 41 children with autism spectrum disorder. **p* < 0.05, v. TD group. ### Predictive model for screening of autism spectrum disorder from peripheral blood mononuclear cells We analyzed transcript levels of *SOD2*, *RORA*, *GPER* and *PGR* in UCBMC, and calculated the ASD index for each participant with ASD (*n* = 41) and matched typically developing participants (*n* = 41); we established the related ROC curve for ASD. We achieved an AUC of 84.1% (SD 2.3%), with a 95% CI of 79.6%–88.7%; the pass/fail cut-off value was defined to 0.5505 when the sensitivity was set to 0.976 and the specificity was set to 0.854 (Figure 7 and Appendix 1, Figure S3), which shows less sensitivity and specificity for ASD screening than when using UCBMC (Table 2). ![Figure 7](http://jpn.ca/https://www.jpn.ca/content/jpn/48/5/E334/F7.medium.gif) [Figure 7](http://jpn.ca/content/48/5/E334/F7) Figure 7 Receiver operating characteristic (ROC) curve of autism spectrum disorder (ASD) index based on transcript levels of superoxide dismutase-2, retinoic acid–related orphan receptor-α, G protein–coupled estrogen receptor-1 and progesterone receptor among typically developing children (*n* = 41) and children with ASD (*n* = 41), showing sensitivity (blue line) and specificity (green line). ## Discussion In this longitudinal follow-up study, we showed that participants with ASD had higher levels of histone methylation on the *RORA* promoter and oxidative stress, as well as altered gene expression in UCBMC, compared with typically developing children; children with DD showed no significant differences. The population-based case–control study showed that, combined, transcript levels of *SOD2*, *RORA*, *GPER* and *PGR* in both UCBMC and PBMC can accurately predict ASD development. It has been reported that prenatal factors, including maternal diabetes,4,14,27,28 progestin levels5,15 and androgen levels,29 contribute to autism-like behaviours in rodent offspring through oxidative stress–mediated epigenetic changes, triggering consistent altered expression of genes such as *SOD2*4 and *RORA*14 in both brain tissues and hematopoietic stem cells.11 Interestingly, mononuclear cells from either cord blood or peripheral blood can inherit the same kinds of epigenetic changes during embryo differentiation11 and show similar epigenetic modifications and altered gene expression during later ASD development in childhood.14,30 This is the rationale for analyzing altered gene expression in either UCBMC or PBMC as a predictive model for ASD screening. We found that children with ASD have increased histone methylation on the *RORA* promoter, as well as increased oxidative stress, as indicated by the potentiated ratio of glutathione to glutathione disulfide and 8-oxo-dG formation, compared with typically developing children. Those with DD showed no significant differences, indicating that the potential driving force within the DD group is not prenatal factor exposure-mediated epigenetic modifications and altered redox balance, and instead may be owing to other factors, such as family history of mental diseases and birth asphyxia.25 It has been reported that prenatal factor exposure results in epigenetic modifications with consistent suppressed gene expression in neural progenitor cells, including ERβ,5 *SOD2*4 and *RORA*,14 contributing to ASD development. Interestingly, prenatal factor–mediated epigenetic changes can also occur in hematopoietic stem cells during embryonic development, and PBMC can inherit the same epigenetic changes during stem cell differentiation, subsequently resulting in a similar gene expression pattern as neurons among offspring;11 this is the rationale for the development of an ASD predictive model based on altered gene expression in both UCBMC and PBMC.30 Our results showed that children with ASD have lower transcript levels of *SOD2*, *RORA* and *GPER*, but have higher *PGR* transcript levels in UCBMC and PBMC than typically developing children. Moreover, the pattern of gene expression is similar to the pattern in neurons from people with ASD, indicating that altered gene expression in either UCBMC or PBMC can be a potential predictive model for ASD screening. In addition, although *GPER* and *PGR* are sex hormone targets, our preliminary results showed that the mRNA levels of *GPER* and *PGR* had no significant differences by sex in both UCBMC and PBMC. Autism spectrum disorder can be triggered by many prenatal factors,31–33 and each factor may mediate a different pattern of gene expression in mononuclear cells. Because of this, expression of a single gene cannot be an accurate predictor of ASD development. Our results showed that either individual genes or a combination of 2 or 3 genes can achieve only an AUC of less than 85 % for prediction of ASD development. In contrast, the combined expression of *SOD2*, *RORA*, *GPER* and *PGR* can identify ASD with an AUC of 87%, with a sensitivity of 1.000 and specificity of 0.854 in UCBMC.34 Altered gene expression in PBMC achieved less accuracy for predicting ASD (AUC 84%), with a sensitivity of 0.976 and specificity of 0.854. Gene expression in PBMC may be altered or interfered with by many other factors, such as redox situations, stressful events or food consumption and metabolism in the human body.35,36 ### Limitations The number of participants with ASD is relatively low for appropriate use of developing a predictive model. Currently, we are conducting an independent study for replication in a different prospective cohort, which is essential for demonstrating reproducibility of a predictive model. Our groups had several significant differences in baseline characteristics, such as maternal age, sex, and preterm birth; further investigation for the potential role of those variables is needed. Cord blood has many varying cell types with different expression and epigenetic patterns, and we determined only the epigenetic changes in a single histone mark on a single gene, rather than global changes. ## Conclusion In this longitudinal follow-up study of 2623 newborns, children with ASD showed significantly higher levels of epigenetic modification on the *RORA* promoter and oxidative stress, as well as altered gene expression of *SOD2*, *RORA*, *GPER* and *PGR* in UCBMC, compared with typically developing children. A predictive model based on altered gene expression of these 4 gene candidates in UCBMC can identify 87% of children with ASD with a sensitivity of 1.000 and specificity of 0.854. We conclude that altered gene expression in umbilical cord blood can predict later development of autism. ## Data availability: All data are available in either the manuscript or supplementary information. ## Footnotes * * These authors contributed equally to this work. * **Competing interests:** None declared. * **Contributors:** Ling Li and Paul Yao conceived and designed the work. Ling Li and Paul Yao acquired the data. All of the authors contributed to data analysis and interpretation. Paul Yao drafted the manuscript. All of the authors revised it critically for important intellectual content, gave final approval of the version to be published and agreed to be accountable for all aspects of the work. * **Funding:** This study was financially supported by The National Natural Science Foundation of China (no. 82060260), The National Natural Science Foundation of Hainan Province (no. 821QN1000), The Excellent Talent Team of Hainan Province (no. QRCBT202121) and Hainan Province Clinical Medical Center (no. QWYH202175). * Received March 24, 2023. * Revision received May 8, 2023. * Revision received May 31, 2023. * Accepted June 11, 2023. This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY-NC-ND 4.0) licence, which permits use, distribution and reproduction in any medium, provided that the original publication is properly cited, the use is noncommercial (i.e., research or educational use), and no modifications or adaptations are made. See: [https://creativecommons.org/licenses/by-nc-nd/4.0/](https://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. Maenner MJ, Shaw KA, Baio J, EdS, Washington A, Patrick 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. [CrossRef](http://jpn.ca/lookup/external-ref?access_num=10.15585/mmwr.ss6904a1&link_type=DOI) [PubMed](http://jpn.ca/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fjpn%2F48%2F5%2FE334.atom) 2. Dunalska A, Rzeszutek M, Debowska Z, et al. Comorbidity of bipolar disorder and autism spectrum disorder — review paper. Psychiatr Pol 2021;55:1421–31. 3. Mutluer T, Aslan Genc H, Ozcan Morey A, et al. Population-based psychiatric comorbidity in children and adolescents with autism spectrum disorder: a meta-analysis. Front Psychiatry 2022;13: 856208. 4. Wang X, Lu J, Xie W, et al. Maternal diabetes induces autism-like behavior by hyperglycemia-mediated persistent oxidative stress and suppression of superoxide dismutase 2. Proc Natl Acad Sci U S A 2019;116:23743–52. [Abstract/FREE Full Text](http://jpn.ca/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE2LzQ3LzIzNzQzIjtzOjQ6ImF0b20iO3M6MTk6Ii9qcG4vNDgvNS9FMzM0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 5. Zou Y, Lu Q, Zheng D, et al. Prenatal levonorgestrel exposure induces autism-like behavior in offspring through ERβ suppression in the amygdala. Mol Autism 2017;8:46. 6. Li L, Li M, Lu J, et al. Prenatal progestin exposure is associated with autism spectrum disorders. Front Psychiatry 2018;9:611. 7. Lin CW, Septyaningtrias DE, Chao HW, et al. A common epigenetic mechanism across different cellular origins underlies systemic immune dysregulation in an idiopathic autism mouse model. Mol Psychiatry 2022;27:3343–54. 8. Shoaib A, Cepeda MS, Murray G, et al. Autism: comorbidities and treatment patterns in the real world, a retrospective cohort study among children, adolescents and adults newly diagnosed with autism. J Autism Dev Disord 2022;52:4311–20. 9. Villarreal-Martinez L, Gonzalez-Martínez G, Saenz-Flores M, et al. Stem cell therapy in the treatment of patients with autism spectrum disorder: a systematic review and meta-analysis. Stem Cell Rev Rep 2022;18:155–64. 10. Frye RE, Vassall S, Kaur G, et al. Emerging biomarkers in autism spectrum disorder: a systematic review. Ann Transl Med 2019;7:792. 11. Lu J, Xiao M, Guo X, et al. Maternal diabetes induces immune dysfunction in autistic offspring through oxidative stress in hematopoietic stem cells. Front Psychiatry 2020;11:576367. 12. Sayad A, Noroozi R, Omrani MD, et al. Retinoic acid-related orphan receptor alpha (RORA) variants are associated with autism spectrum disorder. Metab Brain Dis 2017;32:1595–601. [CrossRef](http://jpn.ca/lookup/external-ref?access_num=10.1007/s11011-017-0049-6&link_type=DOI) 13. Hu VW, Sarachana T, Sherrard RM, et al. Investigation of sex differences in the expression of RORA and its transcriptional targets in the brain as a potential contributor to the sex bias in autism. Mol Autism 2015;6:7. 14. Yu H, Niu Y, Jia G, et al. Maternal diabetes-mediated RORA suppression in mice contributes to autism-like offspring through inhibition of aromatase. Commun Biol 2022;5:51. 15. Huang S, Zeng J, Sun R, et al. Prenatal progestin exposure-mediated oxytocin suppression contributes to social deficits in mouse offspring. Front Endocrinol (Lausanne) 2022;13:840398. 16. Roque C, Mendes-Oliveira J, Duarte-Chendo C, et al. The role of G protein-coupled estrogen receptor 1 on neurological disorders. Front Neuroendocrinol 2019;55:100786. 17. Altun H, Kurutas EB, Sahin N, et al. Decreased levels of G protein-coupled estrogen receptor in children with autism spectrum disorders. Psychiatry Res 2017;257:67–71. [CrossRef](http://jpn.ca/lookup/external-ref?access_num=10.1016/j.psychres.2017.06.008&link_type=DOI) [PubMed](http://jpn.ca/lookup/external-ref?access_num=28734238&link_type=MED&atom=%2Fjpn%2F48%2F5%2FE334.atom) 18. Xie W, Ge X, Li L, et al. Resveratrol ameliorates prenatal progestin exposure-induced autism-like behavior through ERβ activation. Mol Autism 2018;9:43. 19. Elias R, Lord C. Diagnostic stability in individuals with autism spectrum disorder: insights from a longitudinal follow-up study. J Child Psychol Psychiatry 2022;63:973–83. 20. Mutluer T, Karakoc Demirkaya S, Abali O. Assessment of sleep problems and related risk factors observed in Turkish children with autism spectrum disorders. Autism Res 2016;9:536–42. 21. Rausch A, Zhang W, Haak KV, et al. Altered functional connectivity of the amygdaloid input nuclei in adolescents and young adults with autism spectrum disorder: a resting state fMRI study. Mol Autism 2016;7:13. [CrossRef](http://jpn.ca/lookup/external-ref?access_num=10.1186/s13229-015-0060-x&link_type=DOI) 22. Verhoeff ME, Blanken LME, Kocevska D, et al. The bidirectional association between sleep problems and autism spectrum disorder: a population-based cohort study. Mol Autism 2018;9:8. 23. Lord C, Rutter M, Le Couteur A. Autism diagnostic interview-revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders. J Autism Dev Disord 1994;24:659–85. [CrossRef](http://jpn.ca/lookup/external-ref?access_num=10.1007/BF02172145&link_type=DOI) [PubMed](http://jpn.ca/lookup/external-ref?access_num=7814313&link_type=MED&atom=%2Fjpn%2F48%2F5%2FE334.atom) [Web of Science](http://jpn.ca/lookup/external-ref?access_num=A1994PM02100009&link_type=ISI) 24. Lord C, Rutter M, Goode S, et al. Autism diagnostic observation schedule: a standardized observation of communicative and social behavior. J Autism Dev Disord 1989;19:185–212. [CrossRef](http://jpn.ca/lookup/external-ref?access_num=10.1007/BF02211841&link_type=DOI) [PubMed](http://jpn.ca/lookup/external-ref?access_num=2745388&link_type=MED&atom=%2Fjpn%2F48%2F5%2FE334.atom) [Web of Science](http://jpn.ca/lookup/external-ref?access_num=A1989AC89500001&link_type=ISI) 25. Kim DH, Krakowiak P, Meltzer A, et al. Neonatal chemokine markers predict subsequent diagnosis of autism spectrum disorder and delayed development. Brain Behav Immun 2022;100:121–33. 26. Rojas-Muñoz LF, Sanchez-Solano S, Martinez-Rodriguez MC, et al. On-line evaluation and monitoring of security features of an RO-based PUF/TRNG for IoT devices. Sensors (Basel) 2023; 23:4070. 27. Xiao L, Wang M, Zhang W, et al. Maternal diabetes-mediated RORA suppression contributes to gastrointestinal symptoms in autism-like mouse offspring. BMC Neurosci 2022;23:8. 28. Liu J, Liang Y, Jiang X, et al. Maternal diabetes-induced suppression of oxytocin receptor contributes to social deficits in offspring. Front Neurosci 2021;15:634781. 29. Xiang D, Lu J, Wei C, et al. Berberine ameliorates prenatal dihydrotestosterone exposure-induced autism-like behavior by suppression of androgen receptor. Front Cell Neurosci 2020;14:87. 30. Zeng J, Liang Y, Sun R, et al. Hematopoietic stem cell transplantation ameliorates maternal diabetes-mediated gastrointestinal symptoms and autism-like behavior in mouse offspring. Ann N Y Acad Sci 2022;1512:98–113. 31. Modabbernia A, Velthorst E, Reichenberg A. Environmental risk factors for autism: an evidence-based review of systematic reviews and meta-analyses. Mol Autism 2017;8:13. [CrossRef](http://jpn.ca/lookup/external-ref?access_num=10.1186/s13229-017-0121-4&link_type=DOI) [PubMed](http://jpn.ca/lookup/external-ref?access_num=28331572&link_type=MED&atom=%2Fjpn%2F48%2F5%2FE334.atom) 32. 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. [CrossRef](http://jpn.ca/lookup/external-ref?access_num=10.1038/mp.2011.165&link_type=DOI) [PubMed](http://jpn.ca/lookup/external-ref?access_num=22143005&link_type=MED&atom=%2Fjpn%2F48%2F5%2FE334.atom) [Web of Science](http://jpn.ca/lookup/external-ref?access_num=000302173900006&link_type=ISI) 33. Pham C, Symeonides C, O’Hely M, et al. Early life environmental factors associated with autism spectrum disorder symptoms in children at age 2 years: a birth cohort study. Autism 2022; 26:1864–81. 34. Hewitson L, Mathews JA, Devlin M, et al. Blood biomarker discovery for autism spectrum disorder: a proteomic analysis. PLoS One 2021;16:e0246581. 35. Doenyas C. Potential role of epigenetics and redox signaling in the gut-brain communication and the case of autism spectrum disorder. Cell Mol Neurobiol 2022;42:483–7. 36. Park-Cardoso J, Silva A. Insistence on sameness for food space appropriation: An exploratory study on Brazilians with autism (self-) diagnosis in adulthood. Autism 2023;27:938–51.