👤 R J F Loos

To delineate organ-specific and systemic drivers of metabolic dysfunction-associated steatotic liver disease (MASLD), we applied integrative causal inference across clinical, imaging, and proteomic domains in individuals with and without type 2 diabetes (T2D). Bayesian network analyses and complementary two-sample Mendelian randomization were used to quantify causal pathways linking adipose distribution, glycemia, and insulin dynamics with liver fat in the IMI-DIRECT prospective cohort study. Data included frequently sampled metabolic challenge tests, MRI-derived abdominal and hepatic fat content, serological biomarkers, and Olink plasma proteomics from 331 adults with new-onset T2D and 964 adults without diabetes, with harmonized protocols enabling replication. High basal insulin secretion rate (BasalISR), estimated via C-peptide deconvolution, emerged as the primary potential causal driver of liver fat accumulation in both cohorts. BasalISR, a clearance-independent measure of β-cell insulin output distinct from peripheral insulin levels, was independently linked to hepatic steatosis. Visceral adipose tissue exhibited bidirectional associations with liver fat, suggesting a self-reinforcing metabolic loop. Of 446 analyzed proteins, 34 mapped to these metabolic networks (27 in the non-diabetes network, 18 in the T2D network, and 11 shared). Key proteins directly associated with liver fat included GUSB, ALDH1A1, LPL, IGFBP1/2, CTSD, HMOX1, FGF21, AGRP, and ACE2. Sex-stratified analyses identified GUSB in females and LEP in males as the strongest protein predictors of liver fat. BasalISR may better capture early β-cell-driven disturbances contributing to MASLD. These findings outline a multifactorial, sex- and disease stage-specific proteo-metabolic architecture of hepatic steatosis and identify potential biomarkers or therapeutic targets. Show less

The objective of this study is to evaluate obesity-related genetic factors in relation to excess consumption and assess if food cues modify associations. Children (9-12 years) completed a randomized crossover experiment. During two visits, children ate a preload and then snacks ad libitum while watching television, embedded with food or non-food advertisements to assess eating in the absence of hunger (EAH). Primary exposures were obesity-associated genotypes, FTO rs9939609 and MC4R rs571312, and a paediatric-specific polygenic risk score (PRS). Outcomes included consumption of all snacks (total EAH) and gummy candy only (gummy candy EAH). Linear mixed-effects models tested whether genetic exposures related to EAH outcomes. We tested for effect modification by food cues using multiplicative interaction terms. Among 177 children, each FTO risk allele was associated with a 30% increase in gummy candy EAH (p = 0.025) in adjusted models. Food cue exposure exacerbated associations between the FTO variant with gummy candy EAH (p = 0.046). No statistically significant associations were found between MC4R and EAH. The results suggest children with the FTO rs9939609 risk allele may be predisposed to excess consumption of candy and that this association may be exacerbated by food cues. Show less

To delineate organ-specific and systemic drivers of metabolic dysfunction-associated steatotic liver disease (MASLD), we applied integrative causal inference across clinical, imaging, and proteomic domains in individuals with and without type 2 diabetes (T2D). We used Bayesian network analyses to quantify causal pathways linking adipose distribution, glycemia, and insulin dynamics with fatty liver using data from the IMI-DIRECT prospective cohort study. Measurements were made of glucose and insulin dynamics (using frequently-sampled metabolic challenge tests), MRI-derived abdominal and liver fat content, serological biomarkers, and Olink plasma proteomics from 331 adults with new-onset T2D and 964 adults free from diabetes at enrolment. The common protocols used in these two cohorts provided the opportunity for replication analyses to be performed. When the direction of the effect could not be determined with high probability through Bayesian networks, complementary two-sample Mendelian randomization (MR) was employed. High basal insulin secretion rate (BasalISR) was identified as the primary causal driver of liver fat accumulation in both diabetes and non-diabetes. Excess visceral adipose tissue (VAT) was bidirectionally associated with liver fat, indicating a self-reinforcing metabolic loop. Basal insulin clearance (Clinsb) worsened as a consequence of liver fat accumulation to a greater degree before the onset of T2D. Out of 446 analysed proteins, 34 mapped to these metabolic networks and 27 were identified in the non-diabetes network, 18 in the diabetes network, and 11 were common between the two networks. Key proteins directly associated with liver fat included GUSB, ALDH1A1, LPL, IGFBP1/2, CTSD, HMOX1, FGF21, AGRP, and ACE2. Sex-stratified analyses revealed distinct proteomic drivers: GUSB and LEP were most predictive of liver fat in females and males, respectively. Basal insulin hypersecretion is a modifiable, causal driver of MASLD, particularly prior to glycaemic decompensation. Our findings highlight a multifactorial, sex- and disease-stage-specific proteo-metabolic architecture of hepatic steatosis. Proteins such as GUSB, ALDH1A1, LPL, and IGFBPs warrant further investigation as potential biomarkers or therapeutic targets for MASLD prevention and treatment. Show less

