Search

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. The concentration of SARS-CoV-2 changes during an individual's infection, and mutations accumulate as viruses are transmitted between people. Here, the authors use data from Iceland to demonstrate how this information can be exploited at the population-level to determine the phase of the epidemic. [ABSTRACT FROM AUTHOR]

We performed a genome-wide association study on 1,292 individuals with abdominal aortic aneurysms (AAAs) and 30,503 controls from Iceland and The Netherlands, with a follow-up of top markers in up to 3,267 individuals with AAAs and 7,451 controls. The A allele of rs7025486 on 9q33 was found to associate with AAA, with an odds ratio (OR) of 1.21 and P = 4.6 × 10−10. In tests for association with other vascular diseases, we found that rs7025486[A] is associated with early onset myocardial infarction (OR = 1.18, P = 3.1 × 10−5), peripheral arterial disease (OR = 1.14, P = 3.9 × 10−5) and pulmonary embolism (OR = 1.20, P = 0.00030), but not with intracranial aneurysm or ischemic stroke. No association was observed between rs7025486[A] and common risk factors for arterial and venous diseases—that is, smoking, lipid levels, obesity, type 2 diabetes and hypertension. Rs7025486 is located within DAB2IP, which encodes an inhibitor of cell growth and survival. [ABSTRACT FROM AUTHOR]

Chronic kidney disease (CKD) is a worldwide public health problem that is associated with substantial morbidity and mortality. To search for sequence variants that associate with CKD, we conducted a genome-wide association study (GWAS) that included a total of 3,203 Icelandic cases and 38,782 controls. We observed an association between CKD and a variant with 80% population frequency, rs4293393-T, positioned next to the UMOD gene (GeneID: 7369) on chromosome 16p12 (OR = 1.25, P= 4.1×10-10). This gene encodes uromodulin (Tamm-Horsfall protein), the most abundant protein in mammalian urine. The variant also associates significantly with serum creatinine concentration (SCr) in Icelandic subjects (N = 24,635, P= 1.3×10-23) but not in a smaller set of healthy Dutch controls (N = 1,819, P = 0.39). Our findings validate the association between the UMOD variant and both CKD and SCr recently discovered in a large GWAS. In the Icelandic dataset, we demonstrate that the effect on SCr increases substantially with both age (P= 3.0×10-17) and number of comorbid diseases (P = 0.008). The association with CKD is also stronger in the older age groups. These results suggest that the UMOD variant may influence the adaptation of the kidney to age-related risk factors of kidney disease such as hypertension and diabetes. The variant also associates with serum urea (P =1.0×10-6), uric acid (P = 0.0064), and suggestively with gout. In contrast to CKD, the UMOD variant confers protection against kidney stones when studied in 3,617 Icelandic and Dutch kidney stone cases and 43,201 controls (OR = 0.88, P= 5.7×10-5). [ABSTRACT FROM AUTHOR]

Kidney stone disease is a common condition. To search for sequence variants conferring risk of kidney stones, we conducted a genome-wide association study in 3,773 cases and 42,510 controls from Iceland and The Netherlands. We discovered common, synonymous variants in the CLDN14 gene that associate with kidney stones (OR = 1.25 and P = 4.0 × 10−12 for rs219780[C]). Approximately 62% of the general population is homozygous for rs219780[C] and is estimated to have 1.64 times greater risk of developing the disease compared to noncarriers. The CLDN14 gene is expressed in the kidney and regulates paracellular permeability at epithelial tight junctions. The same variants were also found to associate with reduced bone mineral density at the hip (P = 0.00039) and spine (P = 0.0077). [ABSTRACT FROM AUTHOR]

We expanded our genome-wide association study on atrial fibrillation (AF) in Iceland, which previously identified risk variants on 4q25, and tested the most significant associations in samples from Iceland, Norway and the United States. A variant in the ZFHX3 gene on chromosome 16q22, rs7193343-T, associated significantly with AF (odds ratio OR = 1.21, P = 1.4 × 10−10). This variant also associated with ischemic stroke (OR = 1.11, P = 0.00054) and cardioembolic stroke (OR = 1.22, P = 0.00021) in a combined analysis of five stroke samples. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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.). [ABSTRACT FROM AUTHOR]

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

