Search

Background: Bone mineral density influences the risk of osteoporosis later in life and is useful in the evaluation of the risk of fracture. We aimed to identify sequence variants associated with bone mineral density and fracture.Methods: We performed a quantitative trait analysis of data from 5861 Icelandic subjects (the discovery set), testing for an association between 301,019 single-nucleotide polymorphisms (SNPs) and bone mineral density of the hip and lumbar spine. We then tested for an association between 74 SNPs (most of which were implicated in the discovery set) at 32 loci in replication sets of Icelandic, Danish, and Australian subjects (4165, 2269, and 1491 subjects, respectively).Results: Sequence variants in five genomic regions were significantly associated with bone mineral density in the discovery set and were confirmed in the replication sets (combined P values, 1.2x10(-7) to 2.0x10(-21)). Three regions are close to or within genes previously shown to be important to the biologic characteristics of bone: the receptor activator of nuclear factor-kappaB ligand gene (RANKL) (chromosomal location, 13q14), the osteoprotegerin gene (OPG) (8q24), and the estrogen receptor 1 gene (ESR1) (6q25). The two other regions are close to the zinc finger and BTB domain containing 40 gene (ZBTB40) (1p36) and the major histocompatibility complex region (6p21). The 1p36, 8q24, and 6p21 loci were also associated with osteoporotic fractures, as were loci at 18q21, close to the receptor activator of the nuclear factor-kappaB gene (RANK), and loci at 2p16 and 11p11.Conclusions: We have discovered common sequence variants that are consistently associated with bone mineral density and with low-trauma fractures in three populations of European descent. Although these variants alone are not clinically useful in the prediction of risk to the individual person, they provide insight into the biochemical pathways underlying osteoporosis. [ABSTRACT FROM AUTHOR]

Background: The contributions of genetic and environmental risk factors to hidradenitis suppurativa (HS) are both poorly understood., Objective: To identify sequence variants that associate with HS and determine the contribution of environmental risk factors and inflammatory diseases to HS pathogenesis., Methods: A genome-wide association meta-analysis of 4814 HS cases (Denmark: 1977; Iceland: 1266; Finland: 800; UK: 569; and US: 202) and 1.2 million controls, searching for sequence variants associated with HS., Results: We found 8 independent sequence variants associating with HS, 6 common and 2 rare (frequency <1%). Four associations point to candidate causal genes, NCSTN, PSENEN, WNT10A, and TMED10, that all map to the Notch and Wnt/β-catenin signaling pathways, involved in epidermal keratinization., Limitations: Limited racial diversity may prevent identification of sequence variants of particular importance in non-Caucasian populations., Conclusions: These findings demonstrate that genes and pathways involved in epidermal keratinization are the genetic backbone of HS pathology., Competing Interests: Conflicts of interest Dr Andersen reports that his attendance at the Nordic Congress of DermatoVenerology hosted in Copenhagen in 2022 was funded by Eli Lilly. Drs Stefansdottir, Halldorsson, Ferkingstad, Oddsson, Walters, Olafsdottir, Rutsdottir, Saemundsdottir, Melsted, Norddahl, Sulem, Stefansson, Holm, Gudbjartsson, Thorleifsson, Jonsdottir, and Stefansson are employees of deCODE Genetics/Amgen, Inc. Dr Zachariae has been a speaker/consultant for Jansen Cilag, Novartis, Almirall, UCB, Leo Pharma, Galderma, with no relation to the present manuscript. Dr Thomsen has received research support from Abbvie, Janssen, LEO Pharma, Novartis, Sanofi, and UCB; and has been a speaker/consultant for Abbvie, Eli Lilly, Galderma, Janssen, LEO Pharma, Novartis, Pfizer, Sanofi, Symphogen, UCB, and UNION therapeutics, with no relation to the present manuscript. Dr Brodersen's spouse is a principal level employee at LEO Pharma, part of stock incentive program, and currently hold incentive stock options to a value of approx. 7300USD (January 17 2024). Dr Brunak declares that he has ownerships in Intomics A/S, Hoba Therapeutics Aps, Novo Nordisk A/S, Lundbeck A/S, ALK abello A/S, and managing board memberships in Proscion A/S and Intomics A/S. Dr Erikstrup reports unrestricted departmental grants from Novo Nordisk and Abbott Diagnostics, no personal fees. Dr Eidsmo has received speaker honoraria from Novartis, Lilly, Leo Pharma, and Janssen; and has been a consultant for advisory boards and received research grants from Novartis. Dr Saunte has over the past 3 years received honoraria as a consultant for advisory board meetings by AbbVie, Janssen, Sanofi, Leo Pharma, and Novartis and as a speaker and/or received grants from the following companies: AbbVie, Janssen, Novartis, Sanofi, Jamjoom Pharma, UCB, Pfizer, and Leo Pharma. Dr Jemec has received honoraria from AbbVie, Chemocentryx, Coloplast, Incyte, Inflarx, Novartis, Pierre Fabre, and UCB for participation on advisory boards, and grants from AbbVie, AstraZeneca, Inflarx, Janssen-Cilag, Leo Pharma, Novartis, Regeneron, and Sanofi, for participation as an investigator, and received speaker honoraria from AbbVie, Boehringer Ingelheim, Galderma, and MSD. He has also received unrestricted departmental grants from AbbVie, Leo Pharma, and Novartis. Drs Riis, Dinh, Knowlton, Knight, Nadauld, Banasik, Hansen, Hjalgrim, Sørensen, Mikkelsen, Ullum, Nyegaard, Bruun, Ostrowski, Sigurgeirsson, Orvar, and Pedersen have no conflicts of interest to declare., (Copyright © 2024 American Academy of Dermatology, Inc. Published by Elsevier Inc. All rights reserved.)

