Vitamin D and cause-specific vascular disease and mortality: a Mendelian randomisation study involving 99,012 Chinese and 106,911 European adults

Background Randomised control trials and genetic analyses have demonstrated that vitamin D or 25-hydroxyvitamin D (25[OH]D) levels may not play a causal role in the development of cardiovascular disease. However, it is unclear if 25(OH)D is causally associated with cause-specific vascular disease and lipids. Therefore, we examined the causal association of 25(OH)D with myocardial infarction, stroke, ischaemic heart disease, ischaemic stroke, subarachnoid haemorrhage, intracerebral haemorrhage, and lipid levels among both Chinese and Europeans. Methods We used a Mendelian randomisation (MR) design in the China Kadoorie Biobank, the Copenhagen City Heart Study, and the Copenhagen General Population Study. The 25(OH)D-related genetic variants in the CYP2R1 and DCHR7 genes were genotyped in 99,012 Chinese adults and 106,911 Danish adults. Results In Chinese adults, plasma 25(OH)D levels were not significantly associated with cause-specific vascular disease or mortality, with the exception of intracerebral haemorrhage (HR, 1.09 [95% CI, 1.01,1.18] per 25 nmol/L higher plasma 25(OH)D). In Europeans, plasma 25(OH)D levels were inversely associated with all types of vascular diseases and mortality. However, MR analysis did not demonstrate causal associations of genetically increased 25(OH)D levels with cause-specific vascular diseases, or mortality in both Chinese and European adults. In addition, each 25 nmol/L higher 25(OH)D was observationally associated with lower total cholesterol and low-density lipoprotein cholesterol levels, but higher high-density lipoprotein cholesterol levels. Likewise, MR analysis showed that 25(OH)D levels were not causally associated with lipids in both Chinese and European adults after Bonferroni correction. Conclusions We found no evidence to support that genetically increased 25(OH)D was associated with a lower risk of ischaemic stroke, intracerebral haemorrhage, subarachnoid haemorrhage, and lipid levels in both Chinese and European adults. These results suggest that the inverse associations of vitamin D with vascular disease could be the result of confounding. Electronic supplementary material The online version of this article (10.1186/s12916-019-1401-y) contains supplementary material, which is available to authorized users.


Background
Observational studies have shown that low plasma levels of vitamin D are associated with higher risks of cardiovascular disease (CVD) [1][2][3] and mortality from CVD, cancer, and all causes [2,4,5] in European populations. However, randomised intervention trials (RCTs) of vitamin D3 (cholecalciferol) supplementation have not clearly demonstrated any protective effect on CVD [6][7][8][9] or all-cause mortality [10,11]. Findings from Mendelian randomisation (MR) analysis also did not indicate causal association of 25hydroxyvitamin D (25[OH]D), usually used to assess vitamin D status, with increased risk of ischaemic heart disease, myocardial infarction [12], and cardiovascular mortality [13], but showed a possible effect on all-cause and cancer mortality in Europeans [13].
No study has evaluated these associations in the Chinese population, while two studies examined risk of diabetes [14,15]. Importantly, the causal associations between 25(OH)D and cause-specific vascular disease such as stroke, ischaemic stroke, intracerebral haemorrhage, subarachnoid haemorrhage, and lipid levels have not been examined among European and Chinese adults. Whether a high 25(OH)D level is a cause of decreased vascular disease and mortality or simply a consequence of poor health is thus unclear. This question has important public health relevance, as about half of the population have vitamin D deficiency and regularly take vitamin D supplements [16], but without proper scientific evidence for any beneficial effect on many disease outcomes.
With the lack of evidence from intervention trials, a MR approach, which is regarded as a natural RCT, can potentially be used to make causal inference [17][18][19][20][21][22]. MR uses the random assortment of genes during gamete formation as the natural randomisation process. Therefore, genetic variants known to influence 25(OH)D concentrations can be used to examine causality as they partly determine lifetime exposure to 25(OH)D. Genome-wide association studies (GWAS) have identified several single nucleotide polymorphisms (SNPs) that are associated with plasma 25(OH)D concentration in relation to genes involved in the synthesis of 25(OH)D (DHCR7 and CYP2R1), transport of 25(OH)D via vitamin D-binding protein (GC/ DBP), and metabolism of 25(OH)D (CYP24A1) [23,24].
The aims of the present study were to test whether higher 25(OH)D concentrations are observationally associated with cause-specific vascular diseases, or mortality and lipids, and whether such associations are causal among European and Chinese adults.

