👤 Hreinn Stefansson

Obesity is associated with adverse effects on health and quality of life. Improved understanding of its underlying pathophysiology is essential for developing counteractive measures. To search for sequence variants with large effects on BMI, we perform a multi-ancestry meta-analysis of 13 genome-wide association studies on BMI, including data derived from 1,534,555 individuals of European ancestry, 339,657 of Asian ancestry, and 130,968 of African ancestry. We identify an intergenic 262,760 base pair deletion at the MC4R locus that associates with 4.11 kg/m 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

It is not clear if antagonizing the GIP (glucose-dependent insulinotropic polypeptide) receptor (GIPR) for treatment of obesity is likely to increase the risk of fractures, or to lower bone mineral density (BMD) beyond what is expected with rapid weight loss. The objective of this study was to investigate the risk of fracture and BMD of sequence variants in GIPR that reduce the activity of the GIP receptor and have been associated with reduced body mass index (BMI). We analyzed the association of 3 missense variants in GIPR, a common variant, rs1800437 (p.Glu354Gln), and 2 rare variants, rs139215588 (p.Arg190Gln) and rs143430880 (p.Glu288Gly), as well as a burden of predicted loss-of-function (LoF) variants with risk of fracture and with BMD in a large meta-analysis of up to 1.2 million participants. We analyzed associations with fractures at different skeletal sites in the general population: any fractures, hip fractures, vertebral fractures and forearm fractures, and specifically nonvertebral and osteoporotic fractures in postmenopausal women. We also evaluated associations with BMD at the lumbar spine, femoral neck, and total body measured with dual-energy x-ray absorptiometry (DXA), and with BMD estimated from heel ultrasound (eBMD). None of the 3 missense variants in GIPR was significantly associated with increased risk of fractures or with lower BMD. Burden of LoF variants in GIPR was not associated with fractures or with BMD measured with clinically validated DXA, but was associated with eBMD. Missense variants in GIPR, or burden of LoF variants in the gene, are not associated with risk of fractures or with lower BMD. Show less

Endometriosis is a common condition associated with debilitating pelvic pain and infertility. A genome-wide association study meta-analysis, including 60,674 cases and 701,926 controls of European and East Asian descent, identified 42 genome-wide significant loci comprising 49 distinct association signals. Effect sizes were largest for stage 3/4 disease, driven by ovarian endometriosis. Identified signals explained up to 5.01% of disease variance and regulated expression or methylation of genes in endometrium and blood, many of which were associated with pain perception/maintenance (SRP14/BMF, GDAP1, MLLT10, BSN and NGF). We observed significant genetic correlations between endometriosis and 11 pain conditions, including migraine, back and multisite chronic pain (MCP), as well as inflammatory conditions, including asthma and osteoarthritis. Multitrait genetic analyses identified substantial sharing of variants associated with endometriosis and MCP/migraine. Targeted investigations of genetically regulated mechanisms shared between endometriosis and other pain conditions are needed to aid the development of new treatments and facilitate early symptomatic intervention. Show less

Migraine affects over a billion individuals worldwide but its genetic underpinning remains largely unknown. Here, we performed a genome-wide association study of 102,084 migraine cases and 771,257 controls and identified 123 loci, of which 86 are previously unknown. These loci provide an opportunity to evaluate shared and distinct genetic components in the two main migraine subtypes: migraine with aura and migraine without aura. Stratification of the risk loci using 29,679 cases with subtype information indicated three risk variants that seem specific for migraine with aura (in HMOX2, CACNA1A and MPPED2), two that seem specific for migraine without aura (near SPINK2 and near FECH) and nine that increase susceptibility for migraine regardless of subtype. The new risk loci include genes encoding recent migraine-specific drug targets, namely calcitonin gene-related peptide (CALCA/CALCB) and serotonin 1F receptor (HTR1F). Overall, genomic annotations among migraine-associated variants were enriched in both vascular and central nervous system tissue/cell types, supporting unequivocally that neurovascular mechanisms underlie migraine pathophysiology. Show less

