👤 Tim M Frayling

The extent to which genetic variations contribute to interindividual differences in weight loss and metabolic outcomes after bariatric surgery is unknown. Identifying genetic variants that impact surgery outcomes may contribute to clinical decision making. This review evaluates current evidence addressing the association of genetic variants with weight loss and changes in metabolic parameters after bariatric surgery. A search was conducted using Medline, Embase, Scopus, Web of Science, and Cochrane Library. Fifty-two eligible studies were identified. Single nucleotide polymorphisms (SNPs) at ADIPOQ (rs226729, rs1501299, rs3774261, and rs17300539) showed a positive association with postoperative change in measures of glucose homeostasis and lipid profiles (n = 4), but not with weight loss after surgery (n = 6). SNPs at FTO (rs11075986, rs16952482, rs8050136, rs9939609, rs9930506, and rs16945088) (n = 10) and MC4R (rs11152213, rs476828, rs2229616, rs9947255, rs17773430, rs5282087, and rs17782313) (n = 9) were inconsistently associated with weight loss and metabolic improvement. Four studies examining the UCP2 SNP rs660339 reported associations with postsurgical weight loss. In summary, there is limited evidence supporting a role for specific genetic variants in surgical outcomes after bariatric surgery. Most studies have adopted a candidate gene approach, limiting the scope for discovery, suggesting that the absence of compelling evidence is not evidence of absence. Show less

Pancreatic beta-cell glucose sensitivity is the slope of the plasma glucose-insulin secretion relationship and is a key predictor of deteriorating glucose tolerance and development of type 2 diabetes. However, there are no large-scale studies looking at the genetic determinants of beta-cell glucose sensitivity. To understand the genetic determinants of pancreatic beta-cell glucose sensitivity using genome-wide meta-analysis and candidate gene studies. We performed a genome-wide meta-analysis for beta-cell glucose sensitivity in subjects with type 2 diabetes and nondiabetic subjects from 6 independent cohorts (n = 5706). Beta-cell glucose sensitivity was calculated from mixed meal and oral glucose tolerance tests, and its associations between known glycemia-related single nucleotide polymorphisms (SNPs) and genome-wide association study (GWAS) SNPs were estimated using linear regression models. Beta-cell glucose sensitivity was moderately heritable (h2 ranged from 34% to 55%) using SNP and family-based analyses. GWAS meta-analysis identified multiple correlated SNPs in the CDKAL1 gene and GIPR-QPCTL gene loci that reached genome-wide significance, with SNP rs2238691 in GIPR-QPCTL (P value = 2.64 × 10-9) and rs9368219 in the CDKAL1 (P value = 3.15 × 10-9) showing the strongest association with beta-cell glucose sensitivity. These loci surpassed genome-wide significance when the GWAS meta-analysis was repeated after exclusion of the diabetic subjects. After correction for multiple testing, glycemia-associated SNPs in or near the HHEX and IGF2B2 loci were also associated with beta-cell glucose sensitivity. We show that, variation at the GIPR-QPCTL and CDKAL1 loci are key determinants of pancreatic beta-cell glucose sensitivity. Show less

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity. Show less

We screened variants on an exome-focused genotyping array in >300,000 participants (replication in >280,000 participants) and identified 444 independent variants in 250 loci significantly associated with total cholesterol (TC), high-density-lipoprotein cholesterol (HDL-C), low-density-lipoprotein cholesterol (LDL-C), and/or triglycerides (TG). At two loci (JAK2 and A1CF), experimental analysis in mice showed lipid changes consistent with the human data. We also found that: (i) beta-thalassemia trait carriers displayed lower TC and were protected from coronary artery disease (CAD); (ii) excluding the CETP locus, there was not a predictable relationship between plasma HDL-C and risk for age-related macular degeneration; (iii) only some mechanisms of lowering LDL-C appeared to increase risk for type 2 diabetes (T2D); and (iv) TG-lowering alleles involved in hepatic production of TG-rich lipoproteins (TM6SF2 and PNPLA3) tracked with higher liver fat, higher risk for T2D, and lower risk for CAD, whereas TG-lowering alleles involved in peripheral lipolysis (LPL and ANGPTL4) had no effect on liver fat but decreased risks for both T2D and CAD. Show less

