Search

Discovery of cancer risk variants in the sequence of the germline genome can shed light on carcinogenesis. Here we describe gene burden association analyses, aggregating rare missense and loss of function variants, at 22 cancer sites, including 130,991 cancer cases and 733,486 controls from Iceland, Norway and the United Kingdom. We identified four genes associated with increased cancer risk; the pro-apoptotic BIK for prostate cancer, the autophagy involved ATG12 for colorectal cancer, TG for thyroid cancer and CMTR2 for both lung cancer and cutaneous melanoma. Further, we found genes with rare variants that associate with decreased risk of cancer; AURKB for any cancer, irrespective of site, and PPP1R15A for breast cancer, suggesting that inhibition of PPP1R15A may be a preventive strategy for breast cancer. Our findings pinpoint several new cancer risk genes and emphasize autophagy, apoptosis and cell stress response as a focus point for developing new therapeutics., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Background: Signal Transducer and Activator of Transcription 6 (STAT6) is central to Type 2 (T2) inflammation and common non-coding variants at the STAT6 locus associate with various T2 inflammatory traits, including diseases, and its pathway is widely targeted in asthma treatment., Objective: To test the association of a rare missense variant in STAT6, p.L406P, with T2 inflammatory traits, including the risk of asthma and allergic diseases, and to characterize its functional consequences in cell culture., Methods: We tested association of p.L406P with plasma protein levels, white blood cell counts and the risk of asthma and allergic phenotypes. We tested significant associations in other cohorts using a burden test. The effects of p.L406P on STAT6 protein function were examined in cell lines and by comparing CD4 + T-cell responses from carriers and non-carriers of the variant., Results: p.L406P associated with reduced plasma levels of STAT6 and IgE as well as with lower eosinophil and basophil counts in blood. It also protected against asthma, mostly driven by severe T2 high asthma. We showed that p.L406P led to lower IL-4-induced activation in luciferase reporter assays and lower levels of STAT6 in CD4 + T cells. We identified multiple genes with expression that was affected by the p.L406P genotype upon IL-4 treatment of CD4 + T cells; the effect was consistent with a weaker IL-4 response in carriers than non-carriers of p.L406P., Conclusions: We report a partial loss-of-function variant in STAT6, resulting in dampened IL-4 responses and protection from T2 high asthma, implicating STAT6 as an attractive therapeutic target., (Copyright © 2024. Published by Elsevier Inc.)

Human genetic studies of common variants have provided substantial insight into the biological mechanisms that govern ovarian ageing 1 . Here we report analyses of rare protein-coding variants in 106,973 women from the UK Biobank study, implicating genes with effects around five times larger than previously found for common variants (ETAA1, ZNF518A, PNPLA8, PALB2 and SAMHD1). The SAMHD1 association reinforces the link between ovarian ageing and cancer susceptibility 1 , with damaging germline variants being associated with extended reproductive lifespan and increased all-cause cancer risk in both men and women. Protein-truncating variants in ZNF518A are associated with shorter reproductive lifespan-that is, earlier age at menopause (by 5.61 years) and later age at menarche (by 0.56 years). Finally, using 8,089 sequenced trios from the 100,000 Genomes Project (100kGP), we observe that common genetic variants associated with earlier ovarian ageing associate with an increased rate of maternally derived de novo mutations. Although we were unable to replicate the finding in independent samples from the deCODE study, it is consistent with the expected role of DNA damage response genes in maintaining the genetic integrity of germ cells. This study provides evidence of genetic links between age of menopause and cancer risk., (© 2024. The Author(s).)

Parkinson's disease (PD) is a debilitating neurodegenerative disorder and its rising global incidence highlights the need for the identification of modifiable risk factors. In a gene-based burden test of rare variants (8647 PD cases and 777,693 controls) we discovered a novel association between loss-of-function variants in ITSN1 and PD. This association was further supported with burden data from the Neurodegenerative Disease Knowledge Portal and the Accelerating Medicines Partnership Parkinson's Disease Knowledge Platform. Our findings show that Rho GTPases and disruptions in synaptic vesicle transport may be involved in the pathogenesis of PD, pointing to the possibility of novel therapeutic approaches., (© 2024. The Author(s).)

Autoimmune thyroid disease (AITD) is a common autoimmune disease. In a GWAS meta-analysis of 110,945 cases and 1,084,290 controls, 290 sequence variants at 225 loci are associated with AITD. Of these variants, 115 are previously unreported. Multiomics analysis yields 235 candidate genes outside the MHC-region and the findings highlight the importance of genes involved in T-cell regulation. A rare 5'-UTR variant (rs781745126-T, MAF = 0.13% in Iceland) in LAG3 has the largest effect (OR = 3.42, P = 2.2 × 10 -16 ) and generates a novel start codon for an open reading frame upstream of the canonical protein translation initiation site. rs781745126-T reduces mRNA and surface expression of the inhibitory immune checkpoint LAG-3 co-receptor on activated lymphocyte subsets and halves LAG-3 levels in plasma among heterozygotes. All three homozygous carriers of rs781745126-T have AITD, of whom one also has two other T-cell mediated diseases, that is vitiligo and type 1 diabetes. rs781745126-T associates nominally with vitiligo (OR = 5.1, P = 6.5 × 10 -3 ) but not with type 1 diabetes. Thus, the effect of rs781745126-T is akin to drugs that inhibit LAG-3, which unleash immune responses and can have thyroid dysfunction and vitiligo as adverse events. This illustrates how a multiomics approach can reveal potential drug targets and safety concerns., (© 2024. The Author(s).)