Study populations
The China Kadoorie Biobank (CKB) is a prospective study of 512,891 individuals, aged 30-79 years, who were enrolled from 5 urban and 5 rural regions in China during [2004][2005][2006][2007][2008]. Previous studies have reported the details of the study design, participant characteristics, and survey methods [21]. Information on demographic data, lifestyle, use of medications, and medical history were collected by trained health workers by using laptop-based questionnaires. Blood pressure, weight, and height were measured using standard protocols. Body mass index (BMI) was calculated as body weight in kilogrammes divided by height in metres squared. A 10-mL nonfasting venous blood sample was collected. All participants from CKB provided written informed consent. Ethics approval was obtained from the Peking University and University of Oxford ethics committees.
The Copenhagen General Population Study (CGPS) is a prospective cohort study among Danish population initiated in 2003 [25]. The Copenhagen City Heart Study (CCHS) was initiated in 1976-1978, and followed up in 1981-1983, 1991-1994, and 2001-2003 [13]. Body weight and height was measured by standard methods. Blood pressure was measured in millimetres of mercury with a sphygmomanometer or an automated digital blood pressure monitor (Kivex). Leisure time physical activity was considered in four categories. Smoking information was classified as pack-years smoked. Alcohol consumption was based on average weekly consumption of beer, wine, and liquor. Written informed consent was obtained from participants. Ethics approval was obtained from Danish ethical committees.

Biochemical analysis
In CKB, plasma 25(OH)D concentrations were measured using a Beckman Access-2 immunoassay assay on stored EDTA plasma samples in 13,565 individuals. The laboratory participated in the international Vitamin D Quality Assessment Scheme (DEQAS) for 25(OH)D and results had a mean (SD) bias of − 11.8% (7.5) from the target mean value for the present study. Samples of the CCHS were stored at − 20°C (1981)(1982)(1983), and samples from the CGPS were stored at − 80°C (2004)(2005). Colorimetric assays (Mannheim, Boehringer Mannheim, Germany or Konelab, Finland, Espoo) were used to measure plasma high-density lipoprotein (HDL) cholesterol and total cholesterol concentrations. 25(OH)D concentrations were measured in both serum and plasma using a chemiluminescent immunoassay, CLIA (DiaSorin, Stillwater, MN, USA). The inter-and intra-assay coefficients of variation were 11% and 10%, respectively.

SNP selection and genotyping
A total of 95,680 individuals from the CKB were genotyped using a 384-SNP panel (Illumina Golden Gate). Two synthesis SNPs (CYP2R1-rs10741657 and DHCR7-rs12785878), one transport SNP (GC/DBP-rs2282679), and one metabolism SNP (CYP24A1-rs6013897) identified for plasma 25(OH)D concentrations in previously GWAS were included in the SNP panel [20,21]. The genotyping concordance was > 99.9% for 2063 pairs of sample replicates and genotype success rates were 99.9% for each SNP. In the CGPS, genetic variants for 25(OH)D were selected based on a reported GWAS [23,24]. The selected SNPs either influence synthesis of pre-vitamin D from 7dehydrocholesterol in the skin or activate transformation from vitamin D to 25(OH)D in the liver. TaqMan assays were used for genotyping for DHCR7 (rs7944926 and rs11234027) and CYP2R1 (rs10741657 and rs12794714). Sequencing of randomly selected samples verified the genotypes. All genotype call rates were higher than 99% complete.
In the CGPS and CCHS studies, information on diagnoses of vascular diseases such as myocardial infarction (ICD-8: 410; ICD-10: I21-I22) and ischaemic heart disease (ICD-8: 410 to 414, and ICD-10: I20 to I25) were verified based on information on diagnoses entered in the national Danish Patient Registry and hospital admissions. The national Danish Causes of Death Registry records were used to identify the causes of death as well as contributing causes of death [12]. We used underlying causes of death to classify deaths as due to cardiovascular (ICD-8 390-458, ICD-10 I00-I99) and cancer (ICD-8 140-209, ICD-10 C00-C97) [13].

