Effects of endocrine disrupting chemicals and their interactions with genetic risk scores on cardiometabolic traits

Ubiquitous non-persistent endocrine disrupting chemicals (EDCs) have inconsistent associations with cardiometabolic traits. Additionally, large-scale genome-wide association studies (GWASs) have yielded many genetic risk variants for cardiometabolic traits and diseases. This study aimed to investigate the associations between a wide range of EDC exposures (parabens, bisphenols, and phthalates) and 14 cardiometabolic traits and whether these are moderated by their respective genetic risk scores (GRSs). Data were from 1074 participants aged 18 years or older of the Lifelines Cohort Study, a large population-based biobank. GRSs for 14 cardiometabolic traits were calculated based on genome-wide significant common variants from recent GWASs. The concentrations of 15 EDCs in 24-hour urine were measured by isotope dilution liquid chromatography tandem mass spectrometry technology. The main effects of trait-specific GRSs and each of the EDC exposures and their interaction effects on the 14 cardiometabolic traits were examined in multiple linear regression. The present study confirmed significant main effects for all GRSs on their corresponding cardiometabolic trait. Regarding the main effects of EDC exposures, 26 out of 280 EDC-trait tests were significant with explained variances ranging from 0.43 % (MMP- estimated glomerular filtration rate (eGFR)) to 2.37 % (PrP-waist-hip ratio adjusted body mass index (WHRadjBMI)). We confirmed the association of MiBP and MBzP with WHRadjBMI and body mass index (BMI), and showed that parabens, bisphenol F, and many other phthalate metabolites significantly contributed to the variance of WHRadjBMI, BMI, high-density lipoprotein (HDL), eGFR, fasting glucose (FG), and diastolic blood pressure (DBP). Only one association between BMI and bisphenol F was nominally significantly moderated by the GRS explaining 0.36 % of the variance. However, it did not survive multiple testing correction. We showed that non-persistent EDC exposures exerted effects on BMI, WHRadjBMI, HDL, eGFR, FG, and DBP. However no evidence for a modulating role of GRSs was found.

A genome-wide association study of 24-hour urinary excretion of endocrine disrupting chemicals

Ubiquitous exposure to environmental endocrine disrupting chemicals (EDCs) instigates a major public health problem, but much remains unknown on the inter-individual differences in metabolism and excretion of EDCs. To examine this we performed a two-stage genome-wide association study (GWAS) for 24-hour urinary excretions of four parabens, two bisphenols, and nine phthalate metabolites. Results showed five genome-wide significant (p-value < 5x10-8) and replicated single nucleotide polymorphisms (SNPs) representing four independent signals that associated with mono-(2-ethyl-5-carboxypentyl) phthalate (MECPP) and mono-(2-ethyl-5-hydroxyhexyl) phthalate (MEHHP). Three of the four signals were located on chromosome 10 in a locus harboring the cytochrome P450 (CYP) genes CYP2C9, CYP2C58P, and CYP2C19 (rs117529685, pMECPP = 5.38x10-25; rs117033379, pMECPP = 1.96x10-19; rs4918798, pMECPP = 4.01x10-71; rs7895726, pMEHHP = 1.37x10-15, r2 with rs4918798 = 0.93). The other signal was on chromosome 6 close to the solute carrier (SLC) genes SLC17A1, SLC17A3, SLC17A4, and SCGN (rs1359232, pMECPP = 7.6x10-16). These four SNPs explained a substantial part (8.3 % - 9.2 %) of the variance in MECPP in the replication cohort. Bioinformatics analyses supported a likely causal role of CYP2C9 and SLC17A1 in metabolism and excretion of MECPP and MEHHP. Our results provide biological insights into mechanisms of phthalate metabolism and excretion with a likely causal role for CYP2C9 and SLC17A1.

Urinary concentrations of bisphenols and parabens and their association with attention, hyperactivity and impulsivity at adolescence

BACKGROUND

Neurobehavioural disorder diagnoses have been increasing over the last decades, leading to heightened interest in the aetiological factors involved. Endocrine disrupting chemicals, such as parabens and bisphenols, have been suggested as one of those factors. It is unknown whether exposure during adolescence may affect neurobehavioural development.

OBJECTIVE

To determine whether urinary concentrations of parabens and bisphenols are associated with attention and concentration in adolescents, in general and sex-specific.

METHODS

We invited 188 adolescents (13-15 years old) for the follow-up birth cohort-study. Concentrations of five parabens and three bisphenols (BPA; BPF; BPS) were measured in morning urine after overnight fasting, using a validated LC-MS/MS method. Attention and concentration were assessed at the clinic with subtests of the Test of Everyday Attention in Children and the Dutch Attention Deficit Hyperactivity Disorder questionnaire (AVL), the latter being filled in by parents. Linear regression analyses were performed, adjusting for urine creatinine concentrations and potential confounding factors.

RESULTS

101 (54%) adolescents participated (46 girls; 55 boys). Urinary paraben concentrations were higher in girls than in boys. Methylparaben was positively associated with attention in girls (p ≤ .05; B= -2.836; 95%CI= -5.175;-.497), ethylparaben negatively with hyperactivity (p ≤ .05; B= -1.864; 95%CI= -3.587;-.141). Butylparaben was associated with more optimal scores on parent reported attention. Propylparaben was negatively associated with scores on sustained auditory attention in girls (p ≤ .10; B=.444; 95%CI= -.009;.896). Bisphenol concentrations were not associated with scores on attention and concentration after adjusting for confounders.

CONCLUSION

In 13-15-year-old Dutch adolescents, urinary concentrations of methylparaben and ethylparaben were associated with better attention and less hyperactivity, whereas a trend toward significance was found between higher urinary propylparaben concentrations and poorer attention. Bisphenol concentrations were not associated with attention and concentration after adjusting for confounders.

Genome-wide association analyses of physical activity and sedentary behavior provide insights into underlying mechanisms and roles in disease prevention

Although physical activity and sedentary behavior are moderately heritable, little is known about the mechanisms that influence these traits. Combining data for up to 703,901 individuals from 51 studies in a multi-ancestry meta-analysis of genome-wide association studies yields 99 loci that associate with self-reported moderate-to-vigorous intensity physical activity during leisure time (MVPA), leisure screen time (LST) and/or sedentary behavior at work. Loci associated with LST are enriched for genes whose expression in skeletal muscle is altered by resistance training. A missense variant in ACTN3 makes the alpha-actinin-3 filaments more flexible, resulting in lower maximal force in isolated type IIA muscle fibers, and possibly protection from exercise-induced muscle damage. Finally, Mendelian randomization analyses show that beneficial effects of lower LST and higher MVPA on several risk factors and diseases are mediated or confounded by body mass index (BMI). Our results provide insights into physical activity mechanisms and its role in disease prevention.

Epigenome-wide association study of incident type 2 diabetes: a meta-analysis of five prospective European cohorts

AIMS/HYPOTHESIS

Type 2 diabetes is a complex metabolic disease with increasing prevalence worldwide. Improving the prediction of incident type 2 diabetes using epigenetic markers could help tailor prevention efforts to those at the highest risk. The aim of this study was to identify predictive methylation markers for incident type 2 diabetes by combining epigenome-wide association study (EWAS) results from five prospective European cohorts.

METHODS

We conducted a meta-analysis of EWASs in blood collected 7-10 years prior to type 2 diabetes diagnosis. DNA methylation was measured with Illumina Infinium Methylation arrays. A total of 1250 cases and 1950 controls from five longitudinal cohorts were included: Doetinchem, ESTHER, KORA1, KORA2 and EPIC-Norfolk. Associations between DNA methylation and incident type 2 diabetes were examined using robust linear regression with adjustment for potential confounders. Inverse-variance fixed-effects meta-analysis of cohort-level individual CpG EWAS estimates was performed using METAL. The methylGSA R package was used for gene set enrichment analysis. Confirmation of genome-wide significant CpG sites was performed in a cohort of Indian Asians (LOLIPOP, UK).

RESULTS

The meta-analysis identified 76 CpG sites that were differentially methylated in individuals with incident type 2 diabetes compared with control individuals (p values <1.1 × 10-7). Sixty-four out of 76 (84.2%) CpG sites were confirmed by directionally consistent effects and p values <0.05 in an independent cohort of Indian Asians. However, on adjustment for baseline BMI only four CpG sites remained genome-wide significant, and addition of the 76 CpG methylation risk score to a prediction model including established predictors of type 2 diabetes (age, sex, BMI and HbA1c) showed no improvement (AUC 0.757 vs 0.753). Gene set enrichment analysis of the full epigenome-wide results clearly showed enrichment of processes linked to insulin signalling, lipid homeostasis and inflammation.

CONCLUSIONS/INTERPRETATION

By combining results from five European cohorts, and thus significantly increasing study sample size, we identified 76 CpG sites associated with incident type 2 diabetes. Replication of 64 CpGs in an independent cohort of Indian Asians suggests that the association between DNA methylation levels and incident type 2 diabetes is robust and independent of ethnicity. Our data also indicate that BMI partly explains the association between DNA methylation and incident type 2 diabetes. Further studies are required to elucidate the underlying biological mechanisms and to determine potential causal roles of the differentially methylated CpG sites in type 2 diabetes development.