POPDC2 encodes for the Popeye domain-containing protein 2 which has an important role in cardiac pacemaking and conduction, due in part to its cAMP-dependent binding and regulation of TREK-1 potassium channels. Loss of Popdc2 in mice results in sinus pauses and bradycardia and morpholino knockdown of popdc2 in zebrafish results in atrioventricular (AV) block. We identified bi-allelic variants in POPDC2 in 4 families that presented with a phenotypic spectrum consisting of sinus node dysfunction, AV conduction defects and hypertrophic cardiomyopathy. Using homology modelling we show that the identified POPDC2 variants are predicted to diminish the ability of POPDC2 to bind cAMP. In in vitro electrophysiological studies we demonstrated that, while co-expression of wild-type POPDC2 with TREK-1 increased TREK-1 current density, POPDC2 variants found in the patients failed to increase TREK-1 current density. While patient muscle biopsy did not show clear myopathic disease, it showed significant reduction of the expression of both POPDC1 and POPDC2, suggesting that stability and/or membrane trafficking of the POPDC1-POPDC2 complex is impaired by pathogenic variants in any of the two proteins. Single-cell RNA sequencing from human hearts demonstrated that co-expression of POPDC1 and 2 was most prevalent in AV node, AV node pacemaker and AV bundle cells. Sinoatrial node cells expressed POPDC2 abundantly, but expression of POPDC1 was sparse. Together, these results concur with predisposition to AV node disease in humans with loss-of-function variants in POPDC1 and POPDC2 and presence of sinus node disease in POPDC2, but not in POPDC1 related disease in human. Using population-level genetic data of more than 1 million individuals we showed that none of the familial variants were associated with clinical outcomes in heterozygous state, suggesting that heterozygous family members are unlikely to develop clinical manifestations and therefore might not necessitate clinical follow-up. Our findings provide evidence for POPDC2 as the cause of a novel Mendelian autosomal recessive cardiac syndrome, consistent with previous work showing that mice and zebrafish deficient in functional POPDC2 display sinus and AV node dysfunction.

Importance: Recurrent pericarditis is a treatment challenge and often a debilitating condition. Drugs inhibiting interleukin 1 cytokines are a promising new treatment option, but their use is based on scarce biological evidence and clinical trials of modest sizes, and the contributions of innate and adaptive immune processes to the pathophysiology are incompletely understood., Objective: To use human genomics, transcriptomics, and proteomics to shed light on the pathogenesis of pericarditis., Design, Setting, and Participants: This was a meta-analysis of genome-wide association studies of pericarditis from 5 countries. Associations were examined between the pericarditis-associated variants and pericarditis subtypes (including recurrent pericarditis) and secondary phenotypes. To explore mechanisms, associations with messenger RNA expression (cis-eQTL), plasma protein levels (pQTL), and CpG methylation of DNA (ASM-QTL) were assessed. Data from Iceland (deCODE genetics, 1983-2020), Denmark (Copenhagen Hospital Biobank/Danish Blood Donor Study, 1977-2022), the UK (UK Biobank, 1953-2021), the US (Intermountain, 1996-2022), and Finland (FinnGen, 1970-2022) were included. Data were analyzed from September 2022 to August 2023., Exposure: Genotype., Main Outcomes and Measures: Pericarditis., Results: In this genome-wide association study of 4894 individuals with pericarditis (mean [SD] age at diagnosis, 51.4 [17.9] years, 2734 [67.6%] male, excluding the FinnGen cohort), associations were identified with 2 independent common intergenic variants at the interleukin 1 locus on chromosome 2q14. The lead variant was rs12992780 (T) (effect allele frequency [EAF], 31%-40%; odds ratio [OR], 0.83; 95% CI, 0.79-0.87; P = 6.67 × 10-16), downstream of IL1B and the secondary variant rs7575402 (A or T) (EAF, 45%-55%; adjusted OR, 0.89; 95% CI, 0.85-0.93; adjusted P = 9.6 × 10-8). The lead variant rs12992780 had a smaller odds ratio for recurrent pericarditis (0.76) than the acute form (0.86) (P for heterogeneity = .03) and rs7575402 was associated with CpG methylation overlapping binding sites of 4 transcription factors known to regulate interleukin 1 production: PU.1 (encoded by SPI1), STAT1, STAT3, and CCAAT/enhancer-binding protein β (encoded by CEBPB)., Conclusions and Relevance: This study found an association between pericarditis and 2 independent sequence variants at the interleukin 1 gene locus. This finding has the potential to contribute to development of more targeted and personalized therapy of pericarditis with interleukin 1-blocking drugs.

