Percentage of variance in rank transformed NMR explained by the independent SNPs, as well as the known CYP2A6 alleles [CYP2A6*2 (= rs1801272), CYP2A6*9 (= rs28399433), CYP2A6*14 (= rs28399435), CYP2A6*18 (= rs1809810), CYP2A6*21 (= rs6413474), allele naming according to the CYP2A6 Allele Nomenclature (www.cypalleles.ki.se/cyp2a6.htm)) included in our 1000Genomes imputed data set, was estimated by running linear regression models separately in each of the three cohorts, with NMR as the dependent variable and age, sex, BMI, and SNPs as explanatory variables. The variance explained by a specific SNP or SNPs was acquired by subtracting the R2 from the model including only the covariates (age, sex, BMI) from the R2 of the model including also the SNP/SNPs of interest as explanatory variables.