The HERMES (HEart failure Molecular Epidemiology for Therapeutic targetS) consortium aims to identify the genomic and molecular basis of heart failure. The consortium currently includes 51 studies from 11 countries, including 68 157 heart failure cases and 949 888 controls, with data on heart failure events and prognosis. All studies collected biological samples and performed genome-wide genotyping of common genetic variants. The enrolment of subjects into participating studies ranged from 1948 to the present day, and the median follow-up following heart failure diagnosis ranged from 2 to 116 months. Forty-nine of 51 individual studies enrolled participants of both sexes; in these studies, participants with heart failure were predominantly male (34-90%). The mean age at diagnosis or ascertainment across all studies ranged from 54 to 84 years. Based on the aggregate sample, we estimated 80% power to genetic variant associations with risk of heart failure with an odds ratio of ≥1.10 for common variants (allele frequency ≥ 0.05) and ≥1.20 for low-frequency variants (allele frequency 0.01-0.05) at P < 5 × 10 HERMES is a global collaboration aiming to (i) identify the genetic determinants of heart failure; (ii) generate insights into the causal pathways leading to heart failure and enable genetic approaches to target prioritization; and (iii) develop genomic tools for disease stratification and risk prediction. 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

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

Bone mineral density (BMD) is a measure of osteoporosis and is useful in evaluating the risk of fracture. In a genome-wide association study of BMD among 20,100 Icelanders, with follow-up in 10,091 subjects of European and East-Asian descent, we found a new BMD locus that harbours the PTCH1 gene, represented by rs28377268 (freq. 11.4-22.6%) that associates with reduced spine BMD (P=1.0 × 10(-11), β=-0.09). We also identified a new spine BMD signal in RSPO3, rs577721086 (freq. 6.8%), that associates with increased spine BMD (P=6.6 × 10(-10), β=0.14). Importantly, both variants associate with osteoporotic fractures and affect expression of the PTCH1 and RSPO3 genes that is in line with their influence on BMD and known biological function of these genes. Additional new BMD signals were also found at the AXIN1 and SOST loci and a new lead SNP at the EN1 locus. Show less

Sequence variants affecting blood lipids and coronary artery disease (CAD) may enhance understanding of the atherogenicity of lipid fractions. Using a large resource of whole-genome sequence data, we examined rare and low-frequency variants for association with non-HDL cholesterol, HDL cholesterol, LDL cholesterol, and triglycerides in up to 119,146 Icelanders. We discovered 13 variants with large effects (within ANGPTL3, APOB, ABCA1, NR1H3, APOA1, LIPC, CETP, LDLR, and APOC1) and replicated 14 variants. Five variants within PCSK9, APOA1, ANGPTL4, and LDLR associate with CAD (33,090 cases and 236,254 controls). We used genetic risk scores for the lipid fractions to examine their causal relationship with CAD. The non-HDL cholesterol genetic risk score associates most strongly with CAD (P = 2.7 × 10(-28)), and no other genetic risk score associates with CAD after accounting for non-HDL cholesterol. The genetic risk score for non-HDL cholesterol confers CAD risk beyond that of LDL cholesterol (P = 5.5 × 10(-8)), suggesting that targeting atherogenic remnant cholesterol may reduce cardiovascular risk. Show less

The ages of puberty, first sexual intercourse and first birth signify the onset of reproductive ability, behavior and success, respectively. In a genome-wide association study of 125,667 UK Biobank participants, we identify 38 loci associated (P < 5 × 10(-8)) with age at first sexual intercourse. These findings were taken forward in 241,910 men and women from Iceland and 20,187 women from the Women's Genome Health Study. Several of the identified loci also exhibit associations (P < 5 × 10(-8)) with other reproductive and behavioral traits, including age at first birth (variants in or near ESR1 and RBM6-SEMA3F), number of children (CADM2 and ESR1), irritable temperament (MSRA) and risk-taking propensity (CADM2). Mendelian randomization analyses infer causal influences of earlier puberty timing on earlier first sexual intercourse, earlier first birth and lower educational attainment. In turn, likely causal consequences of earlier first sexual intercourse include reproductive, educational, psychiatric and cardiometabolic outcomes. Show less