Many sequence variants have additive effects on blood lipid levels and, through that, on the risk of coronary artery disease (CAD). We show that variants also have non-additive effects and interact to affect lipid levels as well as affecting variance and correlations. Variance and correlation effects are often signatures of epistasis or gene-environmental interactions. These complex effects can translate into CAD risk. For example, Trp154Ter in FUT2 protects against CAD among subjects with the A1 blood group, whereas it associates with greater risk of CAD in others. His48Arg in ADH1B interacts with alcohol consumption to affect lipid levels and CAD. The effect of variants in TM6SF2 on blood lipids is greatest among those who never eat oily fish but absent from those who often do. This work demonstrates that variants that affect variance of quantitative traits can allow for the discovery of epistasis and interactions of variants with the environment., Competing Interests: Declaration of interests The authors affiliated with deCODE genetics/Amgen, Inc. are employed by the company., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Importance: Whether protein risk scores derived from a single plasma sample could be useful for risk assessment for atherosclerotic cardiovascular disease (ASCVD), in conjunction with clinical risk factors and polygenic risk scores, is uncertain., Objective: To develop protein risk scores for ASCVD risk prediction and compare them to clinical risk factors and polygenic risk scores in primary and secondary event populations., Design, Setting, and Participants: The primary analysis was a retrospective study of primary events among 13 540 individuals in Iceland (aged 40-75 years) with proteomics data and no history of major ASCVD events at recruitment (study duration, August 23, 2000 until October 26, 2006; follow-up through 2018). We also analyzed a secondary event population from a randomized, double-blind lipid-lowering clinical trial (2013-2016), consisting of individuals with stable ASCVD receiving statin therapy and for whom proteomic data were available for 6791 individuals., Exposures: Protein risk scores (based on 4963 plasma protein levels and developed in a training set in the primary event population); polygenic risk scores for coronary artery disease and stroke; and clinical risk factors that included age, sex, statin use, hypertension treatment, type 2 diabetes, body mass index, and smoking status at the time of plasma sampling., Main Outcomes and Measures: Outcomes were composites of myocardial infarction, stroke, and coronary heart disease death or cardiovascular death. Performance was evaluated using Cox survival models and measures of discrimination and reclassification that accounted for the competing risk of non-ASCVD death., Results: In the primary event population test set (4018 individuals [59.0% women]; 465 events; median follow-up, 15.8 years), the protein risk score had a hazard ratio (HR) of 1.93 per SD (95% CI, 1.75 to 2.13). Addition of protein risk score and polygenic risk scores significantly increased the C index when added to a clinical risk factor model (C index change, 0.022 [95% CI, 0.007 to 0.038]). Addition of the protein risk score alone to a clinical risk factor model also led to a significantly increased C index (difference, 0.014 [95% CI, 0.002 to 0.028]). Among White individuals in the secondary event population (6307 participants; 432 events; median follow-up, 2.2 years), the protein risk score had an HR of 1.62 per SD (95% CI, 1.48 to 1.79) and significantly increased C index when added to a clinical risk factor model (C index change, 0.026 [95% CI, 0.011 to 0.042]). The protein risk score was significantly associated with major adverse cardiovascular events among individuals of African and Asian ancestries in the secondary event population., Conclusions and Relevance: A protein risk score was significantly associated with ASCVD events in primary and secondary event populations. When added to clinical risk factors, the protein risk score and polygenic risk score both provided statistically significant but modest improvement in discrimination.

Urticaria is a skin disorder characterized by outbreaks of raised pruritic wheals. In order to identify sequence variants associated with urticaria, we performed a meta-analysis of genome-wide association studies for urticaria with a total of 40,694 cases and 1,230,001 controls from Iceland, the UK, Finland, and Japan. We also performed transcriptome- and proteome-wide analyses in Iceland and the UK. We found nine sequence variants at nine loci associating with urticaria. The variants are at genes participating in type 2 immune responses and/or mast cell biology (CBLB, FCER1A, GCSAML, STAT6, TPSD1, ZFPM1), the innate immunity (C4), and NF-κB signaling. The most significant association was observed for the splice-donor variant rs56043070[A] (hg38: chr1:247556467) in GCSAML (MAF = 6.6%, OR = 1.24 (95%CI: 1.20-1.28), P-value = 3.6 × 10 -44 ). We assessed the effects of the variants on transcripts, and levels of proteins relevant to urticaria pathophysiology. Our results emphasize the role of type 2 immune response and mast cell activation in the pathogenesis of urticaria. Our findings may point to an IgE-independent urticaria pathway that could help address unmet clinical need., (© 2023. The Author(s).)

Background: Persistent symptoms are common after SARS-CoV-2 infection but correlation with objective measures is unclear., Methods: We invited all 3098 adults who tested SARS-CoV-2 positive in Iceland before October 2020 to the deCODE Health Study. We compared multiple symptoms and physical measures between 1706 Icelanders with confirmed prior infection (cases) who participated, and 619 contemporary and 13,779 historical controls. Cases participated in the study 5-18 months after infection., Results: Here we report that 41 of 88 symptoms are associated with prior infection, most significantly disturbed smell and taste, memory disturbance, and dyspnea. Measured objectively, cases had poorer smell and taste results, less grip strength, and poorer memory recall. Differences in grip strength and memory recall were small. No other objective measure associated with prior infection including heart rate, blood pressure, postural orthostatic tachycardia, oxygen saturation, exercise tolerance, hearing, and traditional inflammatory, cardiac, liver, and kidney blood biomarkers. There was no evidence of more anxiety or depression among cases. We estimate the prevalence of long Covid to be 7% at a median of 8 months after infection., Conclusions: We confirm that diverse symptoms are common months after SARS-CoV-2 infection but find few differences between cases and controls in objective parameters measured. These discrepancies between symptoms and physical measures suggest a more complicated contribution to symptoms related to prior infection than is captured with conventional tests. Traditional clinical assessment is not expected to be particularly informative in relating symptoms to a past SARS-CoV-2 infection., (© 2023. The Author(s).)

The genetic basis of the human vocal system is largely unknown, as are the sequence variants that give rise to individual differences in voice and speech. Here, we couple data on diversity in the sequence of the genome with voice and vowel acoustics in speech recordings from 12,901 Icelanders. We show how voice pitch and vowel acoustics vary across the life span and correlate with anthropometric, physiological, and cognitive traits. We found that voice pitch and vowel acoustics have a heritable component and discovered correlated common variants in ABCC9 that associate with voice pitch. The ABCC9 variants also associate with adrenal gene expression and cardiovascular traits. By showing that voice and vowel acoustics are influenced by genetics, we have taken important steps toward understanding the genetics and evolution of the human vocal system.