Age at menopause (AOM) has a substantial impact on fertility and disease risk. While many loci with variants that associate with AOM have been identified through genome-wide association studies (GWAS) under an additive model, other genetic models are rarely considered 1 . Here through GWAS meta-analysis under the recessive model of 174,329 postmenopausal women from Iceland, Denmark, the United Kingdom (UK; UK Biobank) and Norway, we study low-frequency variants with a large effect on AOM. We discovered that women homozygous for the stop-gain variant rs117316434 (A) in CCDC201 (p.(Arg162Ter), minor allele frequency ~1%) reached menopause 9 years earlier than other women (P = 1.3 × 10 -15 ). The genotype is present in one in 10,000 northern European women and leads to primary ovarian insufficiency in close to half of them. Consequently, homozygotes have fewer children, and the age at last childbirth is 5 years earlier (P = 3.8 × 10 -5 ). The CCDC201 gene was only found in humans in 2022 and is highly expressed in oocytes. Homozygosity for CCDC201 loss-of-function has a substantial impact on female reproductive health, and homozygotes would benefit from reproductive counseling and treatment for symptoms of early menopause., (© 2024. The Author(s).)

Background: In 2021, the American College of Medical Genetics and Genomics (ACMG) recommended reporting actionable genotypes in 73 genes associated with diseases for which preventive or therapeutic measures are available. Evaluations of the association of actionable genotypes in these genes with life span are currently lacking., Methods: We assessed the prevalence of coding and splice variants in genes on the ACMG Secondary Findings, version 3.0 (ACMG SF v3.0), list in the genomes of 57,933 Icelanders. We assigned pathogenicity to all reviewed variants using reported evidence in the ClinVar database, the frequency of variants, and their associations with disease to create a manually curated set of actionable genotypes (variants). We assessed the relationship between these genotypes and life span and further examined the specific causes of death among carriers., Results: Through manual curation of 4405 sequence variants in the ACMG SF v3.0 genes, we identified 235 actionable genotypes in 53 genes. Of the 57,933 participants, 2306 (4.0%) carried at least one actionable genotype. We found shorter median survival among persons carrying actionable genotypes than among noncarriers. Specifically, we found that carrying an actionable genotype in a cancer gene was associated with survival that was 3 years shorter than that among noncarriers, with causes of death among carriers attributed primarily to cancer-related conditions. Furthermore, we found evidence of association between carrying an actionable genotype in certain genes in the cardiovascular disease group and a reduced life span., Conclusions: On the basis of the ACMG SF v3.0 guidelines, we found that approximately 1 in 25 Icelanders carried an actionable genotype and that carrying such a genotype was associated with a reduced life span. (Funded by deCODE Genetics-Amgen.)., (Copyright © 2023 Massachusetts Medical Society.)

Migraine is a complex neurovascular disease with a range of severity and symptoms, yet mostly studied as one phenotype in genome-wide association studies (GWAS). Here we combine large GWAS datasets from six European populations to study the main migraine subtypes, migraine with aura (MA) and migraine without aura (MO). We identified four new MA-associated variants (in PRRT2, PALMD, ABO and LRRK2) and classified 13 MO-associated variants. Rare variants with large effects highlight three genes. A rare frameshift variant in brain-expressed PRRT2 confers large risk of MA and epilepsy, but not MO. A burden test of rare loss-of-function variants in SCN11A, encoding a neuron-expressed sodium channel with a key role in pain sensation, shows strong protection against migraine. Finally, a rare variant with cis-regulatory effects on KCNK5 confers large protection against migraine and brain aneurysms. Our findings offer new insights with therapeutic potential into the complex biology of migraine and its subtypes., (© 2023. The Author(s).)

