Differential DNA methylation at birth associated with mental disorder in individuals with 22q11.2 deletion syndrome.

Individuals with 22q11.2 deletion syndrome (DS) have an increased risk of comorbid mental disorders including schizophrenia, attention deficit hyperactivity disorder, depression, as well as intellectual disability. Although most 22q11.2 deletion carriers have the long 3-Mb form of the hemizygous deletion, there remains a large variation in the development and progression of psychiatric disorders, which suggests that alternative factors contribute to the pathogenesis. In this study we investigated whether neonatal DNA methylation signatures in individuals with the 22q11.2 deletion associate with mental disorder later in life. DNA methylation was measured genome-wide from neonatal dried blood spots in a cohort of 164 individuals with 22q11.2DS, including 48 individuals diagnosed with a psychiatric disorder. Among several CpG sites with P-value<10-6, we identified cg23546855 (P-value=2.15 × 10-7) mapping to STK32C to be associated with a later psychiatric diagnosis. Pathway analysis of the top findings resulted in the identification of several Gene Ontology pathways to be significantly enriched (P-value<0.05 after Benjamini-Hochberg correction); among them are the following: neurogenesis, neuron development, neuron projection development, astrocyte development, axonogenesis and axon guidance. In addition, we identified differentially methylated CpG sites in LRP2BP (P-value=5.37 × 10-8) to be associated with intellectual disability (F70-79), in TOP1 (P-value=1.86 × 10-7) with behavioral disorders (F90-98), in NOSIP (P-value=5.12 × 10-8) with disorders of psychological development (F80-89) and in SEMA4B (P-value=4.02 × 10-7) with schizophrenia spectrum disorders (F20-29). In conclusion, our study suggests an association of DNA methylation differences at birth with development of mental disorder later in life in 22q11.2DS individuals.

The effect of GWAS identified BMI loci on changes in body weight among middle-aged Danes during a five-year period.

OBJECTIVE: Genome-wide association studies have identified genetic variants associating with BMI, however, it is un-clarified whether the same variants also influence body weight fluctuations. METHODS: Among 3,982 adult individuals that attended both a baseline and a five-year follow-up examination in the Danish Inter99 intervention study, a genetic risk score (GRS) was constructed based on 30 BMI variants to address whether it is associated with body weight changes. Moreover, it was examined whether the effect of lifestyle changes was modulated by the GRS. RESULTS: The GRS associated strongly with baseline body weight, with a per risk allele increase of 0.45 (0.33-0.58) kg (P = 2.7 × 10(-12) ), corresponding to a body weight difference of 3.41 (2.21-4.60) kg comparing the highest (≥ 30 risk alleles) and lowest (≤ 26 risk alleles) risk allele tertile. No association was observed with changes in body weight during the five years. Changes in lifestyle, including physical activity, diet and smoking habits associated strongly with body weight changes, however, no interactions with the GRS was observed. CONCLUSION: The GRS associated with body weight cross-sectionally, but not with changes over a five-year period. Body weight changes were influenced by lifestyle changes, however, independently of the GRS.

Exome sequencing-driven discovery of coding polymorphisms associated with common metabolic phenotypes.

AIMS/HYPOTHESIS: Human complex metabolic traits are in part regulated by genetic determinants. Here we applied exome sequencing to identify novel associations of coding polymorphisms at minor allele frequencies (MAFs) >1% with common metabolic phenotypes. METHODS: The study comprised three stages. We performed medium-depth (8×) whole exome sequencing in 1,000 cases with type 2 diabetes, BMI >27.5 kg/m(2) and hypertension and in 1,000 controls (stage 1). We selected 16,192 polymorphisms nominally associated (p < 0.05) with case-control status, from four selected annotation categories or from loci reported to associate with metabolic traits. These variants were genotyped in 15,989 Danes to search for association with 12 metabolic phenotypes (stage 2). In stage 3, polymorphisms showing potential associations were genotyped in a further 63,896 Europeans. RESULTS: Exome sequencing identified 70,182 polymorphisms with MAF >1%. In stage 2 we identified 51 potential associations with one or more of eight metabolic phenotypes covered by 45 unique polymorphisms. In meta-analyses of stage 2 and stage 3 results, we demonstrated robust associations for coding polymorphisms in CD300LG (fasting HDL-cholesterol: MAF 3.5%, p = 8.5 × 10(-14)), COBLL1 (type 2 diabetes: MAF 12.5%, OR 0.88, p = 1.2 × 10(-11)) and MACF1 (type 2 diabetes: MAF 23.4%, OR 1.10, p = 8.2 × 10(-10)). CONCLUSIONS/INTERPRETATION: We applied exome sequencing as a basis for finding genetic determinants of metabolic traits and show the existence of low-frequency and common coding polymorphisms with impact on common metabolic traits. Based on our study, coding polymorphisms with MAF above 1% do not seem to have particularly high effect sizes on the measured metabolic traits.