Atrial fibrillation (AF) is a prevalent and morbid abnormality of the heart rhythm with a strong genetic component. Here, we meta-analyzed genome and exome sequencing data from 36 studies that included 52,416 AF cases and 277,762 controls. In burden tests of rare coding variation, we identified novel associations between AF and the genes MYBPC3, LMNA, PKP2, FAM189A2 and KDM5B. We further identified associations between AF and rare structural variants owing to deletions in CTNNA3 and duplications of GATA4. We broadly replicated our findings in independent samples from MyCode, deCODE and UK Biobank. Finally, we found that CRISPR knockout of KDM5B in stem-cell-derived atrial cardiomyocytes led to a shortening of the action potential duration and widespread transcriptomic dysregulation of genes relevant to atrial homeostasis and conduction. Our results highlight the contribution of rare coding and structural variants to AF, including genetic links between AF and cardiomyopathies, and expand our understanding of the rare variant architecture for this common arrhythmia. Show less

To broaden our understanding of bradyarrhythmias and conduction disease, we performed common variant genome-wide association analyses in up to 1.3 million individuals and rare variant burden testing in 460,000 individuals for sinus node dysfunction (SND), distal conduction disease (DCD) and pacemaker (PM) implantation. We identified 13, 31 and 21 common variant loci for SND, DCD and PM, respectively. Four well-known loci (SCN5A/SCN10A, CCDC141, TBX20 and CAMK2D) were shared for SND and DCD, while others were more specific for SND or DCD. SND and DCD showed a moderate genetic correlation (r Show less

Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals (P < 5 × 10 Show less

Obesity is a major risk factor for a myriad of diseases, affecting >600 million people worldwide. Genome-wide association studies (GWASs) have identified hundreds of genetic variants that influence body mass index (BMI), a commonly used metric to assess obesity risk. Most variants are non-coding and likely act through regulating genes nearby. Here, we apply multiple computational methods to prioritize the likely causal gene(s) within each of the 536 previously reported GWAS-identified BMI-associated loci. We performed summary-data-based Mendelian randomization (SMR), FINEMAP, DEPICT, MAGMA, transcriptome-wide association studies (TWASs), mutation significance cutoff (MSC), polygenic priority score (PoPS), and the nearest gene strategy. Results of each method were weighted based on their success in identifying genes known to be implicated in obesity, ranking all prioritized genes according to a confidence score (minimum: 0; max: 28). We identified 292 high-scoring genes (≥11) in 264 loci, including genes known to play a role in body weight regulation (e.g., DGKI, ANKRD26, MC4R, LEPR, BDNF, GIPR, AKT3, KAT8, MTOR) and genes related to comorbidities (e.g., FGFR1, ISL1, TFAP2B, PARK2, TCF7L2, GSK3B). For most of the high-scoring genes, however, we found limited or no evidence for a role in obesity, including the top-scoring gene BPTF. Many of the top-scoring genes seem to act through a neuronal regulation of body weight, whereas others affect peripheral pathways, including circadian rhythm, insulin secretion, and glucose and carbohydrate homeostasis. The characterization of these likely causal genes can increase our understanding of the underlying biology and offer avenues to develop therapeutics for weight loss. Show less

To test associations of candidate obesity-related single nucleotide polymorphisms (SNPs) and obesity polygenic risk scores (PRS) with neural reward reactivity to food cues. After consuming a pre-load meal, 9-12-year-old children completed a functional magnetic resonance imaging (fMRI) paradigm with exposure to food and non-food commercials. Genetic exposures included Analyses included 151 children ( Children genetically at risk for obesity, as indicated by the https://clinicaltrials.gov/study/NCT03766191, identifier NCT03766191. Show less

Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals ( Show less

MYBPC3 is the most frequently affected gene in hypertrophic cardiomyopathy (HCM), which is an autosomal-dominant cardiac disease caused by mutations in sarcomeric proteins. Bi-allelic truncating MYBPC3 mutations are associated with severe forms of neonatal cardiomyopathy. We reprogrammed skin fibroblasts from a HCM patient carrying a heterozygous MYBPC3 truncating mutation into human induced pluripotent stem cells (iPSC) and used CRISPR/Cas9 to generate bi-allelic MYBPC3 truncating mutation and isogenic control hiPSC lines. All lines expressed pluripotency markers, had normal karyotype and differentiated into endoderm, ectoderm and cardiomyocytes in vitro. This set of three lines provides a useful tool to study HCM pathomechanisms. Show less

Melanocortin 4 receptor (MC4R) deficiency, caused by mutations in MC4R, is the most common cause of monogenic forms of obesity. However, these mutations have often been identified in small-scale, case-focused studies. Here, we assess the penetrance of previously reported MC4R mutations at a population level. Furthermore, we examine why some carriers of pathogenic mutations remain of normal weight, to gain insight into the mechanisms that control body weight. We identified 59 known obesity-increasing mutations in MC4R from the Human Gene Mutation Database (HGMD) and Clinvar. We assessed their penetrance and effect on obesity (body mass index [BMI] ≥ 30 kg/m2) in >450,000 individuals (age 40-69 years) of the UK Biobank, a population-based cohort study. Of these 59 mutations, only 11 had moderate-to-high penetrance and increased the odds of obesity by more than 2-fold. We subsequently focused on these 11 mutations and examined differences between carriers of normal weight and carriers with obesity. Twenty-eight of the 182 carriers of these 11 mutations were of normal weight. Body composition of carriers of normal weight was similar to noncarriers of normal weight, whereas among individuals with obesity, carriers had a somewhat higher BMI than noncarriers (1.44 ± 0.07 standard deviation scores [SDSs] ± standard error [SE] versus 1.29 ± 0.001, P = 0.03), because of greater lean mass (1.44 ± 0.09 versus 1.15 ± 0.002, P = 0.002). Carriers of normal weight more often reported that, already at age 10 years, their body size was below average or average (72%) compared with carriers with obesity (48%) (P = 0.01). To assess the polygenic contribution to body weight in carriers of normal weight and carriers with obesity, we calculated a genome-wide polygenic risk score for BMI (PRSBMI). The PRSBMI of carriers of normal weight (PRSBMI = -0.64 ± 0.18) was significantly lower than of carriers with obesity (0.40 ± 0.11; P = 1.7 × 10-6), and tended to be lower than that of noncarriers of normal weight (-0.29 ± 0.003; P = 0.05). Among carriers, those with a low PRSBMI (bottom quartile) have an approximately 5-kg/m2 lower BMI (approximately 14 kg of body weight for a 1.7-m-tall person) than those with a high PRS (top quartile). Because the UK Biobank population is healthier than the general population in the United Kingdom, penetrance may have been somewhat underestimated. We showed that large-scale data are needed to validate the impact of mutations observed in small-scale and case-focused studies. Furthermore, we observed that despite the key role of MC4R in obesity, the effects of pathogenic MC4R mutations may be countered, at least in part, by a low polygenic risk potentially representing other innate mechanisms implicated in body weight regulation. Show less

Lean body mass (LM) plays an important role in mobility and metabolic function. We previously identified five loci associated with LM adjusted for fat mass in kilograms. Such an adjustment may reduce the power to identify genetic signals having an association with both lean mass and fat mass. To determine the impact of different fat mass adjustments on genetic architecture of LM and identify additional LM loci. We performed genome-wide association analyses for whole-body LM (20 cohorts of European ancestry with n = 38,292) measured using dual-energy X-ray absorptiometry) or bioelectrical impedance analysis, adjusted for sex, age, age2, and height with or without fat mass adjustments (Model 1 no fat adjustment; Model 2 adjustment for fat mass as a percentage of body mass; Model 3 adjustment for fat mass in kilograms). Seven single-nucleotide polymorphisms (SNPs) in separate loci, including one novel LM locus (TNRC6B), were successfully replicated in an additional 47,227 individuals from 29 cohorts. Based on the strengths of the associations in Model 1 vs Model 3, we divided the LM loci into those with an effect on both lean mass and fat mass in the same direction and refer to those as "sumo wrestler" loci (FTO and MC4R). In contrast, loci with an impact specifically on LM were termed "body builder" loci (VCAN and ADAMTSL3). Using existing available genome-wide association study databases, LM increasing alleles of SNPs in sumo wrestler loci were associated with an adverse metabolic profile, whereas LM increasing alleles of SNPs in "body builder" loci were associated with metabolic protection. In conclusion, we identified one novel LM locus (TNRC6B). Our results suggest that a genetically determined increase in lean mass might exert either harmful or protective effects on metabolic traits, depending on its relation to fat mass. 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

Genome-wide association studies have so far identified 56 loci associated with risk of coronary artery disease (CAD). Many CAD loci show pleiotropy; that is, they are also associated with other diseases or traits. This study sought to systematically test if genetic variants identified for non-CAD diseases/traits also associate with CAD and to undertake a comprehensive analysis of the extent of pleiotropy of all CAD loci. In discovery analyses involving 42,335 CAD cases and 78,240 control subjects we tested the association of 29,383 common (minor allele frequency >5%) single nucleotide polymorphisms available on the exome array, which included a substantial proportion of known or suspected single nucleotide polymorphisms associated with common diseases or traits as of 2011. Suggestive association signals were replicated in an additional 30,533 cases and 42,530 control subjects. To evaluate pleiotropy, we tested CAD loci for association with cardiovascular risk factors (lipid traits, blood pressure phenotypes, body mass index, diabetes, and smoking behavior), as well as with other diseases/traits through interrogation of currently available genome-wide association study catalogs. We identified 6 new loci associated with CAD at genome-wide significance: on 2q37 (KCNJ13-GIGYF2), 6p21 (C2), 11p15 (MRVI1-CTR9), 12q13 (LRP1), 12q24 (SCARB1), and 16q13 (CETP). Risk allele frequencies ranged from 0.15 to 0.86, and odds ratio per copy of the risk allele ranged from 1.04 to 1.09. Of 62 new and known CAD loci, 24 (38.7%) showed statistical association with a traditional cardiovascular risk factor, with some showing multiple associations, and 29 (47%) showed associations at p < 1 × 10 We identified 6 loci associated with CAD at genome-wide significance. Several CAD loci show substantial pleiotropy, which may help us understand the mechanisms by which these loci affect CAD risk. Show less

Genome-wide association studies have identified over 150 loci associated with lipid traits, however, no large-scale studies exist for Hispanics and other minority populations. Additionally, the genetic architecture of lipid-influencing loci remains largely unknown. We performed one of the most racially/ethnically diverse fine-mapping genetic studies of HDL-C, LDL-C, and triglycerides to-date using SNPs on the MetaboChip array on 54,119 individuals: 21,304 African Americans, 19,829 Hispanic Americans, 12,456 Asians, and 530 American Indians. The majority of signals found in these groups generalize to European Americans. While we uncovered signals unique to racial/ethnic populations, we also observed systematically consistent lipid associations across these groups. In African Americans, we identified three novel signals associated with HDL-C (LPL, APOA5, LCAT) and two associated with LDL-C (ABCG8, DHODH). In addition, using this population, we refined the location for 16 out of the 58 known MetaboChip lipid loci. These results can guide tailored screening efforts, reveal population-specific responses to lipid-lowering medications, and aid in the development of new targeted drug therapies. Show less

We performed genome-wide meta-analysis of lipid traits on three samples of Mexican and Mexican American ancestry comprising 4,383 individuals, and followed up significant and highly suggestive associations in three additional Hispanic samples comprising 7,876 individuals. Genome-wide significant signals were observed in or near CELSR2, ZNF259/APOA5, KANK2/DOCK6 and NCAN/MAU2 for total cholesterol, LPL, ABCA1, ZNF259/APOA5, LIPC and CETP for HDL cholesterol, CELSR2, APOB and NCAN/MAU2 for LDL cholesterol, and GCKR, TRIB1, ZNF259/APOA5 and NCAN/MAU2 for triglycerides. Linkage disequilibrium and conditional analyses indicate that signals observed at ABCA1 and LIPC for HDL cholesterol and NCAN/MAU2 for triglycerides are independent of previously reported lead SNP associations. Analyses of lead SNPs from the European Global Lipids Genetics Consortium (GLGC) dataset in our Hispanic samples show remarkable concordance of direction of effects as well as strong correlation in effect sizes. A meta-analysis of the European GLGC and our Hispanic datasets identified five novel regions reaching genome-wide significance: two for total cholesterol (FN1 and SAMM50), two for HDL cholesterol (LOC100996634 and COPB1) and one for LDL cholesterol (LINC00324/CTC1/PFAS). The top meta-analysis signals were found to be enriched for SNPs associated with gene expression in a tissue-specific fashion, suggesting an enrichment of tissue-specific function in lipid-associated loci. Show less

White blood cells play diverse roles in innate and adaptive immunity. Genetic association analyses of phenotypic variation in circulating white blood cell (WBC) counts from large samples of otherwise healthy individuals can provide insights into genes and biologic pathways involved in production, differentiation, or clearance of particular WBC lineages (myeloid, lymphoid) and also potentially inform the genetic basis of autoimmune, allergic, and blood diseases. We performed an exome array-based meta-analysis of total WBC and subtype counts (neutrophils, monocytes, lymphocytes, basophils, and eosinophils) in a multi-ancestry discovery and replication sample of ∼157,622 individuals from 25 studies. We identified 16 common variants (8 of which were coding variants) associated with one or more WBC traits, the majority of which are pleiotropically associated with autoimmune diseases. Based on functional annotation, these loci included genes encoding surface markers of myeloid, lymphoid, or hematopoietic stem cell differentiation (CD69, CD33, CD87), transcription factors regulating lineage specification during hematopoiesis (ASXL1, IRF8, IKZF1, JMJD1C, ETS2-PSMG1), and molecules involved in neutrophil clearance/apoptosis (C10orf54, LTA), adhesion (TNXB), or centrosome and microtubule structure/function (KIF9, TUBD1). Together with recent reports of somatic ASXL1 mutations among individuals with idiopathic cytopenias or clonal hematopoiesis of undetermined significance, the identification of a common regulatory 3' UTR variant of ASXL1 suggests that both germline and somatic ASXL1 mutations contribute to lower blood counts in otherwise asymptomatic individuals. These association results shed light on genetic mechanisms that regulate circulating WBC counts and suggest a prominent shared genetic architecture with inflammatory and autoimmune diseases. Show less

Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome-wide (GW) meta-analysis of predominately regular-type coffee consumption (cups per day) among up to 91,462 coffee consumers of European ancestry with top single-nucleotide polymorphisms (SNPs) followed-up in ~30 062 and 7964 coffee consumers of European and African-American ancestry, respectively. Studies from both stages were combined in a trans-ethnic meta-analysis. Confirmed loci were examined for putative functional and biological relevance. Eight loci, including six novel loci, met GW significance (log10Bayes factor (BF)>5.64) with per-allele effect sizes of 0.03-0.14 cups per day. Six are located in or near genes potentially involved in pharmacokinetics (ABCG2, AHR, POR and CYP1A2) and pharmacodynamics (BDNF and SLC6A4) of caffeine. Two map to GCKR and MLXIPL genes related to metabolic traits but lacking known roles in coffee consumption. Enhancer and promoter histone marks populate the regions of many confirmed loci and several potential regulatory SNPs are highly correlated with the lead SNP of each. SNP alleles near GCKR, MLXIPL, BDNF and CYP1A2 that were associated with higher coffee consumption have previously been associated with smoking initiation, higher adiposity and fasting insulin and glucose but lower blood pressure and favorable lipid, inflammatory and liver enzyme profiles (P<5 × 10(-8)).Our genetic findings among European and African-American adults reinforce the role of caffeine in mediating habitual coffee consumption and may point to molecular mechanisms underlying inter-individual variability in pharmacological and health effects of coffee. Show less

Genome-wide association studies (GWAS) of BMI are mostly undertaken under the assumption that "kg/m(2) " is an index of weight fully adjusted for height, but in general this is not true. The aim here was to assess the contribution of common genetic variation to a adjusted version of that phenotype which appropriately accounts for covariation in height in children. A GWAS of height-adjusted BMI (BMI[x] = weight/height(x) ), calculated to be uncorrelated with height, in 5809 participants (mean age 9.9 years) from the Avon Longitudinal Study of Parents and Children (ALSPAC) was performed. GWAS based on BMI[x] yielded marked differences in genomewide results profile. SNPs in ADCY3 (adenylate cyclase 3) were associated at genome-wide significance level (rs11676272 (0.28 kg/m(3.1) change per allele G (0.19, 0.38), P = 6 × 10(-9) ). In contrast, they showed marginal evidence of association with conventional BMI [rs11676272 (0.25 kg/m(2) (0.15, 0.35), P = 6 × 10(-7) )]. Results were replicated in an independent sample, the Generation R study. Analysis of BMI[x] showed differences to that of conventional BMI. The association signal at ADCY3 appeared to be driven by a missense variant and it was strongly correlated with expression of this gene. Our work highlights the importance of well understood phenotype use (and the danger of convention) in characterising genetic contributions to complex traits. Show less

Plasma triglyceride levels are heritable and are correlated with the risk of coronary heart disease. Sequencing of the protein-coding regions of the human genome (the exome) has the potential to identify rare mutations that have a large effect on phenotype. We sequenced the protein-coding regions of 18,666 genes in each of 3734 participants of European or African ancestry in the Exome Sequencing Project. We conducted tests to determine whether rare mutations in coding sequence, individually or in aggregate within a gene, were associated with plasma triglyceride levels. For mutations associated with triglyceride levels, we subsequently evaluated their association with the risk of coronary heart disease in 110,970 persons. An aggregate of rare mutations in the gene encoding apolipoprotein C3 (APOC3) was associated with lower plasma triglyceride levels. Among the four mutations that drove this result, three were loss-of-function mutations: a nonsense mutation (R19X) and two splice-site mutations (IVS2+1G→A and IVS3+1G→T). The fourth was a missense mutation (A43T). Approximately 1 in 150 persons in the study was a heterozygous carrier of at least one of these four mutations. Triglyceride levels in the carriers were 39% lower than levels in noncarriers (P<1×10(-20)), and circulating levels of APOC3 in carriers were 46% lower than levels in noncarriers (P=8×10(-10)). The risk of coronary heart disease among 498 carriers of any rare APOC3 mutation was 40% lower than the risk among 110,472 noncarriers (odds ratio, 0.60; 95% confidence interval, 0.47 to 0.75; P=4×10(-6)). Rare mutations that disrupt APOC3 function were associated with lower levels of plasma triglycerides and APOC3. Carriers of these mutations were found to have a reduced risk of coronary heart disease. (Funded by the National Heart, Lung, and Blood Institute and others.). Show less

Metabolic syndrome (MetS) has become a health and financial burden worldwide. The MetS definition captures clustering of risk factors that predict higher risk for diabetes mellitus and cardiovascular disease. Our study hypothesis is that additional to genes influencing individual MetS risk factors, genetic variants exist that influence MetS and inflammatory markers forming a predisposing MetS genetic network. To test this hypothesis a staged approach was undertaken. (a) We analyzed 17 metabolic and inflammatory traits in more than 85,500 participants from 14 large epidemiological studies within the Cross Consortia Pleiotropy Group. Individuals classified with MetS (NCEP definition), versus those without, showed on average significantly different levels for most inflammatory markers studied. (b) Paired average correlations between 8 metabolic traits and 9 inflammatory markers from the same studies as above, estimated with two methods, and factor analyses on large simulated data, helped in identifying 8 combinations of traits for follow-up in meta-analyses, out of 130,305 possible combinations between metabolic traits and inflammatory markers studied. (c) We performed correlated meta-analyses for 8 metabolic traits and 6 inflammatory markers by using existing GWAS published genetic summary results, with about 2.5 million SNPs from twelve predominantly largest GWAS consortia. These analyses yielded 130 unique SNPs/genes with pleiotropic associations (a SNP/gene associating at least one metabolic trait and one inflammatory marker). Of them twenty-five variants (seven loci newly reported) are proposed as MetS candidates. They map to genes MACF1, KIAA0754, GCKR, GRB14, COBLL1, LOC646736-IRS1, SLC39A8, NELFE, SKIV2L, STK19, TFAP2B, BAZ1B, BCL7B, TBL2, MLXIPL, LPL, TRIB1, ATXN2, HECTD4, PTPN11, ZNF664, PDXDC1, FTO, MC4R and TOMM40. Based on large data evidence, we conclude that inflammation is a feature of MetS and several gene variants show pleiotropic genetic associations across phenotypes and might explain a part of MetS correlated genetic architecture. These findings warrant further functional investigation. Show less

Recent large-scale genome-wide association studies have identified multiple loci robustly associated with BMI, predominantly in European ancestry (EA) populations. However, associations of these loci with obesity and related traits have not been well described in Chinese Hans. This study aimed to investigate whether BMI-associated loci are, individually and collectively, associated with adiposity-related traits and obesity in Chinese Hans and whether these associations are modified by physical activity (PA). We genotyped 28 BMI-associated single nucleotide polymorphisms (SNPs) in a population-based cohort including 2,894 unrelated Han Chinese. Genetic risk score (GRS), EA and East Asian ancestry (EAA) GRSs were calculated by adding BMI-increasing alleles based on all, EA and EAA identified SNPs, respectively. Interactions of GRS and PA were examined by including the interaction-term in the regression model. Individually, 26 of 28 SNPs showed directionally consistent effects on BMI, and associations of four loci (TMEM18, PCSK1, BDNF and MAP2K5) reached nominal significance (P<0.05). The GRS was associated with increased BMI, trunk fat and body fat percentages; and increased risk of obesity and overweight (all P<0.05). Effect sizes (0.11 vs. 0.17 kg/m2) and explained variance (0.90% vs. 1.45%) of GRS for BMI tended to be lower in Chinese Hans than in Europeans. The EA GRS and EAA GRS were associated with 0.11 and 0.13 kg/m2 higher BMI, respectively. In addition, we found that PA attenuated the effect of the GRS on BMI (Pinteraction = 0.022). Our observations suggest that the combined effect of obesity-susceptibility loci on BMI tended to be lower in Han Chinese than in EA. The overall, EA and EAA GRSs exert similar effects on adiposity traits. Genetic predisposition to increased BMI is attenuated by PA in this population of Han Chinese. Show less

We analysed single nucleotide polymorphisms (SNPs) tagging the genetic variability of six candidate genes (ATF6, FABP1, LPIN2, LPIN3, MLXIPL and MTTP) involved in the regulation of hepatic lipid metabolism, an important regulatory site of energy balance for associations with body mass index (BMI) and changes in weight and waist circumference. We also investigated effect modification by sex and dietary intake. Data of 6,287 individuals participating in the European prospective investigation into cancer and nutrition were included in the analyses. Data on weight and waist circumference were followed up for 6.9 ± 2.5 years. Association of 69 tagSNPs with baseline BMI and annual changes in weight as well as waist circumference were investigated using linear regression analysis. Interactions with sex, GI and intake of carbohydrates, fat as well as saturated, monounsaturated and polyunsaturated fatty acids were examined by including multiplicative SNP-covariate terms into the regression model. Neither baseline BMI nor annual weight or waist circumference changes were significantly associated with variation in the selected genes in the entire study population after correction for multiple testing. One SNP (rs1164) in LPIN2 appeared to be significantly interacting with sex (p = 0.0003) and was associated with greater annual weight gain in men (56.8 ± 23.7 g/year per allele, p = 0.02) than in women (-25.5 ± 19.8 g/year per allele, p = 0.2). With respect to gene-nutrient interaction, we could not detect any significant interactions when accounting for multiple testing. Therefore, out of our six candidate genes, LPIN2 may be considered as a candidate for further studies. Show less

Genetic pleiotropy may contribute to the clustering of obesity and metabolic conditions. We assessed whether genetic variants that are robustly associated with BMI and waist-to-hip ratio (WHR) also influence metabolic and cardiovascular traits, independently of obesity-related traits, in meta-analyses of up to 37,874 individuals from six European population-based studies. We examined associations of 32 BMI and 14 WHR loci, individually and combined in two genetic predisposition scores (GPSs), with glycaemic traits, blood lipids and BP, with and without adjusting for BMI and/or WHR. We observed significant associations of BMI-increasing alleles at five BMI loci with lower levels of 2 h glucose (RBJ [also known as DNAJC27], QPTCL: effect sizes -0.068 and -0.107 SD, respectively), HDL-cholesterol (SLC39A8: -0.065 SD, MTCH2: -0.039 SD), and diastolic BP (SLC39A8: -0.069 SD), and higher and lower levels of LDL- and total cholesterol (QPTCL: 0.041 and 0.042 SDs, respectively, FLJ35779 [also known as POC5]: -0.042 and -0.041 SDs, respectively) (all p < 2.4 × 10(-4)), independent of BMI. The WHR-increasing alleles at two WHR loci were significantly associated with higher proinsulin (GRB14: 0.069 SD) and lower fasting glucose levels (CPEB4: -0.049 SD), independent of BMI and WHR. A higher GPS-BMI was associated with lower systolic BP (-0.005 SD), diastolic BP (-0.006 SD) and 2 h glucose (-0.013 SD), while a higher GPS-WHR was associated with lower HDL-cholesterol (-0.015 SD) and higher triacylglycerol levels (0.014 SD) (all p < 2.9 × 10(-3)), independent of BMI and/or WHR. These pleiotropic effects of obesity-susceptibility loci provide novel insights into mechanisms that link obesity with metabolic abnormalities. Show less

Differences in the interindividual response to dietary intervention could be modified by genetic variation in nutrient-sensitive genes. This study examined single nucleotide polymorphisms (SNPs) in presumed nutrient-sensitive candidate genes for obesity and obesity-related diseases for main and dietary interaction effects on weight, waist circumference, and fat mass regain over 6 mo. In total, 742 participants who had lost ≥ 8% of their initial body weight were randomly assigned to follow 1 of 5 different ad libitum diets with different glycemic indexes and contents of dietary protein. The SNP main and SNP-diet interaction effects were analyzed by using linear regression models, corrected for multiple testing by using Bonferroni correction and evaluated by using quantile-quantile (Q-Q) plots. After correction for multiple testing, none of the SNPs were significantly associated with weight, waist circumference, or fat mass regain. Q-Q plots showed that ALOX5AP rs4769873 showed a higher observed than predicted P value for the association with less waist circumference regain over 6 mo (-3.1 cm/allele; 95% CI: -4.6, -1.6; P/Bonferroni-corrected P = 0.000039/0.076), independently of diet. Additional associations were identified by using Q-Q plots for SNPs in ALOX5AP, TNF, and KCNJ11 for main effects; in LPL and TUB for glycemic index interaction effects on waist circumference regain; in GHRL, CCK, MLXIPL, and LEPR on weight; in PPARC1A, PCK2, ALOX5AP, PYY, and ADRB3 on waist circumference; and in PPARD, FABP1, PLAUR, and LPIN1 on fat mass regain for dietary protein interaction. The observed effects of SNP-diet interactions on weight, waist, and fat mass regain suggest that genetic variation in nutrient-sensitive genes can modify the response to diet. This trial was registered at clinicaltrials.gov as NCT00390637. Show less

Longitudinal growth associations with genetic variants identified for adult BMI may provide insights into the timing of obesity susceptibility. The objective was to explore associations of known BMI loci with measures of body size from birth to adulthood. A total of 2537 individuals from a longitudinal British birth cohort were genotyped for 11 genetic variants robustly associated with adult BMI (in/near FTO, MC4R, TMEM18, GNPDA2, KCTD15, NEGR1, BDNF, ETV5, SEC16B, SH2B1, and MTCH2). We derived an obesity-risk-allele score, comprising the sum of BMI-increasing alleles in each individual, and examined this for an association with birth weight and repeated measures of weight, height, and BMI SD scores (SDS) at 11 time points between ages 2 and 53 y. The obesity-risk-allele score showed borderline significant association with birth weight (0.019 SDS/allele; P = 0.05) and was more clearly associated with higher weight and BMI at all time points between ages 2 and 53 y; the strongest associations with weight occurred at ages 11 and 20 y (both 0.056 SDS/allele). In longitudinal analyses, the score was positively associated with weight gain only between birth and 11 y (0.003 SDS/allele per year; 95% CI: 0.001, 0.004; P = 0.001). The risk-allele score was associated with taller height at 7 y (0.031 SDS/allele; P = 0.002) and greater height gains between 2 and 7 y (0.007 SDS/allele per year; P < 0.001), but not with adult height (P = 0.5). The combined effect of adult obesity susceptibility variants on weight gain was confined to childhood. These variants conferred a faster tempo of height growth that was evident before the pubertal years. Show less

Concentrations of liver enzymes in plasma are widely used as indicators of liver disease. We carried out a genome-wide association study in 61,089 individuals, identifying 42 loci associated with concentrations of liver enzymes in plasma, of which 32 are new associations (P = 10(-8) to P = 10(-190)). We used functional genomic approaches including metabonomic profiling and gene expression analyses to identify probable candidate genes at these regions. We identified 69 candidate genes, including genes involved in biliary transport (ATP8B1 and ABCB11), glucose, carbohydrate and lipid metabolism (FADS1, FADS2, GCKR, JMJD1C, HNF1A, MLXIPL, PNPLA3, PPP1R3B, SLC2A2 and TRIB1), glycoprotein biosynthesis and cell surface glycobiology (ABO, ASGR1, FUT2, GPLD1 and ST3GAL4), inflammation and immunity (CD276, CDH6, GCKR, HNF1A, HPR, ITGA1, RORA and STAT4) and glutathione metabolism (GSTT1, GSTT2 and GGT), as well as several genes of uncertain or unknown function (including ABHD12, EFHD1, EFNA1, EPHA2, MICAL3 and ZNF827). Our results provide new insight into genetic mechanisms and pathways influencing markers of liver function. Show less

Recent genome-wide association studies have identified a number of common variants associated with fasting glucose homeostasis and type 2 diabetes in populations of European origin. This is a replication study to examine whether such associations are also observed in Chinese Hans. We genotyped nine variants in or near MADD, ADRA2A, CRY2, GLIS3, PROX1, FADS1, C2CD4B, IGF1 and IRS1 in a population-based cohort including 3,210 unrelated Chinese Hans from Beijing and Shanghai. We confirmed the associations of GLIS3-rs7034200 with fasting glucose (beta = 0.07 mmol/l, P = 0.03), beta cell function (HOMA-B) (beta = -3.03%, P = 0.009), and type 2 diabetes (OR [95%CI]  = 1.27 [1.09-1.49], P = 0.003) after adjustment for age, sex, region and BMI. The association for type 2 diabetes remained significant after adjusting for other diabetes related risk factors including family history of diabetes, lipid profile, medication information, hypertension and life style factors, while further adjustment for HOMA-B abolished the association. The A-allele of CRY2-rs11605924 was moderately associated with increased risk of combined IFG/type 2 diabetes (OR [95%CI]  = 1.15[1.01-1.30], P = 0.04). SNPs in or near MADD, ADRA2A, PROX1, FADS1, C2CD4B, IGF1, and IRS1 did not exhibit significant associations with type 2 diabetes or related glycemic traits (P≥0.10). In conclusion, our results indicate the associations of GLIS3 locus with type 2 diabetes and impaired fasting glucose in Chinese Hans, partially mediated through impaired beta-cell function. In addition, we also found modest evidence for the association of CRY2-rs11605924 with combined IFG/type 2 diabetes. Show less

To infer the causal association between childhood BMI and age at menarche, we performed a mendelian randomisation analysis using twelve established "BMI-increasing" genetic variants as an instrumental variable (IV) for higher BMI. In 8,156 women of European descent from the EPIC-Norfolk cohort, height was measured at age 39-77 years; age at menarche was self-recalled, as was body weight at age 20 years, and BMI at 20 was calculated as a proxy for childhood BMI. DNA was genotyped for twelve BMI-associated common variants (in/near FTO, MC4R, TMEM18, GNPDA2, KCTD15, NEGR1, BDNF, ETV5, MTCH2, SEC16B, FAIM2 and SH2B1), and for each individual a "BMI-increasing-allele-score" was calculated by summing the number of BMI-increasing alleles across all 12 loci. Using this BMI-increasing-allele-score as an instrumental variable for BMI, each 1 kg/m(2) increase in childhood BMI was predicted to result in a 6.5% (95% CI: 4.6-8.5%) higher absolute risk of early menarche (before age 12 years). While mendelian randomisation analysis is dependent on a number of assumptions, our findings support a causal effect of BMI on early menarche and suggests that increasing prevalence of childhood obesity will lead to similar trends in the prevalence of early menarche. Show less