Understanding the physiological mechanisms by which common variants predispose to type 2 diabetes requires large studies with detailed measures of insulin secretion and sensitivity. Here we performed the largest genome-wide association study of first-phase insulin secretion, as measured by intravenous glucose tolerance tests, using up to 5,567 individuals without diabetes from 10 studies. We aimed to refine the mechanisms of 178 known associations between common variants and glycemic traits and identify new loci. Thirty type 2 diabetes or fasting glucose-raising alleles were associated with a measure of first-phase insulin secretion at Show less

Circulating metabolites associated with insulin sensitivity may represent useful biomarkers, but their causal role in insulin sensitivity and diabetes is less certain. We previously identified novel metabolites correlated with insulin sensitivity measured by the hyperinsulinemic-euglycemic clamp. The top-ranking metabolites were in the glutathione and glycine biosynthesis pathways. We aimed to identify common genetic variants associated with metabolites in these pathways and test their role in insulin sensitivity and type 2 diabetes. With 1,004 nondiabetic individuals from the RISC study, we performed a genome-wide association study (GWAS) of 14 insulin sensitivity-related metabolites and one metabolite ratio. We replicated our results in the Botnia study (n = 342). We assessed the association of these variants with diabetes-related traits in GWAS meta-analyses (GENESIS [including RISC, EUGENE2, and Stanford], MAGIC, and DIAGRAM). We identified four associations with three metabolites-glycine (rs715 at CPS1), serine (rs478093 at PHGDH), and betaine (rs499368 at SLC6A12; rs17823642 at BHMT)-and one association signal with glycine-to-serine ratio (rs1107366 at ALDH1L1). There was no robust evidence for association between these variants and insulin resistance or diabetes. Genetic variants associated with genes in the glycine biosynthesis pathways do not provide consistent evidence for a role of glycine in diabetes-related traits. Show less

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. 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. 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)). 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. Show less

Sex hormone-binding globulin (SHBG) is a glycoprotein responsible for the transport and biologic availability of sex steroid hormones, primarily testosterone and estradiol. SHBG has been associated with chronic diseases including type 2 diabetes (T2D) and with hormone-sensitive cancers such as breast and prostate cancer. We performed a genome-wide association study (GWAS) meta-analysis of 21,791 individuals from 10 epidemiologic studies and validated these findings in 7,046 individuals in an additional six studies. We identified twelve genomic regions (SNPs) associated with circulating SHBG concentrations. Loci near the identified SNPs included SHBG (rs12150660, 17p13.1, p = 1.8 × 10(-106)), PRMT6 (rs17496332, 1p13.3, p = 1.4 × 10(-11)), GCKR (rs780093, 2p23.3, p = 2.2 × 10(-16)), ZBTB10 (rs440837, 8q21.13, p = 3.4 × 10(-09)), JMJD1C (rs7910927, 10q21.3, p = 6.1 × 10(-35)), SLCO1B1 (rs4149056, 12p12.1, p = 1.9 × 10(-08)), NR2F2 (rs8023580, 15q26.2, p = 8.3 × 10(-12)), ZNF652 (rs2411984, 17q21.32, p = 3.5 × 10(-14)), TDGF3 (rs1573036, Xq22.3, p = 4.1 × 10(-14)), LHCGR (rs10454142, 2p16.3, p = 1.3 × 10(-07)), BAIAP2L1 (rs3779195, 7q21.3, p = 2.7 × 10(-08)), and UGT2B15 (rs293428, 4q13.2, p = 5.5 × 10(-06)). These genes encompass multiple biologic pathways, including hepatic function, lipid metabolism, carbohydrate metabolism and T2D, androgen and estrogen receptor function, epigenetic effects, and the biology of sex steroid hormone-responsive cancers including breast and prostate cancer. We found evidence of sex-differentiated genetic influences on SHBG. In a sex-specific GWAS, the loci 4q13.2-UGT2B15 was significant in men only (men p = 2.5 × 10(-08), women p = 0.66, heterogeneity p = 0.003). Additionally, three loci showed strong sex-differentiated effects: 17p13.1-SHBG and Xq22.3-TDGF3 were stronger in men, whereas 8q21.12-ZBTB10 was stronger in women. Conditional analyses identified additional signals at the SHBG gene that together almost double the proportion of variance explained at the locus. Using an independent study of 1,129 individuals, all SNPs identified in the overall or sex-differentiated or conditional analyses explained ~15.6% and ~8.4% of the genetic variation of SHBG concentrations in men and women, respectively. The evidence for sex-differentiated effects and allelic heterogeneity highlight the importance of considering these features when estimating complex trait variance. Show less