Epigenome-wide association study of serum urate reveals insights into urate co-regulation and the SLC2A9 locus

Elevated serum urate levels, a complex trait and major risk factor for incident gout, are correlated with cardiometabolic traits via incompletely understood mechanisms. DNA methylation in whole blood captures genetic and environmental influences and is assessed in transethnic meta-analysis of epigenome-wide association studies (EWAS) of serum urate (discovery, n = 12,474, replication, n = 5522). The 100 replicated, epigenome-wide significant (p < 1.1E-7) CpGs explain 11.6% of the serum urate variance. At SLC2A9, the serum urate locus with the largest effect in genome-wide association studies (GWAS), five CpGs are associated with SLC2A9 gene expression. Four CpGs at SLC2A9 have significant causal effects on serum urate levels and/or gout, and two of these partly mediate the effects of urate-associated GWAS variants. In other genes, including SLC7A11 and PHGDH, 17 urate-associated CpGs are associated with conditions defining metabolic syndrome, suggesting that these CpGs may represent a blood DNA methylation signature of cardiometabolic risk factors. This study demonstrates that EWAS can provide new insights into GWAS loci and the correlation of serum urate with other complex traits.

Meta-analyses identify DNA methylation associated with kidney function and damage

Chronic kidney disease is a major public health burden. Elevated urinary albumin-to-creatinine ratio is a measure of kidney damage, and used to diagnose and stage chronic kidney disease. To extend the knowledge on regulatory mechanisms related to kidney function and disease, we conducted a blood-based epigenome-wide association study for estimated glomerular filtration rate (n = 33,605) and urinary albumin-to-creatinine ratio (n = 15,068) and detected 69 and seven CpG sites where DNA methylation was associated with the respective trait. The majority of these findings showed directionally consistent associations with the respective clinical outcomes chronic kidney disease and moderately increased albuminuria. Associations of DNA methylation with kidney function, such as CpGs at JAZF1, PELI1 and CHD2 were validated in kidney tissue. Methylation at PHRF1, LDB2, CSRNP1 and IRF5 indicated causal effects on kidney function. Enrichment analyses revealed pathways related to hemostasis and blood cell migration for estimated glomerular filtration rate, and immune cell activation and response for urinary albumin-to-creatinineratio-associated CpGs.

Temporal exposure and consistency of endocrine disrupting chemicals in a longitudinal study of individuals with impaired fasting glucose

Endocrine disrupting chemicals (EDCs) include non-persistent exogenous substances such as parabens, bisphenols and phthalates which have been associated with a range of metabolic disorders and disease. It is unclear if exposure remains consistent over time. We investigated change in indicators of EDC exposure between 2009 and 2016 and assessed its consistency between and within individuals over a median follow-up time of 47 months in a sample of Dutch individuals. Of 500 Dutch individuals, two 24 h urine samples were analysed for 5 parabens, 3 bisphenols and 13 metabolites of in total 8 different phthalates. We calculated per-year differences using meta-analysis and assessed temporal correlations between and within individuals using Spearman correlation coefficients, intra-class correlation coefficients (ICC) and kappa-statistics. We found a secular decrease in concentrations of methyl, ethyl, propyl and n-butyl paraben, bisphenol A, and metabolites of di-ethyl phthalate (DEP), di-butyl phthalate (DBP), di-(2-ethyl-hexyl) phthalate (DEHP), and butylbenzyl phthalate (DBzP) which varied from 8 to 96% (ethyl paraben, propyl paraben) between 2009 and 2016. Within-person temporal correlations were highest for parabens (ICC: 0.34 to 0.40) and poorest for bisphenols (ICC: 0.15 to 0.23). For phthalate metabolites, correlations decreased most between time periods (ICC < 48 months: 0.22 to 0.39; ≥48 months: 0.05 to 0.32). When categorizing EDC concentrations, 33-54% of individuals remained in the lowest or highest category and temporal correlations were similar to continuous measurements. Exposure to most EDCs decreased between 2009 and 2016 in a sample of individuals with impaired fasting glucose from the Dutch population. Temporal consistency was generally poor. The inconsistency in disease associations may be influenced by individual-level or temporal variation exhibited by EDCs. Our findings call for the need for repeated measurements of EDCs in observational studies before and during at-risk temporal windows for the disease.

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Endocrine disrupting chemicals during diet-induced weight loss - A post-hoc analysis of the LOWER study

The link between exposure to endocrine disrupting chemicals (EDCs) and the rapid increase in prevalence of obesity has recently been suggested. However, the magnitude and health impact of EDC exposure in at-risk populations remain largely unclear. In this study, we investigated the effect of a dietary intervention driven reduction in adipose tissue on the magnitude of urinary EDC exposure and mobilization, and whether higher EDC exposure leads to impaired weight loss in obese individuals. In this post-hoc analysis of the Lifestyle, OverWeight, Energy Restriction (LOWER) study from the Netherlands, 218 subjects were included. Five parabens, three bisphenols and thirteen metabolites of eight phthalates were measured in 24-h urine using LC-MS/MS, before and after three-months of a calory-restricted weight reduction intervention program. Associations between adiposity-related traits and EDCs were tested using multivariable linear regression and linear mixed effects models. A multiple testing correction based on the false discovery rate (FDR) was applied. After the 3-month intervention, urinary paraben and bisphenol excretions remained similar. Excretions of mono-butyl phthalates and most high-molecular-weight phthalates decreased, whereas mono-ethyl phthalate increased (all FDR<0.05). A reduction in adipose tissue was not associated with higher urinary EDC excretions. Higher baseline EDC excretions were associated with higher post-intervention body-mass index (methyl-, propylparaben), waist circumference (propylparaben, mono-n-butyl phthalate, mono-benzyl phthalate), and body fat percentage (mono-ethyl phthalate, mono-benzyl phthalate). Associations between parabens and body-mass index, and mono-benzyl phthalate and waist circumference and body fat percentage remained after multiple testing correction (all FDR<0.05). In a study of obese participants, we observed a reduction in most phthalates after a weight reduction intervention. A reduction in adipose tissue may not lead to mobilization and successively to higher urinary EDC excretions. Higher baseline paraben and phthalate exposures were associated with reduced weight loss, suggesting obesogenic properties.

Genomic analysis of diet composition finds novel loci and associations with health and lifestyle

We conducted genome-wide association studies (GWAS) of relative intake from the macronutrients fat, protein, carbohydrates, and sugar in over 235,000 individuals of European ancestries. We identified 21 unique, approximately independent lead SNPs. Fourteen lead SNPs are uniquely associated with one macronutrient at genome-wide significance (P < 5 × 10^-8), while five of the 21 lead SNPs reach suggestive significance (P < 1 × 10^-5) for at least one other macronutrient. While the phenotypes are genetically correlated, each phenotype carries a partially unique genetic architecture. Relative protein intake exhibits the strongest relationships with poor health, including positive genetic associations with obesity, type 2 diabetes, and heart disease (r_g ≈ 0.15-0.5). In contrast, relative carbohydrate and sugar intake have negative genetic correlations with waist circumference, waist-hip ratio, and neighborhood deprivation (|r_g| ≈ 0.1-0.3) and positive genetic correlations with physical activity (r_g ≈ 0.1 and 0.2). Relative fat intake has no consistent pattern of genetic correlations with poor health but has a negative genetic correlation with educational attainment (r_g ≈-0.1). Although our analyses do not allow us to draw causal conclusions, we find no evidence of negative health consequences associated with relative carbohydrate, sugar, or fat intake. However, our results are consistent with the hypothesis that relative protein intake plays a role in the etiology of metabolic dysfunction.

An epigenome-wide association study identifies multiple DNA methylation markers of exposure to endocrine disruptors

BACKGROUND

Exposure to environmental endocrine disrupting chemicals (EDCs) may play an important role in the epidemic of metabolic diseases. Epigenetic alterations may functionally link EDCs with gene expression and metabolic traits.

OBJECTIVES

We aimed to evaluate metabolic-related effects of the exposure to endocrine disruptors including five parabens, three bisphenols, and 13 metabolites of nine phthalates as measured in 24-hour urine on epigenome-wide DNA methylation.

METHODS

A blood-based epigenome-wide association study was performed in 622 participants from the Lifelines DEEP cohort using Illumina Infinium HumanMethylation450 methylation data and EDC excretions in 24-hour urine. Out of the 21 EDCs, 13 compounds were detected in >75% of the samples and, together with bisphenol F, were included in these analyses. Furthermore, we explored the putative function of identified methylation markers and their correlations with metabolic traits.

RESULTS

We found 20 differentially methylated cytosine-phosphate-guanines (CpGs) associated with 10 EDCs at suggestive p-value < 1 × 10^-6, of which four, associated with MEHP and MEHHP, were genome-wide significant (Bonferroni-corrected p-value < 1.19 × 10^-7). Nine out of 20 CpGs were significantly associated with at least one of the tested metabolic traits, such as fasting glucose, glycated hemoglobin, blood lipids, and/or blood pressure. 18 out of 20 EDC-associated CpGs were annotated to genes functionally related to metabolic syndrome, hypertension, obesity, type 2 diabetes, insulin resistance and glycemic traits.

CONCLUSIONS