Background: The TNM system is used to assess prognosis after colorectal cancer (CRC) diagnosis. Other prognostic factors reported include histopathological assessments of the tumour, tumour mutations and proteins in the blood. As some of these factors are strongly correlated, it is important to evaluate the independent effects they may have on survival., Methods: Tumour samples from 2162 CRC patients were visually assessed for amount of tumour stroma, severity of lymphocytic infiltrate at the tumour margins and the presence of lymphoid follicles. Somatic mutations in the tumour were assessed for 2134 individuals. Pre-surgical levels of 4963 plasma proteins were measured in 128 individuals. The associations between these features and prognosis were inspected by a Cox Proportional Hazards Model (CPH)., Results: Levels of stroma, lymphocytic infiltration and presence of lymphoid follicles all associate with prognosis, along with high tumour mutation burden, high microsatellite instability and TP53 and BRAF mutations. The somatic mutations are correlated with the histopathology and none of the somatic mutations associate with survival in a multivariate analysis. Amount of stroma and lymphocytic infiltration associate with local invasion of tumours. Elevated levels of two plasma proteins, CA-125 and PPP1R1A, associate with a worse prognosis., Conclusions: Tumour stroma and lymphocytic infiltration variables are strongly associated with prognosis of CRC and capture the prognostic effects of tumour mutation status. CA-125 and PPP1R1A may be useful prognostic biomarkers in CRC., (© 2023. The Author(s).)

High-throughput proteomics platforms measuring thousands of proteins in plasma combined with genomic and phenotypic information have the power to bridge the gap between the genome and diseases. Here we performed association studies of Olink Explore 3072 data generated by the UK Biobank Pharma Proteomics Project 1 on plasma samples from more than 50,000 UK Biobank participants with phenotypic and genotypic data, stratifying on British or Irish, African and South Asian ancestries. We compared the results with those of a SomaScan v4 study on plasma from 36,000 Icelandic people 2 , for 1,514 of whom Olink data were also available. We found modest correlation between the two platforms. Although cis protein quantitative trait loci were detected for a similar absolute number of assays on the two platforms (2,101 on Olink versus 2,120 on SomaScan), the proportion of assays with such supporting evidence for assay performance was higher on the Olink platform (72% versus 43%). A considerable number of proteins had genomic associations that differed between the platforms. We provide examples where differences between platforms may influence conclusions drawn from the integration of protein levels with the study of diseases. We demonstrate how leveraging the diverse ancestries of participants in the UK Biobank helps to detect novel associations and refine genomic location. Our results show the value of the information provided by the two most commonly used high-throughput proteomics platforms and demonstrate the differences between them that at times provides useful complementarity., (© 2023. The Author(s).)

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.

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

Genotypes causing pregnancy loss and perinatal mortality are depleted among living individuals and are therefore difficult to find. To explore genetic causes of recessive lethality, we searched for sequence variants with deficit of homozygosity among 1.52 million individuals from six European populations. In this study, we identified 25 genes harboring protein-altering sequence variants with a strong deficit of homozygosity (10% or less of predicted homozygotes). Sequence variants in 12 of the genes cause Mendelian disease under a recessive mode of inheritance, two under a dominant mode, but variants in the remaining 11 have not been reported to cause disease. Sequence variants with a strong deficit of homozygosity are over-represented among genes essential for growth of human cell lines and genes orthologous to mouse genes known to affect viability. The function of these genes gives insight into the genetics of intrauterine lethality. We also identified 1077 genes with homozygous predicted loss-of-function genotypes not previously described, bringing the total set of genes completely knocked out in humans to 4785., (© 2023. The Author(s).)