Proinsulin is a precursor of mature insulin and C-peptide. Higher circulating proinsulin levels are associated with impaired β-cell function, raised glucose levels, insulin resistance, and type 2 diabetes (T2D). Studies of the insulin processing pathway could provide new insights about T2D pathophysiology. We have conducted a meta-analysis of genome-wide association tests of ∼2.5 million genotyped or imputed single nucleotide polymorphisms (SNPs) and fasting proinsulin levels in 10,701 nondiabetic adults of European ancestry, with follow-up of 23 loci in up to 16,378 individuals, using additive genetic models adjusted for age, sex, fasting insulin, and study-specific covariates. Nine SNPs at eight loci were associated with proinsulin levels (P < 5 × 10(-8)). Two loci (LARP6 and SGSM2) have not been previously related to metabolic traits, one (MADD) has been associated with fasting glucose, one (PCSK1) has been implicated in obesity, and four (TCF7L2, SLC30A8, VPS13C/C2CD4A/B, and ARAP1, formerly CENTD2) increase T2D risk. The proinsulin-raising allele of ARAP1 was associated with a lower fasting glucose (P = 1.7 × 10(-4)), improved β-cell function (P = 1.1 × 10(-5)), and lower risk of T2D (odds ratio 0.88; P = 7.8 × 10(-6)). Notably, PCSK1 encodes the protein prohormone convertase 1/3, the first enzyme in the insulin processing pathway. A genotype score composed of the nine proinsulin-raising alleles was not associated with coronary disease in two large case-control datasets. We have identified nine genetic variants associated with fasting proinsulin. Our findings illuminate the biology underlying glucose homeostasis and T2D development in humans and argue against a direct role of proinsulin in coronary artery disease pathogenesis. Show less

Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and ∼ 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 × 10⁻⁸), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation. Show less

Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studies (n = 15,234 nondiabetic individuals) and a follow-up of 29 independent loci (n = 6,958-30,620). We identify variants at the GIPR locus associated with 2-h glucose level (rs10423928, beta (s.e.m.) = 0.09 (0.01) mmol/l per A allele, P = 2.0 x 10(-15)). The GIPR A-allele carriers also showed decreased insulin secretion (n = 22,492; insulinogenic index, P = 1.0 x 10(-17); ratio of insulin to glucose area under the curve, P = 1.3 x 10(-16)) and diminished incretin effect (n = 804; P = 4.3 x 10(-4)). We also identified variants at ADCY5 (rs2877716, P = 4.2 x 10(-16)), VPS13C (rs17271305, P = 4.1 x 10(-8)), GCKR (rs1260326, P = 7.1 x 10(-11)) and TCF7L2 (rs7903146, P = 4.2 x 10(-10)) associated with 2-h glucose. Of the three newly implicated loci (GIPR, ADCY5 and VPS13C), only ADCY5 was found to be associated with type 2 diabetes in collaborating studies (n = 35,869 cases, 89,798 controls, OR = 1.12, 95% CI 1.09-1.15, P = 4.8 x 10(-18)). Show less

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes. Show less