The identified DNA methylation markers for exposure to the most common EDCs provide suggestive mechanism underlying the contributions of EDCs to metabolic health. Follow-up studies are needed to unravel the causality of EDC-induced methylation changes in metabolic alterations.

Identification, Heritability, and Relation With Gene Expression of Novel DNA Methylation Loci for Blood Pressure

We conducted an epigenome-wide association study meta-analysis on blood pressure (BP) in 4820 individuals of European and African ancestry aged 14 to 69. Genome-wide DNA methylation data from peripheral leukocytes were obtained using the Infinium Human Methylation 450k BeadChip. The epigenome-wide association study meta-analysis identified 39 BP-related CpG sites with P<1×10^-5. In silico replication in the CHARGE consortium of 17 010 individuals validated 16 of these CpG sites. Out of the 16 CpG sites, 13 showed novel association with BP. Conversely, out of the 126 CpG sites identified as being associated (P<1×10^-7) with BP in the CHARGE consortium, 21 were replicated in the current study. Methylation levels of all the 34 CpG sites that were cross-validated by the current study and the CHARGE consortium were heritable and 6 showed association with gene expression. Furthermore, 9 CpG sites also showed association with BP with P<0.05 and consistent direction of the effect in the meta-analysis of the Finnish Twin Cohort (199 twin pairs and 4 singletons; 61% monozygous) and the Netherlands Twin Register (266 twin pairs and 62 singletons; 84% monozygous). Bivariate quantitative genetic modeling of the twin data showed that a majority of the phenotypic correlations between methylation levels of these CpG sites and BP could be explained by shared unique environmental rather than genetic factors, with 100% of the correlations of systolic BP with cg19693031 (TXNIP) and cg00716257 (JDP2) determined by environmental effects acting on both systolic BP and methylation levels.

Development and Interlaboratory Validation of Two Fast UPLC-MS-MS Methods Determining Urinary Bisphenols, Parabens and Phthalates

People are constantly exposed to a wide variety of chemicals. Some of these compounds, such as parabens, bisphenols and phthalates, are known to have endocrine disrupting potencies. Over the years, these endocrine disrupting chemicals (EDCs) have been a rising cause for concern. In this study, we describe setup and validation of two methods to measure EDCs in human urine, using ultra-performance liquid chromatography tandem mass spectrometry. The phenol method determines methyl-, ethyl-, propyl-, n-butyl- and benzylparaben and bisphenol A, F and S. The phthalate method determines in total 13 metabolites of dimethyl, diethyl, diisobutyl, di-n-butyl, di(2-ethylhexyl), butylbenzyl, diiso-nonyl and diisodecyl phthalate. Runtime was 7 and 8 min per sample for phenols and phthalates, respectively. The methods were validated by the National Institute of Standards & Technology (NIST) for 13 compounds. In addition, EDCs were measured in forty 24-h urine samples, of which 12 EDCs were compared with the same samples measured in an established facility (Rigshospitalet, Copenhagen, Denmark). The intra-assay coefficient of variability (CV) was highest at 10% and inter-assay CV was highest at 12%. Recoveries ranged from 86 to 115%. The limit of detection ranged from 0.06 to 0.43 ng/mL. Of 21 compounds, 10 were detected above limit of detection in ≥93% of the samples. Eight compounds were in accordance to NIST reference concentrations. Differences in intercept were found for two compounds whereas slope differed for six compounds between our method and that used in the Danish facility. In conclusion, we set up and validated two high-throughput methods with very short runtime capable of measuring 5 parabens, 3 bisphenols and 13 different metabolites of 8 phthalates. Sensitivity of the phenol method was increased by using ammonium fluoride in the mobile phase.

The effects of bariatric surgery on clinical profile, DNA methylation, and ageing in severely obese patients

BACKGROUND

Severe obesity is a growing, worldwide burden and conventional therapies including radical change of diet and/or increased physical activity have limited results. Bariatric surgery has been proposed as an alternative therapy showing promising results. It leads to substantial weight loss and improvement of comorbidities such as type 2 diabetes. Increased adiposity is associated with changes in epigenetic profile, including DNA methylation. We investigated the effect of bariatric surgery on clinical profile, DNA methylation, and biological age estimated using Horvath's epigenetic clock.

RESULTS

To determine the impact of bariatric surgery and subsequent weight loss on clinical traits, a cohort of 40 severely obese individuals (BMI = 30-73 kg/m²) was examined at the time of surgery and at three follow-up visits, i.e., 3, 6, and 12 months after surgery. The majority of the individuals were women (65%) and the mean age at surgery was 45.1 ± 8.1 years. We observed a significant decrease over time in BMI, fasting glucose, HbA1c, HOMA-IR, insulin, total cholesterol, triglycerides, LDL and free fatty acids levels, and a significant small increase in HDL levels (all p values < 0.05). Epigenome-wide association analysis revealed 4857 differentially methylated CpG sites 12 months after surgery (at Bonferroni-corrected p value < 1.09 × 10^-7). Including BMI change in the model decreased the number of significantly differentially methylated CpG sites by 51%. Gene set enrichment analysis identified overrepresentation of multiple processes including regulation of transcription, RNA metabolic, and biosynthetic processes in the cell. Bariatric surgery in severely obese patients resulted in a decrease in both biological age and epigenetic age acceleration (EAA) (mean = - 0.92, p value = 0.039).

CONCLUSIONS

Our study shows that bariatric surgery leads to substantial BMI decrease and improvement of clinical outcomes observed 12 months after surgery. These changes explained part of the association between bariatric surgery and DNA methylation. We also observed a small, but significant improvement of biological age. These epigenetic changes may be modifiable by environmental lifestyle factors and could be used as potential biomarkers for obesity and in the future for obesity related comorbidities.

Mendelian randomisation analyses find pulmonary factors mediate the effect of height on coronary artery disease

There is evidence that lower height is associated with a higher risk of coronary artery disease (CAD) and increased risk of type 2 diabetes (T2D). It is not clear though whether these associations are causal, direct or mediated by other factors. Here we show that one standard deviation higher genetically determined height (~6.5 cm) is causally associated with a 16% decrease in CAD risk (OR = 0.84, 95% CI 0.80-0.87). This causal association remains after performing sensitivity analyses relaxing pleiotropy assumptions. The causal effect of height on CAD risk is reduced by 1-3% after adjustment for potential mediators (lipids, blood pressure, glycaemic traits, body mass index, socio-economic status). In contrast, our data suggest that lung function (measured by forced expiratory volume [FEV1] and forced vital capacity [FVC]) is a mediator of the effect of height on CAD. We observe no direct causal effect of height on the risk of T2D.

Skin autofluorescence predicts incident type 2 diabetes, cardiovascular disease and mortality in the general population

AIMS/HYPOTHESIS

Earlier studies have shown that skin autofluorescence measured with an AGE reader estimates the accumulation of AGEs in the skin, which increases with ageing and is associated with the metabolic syndrome and type 2 diabetes. In the present study, we examined whether the measurement of skin autofluorescence can predict 4 year risk of incident type 2 diabetes, cardiovascular disease (CVD) and mortality in the general population.

METHODS

For this prospective analysis, we included 72,880 participants of the Dutch Lifelines Cohort Study, who underwent baseline investigations between 2007 and 2013, had validated baseline skin autofluorescence values available and were not known to have diabetes or CVD. Individuals were diagnosed with incident type 2 diabetes by self-report or by a fasting blood glucose ≥7.0 mmol/l or HbA_1c ≥48 mmol/mol (≥6.5%) at follow-up. Participants were diagnosed as having incident CVD (myocardial infarction, coronary interventions, cerebrovascular accident, transient ischaemic attack, intermittent claudication or vascular surgery) by self-report. Mortality was ascertained using the Municipal Personal Records Database.

RESULTS

After a median follow-up of 4 years (range 0.5-10 years), 1056 participants (1.4%) had developed type 2 diabetes, 1258 individuals (1.7%) were diagnosed with CVD, while 928 (1.3%) had died. Baseline skin autofluorescence was elevated in participants with incident type 2 diabetes and/or CVD and in those who had died (all p < 0.001), compared with individuals who survived and remained free of the two diseases. Skin autofluorescence predicted the development of type 2 diabetes, CVD and mortality, independent of several traditional risk factors, such as the metabolic syndrome, glucose and HbA_1c.

CONCLUSIONS/INTERPRETATION

The non-invasive skin autofluorescence measurement is of clinical value for screening for future risk of type 2 diabetes, CVD and mortality, independent of glycaemic measures and the metabolic syndrome.

Dietary patterns and physical activity in the metabolically (un)healthy obese: the Dutch Lifelines cohort study

BACKGROUND

Diversity in the reported prevalence of metabolically healthy obesity (MHO), suggests that modifiable factors may be at play. We evaluated differences in dietary patterns and physical activity between MHO and metabolically unhealthy obesity (MUO).

METHODS

Cross-sectional data of 9270 obese individuals (30-69 years) of the Lifelines Cohort Study was used. MHO was defined as obesity and no metabolic syndrome risk factors and no cardiovascular disease history. MUO was defined as obesity and ≥2 metabolic syndrome risk factors. Sex-specific associations of dietary patterns (identified by principal component analysis) and physical activity with MHO were assessed by multivariable logistic regression (reference group: MUO). Analyses were adjusted for multiple covariates.