Background and Aims: The causal contribution of apolipoprotein B (apoB) particles to coronary artery disease (CAD) is established. We examined whether this atherogenic contribution is better reflected by non-high-density lipoprotein cholesterol (non-HDL-C) or apoB particle concentration., Method and Results: We performed Mendelian randomization (MR) analysis using 235 variants as genetic instruments; testing the relationship between their effects on the exposures, non-HDL-C and apoB, and on the outcome CAD using weighted regression. Variant effect estimates on the exposures came from the UK Biobank (N = 376 336) and on the outcome from a meta-analysis of five CAD datasets (187 451 cases and 793 315 controls). Subsequently, we carried out sensitivity and replication analyses.In univariate MR analysis, both exposures associated with CAD (βnon-HDL-C = 0.40, P = 2.8 × 10-48 and βapoB = 0.38, P = 1.3 × 10-44). Adding effects on non-HDL-C into a model that already included those on apoB significantly improved the genetically predicted CAD effects (P = 3.9 × 10-5), while adding apoB into the model including non-HDL-C did not (P = 0.69). Thirty-five per cent (82/235) of the variants used as genetic instruments had discordant effects on the exposures, associating with non-HDL-C/apoB ratio at P < 2.1 × 10-4 (0.05/235). Fifty-one variants associated at genome-wide significance., Conclusion: Many sequence variants have discordant effects on non-HDL-C and apoB. These variants allowed us to show that the causal mechanism underlying the relationship between apolipoprotein B particles and CAD is more associated with non-HDL-C than apoB particle concentration., Competing Interests: Conflict of interest: The following authors affiliated with deCODE genetics/Amgen, Inc., are employed by the company: A.H., G.Th., A.S., L.S., G.S., V.T., E.B., V.S., S.G., H.H., J.S., P.S., D.O.A., G.T., U.T., D.F.G., K.S, and H.H., (© The Author(s) 2022. Published by Oxford University Press on behalf of the European Society of Cardiology.)

Memory T-cell responses following SARS-CoV-2 infection have been extensively investigated but many studies have been small with a limited range of disease severity. Here we analyze SARS-CoV-2 reactive T-cell responses in 768 convalescent SARS-CoV-2-infected (cases) and 500 uninfected (controls) Icelanders. The T-cell responses are stable three to eight months after SARS-CoV-2 infection, irrespective of disease severity and even those with the mildest symptoms induce broad and persistent T-cell responses. Robust CD4 + T-cell responses are detected against all measured proteins (M, N, S and S1) while the N protein induces strongest CD8 + T-cell responses. CD4 + T-cell responses correlate with disease severity, humoral responses and age, whereas CD8 + T-cell responses correlate with age and functional antibodies. Further, CD8 + T-cell responses associate with several class I HLA alleles. Our results, provide new insight into HLA restriction of CD8 + T-cell immunity and other factors contributing to heterogeneity of T-cell responses following SARS-CoV-2 infection., (© 2022. The Author(s).)

Detailed knowledge of how diversity in the sequence of the human genome affects phenotypic diversity depends on a comprehensive and reliable characterization of both sequences and phenotypic variation. Over the past decade, insights into this relationship have been obtained from whole-exome sequencing or whole-genome sequencing of large cohorts with rich phenotypic data 1,2 . Here we describe the analysis of whole-genome sequencing of 150,119 individuals from the UK Biobank 3 . This constitutes a set of high-quality variants, including 585,040,410 single-nucleotide polymorphisms, representing 7.0% of all possible human single-nucleotide polymorphisms, and 58,707,036 indels. This large set of variants allows us to characterize selection based on sequence variation within a population through a depletion rank score of windows along the genome. Depletion rank analysis shows that coding exons represent a small fraction of regions in the genome subject to strong sequence conservation. We define three cohorts within the UK Biobank: a large British Irish cohort, a smaller African cohort and a South Asian cohort. A haplotype reference panel is provided that allows reliable imputation of most variants carried by three or more sequenced individuals. We identified 895,055 structural variants and 2,536,688 microsatellites, groups of variants typically excluded from large-scale whole-genome sequencing studies. Using this formidable new resource, we provide several examples of trait associations for rare variants with large effects not found previously through studies based on whole-exome sequencing and/or imputation., (© 2022. The Author(s).)

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

Objectives: The spread of SARS-CoV-2 is dependent on several factors, both biological and behavioural. The effectiveness of nonpharmaceutical interventions can be attributed largely to changes in human behaviour, but quantifying this effect remains challenging. Reconstructing the transmission tree of the third wave of SARS-CoV-2 infections in Iceland using contact tracing and viral sequence data from 2522 cases enables us to directly compare the infectiousness of distinct groups of persons., Methods: The transmission tree enables us to model the effect that a given population prevalence of vaccination would have had on the third wave had one of three different vaccination strategies been implemented before that time. This allows us to compare the effectiveness of the strategies in terms of minimizing the number of cases, deaths, critical cases, and severe cases., Results: We found that people diagnosed outside of quarantine (Rˆ=1.31) were 89% more infectious than those diagnosed while in quarantine (Rˆ=0.70) and that infectiousness decreased as a function of time spent in quarantine before diagnosis, with people diagnosed outside of quarantine being 144% more infectious than those diagnosed after ≥3 days in quarantine (Rˆ=0.54). People of working age, 16 to 66 years (Rˆ=1.08), were 46% more infectious than those outside of that age range (Rˆ=0.74)., Discussion: We found that vaccinating the population in order of ascending age or uniformly at random would have prevented more infections per vaccination than vaccinating in order of descending age, without significantly affecting the expected number of deaths, critical cases, or severe cases., (Copyright © 2022 European Society of Clinical Microbiology and Infectious Diseases. Published by Elsevier Ltd. All rights reserved.)

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