Aims: Syncope is a common and clinically challenging condition. In this study, the genetics of syncope were investigated to seek knowledge about its pathophysiology and prognostic implications., Methods and Results: This genome-wide association meta-analysis included 56 071 syncope cases and 890 790 controls from deCODE genetics (Iceland), UK Biobank (United Kingdom), and Copenhagen Hospital Biobank Cardiovascular Study/Danish Blood Donor Study (Denmark), with a follow-up assessment of variants in 22 412 cases and 286 003 controls from Intermountain (Utah, USA) and FinnGen (Finland). The study yielded 18 independent syncope variants, 17 of which were novel. One of the variants, p.Ser140Thr in PTPRN2, affected syncope only when maternally inherited. Another variant associated with a vasovagal reaction during blood donation and five others with heart rate and/or blood pressure regulation, with variable directions of effects. None of the 18 associations could be attributed to cardiovascular or other disorders. Annotation with regard to regulatory elements indicated that the syncope variants were preferentially located in neural-specific regulatory regions. Mendelian randomization analysis supported a causal effect of coronary artery disease on syncope. A polygenic score (PGS) for syncope captured genetic correlation with cardiovascular disorders, diabetes, depression, and shortened lifespan. However, a score based solely on the 18 syncope variants performed similarly to the PGS in detecting syncope risk but did not associate with other disorders., Conclusion: The results demonstrate that syncope has a distinct genetic architecture that implicates neural regulatory processes and a complex relationship with heart rate and blood pressure regulation. A shared genetic background with poor cardiovascular health was observed, supporting the importance of a thorough assessment of individuals presenting with syncope., Competing Interests: Conflicts of interest: The following authors are employees of deCODE genetics/Amgen, Inc.: Hildur M Aegisdottir, Rosa B Thorolfsdottir, Gardar Sveinbjornsson, Olafur A Stefansson, Bjarni Gunnarsson, Vinicius Tragante, Lilja Stefansdottir, Thorgeir E Thorgeirsson, Egil Ferkingstad, Gudmar Thorleifsson, Michael L Frigge, Kristjan E Hjorleifsson, Erna V Ivarsdottir, Anna Helgadottir, Solveig Gretarsdottir, Valgerdur Steinthorsdottir, Asmundur Oddsson, Hannes P Eggertsson, Gisli H Halldorsson, Patrick Sulem, Gudmundur Norddahl, Gudrun Rutsdottir, Gudmundur Thorgeirsson, David O Arnar, Unnur Thorsteinsdottir, Daniel F Gudbjartsson, Hilma Holm, and Kari Stefansson. Karina Banasik and Søren Brunak acknowledge the Novo Nordisk Foundation (grants NNF17OC0027594 and NNF14CC0001)., (© The Author(s) 2023. Published by Oxford University Press on behalf of the European Society of Cardiology. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

The characteristic lobulated nuclear morphology of granulocytes is partially determined by composition of nuclear envelope proteins. Abnormal nuclear morphology is primarily observed as an increased number of hypolobulated immature neutrophils, called band cells, during infection or in rare envelopathies like Pelger-Huët anomaly. To search for sequence variants affecting nuclear morphology of granulocytes, we performed a genome-wide association study using band neutrophil fraction from 88,101 Icelanders. We describe 13 sequence variants affecting band neutrophil fraction at nine loci. Five of the variants are at the Lamin B receptor (LBR) locus, encoding an inner nuclear membrane protein. Mutations in LBR are linked to Pelger-Huët anomaly. In addition, we identify cosegregation of a rare stop-gain sequence variant in LBR and Pelger Huët anomaly in an Icelandic eight generation pedigree, initially reported in 1963. Two of the other loci include genes which, like LBR, play a role in the nuclear membrane function and integrity. These GWAS results highlight the role proteins of the inner nuclear membrane have as important for neutrophil nuclear morphology., (© 2022. The Author(s).)

Predicting the pathogenicity of biallelic missense variants can be challenging. Here, we use a deficit of observed homozygous carriers of missense variants, versus an expected number in a set of 153,054 chip-genotyped Icelanders, to identify potentially pathogenic genotypes. We follow three missense variants with a complete deficit of homozygosity and find that their pathogenic effect in homozygous state ranges from severe childhood disease to early embryonic lethality. One of these variants is in CPSF3, a gene not previously linked to disease. From a set of clinically sequenced Icelanders, and by sequencing archival samples targeted through the Icelandic genealogy, we find four homozygous carriers. Additionally, we find two homozygous carriers of Mexican descent of another missense variant in CPSF3. All six homozygous carriers of missense variants in CPSF3 show severe intellectual disability, seizures, microcephaly, and abnormal muscle tone. Here, we show how the absence of certain homozygous genotypes from a large population set can elucidate causes of previously unexplained recessive diseases and early miscarriage., (© 2022. The Author(s).)

Objective: Biomarkers for diagnosis and progression of osteoarthritis (OA) are lacking. This study was undertaken to identify circulating biomarkers for OA that could predict disease occurrence and/or progression to joint replacement., Methods: Using the SomaScan platform, we measured 4,792 proteins in plasma from 37,278 individuals, of whom 12,178 individuals had OA and 2,524 had undergone joint replacement. We performed a case-control study for identification of potential protein biomarkers for hip, knee, and/or hand OA, and a prospective study for identification of biomarkers for joint replacement., Results: Among the large panel of plasma proteins assessed, cartilage acidic protein 1 (CRTAC1) was the most strongly associated with both OA diagnosis (odds ratio 1.46 [95% confidence interval 1.41-1.52] for knee OA, odds ratio 1.36 [95% confidence interval 1.29-1.43] for hip OA, and odds ratio 1.33 [95% confidence interval 1.26-1.40] for hand OA) and progression to joint replacement (hazard ratio 1.40 [95% confidence interval 1.30-1.51] for knee replacement and hazard ratio 1.31 [95% confidence interval 1.19-1.45] for hip replacement). Patients with OA who were in the highest quintile of risk of joint replacement, based on known risk factors (i.e., age, sex, and body mass index) and plasma CRTAC1 level, were 16 times more likely to undergo knee replacement within 5 years of plasma sample collection than those in the lowest quintile, and 6.5 times more likely to undergo hip replacement. CRTAC1 was not associated with other types of inflammatory arthritis. A specific protein profile was identified in those patients who had undergone joint replacement prior to plasma sample collection., Conclusion: Through a hypothesis-free approach, we identified CRTAC1 in plasma as a novel promising candidate biomarker for OA that is both associated with occurrence of OA and predictive of progression to joint replacement. This biomarker might also be useful in the selection of suitable patients for clinical trial enrollment., (© 2021 The Authors. Arthritis & Rheumatology published by Wiley Periodicals LLC on behalf of American College of Rheumatology.)

Vertigo is the leading symptom of vestibular disorders and a major risk factor for falls. In a genome-wide association study of vertigo (N cases  = 48,072, N controls  = 894,541), we uncovered an association with six common sequence variants in individuals of European ancestry, including missense variants in ZNF91, OTOG, OTOGL, and TECTA, and a cis-eQTL for ARMC9. The association of variants in ZNF91, OTOGL, and OTOP1 was driven by an association with benign paroxysmal positional vertigo. Using previous reports of sequence variants associating with age-related hearing impairment and motion sickness, we found eight additional variants that associate with vertigo. Although disorders of the auditory and the vestibular system may co-occur, none of the six genome-wide significant vertigo variants were associated with hearing loss and only one was associated with age-related hearing impairment. Our results uncovered sequence variants associating with vertigo in a genome-wide association study and implicated genes with known roles in inner ear development, maintenance, and disease., (© 2021. The Author(s).)

Birth weight is a common measure of fetal growth that is associated with a range of health outcomes. It is directly affected by the fetal genome and indirectly by the maternal genome. We performed genome-wide association studies on birth weight in the genomes of the child and parents and further analyzed birth length and ponderal index, yielding a total of 243 fetal growth variants. We clustered those variants based on the effects of transmitted and nontransmitted alleles on birth weight. Out of 141 clustered variants, 22 were consistent with parent-of-origin-specific effects. We further used haplotype-specific polygenic risk scores to directly test the relationship between adult traits and birth weight. Our results indicate that the maternal genome contributes to increased birth weight through blood-glucose-raising alleles while blood-pressure-raising alleles reduce birth weight largely through the fetal genome., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Predicting all-cause mortality risk is challenging and requires extensive medical data. Recently, large-scale proteomics datasets have proven useful for predicting health-related outcomes. Here, we use measurements of levels of 4,684 plasma proteins in 22,913 Icelanders to develop all-cause mortality predictors both for short- and long-term risk. The participants were 18-101 years old with a mean follow up of 13.7 (sd. 4.7) years. During the study period, 7,061 participants died. Our proposed predictor outperformed, in survival prediction, a predictor based on conventional mortality risk factors. We could identify the 5% at highest risk in a group of 60-80 years old, where 88% died within ten years and 5% at the lowest risk where only 1% died. Furthermore, the predicted risk of death correlates with measures of frailty in an independent dataset. Our results show that the plasma proteome can be used to assess general health and estimate the risk of death.

Age-related hearing impairment (ARHI) is the most common sensory disorder in older adults. We conducted a genome-wide association meta-analysis of 121,934 ARHI cases and 591,699 controls from Iceland and the UK. We identified 21 novel sequence variants, of which 13 are rare, under either additive or recessive models. Of special interest are a missense variant in LOXHD1 (MAF = 1.96%) and a tandem duplication in FBF1 covering 4 exons (MAF = 0.22%) associating with ARHI (OR = 3.7 for homozygotes, P = 1.7 × 10 -22 and OR = 4.2 for heterozygotes, P = 5.7 × 10 -27 , respectively). We constructed an ARHI genetic risk score (GRS) using common variants and showed that a common variant GRS can identify individuals at risk comparable to carriers of rare high penetrance variants. Furthermore, we found that ARHI and tinnitus share genetic causes. This study sheds a new light on the genetic architecture of ARHI, through several rare variants in both Mendelian deafness genes and genes not previously linked to hearing.

Aims: The aim of this study was to use human genetics to investigate the pathogenesis of sick sinus syndrome (SSS) and the role of risk factors in its development., Methods and Results: We performed a genome-wide association study of 6469 SSS cases and 1 000 187 controls from deCODE genetics, the Copenhagen Hospital Biobank, UK Biobank, and the HUNT study. Variants at six loci associated with SSS, a reported missense variant in MYH6, known atrial fibrillation (AF)/electrocardiogram variants at PITX2, ZFHX3, TTN/CCDC141, and SCN10A and a low-frequency (MAF = 1.1-1.8%) missense variant, p.Gly62Cys in KRT8 encoding the intermediate filament protein keratin 8. A full genotypic model best described the p.Gly62Cys association (P = 1.6 × 10-20), with an odds ratio (OR) of 1.44 for heterozygotes and a disproportionally large OR of 13.99 for homozygotes. All the SSS variants increased the risk of pacemaker implantation. Their association with AF varied and p.Gly62Cys was the only variant not associating with any other arrhythmia or cardiovascular disease. We tested 17 exposure phenotypes in polygenic score (PGS) and Mendelian randomization analyses. Only two associated with the risk of SSS in Mendelian randomization, AF, and lower heart rate, suggesting causality. Powerful PGS analyses provided convincing evidence against causal associations for body mass index, cholesterol, triglycerides, and type 2 diabetes (P > 0.05)., Conclusion: We report the associations of variants at six loci with SSS, including a missense variant in KRT8 that confers high risk in homozygotes and points to a mechanism specific to SSS development. Mendelian randomization supports a causal role for AF in the development of SSS., (© The Author(s) 2021. Published by Oxford University Press on behalf of the European Society of Cardiology.)

Olfactory receptor (OR) genes in humans form a special class characterized by unusually high DNA sequence diversity, which should give rise to differences in perception and behavior. In the largest genome-wide association study to date based on olfactory testing, we investigated odor perception and naming with smell tasks performed by 9,122 Icelanders, with replication in a separate sample of 2,204 individuals. We discovered an association between a low-frequency missense variant in TAAR5 and reduced intensity rating of fish odor containing trimethylamine (p.Ser95Pro, p combined  = 5.6 × 10 -15 ). We demonstrate that TAAR5 genotype affects aversion to fish odor, reflected by linguistic descriptions of the odor and pleasantness ratings. We also discovered common sequence variants in two canonical olfactory receptor loci that associate with increased intensity and naming of licorice odor (trans-anethole: lead variant p.Lys233Asn in OR6C70, p combined  = 8.8 × 10 -16 and p combined  = 1.4 × 10 -9 ) and enhanced naming of cinnamon (trans-cinnamaldehyde; intergenic variant rs317787-T, p combined  = 5.0 × 10 -17 ). Together, our results show that TAAR5 genotype variation influences human odor responses and highlight that sequence diversity in canonical OR genes can lead to enhanced olfactory ability, in contrast to the view that greater tolerance for mutations in the human OR repertoire leads to diminished function., Competing Interests: Declaration of Interests R.S.G., E.V.I., A.H., N.K.H., G.R., G.A.A., A.S., B.A.J., G.L.N., M.O.U., H. Helgason, B.V.H., M.S.N., V.T., G.S., T.T., A.O., R.P.K., G.B., G. T., I. J., H. Holm, D.F.G., U.T., H.S., P.S., and K.S. are employees of deCODE genetics/Amgen, Inc., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: Little is known about the nature and durability of the humoral immune response to infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)., Methods: We measured antibodies in serum samples from 30,576 persons in Iceland, using six assays (including two pan-immunoglobulin [pan-Ig] assays), and we determined that the appropriate measure of seropositivity was a positive result with both pan-Ig assays. We tested 2102 samples collected from 1237 persons up to 4 months after diagnosis by a quantitative polymerase-chain-reaction (qPCR) assay. We measured antibodies in 4222 quarantined persons who had been exposed to SARS-CoV-2 and in 23,452 persons not known to have been exposed., Results: Of the 1797 persons who had recovered from SARS-CoV-2 infection, 1107 of the 1215 who were tested (91.1%) were seropositive; antiviral antibody titers assayed by two pan-Ig assays increased during 2 months after diagnosis by qPCR and remained on a plateau for the remainder of the study. Of quarantined persons, 2.3% were seropositive; of those with unknown exposure, 0.3% were positive. We estimate that 0.9% of Icelanders were infected with SARS-CoV-2 and that the infection was fatal in 0.3%. We also estimate that 56% of all SARS-CoV-2 infections in Iceland had been diagnosed with qPCR, 14% had occurred in quarantined persons who had not been tested with qPCR (or who had not received a positive result, if tested), and 30% had occurred in persons outside quarantine and not tested with qPCR., Conclusions: Our results indicate that antiviral antibodies against SARS-CoV-2 did not decline within 4 months after diagnosis. We estimate that the risk of death from infection was 0.3% and that 44% of persons infected with SARS-CoV-2 in Iceland were not diagnosed by qPCR., (Copyright © 2020 Massachusetts Medical Society.)