Low plasma levels of carotenoids and tocopherols are associated with increased risk of chronic disease and disability. Because dietary intake of these lipid-soluble antioxidant vitamins is only poorly correlated with plasma levels, we hypothesized that circulating carotenoids (vitamin A-related compounds) and tocopherols (vitamin E-related compounds) are affected by common genetic variation. By conducting a genome-wide association study in a sample of Italians (n = 1190), we identified novel common variants associated with circulating carotenoid levels and known lipid variants associated with alpha-tocopherol levels. Effects were replicated in the Women's Health and Aging Study (n = 615) and in the alpha-Tocopherol, beta-Carotene Cancer Prevention (ATBC) study (n = 2136). In meta-analyses including all three studies, the G allele at rs6564851, near the beta-carotene 15,15'-monooxygenase 1 (BCMO1) gene, was associated with higher beta-carotene (p = 1.6 x 10(-24)) and alpha-carotene (p = 0.0001) levels and lower lycopene (0.003), zeaxanthin (p = 1.3 x 10(-5)), and lutein (p = 7.3 x 10(-15)) levels, with effect sizes ranging from 0.10-0.28 SDs per allele. Interestingly, this genetic variant had no significant effect on plasma retinol (p > 0.05). The SNP rs12272004, in linkage disequilibrium with the S19W variant in the APOA5 gene, was associated with alpha-tocopherol (meta-analysis p = 7.8 x 10(-10)) levels, and this association was substantially weaker when we adjusted for triglyceride levels (p = 0.002). Our findings might shed light on the controversial relationship between lipid-soluble anti-oxidant nutrients and human health. Show less

Common genetic variants influence plasma triacylglycerol, HDL-cholesterol (HDL-C) and glucose levels in cross-sectional studies. However, the longitudinal effects of these established variants have not been studied. Our aim was to examine the longitudinal associations of four such variants in the apolipoprotein A-V (APOA5), lipoprotein lipase (LPL), and glucokinase (GCK) genes with fasting glucose or lipid levels. The individuals analysed were participants in the Busselton Health Survey (n = 4,554). Cross-sectional analyses of family data used the total association test. Longitudinal association analyses of unrelated participant data (n = 2,864) used linear mixed-effects models. The findings of cross-sectional association analyses replicated those of previous studies. We observed associations of the G and C alleles at the APOA5 single nucleotide polymorphisms (SNPs) rs662799 and rs3135506 with raised triacylglycerol levels (p = 0.0003 and p < 0.0001, respectively), the 447X allele at the LPL SNP rs328 with reduced triacylglycerol levels (p = 0.0004) and raised HDL-C levels (p = 0.0004), and the A allele of the GCK SNP rs1799884 with raised fasting glucose level (p = 0.015). Longitudinal association analyses showed that most of these associations did not change in the same individuals over an average follow-up time of 17.4 years, though there was some evidence that the association of the 447X allele of rs328 with raised HDL-C level significantly increased with age (p = 0.01), and that the association of the C allele of rs3135506 with raised triacylglycerol level significantly increased over time (p = 0.0007). The current study suggests that the effects of established gene variants on lipid and glucose traits do not tend to alter with age during adulthood or over time. Show less

Plasma liver-enzyme tests are widely used in the clinic for the diagnosis of liver diseases and for monitoring the response to drug treatment. There is considerable evidence that human genetic variation influences plasma levels of liver enzymes. However, such genetic variation has not been systematically assessed. In the present study, we performed a genome-wide association study of plasma liver-enzyme levels in three populations (total n = 7715) with replication in three additional cohorts (total n = 4704). We identified two loci influencing plasma levels of alanine-aminotransferase (ALT) (CPN1-ERLIN1-CHUK on chromosome 10 and PNPLA3-SAMM50 on chromosome 22), one locus influencing gamma-glutamyl transferase (GGT) levels (HNF1A on chromosome 12), and three loci for alkaline phosphatase (ALP) levels (ALPL on chromosome 1, GPLD1 on chromosome 6, and JMJD1C-REEP3 on chromosome 10). In addition, we confirmed the associations between the GGT1 locus and GGT levels and between the ABO locus and ALP levels. None of the ALP-associated SNPs were associated with other liver tests, suggesting intestine and/or bone specificity. The mechanisms underlying the associations may involve cis- or trans-transcriptional effects (some of the identified variants were associated with mRNA transcription in human liver or lymphoblastoid cells), dysfunction of the encoded proteins (caused by missense variations at the functional domains), or other unknown pathways. These findings may help in the interpretation of liver-enzyme tests and provide candidate genes for liver diseases of viral, metabolic, autoimmune, or toxic origin. The specific associations with ALP levels may point to genes for bone or intestinal diseases. Show less