RESULTS

Among 3442 men and 5828 women, 10.2% and 24.4% had MHO and 56.9% and 35.3% MUO, respectively. We generated four obesity-specific dietary patterns. Two were related to MHO, and in women only. In the highest quartile (Q) of 'bread, potatoes and sweet snacks' pattern, odds ratio (OR) (95% CI) for MHO was 0.52 (0.39-0.70). For the healthier pattern 'fruit, vegetables and fish', an OR of 1.36 (1.09-1.71) in Q3 and 1.55 (1.21-1.97) in Q4 was found for MHO. For physical activity, there was a positive association between moderate physical activity and vigorous physical activity in the highest tertile and MHO in women and men, respectively (OR 1.19 (1.01-1.41) and OR 2.02 (1.50-2.71)).

CONCLUSION

The healthier diet -characterized by 'fruit, vegetables and fish'- and moderate physical activity in women, and vigorous physical activity in men may be related to MHO. The (refined) carbohydrate-rich 'bread, potatoes and sweet snacks' dietary pattern was found to counteract MHO in women.

DNA methylation markers associated with type 2 diabetes, fasting glucose and HbA1c levels: a systematic review and replication in a case-control sample of the Lifelines study

AIMS/HYPOTHESIS

Epigenetic mechanisms may play an important role in the aetiology of type 2 diabetes. Recent epigenome-wide association studies (EWASs) identified several DNA methylation markers associated with type 2 diabetes, fasting glucose and HbA_1c levels. Here we present a systematic review of these studies and attempt to replicate the CpG sites (CpGs) with the most significant associations from these EWASs in a case-control sample of the Lifelines study.

METHODS

We performed a systematic literature search in PubMed and EMBASE for EWASs to test the association between DNA methylation and type 2 diabetes and/or glycaemic traits and reviewed the search results. For replication purposes we selected 100 unique CpGs identified in peripheral blood, pancreas, adipose tissue and liver from 15 EWASs, using study-specific Bonferroni-corrected significance thresholds. Methylation data (Illumina 450K array) in whole blood from 100 type 2 diabetic individuals and 100 control individuals from the Lifelines study were available. Multivariate linear models were used to examine the associations of the specific CpGs with type 2 diabetes and glycaemic traits.

RESULTS

From the 52 CpGs identified in blood and selected for replication, 15 CpGs showed nominally significant associations with type 2 diabetes in the Lifelines sample (p < 0.05). The results for five CpGs (in ABCG1, LOXL2, TXNIP, SLC1A5 and SREBF1) remained significant after a stringent multiple-testing correction (changes in methylation from -3% up to 3.6%, p < 0.0009). All associations were directionally consistent with the original EWAS results. None of the selected CpGs from the tissue-specific EWASs were replicated in our methylation data from whole blood. We were also unable to replicate any of the CpGs associated with HbA_1c levels in the healthy control individuals of our sample, while two CpGs (in ABCG1 and CCDC57) for fasting glucose were replicated at a nominal significance level (p < 0.05).

CONCLUSIONS/INTERPRETATION

A number of differentially methylated CpGs reported to be associated with type 2 diabetes in the EWAS literature were replicated in blood and show promise for clinical use as disease biomarkers. However, more prospective studies are needed to support the robustness of these findings.

Thyroid function and metabolic syndrome in the population-based LifeLines cohort study

BACKGROUND

The metabolic syndrome (MetS) is a combination of unfavourable health factors which includes abdominal obesity, dyslipidaemia, elevated blood pressure and impaired fasting glucose. Earlier studies have reported a relationship between thyroid function and some MetS components or suggested that serum free thyroxine (FT4) or free triiodothyronine (FT3) levels within the normal range were independently associated with insulin resistance. We assessed how thyroid function relates to MetS prevalence in a large population-based study.

METHODS

Data of 26,719 people of western European descent, aged 18-80 years from the Dutch LifeLines Cohort study, all with normal thyroid stimulating hormone (TSH), FT4 and FT3 levels (electrochemiluminescent immunoassay, Roche Modular E170 Analyzer), were available. MetS was defined with the revised National Cholesterol Education Programs Adults Treatment Panel III (NCEP ATP III) criteria. We calculated prevalence of all MetS components according to TSH, FT4 and FT3 quartiles.

RESULTS

At similar TSH levels and age (mean 45 yrs), men had significantly higher levels of FT4, FT3, blood pressure (BP), heart rate, total and LDL-cholesterol, triglycerides (TG), and creatinine, but lower HDL-cholesterol compared to women (all p < 0.001). In total, 11.8% of women and 20.7% of men had MetS. In men, lower FT4 levels were associated with higher prevalence of MetS and all MetS components. In women, lower FT4 quartile was only associated with a higher prevalence of elevated TG, waist circumference, and MetS. However, when corrected for confounding factors like age, BMI, current smoking and alcohol consumption, a significant relationship was found between FT3 and three MetS components in men, and all five components in women. Moreover, the highest quartiles of FT3 and the FT3FT4 ratio predicted a 49% and 67% higher prevalence of MetS in men, and a 62 and 80% higher prevalence in women.

CONCLUSIONS

When corrected for possible confounding factors, higher plasma levels of FT3 are associated with several components of the MetS. Only in men, lower FT4 is related to MetS. In the highest FT3 and FT3FT4 quartiles, there is a 50-80% increased risk of having MetS compared to the lowest quartile. Further studies are needed to assess the possible causality of this relationship.

Distribution of Non-Persistent Endocrine Disruptors in Two Different Regions of the Human Brain

Non-persistent endocrine disrupting chemicals (npEDCs) can affect multiple organs and systems in the body. Whether npEDCs can accumulate in the human brain is largely unknown. The major aim of this pilot study was to examine the presence of environmental phenols and parabens in two distinct brain regions: the hypothalamus and white-matter tissue. In addition, a potential association between these npEDCs concentrations and obesity was investigated. Post-mortem brain material was obtained from 24 individuals, made up of 12 obese and 12 normal-weight subjects (defined as body mass index (BMI) > 30 and BMI < 25 kg/m², respectively). Nine phenols and seven parabens were measured by isotope dilution TurboFlow-LC-MS/MS. In the hypothalamus, seven suspect npEDCs (bisphenol A, triclosan, triclocarban and methyl-, ethyl-, n-propyl-, and benzyl paraben) were detected, while five npEDCs (bisphenol A, benzophenone-3, triclocarban, methyl-, and n-propyl paraben) were found in the white-matter brain tissue. We observed higher levels of methylparaben (MeP) in the hypothalamic tissue of obese subjects as compared to controls (p = 0.008). Our findings indicate that some suspected npEDCs are able to cross the blood-brain barrier. Whether the presence of npEDCs can adversely affect brain function and to which extent the detected concentrations are physiologically relevant needs to be further investigated.

Genotype-covariate interaction effects and the heritability of adult body mass index

Obesity is a worldwide epidemic, with major health and economic costs. Here we estimate heritability for body mass index (BMI) in 172,000 sibling pairs and 150,832 unrelated individuals and explore the contribution of genotype-covariate interaction effects at common SNP loci. We find evidence for genotype-age interaction (likelihood ratio test (LRT) = 73.58, degrees of freedom (df) = 1, P = 4.83 × 10^-18), which contributed 8.1% (1.4% s.e.) to BMI variation. Across eight self-reported lifestyle factors, including diet and exercise, we find genotype-environment interaction only for smoking behavior (LRT = 19.70, P = 5.03 × 10^-5 and LRT = 30.80, P = 1.42 × 10^-8), which contributed 4.0% (0.8% s.e.) to BMI variation. Bayesian association analysis suggests that BMI is highly polygenic, with 75% of the SNP heritability attributable to loci that each explain <0.01% of the phenotypic variance. Our findings imply that substantially larger sample sizes across ages and lifestyles are required to understand the full genetic architecture of BMI.

The association between various smoking behaviors, cotinine biomarkers and skin autofluorescence, a marker for advanced glycation end product accumulation

BACKGROUND

Skin autofluorescence, a biomarker for advanced glycation end products (AGEs) accumulation, has been shown to predict diabetes-related cardiovascular complications and is associated with several environmental and lifestyle factors. In the present study, we examined the association between various smoking behaviors and skin autofluorescence, as well as the association between several cotinine biomarkers and skin autofluorescence, using both epidemiological and metabolomics data.

METHODS

In a cross-sectional study, we evaluated participants from the LifeLines Cohort Study and the Qatar Metabolomics Study on Diabetes (QMDiab). In the LifeLines Cohort Study smoking behavior and secondhand smoking were assessed in 8,905 individuals including 309 individuals (3.5%) with type 2 diabetes. In QMDiab, cotinine biomarkers were measured in saliva, plasma and urine in 364 individuals of whom 188 (51%) had type 2 diabetes. Skin autofluorescence was measured non-invasively in all participants using the AGE Reader.

RESULTS

Skin autofluorescence levels increased with a higher number of hours being exposed to secondhand smoking. Skin autofluorescence levels of former smokers approached levels of never smokers after around 15 years of smoking cessation. Urinary cotinine N-oxide, a biomarker of nicotine exposure, was found to be positively associated with skin autofluorescence in the QMDiab study (p = 0.03).