Common variation in oxidative phosphorylation genes is not a major cause of insulin resistance or type 2 diabetes.

AIMS/HYPOTHESIS: There is substantial evidence that mitochondrial dysfunction is linked to insulin resistance and is present in several tissues relevant to the pathogenesis of type 2 diabetes. Here, we examined whether common variation in genes involved in oxidative phosphorylation (OxPhos) contributes to type 2 diabetes susceptibility or influences diabetes-related metabolic traits. METHODS: OxPhos gene variants (n = 10) that had been nominally associated (p < 0.01) with type 2 diabetes in a recent genome-wide meta-analysis (n = 10,108) were selected for follow-up in 3,599 type 2 diabetic and 4,956 glucose-tolerant Danish individuals. A meta-analysis of these variants was performed in 11,729 type 2 diabetic patients and 43,943 non-diabetic individuals. The impact on OGTT-derived metabolic traits was evaluated in 5,869 treatment-naive individuals from the Danish Inter99 study. RESULTS: The minor alleles of COX10 rs9915302 (p = 0.02) and COX5B rs1466100 (p = 0.005) showed nominal association with type 2 diabetes in our Danish cohort. However, in the meta-analysis, none of the investigated variants showed a robust association with type 2 diabetes after correction for multiple testing. Among the alleles potentially associated with type 2 diabetes, none negatively influenced surrogate markers of insulin sensitivity in non-diabetic participants, while the minor alleles of UQCRC1 rs2228561 and COX10 rs10521253 showed a weak (p < 0.01 to p < 0.05) negative influence on indices of glucose-stimulated insulin secretion. CONCLUSIONS/INTERPRETATION: We cannot rule out the possibility that common variants in or near OxPhos genes may influence beta cell function in non-diabetic individuals. However, our quantitative trait studies and a sufficiently large meta-analysis indicate that common variation in proximity to the examined OxPhos genes is not a major cause of insulin resistance or type 2 diabetes.

Type 2 diabetes risk allele near CENTD2 is associated with decreased glucose-stimulated insulin release.

AIMS/HYPOTHESIS: By combining multiple genome-wide association (GWA) studies and comprehensive replication efforts, 12 novel type 2 diabetes associated loci have recently been discovered. Here we evaluate the effect of lead variants of these loci on estimates of insulin release and insulin resistance derived from an oral glucose tolerance test. METHODS: We examined 12 lead variants in or near HMGA2, CENTD2 (also known as ARAP1), KLF14, PRC1, TP53INP1, ZBED3, ZFAND6, CHCHD9, DUSP9, KCNQ1, BCL11A and HNF1A in 5,722 middle-aged people from the population-based Inter99 sample. RESULTS: Carriers of the major diabetogenic allele of rs1552224 in CENTD2 had increased 30-min plasma glucose values (2.0%, p = 2 × 10(-5)) as well as 4.2% reduced insulin release 30 min after an oral glucose load (p = 0.001). Risk allele carriers also had decreased BIGTT-acute insulin release (AIR), which is a surrogate measure of insulin release where sex, BMI, plasma glucose and serum insulin are integrated (5.3%, p = 8 × 10(-7)). In addition, a decreased corrected insulin response (CIR; 9.9%, p = 3 × 10(-8)) was observed. For rs5945326 near DUSP9 on the X-chromosome we stratified according to sex. Male carriers of the risk allele showed nominally decreased BIGTT-AIR (2.6%, p = 0.01). No associations with intermediate metabolic traits were found in women. For the remaining ten lead variants no consistent associations were demonstrated. CONCLUSIONS/INTERPRETATION: Of the lead variants from 12 novel type 2 diabetes associated loci, CENTD2 significantly associated with increased plasma glucose values and decreased glucose-stimulated insulin release, suggesting that the diabetogenic effect of this locus is mediated through an impaired pancreatic beta cell function.

The diabetogenic VPS13C/C2CD4A/C2CD4B rs7172432 variant impairs glucose-stimulated insulin response in 5,722 non-diabetic Danish individuals.