Introduction: Long QT syndrome (LQTS) is a well-known cause of sudden cardiac death (SCD) in the young. The frequency of LQTS is estimated to be 1:2000. The most common causes of LQTS are variants in KCNQ1 , KCNH2 and SCN5A. The phenotypic effects of such variants has not been well characterized. Hypothesis: The goal of the study was to identify and estimate the frequency of variants causing LQTS in the Icelandic population and to determine if the genotype can predict the severity of the phenotype. Methods: The sequence diversity in KCNQ1 , KCNH2 and SCN5A was analyzed through large-scale whole genome sequencing of 45K Icelanders and imputation into 150K chip-typed. Using a set of 434,000 ECGs from 80,694 individuals, the effect of sequence variants on the QTc-interval was assessed. To explore the effect on outcomes, we analyzed association of the variants with severe cardiac adverse events, including SCD (N=3,133). Medical records of carriers were evaluated to gain more phenotypic information. Results: We identified six rare coding variants in KCNQ1 and two in KCNH2 that have strong effects on the QTc-interval. One of the KCNQ1 variants, carried by one in 3500 Icelanders, has the largest effect on the QTc-interval (P=3.2х10-25, effect=2.2 SD) while the most frequent variant (found in 70%), carried by one in 300 Icelanders, has the smallest effect (P=3.1х10-34, effect=0.84 SD). The variants were tested for association with adverse events. The strongest association with sudden cardiac death was with two of the KCNQ1 variants (P=0.0013 and 0.0051, ORs=4.3 and 3.9) while the most frequent variant showed no association (P=0.97, OR=1.01). We further identified four variants in KCNH2 and two in KCNQ1 , reported as pathogenic by ClinVar, that are too rare (private) to be tested for association. Conclusions: Screening for variants causing LQTS in Iceland and assessing their phenotypic effects reveals new insight into both genetic epidemiology of LQTS and the relationship between genotype and phenotype. These results can be used for further risk assessment of individuals with LQTS. [ABSTRACT FROM AUTHOR]

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

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

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

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 Long-QT syndrome (LQTS) is a cardiac repolarization abnormality that can lead to sudden cardiac death. The most common causes are rare coding variants in the genes KCNQ1 , KCNH2 , and SCN5A . The data on LQTS epidemiology are limited, and information on expressivity and penetrance of pathogenic variants is sparse. Methods and Results We screened for rare coding variants associated with the corrected QT (QTc) interval in Iceland. We explored the frequency of the identified variants, their penetrance, and their association with severe events. Twelve variants were associated with the QTc interval. Five in KCNQ 1, 3 in KCNH2 , 2 in cardiomyopathy genes MYBPC3 and PKP2 , and 2 in genes where coding variants have not been associated with the QTc interval, ISOC1 and MYOM2. The combined carrier frequency of the 8 variants in the previously known LQTS genes was 530 per 100 000 individuals (1:190). p.Tyr315Cys and p.Leu273Phe in KCNQ1 were associated with having a mean QTc interval longer than 500 ms ( P =4.2×10 -7 ; odds ratio [OR], 38.6; P =8.4×10 -10 , OR, 26.5; respectively), and p.Leu273Phe was associated with sudden cardiac death ( P =0.0034; OR, 2.99). p.Val215Met in KCNQ1 was carried by 1 in 280 Icelanders, had a smaller effect on the QTc interval ( P =1.8×10 -44 ; effect, 22.8 ms), and did not associate with severe clinical events. Conclusions The carrier frequency of associating variants in LQTS genes was higher than previous estimates of the prevalence of LQTS. The variants have variable effects on the QTc interval, and carriers of p.Tyr315Cys and p.Leu273Phe have a more severe disease than carriers of p.Val215Met. These data could lead to improved identification, risk stratification, and a more precise clinical approach to those with QTc prolongation.

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

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

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.

Long-read sequencing (LRS) promises to improve the characterization of structural variants (SVs). We generated LRS data from 3,622 Icelanders and identified a median of 22,636 SVs per individual (a median of 13,353 insertions and 9,474 deletions). We discovered a set of 133,886 reliably genotyped SV alleles and imputed them into 166,281 individuals to explore their effects on diseases and other traits. We discovered an association of a rare deletion in PCSK9 with lower low-density lipoprotein (LDL) cholesterol levels, compared to the population average. We also discovered an association of a multiallelic SV in ACAN with height; we found 11 alleles that differed in the number of a 57-bp-motif repeat and observed a linear relationship between the number of repeats carried and height. These results show that SVs can be accurately characterized at the population scale using LRS data in a genome-wide non-targeted approach and demonstrate how SVs impact phenotypes.

