Genetics and epidemiology of mutational barcode-defined clonal hematopoiesis.

Clonal hematopoiesis (CH) arises when a substantial proportion of mature blood cells is derived from a single hematopoietic stem cell lineage. Using whole-genome sequencing of 45,510 Icelandic and 130,709 UK Biobank participants combined with a mutational barcode method, we identified 16,306 people with CH. Prevalence approaches 50% in elderly participants. Smoking demonstrates a dosage-dependent impact on risk of CH. CH associates with several smoking-related diseases. Contrary to published claims, we find no evidence that CH is associated with cardiovascular disease. We provide evidence that CH is driven by genes that are commonly mutated in myeloid neoplasia and implicate several new driver genes. The presence and nature of a driver mutation alters the risk profile for hematological disorders. Nevertheless, most CH cases have no known driver mutations. A CH genome-wide association study identified 25 loci, including 19 not implicated previously in CH. Splicing, protein and expression quantitative trait loci were identified for CD164 and TCL1A.

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.

Genome-wide meta-analysis identifies 93 risk loci and enables risk prediction equivalent to monogenic forms of venous thromboembolism.

We report a genome-wide association study of venous thromboembolism (VTE) incorporating 81,190 cases and 1,419,671 controls sampled from six cohorts. We identify 93 risk loci, of which 62 are previously unreported. Many of the identified risk loci are at genes encoding proteins with functions converging on the coagulation cascade or platelet function. A VTE polygenic risk score (PRS) enabled effective identification of both high- and low-risk individuals. Individuals within the top 0.1% of PRS distribution had a VTE risk similar to homozygous or compound heterozygous carriers of the variants G20210A (c.*97 G > A) in F2 and p.R534Q in F5. We also document that F2 and F5 mutation carriers in the bottom 10% of the PRS distribution had a risk similar to that of the general population. We further show that PRS improved individual risk prediction beyond that of genetic and clinical risk factors. We investigated the extent to which venous and arterial thrombosis share clinical risk factors using Mendelian randomization, finding that some risk factors for arterial thrombosis were directionally concordant with VTE risk (for example, body mass index and smoking) whereas others were discordant (for example, systolic blood pressure and triglyceride levels).

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 predisposition to hypertension is associated with preeclampsia in European and Central Asian women.

Preeclampsia is a serious complication of pregnancy, affecting both maternal and fetal health. In genome-wide association meta-analysis of European and Central Asian mothers, we identify sequence variants that associate with preeclampsia in the maternal genome at ZNF831/20q13 and FTO/16q12. These are previously established variants for blood pressure (BP) and the FTO variant has also been associated with body mass index (BMI). Further analysis of BP variants establishes that variants at MECOM/3q26, FGF5/4q21 and SH2B3/12q24 also associate with preeclampsia through the maternal genome. We further show that a polygenic risk score for hypertension associates with preeclampsia. However, comparison with gestational hypertension indicates that additional factors modify the risk of preeclampsia.

Abdominal aortic aneurysm is associated with a variant in low-density lipoprotein receptor-related protein 1.

Abdominal aortic aneurysm (AAA) is a common cause of morbidity and mortality and has a significant heritability. We carried out a genome-wide association discovery study of 1866 patients with AAA and 5435 controls and replication of promising signals (lead SNP with a p value < 1 × 10(-5)) in 2871 additional cases and 32,687 controls and performed further follow-up in 1491 AAA and 11,060 controls. In the discovery study, nine loci demonstrated association with AAA (p < 1 × 10(-5)). In the replication sample, the lead SNP at one of these loci, rs1466535, located within intron 1 of low-density-lipoprotein receptor-related protein 1 (LRP1) demonstrated significant association (p = 0.0042). We confirmed the association of rs1466535 and AAA in our follow-up study (p = 0.035). In a combined analysis (6228 AAA and 49182 controls), rs1466535 had a consistent effect size and direction in all sample sets (combined p = 4.52 × 10(-10), odds ratio 1.15 [1.10-1.21]). No associations were seen for either rs1466535 or the 12q13.3 locus in independent association studies of coronary artery disease, blood pressure, diabetes, or hyperlipidaemia, suggesting that this locus is specific to AAA. Gene-expression studies demonstrated a trend toward increased LRP1 expression for the rs1466535 CC genotype in arterial tissues; there was a significant (p = 0.029) 1.19-fold (1.04-1.36) increase in LRP1 expression in CC homozygotes compared to TT homozygotes in aortic adventitia. Functional studies demonstrated that rs1466535 might alter a SREBP-1 binding site and influence enhancer activity at the locus. In conclusion, this study has identified a biologically plausible genetic variant associated specifically with AAA, and we suggest that this variant has a possible functional role in LRP1 expression.