Statistical methods
For observational analyses, data were available on 13,565 participants with plasma 25(OH)D levels for observational analyses, and data were available on 82,464 participants with all genotypes. Of these, 3397 had data on both 25(OH)D plasma levels and genotypes in CKB. We included 25,621 participants and 10,271 participants with plasma 25(OH)D levels in CGPS and CCHS, respectively; in genetic analyses, we included 106,911 participants with all genotypes. Of these, 25,332 had both 25(OH)D plasma levels and genotypes in CGPS. Cox proportional hazard regression models were used to analysing the relationship between quartile groups of 25(OH)D levels and risk of disease. We used Schoenfeld residual test to examine the proportional hazard assumption. All results (P > 0.05) showed that the Cox models were acceptable. Multivariable Cox regression analyses included age (years), sex (male or female), smoking status (current smoker), alcohol intake (current drinker), season, region, systolic blood pressure (SBP), physical activity (low activity) in the CKB, and age (years), sex (male or female), season, SBP, physical activity (low activity), smoking status (current smoker), and alcohol intake (current drinker) in the CGPS and CCHS.
Genotype distributions in each study did not differ from the Hardy-Weinberg equilibrium. Linear regression was used to examine the per allele effects of each SNP on plasma 25(OH)D concentrations. Genetic scores were estimated for the two synthesis SNPs and all four SNPs. The selected SNPs are not in linkage disequilibrium (LD). F-statistics were used to estimate the strength of associations of each SNP with 25(OH)D concentrations, and F-statistics > 10 were considered strong. We also used linear regression to assess the relationship of each SNP with plasma lipids and lipoproteins. All genetic analyses were adjusted for age at baseline, sex, season, and region. Analyses of individual SNPs and genetic scores with outcomes were estimated using Cox proportional hazard regression models. The per allele effects of each SNP on plasma 25(OH)D concentrations were expressed as the differences in 25(OH)D concentrations per copy of the 25(OH)D raising allele.
A two-stage least squared regression model using genotype individually and combined allele scores as instruments was used to examine the influence of a 25 nmol/L increase in 25(OH)D concentrations on risk of outcomes. We used the DHCR7/CYP2R1 genetic score as the main instrumental variable. We calculated causal estimates of genetically determined hazard ratios (HRs) by using the Wald-type estimator, which calculates the ratio of the outcome allele score log hazard ratios to the exposure allele score coefficient and then exponentiates the ratio to a hazard ratio. For combination of CVD and lipid data from Asians and Europeans, meta-analyses were conducted. We assessed between-study heterogeneity via Cochrane's Qand I 2 -statistics. A random-effect meta-analysis was used if I 2 > 0.25, otherwise, a fixed-effect model was used [26].
Sensitivity analyses were performed for three-SNP and four-SNP genetic scores to examine the potential effects of pleiotropy using the MR-Egger regression method [27]. We used SAS version 9.2 and R version 3.01. We included 20 diseases or cause-specific mortality as primary outcomes. The threshold of significance was P < 0.0025 (0.05/20) after Bonferroni correction. P ≤ 0.05 but above the threshold of Bonferroni corrected significance was considered as a suggestive causal association.