Iron is essential for many biological functions and iron deficiency and overload have major health implications. We performed a meta-analysis of three genome-wide association studies from Iceland, the UK and Denmark of blood levels of ferritin (N = 246,139), total iron binding capacity (N = 135,430), iron (N = 163,511) and transferrin saturation (N = 131,471). We found 62 independent sequence variants associating with iron homeostasis parameters at 56 loci, including 46 novel loci. Variants at DUOX2, F5, SLC11A2 and TMPRSS6 associate with iron deficiency anemia, while variants at TF, HFE, TFR2 and TMPRSS6 associate with iron overload. A HBS1L-MYB intergenic region variant associates both with increased risk of iron overload and reduced risk of iron deficiency anemia. The DUOX2 missense variant is present in 14% of the population, associates with all iron homeostasis biomarkers, and increases the risk of iron deficiency anemia by 29%. The associations implicate proteins contributing to the main physiological processes involved in iron homeostasis: iron sensing and storage, inflammation, absorption of iron from the gut, iron recycling, erythropoiesis and bleeding/menstruation.

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.

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

A major challenge in genetic association studies is that most associated variants fall in the non-coding part of the human genome. We searched for variants associated with bone mineral density (BMD) after enriching the discovery cohort for loss-of-function (LoF) mutations by sequencing a subset of the Nord-Trøndelag Health Study, followed by imputation in the remaining sample (N = 19,705), and identified ten known BMD loci. However, one previously unreported variant, LoF mutation in MEPE, p.(Lys70IlefsTer26, minor allele frequency [MAF] = 0.8%), was associated with decreased ultradistal forearm BMD (P-value = 2.1 × 10 -18 ), and increased osteoporosis (P-value = 4.2 × 10 -5 ) and fracture risk (P-value = 1.6 × 10 -5 ). The MEPE LoF association with BMD and fractures was further evaluated in 279,435 UK (MAF = 0.05%, heel bone estimated BMD P-value = 1.2 × 10 -16 , any fracture P-value = 0.05) and 375,984 Icelandic samples (MAF = 0.03%, arm BMD P-value = 0.12, forearm fracture P-value = 0.005). Screening for the MEPE LoF mutations before adulthood could potentially prevent osteoporosis and fractures due to the lifelong effect on BMD observed in the study. A key implication for precision medicine is that high-impact functional variants missing from the publicly available cosmopolitan panels could be clinically more relevant than polygenic risk scores.

Aims: To explore whether variability in dietary cholesterol and phytosterol absorption impacts the risk of coronary artery disease (CAD) using as instruments sequence variants in the ABCG5/8 genes, key regulators of intestinal absorption of dietary sterols., Methods and Results: We examined the effects of ABCG5/8 variants on non-high-density lipoprotein (non-HDL) cholesterol (N up to 610 532) and phytosterol levels (N = 3039) and the risk of CAD in Iceland, Denmark, and the UK Biobank (105 490 cases and 844 025 controls). We used genetic scores for non-HDL cholesterol to determine whether ABCG5/8 variants confer greater risk of CAD than predicted by their effect on non-HDL cholesterol. We identified nine rare ABCG5/8 coding variants with substantial impact on non-HDL cholesterol. Carriers have elevated phytosterol levels and are at increased risk of CAD. Consistent with impact on ABCG5/8 transporter function in hepatocytes, eight rare ABCG5/8 variants associate with gallstones. A genetic score of ABCG5/8 variants predicting 1 mmol/L increase in non-HDL cholesterol associates with two-fold increase in CAD risk [odds ratio (OR) = 2.01, 95% confidence interval (CI) 1.75-2.31, P = 9.8 × 10-23] compared with a 54% increase in CAD risk (OR = 1.54, 95% CI 1.49-1.59, P = 1.1 × 10-154) associated with a score of other non-HDL cholesterol variants predicting the same increase in non-HDL cholesterol (P for difference in effects = 2.4 × 10-4)., Conclusions: Genetic variation in cholesterol absorption affects levels of circulating non-HDL cholesterol and risk of CAD. Our results indicate that both dietary cholesterol and phytosterols contribute directly to atherogenesis., (© The Author(s) 2020. Published by Oxford University Press on behalf of the European Society of Cardiology.)

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.