CONCLUSIONS

In the present study, we have demonstrated that secondhand smoking is associated with higher skin autofluorescence levels whereas smoking cessation has a beneficial effect on skin autofluorescence. Finally, urinary cotinine N-oxide might be used as an alternative way for questionnaires to examine the effect of (environmental) tobacco smoking on skin autofluorescence.

Sex, BMI and age differences in metabolic syndrome: the Dutch Lifelines Cohort Study

INTRODUCTION

To evaluate the prevalence of metabolic syndrome (MetS) and its individual components within sex-, body mass index (BMI)- and age combined clusters. In addition, we used the age-adjusted blood pressure thresholds to demonstrate the effect on the prevalence of MetS and elevated blood pressure.

SUBJECTS AND METHODS

Cross-sectional data from 74,531 Western European participants, aged 18-79 years, were used from the Dutch Lifelines Cohort Study. MetS was defined according to the revised NCEP-ATPIII. Age-adjusted blood pressure thresholds were defined as recommended by the eight reports of the Joint National Committee (≥140/90 mmHg for those aged <60 years, and ≥150/90 mmHg for those aged ≥60 years).

RESULTS

19.2% men and 12.1% women had MetS. MetS prevalence increased with BMI and age. Independent of BMI, abdominal obesity dominated MetS prevalence especially in women, while elevated blood pressure was already highly prevalent among young men. Applying age-adjusted blood pressure thresholds resulted in a 0.2-11.9% prevalence drop in MetS and 6.0-36.3% prevalence drop in elevated blood pressure, within the combined sex, BMI and age clusters.

CONCLUSIONS

We observed a gender disparity with age and BMI for the prevalence of MetS and, especially, abdominal obesity and elevated blood pressure. The strict threshold level for elevated blood pressure in the revised NCEP-ATPIII, results in an overestimation of MetS prevalence.

Genome-wide physical activity interactions in adiposity - A meta-analysis of 200,452 adults

Physical activity (PA) may modify the genetic effects that give rise to increased risk of obesity. To identify adiposity loci whose effects are modified by PA, we performed genome-wide interaction meta-analyses of BMI and BMI-adjusted waist circumference and waist-hip ratio from up to 200,452 adults of European (n = 180,423) or other ancestry (n = 20,029). We standardized PA by categorizing it into a dichotomous variable where, on average, 23% of participants were categorized as inactive and 77% as physically active. While we replicate the interaction with PA for the strongest known obesity-risk locus in the FTO gene, of which the effect is attenuated by ~30% in physically active individuals compared to inactive individuals, we do not identify additional loci that are sensitive to PA. In additional genome-wide meta-analyses adjusting for PA and interaction with PA, we identify 11 novel adiposity loci, suggesting that accounting for PA or other environmental factors that contribute to variation in adiposity may facilitate gene discovery.

No Effect of the Thr92Ala Polymorphism of Deiodinase-2 on Thyroid Hormone Parameters, Health-Related Quality of Life, and Cognitive Functioning in a Large Population-Based Cohort Study

INTRODUCTION

The presence of the Thr92Ala polymorphism of deiodinase-2 (D2) has been thought to have several effects. It may influence its enzymatic function, is associated with increased expression of genes involved in oxidative stress in brain tissue, and may predict favorable response to combination levothyroxine (LT4) plus triiodothyronine (T3) therapy. It was hypothesized that homozygous carriers of the D2-92Ala allele have different thyroid hormone parameters, and reduced health-related quality of life (HRQoL) and cognitive functioning.

METHODS

In 12,625 participants from the LifeLines cohort study with genome-wide genetic data available, the effects of the Thr92Ala polymorphism (rs225014) were evaluated in the general population and in 364 people treated with thyroid hormone replacement therapy, the latter mainly because of primary hypothyroidism. In addition to evaluating anthropometric data, medication use, and existence of metabolic syndrome, HRQoL was assessed with the RAND 36-Item Health Survey, and the Ruff Figural Fluency Test was used as a sensitive test for executive functioning. Data on thyrotropin, free thyroxine (fT4), and free T3 (fT3) levels were available in a subset of 4479 participants.

RESULTS

The mean age (±standard deviation) was 53 ± 12 years and the body mass index was 27.0 ± 4.5 kg/m2 in the LT4 users compared with 48 ± 11 years and 26.2 ± 4.1 kg/m2 in participants from the general population. The Ala/Ala genotype of the D2-Thr92Ala polymorphism was present in 11.3% of LT4 users and in 10.7% of the general population. In total, 3742/4479 subjects with thyroid hormone data available had normal TSH (0.4-4.0 mIU/L), and 88% of LT4 users were females. LT4 users had higher fT4, lower fT3, and a lower fT3/fT4 ratio, and female patients had lower scores on the HRQoL domains of physical functioning, vitality, mental health, social functioning, bodily pain, and general health compared with those not using LT4 (p < 0.005). Executive functioning scores, as part of cognitive functioning, were comparable between female LT4 users and the general population. In both groups, the D2-Thr92Ala polymorphism was not associated with differences in TSH, fT4, fT3, the fT3/fT4 ratio, presence of metabolic syndrome or other comorbidities, use of medication, HRQoL, and cognitive functioning.

CONCLUSION

The Thr92Ala polymorphism of D2 was not associated with thyroid parameters, HRQoL, and cognitive functioning in the general population and in participants on thyroid hormone replacement therapy.

Skin autofluorescence, a non-invasive biomarker for advanced glycation end products, is associated with the metabolic syndrome and its individual components

BACKGROUND

The metabolic syndrome (MetS) comprises several cardiometabolic risk factors associated with increased risk for both type 2 diabetes and cardiovascular disease. Skin autofluorescence (SAF), a non-invasive biomarker of advanced glycation end products accumulation, is associated with cardiovascular complications in subjects with diabetes. The aim of the present study was to examine the association between SAF and the presence of MetS as well as its individual components in a general population.

METHODS

For this cross-sectional analysis, we included 78,671 non-diabetic subjects between 18 and 80 years of age who participated in the LifeLines Cohort Study and had SAF measurement obtained non-invasively using the AGE Reader. MetS was defined according to the revised NCEP ATP III criteria. Students unpaired t test was used to test differences between groups. Both logistic and linear regression analyses were performed in order to test associations between the individual MetS components and SAF.

RESULTS

Subjects with MetS had higher SAF (2.07 ± 0.45 arbitrary units, AU) compared to individuals without MetS (1.89 ± 0.42 AU) (p < 0.001). There was a positive association between the number of MetS components and higher SAF Z-scores (p < 0.001). Individuals in the highest SAF tertile had a higher presence of MetS (OR 2.61; 95% CI 2.48-2.75) and some of the individual components compared to subjects in the lowest SAF tertile. After correction for age, gender, creatinine clearance, HbA1c and smoking status, only elevated blood pressure and low HDL cholesterol remained significantly associated with higher SAF (p = 0.002 and p = 0.001 respectively).

CONCLUSION

Skin autofluorescence was associated with the presence of MetS and some of its individual components. In addition, increasing SAF Z-scores were observed with a higher number of MetS components. Prospective studies are needed to establish whether SAF can be used as an (additional) screening tool to predict both cardiovascular disease and type 2 diabetes in high-risk populations.

New Locus for Skin Intrinsic Fluorescence in Type 1 Diabetes Also Associated With Blood and Skin Glycated Proteins

Skin fluorescence (SF) noninvasively measures advanced glycation end products (AGEs) in the skin and is a risk indicator for diabetes complications. N-acetyltransferase 2 (NAT2) is the only known locus influencing SF. We aimed to identify additional genetic loci influencing SF in type 1 diabetes (T1D) through a meta-analysis of genome-wide association studies (N = 1,359) including Diabetes Control and Complications Trial/Epidemiology of Diabetes Interventions and Complications (DCCT/EDIC) and Wisconsin Epidemiologic Study of Diabetic Retinopathy (WESDR). A locus on chromosome 1, rs7533564 (P = 1.9 × 10(-9)), was associated with skin intrinsic fluorescence measured by SCOUT DS (excitation 375 nm, emission 435-655 nm), which remained significant after adjustment for time-weighted HbA1c (P = 1.7 × 10(-8)). rs7533564 was associated with mean HbA1c in meta-analysis (P = 0.0225), mean glycated albumin (P = 0.0029), and glyoxal hydroimidazolones (P = 0.049), an AGE measured in skin biopsy collagen, in DCCT. rs7533564 was not associated with diabetes complications in DCCT/EDIC or with SF in subjects without diabetes (nondiabetic [ND]) (N = 8,721). In conclusion, we identified a new locus associated with SF in T1D subjects that did not show similar effect in ND subjects, suggesting a diabetes-specific effect. This association needs to be investigated in type 2 diabetes.

Correction: The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

Lifestyle and clinical determinants of skin autofluorescence in a population-based cohort study

BACKGROUND

Skin autofluorescence (SAF) is a noninvasive marker of advanced glycation end products (AGEs). In diabetes, higher SAF levels have been positively associated with long-term complications, cardiovascular morbidity and mortality. Because little is known about the factors that influence SAF in nondiabetic individuals, we assessed the association of clinical and lifestyle parameters with SAF as well as their interactions in a large-scale, nondiabetic population and performed the same analysis in a type 2 diabetic subgroup.

