Figure S1: Results of MR-Radial analyses with instruments of IV-M and IV-S. The effect of instrument which only genome-wide associated with 25(OH)D (IV-S) on height was estimated by MR-Radial analysis in a) mixed population, b) males and c) females. Each point represents the SNP effects of 25(OH)D on height. The blue lines represent inverse-variance-weighted (IVW) estimates of the association between a 1-SD increase in 25(OH)D and height, removing the outliers; The effect of variants which genome-wide associated with 25(OH)D and other phenotypes (IV-M) in d) mixed population, e) males and f) females were also analyzed; Table S1: (DOC) Detailed information on the GWASs included in the current study; Table S2, (XLSX) Detailed information on all variants associated with the serum levels of 25(OH)D; Table S3, (XLSX) Detailed information on all variants associated with vitamin D levels; Table S4, (XLSX) Detailed information on all variants associated with the levels of VDBP; Table S5, (DOC) Causal directionality test of primary MR analyses; Table S6, (DOC) Heterogeneity/pleiotropy test results of primary MR analyses; Table S7, (DOC) Genetic effect of 25(OH)D levels on height estimated by the MR-Radial method; Table S8, (DOC) Multiple MR analyses for genetic variants of 25(OH)D and height; Table S9, (DOC) Genetic effect of 25(OH)D levels on height estimation by subinstrument variants MR analysis.