The plasma proteome can help bridge the gap between the genome and diseases. Here we describe genome-wide association studies (GWASs) of plasma protein levels measured with 4,907 aptamers in 35,559 Icelanders. We found 18,084 associations between sequence variants and levels of proteins in plasma (protein quantitative trait loci; pQTL), of which 19% were with rare variants (minor allele frequency (MAF) < 1%). We tested plasma protein levels for association with 373 diseases and other traits and identified 257,490 associations. We integrated pQTL and genetic associations with diseases and other traits and found that 12% of 45,334 lead associations in the GWAS Catalog are with variants in high linkage disequilibrium with pQTL. We identified 938 genes encoding potential drug targets with variants that influence levels of possible biomarkers. Combining proteomics, genomics and transcriptomics, we provide a valuable resource that can be used to improve understanding of disease pathogenesis and to assist with drug discovery and development., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Malignant hyperthermia (MH) susceptibility is a rare life-threatening disorder that occurs upon exposure to a triggering agent. MH is commonly due to protein-altering variants in RYR1 and CACNA1S. The American College of Medical Genetics and Genomics recommends that when pathogenic and likely pathogenic variants in RYR1 and CACNA1S are incidentally found, they should be reported to the carriers. The detection of actionable variants allows the avoidance of exposure to triggering agents during anesthesia. First, we report a 10-year-old Icelandic proband with a suspected MH event, harboring a heterozygous missense variant NM&#95;000540.2:c.6710G>A r.(6710g>a) p.(Cys2237Tyr) in the RYR1 gene that is likely pathogenic. The variant is private to four individuals within a three-generation family and absent from 62,240 whole-genome sequenced (WGS) Icelanders. Haplotype sharing and WGS revealed that the variant occurred as a somatic mosaicism also present in germline of the proband's paternal grandmother. Second, using a set of 62,240 Icelanders with WGS, we assessed the carrier frequency of actionable pathogenic and likely pathogenic variants in RYR1 and CACNA1S. We observed 13 actionable variants in RYR1, based on ClinVar classifications, carried by 43 Icelanders, and no actionable variant in CACNA1S. One in 1450 Icelanders carries an actionable variant for MH. Extensive sequencing allows for better classification and precise dating of variants, and WGS of a large fraction of the population has led to incidental findings of actionable MH genotypes., (© 2021. The Author(s).)

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.

A pressing concern in the SARS-CoV-2 epidemic and other viral outbreaks, is the extent to which the containment measures are halting the viral spread. A straightforward way to assess this is to tally the active cases and the recovered ones throughout the epidemic. Here, we show how epidemic control can be assessed with molecular information during a well characterized epidemic in Iceland. We demonstrate how the viral concentration decreased in those newly diagnosed as the epidemic transitioned from exponential growth phase to containment phase. The viral concentration in the cases identified in population screening decreased faster than in those symptomatic and considered at high risk and that were targeted by the healthcare system. The viral concentration persists in recovering individuals as we found that half of the cases are still positive after two weeks. We demonstrate that accumulation of mutations in SARS-CoV-2 genome can be exploited to track the rate of new viral generations throughout the different phases of the epidemic, where the accumulation of mutations decreases as the transmission rate decreases in the containment phase. Overall, the molecular signatures of SARS-CoV-2 infections contain valuable epidemiological information that can be used to assess the effectiveness of containment measures.