Genome-wide association study identifies a sequence variant within the DAB2IP gene conferring susceptibility to abdominal aortic aneurysm.

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 x 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 x 10(-5)), peripheral arterial disease (OR = 1.14, P = 3.9 x 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.

Sequence variants at CHRNB3-CHRNA6 and CYP2A6 affect smoking behavior.

Smoking is a common risk factor for many diseases. We conducted genome-wide association meta-analyses for the number of cigarettes smoked per day (CPD) in smokers (n = 31,266) and smoking initiation (n = 46,481) using samples from the ENGAGE Consortium. In a second stage, we tested selected SNPs with in silico replication in the Tobacco and Genetics (TAG) and Glaxo Smith Kline (Ox-GSK) consortia cohorts (n = 45,691 smokers) and assessed some of those in a third sample of European ancestry (n = 9,040). Variants in three genomic regions associated with CPD (P < 5 x 10(-8)), including previously identified SNPs at 15q25 represented by rs1051730[A] (effect size = 0.80 CPD, P = 2.4 x 10(-69)), and SNPs at 19q13 and 8p11, represented by rs4105144[C] (effect size = 0.39 CPD, P = 2.2 x 10(-12)) and rs6474412-T (effect size = 0.29 CPD, P = 1.4 x 10(-8)), respectively. Among the genes at the two newly associated loci are genes encoding nicotine-metabolizing enzymes (CYP2A6 and CYP2B6) and nicotinic acetylcholine receptor subunits (CHRNB3 and CHRNA6), all of which have been highlighted in previous studies of smoking and nicotine dependence. Nominal associations with lung cancer were observed at both 8p11 (rs6474412[T], odds ratio (OR) = 1.09, P = 0.04) and 19q13 (rs4105144[C], OR = 1.12, P = 0.0006).

New common variants affecting susceptibility to basal cell carcinoma.

In a follow-up to our previously reported genome-wide association study of cutaneous basal cell carcinoma (BCC), we describe here several new susceptibility variants. SNP rs11170164, encoding a G138E substitution in the keratin 5 (KRT5) gene, affects risk of BCC (OR = 1.35, P = 2.1 x 10(-9)). A variant at 9p21 near CDKN2A and CDKN2B also confers susceptibility to BCC (rs2151280[C]; OR = 1.19, P = 6.9 x 10(-9)), as does rs157935[T] at 7q32 near the imprinted gene KLF14 (OR = 1.23, P = 5.7 x 10(-10)). The effect of rs157935[T] is dependent on the parental origin of the risk allele. None of these variants were found to be associated with melanoma or fair-pigmentation traits. A melanoma- and pigmentation-associated variant in the SLC45A2 gene, L374F, is associated with risk of both BCC and squamous cell carcinoma. Finally, we report conclusive evidence that rs401681[C] in the TERT-CLPTM1L locus confers susceptibility to BCC but protects against melanoma.

A drastic reduction in the life span of cystatin C L68Q carriers due to life-style changes during the last two centuries.

Hereditary cystatin C amyloid angiopathy (HCCAA) is an autosomal dominant disease with high penetrance, manifest by brain hemorrhages in young normotensive adults. In Iceland, this condition is caused by the L68Q mutation in the cystatin C gene, with contemporary carriers reaching an average age of only 30 years. Here, we report, based both on linkage disequilibrium and genealogical evidence, that all known copies of this mutation derive from a common ancestor born roughly 18 generations ago. Intriguingly, the genealogies reveal that obligate L68Q carriers born 1825 to 1900 experienced a drastic reduction in life span, from 65 years to the present-day average. At the same time, a parent-of-origin effect emerged, whereby maternal inheritance of the mutation was associated with a 9 year reduction in life span relative to paternal inheritance. As these trends can be observed in several different extended families, many generations after the mutational event, it seems likely that some environmental factor is responsible, perhaps linked to radical changes in the life-style of Icelanders during this period. A mutation with such radically different phenotypic effects in reaction to normal variation in human life-style not only opens the possibility of preventive strategies for HCCAA, but it may also provide novel insights into the complex relationship between genotype and environment in human disease.

Multiple genetic loci for bone mineral density and fractures.

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

ASIP and TYR pigmentation variants associate with cutaneous melanoma and basal cell carcinoma.

Fair color increases risk of cutaneous melanoma (CM) and basal cell carcinoma (BCC). Recent genome-wide association studies have identified variants affecting hair, eye and skin pigmentation in Europeans. Here, we assess the effect of these variants on risk of CM and BCC in European populations comprising 2,121 individuals with CM, 2,163 individuals with BCC and over 40,000 controls. A haplotype near ASIP, known to affect a similar spectrum of pigmentation traits as MC1R variants, conferred significant risk of CM (odds ratio (OR) = 1.45, P = 1.2 x 10(-9)) and BCC (OR = 1.33, P = 1.2 x 10(-6)). The variant in TYR encoding the R402Q amino acid substitution, previously shown to affect eye color and tanning response, conferred risk of CM (OR = 1.21, P = 2.8 x 10(-7)) and BCC (OR = 1.14, P = 6.1 x 10(-4)). An eye color variant in TYRP1 was associated with risk of CM (OR = 1.15, P = 4.6 x 10(-4)). The association of all three variants is robust with respect to adjustment for the effect of pigmentation.

A common variant on chromosome 9p21 affects the risk of myocardial infarction.

The global endemic of cardiovascular diseases calls for improved risk assessment and treatment. Here, we describe an association between myocardial infarction (MI) and a common sequence variant on chromosome 9p21. This study included a total of 4587 cases and 12,767 controls. The identified variant, adjacent to the tumor suppressor genes CDKN2A and CDKN2B, was associated with the disease with high significance. Approximately 21% of individuals in the population are homozygous for this variant, and their estimated risk of suffering myocardial infarction is 1.64 times as great as that of noncarriers. The corresponding risk is 2.02 times as great for early-onset cases. The population attributable risk is 21% for MI in general and 31% for early-onset cases.