Autoimmune thyroid disease is the most common autoimmune disease and is highly heritable 1 . Here, by using a genome-wide association study of 30,234 cases and 725,172 controls from Iceland and the UK Biobank, we find 99 sequence variants at 93 loci, of which 84 variants are previously unreported 2-7 . A low-frequency (1.36%) intronic variant in FLT3 (rs76428106-C) has the largest effect on risk of autoimmune thyroid disease (odds ratio (OR) = 1.46, P = 2.37 × 10 -24 ). rs76428106-C is also associated with systemic lupus erythematosus (OR = 1.90, P = 6.46 × 10 -4 ), rheumatoid factor and/or anti-CCP-positive rheumatoid arthritis (OR = 1.41, P = 4.31 × 10 -4 ) and coeliac disease (OR = 1.62, P = 1.20 × 10 -4 ). FLT3 encodes fms-related tyrosine kinase 3, a receptor that regulates haematopoietic progenitor and dendritic cells. RNA sequencing revealed that rs76428106-C generates a cryptic splice site, which introduces a stop codon in 30% of transcripts that are predicted to encode a truncated protein, which lacks its tyrosine kinase domains. Each copy of rs76428106-C doubles the plasma levels of the FTL3 ligand. Activating somatic mutations in FLT3 are associated with acute myeloid leukaemia 8 with a poor prognosis and rs76428106-C also predisposes individuals to acute myeloid leukaemia (OR = 1.90, P = 5.40 × 10 -3 ). Thus, a predicted loss-of-function germline mutation in FLT3 causes a reduction in full-length FLT3, with a compensatory increase in the levels of its ligand and an increased disease risk, similar to that of a gain-of-function mutation.