Background: Loss-of-function mutations in the LDL (low-density lipoprotein) receptor gene ( LDLR ) cause elevated levels of LDL cholesterol and premature cardiovascular disease. To date, a gain-of-function mutation in LDLR with a large effect on LDL cholesterol levels has not been described. Here, we searched for sequence variants in LDLR that have a large effect on LDL cholesterol levels., Methods: We analyzed whole-genome sequencing data from 43 202 Icelanders. Single-nucleotide polymorphisms and structural variants including deletions, insertions, and duplications were genotyped using whole-genome sequencing-based data. LDL cholesterol associations were carried out in a sample of >100 000 Icelanders with genetic information (imputed or whole-genome sequencing). Molecular analyses were performed using RNA sequencing and protein expression assays in Epstein-Barr virus-transformed lymphocytes., Results: We discovered a 2.5-kb deletion (del2.5) overlapping the 3' untranslated region of LDLR in 7 heterozygous carriers from a single family. Mean level of LDL cholesterol was 74% lower in del2.5 carriers than in 101 851 noncarriers, a difference of 2.48 mmol/L (96 mg/dL; P =8.4×10 - 8 ). Del2.5 results in production of an alternative mRNA isoform with a truncated 3' untranslated region. The truncation leads to a loss of target sites for microRNAs known to repress translation of LDLR . In Epstein-Barr virus-transformed lymphocytes derived from del2.5 carriers, expression of alternative mRNA isoform was 1.84-fold higher than the wild-type isoform ( P =0.0013), and there was 1.79-fold higher surface expression of the LDL receptor than in noncarriers ( P =0.0086). We did not find a highly penetrant detrimental impact of lifelong very low levels of LDL cholesterol due to del2.5 on health of the carriers., Conclusions: Del2.5 is the first reported gain-of-function mutation in LDLR causing a large reduction in LDL cholesterol. These data point to a role for alternative polyadenylation of LDLR mRNA as a potent regulator of LDL receptor expression in humans.

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

Background: During the current worldwide pandemic, coronavirus disease 2019 (Covid-19) was first diagnosed in Iceland at the end of February. However, data are limited on how SARS-CoV-2, the virus that causes Covid-19, enters and spreads in a population., Methods: We targeted testing to persons living in Iceland who were at high risk for infection (mainly those who were symptomatic, had recently traveled to high-risk countries, or had contact with infected persons). We also carried out population screening using two strategies: issuing an open invitation to 10,797 persons and sending random invitations to 2283 persons. We sequenced SARS-CoV-2 from 643 samples., Results: As of April 4, a total of 1221 of 9199 persons (13.3%) who were recruited for targeted testing had positive results for infection with SARS-CoV-2. Of those tested in the general population, 87 (0.8%) in the open-invitation screening and 13 (0.6%) in the random-population screening tested positive for the virus. In total, 6% of the population was screened. Most persons in the targeted-testing group who received positive tests early in the study had recently traveled internationally, in contrast to those who tested positive later in the study. Children under 10 years of age were less likely to receive a positive result than were persons 10 years of age or older, with percentages of 6.7% and 13.7%, respectively, for targeted testing; in the population screening, no child under 10 years of age had a positive result, as compared with 0.8% of those 10 years of age or older. Fewer females than males received positive results both in targeted testing (11.0% vs. 16.7%) and in population screening (0.6% vs. 0.9%). The haplotypes of the sequenced SARS-CoV-2 viruses were diverse and changed over time. The percentage of infected participants that was determined through population screening remained stable for the 20-day duration of screening., Conclusions: In a population-based study in Iceland, children under 10 years of age and females had a lower incidence of SARS-CoV-2 infection than adolescents or adults and males. The proportion of infected persons identified through population screening did not change substantially during the screening period, which was consistent with a beneficial effect of containment efforts. (Funded by deCODE Genetics-Amgen.)., (Copyright © 2020 Massachusetts Medical Society.)

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.

Background: Alpha-fetoprotein (AFP), cancer antigens 15.3, 19.9, and 125, carcinoembryonic antigen, and alkaline phosphatase (ALP) are widely measured in attempts to detect cancer and to monitor treatment response. However, due to lack of sensitivity and specificity, their utility is debated. The serum levels of these markers are affected by a number of nonmalignant factors, including genotype. Thus, it may be possible to improve both sensitivity and specificity by adjusting test results for genetic effects., Methods: We performed genome-wide association studies of serum levels of AFP ( N = 22,686), carcinoembryonic antigen ( N = 22,309), cancer antigens 15.3 ( N = 7,107), 19.9 ( N = 9,945), and 125 ( N = 9,824), and ALP ( N = 162,774). We also examined the correlations between levels of these biomarkers and the presence of cancer, using data from a nationwide cancer registry., Results: We report a total of 84 associations of 79 sequence variants with levels of the six biomarkers, explaining between 2.3% and 42.3% of the phenotypic variance. Among the 79 variants, 22 are cis (in- or near the gene encoding the biomarker), 18 have minor allele frequency less than 1%, 31 are coding variants, and 7 are associated with gene expression in whole blood. We also find multiple conditions associated with higher biomarker levels., Conclusions: Our results provide insights into the genetic contribution to diversity in concentration of tumor biomarkers in blood., Impact: Genetic correction of biomarker values could improve prediction algorithms and decision-making based on these biomarkers., (©2019 American Association for Cancer Research.)

