- 1Department of Child Health Care, Nanjing Maternity and Child Health Care Hospital, Women’s Hospital of Nanjing Medical University, Nanjing, China
- 2Shanghai Institute of Hematology, State Key Laboratory of Medical Genomics, National Research Center for Translational Medicine at Shanghai, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, China
Autism spectrum disorder (ASD) is a highly genetic heterogeneous neurodevelopmental disorder, which is usually considered a heritable and heterogeneous neurodevelopmental disorder and has caused a great burden to society and families. Emerging roles of ferroptosis have been observed in neurological disorders. This study aimed to construct a diagnostic model based on ferroptosis-related genes (FRGs) to contribute to the early and precise diagnosis of childhood ASD. In the candidate FRGs, we identified 27 differentially expressed genes (DEGs) between ASD patients and typically developing (TD) controls. Four key FRGs were identified using the random forest analysis for further analysis. Utilization of the four gene expression, we constructed a diagnostic model and the AUC value in the training dataset (GSE18123) is 0.7002. We deem that a patient with a score less than 0.9904 is likely to have ASD. Three validation datasets (GSE111176, GSE113834, and GSE28521) were collected and the AUC value is 0.7442, 0.7444, and 0.6474, respectively. A multi-factor regulatory network based on four FRGs indicated that RORA, EAF1, NFYB, miR-4703-3p, and miR-6073 may play a role in the development of ASD. In addition, we found piperaquine may have the potential to be a promising drug for the treatment of ASD. Overall, we constructed a diagnostic model of childhood ASD, which could contribute to the precision diagnosis and timely treatment of childhood ASD.
Introduction
Autism spectrum disorder (ASD) is a set of heterogeneous genetically complex neurodevelopmental disorders, characterized by a deficit in social communication and interaction, and lasting impairments in restricted, stereotyped, and repetitive patterns of behavior or interests (1). The prevalence of ASD continues to increase worldwide since the first epidemiological study, latest study reported an incidence of 1 in 54 children worldwide (2). Males are four times more likely to develop autism than girls (3). Several complex underlying mechanisms were involved in the development of ASD, such as genetic and epigenetic effects, inflammation, oxidative stress, neurotrophic factors, and hypoxic damage (4–6).
Early diagnosis and appropriate intervention for childhood ASD can improve ASD outcomes. However, diagnosis of autism was often delayed until mid-childhood because of complex symptoms, insufficient screening practices, the lack of clinical access, low sensitivity of autism screening instruments, and ignorance of early warning signs (7, 8). In addition, all present medications that have evidence to be beneficial to autism are aimed at related symptoms, rather than directly at autism symptoms (including repetitive behaviors or social communication) (9). With the development of sequencing technologies, abundant transcriptomic analyses improve our understanding of ASD and allow us to develop a novel model contributing to the early precise diagnosis of childhood ASD.
In recent years, ferroptosis is a newly identified iron-dependent form of programmed cell death, primarily caused by the imbalance of oxidation and anti-oxidation in the body (10). Ferroptosis is typified by lipid peroxidation and depends upon the severe lipid peroxidation of intracellular iron accumulation and generation of reactive oxygen species (ROS) (11). A growing body of evidence shows that ferroptosis is strongly implicated in a variety of diseases, including cancers, cardiovascular diseases, kidney damage, and nervous system diseases (12). Ferroptosis is closely linked to the occurrence and development of neurodegenerative diseases, strokes, and brain tumors. It is also involved in the development, maturation, and aging of the nervous system (13, 14). Many studies have found a link between ASD and elevated oxidative stress and ASD is contributed to by oxidative stress in several ways, including protein post-translational changes, abnormal metabolism (e.g., lipid peroxidation), and toxic buildup (e.g., ROS) (15). However, the definite role of ferroptosis in ASD is still unclear. In this study, we aimed to explore the potential association between ferroptosis-related genes (FRGs) and ASD and hope to develop a novel diagnostic model to contribute to the early and precise diagnosis of childhood ASD.
Materials and Methods
Data Acquisition
The transcriptome data and clinical information of ASD and typically developing (TD) samples were downloaded from GSE18123 and GSE111176 datasets. The external validation data were downloaded from GSE113834 and GSE28521. The characteristics of the four datasets are detailed in Supplementary Table 1. Briefly, the sample size of normal control is 82, 126, 12, and 40, respectively (from left to right, GSE18123, GSE111176, GSE113834, and GSE28521); the sample size of ASD is 41, 119, 15, and 39, respectively (from left to right, GSE18123, GSE111176, GSE113834, and GSE28521). Tissues are sourced from blood (GSE18123 and GSE111176) or the brain (GSE113834, and GSE28521). The information on age, gender, time of collection of blood samples, and postmortem parameters was incomplete. The downloaded data were normalized by log2-transformed [log2(x + 1)]. All samples included in this study were accorded with the inclusive criteria: (1) Homo sapiens; (2) ASD and TD samples; (3) Availability of transcriptome data; and (4) Validation data including the expression of key genes included in the model. Three hundred thirty-five FRGs were obtained from the FerrDb database1. The rest FRGs were acquired from three published articles (16–18). After integrating these genes, we at last acquired 376 FRGs for further analysis.
Identification of Differentially Expressed FRGs
The “limma” package from R software was used to identify DEGs between ASD and TD groups. Only genes with P < 0.05 were considered as DEFRGs.
Functional Enrichment Analysis
The “clusterProfiler” R package was used to conduct Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses.
Establishment of a Ferroptosis-Related Diagnostic Model
After obtaining 27 differentially expressed FRGs, random forest analysis was performed to further identify key FRGs. The “RandomForest” R package was used to conduct random forest analysis. Parameters were set to default and genes with Gini values greater than 1 were considered key genes. After taking the intersection of GSE18123 and GSE111176, we obtained four key FRGs. Multiple logistic regression was applied to construct a ferroptosis-related diagnostic model. The diagnostic score was calculated using the equation: Score = AKR1C3 expression*(−0.2837) + CEBPG expression*(−0.7620) + DDIT4 expression*(−0.5943) + LAMP2 expression*(−0.9101) + 21.38.
Evaluation of the Ferroptosis-Related Diagnostic Model
Every individual from the Gene Expression Omnibus (GEO) database was allocated a diagnostic score derived from FRDM. Model specificity and sensitivity were assessed by calculating the area under the curve (AUC) values of ROC curves. The training cohort was GSE18123; the internal validation cohort was GSE111176; the external validation cohort was GSE113834 and GSE28521.
Regulatory Network
The microRNA (miRNA), long non-coding RNA (lncRNA), and transcription factor (TF) that interact with key genes were extracted from microRNA Data Integration Portal (mirDIP), starBase, and transcriptional regulatory relationships unraveled by sentence-based text mining (TRRUST) databases, respectively. Cytoscape software was used to display the multi-factor interaction network of key genes.
Molecular Docking
The protein structure information was downloaded from the PDBbind database2. Compounds’ structure information was downloaded from DrugBank3. AutoDock Vina was used to perform molecular docking, and compounds were filtered according to the binding energy. Screening parameters were set as Affinity < −7. PyMOL software was used to plot a protein-compound binding diagram.
Statistical Analysis
The R software (version 4.1.1)4 was used to perform all statistical analyses. Student t-test or one-way ANOVA was applied to assess differences between groups. P < 0.05 was considered statistically significant.
Results
Identification of Differentially Expressed FRGs and Enrichment Analysis
Forty four FRGs were extracted from the GSE18123 dataset (Figure 1A). One hundred twenty-one FRGs were extracted from the GSE111176 dataset (Figure 1B). After taking the intersection of FRGs from two datasets, we at last obtained 27 FRGs for further analysis (Figure 1C). The top 20 terms of the 27 FRGs’ enrichment analysis were displayed in four parts, including biological process (Figure 2A), cell component (Figure 2B), molecular function (Figure 2C), and KEGG (Figure 2D).
Figure 1. Identification of DEFRGs. Heatmap of DEFRGs from GSE18123 (A) and GSE111176 (B). (C) Venn diagram for DEFRGs.
Figure 2. Enrichment analysis. (A–C) GO enrichment analysis (BP, CC, and MF, respectively). (D) KEGG enrichment analysis.
Random Forest Analysis
For the differentially expressed FRGs, random forest analysis was performed to further identify key FRGs. The top 30 genes with high accuracy and Gini value were displayed in Figures 3A,B. Genes with Gini value >1 were selected and we obtained 28 genes in GSE18123 and 35 genes in GSE111176. Integrating the key genes from two datasets, we at last acquired 4 key FRGs (Supplementary Figure 1).
Figure 3. Random forest analysis. (A) Accuracy value and Gini value of candidate ferroptosis-related genes in GSE18123. (B) Accuracy and Gini values of candidate ferroptosis-related genes in GSE111176.
Expression of Ferroptosis-Related Genes in GSE18123 and GSE111176
AKR1C3, CEBPG, DDIT4, and LAMP2 were all highly expressed in ASD samples in both GSE18123 (Figure 4A) and GSE111176 datasets (Figure 4B). Additionally, in GSE18123, FRGs have positive correlations except for DDIT4 and LAMP2 genes (Supplementary Figure 2A). In GSE111176, AKR1C3 and DDIT4 have a positive correlation. CEBPG was negatively associated with DDIT4; while positively associated with LAMP2. DDIT4 and LAMP2 have a negative correlation (Supplementary Figure 2B).
Figure 4. Expression of FRGs in GSE18123 (A) and GSE111176 (B). *p < 0.05, **p < 0.01, ***p < 0.001.
Construction and Evaluation of a Ferroptosis-Related Diagnostic Model
Based on the four FRGs, we constructed diagnostic models using multiple logistic regression in GSE18123 and GSE111176, respectively. Then, ROC curves were used to evaluate the model specificity and sensitivity. The AUC value for the diagnostic model in GSE18123 is 0.7002; in GSE111176 is 0.7498 (Figure 5A). To further evaluate the reliability of the two models, external datasets were used to perform validations. For the model based on GSE18123, validation datasets including GSE111176, GSE113834, and GSE28521 were adopted and the AUC value is 0.7442, 0.7444, and 0.6474, respectively (Figure 5B). For the model based on GSE111176, validation datasets including GSE18123, GSE113834, and GSE28521 were adopted and the AUC value is 0.6838, 0.6944, and 0.6417, respectively (Figure 5C). Hence, the diagnostic model based on GSE18123 has a better diagnostic performance. The cutoff value of the diagnostic model was determined to be 0.9904 based on the Youden index (19); a diagnostic index ≤0.9904 indicated an ASD, and a diagnostic index >0.4 indicated its absence.
Figure 5. Evaluation of the constructed diagnostic model. (A) ROC curves in development datasets. (B) ROC curves in validation datasets based on the model from GSE18123. (C) ROC curves in validation datasets based on the model from GSE111176.
Multi-Factor Regulatory Network
According to the interaction degree, RORA, EAF1, and NFYB are the only three transcription factors that can interact with CEBPG, DDIT4, and LAMP2. In addition, hsa-miR-4703-3p can regulate LAMP2 and AKR1C3; hsa-miR-6073 can regulate DDIT4 and LAMP2 (Figure 6). These results indicated that TFs and miRNAs may play an important role in the development of ASD.
Figure 6. Multi-factor regulatory network. Orange dots represent lncRNAs, red dots represent transcription factors (TF), green dots represent miRNAs, blue dots represent mRNAs, and the size of dots represents the number of regulatory genes.
Potential Therapeutic Compounds
According to the PDBbind database, AKR1C3, LAMP2, and DDIT4 have available spatial structure information for subsequent analysis. We found only piperaquine could dock to all three FRGs with affinity < −7. The docking conformation of piperaquine to FRGs was displayed in Figure 7A (AKR1C3), Figure 7B (DDIT4), and Figure 7C (LAMP2). Figure 7D showed the compound structure of piperaquine. The docking score was shown in Supplementary Table 2.
Figure 7. Molecular docking. Docking conformation of DB13941 with AKR1C3 (A), DDIT4 (B), LAMP2 (C). The compound structure of DB13941 (D).
Discussion
At present, the ASD diagnosis was confirmed using clinical expert assessment including the Autism Diagnostic Interview-Revised (ADI-R) and the Autism Diagnostic Observation Schedule (ADOS), but only a limited number have been rigorously tested for diagnostic accuracy (20). In this study, based on transcriptome data, we constructed a ferroptosis-related diagnostic model to distinguish ASD patients from the TD population. We used three validation datasets to evaluate the precision of this model and the AUC value is 0.7442, 0.7444, and 0.6474, respectively. Hence, this diagnostic model may have the potential to be applied in the future.
To date, the four genes in this study have not been reported in ASD. CCAAT Enhancer Binding Protein Gamma (CEBPG) is a member of the C/EBP family and has the sequences required for DNA binding and heterodimer formation, but lacks the sequences required for transactivation (21). CEBPG appears as a stress-induced gene in some genome-wide expression studies, which can stimulate promoters of IL-6 and IL-8 in B cell lines (22). Previous studies have implicated CEBPG as an antioxidant regulator that controls redox homeostasis in normal and cancer cells (23). The neurobiology of ASD is considered related to oxidative stress (24) which indicated that CEBPG may be involved in the development of ASD. Aldo-Keto Reductase Family 1 Member C3 (AKR1C3) is known as a hydroxysteroid dehydrogenase, which exists in many normal human tissues at varying levels (25). The mRNA expression and enzyme activity of AKR1C3 in subcortical white matter were significantly higher than in the cerebral cortex, and AKR1C3 activity in adults was higher than in children, but no gender differences were observed (26). In addition, the antioxidant response element (ARE) could bind to the AKR1C3 promoter indicating that oxidative and electrophilic stress might regulate AKR1C3 expression (27). Due to the upregulation of AKR1C3 being closely correlated with various diseases, a large range of research on AKR1C3 inhibitors has been executed in the past few years, which could have the potential to be applied to ASD patients. Lysosomal Associated Membrane Protein 2 (LAMP2) is a major lysosomal membrane glycoprotein, which plays a role in multiple biological processes including antigen presentation, oxidative stress, and regulation of T lymphocyte responses (28). LAMP2 is ubiquitously expressed in the central nervous system and has been reported involved in the development of Parkinson’s disease and Alzheimer’s disease (29, 30). DNA Damage Inducible Transcript 4 (DDIT4), is known as a HIF1A responsive protein that promotes oxidative stress-dependent cell death (31). In addition, numerous studies have suggested that DDIT4 was crucial for optimal T cell proliferation and survival (32). Abnormal immune system regulation is involved in the pathophysiological process of ASD, including T cell-related signaling pathways (33). Overall, all these four genes are associated with oxidative stress or immune response, which provides clues for the exploration of the pathogenesis of ASD.
To elucidate the potential mechanism of the four genes in ASD, we conducted a multi-factor regulation network. We identified that RORA, EAF1, NFYB, miR-4703-3p, and miR-6073 might play a role in the development of ASD. RAR-Related Orphan Receptor A (RORA) is a ligand-dependent nuclear receptor that regulates gene transcription. Recently, studies have identified RORA as a novel candidate gene for ASD, which may be conducive to the known pathophysiology, behaviors, and sex bias of ASD (34). ELL Associated Factor 1 (EAF1) is one of the EAF family members, which plays an essential role in tumor suppression and embryogenesis (35). EAF1 and EAF2 can suppress Wnt/β-catenin signaling to affect neuroectodermal and mesodermal patterning (36). Nuclear Transcription Factor Y Subunit Beta (NFYB) plays a fundamental role in proliferation by binding to and regulating the transcription of numerous cell cycle regulatory genes (37). So far, there have been no reports on miR-4703-3p and miR-6073, and our study provides insights into the potential role of these miRNAs in ASD.
The classical therapeutic drugs for many diseases may have unexpected curative effects in other diseases, which is an approach known as drug repurposing. Piperaquine is a bisquinoline antimalarial drug, which is known as an effective drug to treat malaria, similar to chloroquine (38). We found that piperaquine has good docking scores with all three FRGs, which indicated piperaquine could be a promising drug for the treatment of ASD. The immune complex consists of brain tissue antigens and nerve-specific autoantibodies that can penetrate the blood-brain barrier and damage nerve tissue of children, resulting in cognitive, language development, and social communication disorders of children (39). Previous studies have shown that chloroquine has anti-inflammatory and antiviral effects, which are used to treat many diseases such as systemic lupus erythematosus (SLE), antiphospholipid antibody syndrome (APS), Middle-east respiratory syndrome, and HIV infection (40–42). Virus infection can cause irreversible damage to the central nervous system which produces ASD-related symptoms (43). Hence, piperaquine may play a role in the treatment of ASD through the immune-related pathway. In the future, studies need to be performed to explore the potential effect of piperaquine on ASD patients.
To a certain extent, several limitations of this study should not be ignored. (a) Clinical information was limited due to the retrospective data from the GEO database; (b) the number of samples in this study was relatively small, especially in the GSE113834 and GSE28521. In addition, the data from GSE113834 and GSE28521 was produced from the brain source. Hence, the conclusions generated based on GSE113834 and GSE28521 need to be further validated in the future; (c) ASD is a multifactorial disease, so the investigation of FRGs cannot comprehensively interpret the association of ASD risk.
Collectively, we identified four key FRGs as potential biomarkers and constructed a diagnostic model of childhood ASD, and derived potential criteria for ASD based upon this, that would need future validation. In addition, we found that piperaquine, previously a classic anti-malarial drug, may be worth exploring as a potential therapeutic compound for ASD treatment based on these findings.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
XW and RL designed the work, integrated and analyzed the data, and wrote the manuscript. QH and XC edited and revised the manuscript. All authors approved the manuscript.
Funding
This work was supported by the Science and Technology Development Fund Project of Nanjing Medical University (NMUB2019224).
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2022.886055/full#supplementary-material
Supplementary Figure 1 | Venn diagram of key FRGs.
Supplementary Figure 2 | The correlation between key FRG expression.
Footnotes
- ^ http://www.zhounan.org/ferrdb/index.html
- ^ http://www.pdbbind.org.cn/
- ^ https://www.drugbank.com/
- ^ https://www.r-project.org/
References
1. Orefice LL, Mosko JR, Morency DT, Wells MF, Tasnim A, Mozeika SM, et al. Targeting peripheral somatosensory neurons to improve tactile-related phenotypes in ASD models. Cell. (2019) 178:867–886.e24. doi: 10.1016/j.cell.2019.07.024
2. 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
3. Rubenstein JL, Merzenich MM. Model of autism: increased ratio of excitation/inhibition in key neural systems. Genes Brain Behav. (2003) 2:255–67.
4. Lord C, Brugha TS, Charman T, Cusack J, Dumas G, Frazier T, et al. Autism spectrum disorder. Nat Rev Dis Primers. (2020) 6:5.
5. Chen L, Shi XJ, Liu H, Mao X, Gui LN, Wang H, et al. Oxidative stress marker aberrations in children with autism spectrum disorder: a systematic review and meta-analysis of 87 studies (N = 9109). Transl Psychiatry. (2021) 11:15. doi: 10.1038/s41398-020-01135-3
6. Liu SH, Shi XJ, Fan FC, Cheng Y. Peripheral blood neurotrophic factor levels in children with autism spectrum disorder: a meta-analysis. Sci Rep. (2021) 11:15.
7. Baird G, Cass H, Slonims V. Diagnosis of autism. BMJ. (2003) 327:488–93. doi: 10.1136/bmj.327.7413.488
8. Lord C, Elsabbagh M, Baird G, Veenstra-Vanderweele J. Autism spectrum disorder. Lancet. (2018) 392:508–20.
9. Farmer C, Thurm A, Grant P. Pharmacotherapy for the core symptoms in autistic disorder: current status of the research. Drugs. (2013) 73:303–14. doi: 10.1007/s40265-013-0021-7
10. Chen PH, Wu J, Ding CC, Lin CC, Pan S, Bossa N, et al. Kinome screen of ferroptosis reveals a novel role of ATM in regulating iron metabolism. Cell Death Differ. (2020) 27:1008–22. doi: 10.1038/s41418-019-0393-7
11. Zhou Y, Shen Y, Chen C, Sui X, Yang J, Wang L, et al. The crosstalk between autophagy and ferroptosis: what can we learn to target drug resistance in cancer? Cancer Biol Med. (2019) 16:630–46. doi: 10.20892/j.issn.2095-3941.2019.0158
12. Jiang X, Stockwell BR, Conrad M. Ferroptosis: mechanisms, biology and role in disease. Nat Rev Mol Biol. (2021) 22:266–82. doi: 10.1038/s41580-020-00324-8
13. Zhou Y, Lin W, Rao T, Zheng J, Zhang T, Zhang M, et al. Ferroptosis and its potential role in the nervous system diseases. J Inflamm Res. (2022) 15:1555–74. doi: 10.2147/JIR.S351799
14. Ou M, Jiang Y, Ji Y, Zhou Q, Du Z, Zhu H, et al. Role and mechanism of ferroptosis in neurological diseases. Mol Metab. (2022). doi: 10.1016/j.molmet.2022.101502
15. Liu X, Lin J, Zhang H, Khan NU, Zhang J, Tang X, et al. Oxidative stress in autism spectrum disorder-current progress of mechanisms and biomarkers. Front Psychiatry. (2022) 13:813304. doi: 10.3389/fpsyt.2022.813304
16. Liang JY, Wang DS, Lin HC, Chen XX, Yang H, Zheng Y, et al. A novel ferroptosis-related gene signature for overall survival prediction in patients with hepatocellular carcinoma. Int J Biol Sci. (2020) 16:2430–41. doi: 10.7150/ijbs.45050
17. Zhuo S, Chen Z, Yang Y, Zhang J, Tang J, Yang K. Clinical and biological significances of a ferroptosis-related gene signature in glioma. Front Oncol. (2020) 10:590861. doi: 10.3389/fonc.2020.590861
18. Hong Z, Tang P, Liu B, Ran C, Yuan C, Zhang Y, et al. Ferroptosis-related genes for overall survival prediction in patients with colorectal cancer can be inhibited by gallic acid. Int J Biol Sci. (2021) 17:942–56. doi: 10.7150/ijbs.57164
19. Fluss R, Faraggi D, Reiser B. Estimation of the youden index and its associated cutoff point. Biom J. (2005) 47:458–72. doi: 10.1002/bimj.200410135
20. Hong SJ, Vos de Wael R, Bethlehem RAI, Lariviere S, Paquola C, Valk SL, et al. Atypical functional connectome hierarchy in autism. Nat Commun. (2019) 10:1022. doi: 10.1038/s41467-019-08944-1
21. Cooper C, Henderson A, Artandi S, Avitahl N, Calame K. Ig/EBP (C/EBP gamma) is a transdominant negative inhibitor of C/EBP family transcriptional activators. Nucleic Acids Res. (1995) 23:4371–7. doi: 10.1093/nar/23.21.4371
22. Mullins DN, Crawford EL, Khuder SA, Hernandez DA, Yoon Y, Willey JC. CEBPG transcription factor correlates with antioxidant and DNA repair genes in normal bronchial epithelial cells but not in individuals with bronchogenic carcinoma. BMC Cancer. (2005) 5:141. doi: 10.1186/1471-2407-5-141
23. Huggins CJ, Mayekar MK, Martin N, Saylor KL, Gonit M, Jailwala P, et al. C/EBPγ is a critical regulator of cellular stress response networks through heterodimerization with ATF4. Mol Cell Biol. (2015) 36:693–713. doi: 10.1128/MCB.00911-15
24. 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.
25. Evans K, Duan J, Pritchard T, Jones CD, McDermott L, Gu Z, et al. OBI-3424, a novel AKR1C3-activated prodrug, exhibits potent efficacy against preclinical models of T-all. Clin Cancer Res. (2019) 25:4493–503. doi: 10.1158/1078-0432.CCR-19-0551
26. Steckelbroeck S, Watzka M, Stoffel-Wagner B, Hans VH, Redel L, Clusmann H, et al. Expression of the 17beta-hydroxysteroid dehydrogenase type 5 mRNA in the human brain. Mol Cell Endocrinol. (2001) 171:165–8. doi: 10.1016/s0303-7207(00)00432-9
27. Tebay LE, Robertson H, Durant ST, Vitale SR, Penning TM, Dinkova-Kostova AT, et al. Mechanisms of activation of the transcription factor Nrf2 by redox stressors, nutrient cues, and energy status and the pathways through which it attenuates degenerative disease. Free Radic Biol Med. (2015) 88:108–46. doi: 10.1016/j.freeradbiomed.2015.06.021
28. Valdor R, Mocholi E, Botbol Y, Guerrero-Ros I, Chandra D, Koga H, et al. Chaperone-mediated autophagy regulates T cell responses through targeted degradation of negative regulators of T cell activation. Nat Immunol. (2014) 15:1046–54. doi: 10.1038/ni.3003
29. Boman A, Svensson S, Boxer A, Rojas JC, Seeley WW, Karydas A, et al. Distinct lysosomal network protein profiles in parkinsonian syndrome cerebrospinal fluid. J Parkinsons Dis. (2016) 6:307–15. doi: 10.3233/JPD-150759
30. Armstrong A, Mattsson N, Appelqvist H, Janefjord C, Sandin L, Agholme L, et al. Lysosomal network proteins as potential novel CSF biomarkers for Alzheimer’s disease. Neuromolecular Med. (2014) 16:150–60. doi: 10.1007/s12017-013-8269-3
31. Slattery ML, John EM, Torres-Mejia G, Lundgreen A, Lewinger JP, Stern MC, et al. Angiogenesis genes, dietary oxidative balance and breast cancer risk and progression: the breast cancer health disparities study. Int J Cancer. (2014) 134:629–44. doi: 10.1002/ijc.28377
32. Reuschel EL, Wang J, Shivers DK, Muthumani K, Weiner DB, Ma Z, et al. REDD1 is essential for optimal T cell proliferation and survival. PLoS One. (2015) 10:e0136323. doi: 10.1371/journal.pone.0136323
33. Ashwood P, Krakowiak P, Hertz-Picciotto I, Hansen R, Pessah IN, Van de Water J. Altered T cell responses in children with autism. Brain Behav Immun. (2011) 25:840–9. doi: 10.1016/j.bbi.2010.09.002
34. Sarachana T, Hu VW. Differential recruitment of coregulators to the RORA promoter adds another layer of complexity to gene (dys) regulation by sex hormones in autism. Mol Autism. (2013) 4:39. doi: 10.1186/2040-2392-4-39
35. Pascal LE, Su F, Wang D, Ai J, Song Q, Wang Y, et al. Conditional deletion of eaf1 induces murine prostatic intraepithelial neoplasia in mice. Neoplasi. (2019) 21:752–64. doi: 10.1016/j.neo.2019.05.005
36. Liu JX, Zhang D, Xie X, Ouyang G, Liu X, Sun Y, et al. Eaf1 and Eaf2 negatively regulate canonical Wnt/β-catenin signaling. Development. (2013) 140:1067–78. doi: 10.1242/dev.086157
37. Gurtner A, Manni I, Piaggio G. NF-Y in cancer: impact on cell transformation of a gene essential for proliferation. Biochim Biophys Acta Gene Regul Mech. (2017) 1860:604–16. doi: 10.1016/j.bbagrm.2016.12.005
38. Li R, Ding Z, Jin P, Wu S, Jiang G, Xiang R, et al. Development and validation of a novel prognostic model for acute myeloid leukemia based on immune-related genes. Front Immunol. (2021) 12:639634. doi: 10.3389/fimmu.2021.639634
39. Malani Shukla N, Lotze TE, Muscal E. Inflammatory diseases of the central nervous system. Neurol Clin. (2021) 39:811–28.
40. Bakadia BM, He F, Souho T, Lamboni L, Ullah MW, Boni BO, et al. Prevention and treatment of COVID-19: focus on interferons, chloroquine/hydroxychloroquine, azithromycin, and vaccine. Biomed Pharmacother. (2021) 133:111008. doi: 10.1016/j.biopha.2020.111008
41. Varisli L, Cen O, Vlahopoulos S. Dissecting pharmacological effects of chloroquine in cancer treatment: interference with inflammatory signaling pathways. Immunology. (2020) 159:257–78. doi: 10.1111/imm.13160
42. Davis TM, Hung TY, Sim IK, Karunajeewa HA, Ilett KF. Piperaquine: a resurgent antimalarial drug. Drugs. (2005) 65:75–87. doi: 10.2165/00003495-200565010-00004
Keywords: autism spectrum disorder, diagnosis, ferroptosis, random forest, piperaquine
Citation: Wu X, Li R, Hong Q and Chi X (2022) Development and Validation of a Novel Diagnostic Model for Childhood Autism Spectrum Disorder Based on Ferroptosis-Related Genes. Front. Psychiatry 13:886055. doi: 10.3389/fpsyt.2022.886055
Received: 28 February 2022; Accepted: 26 April 2022;
Published: 12 May 2022.
Edited by:
David Quentin Beversdorf, University of Missouri, United StatesReviewed by:
Yong Cheng, Minzu University of China, ChinaQiyang Li, Southern Medical University, China
Copyright © 2022 Wu, Li, Hong and Chi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xia Chi, chixiananjing@163.com; Qin Hong, rambler_hq@163.com
†These authors have contributed equally to this work