
==== Front
Heliyon
Heliyon
Heliyon
2405-8440
Elsevier

S2405-8440(24)12603-5
10.1016/j.heliyon.2024.e36572
e36572
Research Article
Exploring the Missing link between vitamin D and autism spectrum disorder: Scientific evidence and new perspectives
Gao Tianci a1351468@outlook.com
a
Dang Wenjun 307692456@qq.com
b
Jiang Zhimei 304532709@qq.com
cd
Jiang Yuwei jyuwei@jmsu.edu.cn
cd⁎
a College of Clinical Medicine, Jiamusi University, Hei longJiang Province, China
b Jiamusi College, HeiLongJiang University of Chinese Medicine, Hei longJiang Province, China
c College of Rehabilitation Medicine, Jiamusi University, Hei longJiang Province, China
d Child Neurological Rehabilitation Key Laboratory of Heilongjiang province, China
⁎ Corresponding author. College of Rehabilitation Medicine, Jiamusi University, Hei longJiang Province, China. jyuwei@jmsu.edu.cn
22 8 2024
15 9 2024
22 8 2024
10 17 e3657214 2 2024
30 7 2024
19 8 2024
© 2024 The Authors. Published by Elsevier Ltd.
2024

https://creativecommons.org/licenses/by-nc-nd/4.0/ This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
Aim

This study aims to address the key question of the causal relationship between serum levels of 25-hydroxyvitamin D (vitamin D) and autism spectrum disorders (ASD).

Methods

Publicly available Genome-Wide Association Study (GWAS) datasets were used to conduct the bidirectional Two-sample MR analyses using methods including inverse-variance weighted (IVW), weighted median, MR-Egger regression, simple mode, MR-PRESSO test, Steiger filtering, and weighted mode, followed by BWMR for validation.

Results

The MR analysis indicated that there was no causal relationship between Vitamin D as the exposure and ASD as the outcome in the positive direction of the MR analysis (IVW: OR = 0.984, 95 % CI: 0.821–1.18, P = 0.866). The subsequent BWMR validation stage yielded consistent results (OR = 0.984, 95 % CI 0.829–1.20, P = 0.994). Notably, in the reverse MR analysis with ASD as the exposure and Vitamin D as the outcome, the results suggested that the occurrence of ASD could lead to decreased Vitamin D levels (IVW: OR = 0.976, 95 % CI: 0.961–0.990, P = 0.000855), with BWMR findings in the validation stage confirming the discovery phase (OR = 0.975, 95 % CI: 0.958–0.991, P = 0.00297). For the positive MR analysis, no pleiotropy was detected in the instrumental variables. Similarly, no pleiotropy or heterogeneity was detected in the instrumental variables for the reverse MR analysis. Sensitivity analysis using the leave-one-out approach for both positive and reverse instrumental variables suggested that the MR analysis results were robust.

Conclusion

Through the discovery and validation analysis process, we can confidently assert that there is no causative link between Vitamin D and ASD, and that supplementing Vitamin D is not expected to provide effective improvement for patients with ASD. Our study significantly advances a new perspective in ASD research and has a positive impact on medication guidance for patients with ASD.

Keywords

Autism spectrum disorder
25-Hydroxyvitamin D
Bayes Theorem
Mendelian randomization analysis
Epigenomics
==== Body
pmc1 Introduction

ASD is a neurodevelopmental disorder intricately linked to atypical brain development. Impairments in regions of the brain governing social interaction, communication, and cognitive behavior could contribute to the multifaceted features of ASD [1]. Symptoms vary widely among individuals and typically encompass repetitive behaviors, stereotypical movements, challenges in social engagement, language deficits, narrow interest scopes, and an acute sensitivity to environmental changes [2]. ASD is found worldwide, yet its prevalence differs across regions and countries. The Centers for Disease Control and Prevention (CDC) reported in 2018 that one in 59 children (one in 37 boys and one in 151 girls) is diagnosed with ASD, with an increasing trend observed year by year [3]. Vitamin D, a crucial nutrient, is widely recognized for its significant connection to various mental health disorders. Recent observational studies have revealed lower levels of vitamin D in individuals with ASD compared to the baseline of the general population [4]. Moreover, meta-analyses investigating the link between vitamin D and ASD have reported consistently low levels of vitamin D in those with ASD [5]. Importantly, a prospective study involving pregnant women showed that supplementation of vitamin D during pregnancy does not substantially reduce the risk of developing ASD in offspring [6]. Likewise, another meta-analysis and systematic review suggested that using vitamin D as a treatment does not significantly alleviate the core symptoms of ASD [7].Due to the prevailing uncertainty about the effects of vitamin D supplementation on patients with ASD, a thorough investigation and clearer understanding of the association between vitamin D and ASD is warranted.MR analysis [8] assesses the causality between observed changes in exposure or clinical risk factors and clinical outcomes. When random controlled trials do not provide clear causal relationships, and observational studies suffer from biases due to a multitude of potential confounders or reverse causality, genetic approaches can serve as an effective alternative for assessing causation [9]. This study aims to investigate the relationship between vitamin D and ASD. The findings suggest that vitamin D supplementation does not confer benefits to ASD patients and that ASD may, in fact, contribute to lower levels of vitamin D, providing novel insights for advancing medical treatment strategies for ASD patients.