METHODS

In a cross-sectional study in participants from the LifeLines Cohort Study, extensive clinical and biochemical phenotyping, including SAF measurement, was assessed in 9009 subjects of whom 314 (3·5%) subjects with type 2 diabetes.

RESULTS

Mean SAF was 2·04 ± 0·44 arbitrary units (AU) in nondiabetic individuals and 2·44 ± 0·55 AU in type 2 diabetic subjects (P < 0·0001). Multivariate backward regression analysis showed that in the nondiabetic population, SAF was significantly and independently associated with age, BMI, HbA1c, creatinine clearance, genetic polymorphism in NAT2 (rs4921914), current smoking, pack-years of smoking and coffee consumption. In the type 2 diabetic group, a similar set of factors was associated with SAF, except for coffee consumption.

CONCLUSIONS

In addition to the established literature on type 2 diabetes, we have demonstrated that SAF levels are associated with several clinical and lifestyle factors in the nondiabetic population. These parameters should be taken into consideration when using SAF as a screening or prediction tool for populations at risk for cardiovascular disease and diabetes.

Biomonitoring of human exposures to chlorinated derivatives and structural analogs of bisphenol A

The high reactivity of bisphenol A (BPA) with disinfectant chlorine is evident in the instantaneous formation of chlorinated BPA derivatives (ClxBPA) in various environmental media that show increased estrogen-activity when compared with that of BPA. The documented health risks associated with BPA exposures have led to the gradual market entry of BPA structural analogs, such as bisphenol S (BPS), bisphenol F (BPF), bisphenol B (BPB), etc. A suite of exposure sources to ClxBPA and BPA analogs in the domestic environment is anticipated to drive the nature and range of halogenated BPA derivatives that can form when residual BPA comes in contact with disinfectant in tap water and/or consumer products. The primary objective of this review was to survey all available studies reporting biomonitoring protocols of ClxBPA and structural BPA analogs (BPS, BPF, BPB, etc.) in human matrices. Focus was paid on describing the analytical methodologies practiced for the analysis of ClxBPA and BPA analogs using hyphenated chromatography and mass spectrometry techniques, because current methodologies for human matrices are complex. During the last decade, an increasing number of ecotoxicological, cell-culture and animal-based and human studies dealing with ClxBPA exposure sources and routes of exposure, metabolism and toxicity have been published. Up to date findings indicated the association of ClxBPA with metabolic conditions, such as obesity, lipid accumulation, and type 2 diabetes mellitus, particularly in in-vitro and in-vivo studies. We critically discuss the limitations, research needs and future opportunities linked with the inclusion of ClxBPA and BPA analogs into exposure assessment protocols of relevant epidemiological studies.

Fine mapping the CETP region reveals a common intronic insertion associated to HDL-C

BACKGROUND

Individuals with exceptional longevity and their offspring have significantly larger high-density lipoprotein concentrations (HDL-C) particle sizes due to the increased homozygosity for the I405V variant in the cholesteryl ester transfer protein (CETP) gene. In this study, we investigate the association of CETP and HDL-C further to identify novel, independent CETP variants associated with HDL-C in humans.

METHODS

We performed a meta-analysis of HDL-C within the CETP region using 59,432 individuals imputed with 1000 Genomes data. We performed replication in an independent sample of 47,866 individuals and validation was done by Sanger sequencing.

RESULTS

The meta-analysis of HDL-C within the CETP region identified five independent variants, including an exonic variant and a common intronic insertion. We replicated these 5 variants significantly in an independent sample of 47,866 individuals. Sanger sequencing of the insertion within a single family confirmed segregation of this variant. The strongest reported association between HDL-C and CETP variants, was rs3764261; however, after conditioning on the five novel variants we identified the support for rs3764261 was highly reduced (β_unadjusted=3.179 mg/dl (P value=5.25×10^-509), β_adjusted=0.859 mg/dl (P value=9.51×10^-25)), and this finding suggests that these five novel variants may partly explain the association of CETP with HDL-C. Indeed, three of the five novel variants (rs34065661, rs5817082, rs7499892) are independent of rs3764261.

CONCLUSIONS

The causal variants in CETP that account for the association with HDL-C remain unknown. We used studies imputed to the 1000 Genomes reference panel for fine mapping of the CETP region. We identified and validated five variants within this region that may partly account for the association of the known variant (rs3764261), as well as other sources of genetic contribution to HDL-C.

Genome-wide genetic homogeneity between sexes and populations for human height and body mass index

Sex-specific genetic effects have been proposed to be an important source of variation for human complex traits. Here we use two distinct genome-wide methods to estimate the autosomal genetic correlation (rg) between men and women for human height and body mass index (BMI), using individual-level (n = ∼44 000) and summary-level (n = ∼133 000) data from genome-wide association studies. Results are consistent and show that the between-sex genetic correlation is not significantly different from unity for both traits. In contrast, we find evidence of genetic heterogeneity between sexes for waist-hip ratio (rg = ∼0.7) and between populations for BMI (rg = ∼0.9 between Europe and the USA) but not for height. The lack of evidence for substantial genetic heterogeneity for body size is consistent with empirical findings across traits and species.

Health-Related Quality of Life in Relation to Obesity Grade, Type 2 Diabetes, Metabolic Syndrome and Inflammation

Background

Health-related quality of life (HR-QoL) may be compromised in obese individuals, depending on the presence of other complications. The aim of this study is to assess the effect of obesity-related conditions on HR-QoL. These conditions are i) grade of obesity with and without type 2 diabetes (T2D), ii) metabolic syndrome (MetS), and iii) level of inflammation.

Methods

From the Dutch LifeLines Cohort Study we included 13,686 obese individuals, aged 18–80 years. HR-QoL was measured with the RAND 36-Item Health Survey which encompasses eight health domains. We calculated the percentage of obese individuals with poor HR-QoL, i.e. those scoring below the domain and sex specific cut-off value derived from the normal weight population. Logistic regression analysis was used to calculate the probability of having poor domain scores according to the conditions under study.

Results

Higher grades of obesity and the additional presence of T2D were associated with lower HR-QoL, particularly in the domains physical functioning (men: odds ratios (ORs) 1.48–11.34, P<0.005, and women: ORs 1.66–5.05, P<0.001) and general health (men: ORs 1.44–3.07, P<0.005, and women: ORs 1.36–3.73, P<0.001). A higher percentage of obese individuals with MetS had a poor HR-QoL than those without MetS. Furthermore, we observed a linear trend between inflammation and the percentage of obese individuals with poor scores on the HR-QoL domains. Individuals with MetS were more likely to have poor scores in the domains general health, vitality, social functioning and role limitations due to emotional problems. Obese women with increased inflammation levels were more likely to have poor scores on all domains except role limitations due to emotional problems and mental health.

Conclusions

The impact of obesity on an individual’s quality of life is enhanced by grade of obesity, T2D, MetS and inflammation and are mainly related to reduced physical health. The mental well-being is less often impaired.

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHR_adjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHR_adjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHR_adjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHR_adjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHR_adjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHR_adjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Adult body size and body shape differ substantially between men and women and change over time. More than 100 genetic variants that influence body mass index (measure of body size) or waist-to-hip ratio (measure of body shape) have been identified. While there is evidence that some genetic loci affect body shape differently in men than in women, little is known about whether genetic effects differ in older compared to younger adults, and whether such changes differ between men and women. Therefore, we conducted a systematic genome-wide search, including 114 studies (>320,000 individuals), to specifically identify genetic loci with age- and or sex-dependent effects on body size and shape. We identified 15 loci of which the effect on BMI was different in older compared to younger adults, whereas we found no evidence for loci with different effects in men compared to women. The opposite was seen for body shape as we identified 44 loci of which the effect on waist-to-hip ratio differed between men and women, but no difference between younger and older adults were observed. Our observations may provide new insights into the biology that underlies weight change with age or the sexual dimorphism of body shape.

Genetic and epigenetic regulation of gene expression in fetal and adult human livers

Background

The liver plays a central role in the maintenance of homeostasis and health in general. However, there is substantial inter-individual variation in hepatic gene expression, and although numerous genetic factors have been identified, less is known about the epigenetic factors.

Results

By analyzing the methylomes and transcriptomes of 14 fetal and 181 adult livers, we identified 657 differentially methylated genes with adult-specific expression, these genes were enriched for transcription factor binding sites of HNF1A and HNF4A. We also identified 1,000 genes specific to fetal liver, which were enriched for GATA1, STAT5A, STAT5B and YY1 binding sites. We saw strong liver-specific effects of single nucleotide polymorphisms on both methylation levels (28,447 unique CpG sites (meQTL)) and gene expression levels (526 unique genes (eQTL)), at a false discovery rate (FDR) < 0.05. Of the 526 unique eQTL associated genes, 293 correlated significantly not only with genetic variation but also with methylation levels. The tissue-specificities of these associations were analyzed in muscle, subcutaneous adipose tissue and visceral adipose tissue. We observed that meQTL were more stable between tissues than eQTL and a very strong tissue-specificity for the identified associations between CpG methylation and gene expression.