Hemoglobin is the essential oxygen-carrying molecule in humans and is regulated by cellular iron and oxygen sensing mechanisms. To search for novel variants associated with hemoglobin concentration, we performed genome-wide association studies of hemoglobin concentration using a combined set of 684,122 individuals from Iceland and the UK. Notably, we found seven novel variants, six rare coding and one common, at the ACO1 locus associating with either decreased or increased hemoglobin concentration. Of these variants, the missense Cys506Ser and the stop-gained Lys334Ter mutations are specific to eight and ten generation pedigrees, respectively, and have the two largest effects in the study (Effect Cys506Ser  = -1.61 SD, CI 95  = [-1.98, -1.35]; Effect Lys334Ter  = 0.63 SD, CI 95  = [0.36, 0.91]). We also find Cys506Ser to associate with increased risk of persistent anemia (OR = 17.1, P = 2 × 10 -14 ). The strong bidirectional effects seen in this study implicate ACO1, a known iron sensing molecule, as a major homeostatic regulator of hemoglobin concentration.

Asthma is one of the most common chronic diseases affecting both children and adults. We report a genome-wide association meta-analysis of 69,189 cases and 702,199 controls from Iceland and UK biobank. We find 88 asthma risk variants at 56 loci, 19 previously unreported, and evaluate their effect on other asthma and allergic phenotypes. Of special interest are two low frequency variants associated with protection against asthma; a missense variant in TNFRSF8 and 3' UTR variant in TGFBR1. Functional studies show that the TNFRSF8 variant reduces TNFRSF8 expression both on cell surface and in soluble form, acting as loss of function. eQTL analysis suggests that the TGFBR1 variant acts through gain of function and together with an intronic variant in a downstream gene, SMAD3, points to defective TGFβR1 signaling as one of the biological perturbations increasing asthma risk. Our results increase the number of asthma variants and implicate genes with known role in T cell regulation, inflammation and airway remodeling in asthma pathogenesis.

