Complex effects of sequence variants on lipid levels and coronary artery disease.

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

Large-scale plasma proteomics comparisons through genetics and disease associations.

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 Project1 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 people2, 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.

PSMD11 loss-of-function variants correlate with a neurobehavioral phenotype, obesity, and increased interferon response.

Primary proteasomopathies have recently emerged as a new class of rare early-onset neurodevelopmental disorders (NDDs) caused by pathogenic variants in the PSMB1, PSMC1, PSMC3, or PSMD12 proteasome genes. Proteasomes are large multi-subunit protein complexes that maintain cellular protein homeostasis by clearing ubiquitin-tagged damaged, misfolded, or unnecessary proteins. In this study, we have identified PSMD11 as an additional proteasome gene in which pathogenic variation is associated with an NDD-causing proteasomopathy. PSMD11 loss-of-function variants caused early-onset syndromic intellectual disability and neurodevelopmental delay with recurrent obesity in 10 unrelated children. Our findings demonstrate that the cognitive impairment observed in these individuals could be recapitulated in Drosophila melanogaster with depletion of the PMSD11 ortholog Rpn6, which compromised reversal learning. Our investigations in subject samples further revealed that PSMD11 loss of function resulted in impaired 26S proteasome assembly and the acquisition of a persistent type I interferon (IFN) gene signature, mediated by the integrated stress response (ISR) protein kinase R (PKR). In summary, these data identify PSMD11 as an additional member of the growing family of genes associated with neurodevelopmental proteasomopathies and provide insights into proteasomal biology in human health.

Genetic insights into biological mechanisms governing human ovarian ageing.

Reproductive longevity is essential for fertility and influences healthy ageing in women1,2, but insights into its underlying biological mechanisms and treatments to preserve it are limited. Here we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in about 200,000 women of European ancestry. These common alleles were associated with clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations3. The identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR-associated genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increases fertility and extends reproductive life in mice. Causal inference analyses using the identified genetic variants indicate that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases the risk of hormone-sensitive cancers. These findings provide insight into the mechanisms that govern ovarian ageing, when they act, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease.

Actionable Genotypes and Their Association with Life Span in Iceland.

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

FLT3 stop mutation increases FLT3 ligand level and risk of autoimmune thyroid disease.

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

Genetic predisposition to mosaic Y chromosome loss in blood.

Mosaic loss of chromosome Y (LOY) in circulating white blood cells is the most common form of clonal mosaicism1-5, yet our knowledge of the causes and consequences of this is limited. Here, using a computational approach, we estimate that 20% of the male population represented in the UK Biobank study (n = 205,011) has detectable LOY. We identify 156 autosomal genetic determinants of LOY, which we replicate in 757,114 men of European and Japanese ancestry. These loci highlight genes that are involved in cell-cycle regulation and cancer susceptibility, as well as somatic drivers of tumour growth and targets of cancer therapy. We demonstrate that genetic susceptibility to LOY is associated with non-haematological effects on health in both men and women, which supports the hypothesis that clonal haematopoiesis is a biomarker of genomic instability in other tissues. Single-cell RNA sequencing identifies dysregulated expression of autosomal genes in leukocytes with LOY and provides insights into why clonal expansion of these cells may occur. Collectively, these data highlight the value of studying clonal mosaicism to uncover fundamental mechanisms that underlie cancer and other ageing-related diseases.

Genetic variants associated with syncope implicate neural and autonomic processes.

AIMS: Syncope is a common and clinically challenging condition. In this study, the genetics of syncope were investigated to seek knowledge about its pathophysiology and prognostic implications. METHODS AND RESULTS: This genome-wide association meta-analysis included 56 071 syncope cases and 890 790 controls from deCODE genetics (Iceland), UK Biobank (United Kingdom), and Copenhagen Hospital Biobank Cardiovascular Study/Danish Blood Donor Study (Denmark), with a follow-up assessment of variants in 22 412 cases and 286 003 controls from Intermountain (Utah, USA) and FinnGen (Finland). The study yielded 18 independent syncope variants, 17 of which were novel. One of the variants, p.Ser140Thr in PTPRN2, affected syncope only when maternally inherited. Another variant associated with a vasovagal reaction during blood donation and five others with heart rate and/or blood pressure regulation, with variable directions of effects. None of the 18 associations could be attributed to cardiovascular or other disorders. Annotation with regard to regulatory elements indicated that the syncope variants were preferentially located in neural-specific regulatory regions. Mendelian randomization analysis supported a causal effect of coronary artery disease on syncope. A polygenic score (PGS) for syncope captured genetic correlation with cardiovascular disorders, diabetes, depression, and shortened lifespan. However, a score based solely on the 18 syncope variants performed similarly to the PGS in detecting syncope risk but did not associate with other disorders. CONCLUSION: The results demonstrate that syncope has a distinct genetic architecture that implicates neural regulatory processes and a complex relationship with heart rate and blood pressure regulation. A shared genetic background with poor cardiovascular health was observed, supporting the importance of a thorough assessment of individuals presenting with syncope.

Spread of SARS-CoV-2 in the Icelandic Population.

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

Parental influence on human germline de novo mutations in 1,548 trios from Iceland.