AIMS/HYPOTHESIS: A genome-wide association study in the Japanese population reported two genome-wide significant loci associated with type 2 diabetes of which the VPS13C/C2CD4A/C2CD4B locus was replicated in Europeans. We looked for potential associations between the diabetogenic VPS13C/C2CD4A/C2CD4B rs7172432 variant and diabetes-related intermediary traits. METHODS: We genotyped the rs7172432 variant in the population-based Inter99 cohort (n = 6,784) and analysed quantitative diabetes-related traits in 5,722 non-diabetic participants who all were examined by an OGTT. RESULTS: The diabetes-associated A allele was associated with 0.60 cm higher waist circumference (p = 0.004), 0.037 mmol/l higher fasting plasma glucose (p = 4 × 10(-5)) and 0.11 mmol/l higher plasma glucose at 30 min during an OGTT (p = 4 × 10(-4)). In analyses adjusted for concomitant insulin sensitivity levels the diabetogenic allele was associated with a lower acute glucose-stimulated insulin response (GSIR) as estimated by 30 min serum insulin (ß = -0.039, p = 2 × 10(-7)), insulinogenic index (ß = -0.057, p = 1 × 10(-8)) and BIGTT-acute insulin release (ß = -0.041, p = 9 × 10(-9)). As rs7172432 is situated in a region previously associated with glycaemic traits, we tested linkage disequilibrium (LD) with the reported regional lead single-nucleotide polymorphisms for fasting (rs11071657) and 2 h plasma glucose (rs17271305), and performed conditional analyses of rs7172432. Rs7172432 showed moderate LD with rs11071657 and rs17271305 (R (2) < 0.34) and we found strong association by almost unchanged effect sizes of rs7172432 with plasma glucose and estimates of GSIR in analyses conditional on rs11071657 and rs17271305. CONCLUSIONS/INTERPRETATION: The diabetogenic VPS13C/C2CD4A/C2CD4B rs7172432 A allele associates with GSIR in non-diabetic individuals from the general population, suggesting an impaired beta cell function as an intermediary diabetes-related trait.

Combined analysis of 19 common validated type 2 diabetes susceptibility gene variants shows moderate discriminative value and no evidence of gene-gene interaction.