Background: Lipoprotein(a) [Lp(a)] is a causal risk factor for cardiovascular diseases that has no established therapy. The attribute of Lp(a) that affects cardiovascular risk is not established. Low levels of Lp(a) have been associated with type 2 diabetes (T2D)., Objectives: This study investigated whether cardiovascular risk is conferred by Lp(a) molar concentration or apolipoprotein(a) [apo(a)] size, and whether the relationship between Lp(a) and T2D risk is causal., Methods: This was a case-control study of 143,087 Icelanders with genetic information, including 17,715 with coronary artery disease (CAD) and 8,734 with T2D. This study used measured and genetically imputed Lp(a) molar concentration, kringle IV type 2 (KIV-2) repeats (which determine apo(a) size), and a splice variant in LPA associated with small apo(a) but low Lp(a) molar concentration to disentangle the relationship between Lp(a) and cardiovascular risk. Loss-of-function homozygotes and other subjects genetically predicted to have low Lp(a) levels were evaluated to assess the relationship between Lp(a) and T2D., Results: Lp(a) molar concentration was associated dose-dependently with CAD risk, peripheral artery disease, aortic valve stenosis, heart failure, and lifespan. Lp(a) molar concentration fully explained the Lp(a) association with CAD, and there was no residual association with apo(a) size. Homozygous carriers of loss-of-function mutations had little or no Lp(a) and increased the risk of T2D., Conclusions: Molar concentration is the attribute of Lp(a) that affects risk of cardiovascular diseases. Low Lp(a) concentration (bottom 10%) increases T2D risk. Pharmacologic reduction of Lp(a) concentration in the 20% of individuals with the greatest concentration down to the population median is predicted to decrease CAD risk without increasing T2D risk., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

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.

Background: Most pathogenic mutations in the BRCA2 gene carry a high risk of hereditary breast and ovarian cancer (HBOC). However, a stop-gain mutation, K3326* (rs11571833), confers risk of lung cancer and cancers of the upper-aero-digestive tract but only a modest risk of breast or ovarian cancer. The Icelandic population provides an opportunity for comprehensive characterization of the cancer risk profiles of K3326* and HBOC mutations because a single mutation, BRCA2 999del5, is responsible for almost all BRCA2-related HBOC in the population., Methods: Genotype information on 43 641 cancer patients and 370 971 control subjects from Iceland, the Netherlands, and the United States was used to assess the cancer risk profiles of K3326* and BRCA2 999del5. BRCA2 expression was assessed using RNAseq data from blood (n = 2233), as well as 52 tissues reported in the GTEx database., Results: The cancer risks associated with K3326* are fundamentally different from those associated with 999del5. We report for the first time an association between K3326* and small cell lung cancer (odds ratio [OR] = 2.06, 95% confidence interval [CI] = 1.35 to 3.16) and squamous cell carcinoma of the skin (OR = 1.69, 95% CI = 1.26 to 2.26). Individuals homozygous for K3326* reach old age and have children. Unlike BRCA2 999del5, the K3326* allele does not affect the level of BRCA2 transcripts, and the allele is expressed to the same extent as the wild-type allele., Conclusions: K3326* associates primarily with cancers that have strong environmental genotoxic risk factors. Expression of the K3326* allele suggests that a variant protein may be made that retains the DNA repair capabilities important to hormone-responsive tissues but may be less efficient in responding to genotoxic stress.

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.

Background: Rare missense mutations in the gene encoding coatomer subunit alpha (COPA) have recently been shown to cause autoimmune interstitial lung, joint and kidney disease, also known as COPA syndrome, under a dominant mode of inheritance., Case Presentation: Here we describe an Icelandic family with three affected individuals over two generations with a rare clinical presentation of lung and joint disease and a histological diagnosis of follicular bronchiolitis. We performed whole-genome sequencing (WGS) of the three affected as well as three unaffected members of the family, and searched for rare genotypes associated with disease using 30,067 sequenced Icelanders as a reference population. We assessed all coding and splicing variants, prioritizing variants in genes known to cause interstitial lung disease. We detected a heterozygous missense mutation, p.Glu241Lys, in the COPA gene, private to the affected family members. The mutation occurred de novo in the paternal germline of the index case and was absent from 30,067 Icelandic genomes and 141,353 individuals from the genome Aggregation Database (gnomAD). The mutation occurs within the conserved and functionally important WD40 domain of the COPA protein., Conclusions: This is the second report of the p.Glu241Lys mutation in COPA, indicating the recurrent nature of the mutation. The mutation was reported to co-segregate with COPA syndrome in a large family from the USA with five affected members, and classified as pathogenic. The two separate occurrences of the p.Glu241Lys mutation in cases and its absence from a large number of sequenced genomes confirms its role in the pathogenesis of the COPA syndrome.

