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.

The sequences of 150,119 genomes in the UK Biobank.

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

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.

A partial loss-of-function variant in STAT6 protects against T2 asthma.

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

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

Genome-wide investigation of persistence to treatment with methotrexate in early rheumatoid arthritis.

OBJECTIVES: To investigate the influence of genetic factors on persistence to treatment of early rheumatoid arthritis (RA) with methotrexate (MTX) monotherapy. METHODS: We conducted a genome-wide association study (GWAS) in a sample of 3902 Swedish early RA patients initiating MTX in DMARD-monotherapy as their first ever DMARD. The outcome, short- and long-term persistence to this treatment, was defined as remaining on MTX at one and at three years, respectively, with no additional DMARDs added. As genetic predictors, we investigated individual SNPs, and a polygenic risk score (PRS) based on SNPs associated with RA risk. The SNP-based heritability of persistence was estimated overall and by RA serostatus. RESULTS: No individual SNP reached genome-wide significance (p < 5e-8), neither for persistence at one nor at three years. The RA PRS was not significantly associated with persistence at one (RR = 0.98 (0.96-1.01)) nor three years (RR = 0.96 (0.93-1.00)). The heritability for persistence was estimated to be 0.45 (0.15-0.75) at one year and 0.14 (0-0.40) at three years. Results in seropositive RA were comparable to those in the analysis of RA overall, while heritability estimates and PRS RRs were attenuated towards the null in seronegative RA. CONCLUSIONS: Despite being the largest GWAS on an MTX treatment outcome to date, no genome-wide significant associations were detected. The modest heritability observed, coupled with the broad spread of suggestively associated loci, indicate that genetic influence is of polygenic nature. Nevertheless, persistence to MTX monotherapy was lower in patients with a greater genetic disposition, per the PRS, towards RA.

Short Vi-polysaccharide abrogates T-independent immune response and hyporesponsiveness elicited by long Vi-CRM197 conjugate vaccine.

Polysaccharide-protein conjugates have been developed to overcome the T-independent response, hyporesponsiveness to repeated vaccination, and poor immunogenicity in infants of polysaccharides. To address the impact of polysaccharide length, typhoid conjugates made with short- and long-chain fractions of Vi polysaccharide with average sizes of 9.5, 22.8, 42.7, 82.0, and 165 kDa were compared. Long-chain-conjugated Vi (165 kDa) induced a response in both wild-type and T cell-deficient mice, suggesting that it maintains a T-independent response. In marked contrast, short-chain Vi (9.5 to 42.7 kDa) conjugates induced a response in wild-type mice but not in T cell-deficient mice, suggesting that the response is dependent on T cell help. Mechanistically, this was explained in neonatal mice, in which long-chain, but not short-chain, Vi conjugate induced late apoptosis of Vi-specific B cells in spleen and early depletion of Vi-specific B cells in bone marrow, resulting in hyporesponsiveness and lack of long-term persistence of Vi-specific IgG in serum and IgG+ antibody-secreting cells in bone marrow. We conclude that while conjugation of long-chain Vi generates T-dependent antigens, the conjugates also retain T-independent properties, leading to detrimental effects on immune responses. The data reported here may explain some inconsistencies observed in clinical trials and help guide the design of effective conjugate vaccines.

Polygenic risk score for ACE-inhibitor-associated cough based on the discovery of new genetic loci.

AIMS: To search for sequence variants associated with ACEi discontinuation and to test their association with ACEi-associated adverse drug reactions (ADRs). METHODS AND RESULTS: A genome-wide association study (GWAS) on ACEi discontinuation was conducted, including 33 959 ACEi-discontinuers and 44 041 controls. Cases were defined as persons who switched from an ACEi treatment to an angiotensin receptor blocker. Controls were defined as persons who continued ACEi treatment for at least 1 year. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were computed for ACEi discontinuation risk by mixed model regression analysis. Summary statistics from the individual cohorts were meta-analyzed with a fixed-effects model. To test for association with specific ACEi-associated ADRs, any genome-wide significant (P < 5 × 10-8) ACEi discontinuation variants was tested for association with ACEi-associated cough and angioedema. A polygenetic risk score (PRS) based on ACEi discontinuation GWAS data was constructed and tested for association with ACEi-associated cough and angioedema in two population-based samples. In total, seven genetic genome-wide loci were identified, of which six were previously unreported. The strongest association with ACEi discontinuation was at 20q13.3 (NTSR1; OR: 1.21; 95% CI: 1.17-1.24; P = 2.1 × 10-34). Five of seven lead variants were associated with ACEi-associated cough, whereas none were associated with ACEi-associated angioedema. The ACEi discontinuation PRS was associated with ACEi-associated cough in a dose-response manner but not with ACEi-associated angioedema. ACEi discontinuation was genetically correlated with important causes for cough, including gastro-esophageal reflux disease, allergic rhinitis, hay fever, and asthma, which indicates partly shared genetic underpinning between these traits. CONCLUSION: This study showed the advantage of using prescription patterns to discover genetic links with ADRs. In total, seven genetic loci that associated with ACEi discontinuation were identified. There was evidence of a strong association between our ADR phenotype and ACEi-associated cough. Taken together, these findings increase insight into the pathophysiological processes that underlie ACEi-associated ADRs.

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.

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.

A rare IL33 loss-of-function mutation reduces blood eosinophil counts and protects from asthma.