Conclusions

Our analyses generated a comprehensive resource of factors involved in the regulation of hepatic gene expression, and allowed us to estimate the proportion of variation in gene expression that could be attributed to genetic and epigenetic variation, both crucial to understanding differences in drug response and the etiology of liver diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-860) contains supplementary material, which is available to authorized users.

GWAS identifies an NAT2 acetylator status tag single nucleotide polymorphism to be a major locus for skin fluorescence

AIMS/HYPOTHESIS

Skin fluorescence (SF) is a non-invasive marker of AGEs and is associated with the long-term complications of diabetes. SF increases with age and is also greater among individuals with diabetes. A familial correlation of SF suggests that genetics may play a role. We therefore performed parallel genome-wide association studies of SF in two cohorts.

METHODS

Cohort 1 included 1,082 participants, 35-67 years of age with type 1 diabetes. Cohort 2 included 8,721 participants without diabetes, aged 18-90 years.

RESULTS

rs1495741 was significantly associated with SF in Cohort 1 (p < 6 × 10(-10)), which is known to tag the NAT2 acetylator phenotype. The fast acetylator genotype was associated with lower SF, explaining up to 15% of the variance. In Cohort 2, the top signal associated with SF (p = 8.3 × 10(-42)) was rs4921914, also in NAT2, 440 bases upstream of rs1495741 (linkage disequilibrium r (2) = 1.0 for rs4921914 with rs1495741). We replicated these results in two additional cohorts, one with and one without type 1 diabetes. Finally, to understand which compounds are contributing to the NAT2-SF signal, we examined 11 compounds assayed from skin biopsies (n = 198): the fast acetylator genotype was associated with lower levels of the AGEs hydroimidazolones of glyoxal (p = 0.017).

CONCLUSIONS/INTERPRETATION

We identified a robust association between NAT2 and SF in people with and without diabetes. Our findings provide proof of principle that genetic variation contributes to interindividual SF and that NAT2 acetylation status plays a major role.

Combined effects of smoking and alcohol on metabolic syndrome: the LifeLines cohort study

Introduction

The development of metabolic syndrome (MetS) is influenced by environmental factors such as smoking and alcohol consumption. We determined the combined effects of smoking and alcohol on MetS and its individual components.

Methods

64,046 participants aged 18–80 years from the LifeLines Cohort study were categorized into three body mass index (BMI) classes (BMI<25, normal weight; BMI 25–30, overweight; BMI≥30 kg/m², obese). MetS was defined according to the revised criteria of the National Cholesterol Education Program’s Adult Treatment Panel III (NCEP ATP III). Within each BMI class and smoking subgroup (non-smoker, former smoker, <20 and ≥20 g tobacco/day), the cross-sectional association between alcohol and individual MetS components was tested using regression analysis.

Results

Prevalence of MetS varied greatly between the different smoking-alcohol subgroups (1.7–71.1%). HDL cholesterol levels in all alcohol drinkers were higher than in non-drinkers (0.02 to 0.29 mmol/L, P values<0.001). HDL cholesterol levels were lower when they were also a former or current smoker (<20 and ≥20 g tobacco/day). Consumption of ≤1 drink/day indicated a trend towards lower triglyceride levels (non-significant). Concurrent use alcohol (>1 drink/day) and tobacco showed higher triglycerides levels. Up to 2 drinks/day was associated with a smaller waist circumference in overweight and obese individuals. Consumption of >2 drinks/day increased blood pressure, with the strongest associations found for heavy smokers. The overall metabolic profile of wine drinkers was better than that of non-drinkers or drinkers of beer or spirits/mixed drinks.

Conclusion

Light alcohol consumption may moderate the negative associations of smoking with MetS. Our results suggest that the lifestyle advice that emphasizes smoking cessation and the restriction of alcohol consumption to a maximum of 1 drink/day, is a good approach to reduce the prevalence of MetS.

The prevalence of metabolic syndrome and metabolically healthy obesity in Europe: a collaborative analysis of ten large cohort studies

BACKGROUND

Not all obese subjects have an adverse metabolic profile predisposing them to developing type 2 diabetes or cardiovascular disease. The BioSHaRE-EU Healthy Obese Project aims to gain insights into the consequences of (healthy) obesity using data on risk factors and phenotypes across several large-scale cohort studies. Aim of this study was to describe the prevalence of obesity, metabolic syndrome (MetS) and metabolically healthy obesity (MHO) in ten participating studies.

METHODS

Ten different cohorts in seven countries were combined, using data transformed into a harmonized format. All participants were of European origin, with age 18-80 years. They had participated in a clinical examination for anthropometric and blood pressure measurements. Blood samples had been drawn for analysis of lipids and glucose. Presence of MetS was assessed in those with obesity (BMI ≥ 30 kg/m2) based on the 2001 NCEP ATP III criteria, as well as an adapted set of less strict criteria. MHO was defined as obesity, having none of the MetS components, and no previous diagnosis of cardiovascular disease.

RESULTS

Data for 163,517 individuals were available; 17% were obese (11,465 men and 16,612 women). The prevalence of obesity varied from 11.6% in the Italian CHRIS cohort to 26.3% in the German KORA cohort. The age-standardized percentage of obese subjects with MetS ranged in women from 24% in CHRIS to 65% in the Finnish Health2000 cohort, and in men from 43% in CHRIS to 78% in the Finnish DILGOM cohort, with elevated blood pressure the most frequently occurring factor contributing to the prevalence of the metabolic syndrome. The age-standardized prevalence of MHO varied in women from 7% in Health2000 to 28% in NCDS, and in men from 2% in DILGOM to 19% in CHRIS. MHO was more prevalent in women than in men, and decreased with age in both sexes.

CONCLUSIONS

Through a rigorous harmonization process, the BioSHaRE-EU consortium was able to compare key characteristics defining the metabolically healthy obese phenotype across ten cohort studies. There is considerable variability in the prevalence of healthy obesity across the different European populations studied, even when unified criteria were used to classify this phenotype.

Associations between smoking, components of metabolic syndrome and lipoprotein particle size

BACKGROUND

The clustering of metabolic and cardiovascular risk factors is known as metabolic syndrome (MetS). The risk of having MetS is strongly associated with increased adiposity and can be further modified by smoking behavior. Apolipoproteins (apo) associated with low-density lipoprotein-cholesterol (LDL-C) and high-density lipoprotein-cholesterol (HDL-C) may be altered in MetS. This study aimed to examine the association between smoking and the following parameters: MetS and its components, levels of apolipoproteins and estimated lipoprotein particle size, separately for men and women, and in different body mass index (BMI) classes.

METHODS

We included 24,389 men and 35,078 women aged between 18 and 80 years who participated in the LifeLines Cohort Study between December 2006 and January 2012; 5,685 men and 6,989 women were current smokers. Participants were categorized into three different body mass index (BMI) classes (BMI <25; BMI 25 to 30; BMI ≥30 kg/m²). MetS was defined according to the National Cholesterol Education Program's Adult Treatment Panel III (NCEP:ATPIII) criteria. Blood pressure, anthropometric and lipid measurements were rigorously standardized, and the large sample size enabled a powerful estimate of quantitative changes. The association between smoking and the individual MetS components, and apoA1 and apoB, was tested with linear regression. Logistic regression was used to examine the effect of smoking and daily tobacco smoked on risk of having MetS. All models were age adjusted and stratified by sex and BMI class.

RESULTS

Prevalence of MetS increased with higher BMI levels. A total of 64% of obese men and 42% of obese women had MetS. Current smoking was associated with a higher risk of MetS in both sexes and all BMI classes (odds ratio 1.7 to 2.4 for men, 1.8 to 2.3 for women, all P values <0.001). Current smokers had lower levels of HDL cholesterol and apoA1, higher levels of triglycerides and apoB, and higher waist circumference than non-smokers (all P <0.001). Smoking had no consistent association with blood pressure or fasting blood glucose. In all BMI classes, we found a dose-dependent association of daily tobacco consumption with MetS prevalence as well as with lower levels of HDL cholesterol, higher triglyceride levels and lower ratios of HDL cholesterol/apoA1 and, only in those with BMI <30, LDL cholesterol/apoB (all P <0.001).

CONCLUSIONS

Smoking is associated with an increased prevalence of MetS, independent of sex and BMI class. This increased risk is mainly related to lower HDL cholesterol, and higher triglycerides and waist circumference. In addition, smoking was associated with unfavorable changes in apoA1 and apoB, and in lipoprotein particle size. Please see related commentary: http://www.biomedcentral.com/1741-7015/11/196.

Gene-Lifestyle Interactions in Obesity

Obesity is a complex multifaceted disease resulting from interactions between genetics and lifestyle. The proportion of phenotypic variance ascribed to genetic variance is 0.4 to 0.7 for obesity and recent years have seen considerable success in identifying disease-susceptibility variants. Although with the advent of genome-wide association studies the list of genetic variants predisposing to obesity has significantly increased the identified variants only explain a fraction of disease heritability. Studies of gene-environment interactions can provide more insight into the biological mechanisms involved in obesity despite the challenges associated with such designs. Epigenetic changes that affect gene function without DNA sequence modifications may be a key factor explaining interindividual differences in obesity, with both genetic and environmental factors influencing the epigenome. Disentangling the relative contributions of genetic, environmental and epigenetic marks to the establishment of obesity is a major challenge given the complex interplay between these determinants.