AIMS/HYPOTHESIS: The list of validated type 2 diabetes susceptibility variants has recently been expanded from three to 19. The variants identified are common and have low penetrance in the general population. The aim of the study is to investigate the combined effect of the 19 variants by applying receiver operating characteristics (ROC) to demonstrate the discriminatory value between glucose-tolerant individuals and type 2 diabetes patients in a cross-sectional population of Danes. METHODS: The 19 variants were genotyped in three study populations: the population-based Inter99 study; the ADDITION study; and additional type 2 diabetic patients and glucose-tolerant individuals. The case-control studies involved 4,093 type 2 diabetic patients and 5,302 glucose-tolerant individuals. RESULTS: Single-variant analyses demonstrated allelic odds ratios ranging from 1.04 (95% CI 0.98-1.11) to 1.33 (95% CI 1.22-1.45). When combining the 19 variants, subgroups with extreme risk profiles showed a threefold difference in the risk of type 2 diabetes (lower 10% carriers with < or =15 risk alleles vs upper 10% carriers with > or =22 risk alleles, OR 2.93 (95% CI 2.38-3.62, p = 1.6 x 10(-25)). We calculated the area under a ROC curve to estimate the discrimination rate between glucose-tolerant individuals and type 2 diabetes patients based on the 19 variants. We found an area under the ROC curve of 0.60. Two-way gene-gene interaction showed few nominal interaction effects. CONCLUSIONS/INTERPRETATION: Combined analysis of the 19 validated variants enables detection of subgroups at substantially increased risk of type 2 diabetes; however, the discrimination between glucose-tolerant and type 2 diabetes individuals is still too inaccurate to achieve clinical value.

Impact of polymorphisms in WFS1 on prediabetic phenotypes in a population-based sample of middle-aged people with normal and abnormal glucose regulation.

AIM/HYPOTHESIS: Recently, variants in WFS1 have been shown to be associated with type 2 diabetes. We aimed to examine metabolic risk phenotypes of WFS1 variants in glucose-tolerant people and in individuals with abnormal glucose regulation. METHODS: The type 2 diabetes-associated WFS1 variant rs734312 (His611Arg) was studied in the population-based Inter99 cohort involving 4,568 glucose-tolerant individuals and 1,471 individuals with treatment-naive abnormal glucose regulation, and in an additional 3,733 treated type 2 diabetes patients. RESULTS: The WFS1 rs734312 showed a borderline significant association with type 2 diabetes with directions and relative risks consistent with previous reports. In individuals with abnormal glucose regulation, the diabetogenic risk A allele of rs734312 was associated in an allele-dependent manner with a decrease in insulinogenic index (p = 0.025) and decreased 30-min serum insulin levels (p = 0.047) after an oral glucose load. In glucose-tolerant individuals the same allele was associated with increased fasting serum insulin concentration (p = 0.019) and homeostasis model assessment of insulin resistance (HOMA-IR; p = 0.026). To study the complex interaction of WFS1 rs734312 on insulin release and insulin resistance we introduced Hotelling's T (2) test. Assuming bivariate normal distribution, we constructed standard error ellipses of the insulinogenic index and HOMA-IR when stratified according to glucose tolerance status around the means of each WFS1 rs734312 genotype level. The interaction term between individuals with normal glucose tolerance and abnormal glucose regulation on the insulinogenic index and HOMA-IR was significantly associated with the traits (p = 0.0017). CONCLUSIONS/INTERPRETATION: Type 2 diabetes-associated risk alleles of WFS1 are associated with estimates of a decreased pancreatic beta cell function among middle-aged individuals with abnormal glucose regulation.

The GCKR rs780094 polymorphism is associated with elevated fasting serum triacylglycerol, reduced fasting and OGTT-related insulinaemia, and reduced risk of type 2 diabetes.

AIMS/HYPOTHESIS: Recent genome-wide association studies have suggested that a polymorphism in GCKR, the gene encoding the glucokinase regulatory protein, is involved in triacylglycerol regulation. Our aim was to examine in large-scale studies the common GCKR rs780094 polymorphism in relation to metabolic traits (mainly fasting hypertriacylglycerolaemia) and traits related to pancreatic beta cell function. METHODS: The polymorphism was genotyped in 16,853 Danes using Taqman allelic discrimination. Association was analysed in case-control studies and quantitative trait analyses. We also analysed the possible interactive effect between the GCK -30G>A polymorphism and the GCKR rs780094 variant on metabolic traits. RESULTS: The minor GCKR A-allele of rs780094 is associated with an increased level of fasting serum triacylglycerol (p = 6 x 10(-14)), impaired fasting (p = 0.001) and OGTT-related insulin release (p = 3 x 10(-6)), reduced homeostasis model assessment of insulin resistance (p = 0.0004), WHO-defined dyslipidaemia (p = 6 x 10(-9)) and a modestly decreased risk of type 2 diabetes (p = 0.01). Significantly increased fasting serum insulin concentrations were demonstrated when analysing the GCK -30A and GCKR rs780094 G-alleles in an additive model. CONCLUSIONS/INTERPRETATION: The GCKR rs780094 polymorphism, or another variant with which it is in tight linkage disequilibrium, is likely to increase glucokinase regulatory protein activity to induce improved glycaemic regulation at the expense of hypertriacylglycerolaemia as reflected in the present study of 16,853 Danes. We also suggest an additive effect of GCK and GCKR risk alleles on [corrected] serum insulin release.

A -243A-->G polymorphism upstream of the gene encoding GAD65 associates with lower levels of body mass index and glycaemia in a population-based sample of 5857 middle-aged White subjects.

AIMS: The glutamate decarboxylase gene (GAD2) encodes GAD65, an enzyme catalysing the production of the gamma-aminobutyric acid (GABA) which interacts with neuropeptide Y to stimulate food intake. It has been suggested that in pancreatic islets, GABA serves as a functional regulator of pancreatic hormone release. Conflicting results have been reported concerning the potential impact of GAD2 variation on estimates of energy metabolism. The aim of this study was to elucidate potential associations between the GAD2-243A-->G polymorphism and levels of body mass index (BMI) and estimates of glycaemia. METHODS: Using high-throughput chip-based matrix-assisted laser desorption/ionization time-of-flight mass spectrometry, the GAD2-243A-->G (rs2236418) polymorphism was genotyped in a population-based sample (Inter99) of 5857 middle-aged, unrelated Danish White subjects. RESULTS: The G-allele was associated with modestly lower BMI (P = 0.01). In a case-control study of obesity, the G-allele frequency in 2582 participants with BMI < 25 kg/m2 was 19.5% (18.4-20.6) compared with 17.1% (15.5-18.8) in 968 participants having BMI > or = 30 kg/m2 (P = 0.03), odds ratio 0.9 (0.7-1.0). Of the 5857 subjects, GG carriers had lower fasting plasma glucose levels (mmol/l) [AA (n = 3859) 5.6 +/- 0.8; AG (n = 1792) 5.5 +/- 0.8; GG (n = 206) 5.5 +/- 0.8, P = 0.008] and lower 30-min oral glucose tolerance test (OGTT)-related plasma glucose levels (AA 8.7 +/- 1.9; AG 8.6 +/- 1.9; GG 8.6 +/- 2.0, P = 0.04), adjusted for sex, age and BMI. Analysing subjects who were both normoglycaemic and glucose tolerant (n = 4431) GG carriers still had lower fasting plasma glucose concentrations: AA (n = 2895) 5.3 +/- 0.4; AG (n = 1383) 5.3 +/- 0.4; GG (n = 153) 5.2 +/- 0.4 (P = 9.10(-5)). CONCLUSION: The present study suggests that the GAD2-243A-->G polymorphism in a population of middle-aged White people associates with a modest reduction in BMI and fasting and OGTT-related plasma glucose levels.

Genetic analysis of the estrogen-related receptor alpha and studies of association with obesity and type 2 diabetes.

BACKGROUND: The estrogen-related receptor alpha (ERRalpha or NR3B1) is a transcription factor from the nuclear receptor super-family, group III. The gene encoding ERRalpha (ESRRA) is located on chromosome 11q13, a region showing genetic linkage to body mass index and fat percentage. Through interaction with the peroxisome proliferator-activated receptor-gamma coactivator-1alpha (PGC-1alpha), ERRalpha regulates key enzymes involved in the beta-oxidation of fatty acids. RESULTS: By screening 48 overweight or obese subjects for variants in the exons, exon-intron boundaries and 1000 base pairs (bp) of the promoter region of ESRRA using bi-directional nucleotide sequencing, we identified seven variants. Four rare variants had minor allele frequencies (MAF) below 1%: Pro369Pro, Gly406Asp, 3'UTR+418G>A, 3'UTR+505C>A. Two single-nucleotide polymorphisms, Pro116Pro and IVS6+65C>T (MAF 15%), were in complete linkage disequilibrium (LD) (r (2)=1). We also confirmed the presence of a reported 23 bp microsatellite repeat (ESRRA23). The Pro116Pro and ESRRA23 variants were not associated with obesity, type 2 diabetes or related phenotypes in a large population-based study of 6365 Danish whites. The two variants were examined for interactions with variants in the peroxisome proliferator-activated receptor-gamma coactivator-1alpha and -beta; however, no evidence of epistatic effects between the variants was demonstrated. CONCLUSION: The ESRRA23 and Pro116Pro variants of the gene encoding ERRalpha are not associated with obesity, type 2 diabetes or related quantitative traits in the examined Danish whites.

Large-scale study of the -232C > G polymorphism of PCK1 in Type 2 diabetes.

AIMS: Phosphoenolpyruvate carboxykinase (PEPCK) is a catalyst of the rate-limiting step in the gluconeogenic pathway and is regulated at the transcriptional level predominantly by insulin, glucocorticoids, glucagon, and cAMP. The -232C > G polymorphism in the gene encoding PEPCK (PCK1) is reported to associate with Type 2 diabetes in Canadian Caucasians and Oji-Cree populations. We have estimated the impact of the PCK1-232C > G polymorphism in a relatively large-scale case-control study of Type 2 diabetes and in association studies of common metabolic phenotypes. Interaction studies of the PCK1-232C > G polymorphism with variants in the genes encoding peroxisome proliferator-activated receptor-gamma co-activator (PGC)-1alpha and hepatic nuclear factor (HNF)-4alpha were also performed. METHODS: PCK1-232C > G was genotyped in a total of 7467 Danish white subjects using TaqMan allelic discrimination. A case-control study of Type 2 diabetes was performed using 6057 of the participants, and quantitative trait studies of metabolic variables were carried out in a subgroup of 5718 non-diabetic subjects. Additionally, variants in PGC-1alpha (Gly482Ser) and HNF-4alpha (Thr130Ile, Val255Met, and rs1884614) were investigated for epistatic interaction with the PCK1-232C > G polymorphism. RESULTS: In the case-control study of Type 2 diabetes of 1377 Type 2 diabetic patients and 4680 normoglycaemic and normal glucose-tolerant subjects we found no association of the PCK1-232C > G polymorphism with diabetes. In addition, the variant was not associated with age of clinical onset of Type 2 diabetes. In the study of 5718 non-diabetic subjects, we found no relationships of quantitative metabolic traits with the PCK1-232C > G polymorphism. We failed to demonstrate any convincing epistatic effects of the variants in the genes encoding PGC-1alpha and HNF-4alpha with the PCK1-232C > G polymorphism. CONCLUSIONS: The PCK1-232C > G polymorphism is not a major contributor to the pathogenesis of Type 2 diabetes in the Danish population.