Primary analysis was performed using Generalized Summary Mendelian Randomisation
40 (link) (GSMR), implemented in GCTA (v1.92.0). Compared with other summary-statistic based estimators, GSMR is reported to have greater statistical power because GSMR uses genome-wide data to account for sampling variance in SNP-exposure and SNP-outcome estimates
40 (link).
Genome-wide summary statistics for seven metabolic traits and two cardiovascular outcomes were selected as follows: BMI; meta-analysis of GIANT with UK Biobank
65 (link); for waist-to-hip ratio adjusted for BMI; GIANT consortium
66 (link), for coronary artery disease CARDIOGRAM plus C4D meta-analysis
67 (link); for stroke the MEGASTROKE consortium
68 (link); for type 2 diabetes, the DIAGRAM consortium
69 (link); for fasting glucose the ENGAGE consortium
70 (link); for HDL cholesterol, LDL cholesterol and triglycerides, the GLGC data sets were used
71 (link).
To achieve adequate variance explained and therefore statistical power, multiple SNPs were included as instrumental variables for all traits. All variants associated with dental disease or metabolic traits (
P < 5 × 10
−8) after LD clumping using reference data from the cohorts arm of the UK10K project (
https://www.uk10k.org/data.html) to produce index SNPs (r
2 threshold, 0.01) were included as potential instrumental variables. No attempt was made to manually screen variants for possible undesirable pleiotropic effects, with filtering instead performed as part of the HEIDI-outlier procedure.
HEIDI-outlier filtering, is an extension of the heterogeneity in dependent instruments method
72 (link). This analysis projects a plausible distribution of causal effect (βxy) estimated from a non-outlying genetic instrumental variable, and tests whether other single nucleotide polymorphisms (SNPs) have values of βxy compatible with this estimate, on the assumption that pleiotropic variants will have outlying values
40 (link). Following default criteria, potentially outlying variants with (
P < 0.01 for pleiotropy) were removed and βxy was estimated from the remaining instruments.
In the primary analysis, variants in the HLA region (chr6: 25–35 Mb) were removed for estimation of casual effects of DMFS/dentures, except for a single-lead variant. For estimation of causal effects of other traits on DMFS/dentures, full genome-wide data were used.
Methods used for sensitivity analysis are described in Supplementary Note
4.
Shungin D., Haworth S., Divaris K., Agler C.S., Kamatani Y., Keun Lee M., Grinde K., Hindy G., Alaraudanjoki V., Pesonen P., Teumer A., Holtfreter B., Sakaue S., Hirata J., Yu Y.H., Ridker P.M., Giulianini F., Chasman D.I., Magnusson P.K., Sudo T., Okada Y., Völker U., Kocher T., Anttonen V., Laitala M.L., Orho-Melander M., Sofer T., Shaffer J.R., Vieira A., Marazita M.L., Kubo M., Furuichi Y., North K.E., Offenbacher S., Ingelsson E., Franks P.W., Timpson N.J, & Johansson I. (2019). Genome-wide analysis of dental caries and periodontitis combining clinical and self-reported data. Nature Communications, 10, 2773.