Pelvic organ prolapse (POP) is a downward descent of one or more of the pelvic organs, resulting in a protrusion of the vaginal wall and/or uterus. We performed a genome-wide association study of POP using data from Iceland and the UK Biobank, a total of 15,010 cases with hospital-based diagnosis code and 340,734 female controls, and found eight sequence variants at seven loci associating with POP (P < 5 × 10 -8 ); seven common (minor allele frequency >5%) and one with minor allele frequency of 4.87%. Some of the variants associating with POP also associated with traits of similar pathophysiology. Of these, rs3820282, which may alter the estrogen-based regulation of WNT4, also associates with leiomyoma of uterus, gestational duration and endometriosis. Rs3791675 at EFEMP1, a gene involved in connective tissue homeostasis, also associates with hernias and carpal tunnel syndrome. Our results highlight the role of connective tissue metabolism and estrogen exposure in the etiology of POP.

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

Features of the QRS complex of the electrocardiogram, reflecting ventricular depolarisation, associate with various physiologic functions and several pathologic conditions. We test 32.5 million variants for association with ten measures of the QRS complex in 12 leads, using 405,732 electrocardiograms from 81,192 Icelanders. We identify 190 associations at 130 loci, the majority of which have not been reported before, including associations with 21 rare or low-frequency coding variants. Assessment of genes expressed in the heart yields an additional 13 rare QRS coding variants at 12 loci. We find 51 unreported associations between the QRS variants and echocardiographic traits and cardiovascular diseases, including atrial fibrillation, complete AV block, heart failure and supraventricular tachycardia. We demonstrate the advantage of in-depth analysis of the QRS complex in conjunction with other cardiovascular phenotypes to enhance our understanding of the genetic basis of myocardial mass, cardiac conduction and disease.

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

Benign prostatic hyperplasia and associated lower urinary tract symptoms (BPH/LUTS) are common conditions affecting the majority of elderly males. Here we report the results of a genome-wide association study of symptomatic BPH/LUTS in 20,621 patients and 280,541 controls of European ancestry, from Iceland and the UK. We discovered 23 genome-wide significant variants, located at 14 loci. There is little or no overlap between the BPH/LUTS variants and published prostate cancer risk variants. However, 15 of the variants reported here also associate with serum levels of prostate specific antigen (PSA) (at a Bonferroni corrected P < 0.0022). Furthermore, there is a strong genetic correlation, r g  = 0.77 (P = 2.6 × 10 -11 ), between PSA and BPH/LUTS, and one standard deviation increase in a polygenic risk score (PRS) for BPH/LUTS increases PSA levels by 12.9% (P = 1.6×10 -55 ). These results shed a light on the genetic background of BPH/LUTS and its substantial influence on PSA levels.

Imprinting is the preferential expression of one parental allele over the other. It is controlled primarily through differential methylation of cytosine at CpG dinucleotides. Here we combine 285 methylomes and 11,617 transcriptomes from peripheral blood samples with parent-of-origin phased haplotypes, to produce a new map of imprinted methylation and gene expression patterns across the human genome. We demonstrate how imprinted methylation is a continuous rather than a binary characteristic. We describe at high resolution the parent-of-origin methylation pattern at the 15q11.2 Prader-Willi/Angelman syndrome locus, with nearly confluent stochastic paternal methylation punctuated by 'spikes' of maternal methylation. We find examples of polymorphic imprinted methylation unrelated (at VTRNA2-1 and PARD6G) or related (at CHRNE) to nearby SNP genotypes. We observe RNA isoform-specific imprinted expression patterns suggestive of a methylation-sensitive transcriptional elongation block. Finally, we gain new insights into parent-of-origin-specific effects on phenotypes at the DLK1/MEG3 and GNAS loci.

Background: Dilated cardiomyopathy (DCM) is an important cause of heart failure. Variants in >50 genes have been reported to cause DCM, but causative variants have been found in less than half of familial cases. Variants causing DCM in Iceland have not been reported before., Methods: We performed a genome-wide association study on DCM based on whole genome sequencing. We tested the association of 32.5 million sequence variants in 424 cases and 337 689 population controls in Iceland., Results: We identified 2 DCM variants in established cardiomyopathy genes, a missense variant p.Phe145Leu in NKX2-5 carried by 1 in 7100 Icelanders ( P=7.0×10 -12 ) and a frameshift variant p.Phe1626Serfs*40 in FLNC carried by 1 in 3600 Icelanders ( P=2.1×10 -10 ). Both variants associate with heart failure and sudden cardiac death. Additionally, p.Phe145Leu in NKX2-5 associates with high degree atrioventricular block and atrial septal defect ( P<1.4×10 -4 ). The penetrance of serious heart disease among carriers of the NKX2-5 variant is high and higher than that of the FLNC variant., Conclusions: Two rare variants in NKX2-5 and FLNC, carried by 1 in 2400 Icelanders, cause familial DCM in Iceland. These genes have recently been associated with DCM. Given the serious consequences of these variants, we suggest screening for them in individuals with DCM and their family members, with subsequent monitoring of carriers, offering early intervention.