2 Materials and methods

2.1 Data source

Relevant GWAS datasets were collected via the website (https://gwas.mrcieu.ac.uk). The vitamin D GWAS dataset (ID: ebi-a-GCST90000615) originates from a GWAS meta-analysis, which includes 417,580 participants of European ancestry [10]. The ASD GWAS dataset (ID: ieu-a-1185) comes from another study from the iPSYCH consortium [11], downloaded from the iPSYCH-PGC ASD portal (http://www.med.unc.edu/pgc), comprising 46,351 individuals of European ancestry and more than 8,643,988 single nucleotide polymorphisms (SNPs). All of the above data come from publicly available databases.

2.2 Instrumental variable selection

Using SNPs as instrumental variables, to avoid bias caused by linkage disequilibrium (LD) [12], we conducted LD clumping analysis with PLINK2, utilizing an r2 of 0.001 and a 10,000 kb window threshold on the relevant SNPs [13]. It is noteworthy that these thresholds are consistent with the widely accepted standards in MR analysis. To minimize the bias due to weak instrumental variables, we used the following formulas to calculate the R2 and F-statistic: R2 = (2 *β^2 * EAF * (1 - EAF))/(2 *β^2 * EAF * (1 - EAF) + 2 * (se^2 * N * EAF * (1 - EAF))); F = (N - K - 1) *R2/(K * (1 - R2)) [14]. In these formulas, β is the effect size of the selected SNP, EAF is the effect allele frequency of the selected SNP, SE is the standard error of the effect size of the selected SNP, R2 is the exposure level of the selected SNP, N is the total number of participants exposed, and K is the number of SNPs filtered out [15]. To ensure the robustness of our analysis, we removed any instrumental variables with an F-value lower than 10. Additionally, we excluded SNPs associated with outcomes with p < 5.0 * 10^-8 and SNPs with p < 5.0 * 10^-6 [16]. Significant SNPs associated with the exposure factor were extracted from the outcome GWAS dataset, and the final instrumental variables were documented, together with information such as effect allele, allele effect size (β), standard error (SE), and p-value.

2.3 Study design

2.3.1 Statistical methods

Data statistical analysis was performed using R version 4.3.2 and packages such as BWMR, TwoSampleMR, MR-PRESSO, Mendelian Randomization, and others. Summary statistics for the coordinated adjustment of related exposure factors and clinical outcome datasets were processed to align the effects of SNPs on the exposure factors with their effects on clinical outcomes in relation to the same alleles.

For the bidirectional Two-sample MR analysis to infer causal relationships, methods including IVW, weighted median, MR-Egger, simple mode, and weighted mode were employed, with IVW as the primary MR analytical method [17]. IVW combines the MR effect estimates of each SNP to produce an overall weighted estimate of the potential causal effect; the results of IVW are most reliable when the instrumental variables do not have horizontal pleiotropy [18]. Even when up to 50 % of the information from the exposure or outcome instrumental variables comes from genetic variants of invalid tools, the weighted median method can still achieve a consistent estimate of causal effects [19]. Potential pleiotropy associated with the selected instrumental variables was assessed using the MR-Egger intercept and MR-PRESSO global test. If pleiotropy was detected, SNPs with values less than 1.00 were excluded as outliers [20], and the Steiger test was used to determine the direction of the causal inference. For accuracy, Leave-One-Out analysis was conducted to identify outlier SNPs that may affect the precision of the results. Heterogeneity was quantified with Cochran's Q [21]test, with P < 0.05 indicating heterogeneity, though negligible heterogeneity can be ignored. If heterogeneity exists among the instrumental variables, the IVW random effects model was used to estimate the causal effect. The level of significance α was set at 0.05, with P < 0.05 considered statistically significant [22].(see Fig. 1)Fig. 1 The flowchart is designed to highlight the selection of instrumental variables, labeled "instrumental variables," which must remain unaffected by any confounders that influence the outcome (ASD). Additionally, these instrumental variables cannot have a direct effect on ASD, thereby ensuring compliance with Mendelian randomization requirements. They may only influence the outcome (ASD) through the exposure factor (Vitamin D). The aforementioned phase of result discovery is based on the verification of a non-causal relationship between Vitamin D and ASD using BWMR.

Fig. 1

The basic steps of BWMR.Eq (1) [23].Eq(1) p(γ^,Γ^,β,π1,γ,w|σX2,σY2,τ2,σ2,σ02,a0)=1AN(β|0,σ02)Pr(π1|a0)∏j=1NN(γj|0,σ2)∏j=1NPr(wj|π1)∏j=1NN(γ^j|γj,σXj2)∏j=1NN(Γ^j|βγj,σYj2+τ2)wj

include: utilizing the principles of MR to estimate the causal association between exposure and outcome using instrumental variables. By establishing a Bayesian model, causal inference is enhanced. Bayesian inference techniques, such as Markov chain Monte Carlo (MCMC) methods, are employed to obtain point estimates and confidence intervals of parameters from the posterior distribution. The advantages of BWMR methods include the ability to handle issues such as genetic correlation, measurement error, and confounding factors, while providing a probability distribution for causal inference. It can further enhance the reliability and accuracy of causal inference. Refer to https://github.com/jiazhao97/BWMR for detailed code examples and instructions on using the R package.

3 Results

3.1 Forward MR analysis (causal effect of vitamin D on ASD)

For the causal effect of vitamin D as the exposure and ASD as the outcome, with the significance threshold set at P < 5.0 × 10^-8, 79 SNPs were ultimately selected as instrumental variables, all with F-values >10. The MR-Egger regression intercept was near zero (Egger's intercept = −0.00751, P = 0.0843), indicating no horizontal pleiotropy in the instrumental variables and no significant impact on the MR analysis results. Directionality was confirmed correct by the Steiger test (P = 3.14 × 10^-109). As there was no horizontal pleiotropy in the instrumental variables, the IVW random effects model was used as the main analytical method for the MR analysis, with IVW: OR = 0.984, 95 % CI: 0.82–1.18, P = 0.866, suggesting no causal relationship between vitamin D and ASD. The BWMR results in the validation stage also indicated no causal link (OR = 0.984, 95 % CI: 0.829–1.20, P = 0.994) with Cochran's Q = 149, Cochran's Q degrees of freedom (Qdf) = 78, (P = 2.34 × 10^-6), indicating only weak heterogeneity between the instrumental variables. Sensitivity analysis was performed using leave-one-out, where each SNP was removed in turn, comparing the causal effect of the remaining SNPs with the MR analysis results for all SNPs, indicating a reliable MR analysis outcome. The detailed analysis results are shown in (Fig. 2a and b, and Table).Fig. 2 Fig. 2a Forest plots from five methods show that a p-value less than 0.05 denotes statistical significance. An odds ratio (OR) greater than 1 indicates a risk factor, while an OR less than 1 suggests a protective factor.

Fig. 2b Leave-one-out method, forest plots, and funnel plot scatter diagrams are used to assess the robustness of the results.

Fig. 2

Table 1 He results indicate that the number of SNPs (single nucleotide polymorphisms) and their directional p-values being less than 0.005 signify statistical significance. A Steiger p-value of less than 0.005 also holds statistical meaning, suggesting the correct causal direction is indicated by correct_causal_direction.

Table Results. Table. The results show that the number of SNPs and their directional p-values less than 0.005 indicate statistical significance.

Table 1TwoSampleMR	
id.exposure	id.outcome	outcome	exposure	steiger p	correct_causal_direction	nsnp	AVERAGE F	
ebi-a-GCST90000615	ieu-a-1185	ASD	VitaminD	3.13 × 10*-109	TURE	79	142.452	
id.exposure	id.outcome	outcome	exposure	steiger p	correct_causal_direction	nsnp	AVERAGE F	
ieu-a-1185	ebi-a-GCST90000615	VitaminD	ASD	5.46 × 10*-63	TURE	14	26.572	
BWMR	
id.exposure	id.outcome	outcome	exposure	steiger p	correct_causal_direction	nsnp	AVERAGE F	
ebi-a-GCST90000615	ieu-a-1185	ASD	VitaminD	3.13 × 10*-109	TURE	79	142.452	
id.exposure	id.outcome	outcome	exposure	steiger p	correct_causal_direction	nsnp	AVERAGE F	
ieu-a-1185	ebi-a-GCST90000615	VitaminD	ASD	5.46 × 10*-63	TURE	14	26.572	

3.2 Reverse MR analysis (causal effect of ASD on vitamin D)

Considering ASD as the exposure and vitamin D as the outcome, the study adopted a more liberal significance threshold (P < 5.0 × 10^-6) [23] to include more SNPs. Ultimately, 14 SNPs were selected as instrumental variables, all with F-values >10. Cochran's Q = 14.2, Q degrees of freedom (Qdf) = 13, (P = 0.359) indicated no heterogeneity among the instrumental variables. The MR-PRESSO test showed no outliers (P = 0.398), and the MR-Egger regression intercept was near zero (Egger's intercept = −0.000688, P = 0.748), which indicates no pleiotropy. The directionality was confirmed by the Steiger test (P = 5.46 × 10^-63). Therefore, the IVW random effects method was used, and the IVW results (IVW: OR = 0.976, 95 % CI: 0.96–0.99, P = 0.000855) indicated that the occurrence of ASD could lead to a reduction in vitamin D levels. The BWMR validation phase (OR = 0.975, 95 % CI: 0.958–0.991, P = 0.00297) also supported the presence of the above-mentioned causal relationship. Sensitivity analysis using leave-one-out was performed; SNPs were removed individually, and the causal effect of the remaining SNPs was compared with the MR analysis results for all SNPs, showing the MR analysis to be reliable. Detailed analysis results are provided in (Fig. 3a and b, and Table).Fig. 3 Fig. 3a Forest plots from five methods show that a p-value less than 0.05 denotes statistical significance. An odds ratio (OR) greater than 1 indicates a risk factor, while an OR less than 1 suggests a protective factor. Fig. 3b Leave-one-out method, forest plots, and funnel plot scatter diagrams are used to assess the robustness of the results.

Fig. 3

4 Discussion

ASD is a complex neurodevelopmental disorder with a multifaceted etiology, involving various genetic and environmental factors [24] Although current research has not yet determined the specific mechanisms that cause ASD, our study has meticulously confirmed that there is no direct link between vitamin D and ASD. However, in a reverse analysis where ASD leads to a decrease in vitamin D levels, our findings lend some support to the rationale for choosing vitamin D as an intervention for ASD. This significantly bolsters our confidence in our procedural design and conclusions.

From the molecular perspective, ASD may be a highly complex brain neurosystem disorder where the causes may involve multiple interactions between genetic regulation and environmental factors [25]. Vitamin D is widely considered to play a crucial role in the development and maintenance of the nervous system. Its main roles may include regulation of gene expression, immune response, inflammation, and its impact on the development of the nervous system (vitamin D's involvement in the brain's neural construction may affect early-stage ASD) [[26], [27], [28]].However, to date, the role and regulatory mechanisms between vitamin D and ASD remain contested.

Our focus will be on the potential impacts of vitamin D on neurodevelopment, immune responses, and inflammatory processes, as well as gene expression regulation, aiming to delve deeply into whether there are substantive connections between these biological processes.

First, the limitations of neuroprotection: Although vitamin D has some neuroprotective effects [29], what might be needed for treating ASD patients is a more complex and precise neural system regulation, including neurotransmitter levels, synaptic function, and neural circuitry shaping [30], where potential neural system regulation likely does not wholly depend on the biological role played by vitamin D [31]. Secondly, the dosage and timing: Vitamin D is indeed more capable of crossing the blood-brain barrier (BBB) compared to water-soluble substances [32]. This is because the BBB has a certain permeability to lipid-soluble molecules, which can pass through the lipid bilayer of cell membranes. Vitamin D [33], as a fat-soluble hormone, can pass into the brain relatively freely through the BBB. However, even if vitamin D can relatively easily cross the BBB, this does not mean that it can freely participate in regulatory mechanisms or exert biological effects in the brain [34]. Overall levels of vitamin D in the body, the state of the BBB, the residual amount of vitamin D after crossing the BBB, and its metabolism and utilization in the brain will all affect whether vitamin D can exert real effects on brain cells [[35], [36], [37]].

The limitations of immunity and inflammation: Firstly, the impact of neuroinflammatory pathways [38]: Some studies have found signs of neuroinflammation in ASD patients [39]. Although vitamin D can reduce inflammatory responses by lowering the production of pro-inflammatory cytokines [40], the root cause of neuroinflammation in ASD may not be merely due to immune system abnormalities but could involve other contributory factors [41,42]. Secondly, the contribution of glial cells: Functions of microglia and astrocytes, which regulate the integration of brain neuronal signals and participate in maintaining the neuronal environment in the brain, may be impaired in ASD [43]. However, these changes are not merely related to the singular immune response involving vitamin D but may also involve broader regulation of the brain's internal environment homeostasis [[44], [45], [46]].

The complexity of gene expression and regulation: Firstly, vitamin D can affect gene expression through its receptor (Vitamin D receptor, VDR) [47] within the cell nucleus. However, ASD may involve a vast number of genes with extremely complex mechanisms of interaction, and a singular vitamin D intervention is unlikely to broadly change this complex gene regulatory network. Secondly, the temporality of neurodevelopment: The neurodevelopmental abnormalities in ASD likely take place early in life, which is a precisely regulated process involving temporal gene expression changes [48,49]. Vitamin D's gene regulatory behavior may not precisely match the specific timing requirements associated with neurodevelopment [50]. Lastly, it is worth noting the limitations of epigenetic regulation: Besides changes in the genes themselves, the development of ASD may also relate to epigenetic factors, such as DNA methylation and histone modification. These modifications do not alter the DNA sequence but affect the regulation of gene expression, and vitamin D's role is likely limited here [[51], [52], [53]].

In our reverse analysis, we accidently discovered that the occurrence of ASD might lead to a decrease in vitamin D levels in the body, a finding that caught our high attention. Therefore, we explored its causes from both biological and environmental perspectives. First, at the genetic level, ASD may influence the gene regulation of the VDR and associated enzymes, such as the 1-α-hydroxylase which is responsible for vitamin D activation. Variations or abnormal expressions of these genes might interfere with the normal metabolism of vitamin D [54], leading to a decreased level in the bloodstream.Next is the regulation of Vitamin D Binding Protein (VDBP) [55]. Since VDBP is a crucial protein that transports vitamin D, if VDBP expression or function is affected in patients with ASD, the level of available vitamin D in the blood may be reduced correspondingly.

Furthermore, the inflammatory pathway is also a factor to consider. Studies have shown that ASD may be related to the body's state of inflammation [56]. Considering that vitamin D plays an essential role in modulating the immune response, it may be consumed during the inflammatory process, which might mean that higher levels of inflammation could reduce the levels of vitamin D in the blood [41].

Lastly, environmental factors were considered. For example, children and adults with ASD might spend less time engaging in outdoor activities [57]. Since vitamin D is primarily produced in the skin through exposure to sunlight [58], the reduction in sunlight exposure can directly lead to decreased production of vitamin D in the body [59]. Additionally, problematic behaviors based on restrictive and repetitive actions, such as particular dietary preferences and aversions, could lead to an imbalance in nutritional intake, including insufficient vitamin D [60,61].

Our research methodology offers several distinctive advantages. This study represents the inaugural application of the bidirectional MR framework for investigative purposes, with the BWMR technique providing subsequent validation. The latter is acclaimed for its exceptional adaptability in addressing numerous interrelated genetic variations, outcomes, or the potency of genetic instrumental variables [62]. In instances of genetic pleiotropy or when the instrumental variables demonstrate weak associations, BWMR delivers estimations of enhanced robustness by amalgamating diverse information sources and types [63]. BWMR also effectively ameliorates the minor heterogeneity detected in our investigation. The Bayesian framework furnishes a more solid evaluation in the face of result heterogeneity [63]. Moreover, the utilization of GWAS data sourced from an extensive sample size and recent collection substantially corroborates the credibility of our findings.

Our investigation reveals certain limitations. The GWAS dataset employed originates from a European cohort and may reflect ethnic biases, rendering the extrapolation of the positive findings to other ethnicities uncertain. Furthermore, the precision of the instrumental variables is contingent upon the sample size of the GWAS; consequently, an expanded GWAS is imperative to uncover a broader array of genetic variations pertinent to MR analysis. Collectively, our findings indicate no causal link between vitamin D levels and the incidence of ASD. It may be essential to consider reverse causality; for instance, the development of ASD could precipitate diminished vitamin D concentrations. These insights confer an elevated tier of evidence in clinical research, markedly diminishing expenditures and human resource investments. They also introduce an innovative angle to ASD research, inform the development of therapeutic strategies by clinicians for ASD, and alleviate the necessity for redundant trials by clinical investigators.

5 Conclusion

In summary, our research found that the levels of Vitamin D did not have a causal intervention effect on autism, but surprisingly, the occurrence of autism actually led to a decrease in the levels of Vitamin D in the body.

Ethics approval and consent to participate

I would like to confirm that the data used in this study is publicly available and can be reused without restrictions under an open license. All datasets in this research have been downloaded from public databases. These public databases permit researchers to download and analyze public datasets for scientific purposes, thereby negating the need for ethical approval.

Code availability

The codes used in this study can be found at:

TwoSampleMR:https://github.com/MRCIEU/TwoSampleMR.

MR-PRESSO:https://github.com/rondolab/MR-PRESSO.

BWMR:https://github.com/jiazhao97/BWMR.

Funding

Jiamusi University Doctoral Startup Project (JMSUBZ2020-07); Jiamusi University 2023 Young Innovative Talent Development Support Program Project (JMSUQP23039); Heilongjiang Youth Science and Technology Talent Support Program Project (2022QNTJ016); Heilongjiang Province Provincial Higher Education Institutions Basic Research Operating Expenses Team Project(2022-KYYWF-0653 );10.13039/501100014902 Jiamusi University "East Pole" Academic Children's Intelligent Rehabilitation Team(DJXSTD202413 ).

CRediT authorship contribution statement

Tianci Gao: Writing – review & editing, Writing – original draft, Visualization, Validation, Supervision, Software, Resources, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization. Wenjun Dang: Writing – review & editing, Writing – original draft, Validation, Methodology, Funding acquisition, Formal analysis, Conceptualization. Zhimei Jiang: Visualization, Supervision, Software, Project administration, Funding acquisition, Formal analysis. Yuwei Jiang: Writing – original draft, Visualization, Software, Resources, Project administration, Methodology, Investigation, Funding acquisition, Formal analysis, Data curation, Conceptualization.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
==== Refs
References

1 Liang X. Age-related differences in accelerometer-assessed physical activity and sleep parameters among children and adolescents with and without autism spectrum disorder: a meta-analysis JAMA Netw. Open 6 10 2023 e2336129
2 Oliver L.D. Social cognitive performance in schizophrenia spectrum disorders compared with autism spectrum disorder: a systematic review, meta-analysis, and meta-regression JAMA Psychiatr. 78 3 2021 281 292
3 Bakthavachalu P. Kannan S.M. Qoronfleh M.W. Food color and autism: a meta-analysis Adv Neurobiol 24 2020 481 504 32006369
4 Tsujiguchi H. Prospective relationship between autistic traits and nutrient intakes among Japanese children: results of the Shika study Autism 27 2 2023 389 401 35652325
5 Jiang Y. Omega-3 polyunsaturated fatty acids and/or vitamin D in autism spectrum disorders: a systematic review Front Psychiatry 14 2023 1238973
6 Madley-Dowd P. Maternal vitamin D during pregnancy and offspring autism and autism-associated traits: a prospective cohort study Mol. Autism. 13 1 2022 44 36371219
7 Li B. The effect of vitamin D supplementation in treatment of children with autism spectrum disorder: a systematic review and meta-analysis of randomized controlled trials Nutr. Neurosci. 25 4 2022 835 845 32893747
8 Cui R. Correspondence on 'Systemic evaluation of the relationship between psoriasis, psoriatic arthritis and osteoporosis: observational and Mendelian randomisation study' Ann. Rheum. Dis. 82 1 2023 e13 33144297
9 Liu K.Q. Response to: 'Correspondence on 'Systemic evaluation of the relationship between psoriasis, psoriatic arthritis and osteoporosis: observational and Mendelian randomisation study'' by Cui et al Ann. Rheum. Dis. 82 1 2023 e14 33144296
10 Revez J.A. Genome-wide association study identifies 143 loci associated with 25 hydroxyvitamin D concentration Nat. Commun. 11 1 2020 1647 32242144
11 Pedersen E.M. ADuLT: an efficient and robust time-to-event GWAS Nat. Commun. 14 1 2023 5553 37689771
12 Forstner A.J. Genome-wide association study of panic disorder reveals genetic overlap with neuroticism and depression Mol Psychiatry 26 8 2021 4179 4190 31712720
13 Brito Nunes C. Mendelian randomization study of maternal coffee consumption and its influence on birthweight, stillbirth, miscarriage, gestational age and pre-term birth Int. J. Epidemiol. 52 1 2023 165 177 35679582
14 Ren Q. Major depression disorder may causally associate with the increased breast cancer risk: evidence from two-sample mendelian randomization analyses Cancer Med. 12 2 2023 1984 1996 35852181
15 Spiga F. Tools for assessing quality and risk of bias in Mendelian randomization studies: a systematic review Int. J. Epidemiol. 52 1 2023 227 249 35900265
16 Geurts S. Heart rate variability and atrial fibrillation in the general population: a longitudinal and Mendelian randomization study Clin. Res. Cardiol. 112 6 2023 747 758 35962833
17 Huang J. Association between cannabis use with urological cancers: a population-based cohort study and a mendelian randomization study in the UK biobank Cancer Med. 12 3 2023 3468 3476 35975633
18 Baranova A. Cao H. Zhang F. Causal effect of COVID-19 on Alzheimer's disease: a Mendelian randomization study J. Med. Virol. 95 1 2023 e28107
19 Yuan S. Morning chronotype and digestive tract cancers: mendelian randomization study Int. J. Cancer 152 4 2023 697 704 36093575
20 Rukh G. de Ruijter M. Schiöth H.B. Effect of worry, depression, and sensitivity to environmental stress owing to neurotic personality on risk of cardiovascular disease: a Mendelian randomization study J. Pers. 91 3 2023 856 867 36165189
21 Liu N. Non-alcoholic fatty liver disease and complications in type 1 and type 2 diabetes: a Mendelian randomization study Diabetes Obes Metab 25 2 2023 365 376 36181433
22 Decina C.S. Investigating a possible causal relationship between maternal serum urate concentrations and offspring birthweight: a Mendelian randomization study Int. J. Epidemiol. 52 1 2023 178 189 36191079
23 Zhang J. Causal relationship between coffee intake and neurological diseases: a Mendelian randomization study Eur. J. Clin. Nutr. 78 2 2024 Feb 114 119 37845420
24 Russo A.M. Social approach, anxiety, and altered tryptophan hydroxylase 2 activity in juvenile BALB/c and C57BL/6J mice Behav. Brain Res. 359 2019 918 926 29935278
25 Cattane N. Richetto J. Cattaneo A. Prenatal exposure to environmental insults and enhanced risk of developing Schizophrenia and Autism Spectrum Disorder: focus on biological pathways and epigenetic mechanisms Neurosci. Biobehav. Rev. 117 2020 253 278 29981347
26 Filgueiras M.S. Vitamin D status, oxidative stress, and inflammation in children and adolescents: a systematic review Crit. Rev. Food Sci. Nutr. 60 4 2020 660 669 30596263
27 Khairy E.Y. Attia M.M. Protective effects of vitamin D on neurophysiologic alterations in brain aging: role of brain-derived neurotrophic factor (BDNF) Nutr. Neurosci. 24 8 2021 650 659 31524100
28 Liu H. Molecular regulatory mechanism of LILRB4 in the immune response Cent. Eur. J. Immunol. 48 1 2023 43 47 37206591
29 Iacopetta K. Are the protective benefits of vitamin D in neurodegenerative disease dependent on route of administration? A systematic review Nutr. Neurosci. 23 4 2020 251 280 29985117
30 AlJohri R. AlOkail M. Haq S.H. Neuroprotective role of vitamin D in primary neuronal cortical culture eNeurologicalSci 14 2019 43 48 30619951
31 Trojsi F. Vitamin D supplementation has no effects on progression of motor dysfunction in amyotrophic lateral sclerosis (ALS) Eur. J. Clin. Nutr. 74 1 2020 167 175 31197218
32 Riccio P. Rossano R. Diet, gut microbiota, and vitamins D + A in multiple sclerosis Neurotherapeutics 15 1 2018 75 91 29067566
33 Sayeed I. Vitamin D deficiency increases blood-brain barrier dysfunction after ischemic stroke in male rats Exp. Neurol. 312 2019 63 71 30502340
34 Alam C. Upregulation of reduced folate carrier by vitamin D enhances brain folate uptake in mice lacking folate receptor alpha Proc Natl Acad Sci U S A 116 35 2019 17531 17540 31405972
35 Uher T. Neurofilament levels are associated with blood-brain barrier integrity, lymphocyte extravasation, and risk factors following the first demyelinating event in multiple sclerosis Mult. Scler. 27 2 2021 220 231 32255388
36 Yang J. Vitamin D3 supplement attenuates blood-brain barrier disruption and cognitive impairments in a rat model of traumatic brain injury NeuroMolecular Med. 23 4 2021 491 499 33616826
37 Zelzer S. A novel method for the determination of vitamin D metabolites assessed at the blood-cerebrospinal fluid barrier Biomolecules 11 9 2021
38 Alrefaie Z. Uncoupling proteins: are they involved in vitamin D3 protective effect against high-fat diet-induced cardiac apoptosis in rats? Arch. Physiol. Biochem. 128 2 2022 438 446 31794287
39 Liang Y. Vitamin D deficiency worsens maternal diabetes induced neurodevelopmental disorder by potentiating hyperglycemia-mediated epigenetic changes Ann. N. Y. Acad. Sci. 1491 1 2021 74 88 33305416
40 Esnafoglu E. Subaşı B. Association of low 25-OH-vitamin D levels and peripheral inflammatory markers in patients with autism spectrum disorder: vitamin D and inflammation in Autism Psychiatry Res 316 2022 114735
41 Harutyunyan A.A. Harutyunyan H.A. Yenkoyan K.B. Novel probable glance at inflammatory scenario development in autistic pathology Front Psychiatry 12 2021 788779
42 Tavasoli S. Taheri M. Vitamin D and calcium kidney stones: a review and a proposal Int. Urol. Nephrol. 51 1 2019 101 111 30136085
43 Lee K.M. The application of human pluripotent stem cells to model the neuronal and glial components of neurodevelopmental disorders Mol Psychiatry 25 2 2020 368 378 31455859
44 Tomova A. Specificity of gut microbiota in children with autism spectrum disorder in Slovakia and its correlation with astrocytes activity marker and specific behavioural patterns Physiol. Behav. 214 2020 112745
45 Sciara A.N. Neuroinflammatory gene expression alterations in anterior cingulate cortical white and gray matter of males with autism spectrum disorder Autism Res. 13 6 2020 870 884 32129578
46 Liao X. Postmortem studies of neuroinflammation in autism spectrum disorder: a systematic review Mol. Neurobiol. 57 8 2020 3424 3438 32529489
47 Ustun N. Association of vitamin D receptor gene FokI and TaqI polymorphisms and risk of RDS J. Matern. Fetal Neonatal Med. 33 21 2020 3640 3646 30760068
48 Ali A. Cui X. Eyles D. Developmental vitamin D deficiency and autism: putative pathogenic mechanisms J. Steroid Biochem. Mol. Biol. 175 2018 108 118 28027915
49 Ganguly A. Vitamin D, the placenta and early pregnancy: effects on trophoblast function J. Endocrinol. 236 2 2018 R93 r103 29109081
50 Lisi G. Maternal vitamin D and its role in determining fetal origins of mental health Curr Pharm Des 26 21 2020 2497 2509 32370709
51 Noroozi R. Identification of miRNA-mRNA network in autism spectrum disorder using a bioinformatics method J. Mol. Neurosci. 71 4 2021 761 766 32875540
52 Boulting G.L. Activity-dependent regulome of human GABAergic neurons reveals new patterns of gene regulation and neurological disease heritability Nat. Neurosci. 24 3 2021 437 448 33542524
53 Abdul F. Disruption of circadian rhythm and risk of autism spectrum disorder: role of immune-inflammatory, oxidative stress, metabolic and neurotransmitter pathways Rev. Neurosci. 33 1 2022 93 109 34047147
54 Ogbu D. Xia E. Sun J. Gut instincts: vitamin D/vitamin D receptor and microbiome in neurodevelopment disorders Open Biol 10 7 2020 200063
55 Greilberger J. Herwig R. Vitamin D - deglycosylated vitamin D binding protein dimer: positive synergistic effects on recognition, activation, phagocytosis and oxidative stress on macrophages Clin. Lab. 66 1 2020
56 Marballi K. MacDonald J.L. Proteomic and transcriptional changes associated with MeCP2 dysfunction reveal nodes for therapeutic intervention in Rett syndrome Neurochem. Int. 148 2021 105076
57 Mpoulimari I. Zintzaras E. Synthesis of genetic association studies on autism spectrum disorders using a genetic model-free approach Psychiatr. Genet. 32 3 2022 91 104 35353796
58 Ghafouri-Fard S. Expression analysis of VDR-related LncRNAs in autism spectrum disorder J. Mol. Neurosci. 71 7 2021 1403 1409 34125396
59 Grübler M.R. Effects of Vitamin D3 on asymmetric- and symmetric dimethylarginine in arterial hypertension J. Steroid Biochem. Mol. Biol. 175 2018 157 163 28027911
60 Muskens J. Vitamin D status in children with a psychiatric diagnosis, autism spectrum disorders, or internalizing disorders Front Psychiatry 13 2022 958556
61 Bolognesi E. GC1f vitamin D binding protein isoform as a marker of severity in autism spectrum disorders Nutrients 14 23 2022
62 Liu A. Prioritization of risk genes in multiple sclerosis by a refined Bayesian framework followed by tissue-specificity and cell type feature assessment BMC Genom. 23 Suppl 4 2022 362
63 Deng Y.T. Identifying causal genes for depression via integration of the proteome and transcriptome from brain and blood Mol Psychiatry 27 6 2022 2849 2857 35296807
