- 1Department of Rehabilitation Medicine, Xiangyang No.1 People’s Hospital, Hubei University of Medicine, Xiangyang, China
- 2Department of Psychiatric Medicine, Xiangyang Central Hospital, Hubei University of Arts and Science, Xiangyang, China
- 3Department of Clinical Laboratory, Xiangyang No.1 People’s Hospital, Hubei University of Medicine, Xiangyang, China
Background: Observational studies have reported a strong association between autistic spectrum disorder (ASD) and intestinal metabolites. However, it is unclear whether this correlation is causally or violated by confounding or backward causality. Therefore, this study explored the potential causal relationship between intestinal metabolites and dependent metabolites on ASD.
Methods: We used a two-sample Mendelian random analysis and selected variants closely related to intestinal flora-dependent metabolites as instrumental variables. MR-Egger, inverse variance weighted (IVW), MR-PRESSO, maximum likelihood, and weighted median were performed to reveal their causal relationships. Ten metabolites were studied, which included trimethylamine-N-oxide, betaine, carnitine, choline, glutamate, kynurenine, phenylalanine, serotonin, tryptophan, and tyrosine. Sensitivity tests were also performed to evaluate the robustness of the MR study.
Results: The IVW method revealed that serotonin may increase the ASD risk (OR 1.060, 95% CI: 1.006–1.118), while choline could decrease the ASD risk (OR 0.925, 95% CI: 0.868–0.988). However, no definite causality was observed between other intestinal metabolites (e.g., trimethylamine-N-oxide, betaine, and carnitine) with ASD. Additionally, neither the funnel plot nor the MR-Egger test showed horizontal pleiotropy, and the MR-PRESSO test found no outliers. Cochran’s Q test showed no significant heterogeneity among the studies, suggesting the robustness of the study.
Conclusion: Our study found potential causality from intestinal metabolites on ASD. Clinicians are encouraged to offer preventive measures to such populations.
1. Introduction
In the groundbreaking clinical description based on 11 boys with emotional contact disorder, Kanner first identified autistic spectrum disorder (ASD) as a coherence disorder in 1943 (1). With the developing understanding of cognitive and behavioral disorders as organic and cerebral pathologies, the concept of this disease in psychiatry has changed significantly. In the 1960s and 1970s, ASD was considered a form of psychosis similar to childhood schizophrenia. Then, it was primarily related to the parenting style (2). Epidemiological studies have shown that the prevalence of ASD was about 4 out of 10,000 children. However, with the passage of time and the update of diagnostic criteria, this estimate had increased from 1/1,000 in 1988 to 1/150–1/200 in 2002 (3). The prevalence was 16.8 per 1,000 children aged 8 years (4). Additionally, with the continuous development of technology, understanding of the mechanism of ASD had broadened. For example, some rare genetic or biomedical diseases, such as phenylketonuria, may lead to ASD (5). Although those various potential interventions have been supported by quality research, only a few pharmacological interventions have been approved by FDA for ASD. This mainly arose due to the absence of adequate randomized controlled studies. Moreover, most of these drugs target maladaptive behavior instead of core symptoms. Finding modifiable risk factors for ASD is of great significance.
The Connectivome Theory proposed by Leonardo et al. shows that ASD is also related to many body devices such as the immune system, motor system and brain-gut axis (6). As a unique ecosystem, the intestinal flora is rich in numerous species and widely participate in various physiological activities essential to human health. Accumulating evidence from observational studies or animal studies suggested that intestinal flora was correlated with various human diseases, including various cancers, cardiovascular diseases, autoimmune diseases, and neuropsychiatric diseases. Intestinal flora could affect brain activity such as behavior by regulating the intestinal metabolites. The research of intestinal flora regulating brain and behavior involves many aspects, such as the intestinal nervous system (7), neuroimaging (8), the interaction between intestinal flora and host (9, 10), and intestinal flora-intestinal brain axis (11–13). These studies demonstrated that intestinal flora affects the brain and associated behavior through the intestinal-brain axis. Additionally, it has been reported that intestinal flora is closely related to ASD (14). However, previous associations were suggested by observational studies. Since intestinal flora could be altered with diet behavior change or diseases, it remained unclear whether the association was causal or not. Besides, whether the association was violated by reverse causality or confounding remained undetermined.
As an emerging technique, Mendelian randomization uses the genetic variation of risk factors to evaluate the causality between risk factors and specific diseases (13, 15, 16). There were large amounts of MR analyses successfully revealing the potential risk factors for many neuropsychiatric diseases (17, 18). Mendelian randomization is conceptually similar to a prospective randomized controlled trial (RCT), while the Mendelian randomization method can be conducted retrospectively. In addition, since all genetic variations are determined during pregnancy and occur before the onset of the disease, Mendelian randomization can avoid non-differential measurement error and potential deviation. With large sample sizes in GWAS data booming, MR is a suitable statistical analysis to reveal the causal effect of human intestinal metabolites on ASD.
Herein, this study discusses the potential causality between internal metabolites and ASD through Mendelian randomization analysis. Our findings can fill the gaps between human intestinal metabolites and ASD and shed light on ASD prevention through diet behavior management or drug administration.
2. Materials and methods
2.1. Ethical approval, data availability, and report
All data used in this study were acquired from public databases, which had been approved by their respective committees. Thus, this study requires no further ethical approval.
2.2. Data source
A total of 10 intestinal metabolites were noted in our study, including trimethylamine-N-oxide, betaine, carnitine, choline, glutamate, kynurenine, phenylalanine, serotonin, tryptophan, and tyrosine. These intestinal metabolites play important roles in maintaining healthy neuropsychiatric function (19, 20), thus their causal effects on ASD are most concerning. Data about the 10 intestinal metabolites were acquired from the results of a genome-wide association study (GWAS) based on the pooled data from the human metabolome conducted among 2,076 European participants in the Framingham Heart Study (21).
The data on ASD was obtained from the dataset published by the cross-disease group of the psychiatric genomics federation. The sample size was 61,220, including 33,332 patients with ASD and 27,888 control population. All participants in the above two GWAS were of European ancestry.
The instrument variables were identified in a GWAS of human intestinal metabolites with the following standard: (1) to obtain more SNPs as instrument variables, we selected SNPs correlated with the risk factors at a relatively relaxed statistical significance threshold (P < 5 × 10–5), (2) to remove SNPs at linkage disequilibrium R2 > 0.001 based on European ancestry reference data from the 1,000 Genomes project, and (3) to calculate the F statistic for each SNP and discard weak instrumental SNPs (F < 10). The instrument variables were used to investigate the causal effect of metabolome on ASD. That is, if a causal effect was found from instrument variables on ASD, the causal effect from human intestinal metabolites on ASD can be concluded.
2.3. Statistical analysis
Two-sample MR package of R software was used for Mendelian randomization analysis. In this study, we applied several MR approaches, including the inverse-variance weighted (IVW) method, MR-Egger, MR-PRESSO, maximum likelihood, and weighted median. The traditional IVW method was used as the main finding to estimate the association between intestinal metabolites and ASD. However, since IVW estimation may be affected by invalid instrument deviation or pleiotropy, the effectiveness and robustness of IVW results were tested using other methods, including MR-Egger, MR-PRESSO, maximum likelihood, and weighted median. To reduce the deviation caused by horizontal pleiotropy (affecting the results through causal pathways rather than exposure), MR-PRESSO (pleiotropy residuals and outliers) was used to detect extensive horizontal pleiotropy in all results. Outliers identified by MR-PRESSO would be discarded and MR analyses would be re-performed. Finally, to evaluate the robustness of the results, heterogeneity and pleiotropy tests were carried out for the statistically significant results, including the MR-Egger intercept test and modified Cochran’s Q statistics. A p-value of > 0.05 is expected in heterogeneity and pleiotropy analyses to support our study is unbiased. Funnel plots and leave-one-out analyses were also performed for significant MR results. Symmetrical funnel plots would be expected, which can reflect the absence of outliers. The perfect results of leave-one-out analyses would be for all estimates to remain similar and significant after each SNP is removed.
All statistical tests were conducted in a two-tailed manner. Since the results were binary variables, the estimated value of the effects was further converted into an odds ratio (OR) to demonstrate the correlation of metabolites with ASD. To account for multiple testing in our study, the Bonferroni’s-corrected significance level of P < 5 × 10–3 (0.05 divided by 10 risk factors) was used. The p-value between 5 × 10–3 and 0.05 was considered a potential association.
3. Results
3.1. Determination of instrumental variables and their weak bias
In total, 2,117 SNPs related to intestinal metabolites (77 associated with beta-hydroxybutyric acid, 264 associated with betaine, 159 associated with carnitine, 105 associated with choline, 133 associated with glutamate, 257 associated with kynurenine, 161 associated with phenylalanine, 52 associated with propionic acid, 216 associated with serotonin, 293 associated with trimethylamine-N-oxide, 170 associated with tryptophan, and 230 associated with tyrosine) were extracted from GWAS and used as IVs (R2 < 0.001, P < 5 × 10–5). The F statistics of involved SNPs were greater than 10, indicating that the instrument variables used in our study are strong. Therefore, no weak bias was detected in the results, and the conclusions of this study were reliable (Supplementary Table 1).
3.2. Causal effects from intestinal metabolites on ASD risk
Results of IVW supported that serotonin could increase the ASD risk (OR = 1.060, 95% CI: 1.006–1.118). Meanwhile, the results of IVW showed that choline might reduce the ASD risk (OR = 0.925, 95% CI: 0.868–0.988). Similar results were observed in other MR analyses (Table 1). Both maximum likelihood and MR-PRESSO present similar and significant estimates, while MR-Egger and weighted median present similar but insignificant estimates. No evidence of causal effects from other intestinal metabolites on ASD was found.
3.3. Sensitivity analysis
To confirm the stability of the above results, MR-Egger intercept analyses and MR-PRESSO tests were performed on the involved SNP variants. No potential horizontal pleiotropy was detected (P > 0.05), and the funnel chart did not show any deviation in this study (Supplementary Figures 11–20). The revised Cochran’s Q statistics showed no significant heterogeneity in the effect of involved SNP (P = 0.90). Additionally, sensitivity analysis was conducted to verify the impact of each SNP site on the overall causality. When a single SNP was systematically removed and MR analysis was repeated, no significant difference was observed (Supplementary Figures 1–10), proving that the estimated effect cannot be explained by any single genetic tool.
4. Discussion
This study is the first Mendelian randomization study to systematically evaluate the causality between intestinal metabolites and ASD. Our results supported the idea that serotonin may be related to increased ASD risk, while choline may reduce ASD risk. The relationship between the two intestinal metabolites and ASD was further confirmed by sensitivity analysis. In our study, we found no evidence of causal effects from other metabolites on ASD.
Central nervous system development disorders have long been associated with intestinal metabolites, creating the so-called “gut-brain axis.” ASD is one of the most frequent neurodevelopmental disorders associated with dysregulation of the microbiota-gut-brain axis. Metabolism is one of the major drivers of microbiome-host interactions. Therefore, metabolites were hypothesized to affect the risk of ASD. Previous observational studies have suggested associations between metabolites and ASD. However, it is unclear whether intestinal metabolites cause ASD or whether ASD causes changes in intestinal metabolites. Herein, we applied two-sample MR to investigate their causal effects.
The study result was in concurrence with previous studies (i.e., serotonin may increase ASD risk). A transmission/disequilibrium test study conducted on 184 cases showed that serotonin was not related to ASD risk. Another meta-analysis, including 22 observational studies, showed that serotonin might promote ASD risk (OR = 4.6, 95% CI = 3.1–5.2) (22). Overall, the exact mechanism by which serotonin causes ASD remains unclear. However, few studies have shown that the onset of ASD may be related to the function of serotonin. For instance, serotonin can regulate neuron differentiation and migration, axon growth, myelin, and synapse formation (23). Additionally, serotonin connects the intestine with the central nervous system. It also affects neurogenesis, guiding the development and survival of post-mitotic neurons [dopamine, γ-Aminobutyric acid (GABA), and calcitonin gene-related peptide-secreting neurons] (24). However, in early childhood, the blood-brain barrier is often not completed until the age of 2 years (25). Therefore, during the fetal and neonatal period, 5-HT with a high level in the blood can enter the developing brain. These high levels may inhibit 5-HT-secreting neurons through the negative feedback mechanism and lead to a decreased serotonin concentration. Azmita et al. conducted an immunocytochemical analysis, in which they found an increase of serotoninergic axons in the cerebral cortex of patients with ASD, which indirectly confirms the above hypothesis. Nevertheless, these theories cannot explain the serotonin-mediated ASD mechanism (25). Additionally, other intestinal metabolites (e.g., choline) may reduce the ASD risk, as confirmed in other studies (26); however, some reports have suggested that choline increases ASD risk (27). Trimethylamine-N-oxide, an important molecule affecting human health, was thought to be a waste product of choline metabolism without action in our organism (28). A previous study reported that TMAO was found to increase the risk of ASD (28). However, evidence in our study supports the finding that choline reduces the risk of ASD but trimethylamine-N-oxide has no effect on ASD. Therefore, there might exist another pathway by which choline acts on ASD, and the specific mechanism remains to be explored.
In MR analyses, three assumptions should be satisfied to derive robust estimates. First, the instrument variables should be strongly related to exposure. In our study, to obtain more SNPs as instrument variables, we used a relaxed P-value. But the F statistics had been calculated to support the presence of strong instrument variables. Second, the instrument variables should not affect the outcomes through other pathways rather than exposure. This assumption is easy to fulfill for the randomly assigned SNPs during pregnancy. Because there were no well-established risk factors for ASD, we did not perform additional analyses to reveal the association between instrument variables and potential risk factors. Third, the instrument variables should not correlate with the outcomes, and we have checked that each SNP of the instrument variables in our study was not correlated with ASD.
In our study, we used IVW estimates as the main findings and other MR approaches as complementary analyses. For the significant estimates identified by IVW, both maximum likelihood and MR-PRESSO also derive significant estimates, but insignificant results are found in MR-Egger and the weighted median. Because MR-Egger suggests that all SNPs used in the study are invalid and the weighted median suggests that half SNPs are invalid. This basic theory causes these two approaches more difficult to derive significant estimates than the other three approaches. We have searched the literature and found support from other studies to believe in the IVW estimates and draw a significant conclusion. In addition, we did not use multiple tests corrected P-value in our study. Our purpose is to discover any potential risk factors for ASD, and Bonferroni’s or FDR correction would reduce our opportunity to do so. With MR analyses and sensitivity analyses, we identified that choline and serotonin were correlated with ASD risk. Further studies should be performed to discover the underlying mechanisms.
Taking the results that serotonin and choline are correlated with the risk of ASD robust, we should carry out some efficient protocols to adjust these two metabolites. Antibiotic administration is one useful method that can prevent the absorption of these two metabolites from the gastric and intestinal tracts. There were some non-absorbable antibiotics (29), for example, vancomycin, neomycin, and bacitracin. These non-absorbable antibiotics offer the advantage of killing the microbiota in the gut while avoiding the disadvantage of causing reverse effects on the human body after entering systemic circulation. Another efficient method is fecal microbiota transplant, which involves the transfer of intestinal microbiota from one human being to another. The fecal microbiota transplant is becoming more and more popular to prevent or treat many diseases. The procedure of fecal microbiota transplant is mature, and some studies have shown successful transfer of microbiota (30). All these methods might be useful to alter the concentration of serotonin and choline, and thus prevent the occurrence of ASD.
There were several strengths and novelties to our study. First, this is the first MR study to investigate the causal relationship between 10 human intestinal metabolites and the risk of an ASD. The MR framework is superior to most observational studies in addressing confounding risk factors, which makes our findings robust. Notably, the results of this study represent the relationship between intestinal metabolites and the lifetime risk of ASD, as genetic variation does not change during a lifetime. Meanwhile, since alleles are randomly classified and determined during pregnancy, this study did not observe deviation due to confusion and backward causality. Second, our findings are robust because of the large sample size. The larger the sample size, the more robust the MR findings. The GWAS of ASD enrolled more than 60,000 participants, tremendously larger than most observational studies.
The study had some limitations. First, this study included only people in Europe. So, it may not represent another ethnicity, such as people in Asia and Africa. Nevertheless, it can effectively reduce deviation caused by population stratification. Second, subgroup analysis could not be performed due to the lack of detailed clinical information. Third, this study showed the potential causality between intestinal metabolites and ASD. However, more extensive research is required to confirm the definite causality between intestinal metabolites and ASD.
5. Conclusion
Mendelian randomization analysis revealed that serotonin could increase the ASD risk. Therefore, managing children’s diet should be an essential component to reduce ASD risk. Additionally, due to the genetic differences among different nationalities, countries, and regions, it is necessary to explore this relationship in wider and more diverse populations.
Data availability statement
The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
DL and DB: study design. HL: data collection and data analysis. QW and XD: data interpretation. DL and XF: drafting manuscript. All authors take responsibility for the integrity of the data analysis and approved the final version of the manuscript.
Acknowledgments
We thank the Framingham Heart Study, the cross-disease group of the psychiatric genomics federation, for providing GWAS data.
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.1034214/full#supplementary-material
References
3. Autism and Developmental Disabilities Monitoring Network Surveillance Year 2002 Principal Investigators, Centers for Disease Control and Prevention. Prevalence of autism spectrum disorders-autism and developmental disabilities monitoring network, 14 sites, United States, 2002. MMWR Surveill Summ. (2007) 56:12–28.
4. Baio J, Wiggins L, Christensen DL, Maenner MJ, Daniels J, Warren Z, et al. Prevalence of autism spectrum disorder among children aged 8 years - autism and developmental disabilities monitoring network, 11 sites, United States, 2014. MMWR Surveill Summ. (2018) 67:1–23. doi: 10.15585/mmwr.ss6706a1
5. Geschwind DH. Advances in autism. Annu Rev Med. (2009) 60:367–80. doi: 10.1146/annurev.med.60.053107.121225
6. Zoccante L, Ciceri ML, Gozzi LA, Gennaro GD, Zerman N. The “Connectivome Theory”: a new model to understand autism spectrum disorders. Front Psychiatry. (2022) 12:794516. doi: 10.3389/fpsyt.2021.794516
7. Furness JB. The enteric nervous system and neurogastroenterology. Nat Rev Gastroenterol Hepatol. (2012) 9:286–94. doi: 10.1038/nrgastro.2012.32
8. Mayer EA, Aziz Q, Coen S, Kern M, Labus JS, Lane R, et al. Brain imaging approaches to the study of functional GI disorders: a Rome working team report. Neurogastroenterol Motil. (2009) 21:579–96. doi: 10.1111/j.1365-2982.2009.01304.x
9. Round JL, Mazmanian SK. The gut microbiota shapes intestinal immune responses during health and disease. Nat Rev Immunol. (2009) 9:313–23. doi: 10.1038/nri2515
10. Raybould HE. Gut microbiota, epithelial function, and derangements in obesity. J Physiol. (2012) 590:441–6. doi: 10.1113/jphysiol.2011.222133
11. Collins SM, Bercik P. The relationship between intestinal microbiota and the central nervous system in normal gastrointestinal function and disease. Gastroenterology. (2009) 136:2003–14. doi: 10.1053/j.gastro.2009.01.075
12. Collins SM, Surette M, Bercik P. The interplay between the intestinal microbiota and the brain. Nat Rev Microbiol. (2012) 10:735–42. doi: 10.1038/nrmicro2876
13. Bercik P. The microbiota-gut-brain axis: learning from intestinal bacteria? Gut. (2011) 60:288–9. doi: 10.1136/gut.2010.226779
14. Horn J, Mayer DE, Chen S, Mayer EA. Role of diet and its effects on the gut microbiome in the pathophysiology of mental disorders. Transl Psychiatry. (2022) 12:164. doi: 10.1038/s41398-022-01922-0
16. Nievergelt C, Maihofer A, Klengel T, Atkinson EG, Chen CY, Choi KW, et al. International meta-analysis of PTSD genome-wide association studies identifies sex- and ancestry-specific genetic risk loci. Nat Commun. (2019) 10:4558. doi: 10.1038/s41467-019-12576-w
17. Polimanti R, Ratanatharathorn A, Maihofer AX, Choi KW, Stein MB, Morey RA, et al. Association of economic status and educational attainment with posttraumatic stress disorder: a Mendelian randomization study. JAMA Netw Open. (2019) 2:e193447. doi: 10.1001/jamanetworkopen.2019.3447
18. Colodro-Conde L, Couvy-Duchesne B, Whitfield JB, Streit F, Gordon S, Kemper KE, et al. Association between population density and genetic risk for schizophrenia. JAMA Psychiatry. (2018) 75:901–10. doi: 10.1001/jamapsychiatry.2018.1581
19. Sabokdast M, Habibi-Rezaei M, Moosavi-Movahedi AA, Ferdousi M, Azimzadeh-Irani E, Poursasan N. Protection by beta-Hydroxybutyric acid against insulin glycation, lipid peroxidation and microglial cell apoptosis. Daru. (2015) 23:42. doi: 10.1186/s40199-015-0126-5
20. Zalar B, Haslberger A, Peterlin B. The role of microbiota in depression - a brief review. Psychiatr Danub. (2018) 30:136–41. doi: 10.24869/psyd.2018.136
21. Rhee EP, Ho JE, Chen MH, Shen D, Cheng S, Larson MG, et al. A genome-wide association study of the human metabolome in a community-based cohort. Cell Metab. (2013) 18:130–43. doi: 10.1016/j.cmet.2013.06.013
22. Gabriele S, Sacco R, Persico A. Blood serotonin levels in autism spectrum disorder: a systematic review and meta-analysis. Eur Neuropsychopharmacol. (2014) 24:919–29. doi: 10.1016/j.euroneuro.2014.02.004
23. Gaspar P, Cases O, Maroteaux L. The developmental role of serotonin: news from mouse molecular genetics. Nat Rev Neurosci. (2003) 4:1002–12. doi: 10.1038/nrn1256
24. Li ZS, Pham TD, Tamir H, Chen JJ, Gershon MD. Enteric dopaminergic neurons: definition, developmental lineage, and effects of extrinsic denervation. J Neurosci. (2004) 24:1330–9. doi: 10.1523/JNEUROSCI.3982-03.2004
25. Jaiswal P, Mohanakumar K, Rajamma U. Serotonin mediated immunoregulation and neural functions: complicity in the aetiology of autism spectrum disorders. Neurosci Biobehav Rev. (2015) 55:413–31. doi: 10.1016/j.neubiorev.2015.05.013
26. Gabis LV, Ben-Hur R, Shefer S, Jokel A, Shalom DB. Improvement of language in children with autism with combined donepezil and choline treatment. J Mol Neurosci. (2019) 69:224–34. doi: 10.1007/s12031-019-01351-7
27. Peralta-Marzal L, Prince N, Bajic D, Roussin L, Naudon L, Rabot S, et al. The impact of gut microbiota-derived metabolites in autism spectrum disorders. Int J Mol Sci. (2021) 22:10052. doi: 10.3390/ijms221810052
28. Quan L, Yi J, Zhao Y, Zhang F, Shi XT, Feng Z, et al. Plasma trimethylamine-N-oxide, a gut microbe-generated phosphatidylcholine metabolite, is associated with autism spectrum disorders. Neurotoxicology. (2020) 76:93–8. doi: 10.1016/j.neuro.2019.10.012
29. Cryan JF, O’Riordan KJ, Cowan CSM, Sandhu KV, Bastiaanssen TFS, Boehme M, et al. The microbiota-gut-brain axis. Physiol Rev. (2019) 99:1877–2013. doi: 10.1152/physrev.00018.2018
Keywords: autistic spectrum disorder, intestinal metabolites, Mendelian randomization, serotonin, gut-brain axis
Citation: Liu D, Bu D, Li H, Wang Q, Ding X and Fang X (2023) Intestinal metabolites and the risk of autistic spectrum disorder: A two-sample Mendelian randomization study. Front. Psychiatry 13:1034214. doi: 10.3389/fpsyt.2022.1034214
Received: 01 September 2022; Accepted: 22 December 2022;
Published: 12 January 2023.
Edited by:
Leonardo Zoccante, Integrated University Hospital Verona, ItalyReviewed by:
Xiong Chen, Guangzhou Medical University, ChinaHan Wang, The Affiliated Hospital of Southwest Jiaotong University, China
Copyright © 2023 Liu, Bu, Li, Wang, Ding and Fang. 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: Xiaolu Fang, MTUyNzE5ODYxOThAMTYzLmNvbQ==
†These authors have contributed equally to this work