👤 Alan R Shuldiner

Rare cases of monogenic obesity, which may respond to specific therapeutics, can remain undetected in populations in which polygenic obesity is prevalent. This study examined rare DNA variation in established monogenic obesity genes within a community using whole-exome sequence data from 6803 longitudinally studied individuals. Exome data across 15 monogenic obesity genes were analyzed for nonsynonymous variants observed in any child with a maximum BMI z score > 2 (N = 279) but not observed in a child with a maximum BMI z score ≤ 0 (n = 1542) or that occurred in adults in the top 5th percentile of BMI (n = 263) but not in adults below the median BMI (n = 2629). Variants were then functionally analyzed using luciferase assays. The comparisons between cases of obesity and controls identified eight missense variants in six genes: DYRK1B, KSR2, MC4R, NTRK2, PCSK1, and SIM1. Among these, MC4R p.A303P and p.R165G were previously shown to impair MC4R function. Functional analyses of the remaining six variants suggest that KSR2 p.I402F and p.T193I and NTRK2 p.S249Y alter protein function. In addition to MC4R, rare missense variants in KSR2 and NTRK2 may potentially explain the severe obesity observed for the carriers. Show less

Identifying the genetic determinants of inter-individual variation in lipid species (lipidome) may provide deeper understanding and additional insight into the mechanistic effect of complex lipidomic pathways in CVD risk and progression beyond simple traditional lipids. Previous studies have been largely population based and thus only powered to discover associations with common genetic variants. Founder populations represent a powerful resource to accelerate discovery of previously unknown biology associated with rare population alleles that have risen to higher frequency due to genetic drift. We performed a genome-wide association scan of 355 lipid species in 650 individuals from the Amish founder population including 127 lipid species not previously tested. To the best of our knowledge, we report for the first time the lipid species associated with two rare-population but Amish-enriched lipid variants: APOB_rs5742904 and APOC3_rs76353203. We also identified novel associations for 3 rare-population Amish-enriched loci with several sphingolipids and with proposed potential functional/causal variant in each locus including GLTPD2_rs536055318, CERS5_rs771033566, and AKNA_rs531892793. We replicated 7 previously known common loci including novel associations with two sterols: androstenediol with UGT locus and estriol with SLC22A8/A24 locus. Our results show the double power of founder populations and detailed lipidome to discover novel trait-associated variants. Show less