IL-33 is a tissue-derived cytokine that induces and amplifies eosinophilic inflammation and has emerged as a promising new drug target for asthma and allergic disease. Common variants at IL33 and IL1RL1, encoding the IL-33 receptor ST2, associate with eosinophil counts and asthma. Through whole-genome sequencing and imputation into the Icelandic population, we found a rare variant in IL33 (NM_001199640:exon7:c.487-1G>C (rs146597587-C), allele frequency = 0.65%) that disrupts a canonical splice acceptor site before the last coding exon. It is also found at low frequency in European populations. rs146597587-C associates with lower eosinophil counts (ß = -0.21 SD, P = 2.5×10-16, N = 103,104), and reduced risk of asthma in Europeans (OR = 0.47; 95%CI: 0.32, 0.70, P = 1.8×10-4, N cases = 6,465, N controls = 302,977). Heterozygotes have about 40% lower total IL33 mRNA expression than non-carriers and allele-specific analysis based on RNA sequencing and phased genotypes shows that only 20% of the total expression is from the mutated chromosome. In half of those transcripts the mutation causes retention of the last intron, predicted to result in a premature stop codon that leads to truncation of 66 amino acids. The truncated IL-33 has normal intracellular localization but neither binds IL-33R/ST2 nor activates ST2-expressing cells. Together these data demonstrate that rs146597587-C is a loss of function mutation and support the hypothesis that IL-33 haploinsufficiency protects against asthma.

Variant ASGR1 Associated with a Reduced Risk of Coronary Artery Disease.

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

Clonal hematopoiesis, with and without candidate driver mutations, is common in the elderly.

Clonal hematopoiesis (CH) arises when a substantial proportion of mature blood cells is derived from a single dominant hematopoietic stem cell lineage. Somatic mutations in candidate driver (CD) genes are thought to be responsible for at least some cases of CH. Using whole-genome sequencing of 11 262 Icelanders, we found 1403 cases of CH by using barcodes of mosaic somatic mutations in peripheral blood, whether or not they have a mutation in a CD gene. We find that CH is very common in the elderly, trending toward inevitability. We show that somatic mutations in TET2, DNMT3A, ASXL1, and PPM1D are associated with CH at high significance. However, known CD mutations were evident in only a fraction of CH cases. Nevertheless, the highly prevalent CH we detect associates with increased mortality rates, risk for hematological malignancy, smoking behavior, telomere length, Y-chromosome loss, and other phenotypic characteristics. Modeling suggests some CH cases could arise in the absence of CD mutations as a result of neutral drift acting on a small population of active hematopoietic stem cells. Finally, we find a germline deletion in intron 3 of the telomerase reverse transcriptase (TERT) gene that predisposes to CH (rs34002450; P = 7.4 × 10-12; odds ratio, 1.37).

A rare missense mutation in MYH6 associates with non-syndromic coarctation of the aorta.

Aims: Coarctation of the aorta (CoA) accounts for 4-8% of congenital heart defects (CHDs) and confers substantial morbidity despite treatment. It is increasingly recognized as a highly heritable condition. The aim of the study was to search for sequence variants that affect the risk of CoA. Methods and results: We performed a genome-wide association study of CoA among Icelanders (120 cases and 355 166 controls) based on imputed variants identified through whole-genome sequencing. We found association with a rare (frequency = 0.34%) missense mutation p.Arg721Trp in MYH6 (odds ratio = 44.2, P = 5.0 × 10-22), encoding the alpha-heavy chain subunit of cardiac myosin, an essential sarcomere protein. Approximately 20% of individuals with CoA in Iceland carry this mutation. We show that p.Arg721Trp also associates with other CHDs, in particular bicuspid aortic valve. We have previously reported broad effects of p.Arg721Trp on cardiac electrical function and strong association with sick sinus syndrome and atrial fibrillation. Conclusion: Through a population approach, we found that a rare missense mutation p.Arg721Trp in the sarcomere gene MYH6 has a strong effect on the risk of CoA and explains a substantial fraction of the Icelanders with CoA. This is the first mutation associated with non-familial or sporadic form of CoA at a population level. The p.Arg721Trp in MYH6 causes a cardiac syndrome with highly variable expressivity and emphasizes the importance of sarcomere integrity for cardiac development and function.

Causal Effect of Genetic Variants Associated With Body Mass Index on Multiple Sclerosis Susceptibility.

Multiple sclerosis (MS) is an autoimmune disease with both genetic and environmental risk factors. Recent studies indicate that childhood and adolescent obesity double the risk of MS, but this association may reflect unmeasured confounders rather than causal effects of obesity. We used separate-sample Mendelian randomization to estimate the causal effect of body mass index (BMI) on susceptibility to MS. Using data from non-Hispanic white members of the Kaiser Permanente Medical Care Plan of Northern California (KPNC) (2006-2014; 1,104 cases of MS and 10,536 controls) and a replication data set from Sweden (the Epidemiological Investigation of MS (EIMS) and the Genes and Environment in MS (GEMS) studies, 2005-2013; 5,133 MS cases and 4,718 controls), we constructed a weighted genetic risk score using 97 variants previously established to predict BMI. Results were adjusted for birth year, sex, education, smoking status, ancestry, and genetic predictors of MS. Estimates in KPNC and Swedish data sets suggested that higher genetically induced BMI predicted greater susceptibility to MS (odds ratio = 1.13, 95% confidence interval: 1.04, 1.22 for the KPNC sample; odds ratio = 1.09, 95% confidence interval: 1.03, 1.15 for the Swedish sample). Although the mechanism remains unclear, to our knowledge, these findings support a causal effect of increased BMI on susceptibility to MS for the first time, and they suggest a role for inflammatory pathways that characterize both obesity and the MS disease process.