Abstract
Many biomarkers are associated with type 2 diabetes (T2D) risk in epidemiological observations. The aim of this study was to identify and summarize current evidence for causal effects of biomarkers on T2D. A systematic literature search in PubMed and EMBASE (until April 2015) was done to identify Mendelian randomization studies that examined potential causal effects of biomarkers on T2D. To replicate the findings of identified studies, data from two large-scale, genome-wide association studies (GWAS) were used: DIAbetes Genetics Replication And Meta-analysis (DIAGRAMv3) for T2D and the Meta-Analyses of Glucose and Insulin-related traits Consortium (MAGIC) for glycaemic traits. GWAS summary statistics were extracted for the same genetic variants (or proxy variants), which were used in the original Mendelian randomization studies. Of the 21 biomarkers (from 28 studies), ten have been reported to be causally associated with T2D in Mendelian randomization. Most biomarkers were investigated in a single cohort study or population. Of the ten biomarkers that were identified, nominally significant associations with T2D or glycaemic traits were reached for those genetic variants related to bilirubin, pro-B-type natriuretic peptide, delta-6 desaturase and dimethylglycine based on the summary data from DIAGRAMv3 or MAGIC. Several Mendelian randomization studies investigated the nature of associations of biomarkers with T2D. However, there were only a few biomarkers that may have causal effects on T2D. Further research is needed to broadly evaluate the causal effects of multiple biomarkers on T2D and glycaemic traits using data from large-scale cohorts or GWAS including many different genetic variants.
Introduction
Over the past decade, interest in studying biological markers (biomarkers) for type 2 diabetes (T2D) has increased intensely. This happened because multiple pathobiological processes may contribute to the disease progression, which provides an opportunity to introduce preventive and therapeutic interventions for T2D (1). In clinical practice, such biomarkers (e.g., glucose and glycated haemoglobin tests) are widely used for the diagnosis of diabetes or for the monitoring of therapeutic intervention (2, 3). Targeted intervention at the biomarker level would be useful where there is evidence for a causal relationship between an exposure (like a biomarker) and T2D (4).
Traditional epidemiological studies lack sufficient information to fill the evidence gap due to unmeasured confounding or reverse causality (4, 5, 6, 7). It has been successfully shown that a complementary analysis of genetic data, termed ‘Mendelian randomization,’ has additive value to infer a causal association (4, 5, 6, 7, 8). The main assumption for Mendelian randomization is that the genetic variants do not change over time and are inherited randomly (based on Mendel's laws). In other words, the genetic variants as proxy measures for exposures (e.g., biomarkers) are essentially considered free from confounding and reverse causation. Therefore, the analysis of integrated observational-genetic data is considered similar to that of the randomized trials (9, 10). In Mendelian randomization, if there is a causal association between a biomarker and T2D, the genetic variant(s) influencing the biomarker and the outcome of interest should be associated (5, 6, 8). In the current study, evidence for causal associations between biomarkers and the risk of T2D was updated via a systematic literature search to identify Mendelian randomization studies. Next, summary data from the two largest genome-wide association studies (GWAS) for T2D or glycaemic traits (10, 11) were used to examine the effect estimates for each genetic variant compared with that of the identified studies.
Methods
Search strategy for candidate biomarkers
PubMed and EMBASE were searched to identify Mendelian randomization studies examining the associations between biomarkers and T2D until April 2015. The overview of this systematic literature search was conducted in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-analysis (PRISMA) guidelines, when applicable (12). A manual search was also done for the references of included articles to identify other relevant studies. Because a review of published original studies was performed, the Declaration of Helsinki items related to ‘approval of medical ethics committee’ and ‘permission acquisition’ are not applicable to the current study.
Selection criteria
Studies were included if they formally quantified a causal association between one or more biomarkers (as main exposures) and T2D (as the main outcome); used data on biomarker-associated genetic variant; and classified/defined the exposure as a biomarker that has been objectively measured in serum, plasma or urine. The following MeSH keywords or related terms were used: Diabetes mellitus, Type 2[Mesh], Type 2 Diabetes[tiab], Diabetes[ti]; and Mendelian randomization, Mendelian randomisation, instrumental variable, genetic risk score or causal association.
Data extraction and quality assessment
A primary plan was made to extract necessary data from the full text of the original studies or to contact the corresponding author(s) when appropriate. Fig. 1 depicts the workflow of the literature search. T2D was determined as the main outcome if one or more of the following conditions were fulfilled: a physician diagnosed T2D as indicated by a self-report or in a primary-care database; fasting plasma glucose ≥7.0 mmol/l, a random sample plasma glucose ≥11.1 mmol/l; or the initiation of glucose-lowering medication as retrieved from a pharmacy registry or hospital records (11).
Statistical analysis
All genetic variants that affect identified biomarkers were obtained from the original Mendelian randomization studies. To replicate the nature of the relationship between each biomarker and the outcome (i.e., T2D or glycaemic traits), a genetic approach using outcome-association data for the biomarker-related genetic variants was applied (5, 6, 8). In brief, to determine whether the same genetic variants (or a suitable proxy variant) were associated with T2D (or glycaemic traits), the corresponding summary association statistics from GWAS for T2D (DIAbetes Genetics Replication And Meta-analysis (DIAGRAMv3)) (13) and glycaemic traits (Meta-Analyses of Glucose and Insulin-related traits Consortium (MAGIC)) (14) were extracted. DIAGRAMv3 is a meta-analysis of multiple GWAS with a total number of 12 171 diabetes cases and 56 862 controls of European descent (13). Previously, details regarding the use of GWAS data for Mendelian randomization were described (5, 6, 10, 11). These selected single-test association analyses in the GWAS data are equivalent to that of individual-level data analysis (5, 6, 11).
Extracted data were tabulated on Excel spreadsheets. All statistical analyses were conducted using Excel or Stata/SE version 13.1 for Windows (http://cran.r-project.org/). A two-sided P value <0.05 was considered nominally significant.
Results
Literature search and study characteristics
After scanning 812 titles and selected abstracts, 33 articles were selected for full text review (11, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46). The literature search process is shown in Fig. 1. Five studies were excluded as they did not measure a specific biomarker (n=1) or investigated a measure of insulin resistance as the main exposure (n=4). The characteristics of the 28 studies are summarized in Table 1.The studies were performed in different populations; most were cohort studies including middle-aged adults in the USA or Europe and were published between 2008 and 2015. Four studies were conducted only in Asian populations including Chinese or Taiwanese. In 12 studies, data from multiple GWAS were combined to examine the associations between genotypes and T2D. Six of these studies used data from DIAGRAMv2 (n=3) (47) or DIAGRAMv3 (n=3) (13), and the other six used at least two cohorts with genome-wide genotyping. DIAGRAMv2 was a meta-analysis of eight GWAS comprising 8130 individuals with T2D and 38 987 controls of European descent (47). All Mendelian randomization studies investigated only one biomarker as an exposure in relation to T2D. Some studies also examined the causal associations of a single biomarker with several other outcomes such as cardiovascular disease (18, 24), rheumatoid arthritis (18) or osteoporosis (24). In the final sample, studies included up to 28 144 T2D cases and 76 344 T2D controls or total participants.
Mendelian randomization studies for biomarkers and type 2 diabetes.
Biomarkers | No. of study | Study | Country/year | Cases/total or controls | Genetic variant(s)/locus | Causal estimate | Type of associations |
---|---|---|---|---|---|---|---|
Adiponectin | 3 (19, 35, 42) | DIAGRAMv2, ARIC, FUSION, WTCCC, EPIC-InterAct | US, Europe/2013 | 15 960/64 731 | rs17300539, rs1736653, rs3774261, rs3821799/ADIPOQ | Null | |
BHS, CUDAS, FDS | Australia/2013 | 967/2355 | rs12637534, rs16861209, rs17366568/ADIPOQ | OR, 0.86 | Per allele | ||
DIAGRAMv2 | US, Europe/2012 | –/22 044 | rs3001032/LYPLAL1, rs1108842/GNL3, rs1597466/TSC22D2, rs6810075/ADIPOQ, rs998584/VEGFA, rs2980879/TRIB1, rs7955516/PDE3A, rs601339/GPR109A, rs7133378/DNAH10, rs2925979/CMIP, rs12922394/CDH13, rs731839/PEPD | β, −0.30 | Average additive effect of adiponectin- raising alleles | ||
Beta-carotene | 1 (33) | WTCCC, DGI, FUSION | USA, Europe/2009 | 4549/5579 | rs6564851/BCMO1 | Null | |
Bilirubin | 1 (11) | PREVEND | Netherlands/2015 | 210/3381 | rs6742078/UGT1A1 | OR, 0.58 | 1-s.d. |
N-terminal pro B-type natriuretic peptide | 1 (37) | DIAGRAMv2, Norfolk Diabetes, Cambridgeshire, ADDITION-Ely, French, Swiss 1, Swiss 2, Austrian, MONICA, INCA, UK | US, Europe/2011 | 23 382/57 898 | rs198389/NPPB | OR, 0.96 | C allele |
C-reactive protein | 1 (16) | Whitehall II | UK/2008 | 354/5274 | rs1130864, rs1205, rs3093077/ CRP | Null | |
Delta-6 desaturase | 1 (30) | EPIC-Potsdam | Germany/2011 | 673/2724 | rs174546/FADS1 | RR, 0.64 | TT genotype |
Dimethylglycine | 1 (32) | MDC-CC, MDC, MPP | Sweden/2015 | 4201/33 898 | rs2431332/ DMGDH | OR, 1.1 | A allele |
Ferritin/TMPRSS6 | 1 (22) | Beijing | China/2012 | 272/1574 | rs855791, rs4820268/ TMPRSS6 | OR, 0.79 | AG haplotype |
Fetuin-a | 1 (27) | CHS | US/2013 | 259/3093 | rs4917, rs2248690/AHSG | Null | |
Homocysteine | 2 (24, 25) | Hangzhou | China/2014 | 774/500 | rs1801131/MTHFR, rs4646356/PEMT | OR, 1.93, 1.52 | CC, CT+TT |
– | Europe, Asia, Africa/2013 | 4011/4303 | rs12134663, rs1801133/MTHFR | OR, 1.29 | 5 μmol/l | ||
IL-1Ra | 1 (26) | DIAGRAMv3 and EPIC-InterAct | US, Europe/2015 | 18 715/61 692 | rs6743376, rs1542176/IL1RN | Null | |
Leukocyte telomere length | 1 (46) | WHI-OS | US/2012 | 1675/2382 | rs34368910/TRF1, rs4888444/TRF2, rs4975605/TERT, rs938886, rs2228041,rs12880583/TPP1 | Null | |
Lipoprotein(a) | 2 (29, 44) | Danish general population | Denmark/2014 | 2157/28 567 | rs10455872/LPA | Null | |
DIAGRAMv3,EPIC-Norfolk | USA, Europe/2014 | 10 088/68 346 | rs10455872/LPA | Null | |||
Macrophage migration inhibitory factor | 1 (23) | MONICA/KORA Augsburg | Germany/2008 | 502/1632 | rs1007888/MIF | HR, 1.74 | C allele |
miRNAs | 1 (41) | Harbin | China/2015 | 995/967 | rs895819/miR-27a, rs531564/miR-124a, rs2910164/miR-146a | Null | |
Resistin | 1 (18) | CVDFACTS | Taiwan/2014 US, Europe/2010 | 230/3400 | rs3745367,rs1423096/RETN | OR, 1.38 | GG |
Sex hormone-binding globulin | 2 (21, 34) | WTCCC, Dundee, EFS-Y2TD, Danish, DGI, MCC, FUSION 1-2, KORA, Cambridgeshire, ADDITION-Ely, NDCCS, DeCODE, METSIM, DIAGEN | 27 657/58 481 | rs1799941/SHBG | OR, 0.92 | G allele | |
WHS | USA/2009 | 359/359 | rs6259/SHBG | OR, 0.3 | 1-s.d. | ||
Triglycerides | 1 (20) | Go-DARTS | UK/2011 | 5637/6860 | rs2954029/ TRIB1, rs714052/MLXIPL, rs7557067/APOB, rs17216525/NCAN, rs10889353/ANGPTL3, rs7679/PLTP, rs7819412/XKR6-AMAC1L2, rs328/LPL, rs3135506/APOA5, rs662799/APOA5 | Null | |
Uric acid | 1 (36) | Cambridgeshire, ADDITION-Ely, Norfolk Diabetes | UK/2011 | 7504/8560 | rs12129861/PDZK1, rs742132/LRRC16A, rs505802/SLC22A12, rs12356193/SLC16A9, rs17300741/SLC22A11, rs1165151/SLC17A1, rs2231142/ABCG2, rs734553/ SLC2A9 | Null | |
Vitamin D | 3 (17, 28, 45) | DIAGRAMv3, ADDITION-Ely, Norfolk Diabetes, Cambridgeshire, EPIC-InterAct | US, Europe/2015 | 281 44/76 344 | rs12785878/DHCR7, rs10741657/CYP2R1, rs4588/DBP, rs17217119/CYP24A1 | Null | |
EPIC-Potsdam | Germany/2013 | 1572/2121 | rs2282679, rs1155563/GC, rs12785878, rs3829251/DHCR7, rs10741657/CYP2R1, rs6013897/CYP24A1, rs6599638/Ch10orf88, rs10877012/CYP27B1 | Null | |||
Tromsø | Norway/2012 | 1092/9528 | rs2298850/DBP GC, rs10741657/CYP2R1, rs3794060/NADSYN1, rs6013897/CYP24A1 | Null | |||
Vitamin D binding protein | 1 (31) | CaMos | Canada/2014 | 201/2254 | rs2282679/DBP GC | Null |
RR, relative risk; IL1-Ra, interleukin 1 receptor antagonist; DIAGRAM, DIAbetes Genetics Replication And Meta-analysis; EPIC-InterAct, European Prospective Investigation into Cancer and Nutrition (EPIC)-InterAct; ARIC, Atherosclerosis Risk In Communities; WTCCC, Welcome Trust Case Control Cohort; FUSION, Finland-US Investigation of NIDDM genetics; BHS, Busselton Population Health Survey; CUDAS, Carotid Ultrasound Disease Assessment Study; FDS, Fremantle Diabetes Study; PREVEND, Prevention of Renal and Vascular ENd-stage Disease; iNCA, iNsuffisance CAdiaque; MONICA, Multinational MONItoring of trends and determinants in CArdiovascular disease; MDC, Malmö Diet and Cancer Study; MDC-CC, Malmö Diet and Cancer Cardiovascular Cohort; MPP, Malmö Preventive Project; CVDFACTS, CardioVascular Disease risk FACtors Two-township Study; KORA, Kooperative Gesundheitsforschung in der Region Augsburg; EFS-YT2D, Exeter Family Study-Young-onset T2D; NDCCS, Norfolk Diabetes Case Control Study; METSIM, METabolic Syndrome In Men; DIAGEN, DIAbetes GENetic study; CHS, Cardiovascular Health Study; WHI-OS, Women's Health Initiative Observational Study Cohort; ADDITION, Anglo-Danish-Dutch Study of Intensive Treatment In People with Screen Detected Diabetes in Primary Care; Go-DARTS, Genetics of Diabetes Audit and Research in Tayside Scotland; CaMos, Canadian Multicentre Osteoporosis Study.
Biomarkers and Mendelian randomization
From 28 studies, data were retrieved on causal associations of biomarkers with T2D. In these studies, 21 unique biomarkers that were investigated at least once (16 biomarkers), twice (three biomarkers) and three times (two biomarkers) were identified (Table 1); 19 biomarkers were investigated in independent studies. However, the two biomarkers with three Mendelian randomization studies were investigated in combined studies where the same sample, e.g., DIAGRAMv2 for adiponectin (19, 43), or a part of whole cohort, e.g., European Prospective Investigation into Cancer and Nutrition (EPIC)-Potsdam for vitamin D (17, 45), were used to make total cases/controls. Eleven studies used one genetic variant as a single instrumental variable in Mendelian randomization. Eight studies used at least two independent genetic variants in the same locus as instrumental variables. The rest of the nine studies used multiple genetic variants in different loci or created a multi-locus genetic risk score for each biomarker.
In main or sensitivity analyses, four studies (11, 26, 44, 45) examined causal associations of bilirubin, lipoprotein (a), vitamin D, interleukin-1 receptor antagonist (with corresponding biomarker-associated genotypes) with T2D using the DIAGRAMv3 summary data. In Table 1, the value for each causal estimate can be interpreted as odds ratio (OR), hazard ratio (HR) or log OR (β coefficient) for diabetes per one unit change in genetically determined biomarkers or biomarker-associated genotypes. Evidence of causal association was reported for adiponectin, bilirubin, N-terminal pro B-type natriuretic peptide (NT proBNP), delta-6 desaturase (D6D), dimethylglycine, ferritin (transmembrane protease serine 6), homocysteine, macrophage migration inhibitory factor (MIF), sex hormone-binding protein (SHBG) and resistin (Table 1). For adiponectin, causal estimates were calculated as an OR of 0.86 per allele (35) or a β coefficient of 0.3 for adiponectin multi-locus genotypic risk score (19). Causal estimate for bilirubin was reported as an OR of 0.58 per 1-s.d. genetically increased log-transformed bilirubin (11). For NT proBNP (37), D6D (30), dimethylglycine (32), ferritin (22), MIF (23) and resistin (18), the expected ORs for T2D were 0.96, 0.64, 1.1, 0.79, 1.74 and 1.38 per each copy of risk alleles or genotypes, respectively. For homocysteine, causal estimates were reported as an OR of 1.29/5 μmol/l genetically increased homocysteine (24) or ORs of 1.93 and 1.52 for CC and CC/TT genotypes (25) respectively. For SHBG, the expected ORs for T2D were 0.92 per each copy of the SHBG lowering allele (34) or 0.3/1-s.d. genetically increased natural log transformed SHBG (21).
Using the GWAS data from DIAGRAMv3 or MAGIC (13, 14), associations of the genetic variants influencing these biomarkers with T2D and glycaemic traits were tested. A nominally significance (P<0.05) was reached for the genetic variants that affect bilirubin (P=0.03 for glucose and Homeostatic model assessment (HOMA)-insulin resistance (IR)), NT proBNP (P=0.03 for T2D), D6D activity (P=0.003 for T2D; P=2.7×10−8 for glucose) and dimethylglycine (P=0.004 for glucose) in relation to T2D or glycaemic traits (Table 2). For adiponectin, a nominally significant association with T2D or glycaemic traits was observed for seven out of 19 genetic variants. All of these seven variants, except rs12637534, were non-ADIPOQ adiponectin genetic variants. In line with previous Mendelian randomization studies, the overall effect of non-ADIPOQ variants on T2D or glycaemic traits together with null association of ADIPOQ variants is compatible with pleiotropic effects of adiponectin on T2D (42).
Nominally significance for each biomarker-associated genetic variants with glucose, HOMA-IR and the risk of type 2 diabetes.
P value* | |||||
---|---|---|---|---|---|
Biomarkers | Genetic variants | T2DM13 | Glucose14 | HOMA-IR14 | Nominally significance |
Adiponectin | rs1736653 | – | – | – | NA |
rs2980879 | 0.96 | 0.96 | 0.93 | Null | |
rs16861209 | 0.89 | 0.85 | 0.51 | Null | |
rs601339 | 0.89 | 0.83 | 0.03 | Yes | |
rs17300539 | 0.87 | 0.81 | 0.15 | Null | |
rs7955516 | 0.86 | 0.84 | 0.83 | Null | |
rs12922394 | 0.54 | 0.07 | 0.19 | Null | |
rs3821799 | 0.51 | 0.19 | 0.80 | Null | |
rs731839 | 0.41 | 0.34 | 0.01 | Yes | |
rs3774261 | 0.32 | 0.44 | 0.76 | Null | |
rs17366568 | 0.31 | 0.97 | 0.40 | Null | |
rs7133378 | 0.16 | 0.21 | 0.33 | Null | |
rs12637534 | 0.15 | 0.38 | 0.00 | Yes | |
rs6810075 | 0.1 | 0.77 | 0.52 | Null | |
rs998584 | 0.06 | 0.81 | 0.68 | Null | |
rs1108842 | 0.032 | 0.47 | 0.05 | Yes | |
rs3001032 | 0.029 | 0.27 | 0.00 | Yes | |
rs1597466 | 0.009 | 0.06 | 0.47 | Yes | |
rs2925979 | 0.002 | 0.05 | 0.18 | Yes | |
Beta-carotene | rs6564851 | 0.22 | 0.39 | 0.27 | Null |
Bilirubin | rs6742078 | 0.1 | 0.02 | 0.03 | Yes |
NT proBNP | rs198389 | 0.03 | 0.37 | 0.39 | Yes |
CRP | rs1205 | 0.66 | 0.68 | 0.64 | Null |
rs3093077 | 0.048 | 0.95 | 0.39 | Yes | |
rs1130864 | – | 0.12 | 0.77 | Null | |
D6-desaturase | rs174546 | 0.003 | 2.7×10−8 | 0.89 | Yes |
Dimethylglycine | rs2431332 | 0.51 | 0.004 | 0.29 | Yes |
Ferritin/TMPRSS6 | rs855791 | 1 | 0.09 | 0.25 | Null |
rs4820268 | 0.73 | 0.15 | 0.13 | Null | |
Fetuin-a | rs4917 | 0.08 | 0.11 | 0.67 | Null |
rs2248690 | 0.055 | 0.96 | 0.77 | Null | |
Homocysteine | rs4646356 | 0.88 | 0.43 | 0.35 | Null |
rs12134663 | 0.78 | 0.69 | 0.62 | Null | |
rs1801131 | 0.22 | 0.61 | 0.24 | Null | |
rs1801133 | 0.10 | 0.46 | 0.43 | Null | |
IL-1Ra | rs1542176 | 0.64 | 0.008 | 0.71 | Yes |
rs6743376 | – | 0.91 | 0.46 | Null | |
LTL | rs12880583 | – | – | – | NA |
rs2228041 | – | – | – | NA | |
rs34368910 | – | – | – | NA | |
rs938886 | 0.88 | 0.042 | 0.74 | Yes | |
rs4975605 | 0.66 | 0.45 | 0.61 | Null | |
rs4888444 | 0.09 | 0.91 | 0.34 | Null | |
Lp(a) | rs10455872 | 0.36 | 0.42 | 0.32 | Null |
MIF | rs1007888 | 0.39 | 0.22 | 0.94 | Null |
miRNAs | rs12610873 | 0.033 | 0.80 | 0.71 | Yes |
rs531564 | 0.66 | 0.53 | 0.97 | Null | |
rs2910164 | 0.52 | 0.01 | 0.36 | Yes | |
Resistin | rs3745367 | 0.77 | 0.99 | 0.58 | Null |
rs1423096 | 0.72 | 0.81 | 0.34 | Null | |
SHBG | rs12150660 | 0.18 | 0.52 | 0.24 | Null |
rs6259 | 0.94 | 0.42 | 0.35 | Null | |
Triglycerides | rs12285095 | 0.41 | 0.97 | 0.62 | Null |
rs7557067 | 0.9 | 0.78 | 0.15 | Null | |
rs7679 | 0.88 | 0.90 | 0.29 | Null | |
rs2954029 | 0.72 | 0.64 | 0.87 | Null | |
rs662799 | 0.28 | 0.04 | 0.55 | Yes | |
rs714052 | 0.28 | 0.81 | 0.80 | Null | |
rs10889353 | 0.19 | 0.42 | 0.21 | Null | |
rs7819412 | 0.08 | 0.80 | 0.01 | Yes | |
rs328 | 0.03 | 0.05 | 0.04 | Yes | |
rs17216525 | 0.00 | 0.98 | 0.81 | Yes | |
Uric acid | rs1165151 | 0.74 | 0.42 | 0.86 | Null |
rs2231142 | 0.73 | 0.83 | 0.36 | Null | |
rs505802 | 0.62 | 0.60 | 0.62 | Null | |
rs734553 | 0.57 | 0.20 | 0.72 | Null | |
rs742132 | 0.25 | 0.50 | 0.91 | Null | |
rs12356193 | 0.14 | 0.95 | 0.47 | Null | |
rs12129861 | 0.08 | 0.94 | 0.50 | Null | |
rs17300741 | 0.03 | 0.34 | 0.36 | Yes | |
Vitamin-D | rs10877012 | – | – | – | NA |
rs3755967 | 0.74 | 0.85 | 0.19 | Null | |
rs1155563 | 0.93 | 0.77 | 0.06 | Null | |
rs2298850 | 0.66 | 0.89 | 0.15 | Null | |
rs3829251 | 0.54 | 0.09 | 0.26 | Null | |
rs3794060 | 0.34 | 0.28 | 0.10 | Null | |
rs6599638 | 0.30 | 0.17 | 0.43 | Null | |
rs10741657 | 0.23 | 0.37 | 0.08 | Null | |
rs12785878 | 0.14 | 0.32 | 0.17 | Null | |
rs17217119 | 0.06 | 0.36 | 0.21 | Null | |
rs6013897 | 0.06 | 0.33 | 0.12 | Null | |
Vitamin D-BP | rs2282679 | 0.76 | 1.00 | 0.12 | Null |
*P values for each of the biomarker variants were extracted from publicly available meta-analyses of genome-wide association studies (13, 14). NA, not applicable; CRP, C-reactive protein; BNP, brain natriuretic peptide; interleukin 1 receptor antagonist (IL1-Ra); Lp(a), Lipoprotein(a); LTL, leukocyte telomere length; MIF, macrophage migration inhibitory factor; NT proBNP, N-terminal pro B-type natriuretic peptide; SHBG, sex hormone binding globulin, Vitamin D-BP, Vitamin D binding protein.
Discussion
This literature search of Mendelian randomization studies shows here that ten out of 21 identified biomarkers were reported to be causally associated with T2D. In particular, the presence of potential causal associations (defined as nominally significant) between the biomarker-related variants and T2D and/or glycaemic traits can be confirmed for four biomarkers using the publically available GWAS data. The inconsistency between the identified studies and the summary data from DIAGRAMv3 or MAGIC can be to some extent explained by different design and populations applying across studies, the possibility of false positive associations, selection bias, the heterogeneity in T2D, the use of varied sources to ascertain T2D cases, differences in data quality control and pre-analysis preparations (6, 11).
Taken together, these findings support that the oxidative stress system (bilirubin or metabolites of bilirubin), the brain natriuretic peptide (BNP) hormone system, D6D activity and dimethylglycine may contribute to the development of diabetes or insulin resistance through secondary effects (i.e., pleiotropy) or direct mechanisms (11, 37, 42). Bilirubin that is the major end-product of heme catabolism has antioxidant properties and may compensate the oxidative stress (11, 48). Oxidative stress has been shown as an important factor in the pathophysiology of diabetes (11). At the cellular level, bilirubin can be oxidized to its precursor, biliverdin, to detoxify the excess of oxidants. Biliverdin is rapidly recycled to bilirubin via the action of biliverdin reductase, generating a physiologic cytoprotective cycle in several tissues (49). The underlying mechanism of a protective role of BNP in the aetiology of T2D is unknown in humans (37). In mice, overexpressed BNP signalling cascade can protect against diet-induced insulin resistance and obesity through promoting muscle mitochondrial biogenesis and fat oxidation (37, 50). The biological mechanisms of the relationship between D6D activity and T2D are not well understood (30). Although data of human experimental studies are scarce, observational studies have shown that D6D activity or lifestyle-induced changes in D6D activity was associated with insulin resistance (30, 51, 52). Because D6D catalyses the synthesis of fatty acids, one can speculate that the link between D6D activity and T2D is likely to be mediated by changes in fatty acid composition, which in turn may affect insulin signalling and receptor-binding affinities (30). Dimethylglycine is metabolized to glycine by dimethylglycine dehydrogenase (DMGDH) in mammals (32). Accordingly, a recent GWAS identified that the DMGDH genetic variants were strongly associated with blood-based dimethylglycine (53). Epidemiological studies have observed an inverse association between the precursor of dimethylglycine, betaine, and metabolic risk factors (54) but a positive association of elevated glycine with increased insulin sensitivity (32). In humans, cardiometabolic effects of the inhibition of DMHDH or the supplementation of dimethylglycine have not been investigated (32). However, experimental animal studies have suggested a protective role of dimethylglycine in glucose metabolism through a reduction in DMGDH function (32, 55, 56).
In the post-omics era, epidemiological studies basically suggest that the levels of a given biomarker differ between patients with T2D (or the individuals at high risk for developing diabetes) and individuals without diabetes (2, 57). If the biomarker is not causally related to the disease outcome, the process of developing diabetes may cause the increase or decrease in the levels of the biomarker, as one of the T2D consequences, called reverse causality (4, 26, 37). Unmeasured confounding or measured confounding factors with errors (for example, physical activity by self-report) is another explanation for the observed associations between most biomarkers and T2D. In this context, the use of genetic data (as in Mendelian randomization) can enhance the likelihood that the association between a biomarker and T2D is causal or not (4, 36, 37). The potential role of biomarkers in the development of T2D and the trajectories of glycaemic traits using longitudinal analysis remains to be further confirmed (6, 11). The latter analytical approach can provide insight into the potential value of biomarkers, which indicates pathobiological signals of metabolic changes in the aetiology of T2D. T2D is influenced by the interactions of multiple genes or a gene may have been mapped to multiple biomarkers rather than the biomarker of interest (4, 13, 14). Thus, an extensive knowledge of gene function and biological processes and that the genome interacts with environmental factors is needed to better understand how genetic variations in the human genome contributes to lifelong risk of T2D (4, 11, 13, 14, 57).
In this review, most studies only investigated a single biomarker-diabetes association, and statistically significant associations are reported more often. Here, publication bias should be considered. Other limitations include the lack of complementary analyses in the genetic associations for glycaemic traits and that set of genetic variants or large-scale GWAS for biomarkers were scarce. Similarly, the Mendelian randomization approach using GWAS summary statistics can be extended as a secondary analysis of several datasets that have data on the biomarker-associated variants (or their proxies) for several diseases or traits. This multidisciplinary research requires that a large group of consortia are contacted to obtain summary association statistics from GWAS consortia for outcomes of interest. Moreover, for the biomarkers linked to T2D, underlying biological mechanisms remain unknown. To uncover the underlying mechanisms, one needs to perform a complementary strand of experimental research. Before that, an in silico functional gene network (pathway) analysis can be used to speculate on the possible biological mechanisms of biomarker-disease associations (58, 59). Finally, Mendelian randomization cannot completely control for the possibility of developmental compensation, called canalization, confounding and pleiotropic effects (5, 6, 8, 11, 37).
In conclusion, this is the first study that updates evidence for causal associations between biomarkers and T2D to date. Several Mendelian randomization studies investigated the nature of associations of biomarkers with T2D. Most biomarkers were investigated in a single cohort study or population. However, there were only a few biomarkers that may have causal effects on T2D. Further research is warranted to broadly evaluate the causal effects of multiple biomarkers on T2D and glycaemic traits using data from large-scale cohorts or combined GWAS including many different genetic variants. This genetic approach may advance our understanding of the causes of T2D and potentially enable us to explore novel targets for the prevention and treatment of diabetes.
Declaration of interest
The author declares that there is no conflict of interest that could be perceived as prejudicing the impartiality of the research reported.
Funding
This work was supported by The Netherlands Organization for Scientific Research project (NWO) and the Medical Research Council UK (grant number MC_UU_12015/1). A Abbasi is supported by a Rubicon grant from the NWO (Project no. 825.13.004).
References
- 1↑
Knowler WC, Barrett-Connor E, Fowler SE, Hamman RF, Lachin JM, Walker EA, Nathan DM, Diabetes Prevention Program Research G. Reduction in the incidence of type 2 diabetes with lifestyle intervention or metformin. New England Journal of Medicine 2002 346 393–403. (doi:10.1056/NEJMoa012512).
- 2
Abbasi A, Stolk RP, Bakker SJ. Identification of relevant biomarkers for type 2 diabetes. Lancet. Diabetes & Endocrinology 2014 2 106–107. (doi:10.1016/S2213-8587(14)70005-1).
- 3
Lyons TJ, Basu A. Biomarkers in diabetes: hemoglobin A1c, vascular and tissue markers. Translational Research 2012 159 303–312. (doi:10.1016/j.trsl.2012.01.009).
- 4↑
Sandhu MS, Debenham SL, Barroso I, Loos RJ. Mendelian randomisation studies of type 2 diabetes: future prospects. Diabetologia 2008 51 211–213. (doi:10.1007/s00125-007-0903-x).
- 5
Nelson CP, Hamby SE, Saleheen D, Hopewell JC, Zeng L, Assimes TL, Kanoni S, Willenborg C, Burgess S, Amouyel P et al.. Genetically Determined Height and Coronary Artery Disease. New England Journal of Medicine 2015 372 1608–1618. (doi:10.1056/NEJMoa1404881).
- 6
VanderWeele TJ, Tchetgen Tchetgen EJ, Cornelis M, Kraft P. Methodological challenges in mendelian randomization. Epidemiology 2014 25 427–435. (doi:10.1097/EDE.0000000000000081).
- 7
Nikolopoulos GK, Bagos PG, Tsangaris I, Tsiara CG, Kopterides P, Vaiopoulos A, Kapsimali V, Bonovas S, Tsantes AE. The association between plasminogen activator inhibitor type 1 (PAI-1) levels, PAI-1 4G/5G polymorphism, and myocardial infarction: a Mendelian randomization meta-analysis. Clinical Chemistry and Laboratory Medicine 2014 52 937–950. (doi:10.1515/cclm-2013-1124).
- 8
Haase CL, Tybjaerg-Hansen A, Qayyum AA, Schou J, Nordestgaard BG, Frikke-Schmidt R. LCAT, HDL cholesterol and ischemic cardiovascular disease: a Mendelian randomization study of HDL cholesterol in 54,500 individuals. Journal of Clinical Endocrinology and Metabolism 2012 97 E248–E256. (doi:10.1210/jc.2011-1846).
- 9
Smith GD, Ebrahim S. Mendelian randomization: prospects, potentials, and limitations. International Journal of Epidemiology 2004 33 30–42. (doi:10.1093/ije/dyh132).
- 10
Burgess S, Scott RA, Timpson NJ, Davey Smith G, Thompson SG, EPIC- InterAct Consortium. Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. European Journal of Epidemiology 2015 30 543–552. (doi:10.1007/s10654-015-0011-z).
- 11↑
Abbasi A, Deetman PE, Corpeleijn E, Gansevoort RT, Gans RO, Hillege HL, van der Harst P, Stolk RP, Navis G, Alizadeh BZ et al.. Bilirubin as a potential causal factor in type 2 diabetes risk: a mendelian randomization study. Diabetes 2015 64 1459–1469. (doi:10.2337/db14-0228).
- 12↑
Liberati A, Altman DG, Tetzlaff J, Mulrow C, Gotzsche PC, Ioannidis JP, Clarke M, Devereaux PJ, Kleijnen J, Moher D et al.. The PRISMA statement for reporting systematic reviews and meta-analyses of studies that evaluate healthcare interventions: explanation and elaboration. British Medical Journal 2009 339 b2700. (doi:10.1136/bmj.b2700).
- 13↑
Morris AP, Voight BF, Teslovich TM, Ferreira T, Segre AV, Steinthorsdottir V, Strawbridge RJ, Khan H, Grallert H, Mahajan A et al.. Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nature Genetics 2012 44 981–990. (doi:10.1038/ng.2383).
- 14↑
Dupuis J, Langenberg C, Prokopenko I, Saxena R, Soranzo N, Jackson AU, Wheeler E, Glazer NL, Bouatia-Naji N, Gloyn AL et al.. New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nature Genetics 2010 42 105–116. (doi:10.1038/ng.520).
- 15
Besseling J, Kastelein JJ, Defesche JC, Hutten BA, Hovingh GK. Association between familial hypercholesterolemia and prevalence of type 2 diabetes mellitus. Journal of the American Medical Association 2015 313 1029–1036. (doi:10.1001/jama.2015.1206).
- 16↑
Brunner EJ, Kivimaki M, Witte DR, Lawlor DA, Davey Smith G, Cooper JA, Miller M, Lowe GD, Rumley A, Casas JP et al.. Inflammation, insulin resistance, and diabetes--Mendelian randomization using CRP haplotypes points upstream. PLoS Medicine 2008 5 e155. (doi:10.1371/journal.pmed.0050155).
- 17
Buijsse B, Boeing H, Hirche F, Weikert C, Schulze MB, Gottschald M, Kuhn T, Katzke VA, Teucher B, Dierkes J et al.. Plasma 25-hydroxyvitamin D and its genetic determinants in relation to incident type 2 diabetes: a prospective case-cohort study. European Journal of Epidemiology 2013 28 743–752. (doi:10.1007/s10654-013-9844-5).
- 18↑
Chung CM, Lin TH, Chen JW, Leu HB, Yin WH, Ho HY, Sheu SH, Tsai WC, Chen JH, Lin SJ et al.. Common quantitative trait locus downstream of RETN gene identified by genome-wide association study is associated with risk of type 2 diabetes mellitus in Han Chinese: a Mendelian randomization effect. Diabetes/Metabolism Research and Reviews 2014 30 232–240. (doi:10.1002/dmrr.2481).
- 19↑
Dastani Z, Hivert MF, Timpson N, Perry JR, Yuan X, Scott RA, Henneman P, Heid IM, Kizer JR, Lyytikainen LP et al.. Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals. PLoS Genetics 2012 8 e1002607. (doi:10.1371/journal.pgen.1002607).
- 20↑
De Silva NM, Freathy RM, Palmer TM, Donnelly LA, Luan J, Gaunt T, Langenberg C, Weedon MN, Shields B, Knight BA et al.. Mendelian randomization studies do not support a role for raised circulating triglyceride levels influencing type 2 diabetes, glucose levels, or insulin resistance. Diabetes 2011 60 1008–1018. (doi:10.2337/db10-1317).
- 21↑
Ding EL, Song Y, Manson JE, Hunter DJ, Lee CC, Rifai N, Buring JE, Gaziano JM, Liu S. Sex hormone-binding globulin and risk of type 2 diabetes in women and men. New England Journal of Medicine 2009 361 1152–1163. (doi:10.1056/NEJMoa0804381).
- 22↑
Gan W, Guan Y, Wu Q, An P, Zhu J, Lu L, Jing L, Yu Y, Ruan S, Xie D et al.. Association of TMPRSS6 polymorphisms with ferritin, hemoglobin, and type 2 diabetes risk in a Chinese Han population. American Journal of Clinical Nutrition 2012 95 626–632. (doi:10.3945/ajcn.111.025684).
- 23↑
Herder C, Klopp N, Baumert J, Muller M, Khuseyinova N, Meisinger C, Martin S, Illig T, Koenig W, Thorand B. Effect of macrophage migration inhibitory factor (MIF) gene variants and MIF serum concentrations on the risk of type 2 diabetes: results from the MONICA/KORA Augsburg Case-Cohort Study, 1984-2002. Diabetologia 2008 51 276–284. (doi:10.1007/s00125-007-0800-3).
- 24↑
Huang T, Ren J, Huang J, Li D. Association of homocysteine with type 2 diabetes: a meta-analysis implementing Mendelian randomization approach. BMC Genomics 2013 14 867. (doi:10.1186/1471-2164-14-867).
- 25↑
Huang T, Sun J, Chen Y, Xie H, Xu D, Li D. Associations of common variants in methionine metabolism pathway genes with plasma homocysteine and the risk of type 2 diabetes in Han Chinese. Journal of Nutrigenetics and Nutrigenomics 2014 7 63–74. (doi:10.1159/000365007).
- 26↑
Interleukin 1 Genetics C Cardiometabolic effects of genetic upregulation of the interleukin 1 receptor antagonist: a Mendelian randomisation analysis. Lancet. Diabetes & Endocrinology 2015 3 243–253. (doi:10.1016/S2213-8587(15)00034-0).
- 27↑
Jensen MK, Bartz TM, Djousse L, Kizer JR, Zieman SJ, Rimm EB, Siscovick DS, Psaty BM, Ix JH, Mukamal KJ. Genetically elevated fetuin-A levels, fasting glucose levels, and risk of type 2 diabetes: the cardiovascular health study. Diabetes Care 2013 36 3121–3127. (doi:10.2337/dc12-2323).
- 28
Jorde R, Schirmer H, Wilsgaard T, Joakimsen RM, Mathiesen EB, Njolstad I, Lochen ML, Figenschau Y, Berg JP, Svartberg J et al.. Polymorphisms related to the serum 25-hydroxyvitamin D level and risk of myocardial infarction, diabetes, cancer and mortality. The Tromso Study. PLoS ONE 2012 7 e37295. (doi:10.1371/journal.pone.0037295).
- 29
Kamstrup PR, Nordestgaard BG. Lipoprotein(a) concentrations, isoform size, and risk of type 2 diabetes: a Mendelian randomisation study. Lancet. Diabetes & Endocrinology 2013 1 220–227. (doi:10.1016/S2213-8587(13)70064-0).
- 30↑
Kroger J, Zietemann V, Enzenbach C, Weikert C, Jansen EH, Doring F, Joost HG, Boeing H, Schulze MB. Erythrocyte membrane phospholipid fatty acids, desaturase activity, and dietary fatty acids in relation to risk of type 2 diabetes in the European Prospective Investigation into Cancer and Nutrition (EPIC)-Potsdam Study. American Journal of Clinical Nutrition 2011 93 127–142. (doi:10.3945/ajcn.110.005447).
- 31↑
Leong A, Rehman W, Dastani Z, Greenwood C, Timpson N, Langsetmo L, Berger C, Metastroke, Fu L, Wong BY et al.. The causal effect of vitamin D binding protein (DBP) levels on calcemic and cardiometabolic diseases: a Mendelian randomization study. PLoS Medicine 2014 11 e1001751. (doi:10.1371/journal.pmed.1001751).
- 32↑
Magnusson M, Wang TJ, Clish C, Engstrom G, Nilsson P, Gerszten RE, Melander O. Dimethylglycine deficiency and the development of diabetes mellitus. Diabetes 2015 64 3010–3016. (doi:10.2337/db14-1863).
- 33↑
Perry JR, Ferrucci L, Bandinelli S, Guralnik J, Semba RD, Rice N, Melzer D, Consortium D, Saxena R, Scott LJ et al.. Circulating beta-carotene levels and type 2 diabetes-cause or effect? Diabetologia 2009 52 2117–2121. (doi:10.1007/s00125-009-1475-8).
- 34↑
Perry JR, Weedon MN, Langenberg C, Jackson AU, Lyssenko V, Sparso T, Thorleifsson G, Grallert H, Ferrucci L, Maggio M et al.. Genetic evidence that raised sex hormone binding globulin (SHBG) levels reduce the risk of type 2 diabetes. Human Molecular Genetics 2010 19 535–544. (doi:10.1093/hmg/ddp522).
- 35↑
Peters KE, Beilby J, Cadby G, Warrington NM, Bruce DG, Davis WA, Davis TM, Wiltshire S, Knuiman M, McQuillan BM et al.. A comprehensive investigation of variants in genes encoding adiponectin (ADIPOQ) and its receptors (ADIPOR1/R2), and their association with serum adiponectin, type 2 diabetes, insulin resistance and the metabolic syndrome. BMC Medical Genetics 2013 14 15. (doi:10.1186/1471-2350-14-15).
- 36↑
Pfister R, Barnes D, Luben R, Forouhi NG, Bochud M, Khaw KT, Wareham NJ, Langenberg C. No evidence for a causal link between uric acid and type 2 diabetes: a Mendelian randomisation approach. Diabetologia 2011 54 2561–2569. (doi:10.1007/s00125-011-2235-0).
- 37↑
Pfister R, Sharp S, Luben R, Welsh P, Barroso I, Salomaa V, Meirhaeghe A, Khaw KT, Sattar N, Langenberg C et al.. Mendelian randomization study of B-type natriuretic peptide and type 2 diabetes: evidence of causal association from population studies. PLoS Medicine 2011 8 e1001112. (doi:10.1371/journal.pmed.1001112).
- 38
Song Y, Yeung E, Liu A, Vanderweele TJ, Chen L, Lu C, Liu C, Schisterman EF, Ning Y, Zhang C. Pancreatic β-cell function and type 2 diabetes risk: quantify the causal effect using a Mendelian randomization approach based on meta-analyses. Human Molecular Genetics 2012 21 5010–5018. (doi:10.1093/hmg/dds339).
- 39
Strawbridge RJ, Dupuis J, Prokopenko I, Barker A, Ahlqvist E, Rybin D, Petrie JR, Travers ME, Bouatia-Naji N, Dimas AS et al.. Genome-wide association identifies nine common variants associated with fasting proinsulin levels and provides new insights into the pathophysiology of type 2 diabetes. Diabetes 2011 60 2624–2634. (doi:10.2337/db11-0415).
- 40
Trombetta M, Bonetti S, Boselli ML, Miccoli R, Trabetti E, Malerba G, Pignatti PF, Bonora E, Del Prato S, Bonadonna RC. PPARG2 Pro12Ala and ADAMTS9 rs4607103 as "insulin resistance loci" and "insulin secretion loci" in Italian individuals. The GENFIEV study and the Verona Newly Diagnosed Type 2 Diabetes Study (VNDS) 4. Acta Diabetologica 2013 50 401–408. (doi:10.1007/s00592-012-0443-9).
- 41↑
Wang TT, Chen YJ, Sun LL, Zhang SJ, Zhou ZY, Qiao H. Affection of single-nucleotide polymorphisms in miR-27a, miR-124a, and miR-146a on susceptibility to type 2 diabetes mellitus in Chinese Han people. Chinese Medical Journal 2015 128 533–539. (doi:10.4103/0366-6999.151112).
- 42↑
Yaghootkar H, Lamina C, Scott RA, Dastani Z, Hivert MF, Warren LL, Stancakova A, Buxbaum SG, Lyytikainen LP, Henneman P et al.. Mendelian randomization studies do not support a causal role for reduced circulating adiponectin levels in insulin resistance and type 2 diabetes. Diabetes 2013 62 3589–3598. (doi:10.2337/db13-0128).
- 43
Yaghootkar H, Scott RA, White CC, Zhang W, Speliotes E, Munroe PB, Ehret GB, Bis JC, Fox CS, Walker M et al.. Genetic evidence for a normal-weight "metabolically obese" phenotype linking insulin resistance, hypertension, coronary artery disease, and type 2 diabetes. Diabetes 2014 63 4369–4377. (doi:10.2337/db14-0318).
- 44
Ye Z, Haycock PC, Gurdasani D, Pomilla C, Boekholdt SM, Tsimikas S, Khaw KT, Wareham NJ, Sandhu MS, Forouhi NG. The association between circulating lipoprotein(a) and type 2 diabetes: is it causal? Diabetes 2014 63 332–342. (doi:10.2337/db13-1144).
- 45
Ye Z, Sharp SJ, Burgess S, Scott RA, Imamura F, InterAct C, Langenberg C, Wareham NJ, Forouhi NG. Association between circulating 25-hydroxyvitamin D and incident type 2 diabetes: a mendelian randomisation study. Lancet. Diabetes & Endocrinology 2015 3 35–42. (doi:10.1016/S2213-8587(14)70184-6).
- 46↑
You NC, Chen BH, Song Y, Lu X, Chen Y, Manson JE, Kang M, Howard BV, Margolis KL, Curb JD et al.. A prospective study of leukocyte telomere length and risk of type 2 diabetes in postmenopausal women. Diabetes 2012 61 2998–3004. (doi:10.2337/db12-0241).
- 47↑
Voight BF, Scott LJ, Steinthorsdottir V, Morris AP, Dina C, Welch RP, Zeggini E, Huth C, Aulchenko YS, Thorleifsson G et al.. Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis. Nature Genetics 2010 42 579–589. (doi:10.1038/ng.609).
- 48
Van Campenhout A, Van Campenhout C, Lagrou AR, Abrams P, Moorkens G, Van Gaal L, Manuel-y-Keenoy B. Impact of diabetes mellitus on the relationships between iron-, inflammatory- and oxidative stress status. Diabetes/Metabolism Research and Reviews 2006 22 444–454. (doi:10.1002/dmrr.635).
- 49↑
Sedlak TW, Snyder SH. Messenger molecules and cell death: therapeutic implications. Journal of the American Medical Association 2006 295 81–89. (doi:10.1001/jama.295.1.81).
- 50
Miyashita K, Itoh H, Tsujimoto H, Tamura N, Fukunaga Y, Sone M, Yamahara K, Taura D, Inuzuka M, Sonoyama T et al.. Natriuretic peptides/cGMP/cGMP-dependent protein kinase cascades promote muscle mitochondrial biogenesis and prevent obesity. Diabetes 2009 58 2880–2892. (doi:10.2337/db09-0393).
- 51
Corpeleijn E, Feskens EJ, Jansen EH, Mensink M, Saris WH, de Bruin TW, Blaak EE. Improvements in glucose tolerance and insulin sensitivity after lifestyle intervention are related to changes in serum fatty acid profile and desaturase activities: the SLIM study. Diabetologia 2006 49 2392–2401. (doi:10.1007/s00125-006-0383-4).
- 52
Vessby B, Gustafsson IB, Tengblad S, Boberg M, Andersson A. Desaturation and elongation of Fatty acids and insulin action. Annals of the New York Academy of Sciences 2002 967 183–195. (doi:10.1111/j.1749-6632.2002.tb04275.x).
- 53↑
Rhee EP, Ho JE, Chen MH, Shen D, Cheng S, Larson MG, Ghorbani A, Shi X, Helenius IT, O'Donnell CJ et al.. A genome-wide association study of the human metabolome in a community-based cohort. Cell Metabolism 2013 18 130–143. (doi:10.1016/j.cmet.2013.06.013).
- 54↑
Konstantinova SV, Tell GS, Vollset SE, Nygard O, Bleie O, Ueland PM. Divergent associations of plasma choline and betaine with components of metabolic syndrome in middle age and elderly men and women. Journal of Nutrition 2008 138 914–920.
- 55
Ho JE, Larson MG, Vasan RS, Ghorbani A, Cheng S, Rhee EP, Florez JC, Clish CB, Gerszten RE, Wang TJ. Metabolite profiles during oral glucose challenge. Diabetes 2013 62 2689–2698. (doi:10.2337/db12-0754).
- 56
Sheikh K, Camejo G, Lanne B, Halvarsson T, Landergren MR, Oakes ND. Beyond lipids, pharmacological PPARalpha activation has important effects on amino acid metabolism as studied in the rat. American Journal of Physiology. Endocrinology and Metabolism 2007 292 E1157–E1165. (doi:10.1152/ajpendo.00254.2006).
- 57
Ahluwalia TS, Allin KH, Sandholt CH, Sparso TH, Jorgensen ME, Rowe M, Christensen C, Brandslund I, Lauritzen T, Linneberg A et al.. Discovery of coding genetic variants influencing diabetes-related serum biomarkers and their impact on risk of type 2 diabetes. Journal of Clinical Endocrinology and Metabolism 2015 100 E664–E671. (doi:10.1210/jc.2014-3677).
- 58
Kraja AT, Chasman DI, North KE, Reiner AP, Yanek LR, Kilpelainen TO, Smith JA, Dehghan A, Dupuis J, Johnson AD et al.. Pleiotropic genes for metabolic syndrome and inflammation. Molecular Genetics and Metabolism 2014 112 317–338. (doi:10.1016/j.ymgme.2014.04.007).
- 59
Vaez A, Jansen R, Prins BP, Hottenga JJ, de Geus EJ, Boomsma DI, Penninx BW, Nolte IM, Snieder H, Alizadeh BZ. In Silico Post Genome-Wide Association Studies Analysis of C-Reactive Protein Loci Suggests an Important Role for Interferons. Circulation. Cardiovascular Genetics 2015 8 487–497. (doi:10.1161/CIRCGENETICS.114.000714).