Large-scale human exome sequencing can identify rare protein-coding variants with a large impact on complex traits such as body adiposity. We sequenced the exomes of 645,626 individuals from the United Kingdom, the United States, and Mexico and estimated associations of rare coding variants with body mass index (BMI). We identified 16 genes with an exome-wide significant association with BMI, including those encoding five brain-expressed G protein-coupled receptors ( Show less

Angiopoietin-like 4 (ANGPTL4) is an endogenous inhibitor of lipoprotein lipase that modulates lipid levels, coronary atherosclerosis risk, and nutrient partitioning. We hypothesize that loss of ANGPTL4 function might improve glucose homeostasis and decrease risk of type 2 diabetes (T2D). We investigate protein-altering variants in ANGPTL4 among 58,124 participants in the DiscovEHR human genetics study, with follow-up studies in 82,766 T2D cases and 498,761 controls. Carriers of p.E40K, a variant that abolishes ANGPTL4 ability to inhibit lipoprotein lipase, have lower odds of T2D (odds ratio 0.89, 95% confidence interval 0.85-0.92, p = 6.3 × 10 Show less

Lipoprotein (a) [Lp(a)] is an independent risk factor for atherosclerosis-related events that is under strong genetic control (heritability = 0.68-0.98). However, causal mutations and functional validation of biological pathways modulating Lp(a) metabolism are lacking. We performed a genome-wide association scan to identify genetic variants associated with Lp(a)-cholesterol levels in the Old Order Amish. We confirmed a previously known locus on chromosome 6q25-26 and found Lp(a) levels also to be significantly associated with a SNP near the APOA5-APOA4-APOC3-APOA1 gene cluster on chromosome 11q23 linked in the Amish to the APOC3 R19X null mutation. On 6q locus, we detected associations of Lp(a)-cholesterol with 118 common variants (P = 5 × 10(-8) to 3.91 × 10(-19)) spanning a ∼5.3 Mb region that included the LPA gene. To further elucidate variation within LPA, we sequenced LPA and identified two variants most strongly associated with Lp(a)-cholesterol, rs3798220 (P = 1.07 × 10(-14)) and rs10455872 (P = 1.85 × 10(-12)). We also measured copy numbers of kringle IV-2 (KIV-2) in LPA using qPCR. KIV-2 numbers were significantly associated with Lp(a)-cholesterol (P = 2.28 × 10(-9)). Conditional analyses revealed that rs3798220 and rs10455872 were associated with Lp(a)-cholesterol levels independent of each other and KIV-2 copy number. Furthermore, we determined for the first time that levels of LPA mRNA were higher in the carriers than non-carriers of rs10455872 (P = 0.0001) and were not different between carriers and non-carriers of rs3798220. Protein levels of apo(a) were higher in the carriers than non-carriers of both rs10455872 and rs3798220. In summary, we identified multiple independent genetic determinants for Lp(a)-cholesterol. These findings provide new insights into Lp(a) regulation. Show less

The triglyceride (TG) response to a high-fat meal (postprandial lipemia, PPL) affects cardiovascular disease risk and is influenced by genes and environment. Genes involved in lipid metabolism have dominated genetic studies of PPL TG response. We sought to elucidate common genetic variants through a genome-wide association (GWA) study in the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN). The GOLDN GWAS discovery sample consisted of 872 participants within families of European ancestry. Genotypes for 2,543,887 variants were measured or imputed from HapMap. Replication of our top results was performed in the Heredity and Phenotype Intervention (HAPI) Heart Study (n = 843). PPL TG response phenotypes were constructed from plasma TG measured at baseline (fasting, 0 hour), 3.5 and 6 hours after a high-fat meal, using a random coefficient regression model. Association analyses were adjusted for covariates and principal components, as necessary, in a linear mixed model using the kinship matrix; additional models further adjusted for fasting TG were also performed. Meta-analysis of the discovery and replication studies (n = 1715) was performed on the top SNPs from GOLDN. GOLDN revealed 111 suggestive (p < 1E-05) associations, with two SNPs meeting GWA significance level (p < 5E-08). Of the two significant SNPs, rs964184 demonstrated evidence of replication (p = 1.20E-03) in the HAPI Heart Study and in a joint analysis, was GWA significant (p = 1.26E-09). Rs964184 has been associated with fasting lipids (TG and HDL) and is near ZPR1 (formerly ZNF259), close to the APOA1/C3/A4/A5 cluster. This association was attenuated upon additional adjustment for fasting TG. This is the first report of a genome-wide significant association with replication for a novel phenotype, namely PPL TG response. Future investigation into response phenotypes is warranted using pathway analyses, or newer genetic technologies such as metabolomics. Show less

Common genetic variants have been recently associated with fasting glucose and insulin levels in white populations. Whether these associations replicate in pre-diabetes is not known. We extended these findings to the Diabetes Prevention Program, a clinical trial in which participants at high risk for diabetes were randomized to placebo, lifestyle modification or metformin for diabetes prevention. We genotyped previously reported polymorphisms (or their proxies) in/near G6PC2, MTNR1B, GCK, DGKB, GCKR, ADCY5, MADD, CRY2, ADRA2A, FADS1, PROX1, SLC2A2, GLIS3, C2CD4B, IGF1, and IRS1 in 3,548 Diabetes Prevention Program participants. We analyzed variants for association with baseline glycemic traits, incident diabetes and their interaction with response to metformin or lifestyle intervention. We replicated associations with fasting glucose at MTNR1B (P<0.001), G6PC2 (P = 0.002) and GCKR (P = 0.001). We noted impaired β-cell function in carriers of glucose-raising alleles at MTNR1B (P<0.001), and an increase in the insulinogenic index for the glucose-raising allele at G6PC2 (P<0.001). The association of MTNR1B with fasting glucose and impaired β-cell function persisted at 1 year despite adjustment for the baseline trait, indicating a sustained deleterious effect at this locus. We also replicated the association of MADD with fasting proinsulin levels (P<0.001). We detected no significant impact of these variants on diabetes incidence or interaction with preventive interventions. The association of several polymorphisms with quantitative glycemic traits is replicated in a cohort of high-risk persons. These variants do not have a detectable impact on diabetes incidence or response to metformin or lifestyle modification in the Diabetes Prevention Program. Show less

Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and stratified for key CKD risk factors. We uncovered 6 new loci in association with estimated glomerular filtration rate (eGFR), the primary clinical measure of CKD, in or near MPPED2, DDX1, SLC47A1, CDK12, CASP9, and INO80. Morpholino knockdown of mpped2 and casp9 in zebrafish embryos revealed podocyte and tubular abnormalities with altered dextran clearance, suggesting a role for these genes in renal function. By providing new insights into genes that regulate renal function, these results could further our understanding of the pathogenesis of CKD. Show less

C-reactive protein (CRP) is a heritable marker of chronic inflammation that is strongly associated with cardiovascular disease. We sought to identify genetic variants that are associated with CRP levels. We performed a genome-wide association analysis of CRP in 66 185 participants from 15 population-based studies. We sought replication for the genome-wide significant and suggestive loci in a replication panel comprising 16 540 individuals from 10 independent studies. We found 18 genome-wide significant loci, and we provided evidence of replication for 8 of them. Our results confirm 7 previously known loci and introduce 11 novel loci that are implicated in pathways related to the metabolic syndrome (APOC1, HNF1A, LEPR, GCKR, HNF4A, and PTPN2) or the immune system (CRP, IL6R, NLRP3, IL1F10, and IRF1) or that reside in regions previously not known to play a role in chronic inflammation (PPP1R3B, SALL1, PABPC4, ASCL1, RORA, and BCL7B). We found a significant interaction of body mass index with LEPR (P<2.9×10(-6)). A weighted genetic risk score that was developed to summarize the effect of risk alleles was strongly associated with CRP levels and explained ≈5% of the trait variance; however, there was no evidence for these genetic variants explaining the association of CRP with coronary heart disease. We identified 18 loci that were associated with CRP levels. Our study highlights immune response and metabolic regulatory pathways involved in the regulation of chronic inflammation. 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

Chronic kidney disease (CKD) is a significant public health problem, and recent genetic studies have identified common CKD susceptibility variants. The CKDGen consortium performed a meta-analysis of genome-wide association data in 67,093 individuals of European ancestry from 20 predominantly population-based studies in order to identify new susceptibility loci for reduced renal function as estimated by serum creatinine (eGFRcrea), serum cystatin c (eGFRcys) and CKD (eGFRcrea < 60 ml/min/1.73 m(2); n = 5,807 individuals with CKD (cases)). Follow-up of the 23 new genome-wide-significant loci (P < 5 x 10(-8)) in 22,982 replication samples identified 13 new loci affecting renal function and CKD (in or near LASS2, GCKR, ALMS1, TFDP2, DAB2, SLC34A1, VEGFA, PRKAG2, PIP5K1B, ATXN2, DACH1, UBE2Q2 and SLC7A9) and 7 loci suspected to affect creatinine production and secretion (CPS1, SLC22A2, TMEM60, WDR37, SLC6A13, WDR72 and BCAS3). These results further our understanding of the biologic mechanisms of kidney function by identifying loci that potentially influence nephrogenesis, podocyte function, angiogenesis, solute transport and metabolic functions of the kidney. 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

To identify loci for age at menarche, we performed a meta-analysis of 32 genome-wide association studies in 87,802 women of European descent, with replication in up to 14,731 women. In addition to the known loci at LIN28B (P = 5.4 × 10⁻⁶⁰) and 9q31.2 (P = 2.2 × 10⁻³³), we identified 30 new menarche loci (all P < 5 × 10⁻⁸) and found suggestive evidence for a further 10 loci (P < 1.9 × 10⁻⁶). The new loci included four previously associated with body mass index (in or near FTO, SEC16B, TRA2B and TMEM18), three in or near other genes implicated in energy homeostasis (BSX, CRTC1 and MCHR2) and three in or near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and gene-set enrichment pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing. Show less

Central abdominal fat is a strong risk factor for diabetes and cardiovascular disease. To identify common variants influencing central abdominal fat, we conducted a two-stage genome-wide association analysis for waist circumference (WC). In total, three loci reached genome-wide significance. In stage 1, 31,373 individuals of Caucasian descent from eight cohort studies confirmed the role of FTO and MC4R and identified one novel locus associated with WC in the neurexin 3 gene [NRXN3 (rs10146997, p = 6.4x10(-7))]. The association with NRXN3 was confirmed in stage 2 by combining stage 1 results with those from 38,641 participants in the GIANT consortium (p = 0.009 in GIANT only, p = 5.3x10(-8) for combined analysis, n = 70,014). Mean WC increase per copy of the G allele was 0.0498 z-score units (0.65 cm). This SNP was also associated with body mass index (BMI) [p = 7.4x10(-6), 0.024 z-score units (0.10 kg/m(2)) per copy of the G allele] and the risk of obesity (odds ratio 1.13, 95% CI 1.07-1.19; p = 3.2x10(-5) per copy of the G allele). The NRXN3 gene has been previously implicated in addiction and reward behavior, lending further evidence that common forms of obesity may be a central nervous system-mediated disorder. Our findings establish that common variants in NRXN3 are associated with WC, BMI, and obesity. Show less

To identify genetic variants influencing plasma lipid concentrations, we first used genotype imputation and meta-analysis to combine three genome-wide scans totaling 8,816 individuals and comprising 6,068 individuals specific to our study (1,874 individuals from the FUSION study of type 2 diabetes and 4,184 individuals from the SardiNIA study of aging-associated variables) and 2,758 individuals from the Diabetes Genetics Initiative, reported in a companion study in this issue. We subsequently examined promising signals in 11,569 additional individuals. Overall, we identify strongly associated variants in eleven loci previously implicated in lipid metabolism (ABCA1, the APOA5-APOA4-APOC3-APOA1 and APOE-APOC clusters, APOB, CETP, GCKR, LDLR, LPL, LIPC, LIPG and PCSK9) and also in several newly identified loci (near MVK-MMAB and GALNT2, with variants primarily associated with high-density lipoprotein (HDL) cholesterol; near SORT1, with variants primarily associated with low-density lipoprotein (LDL) cholesterol; near TRIB1, MLXIPL and ANGPTL3, with variants primarily associated with triglycerides; and a locus encompassing several genes near NCAN, with variants strongly associated with both triglycerides and LDL cholesterol). Notably, the 11 independent variants associated with increased LDL cholesterol concentrations in our study also showed increased frequency in a sample of coronary artery disease cases versus controls. Show less

Apolipoprotein C-III (apoC-III) inhibits triglyceride hydrolysis and has been implicated in coronary artery disease. Through a genome-wide association study, we have found that about 5% of the Lancaster Amish are heterozygous carriers of a null mutation (R19X) in the gene encoding apoC-III (APOC3) and, as a result, express half the amount of apoC-III present in noncarriers. Mutation carriers compared with noncarriers had lower fasting and postprandial serum triglycerides, higher levels of HDL-cholesterol and lower levels of LDL-cholesterol. Subclinical atherosclerosis, as measured by coronary artery calcification, was less common in carriers than noncarriers, which suggests that lifelong deficiency of apoC-III has a cardioprotective effect. Show less

Alteration of single genes involved in nutrient and lipoprotein metabolism increases longevity in several animal models. Because exceptional longevity in humans is familial, it is likely that polymorphisms in genes favorably influence certain phenotypes and increase the likelihood of exceptional longevity. A group of Ashkenazi Jewish centenarians (n = 213), their offspring (n = 216), and an age-matched Ashkenazi control group (n = 258) were genotyped for 66 polymorphisms in 36 candidate genes related to cardiovascular disease (CVD). These genes were tested for association with serum lipoprotein levels and particle sizes, apolipoprotein A1, B, and C-3 levels and with outcomes of hypertension, insulin resistance, and mortality. The prevalence of homozygosity for the -641C allele in the APOC3 promoter (rs2542052) was higher in centenarians (25%) and their offspring (20%) than in controls (10%) (p = 0.0001 and p = 0.001, respectively). This genotype was associated with significantly lower serum levels of APOC3 and a favorable pattern of lipoprotein levels and sizes. We found a lower prevalence of hypertension and greater insulin sensitivity in the -641C homozygotes, suggesting a protective effect against CVD and the metabolic syndrome. Finally, in a prospectively studied cohort, a significant survival advantage was demonstrated in those with the favorable -641C homozygote (p < 0.0001). Homozygosity for the APOC3 -641C allele is associated with a favorable lipoprotein profile, cardiovascular health, insulin sensitivity, and longevity. Because modulation of lipoproteins is also seen in genetically altered longevity models, it may be a common pathway influencing lifespan from nematodes to humans. Show less