Aims: Scavenger receptor Class B Type 1 (SR-BI) is a major receptor for high-density lipoprotein (HDL) that promotes hepatic uptake of cholesterol from HDL. A rare mutation p.P376L, in the gene encoding SR-BI, SCARB1, was recently reported to associate with elevated HDL cholesterol (HDL-C) and increased risk of coronary artery disease (CAD), suggesting that increased HDL-C caused by SR-BI impairment might be an independent marker of cardiovascular risk. We tested the hypothesis that alleles in or close to SCARB1 that associate with elevated levels of HDL-C also associate with increased risk of CAD in the relatively homogeneous population of Iceland., Methods and Results: Using a large resource of whole-genome sequenced Icelanders, we identified thirteen SCARB1 coding mutations that we examined for association with HDL-C (n = 136 672). Three rare SCARB1 mutations, encoding p.G319V, p.V111M, and p.V32M (combined allelic frequency = 0.2%) associate with elevated levels of HDL-C (p.G319V: β = 11.1 mg/dL, P = 8.0 × 10-7; p.V111M: β = 8.3 mg/dL, P = 1.1 × 10-6; p.V32M: β = 10.2 mg/dL, P = 8.1 × 10-4). These mutations do not associate with CAD (36 886 cases/306 268 controls) (odds ratio = 0.90, 95% confidence interval 0.67-1.22, P = 0.49), despite effects on HDL-C comparable to that reported for p.P376L, both in terms of direction and magnitude. Furthermore, HDL-C raising alleles of three common SCARB1 non-coding variants, including one previously unreported (rs61941676-C: β = 1.25 mg/dL, P = 1.7 × 10-18), and of one low frequency coding variant (p.V135I) that independently associate with higher HDL-C, do not confer increased risk of CAD., Conclusion: Elevated HDL-C due to genetically compromised SR-BI function is not a marker of CAD risk.

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

Background: Several sequence variants are known to have effects on serum levels of non-high-density lipoprotein (HDL) cholesterol that alter the risk of coronary artery disease., Methods: We sequenced the genomes of 2636 Icelanders and found variants that we then imputed into the genomes of approximately 398,000 Icelanders. We tested for association between these imputed variants and non-HDL cholesterol levels in 119,146 samples. We then performed replication testing in two populations of European descent. We assessed the effects of an implicated loss-of-function variant on the risk of coronary artery disease in 42,524 case patients and 249,414 controls from five European ancestry populations. An augmented set of genomes was screened for additional loss-of-function variants in a target gene. We evaluated the effect of an implicated variant on protein stability., Results: We found a rare noncoding 12-base-pair (bp) deletion (del12) in intron 4 of ASGR1, which encodes a subunit of the asialoglycoprotein receptor, a lectin that plays a role in the homeostasis of circulating glycoproteins. The del12 mutation activates a cryptic splice site, leading to a frameshift mutation and a premature stop codon that renders a truncated protein prone to degradation. Heterozygous carriers of the mutation (1 in 120 persons in our study population) had a lower level of non-HDL cholesterol than noncarriers, a difference of 15.3 mg per deciliter (0.40 mmol per liter) (P=1.0×10(-16)), and a lower risk of coronary artery disease (by 34%; 95% confidence interval, 21 to 45; P=4.0×10(-6)). In a larger set of sequenced samples from Icelanders, we found another loss-of-function ASGR1 variant (p.W158X, carried by 1 in 1850 persons) that was also associated with lower levels of non-HDL cholesterol (P=1.8×10(-3))., Conclusions: ASGR1 haploinsufficiency was associated with reduced levels of non-HDL cholesterol and a reduced risk of coronary artery disease. (Funded by the National Institutes of Health and others.).