The characterization of mutational processes that generate sequence diversity in the human genome is of paramount importance both to medical genetics and to evolutionary studies. To understand how the age and sex of transmitting parents affect de novo mutations, here we sequence 1,548 Icelanders, their parents, and, for a subset of 225, at least one child, to 35× genome-wide coverage. We find 108,778 de novo mutations, both single nucleotide polymorphisms and indels, and determine the parent of origin of 42,961. The number of de novo mutations from mothers increases by 0.37 per year of age (95% CI 0.32-0.43), a quarter of the 1.51 per year from fathers (95% CI 1.45-1.57). The number of clustered mutations increases faster with the mother's age than with the father's, and the genomic span of maternal de novo mutation clusters is greater than that of paternal ones. The types of de novo mutation from mothers change substantially with age, with a 0.26% (95% CI 0.19-0.33%) decrease in cytosine-phosphate-guanine to thymine-phosphate-guanine (CpG>TpG) de novo mutations and a 0.33% (95% CI 0.28-0.38%) increase in C>G de novo mutations per year, respectively. Remarkably, these age-related changes are not distributed uniformly across the genome. A striking example is a 20 megabase region on chromosome 8p, with a maternal C>G mutation rate that is up to 50-fold greater than the rest of the genome. The age-related accumulation of maternal non-crossover gene conversions also mostly occurs within these regions. Increased sequence diversity and linkage disequilibrium of C>G variants within regions affected by excess maternal mutations indicate that the underlying mutational process has persisted in humans for thousands of years. Moreover, the regional excess of C>G variation in humans is largely shared by chimpanzees, less by gorillas, and is almost absent from orangutans. This demonstrates that sequence diversity in humans results from evolving interactions between age, sex, mutation type, and genomic location.

Evaluation of Large-Scale Proteomics for Prediction of Cardiovascular Events.

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.

Multiomics analysis of rheumatoid arthritis yields sequence variants that have large effects on risk of the seropositive subset.

OBJECTIVES: To find causal genes for rheumatoid arthritis (RA) and its seropositive (RF and/or ACPA positive) and seronegative subsets. METHODS: We performed a genome-wide association study (GWAS) of 31 313 RA cases (68% seropositive) and ~1 million controls from Northwestern Europe. We searched for causal genes outside the HLA-locus through effect on coding, mRNA expression in several tissues and/or levels of plasma proteins (SomaScan) and did network analysis (Qiagen). RESULTS: We found 25 sequence variants for RA overall, 33 for seropositive and 2 for seronegative RA, altogether 37 sequence variants at 34 non-HLA loci, of which 15 are novel. Genomic, transcriptomic and proteomic analysis of these yielded 25 causal genes in seropositive RA and additional two overall. Most encode proteins in the network of interferon-alpha/beta and IL-12/23 that signal through the JAK/STAT-pathway. Highlighting those with largest effect on seropositive RA, a rare missense variant in STAT4 (rs140675301-A) that is independent of reported non-coding STAT4-variants, increases the risk of seropositive RA 2.27-fold (p=2.1×10-9), more than the rs2476601-A missense variant in PTPN22 (OR=1.59, p=1.3×10-160). STAT4 rs140675301-A replaces hydrophilic glutamic acid with hydrophobic valine (Glu128Val) in a conserved, surface-exposed loop. A stop-mutation (rs76428106-C) in FLT3 increases seropositive RA risk (OR=1.35, p=6.6×10-11). Independent missense variants in TYK2 (rs34536443-C, rs12720356-C, rs35018800-A, latter two novel) associate with decreased risk of seropositive RA (ORs=0.63-0.87, p=10-9-10-27) and decreased plasma levels of interferon-alpha/beta receptor 1 that signals through TYK2/JAK1/STAT4. CONCLUSION: Sequence variants pointing to causal genes in the JAK/STAT pathway have largest effect on seropositive RA, while associations with seronegative RA remain scarce.

Large-Scale Screening for Monogenic and Clinically Defined Familial Hypercholesterolemia in Iceland.

Objective: Familial hypercholesterolemia (FH) is traditionally defined as a monogenic disease characterized by severely elevated LDL-C (low-density lipoprotein cholesterol) levels. In practice, FH is commonly a clinical diagnosis without confirmation of a causative mutation. In this study, we sought to characterize and compare monogenic and clinically defined FH in a large sample of Icelanders. Approach and Results: We whole-genome sequenced 49 962 Icelanders and imputed the identified variants into an overall sample of 166 281 chip-genotyped Icelanders. We identified 20 FH mutations in LDLR, APOB, and PCSK9 with combined prevalence of 1 in 836. Monogenic FH was associated with severely elevated LDL-C levels and increased risk of premature coronary disease, aortic valve stenosis, and high burden of coronary atherosclerosis. We used a modified version of the Dutch Lipid Clinic Network criteria to screen for the clinical FH phenotype among living adult participants (N=79 058). Clinical FH was found in 2.2% of participants, of whom only 5.2% had monogenic FH. Mutation-negative clinical FH has a strong polygenic basis. Both individuals with monogenic FH and individuals with mutation-negative clinical FH were markedly undertreated with cholesterol-lowering medications and only a minority attained an LDL-C target of <2.6 mmol/L (<100 mg/dL; 11.0% and 24.9%, respectively) or <1.8 mmol/L (<70 mg/dL; 0.0% and 5.2%, respectively), as recommended for primary prevention by European Society of Cardiology/European Atherosclerosis Society cholesterol guidelines. Conclusions: Clinically defined FH is a relatively common phenotype that is explained by monogenic FH in only a minority of cases. Both monogenic and clinical FH confer high cardiovascular risk but are markedly undertreated.

Genetic insight into sick sinus syndrome.

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

Genetic insight into sick sinus syndrome.

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

Sequence variants associating with urinary biomarkers.

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.

Genetic variability in the absorption of dietary sterols affects the risk of coronary artery disease.

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.