Characteristics of the CKB, CGPS, and CCHS populations
Among the 82,464 CKB individuals and 106,911 CGPS individuals in the genetic study, the mean (SD) of age was 51.4 (10.6) and 58.0 (13.1) years, 61% and 55% were women, and the mean (SD) of BMI was 23.7 (3.4) and 26.1  (Table 1).

Observational associations
In CKB, plasma 25(OH)D levels were not statistically associated with major vascular events (multivariable adjusted HR per 25 nmol/L higher plasma 25 (Table 2). Likewise, plasma 25(OH)D levels were not significantly associated with cancer, respiratory diseases, infection-cause mortality, or all-cause mortality after adjustment for confounding factors (Table 3). In addition, each 25 nmol/ L higher plasma 25(OH)D concentration was associated with a 1.52 mg/dL higher apoA, 0.04 mg/dL higher apoB, and 0.03 mmol/L (1 mg/dL) higher HDL cholesterol, but a 0.05 mmol/L (2 mg/dL) lower total cholesterol, 0.03 mmol/L (1 mg/dL) lower LDL cholesterol, and a 0.16 mmol/L (14 mg/dL) lower triglycerides (Table 4).
In CGPS and CCHS, each 25 nmol/L higher plasma 25(OH)D concentration was associated with about a 10%  (Table 3).  Table S1). In CGPS, genetic variants and genetic score were also significantly associated with plasma 25(OH)D concentrations ( Fig. 1 and Additional file 1: Table S1). None of the SNPs was associated with differences in confounders such as age, gender, season, region, smoking, drinking, physical activity, blood pressure, and BMI in CKB or in CGPS (Additional file 1: Tables S2 & S3).
Causal associations with cause-specific vascular disease or mortality In both CKB and CGPS, there was no significant genetic association of the DHCR7/CYP2R1 risk score with cause-specific vascular disease (Additional file 1: Tables S4 and S5) or mortality (Additional file 1: Table S6 and S7) Multivariable cox proportional hazard regression models were used to examine the association between quartile groups of 25(OH)D concentrations and risk of vascular diseases All values are adjusted for age (years), sex (male or female), smoking status (current smoker), alcohol intake (current drinker), season, region, SBP, and physical activity (low activity) in CKB All values are adjusted for age (years), sex (male or female), latitude, season, SBP, physical activity (METs, h/day), smoking status (never smoker, occasional smoker, former smoker, or regular smoker), and alcohol intake (non-drinker, occasional drinker, former drinker, or regular drinker) in CGPS and CCHS after Bonferroni correction. Instrumental variable analysis did not show any significant causal association of genetically predicted 25(OH)D with any cause-specific vascular disease in CKB and CGPS (Fig. 2). Furthermore, genetically predicted 25(OH)D was not significantly associated with cause-specific mortality in CKB and CGPS (Fig. 3). However, the instrumental variable analysis indicated a marginal association for all-cause mortality per a 25 nmol/ L higher in genetically predicted 25(OH)D concentrations (0.98 [0.96,1.00]) (Fig. 3). Sensitivity analyses yielded similar findings for genetically instrumented differences in 25(OH)D concentrations using the individual SNP or four-SNP score for each disease (Additional file 1: Table  S4 to S7). In addition, instrumental variable analyses showed significant associations with total and LDL cholesterol in CKB ( Fig. 4 and Additional file 1: Table S8). Each 25 nmol/l genetic increase in 25(OH)D concentrations was associated with 0.058 mmol/L (0.22 mg/dL) and 0.034 mmol/L (0.13 mg/dL) lower total and LDL cholesterol in Multivariable cox proportional hazard regression models were used to examine the association between quartile groups of 25(OH)D concentrations and risk of mortality All values are adjusted for age (years), sex (male or female), smoking status (current smoker), alcohol intake (current drinker), season, region, SBP, and physical activity (low activity) in CKB All values are adjusted for age (years), sex (male or female), latitude, season, SBP, physical activity (METs, h/day), smoking status (never smoker, occasional smoker, former smoker, or regular smoker), and alcohol intake (non-drinker, occasional drinker, former drinker, or regular drinker) in CGPS and CCHS CKB. However, we did not observe causal associations for lipids and lipoproteins in CGPS and pooled cohorts ( Fig. 4 and Additional file 1: Table S9).