Common variants in the type 2 diabetes KCNQ1 gene are associated with impairments in insulin secretion during hyperglycaemic glucose clamp

BACKGROUND

Genome-wide association studies in Japanese populations recently identified common variants in the KCNQ1 gene to be associated with type 2 diabetes. We examined the association of these variants within KCNQ1 with type 2 diabetes in a Dutch population, investigated their effects on insulin secretion and metabolic traits and on the risk of developing complications in type 2 diabetes patients.

METHODOLOGY

The KCNQ1 variants rs151290, rs2237892, and rs2237895 were genotyped in a total of 4620 type 2 diabetes patients and 5285 healthy controls from the Netherlands. Data on macrovascular complications, nephropathy and retinopathy were available in a subset of diabetic patients. Association between genotype and insulin secretion/action was assessed in the additional sample of 335 individuals who underwent a hyperglycaemic clamp.

PRINCIPAL FINDINGS

We found that all the genotyped KCNQ1 variants were significantly associated with type 2 diabetes in our Dutch population, and the association of rs151290 was the strongest (OR 1.20, 95% CI 1.07-1.35, p = 0.002). The risk C-allele of rs151290 was nominally associated with reduced first-phase glucose-stimulated insulin secretion, while the non-risk T-allele of rs2237892 was significantly correlated with increased second-phase glucose-stimulated insulin secretion (p = 0.025 and 0.0016, respectively). In addition, the risk C-allele of rs2237892 was associated with higher LDL and total cholesterol levels (p = 0.015 and 0.003, respectively). We found no evidence for an association of KCNQ1 with diabetic complications.

CONCLUSIONS

Common variants in the KCNQ1 gene are associated with type 2 diabetes in a Dutch population, which can be explained at least in part by an effect on insulin secretion. Furthermore, our data suggest that KCNQ1 is also associated with lipid metabolism.

A genome-wide association search for type 2 diabetes genes in African Americans

African Americans are disproportionately affected by type 2 diabetes (T2DM) yet few studies have examined T2DM using genome-wide association approaches in this ethnicity. The aim of this study was to identify genes associated with T2DM in the African American population. We performed a Genome Wide Association Study (GWAS) using the Affymetrix 6.0 array in 965 African-American cases with T2DM and end-stage renal disease (T2DM-ESRD) and 1029 population-based controls. The most significant SNPs (n = 550 independent loci) were genotyped in a replication cohort and 122 SNPs (n = 98 independent loci) were further tested through genotyping three additional validation cohorts followed by meta-analysis in all five cohorts totaling 3,132 cases and 3,317 controls. Twelve SNPs had evidence of association in the GWAS (P<0.0071), were directionally consistent in the Replication cohort and were associated with T2DM in subjects without nephropathy (P<0.05). Meta-analysis in all cases and controls revealed a single SNP reaching genome-wide significance (P<2.5×10(-8)). SNP rs7560163 (P = 7.0×10(-9), OR (95% CI) = 0.75 (0.67-0.84)) is located intergenically between RND3 and RBM43. Four additional loci (rs7542900, rs4659485, rs2722769 and rs7107217) were associated with T2DM (P<0.05) and reached more nominal levels of significance (P<2.5×10(-5)) in the overall analysis and may represent novel loci that contribute to T2DM. We have identified novel T2DM-susceptibility variants in the African-American population. Notably, T2DM risk was associated with the major allele and implies an interesting genetic architecture in this population. These results suggest that multiple loci underlie T2DM susceptibility in the African-American population and that these loci are distinct from those identified in other ethnic populations.

Physical activity attenuates the influence of FTO variants on obesity risk: a meta-analysis of 218,166 adults and 19,268 children

BACKGROUND

The FTO gene harbors the strongest known susceptibility locus for obesity. While many individual studies have suggested that physical activity (PA) may attenuate the effect of FTO on obesity risk, other studies have not been able to confirm this interaction. To confirm or refute unambiguously whether PA attenuates the association of FTO with obesity risk, we meta-analyzed data from 45 studies of adults (n = 218,166) and nine studies of children and adolescents (n = 19,268).

METHODS AND FINDINGS

All studies identified to have data on the FTO rs9939609 variant (or any proxy [r(2)>0.8]) and PA were invited to participate, regardless of ethnicity or age of the participants. PA was standardized by categorizing it into a dichotomous variable (physically inactive versus active) in each study. Overall, 25% of adults and 13% of children were categorized as inactive. Interaction analyses were performed within each study by including the FTO×PA interaction term in an additive model, adjusting for age and sex. Subsequently, random effects meta-analysis was used to pool the interaction terms. In adults, the minor (A-) allele of rs9939609 increased the odds of obesity by 1.23-fold/allele (95% CI 1.20-1.26), but PA attenuated this effect (p(interaction)  = 0.001). More specifically, the minor allele of rs9939609 increased the odds of obesity less in the physically active group (odds ratio  = 1.22/allele, 95% CI 1.19-1.25) than in the inactive group (odds ratio  = 1.30/allele, 95% CI 1.24-1.36). No such interaction was found in children and adolescents.

CONCLUSIONS

The association of the FTO risk allele with the odds of obesity is attenuated by 27% in physically active adults, highlighting the importance of PA in particular in those genetically predisposed to obesity.

Twelve type 2 diabetes susceptibility loci identified through large-scale association analysis

By combining genome-wide association data from 8,130 individuals with type 2 diabetes (T2D) and 38,987 controls of European descent and following up previously unidentified meta-analysis signals in a further 34,412 cases and 59,925 controls, we identified 12 new T2D association signals with combined P<5x10(-8). These include a second independent signal at the KCNQ1 locus; the first report, to our knowledge, of an X-chromosomal association (near DUSP9); and a further instance of overlap between loci implicated in monogenic and multifactorial forms of diabetes (at HNF1A). The identified loci affect both beta-cell function and insulin action, and, overall, T2D association signals show evidence of enrichment for genes involved in cell cycle regulation. We also show that a high proportion of T2D susceptibility loci harbor independent association signals influencing apparently unrelated complex traits.

Obesity genes identified in genome-wide association studies are associated with adiposity measures and potentially with nutrient-specific food preference

BACKGROUND

New genetic loci, most of which are expressed in the brain, have recently been reported to contribute to the development of obesity. The brain, especially the hypothalamus, is strongly involved in regulating weight and food intake.

OBJECTIVES

We investigated whether the recently reported obesity loci are associated with measures of abdominal adiposity and whether these variants affect dietary energy or macronutrient intake.

DESIGN

We studied 1700 female Dutch participants in the European Prospective Investigation into Cancer and Nutrition (EPIC). Their anthropometric measurements and intake of macronutrients were available. Genotyping was performed by using KASPar chemistry. A linear regression model, with an assumption of an additive effect, was used to analyze the association between genotypes of 12 single nucleotide polymorphisms (SNPs) and adiposity measures and dietary intake.

RESULTS

Seven SNPs were associated (P < 0.05) with weight, body mass index (BMI), and waist circumference (unadjusted for BMI). They were in or near to 6 loci: FTO, MC4R, KCTD15, MTCH2, NEGR1, and BDNF. Five SNPs were associated with dietary intake (P < 0.05) and were in or near 5 loci: SH2B1 (particularly with increased fat), KCTD15 (particularly with carbohydrate intake), MTCH2, NEGR1, and BDNF.

CONCLUSIONS

We confirmed some of the findings for the newly identified obesity loci that are associated with general adiposity in a healthy Dutch female population. Our results suggest that these loci are not specifically associated with abdominal adiposity but more generally with obesity. We also found that some of the SNPs were associated with macronutrient-specific food intake.

Upstream transcription factor 1 (USF1) in risk of type 2 diabetes: association study in 2000 Dutch Caucasians

Type 2 diabetes shares substantial genetic and phenotypic overlap with familial combined hyperlipidemia. Upstream stimulatory factor 1 (USF1), a well-established susceptibility gene for familial combined hyperlipidemia, is postulated to be such a shared genetic determinant. We evaluated two established variants in familial combined hyperlipidemia (rs2073658 and rs3737787) for association with type 2 diabetes in two Dutch case-control samples (N=2011). The first case-control sample comprised 501 subjects with type 2 diabetes from the Breda cohort and 920 healthy blood bank donors of Dutch Caucasian origin. The second case-control sample included 211 subjects with type 2 diabetes, and 379 normoglycemic controls. SNP rs2073658 and SNP rs3737787 were in perfect linkage disequilibrium. In the first case-control sample, prevalence of the major allele was higher in patients than in controls (75% versus 71%, OR=1.25, p=0.018). A similar effect-size and -direction was observed in the second case-control sample (76% versus 72%, OR=1.22, p=0.16). A combined analysis strengthened the evidence for association (OR=1.23, p=0.006). Notably, the increased risk for type 2 diabetes could be ascribed to the major allele, and its high frequency translated to a substantial population attributable risk of 14.5%. In conclusion, the major allele of rs2073658 in the USF1 gene is associated with a modestly increased risk to develop type 2 diabetes in Dutch Caucasians, with considerable impact at the population level.