The consensus approach to genome-wide association studies (GWAS) has been to assign equal prior probability of association to all sequence variants tested. However, some sequence variants, such as loss-of-function and missense variants, are more likely than others to affect protein function and are therefore more likely to be causative. Using data from whole-genome sequencing of 2,636 Icelanders and the association results for 96 quantitative and 123 binary phenotypes, we estimated the enrichment of association signals by sequence annotation. We propose a weighted Bonferroni adjustment that controls for the family-wise error rate (FWER), using as weights the enrichment of sequence annotations among association signals. We show that this weighted adjustment increases the power to detect association over the standard Bonferroni correction. We use the enrichment of associations by sequence annotation we have estimated in Iceland to derive significance thresholds for other populations with different numbers and combinations of sequence variants.

Creatine kinase (CK) and lactate dehydrogenase (LDH) are widely used markers of tissue damage. To search for sequence variants influencing serum levels of CK and LDH, 28.3 million sequence variants identified through whole-genome sequencing of 2,636 Icelanders were imputed into 63,159 and 98,585 people with CK and LDH measurements, respectively. Here we describe 13 variants associating with serum CK and 16 with LDH levels, including four that associate with both. Among those, 15 are non-synonymous variants and 12 have a minor allele frequency below 5%. We report sequence variants in genes encoding the enzymes being measured (CKM and LDHA), as well as in genes linked to muscular (ANO5) and immune/inflammatory function (CD163/CD163L1, CSF1, CFH, HLA-DQB1, LILRB5, NINJ1 and STAB1). A number of the genes are linked to the mononuclear/phagocyte system and clearance of enzymes from the serum. This highlights the variety in the sources of normal diversity in serum levels of enzymes.

Kidney stone disease is a complex disorder with a strong genetic component. We conducted a genome-wide association study of 28.3 million sequence variants detected through whole-genome sequencing of 2,636 Icelanders that were imputed into 5,419 kidney stone cases, including 2,172 cases with a history of recurrent kidney stones, and 279,870 controls. We identify sequence variants associating with kidney stones at ALPL (rs1256328[T], odds ratio (OR)=1.21, P=5.8 × 10(-10)) and a suggestive association at CASR (rs7627468[A], OR=1.16, P=2.0 × 10(-8)). Focusing our analysis on coding sequence variants in 63 genes with preferential kidney expression we identify two rare missense variants SLC34A1 p.Tyr489Cys (OR=2.38, P=2.8 × 10(-5)) and TRPV5 p.Leu530Arg (OR=3.62, P=4.1 × 10(-5)) associating with recurrent kidney stones. We also observe associations of the identified kidney stone variants with biochemical traits in a large population set, indicating potential biological mechanism.

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.