Discussion
The novel findings of this MR study involving 99,012 Chinese adults and 106,911 European adults suggest that genetically higher plasma 25(OH)D levels were not associated with lower risk of cause-specific vascular disease such as myocardial infarction, ischaemic heart disease, ischaemic stroke, stroke, intracerebral haemorrhage, and subarachnoid haemorrhage among both Chinese and European adults. Likewise, 25(OH)D may not influence lipid levels. Our results suggest that the inverse associations of vitamin D with cause-specific vascular disease and lipids could be the results of confounding.
Using an MR approach free from confounding and reverse causation, our study for the first time showed that genetically predicted 25(OH)D was not associated with intracerebral haemorrhage, subarachnoid haemorrhage, stroke, and ischaemic stroke, thus suggesting that vitamin D may not play a causal role in development of cause-specific vascular disease in Chinese or Europeans.
The results suggest that individuals should be cautious about long-term supplementation with high doses of Linear regression was also used to assess the associations of 25(OH)D with lipids All values are adjusted for age (years), sex (male or female), smoking status (current smoker), alcohol intake (current drinker), season, SBP, and physical activity (low activity) in CKB All values are adjusted for age (years), sex (male or female), latitude, season, region, SBP, physical activity (METs, h/day), smoking status (never smoker, occasional smoker, former smoker, or regular smoker), and alcohol intake (non-drinker, occasional drinker, former drinker, or regular drinker) in CGPS and CCHS *Data only available in CGPS vitamin D3 [10]. The results are consistent with recent MR studies which did not show a causal effect of circulating vitamin D levels on risk of ischaemic heart disease or myocardial infarction [12], or cardiovascular mortality [13]. Findings from RCTs that do not show a protective effect of vitamin D3 supplementation on cardiovascular outcomes [9][10][11], supported our findings. However, previous RCTs did not specifically focus on vitamin Ddeficient population. Thus, available evidence is not enough to support cardiovascular benefits or hazards of the commonly used vitamin D doses, but whether higher doses of vitamin D are beneficial is uncertain, or whether vitamin D supplementation has vascular effects in individuals with vitamin D deficiency remains unclear [9]. It worth noting that there are two major forms of vitamin D: vitamin D2 (ergocalciferol) and vitamin D3 (cholecalciferol). A meta-analysis demonstrated that vitamin D3 is more effective than vitamin D2 in maintaining circulating concentrations of 25(OH)D [28]. Previous randomised double-blind placebo-controlled trial examining the effects of vitamin D2 or D3 supplementation on cardiometabolic risk showed that both vitamin D2 or D3 did not influence blood pressure and cardiovascular risk, but decreased ApoB concentration and pulse wave velocity. Only vitamin D2 decreased total cholesterol [29]. However, so far, no RCT has investigated the benefits of supplemental vitamin D2 or D3 in preventing vascular diseases in Chinese adults. Whether vitamin D2 or D3 supplementation will be beneficial in patients with CVD or high-risk older patients is unknown and warrants further investigation. Our findings could be explained by previous evidence that low 25(OH)D concentration is a marker of an unhealthy lifestyle or poor health [30], which has been associated with increased risk of vascular disease. Furthermore, previous MR studies suggest that cardiovascular factors [31], such as high remnant cholesterol, LDL cholesterol [32], and obesity [33], decrease plasma concentrations of 25(OH)D and increase vascular risk. That means low Fig. 2 Instrumental variable estimates for vascular diseases. Analyses of two-SNP score with vascular diseases were estimated using cox proportional hazard regression models. We used a two-SNP score as instrument to estimate the influence of a 25 nmol/L increase in 25(OH)D concentrations on risk of vascular diseases. We calculated instrumental variable estimates of genetically determined hazard ratios by using the Wald-type estimator, which involves taking the ratio of the gene-outcome log hazard ratios to the gene-exposure coefficient and then exponentiating to express it as a hazard ratio. Two-SNP score was calculated based on DHCR7 + CYP2R1: rs12785878 + rs10741657 in CKB and DHCR7 + CYP2R1: rs7944926 + rs10741657 in CGPS. The r 2 between rs12785878 and rs7944926 is 0.87. All values are adjusted for age, sex, and season and stratified by region plasma 25(OH)D is a consequence of CVD rather than a causal factor. This could explain why we did not observe a causal association in MR analyses, despite strong observational associations in Europeans [4,5]. In addition, a previous study showed that elevated parathyroid hormone levels identify a population with higher cardiovascular risk. Parathyroid hormone correlates to cardiovascular risk, but not vitamin D [34]. The potential mechanism is that parathyroid hormone induces oxidative stress and affects endothelial function [35]. Therefore, we have to acknowledge that the vitamin D-dependent parathyroid hormone, rather than vitamin D itself, might be a better predictor of cardiovascular risk. However, we have to acknowledge that the possibility of a non-linear association of vitamin D with cardiovascular outcomes cannot be excluded as it was not tested in the present MR study.
Although our MR analysis did not show causal effects of 25(OH)D on all-cause and cause-specific mortality in Chinese, results from Europeans demonstrated that 25(OH)D might be causally related to all-cause mortality and cancer mortality. Our findings are consistent with a European MR study, which showed that genetically reduced vitamin D concentrations were associated with increased all-cause, cancer mortality, but not with cardiovascular mortality [13]. Our results in Europeans are also supported by a meta-analysis of 56 RCTs mainly conducted in Europeans, which showed that vitamin D decreased allcause and cancer mortality [11]. However, when different forms of vitamin D were assessed, only vitamin D3 (cholecalciferol) decreased mortality. Vitamin D2 (ergocalciferol) may even increase mortality [11]. It should be noted that recent European MR studies found evidence that genetically determined vitamin D concentrations were inversely associated with type 2 diabetes [36], and blood pressure and risk of hypertension [37]. These results could at least partly explain the reduced mortality among Europeans observed by us. However, so far, neither RCTs nor MR studies have investigated the causal effect of vitamin D on all-cause and cause-specific mortality in Chinese adults. Our present MR analysis for the first time documented Fig. 3 Instrumental variable estimates for mortality. Analyses of two-SNP score with mortality were estimated using cox proportional hazard regression models. We used a two-SNP score as instrument to estimate the influence of a 25 nmol/L increase in 25(OH)D concentrations on risk of mortality. We calculated instrumental variable estimates of genetically determined hazard ratios by using the Wald-type estimator, which involves taking the ratio of the gene-outcome log hazard ratios to the gene-exposure coefficient and then exponentiating to express it as a hazard ratio. Number of individuals in CKB and CGPS are 99,012 and 106,911, respectively. Two-SNP score was calculated based on DHCR7 + CYP2R1: rs12785878 + rs10741657 in CKB and DHCR7 + CYP2R1: rs7944926 + rs10741657 in CGPS. The r 2 between rs12785878 and rs7944926 is 0.87. All values are adjusted for age, sex, and season and stratified by region that genetically elevated 25(OH)D was not associated with cause-specific vascular disease or mortality. These findings suggest possible ethnic differences in relationships between 25(OH)D levels and cancer mortality. A previous study showed that the effect of vitamin D3 supplementation on advanced colorectal adenomas varied according to vitamin D receptor genotypes [38]. Among the rs7968585AA genotype individuals (26%), supplemental vitamin D3 reduced risk by 64%, but increased risk by 41% among individuals with rs7968585 G alleles (74%) [38]. Therefore, whether benefits from vitamin D supplementation for the prevention of cancer vary according to vitamin D receptor genotype is unknown and warrants further investigation.
Previous RCTs have shown inconsistent results for the effects of vitamin D on plasma lipid levels [39][40][41][42][43][44]. However, it is worth noting that none of the RCTs were specifically conducted for evaluating the effect of supplemental vitamin D on lipids and lipoproteins. Interestingly, the largest trial (The WHI-Calcium+Vitamin D trial) conducted in 36,282 women showed that patients receiving combined vitamin D and omega-3 fatty acid supplementation had decreased serum very low-density lipoprotein cholesterol and serum triglycerides [42] while 400 IU vitamin D3 plus 1000 mg elemental calcium carbonate daily showed a decrease in LDL cholesterol after 6 weeks of intervention [43]. Consistently, our MR analysis in Chinese adults for the first time showed that each 25 nmol/l genetic increase in 25(OH)D concentration was associated with 0.058 mmol/L(0.22 mg/dL) and 0.034 mmol/ L(0.13 mg/dL) lower total and LDL cholesterol, thus suggesting that 25(OH)D may play a causal role in lipid metabolism at least in Chinese adults. However, the potential mechanism warrants further investigation.
In the present MR analysis, we considered three assumptions to ensure the validity of our analyses [17,18,27]. First, our data demonstrated that the individual SNP and genetic score were associated with 25(OH)D (F statistic > 15), suggesting strong instrumental variables. Second, the instruments should be independent of measured and unmeasured confounders. Indeed, we found no association of genetic variants or score with measured confounders, thus satisfying an important condition for a valid MR experiment. Third, the genetic score should have no pleiotropic effect on the outcomes, which is an important consideration in MR studies. The third assumption may be violated by pleiotropy [45]. The SNPs encoding enzymes involved in the synthesis of 25(OH)D are believed to be more specific for 25(OH)D, but SNPs involved in the transport of 25(OH)D may have pleiotropic effects in addition to their effects on plasma 25(OH)D concentrations alone [31,46,47]. However, when we excluded SNPs involved in the transport of 25(OH)D, the association remained unchanged.
Several strengths of the present study merit consideration. This is the first MR analysis to examine the causal association of 25(OH)D with intracerebral haemorrhage and subarachnoid haemorrhage. This large MR study had a large number of disease events and high-quality 25(OH)D measurements by well-trained technicians and an appropriate design allowed us to gain sufficient power for estimation of causal effects. Another strength of our study was detailed information on known confounding factors, which allowed us to examine the three MR assumptions. Finally, we used standardised analytic methods to compare the differences in causal effects between Chinese and Europeans.

Conclusion
In summary, this large-scale analysis provides reliable evidence that higher 25(OH)D may not play a causal role in developing cause-specific vascular disease in Chinese or Europeans. MR analysis further demonstrated that high vitamin D may not influence lipids.

Additional file
Additional file 1: Supplemental tables providing additional information on genetic association with 25(OH)D, lipids and CVD. (DOCX 52 kb) Fig. 4 Instrumental variable estimates for lipids. Linear regression was also used to assess the associations of a two-SNP score with lipids. We used a two-SNP score as instrument to estimate the influence of a 25 nmol/L increase in 25(OH)D concentrations on lipids. We calculated instrumental variable estimates of genetically determined hazard ratios by using the Wald-type estimator, which involves taking the ratio of the gene-outcome coefficient to the gene-exposure coefficient. Two-SNP score was calculated based on DHCR7 + CYP2R1: rs12785878 + rs10741657 in CKB and DHCR7 + CYP2R1: rs7944926 + rs10741657 in CGPS. The r 2 between rs12785878 and rs7944926 is 0.87. All values are adjusted for age, sex, and season and stratified by region