Vitamin D is a steroid hormone that regulates the metabolism of Ca, P and bone and promotes gene expression through endocrine and autocrine mechanisms(Reference Heaney1). 25-Hydroxyvitamin D is produced by the metabolism of vitamin D in the liver, and it is the main circulating form of vitamin D. Therefore, 25-hydroxyvitamin D can be used as a clinical indicator of vitamin D levels(Reference Bouillon, Antonio and Olarte2). There is growing evidence that vitamin D may play a role in the development of a number of diseases, including cancer, diabetes, autoimmune diseases and infectious diseases(Reference Manson, Cook and Lee3–Reference Sîrbe, Rednic and Grama5). The results of these studies indicated that vitamin D plays a significant role in Ca homeostasis, immunomodulation and anti-inflammation.
Some publicly available studies have reported an association between vitamin D and the development of autoimmune thyroid disease, but the results have been inconclusive(Reference Kivity, Agmon-Levin and Zisappl6–Reference Hahn, Cook and Alexander10). For example, a study conducted by Kivity et al. revealed that the prevalence of vitamin D deficiency was considerably greater among 50 patients with autoimmune thyroid disease than among 98 healthy individuals (72 % v. 30·6 %; P < 0·001)(Reference Kivity, Agmon-Levin and Zisappl6). In a recent 5-year randomised controlled trial conducted in the USA, 25 871 participants were followed, and a 22 % reduction in autoimmune disease incidence was observed after vitamin D supplementation after 5 years(Reference Hahn, Cook and Alexander10). The association between vitamin D and autoimmune diseases may be related to variations in the vitamin D receptor gene, which is involved in vitamin D function(Reference Feng, Li and Chen11,Reference Inoue, Watanabe and Ishido12) . Nevertheless, a study conducted by D’Aurizio et al. revealed that there was no statistically significant difference in vitamin D levels between patients with autoimmune thyroid disease and the general population(Reference D’Aurizio, Villalta and Metus13). Furthermore, vitamin D has been reported to reduce the incidence of thyroid cancer. Vitamin D can play a role through the following mechanisms: increasing apoptosis, arresting the cell cycle, inhibiting proliferation and differentiation, decreasing the inflammatory response and decreasing aggressiveness(Reference Hansen, Binderup and Hamberg14–Reference Liu, Zhang and Xu16). The parathyroid gland is a direct target of vitamin D. Parathyroid cells express both the vitamin D receptor and 1-α-hydroxylase. The existing literature indicates that vitamin D metabolites exert an influence on parathyroid hormone secretion and might act to prevent parathyroid cell proliferation(Reference Bienaimé, Prié and Friedlander17). The majority of published studies have focused on the impact of vitamin D on autoimmune thyroid disease and thyroid cancer. Few studies have explored the effects of vitamin D on other thyroid- and parathyroid-related diseases.
Data on vitamin D and thyroid- and parathyroid-related diseases derived from observational studies will inevitably be affected by factors such as sample size, ethnicity and other confounding variables, making causal inference challenging. Mendelian randomisation (MR) is a highly effective methodology for investigating the causal relationship between exposure and disease since disease status typically does not alter the germline DNA sequences(Reference Nattel18). This is accomplished by utilising genetic variation as an instrumental variable (IV)(Reference O’Donnell and Sabatine19). Genetic variants are randomly allocated during meiosis, much like a random assignment in a randomised controlled trial; additionally, the genetic variants undergo minimal changes throughout an individual’s lifetime(Reference Zheng, Baird and Borges20), which minimises unmeasured confounding factors and biases caused by reverse causation(Reference Sanderson, Glymour and Holmes21). Previous MR studies have demonstrated that there is no causal relationship between serum vitamin D levels and the development of Graves’ disease or thyroid cancer(Reference Yu, Yang and Wu22,Reference Shen, Zhang and Jiang23) . The causal relationships between vitamin D and other thyroid- and parathyroid-related diseases are unclear. It is not clear what role vitamin D plays in the development of thyroid- and parathyroid-related diseases, whether it is a cause of the disease, a side effect of treatment or a consequence of the disease(Reference Vieira, Rodrigues and Paiva24). The MR studies can clarify the relationship between vitamin D and thyroid- and parathyroid-related diseases from a gene- and bidirectional perspective, which can facilitate the optimisation of disease prevention and management strategies.
Therefore, we performed a two-sample bidirectional MR to comprehensively evaluate the causal relationships between vitamin D levels and diseases related to the thyroid and parathyroid glands using an extensive genome-wide association study (GWAS, which is a strategy that is used in the analysis of complex traits. GWAS involves the scanning of the genome for millions of SNP molecular markers, with the aim of identifying genotypic and phenotypic correlations that affect these traits(Reference Flint25)).
Methods
Data availability and ethics statement
This study involved a secondary examination of publicly available data, and it did not involve new human or animal research. All utilised GWAS datasets that were used are openly accessible, negating the need for ethical approval or informed consent. This study followed the Strengthening the Reporting of Observational Studies in Epidemiology Using Mendelian Randomisation reporting guidelines(Reference Skrivankova, Richmond and Woolf26).
Study design
To examine the possible causal relationship between vitamin D levels and thyroid- and parathyroid-related diseases, we conducted a two-way, bidirectional MR analysis. The use of MR necessitates satisfying three fundamental assumptions(Reference König and Greco27): (i) association, where the chosen IV must be highly connected to the exposure; (ii) exclusivity, where the IV can solely influence the outcome via the exposure; and (iii) independence, where confounding factors cannot impact the effect of exposure on the outcome. Fulfilling these three criteria allows us to estimate the exposure-outcome relationship, utilising the obtained IV(Reference Zheng, Baird and Borges20). The framework is presented in Fig. 1.
Exposure and outcome data sources
Summary statistics on vitamin D levels were extracted from the UK Biobank(Reference Collins28), comprising phenotypic, genotypic and clinical data for 417 580 individuals of European descent (age range 40–69 years), with ID number ‘ebi-a-GCST90000614’. Since vitamin D from any source is quickly transformed into 25(OH)D, which has a long half-life of approximately 2 weeks and is a crucial precursor of active hormones, it is generally accepted that the serum 25(OH)D levels represent the levels of vitamin D that are stored in the body and the serum 25(OH)D levels are the optimal indicator for estimating vitamin D nutritional intake(Reference Bouillon, Antonio and Olarte2). In this study, the serum 25(OH)D concentration (nmol/l) was measured using the LIAISON XL 25(OH)D assay (DiaSorin). The median, mean and interquartile range of serum 25(OH)D levels were 47·9, 49·6 and 33·5–63·2 nmol/l, respectively. Serum 25(OH)D concentrations below 25 nmol/l were considered to indicate vitamin D deficiency(Reference Revez, Lin and Qiao29).
To eliminate potential bias caused by overlapping samples in terms of both exposure and outcomes, we collected summary statistics on conditions involving the thyroid (such as autoimmune diseases, thyroid cancer, thyroiditis and thyroid nodules) and parathyroid glands from the FinnGen Biobank. FinnGen Biobank conducted a genomics and personalised medicine research initiative that involved by examining genomic and health data taken from 500 000 Finnish biobanks to understand the genetic basis of disease(30). Specific GWAS summary data are provided in Table 1. To prevent differences in ethnicity from impacting pleiotropy, the biobank included only study participants of European origin. Using the summary data from GWAS, an independent study employing two-sample bidirectional MR was conducted to investigate the relationship between serum vitamin D levels and thyroid- and parathyroid-related diseases.
Mendelian randomisation analysis
Selection of instrumental variables
Instrumental variables were determined through various processes. First, to meet the association assumption, a threshold of P < 5 × 10–8 was established, and the threshold was relaxed to 5 × 10–6 or 1 × 10–5 when adequate SNP, which primarily DNA sequence polymorphisms that result from variation in a single nucleotide at the genomic level(31), were not available for analysis at the P < 5 × 10–8 threshold. Second, to minimise the occurrence of multiple results due to linkage disequilibrium, linkage disequilibrium analyses were conducted (r 2 < 0·001, kb = 10 000) in accordance with the required independence assumptions. The linkage disequilibrium levels were then obtained from the European samples of the 1000 Genomes Project(Reference Abecasis, Altshuler and Auton32). To ensure that the impact of selected SNP on both exposure and outcome aligned with the corresponding alleles (an allele is one of two or more versions of the DNA sequence at a given genomic location(33)), we eliminated the palindromic structure and utilised surrogate SNP when relevant SNP were absent from the GWAS dataset for outcome. Confounders not associated with SNP (including other types of thyroid disease, other vitamin levels, smoking, alcohol consumption, obesity, etc.) were manually removed by PhenoScanner(Reference Kamat, Blackshaw and Young34). The direction of causality between exposure and outcome was evaluated using Steiger filtering. If the IV met the criteria, the instrument’s direction was ‘TRUE’, and SNP with a direction of ‘FALSE’ were excluded(Reference Hemani, Tilling and Davey Smith35). To assess the strength of the IV, we calculated the F-statistic (F reflects the bias of an IV or a set of IV(Reference Pierce, Ahsan and Vanderweele36)), using the formula(Reference Burgess and Thompson37) $\left( {{{n - k - 1} \over k}} \right)\left( {{{{R^2}} \over {1 - {R^2}}}} \right)$ , with R 2 representing the proportion of variance explained by the IV, n representing sample size and k representing the number of SNP(Reference Papadimitriou, Dimou and Tsilidis38). An F-statistic greater than 10 indicates that the IV may be resistant to the effects of weak instrumental bias. This analysis adheres to the conventional academic structure and employs clear, objective language with precise technical terms(Reference Burgess and Thompson39).
Statistical analysis
All the MR analyses were conducted using R version 4.3.1 (The R Foundation for Statistical Computing). The estimation of causal relationships was conducted using the ‘TwoSampleMR’ (version 0.5.7), ‘MR-PRESSO’ (version 1.0), ‘ggplot2’ (version 3.4.3), ‘plyr’ (version 1.8.8) and ‘phenoscanner’ (version 1.0) packages. The MR methods applied included the inverse variance weighted (IVW), weighted median and MR-Egger methods. The IVW method assumed that all IV met the validity criteria to obtain unbiased estimates(Reference Burgess, Butterworth and Thompson40); specifically, it assumed that all SNP were not related to the pleiotropic effect of exposure (known as the InSIDE hypothesis)(Reference Burgess and Thompson41), which had the highest test efficacy. When the weighted median method assumed that only more than half of the IV were unbiased(Reference Bowden, Davey Smith and Haycock42), the IVW method was the primary reference for obtaining results. The weighted median and MR-Egger methods were used to complement the IVW analysis. Causal effects were assessed as OR, indicating an increased risk of outcome for each increase in the exposure log ratio. Our sensitivity analyses used four main methods, including Cochran’s Q test, MR-Egger intercept analysis, MR-PRESSO and the leave-one-out sensitivity test. Horizontal pleiotropy was evaluated through the MR-Egger test’s intercept, with a significance threshold of P < 0·05 denoting the existence of horizontal pleiotropy and the intercept term’s value distance from 0 indicating the magnitude of horizontal pleiotropy(Reference Bowden, Davey Smith and Burgess43). Heterogeneity was evaluated using the Cochran Q statistic, whereby a P value of less than 0·05 indicated the existence of heterogeneity(Reference Hemani, Zheng and Elsworth44). When heterogeneity was detected, we utilised MR-PRESSO to identify potential outliers and then removed them before re-evaluating causality with the remaining SNP(Reference Ong and MacGregor45). The ‘leave-one-out’ test was utilised to evaluate the impact of a single SNP on the analysis, and this test enhanced the robustness of the outcomes(Reference Zheng, Baird and Borges20). We performed RadialMR analysis using modified second-order weights to identify outliers(Reference Bowden, Spiller and Del Greco46). Radial plots offered advantages in identifying peripheral studies, detecting small study biases and identifying outliers more directly than traditional scatter plots.
The outcomes of the MR analyses are presented as scatter plots, forest plots, ‘leave-one-out’ plots and funnel plots. Each point in the scatter plots represents an SNP, thereby demonstrating the association of that SNP with exposure and outcome. The forest plot comprised horizontal lines and points, with each line representing the effect size of an SNP and its 95 % CI. Given the lack of robustness of the results for individual SNP, it was necessary to combine them, and this is represented by the bottom red line (All − IVW). The leave-one-out forest plot was used to calculate the meta-effect of the remaining SNP after removing each SNP one by one. If all the error lines were consistent to the right or left of 0, the results were deemed reliable. The funnel plot was generated to determine whether the points situated on either side of the IVW line were approximately symmetrical. The presence of any outlying points indicated the potential for outliers, which could be removed, and the analysis process was repeated.
To consider multiple testing, we used a conservative approach and applied a Bonferroni-corrected significance level of 0·05 divided by 13 (1 exposure × 13 outcomes, i.e. 0·0038)(Reference Curtin and Schulz47). P < 0·0038 was considered to indicate a significant association. A P < 0·05 but above the Bonferroni-corrected significance threshold was considered to indicate a potential association.
Results
Causal association of serum vitamin D levels on thyroid- and parathyroid-related diseases
Mendelian randomisation analysis of serum vitamin D levels on thyroid- and parathyroid-related diseases
Forty-eight vitamin D alleles were linked to benign parathyroid adenoma. The F-statistic for assessing the relationship of serum vitamin D levels to benign parathyroid adenoma ranged from 25 to 1220, indicating the absence of weak IV (online Supplementary Table S1). After applying the Bonferroni correction, MR estimation through the IVW method revealed that serum vitamin D levels was suggestively decreased the risk of benign parathyroid adenoma (OR = 0·244; 95 % CI 0·074, 0·802; P = 0·0202). The MR-Egger and weighted median OR were also less than 1, and the P value was less than 0·05, further validating our findings (Fig. 2).
Based on the IVW method, the MR estimation did not reveal a causal association of vitamin D with other thyroid- and parathyroid-related diseases (all P values were greater than 0·05) (online Supplementary Table S2). The F-statistics of the serum vitamin D levels for various thyroid- and parathyroid-related diseases were all greater than 10.
Sensitivity analysis and visualisation of results of serum vitamin D levels on thyroid- and parathyroid-related diseases
The selected SNP indicated no heterogeneity or pleiotropy (Cochran’s Q P value, MR-Egger intercept P value and MR-PRESSO test P value were all greater than 0·05) (Table 2). The leave-one-out analysis indicated that none of the SNP exerted a dominant influence on the effect of serum vitamin D levels on benign parathyroid adenomas (Fig. 3(c)). RadialMR detected two outliers (online Supplementary Fig. S1) After excluding the outliers that were detected by the RadialMR analysis, the MR estimates did not change significantly (OR = 0·254; 95 % CI 0·076, 0·847; P = 0·026). The funnel plot displayed an even distribution of points, indicating that causal associations were less susceptible to potential bias (Fig. 3(d)).
Causal association of thyroid- and parathyroid-related diseases on serum vitamin D levels
Mendelian randomisation analysis of thyroid- and parathyroid-related diseases on serum vitamin D levels
Five thyroid cancer alleles, three chronic thyroiditis alleles and nineteen thyroid nodule alleles were associated with serum vitamin D levels. The F-statistics were greater than 10, indicating that there were strong associations of thyroid cancer, chronic thyroiditis and thyroid nodules with serum levels of vitamin D. Additionally, the results indicated a low probability of being influenced by weak IV (online Supplementary Table S3–5) After applying the Bonferroni correction, MR estimation using the IVW approach suggested genetically predicted risk of thyroid cancer was suggestively increased the risk of elevated vitamin D (OR = 1·007; 95 % CI 1·010, 1·013; P = 0·0284), and chronic thyroiditis was significantly increased the risk of elevated vitamin D (OR = 1·007; 95 % CI 1·002, 1·011; P = 0. 0030); thyroid nodules were significantly decreased the vitamin D levels (OR = 0·991; 95 % CI 0·985, 0·997; P = 0·0034). Both the MR-Egger and weighted median OR and P values were consistent with the results of IVW results, further confirming our findings (Fig. 4).
Based on the IVW method, the MR estimation did not reveal a causal effect of other thyroid and parathyroid diseases on serum vitamin D levels (all P values greater than 0·05) (online Supplementary Table S2). The F-statistics of various thyroid- and parathyroid-related diseases were all greater than 10.
Sensitivity analysis and visualisation of results of thyroid- and parathyroid-related diseases on serum vitamin D levels
No heterogeneity or pleiotropy was observed among the chosen SNP (Cochran’s Q P values, MR-Egger intercept P values and MR-PRESSO test P values, except for the SNP linked to chronic thyroiditis, which was inadequate for MR-PRESSO) (Table 2). Furthermore, the leave-one-out method indicated that none of the SNP played a dominant role in the causal associations of thyroid cancer, chronic thyroiditis or thyroid nodules on serum vitamin D levels (Fig. 5–7(c)). RadialMR did not detect an outlier (online Supplementary Fig. S2–4). The funnel plot displayed an even distribution of points, indicating that causal associations are less susceptible to potential bias (Fig. 5–7(d)).
Discussion
To the best of our knowledge, this is the first comprehensive study using a bidirectional MR approach to investigate the relationship between serum vitamin D levels and thyroid- and parathyroid-related diseases. Our genetic analysis revealed an inverse suggestive association of vitamin D levels on parathyroid tumour development. However, our reverse MR analysis revealed a genetically suggestive causal association of thyroid cancer on serum vitamin D levels and significant causal associations of chronic thyroiditis and thyroid nodules on serum vitamin D levels. Additionally, we did not observe the causal effect between serum vitamin D levels and other thyroid and parathyroid diseases. These findings highlighted the role of vitamin D levels in potential disease intervention, prognosis and monitoring goals.
Our research indicated that vitamin D has a potential role in benign parathyroid adenoma development. Active vitamin D has been demonstrated to hinder parathyroid cell growth both in vitro and in vivo (Reference Nygren, Larsson and Johansson48,Reference Cantley, Russell and Lettieri49) . Furthermore, it can impede cell cycle progression and contribute to multiple cellular pathways that are involved in tumour formation(Reference Buchwald, Westin and Akerström50,Reference Christakos, Dhawan and Verstuyf51) . A meta-analysis revealed a negative correlation between the 25(OH)D concentration and total cancer incidence and mortality(Reference Yin, Ordóñez-Mena and Chen52). Additionally, the results of a recent animal study indicated that the growth of parathyroid tumours was accelerated by vitamin D deficiency(Reference Costa-Guda, Corrado and Bellizzi53). These studies were consistent with our findings. Our research also indicated that thyroid cancer had a suggestive association on vitamin D levels. Most previous studies have focused on the low serum concentration of vitamin D in patients with thyroid cancer(Reference Zhao, Wang and Zhang54). However, as multiple recent studies have suggested the anti-tumour effects of vitamin D(Reference Buchwald, Westin and Akerström50,Reference Christakos, Dhawan and Verstuyf51,Reference Carlberg and Muñoz55) , further studies are necessary to validate these findings.
The results of our study indicated that thyroid nodules and chronic thyroiditis might have a causal association on vitamin D. Patients with thyroid nodules had serum vitamin D levels that were lower than the normal range(Reference Aboelnaga, Elshafei and Elsayed56,Reference Du, Liu and Zhao57) . Kim et al. discovered that patients with larger nodules had lower vitamin D levels and that there was a negative correlation between vitamin D levels and nodule diameter(Reference Kim, Kim and Kim58). There are several potential explanations for the causal effect of thyroid nodules and chronic thyroiditis on vitamin D levels. A review of clinical data indicated that vitamin D deficiency was a common occurrence in patients with thyroid disease. Research has indicated that ethnicity, geography and seasonal factors might be associated with vitamin D deficiency, but the precise effect of each factor remains unclear(Reference Diffey59,Reference Datta, Pal and De60) . Furthermore, thyroid-stimulating hormone levels have been proposed to be a risk factor for thyroid nodules(Reference Fiore and Vitti61). Additionally, vitamin D deficiency has been linked to an elevated risk of impaired thyroid hormone sensitivity(Reference Zhou, Wang and Su62), which suggests that there might be an association between thyroid-stimulating hormone, thyroid nodules and vitamin D. However, the underlying mechanisms of the effects of thyroid nodules and chronic thyroiditis on vitamin D remain unclear and require further investigation.
In conclusion, the complex mechanisms that link vitamin D to thyroid- and parathyroid-related diseases necessitate further investigation. Observational studies suffer from reverse causality and are limited by confounding factors, thereby limiting their ability to identify aetiologic explanations(Reference Lawlor, Davey Smith and Kundu63). The present study used MR to minimise confounding factors and establish a causal association. The strengths of the study included the use of MR, which overcomes the limitations of observational studies in terms of confounders and causality associations. Additionally, all F-statistics were greater than 10, indicating a reduced susceptibility to weak instrument bias.
However, our study had several limitations. First, the data from our study were obtained from summary GWAS data, and specific information necessary for further analysing age, sex and time of blood collection among the study population is lacking. Second, the data in the FinnGen database were derived from the primary diagnosis of the study population. However, patients might have additional comorbidities, which could lead to biased results. Third, because the number of SNP included in some diseases was insufficient for MR analysis, the thresholds for selecting SNP were appropriately lowered. Additionally, some of the available SNP for the exposure-disease associations were low, which may have affected the study conclusions. Fourth, while we have employed various approaches to minimise pleiotropy, potential unidentified pathways and confounders between exposure and outcome might still lead to inaccuracies in our findings. Fifth, the study subjects were primarily of European origin (the study’s participants were primarily from FinnGen and the UK database), which limits the generalisability of the study findings to the broader European population, and there may be genetic variations between different races. Whether the findings from this research can be extended to other racial groups remains uncertain, and additional corroboration of the outcomes is necessary in the future.
Conclusion
Our study confirmed that chronic thyroiditis and thyroid nodules impact the vitamin D levels, although the underlying mechanisms require further investigation. However, the relationships of serum vitamin D levels on benign parathyroid adenoma and the relationship of thyroid cancer on vitamin D levels still need to be examined in larger multicentre GWAS that include more SNP to validate or reassess our conclusions. Moreover, our study did not reveal a causal relationship between vitamin D and other thyroid- and parathyroid-related diseases.
Acknowledgements
Our thanks go out to the original GWAS participants and investigators, as well as to the FinnGen and UK Biobank research group for sharing and managing the summary statistics.
This study was supported by the Natural Science Foundation of Fujian, China (nos. 2022J011004, 2023J011207 and 2021J01397) and the Fujian provincial health technology project (no. 2022GGA010).
L. R. Z. was responsible for the conceptualisation, data curation, formal analysis, investigation, methodology, resources, software, writing – original draft, visualisation and validation. C. T. H., M. X. L., W. H. L. W and H. T. L. were responsible for the conceptualisation and writing – review and editing. X. X. W., J. Q. C. and H. S. were responsible for writing – review and editing and funding acquisition. All authors have read and approved the final manuscript.
The authors have no potential, perceived or real conflicts of interest to disclose.
Information on the datasets covered in this manuscript can be found in https://gwas.mrcieu.ac./datasets/ebi-a-GCST90000614/uk and r9.risteys.finngen.fi.
Supplementary material
For supplementary material/s referred to in this article, please visit https://doi.org/10.1017/S0007114524001843.