The geographic isolation and homogeneous population of Iceland are ideally suited to ascertain clinical and genetic characteristics of hypertrophic cardiomyopathy (HCM) at the population level. Medical records and cardiac imaging studies obtained between 1997 and 2010 were reviewed to identify Icelandic patients with HCM. Surviving patients were recruited for clinical and genetic studies. A previously identified Icelandic mutation, MYBPC3 c.927-2A>G, was genotyped, and mutation-negative samples were sequenced for HCM genes and other hypertrophic genes. Record review identified 180 patients with HCM. Genetic analyses of 151 patients defined pathogenic mutations in 101 (67%), including MYBPC3 c.927-2A>G (88 patients, 58%), 4 other MYBPC3 or MYH7 mutations (5 patients, 3.3%), and 2 GLA mutations (8 patients, 5.3%). Haplotype and genetic genealogical data defined MYBPC3 c.927-2A>G as a founder mutation, introduced into the Icelandic population in the 15th century, with a current population prevalence of 0.36%. MYBPC3 c.927-2A>G mutation carriers exhibited phenotypic diversity but were younger at diagnosis (42 versus 49 years; P=0.001) and sustained more adverse events (15% versus 2%; P=0.02) than mutation-negative patients. All-cause mortality for patients with HCM was similar to that of an age-matched Icelandic population (hazard ratio, 0.98; P=0.9). HCM-related mortality (0.78%/y) occurred at a mean age of 68 compared with 81 years for non-HCM-related mortality (P=0.02). A founder MYBPC3 mutation that arose >550 years ago is the predominant cause of HCM in Iceland. The MYBPC3 c.927-2A>G mutation is associated with low adverse event rates but earlier cardiovascular mortality, illustrating the impact of genotype on outcomes in HCM. 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

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

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

We identified a marker in LINGO1 showing genome-wide significant association (P = 1.2 x 10(-9), odds ratio = 1.55) with essential tremor. LINGO1 has potent, negative regulatory influences on neuronal survival and is also important in regulating both central-nervous-system axon regeneration and oligodendrocyte maturation. Increased axon integrity observed in Lingo1 mouse [corrected] knockout models highlights the potential role of LINGO1 in the pathophysiology of ET [corrected] Show less

Deletions within the neurexin 1 gene (NRXN1; 2p16.3) are associated with autism and have also been reported in two families with schizophrenia. We examined NRXN1, and the closely related NRXN2 and NRXN3 genes, for copy number variants (CNVs) in 2977 schizophrenia patients and 33 746 controls from seven European populations (Iceland, Finland, Norway, Germany, The Netherlands, Italy and UK) using microarray data. We found 66 deletions and 5 duplications in NRXN1, including a de novo deletion: 12 deletions and 2 duplications occurred in schizophrenia cases (0.47%) compared to 49 and 3 (0.15%) in controls. There was no common breakpoint and the CNVs varied from 18 to 420 kb. No CNVs were found in NRXN2 or NRXN3. We performed a Cochran-Mantel-Haenszel exact test to estimate association between all CNVs and schizophrenia (P = 0.13; OR = 1.73; 95% CI 0.81-3.50). Because the penetrance of NRXN1 CNVs may vary according to the level of functional impact on the gene, we next restricted the association analysis to CNVs that disrupt exons (0.24% of cases and 0.015% of controls). These were significantly associated with a high odds ratio (P = 0.0027; OR 8.97, 95% CI 1.8-51.9). We conclude that NRXN1 deletions affecting exons confer risk of schizophrenia. Show less

Protein phosphatase 6 (PP6) is an essential Ser/Thr phosphatase conserved among eukaryotes. The Saccharomyces cerevisiae homologue of PP6 called Sit4 depends on association with SAPS domain subunits. This study used a human SAPS domain subunit FLAG-PP6R1 to identify endogenous interacting proteins. Mass spectrometry identified coprecipitating proteins as PP6 catalytic subunit and three ankyrin repeat proteins (Ankrd28, Ankrd44, and Ankrd52). These proteins have extensive sequence identity to one another but segregate into separate branches on a phylogenetic tree for vertebrate species, suggesting individual biological functions. Tagged Ankrd28 coprecipitated with PP6, not with PP2A or PP4, and with SAPS domain subunits PP6R1 and PP6R3. Tagged PP6 coprecipitated endogenous SAPS domain subunits and Ankrd28. The C-terminal region of PP6R1 was sufficient to coprecipitate Ankrd28, but not PP6, demonstrating that PP6R1 acts as a scaffold with separate regions for binding to PP6 and to Ankrd28. Endogenous PP6 holoenzymes with PP6R1 and PP6R3 subunits were resolved by DEAE chromatography and eluted together with Ankrd28 at Mr > 440 kDa from Superose 12. Knockdown of PP6R1 or Ankrd28, but not PP6R3, produced equivalent enhancement of IkappaBepsilon degradation in response to TNFalpha. The results suggest that PP6 functions as a heterotrimer, composed of the PP6 catalytic subunit bound to a SAPS domain scaffold subunit that associates with Ankrd28. We propose that the SAPS and ankyrin repeat regulatory subunits determine the function and specificity of PP6. Show less