Chronic kidney disease (CKD) is a complex disorder with a strong genetic component. A number of common sequence variants have been found to associate with serum creatinine (SCr), estimated glomerular filtration rate (eGFR) and/or CKD. We imputed 24 million single-nucleotide polymorphisms and insertions/deletions identified by whole-genome sequencing of 2230 Icelanders into 81 656 chip-typed individuals and 112 630 relatives of genotyped individuals over the age of 18 with SCr measurements. The large set of sequenced individuals allowed accurate imputation of variants to a minor allele frequency (MAF) of 0.1%. We tested the imputed variants for association with SCr. In addition to replicating established loci, we discovered missense and loss-of-function variants associating with SCr in three solute carriers (SLC6A19, SLC25A45 and SLC47A1) and two E3 ubiquitin ligases (RNF186 and RNF128). All the variants are within coding sequences and all but one are rare (MAF <2%) with SCr effects between 0.085 and 0.129 standard deviations. These rare variants have a larger effect on SCr than previously reported common variants, explaining 0.5% of the variability of SCr in Icelanders in addition to the 1% already accounted for. We tested the five variants associating with SCr for association with CKD in an Icelandic sample of 15 594 cases and 291 428 controls. Three of the variants also associated with CKD. These variants may either affect kidney function or creatinine synthesis and excretion. Of note were four mutations in SLC6A19 that associate with reduced SCr, three of which have been shown to cause Hartnup disease., (© The Author 2014. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

To search for sequence variants conferring risk of nonmedullary thyroid cancer, we focused our analysis on 22 SNPs with a P < 5 × 10(-8) in a genome-wide association study on levels of thyroid stimulating hormone (TSH) in 27,758 Icelanders. Of those, rs965513 has previously been shown to associate with thyroid cancer. The remaining 21 SNPs were genotyped in 561 Icelandic individuals with thyroid cancer (cases) and up to 40,013 controls. Variants suggestively associated with thyroid cancer (P < 0.05) were genotyped in an additional 595 non-Icelandic cases and 2,604 controls. After combining the results, three variants were shown to associate with thyroid cancer: rs966423 on 2q35 (OR = 1.34; P(combined) = 1.3 × 10(-9)), rs2439302 on 8p12 (OR = 1.36; P(combined) = 2.0 × 10(-9)) and rs116909374 on 14q13.3 (OR = 2.09; P(combined) = 4.6 × 10(-11)), a region previously reported to contain an uncorrelated variant conferring risk of thyroid cancer. A strong association (P = 9.1 × 10(-91)) was observed between rs2439302 on 8p12 and expression of NRG1, which encodes the signaling protein neuregulin 1, in blood.

We tested 16 million SNPs, identified through whole-genome sequencing of 457 Icelanders, for association with gout and serum uric acid levels. Genotypes were imputed into 41,675 chip-genotyped Icelanders and their relatives, for effective sample sizes of 968 individuals with gout and 15,506 individuals for whom serum uric acid measurements were available. We identified a low-frequency missense variant (c.1580C>G) in ALDH16A1 associated with gout (OR = 3.12, P = 1.5 × 10(-16), at-risk allele frequency = 0.019) and serum uric acid levels (effect = 0.36 s.d., P = 4.5 × 10(-21)). We confirmed the association with gout by performing Sanger sequencing on 6,017 Icelanders. The association with gout was stronger in males relative to females. We also found a second variant on chromosome 1 associated with gout (OR = 1.92, P = 0.046, at-risk allele frequency = 0.986) and serum uric acid levels (effect = 0.48 s.d., P = 4.5 × 10(-16)). This variant is close to a common variant previously associated with serum uric acid levels. This work illustrates how whole-genome sequencing data allow the detection of associations between low-frequency variants and complex traits.

Through complementary application of SNP genotyping, whole-genome sequencing and imputation in 38,384 Icelanders, we have discovered a previously unidentified sick sinus syndrome susceptibility gene, MYH6, encoding the alpha heavy chain subunit of cardiac myosin. A missense variant in this gene, c.2161C>T, results in the conceptual amino acid substitution p.Arg721Trp, has an allelic frequency of 0.38% in Icelanders and associates with sick sinus syndrome with an odds ratio = 12.53 and P = 1.5 × 10⁻²⁹. We show that the lifetime risk of being diagnosed with sick sinus syndrome is around 6% for non-carriers of c.2161C>T but is approximately 50% for carriers of the c.2161C>T variant.

Electrocardiographic measures are indicative of the function of the cardiac conduction system. To search for sequence variants that modulate heart rate, PR interval and QRS duration in individuals of European descent, we performed a genome-wide association study in approximately 10,000 individuals and followed up the top signals in an additional approximately 10,000 individuals. We identified several genome-wide significant associations (with P < 1.6 x 10(-7)). We identified one locus for heart rate (MYH6), four for PR interval (TBX5, SCN10A, CAV1 and ARHGAP24) and four for QRS duration (TBX5, SCN10A, 6p21 and 10q21). We tested for association between these loci and subjects with selected arrhythmias in Icelandic and Norwegian case-control sample sets. We observed correlations between TBX5 and CAV1 and atrial fibrillation (P = 4.0 x 10(-5) and P = 0.00032, respectively), between TBX5 and advanced atrioventricular block (P = 0.0067), and between SCN10A and pacemaker implantation (P = 0.0029). We also replicated previously described associations with the QT interval.

Earlier menarche correlates with shorter adult height and higher childhood body fat. We conducted a genome-wide association study of age at menarche (AAM) on 15,297 Icelandic women. Combined analysis with replication sets from Iceland, Denmark and the Netherlands (N = 10,040) yielded a significant association between rs314280[T] on 6q21, near the LIN28B gene, and AAM (effect = 1.2 months later per allele; P = 1.8 × 10(-14)). A second SNP within the same linkage disequilibrium (LD) block, rs314277, splits rs314280[T] into two haplotypes with different effects (0.9 months and 1.9 months per allele). These variants have been associated with greater adult height. The association with adult height did not account for the association with AAM or vice versa. Other variants, previously associated with height, did not associate significantly with AAM. Given the link between body fat and AAM, we also assessed 11 variants recently associated with higher body mass index (BMI) and 5 of those associated with earlier AAM.