Bone area is one measure of bone size that is easily derived from dual-energy X-ray absorptiometry (DXA) scans. In a GWA study of DXA bone area of the hip and lumbar spine (N ≥ 28,954), we find thirteen independent association signals at twelve loci that replicate in samples of European and East Asian descent (N = 13,608 - 21,277). Eight DXA area loci associate with osteoarthritis, including rs143384 in GDF5 and a missense variant in COL11A1 (rs3753841). The strongest DXA area association is with rs11614913[T] in the microRNA MIR196A2 gene that associates with lumbar spine area (P = 2.3 × 10 -42 , β = -0.090) and confers risk of hip fracture (P = 1.0 × 10 -8 , OR = 1.11). We demonstrate that the risk allele is less efficient in repressing miR-196a-5p target genes. We also show that the DXA area measure contributes to the risk of hip fracture independent of bone density.

Nerve conduction (NC) studies generate measures of peripheral nerve function that can reveal underlying pathology due to axonal loss, demyelination or both. We perform a genome-wide association study of sural NC amplitude and velocity in 7045 Icelanders and find a low-frequency splice-donor variant in PRPH (c.996+1G>A; MAF = 1.32%) associating with decreased NC amplitude but not velocity. PRPH encodes peripherin, an intermediate filament (IF) protein involved in cytoskeletal development and maintenance of neurons. Through RNA and protein studies, we show that the variant leads to loss-of-function (LoF), as when over-expressed in a cell line devoid of other IFs, it does not allow formation of the normal filamentous structure of peripherin, yielding instead punctate protein inclusions. Recall of carriers for neurological assessment confirms that from an early age, homozygotes have significantly lower sural NC amplitude than non-carriers and are at risk of a mild, early-onset, sensory-negative, axonal polyneuropathy.

Urine dipstick tests are widely used in routine medical care to diagnose kidney and urinary tract and metabolic diseases. Several environmental factors are known to affect the test results, whereas the effects of genetic diversity are largely unknown. We tested 32.5 million sequence variants for association with urinary biomarkers in a set of 150 274 Icelanders with urine dipstick measurements. We detected 20 association signals, of which 14 are novel, associating with at least one of five clinical entities defined by the urine dipstick: glucosuria, ketonuria, proteinuria, hematuria and urine pH. These include three independent glucosuria variants at SLC5A2, the gene encoding the sodium-dependent glucose transporter (SGLT2), a protein targeted pharmacologically to increase urinary glucose excretion in the treatment of diabetes. Two variants associating with proteinuria are in LRP2 and CUBN, encoding the co-transporters megalin and cubilin, respectively, that mediate proximal tubule protein uptake. One of the hematuria-associated variants is a rare, previously unreported 2.5 kb exonic deletion in COL4A3. Of the four signals associated with urine pH, we note that the pH-increasing alleles of two variants (POU2AF1, WDR72) associate significantly with increased risk of kidney stones. Our results reveal that genetic factors affect variability in urinary biomarkers, in both a disease dependent and independent context., (© The Author(s) 2018. Published by Oxford University Press.)

The corneal endothelium is vital for transparency and proper hydration of the cornea. Here, we conduct a genome-wide association study of corneal endothelial cell density (cells/mm 2 ), coefficient of cell size variation (CV), percentage of hexagonal cells (HEX) and central corneal thickness (CCT) in 6,125 Icelanders and find associations at 10 loci, including 7 novel. We assess the effects of these variants on various ocular biomechanics such as corneal hysteresis (CH), as well as eye diseases such as glaucoma and corneal dystrophies. Most notably, an intergenic variant close to ANAPC1 (rs78658973[A], frequency = 28.3%) strongly associates with decreased cell density and accounts for 24% of the population variance in cell density (β = -0.77 SD, P = 1.8 × 10 -314 ) and associates with increased CH (β = 0.19 SD, P = 2.6 × 10 -19 ) without affecting risk of corneal diseases and glaucoma. Our findings indicate that despite correlations between cell density and eye diseases, low cell density does not increase the risk of disease.