The APOA5 gene variants, -1131T>C and S19W, are associated with altered triglyceride concentrations in studies of subjects of Caucasian and East Asian descent. There are few studies of these variants in South Asians. We investigated whether the two APOA5 variants also show similar association with various lipid parameters in Indian population as in the UK white subjects. We genotyped 557 Indian adults from Pune, India, and 237 UK white adults for -1131T>C and S19W variants in the APOA5 gene, compared their allelic and genotype frequency and determined their association with fasting serum triglycerides, total cholesterol, HDL and LDL cholesterol levels using univariate general linear analysis. APOC3 SstI polymorphism was also analyzed in 175 Pune Indian subjects for analysis of linkage disequilibrium with the APOA5 variants. The APOA5 -1131C allele was more prevalent in Indians from Pune (Pune Indians) compared to UK white subjects (allele frequency 20% vs. 4%, p = 0.00001), whereas the 19W allele was less prevalent (3% vs. 6% p = 0.0015). Patterns of linkage disequilibrium between the two variants were similar between the two populations and confirmed that they occur on two different haplotypes. In Pune Indians, the presence of -1131C allele and the 19W allele was associated with a 19% and 15% increase respectively in triglyceride concentrations although only -1131C was significant (p = 0.0003). This effect size was similar to that seen in the UK white subjects. Analysis of the APOC3 SstI polymorphism in 175 Pune Indian subjects showed that this variant is not in appreciable linkage disequilibrium with the APOA5 -1131T>C variant (r2 = 0.07). This is the first study to look at the role of APOA5 in Asian Indian subjects that reside in India. The -1131C allele is more prevalent and the 19W allele is less prevalent in Pune Indians compared to UK Caucasians. We confirm that the APOA5 variants are associated with triglyceride levels independent of ethnicity and that this association is similar in magnitude in Asian Indians and Caucasians. The -1131C allele is present in 36% of the Pune Indian population making it a powerful marker for looking at the role of elevated triglycerides in important conditions such as pancreatitis, diabetes and coronary heart disease. Show less

Apolipoprotein AV (ApoAV) gene variant, -1131T>C, is associated with increased triglyceride concentrations in all ethnic groups studied. An MseI based RFLP analysis is the most commonly used method for genotyping this SNP. We genotyped a large cohort comprising 1185 Asian Indians and 173 UK Caucasians for -1131T>C using an ARMS-PCR based tetra-primer method. For quality control, we re-genotyped approximately 10% random samples from this cohort utilizing the MseI RFLP, which showed a 2.9% (3/102) genotyping error rate between the two methods. To investigate further, we sequenced the 900 bp region around the -1131T>C polymorphism in 25 Asian Indians and 15 UK Caucasians and found a number of polymorphisms including the -987C>T polymorphism. Further analysis of the -987C>T SNP showed a higher rare allele frequency of 0.23 in Asian Indians (n = 158) compared to 0.09 in the UK Caucasians (n = 157). This SNP is located 4 bp from the 3' end of the RFLP forward primer and is in weak linkage disequilibrium with -1131T>C variant (r2 = 0.084 and D' = 1). Repeated RFLP analysis of seven subjects heterozygous for -987C>T (seven times), showed discordant results with the sequence at -1131T>C SNP nearly one third (15/49) of the time. We conclude that presence of -987C>T polymorphism in the forward primer of the MseI RFLP assay may lead to allelic drop-out and generate unforeseen errors in genotyping the -1131T>C polymorphism. Our results also emphasise the need for careful quality control in all molecular genetic studies, particularly while transferring genotyping methods between various ethnic groups. Show less