Background: Epileptic encephalopathies are a group of childhood epilepsies that display high phenotypic and genetic heterogeneity. The recent, extensive use of next-generation sequencing has identified a large number of genes in epileptic encephalopathies, including UBA5 in which biallelic mutations were first described as pathogenic in 2016 (Colin E et al., Am J Hum Genet 99(3):695-703, 2016. Muona M et al., Am J Hum Genet 99(3):683-694, 2016). UBA5 encodes an activating enzyme for a post-translational modification mechanism known as ufmylation, and is the first gene from the ufmylation pathway that is linked to disease., Case Presentation: We sequenced the genomes of two sisters with early-onset epileptic encephalopathy along with their unaffected parents in an attempt to find a genetic cause for their condition. The sisters, born in 2004 and 2006, presented with infantile spasms at six months of age, which later progressed to recurrent, treatment-resistant seizures. We detected a compound heterozygous genotype in UBA5 in the sisters, a genotype not seen elsewhere in an Icelandic reference set of 30,067 individuals nor in public databases. One of the mutations, c.684G > A, is a paternally inherited exonic splicing mutation, occuring at the last nucleotide of exon 7 of UBA5. The mutation is predicted to disrupt the splice site, resulting in loss-of-function of one allele of UBA5. The second mutation is a maternally inherited missense mutation, p.Ala371Thr, previously reported as pathogenic when in compound heterozygosity with a loss-of-function mutation in UBA5 and is believed to produce a hypomorphic allele. Supportive of this, we have identified three adult Icelanders homozygous for the p.Ala371Thr mutation who show no signs of neurological disease., Conclusions: We describe compound heterozygous mutations in the UBA5 gene in two sisters with early-onset epileptic encephalopathy. To our knowledge, this is the first description of mutations in UBA5 since the initial discovery that pathogenic biallelic variants in the gene cause early-onset epileptic encephalopathy. We further provide confirmatory evidence that p.Ala371Thr is a hypomorphic mutation, by presenting three adult homozygotes who show no signs of neurological disease.

Aims: Atrial fibrillation (AF) is the most common sustained cardiac arrhythmia in man, causing substantial morbidity and mortality with a major worldwide public health impact. It is increasingly recognized as a highly heritable condition. This study aimed to determine genetic risk factors for early-onset AF., Methods and Results: We sequenced the whole genomes of 8453 Icelanders and imputed genotypes of the 25.5 million sequence variants we discovered into 1799 Icelanders with early-onset AF (diagnosed before 60 years of age) and 337 453 controls. Each sequence variant was tested for association based on multiplicative and recessive inheritance models. We discovered a rare frameshift deletion in the myosin MYL4 gene (c.234delC) that associates with early-onset AF under a recessive mode of inheritance (allelic frequency = 0.58%). We found eight homozygous carriers of the mutation, all of whom had early-onset AF. Six of the homozygotes were diagnosed by the age of 30 and the remaining two in their 50s. Three of the homozygotes had received pacemaker implantations due to sick sinus syndrome, three had suffered an ischemic stroke, and one suffered sudden cardiac death., Conclusions: Through a population approach we found a loss of function mutation in the myosin gene MYL4 that, in the homozygous state, is completely penetrant for early-onset AF. The finding may provide novel mechanistic insight into the pathophysiology of this complex arrhythmia., (Published on behalf of the European Society of Cardiology. All rights reserved. © The Author 2016. For Permissions, please email: journals.permissions@oup.com.)

Here we describe the insights gained from sequencing the whole genomes of 2,636 Icelanders to a median depth of 20×. We found 20 million SNPs and 1.5 million insertions-deletions (indels). We describe the density and frequency spectra of sequence variants in relation to their functional annotation, gene position, pathway and conservation score. We demonstrate an excess of homozygosity and rare protein-coding variants in Iceland. We imputed these variants into 104,220 individuals down to a minor allele frequency of 0.1% and found a recessive frameshift mutation in MYL4 that causes early-onset atrial fibrillation, several mutations in ABCB4 that increase risk of liver diseases and an intronic variant in GNAS associating with increased thyroid-stimulating hormone levels when maternally inherited. These data provide a study design that can be used to determine how variation in the sequence of the human genome gives rise to human diversity.