Mutations in genes encoding subunits of the phagocyte NADPH oxidase complex are recognized to cause chronic granulomatous disease (CGD), a severe primary immunodeficiency. Here we describe how deficiency of CYBC1, a previously uncharacterized protein in humans (C17orf62), leads to reduced expression of NADPH oxidase's main subunit (gp91 phox ) and results in CGD. Analyzing two brothers diagnosed with CGD we identify a homozygous loss-of-function mutation, p.Tyr2Ter, in CYBC1. Imputation of p.Tyr2Ter into 155K chip-genotyped Icelanders reveals six additional homozygotes, all with signs of CGD, manifesting as colitis, rare infections, or a severely impaired PMA-induced neutrophil oxidative burst. Homozygosity for p.Tyr2Ter consequently associates with inflammatory bowel disease (IBD) in Iceland (P = 8.3 × 10 -8 ; OR = 67.6), as well as reduced height (P = 3.3 × 10 -4 ; -8.5 cm). Overall, we find that CYBC1 deficiency results in CGD characterized by colitis and a distinct profile of infections indicative of macrophage dysfunction.

Most sequence variants identified hitherto in genome-wide association studies (GWAS) of atrial fibrillation are common, non-coding variants associated with risk through unknown mechanisms. We performed a meta-analysis of GWAS of atrial fibrillation among 29,502 cases and 767,760 controls from Iceland and the UK Biobank with follow-up in samples from Norway and the US, focusing on low-frequency coding and splice variants aiming to identify causal genes. We observe associations with one missense (OR = 1.20) and one splice-donor variant (OR = 1.50) in RPL3L , the first ribosomal gene implicated in atrial fibrillation to our knowledge. Analysis of 167 RNA samples from the right atrium reveals that the splice-donor variant in RPL3L results in exon skipping. We also observe an association with a missense variant in MYZAP (OR = 1.38), encoding a component of the intercalated discs of cardiomyocytes. Both discoveries emphasize the close relationship between the mechanical and electrical function of the heart., Competing Interests: The following authors affiliated with deCODE genetics/Amgen, Inc. are employed by the company: R.B.T., G.S., P.S., S.J., G.H.H., P.M., E.V.I., O.B.D., R.P.K., G.T., A.H., S.G., G.N., S.R., V.T., B.V.H., I.J., D.O.A., U.T., D.F.G., H.H., and K.S. T.R.P. has received consulting and/or speakers honoraria from Amgen, Sanofi, and Merck (all minor). M.S.S. has received research grant support through Brigham and Women’s Hospital from Amgen, AstraZeneca, Daiichi-Sankyo, Eisai, GlaxoSmithKline, Intarcia, Janssen Research and Development, MedImmune, Merck, Novartis, Pfizer, Poxel, and Takeda and consulting honoraria from Amgen, CVS Caremark, Esperion, Intarcia, Ionis, Janssen Research and Development, MedImmune, and Merck. B.T., A.S.V., J.B.N., J.T.S., O.L.H., F.W.A., D.M.R., D.D., C.J.W., M.-L.L., and K.H. have no relationship with industry to disclose.

The cytokine erythropoietin (EPO), signalling through the EPO receptor (EPO-R), is essential for the formation of red blood cells. We performed a genome-wide association study (GWAS) testing 32.5 million sequence variants for association with serum EPO levels in a set of 4187 individuals. We detect an association between a rare and well imputed stop-gained variant rs370865377[A] (p.Gln82Ter) in EPOR , carried by 1 in 550 Icelanders, and increased serum EPO levels (MAF = 0.09%, Effect = 1.47 SD, P  = 3.3 × 10 -7 ). We validated these findings by measuring serum EPO levels in 34 additional pairs of carriers and matched controls and found carriers to have 3.23-fold higher EPO levels than controls ( P  = 1.7 × 10 -6 ; P combined  = 1.6 × 10 -11 ). In contrast to previously reported EPOR mutations, p.Gln82Ter does not associate with haemoglobin levels (Effect = -0.045 SD, P  = 0.32, N  = 273,160), probably due to a compensatory EPO upregulation in response to EPO-R hypo-responsiveness., Competing Interests: Authors affiliated with deCODE genetics/Amgen declare competing financial interests as employees. The remaining authors declare no competing financial interests. All authors declare no non-financial competing interests.

Sequence variants that affect mean fasting glucose levels do not necessarily affect risk for type 2 diabetes (T2D). We assessed the effects of 36 reported glucose-associated sequence variants on between- and within-subject variance in fasting glucose levels in 69,142 Icelanders. The variant in TCF7L2 that increases fasting glucose levels increases between-subject variance (5.7% per allele, P = 4.2 × 10 -10 ), whereas variants in GCK and G6PC2 that increase fasting glucose levels decrease between-subject variance (7.5% per allele, P = 4.9 × 10 -11 and 7.3% per allele, P = 7.5 × 10 -18 , respectively). Variants that increase mean and between-subject variance in fasting glucose levels tend to increase T2D risk, whereas those that increase the mean but reduce variance do not (r 2 = 0.61). The variants that increase between-subject variance increase fasting glucose heritability estimates. Intuitively, our results show that increasing the mean and variance of glucose levels is more likely to cause pathologically high glucose levels than increase in the mean offset by a decrease in variance.