Meta-analysis of ACE inhibitor-induced angioedema identifies novel risk locus

BACKGROUND

Angioedema is a rare but potentially life-threatening adverse drug reaction in patients receiving angiotensin-converting enzyme inhibitors (ACEis). Research suggests that susceptibility to ACEi-induced angioedema (ACEi-AE) involves both genetic and nongenetic risk factors. Genome- and exome-wide studies of ACEi-AE have identified the first genetic risk loci. However, understanding of the underlying pathophysiology remains limited.

OBJECTIVE

We sought to identify further genetic factors of ACEi-AE to eventually gain a deeper understanding of its pathophysiology.

METHODS

By combining data from 8 cohorts, a genome-wide association study meta-analysis was performed in more than 1000 European patients with ACEi-AE. Secondary bioinformatic analyses were conducted to fine-map associated loci, identify relevant genes and pathways, and assess the genetic overlap between ACEi-AE and other traits. Finally, an exploratory cross-ancestry analysis was performed to assess shared genetic factors in European and African-American patients with ACEi-AE.

RESULTS

Three genome-wide significant risk loci were identified. One of these, located on chromosome 20q11.22, has not been implicated previously in ACEi-AE. Integrative secondary analyses highlighted previously reported genes (BDKRB2 [bradykinin receptor B2] and F5 [coagulation factor 5]) as well as biologically plausible novel candidate genes (PROCR [protein C receptor] and EDEM2 [endoplasmic reticulum degradation enhancing alpha-mannosidase like protein 2]). Lead variants at the risk loci were found with similar effect sizes and directions in an African-American cohort.

CONCLUSIONS

The present results contributed to a deeper understanding of the pathophysiology of ACEi-AE by (1) providing further evidence for the involvement of bradykinin signaling and coagulation pathways and (2) suggesting, for the first time, the involvement of the fibrinolysis pathway in this adverse drug reaction. An exploratory cross-ancestry comparison implicated the relevance of the associated risk loci across diverse ancestries.

Whole genome case-control study of central nervous system toxicity due to antimicrobial drugs

A genetic predisposition to central nervous system (CNS) toxicity induced by antimicrobial drugs (antibiotics, antivirals, antifungals, and antiparasitic drugs) has been suspected. Whole genome sequencing of 66 cases and 833 controls was performed to investigate whether antimicrobial drug-induced CNS toxicity was associated with genetic variation. The primary objective was to test whether antimicrobial-induced CNS toxicity was associated with seventeen efflux transporters at the blood-brain barrier. In this study, variants or structural elements in efflux transporters were not significantly associated with CNS toxicity. Secondary objectives were to test whether antimicrobial-induced CNS toxicity was associated with genes over the whole genome, with HLA, or with structural genetic variation. Uncommon variants in and close to three genes were significantly associated with CNS toxicity according to a sequence kernel association test combined with an optimal unified test (SKAT-O). These genes were LCP1 (q = 0.013), RETSAT (q = 0.013) and SFMBT2 (q = 0.035). Two variants were driving the LCP1 association: rs6561297 (p = 1.15x10-6, OR: 4.60 [95% CI: 2.51-8.46]) and the regulatory variant rs10492451 (p = 1.15x10-6, OR: 4.60 [95% CI: 2.51-8.46]). No common genetic variant, HLA-type or structural variation was associated with CNS toxicity. In conclusion, CNS toxicity due to antimicrobial drugs was associated with uncommon variants in LCP1, RETSAT and SFMBT2.

No link between type I interferon autoantibody positivity and adverse reactions to COVID-19 vaccines

Type I interferons act as gatekeepers against viral infection, and autoantibodies that neutralize these signaling molecules have been associated with COVID-19 severity and adverse reactions to the live-attenuated yellow fever vaccine. On this background, we sought to examine whether autoantibodies against type I interferons were associated with adverse events following COVID-19 vaccination. Our nationwide analysis suggests that type I interferon autoantibodies were not associated with adverse events after mRNA or viral-vector COVID-19 vaccines.

Identification of risk factors for adverse drug reactions in a pharmacovigilance database

INTRODUCTION

In addition to identifying new safety signals, pharmacovigilance databases could be used to identify potential risk factors for adverse drug reactions (ADRs).

OBJECTIVE

To evaluate whether data mining in a pharmacovigilance database can be used to identify known and possible novel risk factors for ADRs, for use in pharmacovigilance practice.

METHOD

Exploratory data mining was performed within the Swedish national database of spontaneously reported ADRs. Bleeding associated with direct oral anticoagulants (DOACs)-rivaroxaban, apixaban, edoxaban, and dabigatran-was used as a test model. We compared demographics, drug treatment, and clinical features between cases with bleeding (N = 965) and controls who had experienced other serious ADRs to DOACs (N = 511). Statistical analysis was performed by unadjusted and age adjusted logistic regression models, and the random forest based machine-learning method Boruta.

RESULTS

In the logistic regression, 13 factors were significantly more common among cases of bleeding compared with controls. Eleven were labelled or previously proposed risk factors. Cardiac arrhythmia (e.g., atrial fibrillation), hypertension, mental impairment disorders (e.g., dementia), renal and urinary tract procedures, gastrointestinal ulceration and perforation, and interacting drugs remained significant after adjustment for age. In the Boruta analysis, high age, arrhythmia, hypertension, cardiac failure, thromboembolism, and pharmacodynamically interacting drugs had a larger than random association with the outcome. High age, cardiac arrhythmia, hypertension, cardiac failure, and pharmacodynamically interacting drugs had odds ratios for bleeding above one, while thromboembolism had an odds ratio below one.

CONCLUSIONS

We demonstrated that data mining within a pharmacovigilance database identifies known risk factors for DOAC bleeding, and potential risk factors such as dementia and atrial fibrillation. We propose that the method could be used in pharmacovigilance for identification of potential ADR risk factors that merit further evaluation.

Narcolepsy risk loci outline role of T cell autoimmunity and infectious triggers in narcolepsy

Narcolepsy type 1 (NT1) is caused by a loss of hypocretin/orexin transmission. Risk factors include pandemic 2009 H1N1 influenza A infection and immunization with Pandemrix®. Here, we dissect disease mechanisms and interactions with environmental triggers in a multi-ethnic sample of 6,073 cases and 84,856 controls. We fine-mapped GWAS signals within HLA (DQ0602, DQB1*03:01 and DPB1*04:02) and discovered seven novel associations (CD207, NAB1, IKZF4-ERBB3, CTSC, DENND1B, SIRPG, PRF1). Significant signals at TRA and DQB1*06:02 loci were found in 245 vaccination-related cases, who also shared polygenic risk. T cell receptor associations in NT1 modulated TRAJ*24, TRAJ*28 and TRBV*4-2 chain-usage. Partitioned heritability and immune cell enrichment analyses found genetic signals to be driven by dendritic and helper T cells. Lastly comorbidity analysis using data from FinnGen, suggests shared effects between NT1 and other autoimmune diseases. NT1 genetic variants shape autoimmunity and response to environmental triggers, including influenza A infection and immunization with Pandemrix®.

Genetic determinants of apixaban plasma levels and their relationship to bleeding and thromboembolic events

Apixaban is a direct oral anticoagulant, a factor Xa inhibitor, used for the prevention of ischemic stroke in patients with atrial fibrillation. Despite using recommended dosing a few patients might still experience bleeding or lack of efficacy that might be related to inappropriate drug exposure. We conducted a genome-wide association study using data from 1,325 participants in the pivotal phase three trial of apixaban with the aim to identify genetic factors affecting the pharmacokinetics of apixaban. A candidate gene analysis was also performed for pre-specified variants in ABCB1, ABCG2, CYP3A4, CYP3A5, and SULT1A1, with a subsequent analysis of all available polymorphisms within the candidate genes. Significant findings were further evaluated to assess a potential association with clinical outcome such as bleeding or thromboembolic events. No variant was consistently associated with an altered apixaban exposure on a genome-wide level. The candidate gene analyses showed a statistically significant association with a well-known variant in the drug transporter gene ABCG2 (c.421G > T, rs2231142). Patients carrying this variant had a higher exposure to apixaban [area under the curve (AUC), beta = 151 (95% CI 59–243), p = 0.001]. On average, heterozygotes displayed a 5% increase of AUC and homozygotes a 17% increase of AUC, compared with homozygotes for the wild-type allele. Bleeding or thromboembolic events were not significantly associated with ABCG2 rs2231142. This large genome-wide study demonstrates that genetic variation in the drug transporter gene ABCG2 is associated with the pharmacokinetics of apixaban. However, the influence of this finding on drug exposure was small, and further studies are needed to better understand whether it is of relevance for ischemic and bleeding events.

Genome-wide association study of liver enzyme elevation in an extended cohort of rheumatoid arthritis patients starting low-dose methotrexate

Aim: A follow-up genome-wide association study (GWAS) in an extended cohort of rheumatoid arthritis (RA) patients starting low-dose methotrexate (MTX) treatment was performed to identify further genetic variants associated with alanine aminotransferase (ALT) elevation. Patients & methods: A GWAS was performed on 346 RA patients. Two outcomes within the first 6 months of MTX treatment were assessed: ALT >1.5-times the upper level of normal (ULN) and maximum level of ALT. Results: SPATA9 (rs72783407) was significantly associated with maximum level of ALT (p = 2.58 × 10^-8) and PLCG2 (rs60427389) was tentatively associated with ALT >1.5 × ULN. Conclusion: Associations with SNPs in genes related to male fertility (SPATA9) and inflammatory processes (PLCG2) were identified.

Molecular Genetic Screening in Patients With ACE Inhibitor/Angiotensin Receptor Blocker-Induced Angioedema to Explore the Role of Hereditary Angioedema Genes

Angioedema is a relatively rare but potentially life-threatening adverse reaction to angiotensin-converting enzyme inhibitors (ACEi) and angiotensin receptor blockers (ARBs). As with hereditary forms of angioedema (HAE), this adverse reaction is mediated by bradykinin. Research suggests that ACEi/ARB-induced angioedema has a multifactorial etiology. In addition, recent case reports suggest that some ACEi/ARB-induced angioedema patients may carry pathogenic HAE variants. The aim of the present study was to investigate the possible association between ACEi/ARB-induced angioedema and HAE genes via systematic molecular genetic screening in a large cohort of ACEi/ARB-induced angioedema cases. Targeted re-sequencing of five HAE-associated genes (SERPING1, F12, PLG, ANGPT1, and KNG1) was performed in 212 ACEi/ARB-induced angioedema patients recruited in Germany/Austria, Sweden, and Denmark, and in 352 controls from a German cohort. Among patients, none of the identified variants represented a known pathogenic variant for HAE. Moreover, no significant association with ACEi/ARB-induced angioedema was found for any of the identified common [minor allele frequency (MAF) >5%] or rare (MAF < 5%) variants. However, several non-significant trends suggestive of possible protective effects were observed. The lowest p-value for an individual variant was found in PLG (rs4252129, p.R523W, p = 0.057, p.adjust > 0.999, Fisher’s exact test). Variant p.R523W was found exclusively in controls and has previously been associated with decreased levels of plasminogen, a precursor of plasmin which is part of a pathway directly involved in bradykinin production. In addition, rare, potentially functional variants (MAF < 5%, Phred-scaled combined annotation dependent depletion score >10) showed a nominally significant enrichment in controls both: 1) across all five genes; and 2) in the F12 gene alone. However, these results did not withstand correction for multiple testing. In conclusion, our results suggest that HAE-associated mutations are, at best, a rare cause of ACEi/ARB-induced angioedema. Furthermore, we were unable to identify a significant association between ACEi/ARB-induced angioedema and other variants in the investigated genes. Further studies with larger sample sizes are warranted to draw more definite conclusions concerning variants with limited effect sizes, including protective variants.

The Clinical Pharmacogenetics Implementation Consortium (CPIC) guideline for SLCO1B1, ABCG2 , and CYP2C9 and statin‐associated musculoskeletal symptoms

Statins reduce cholesterol, prevent cardiovascular disease, and are among the most commonly prescribed medications in the world. Statin-associated musculoskeletal symptoms (SAMS) impact statin adherence and ultimately can impede the long-term effectiveness of statin therapy. There are several identified pharmacogenetic variants that impact statin disposition and adverse events during statin therapy. SLCO1B1 encodes a transporter (SLCO1B1; alternative names include OATP1B1 or OATP-C) that facilitates the hepatic uptake of all statins. ABCG2 encodes an efflux transporter (BCRP) that modulates the absorption and disposition of rosuvastatin. CYP2C9 encodes a phase I drug metabolizing enzyme responsible for the oxidation of some statins. Genetic variation in each of these genes alters systemic exposure to statins (i.e., simvastatin, rosuvastatin, pravastatin, pitavastatin, atorvastatin, fluvastatin, lovastatin), which can increase the risk for SAMS. We summarize the literature supporting these associations and provide therapeutic recommendations for statins based on SLCO1B1, ABCG2, and CYP2C9 genotype with the goal of improving the overall safety, adherence, and effectiveness of statin therapy. This document replaces the 2012 and 2014 Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines for SLCO1B1 and simvastatin-induced myopathy.

HLA-variants associated with azathioprine-induced pancreatitis in patients with Crohn's disease

The immunosuppressant drug azathioprine is associated with a 4% risk of acute pancreatitis in patients with inflammatory bowel disease (IBD). Studies have demonstrated an increased risk in carriers of HLA‐DQA1*02:01 and HLA‐DRB1*07:01. We investigated whether these human leukocyte antigen (HLA) types were associated with azathioprine‐induced pancreatitis also in Swedish patients with IBD, and whether the type of disease affected the association. Nineteen individuals with IBD who developed acute pancreatitis after initiation of azathioprine were genotyped and compared with a population control cohort (n = 4891) and a control group matched for disease (n = 81). HLA‐DQA1*02:01 and HLA‐DRB1*07:01 were in full linkage disequilibrium, and were significantly associated with acute pancreatitis both when cases were compared with population controls (OR 3.97 [95% CI 1.57–9.97], p = 0.0035) and matched controls (OR 3.55 [95% CI 1.23–10.98], p = 0.0275). In a disease‐specific analysis, the correlation was positive in patients with Crohn's disease versus matched controls (OR 9.27 [95% CI 1.86–46.19], p = 0.0066), but not in those with ulcerative colitis versus matched controls (OR 0.69 [95% CI 0.07–6.74], p = 0.749). In patients with Crohn's disease, we estimated the conditional risk of carriers of HLA‐DQA1*02:01‐HLA‐DRB1*07:01 to 7.3%, and the conditional risk of a non‐carrier to 2.2%. We conclude that HLA‐DQA1*02:01‐HLA‐DRB1*07:01 is a marker for increased risk of acute pancreatitis in individuals of Swedish genetic origin, treated with azathioprine for Crohn's disease.

Genome-wide association study of liver enzyme elevation in rheumatoid arthritis patients starting methotrexate

Aim: To identify novel genetic variants predisposing to elevation of Alanine aminotransferase (ALT) in rheumatoid arthritis (RA) patients after initiation of methotrexate (MTX) treatment. Patients & methods: We performed genome-wide association studies in 198 RA patients starting MTX. Outcomes were maximum level of ALT and ALT >1.5-times the upper level of normal within the first 6 months of treatment. Results: RAVER2 (rs72675408) was significantly associated with maximum level of ALT (p = 4.36 × 10^-8). This variant is in linkage disequilibrium with rs72675451, which is associated with differential expression of JAK1 and RAVER2. Conclusion: We found an association between ALT elevation and genetic variants that may regulate the expression of JAK1 and RAVER2. JAK1 encodes a janus kinase involved in the pathogenesis of RA.

Technological readiness and implementation of genomic-driven precision medicine for complex diseases

The fields of human genetics and genomics have generated considerable knowledge about the mechanistic basis of many diseases. Genomic approaches to diagnosis, prognostication, prevention and treatment - genomic-driven precision medicine (GDPM) - may help optimize medical practice. Here, we provide a comprehensive review of GDPM of complex diseases across major medical specialties. We focus on technological readiness: how rapidly a test can be implemented into health care. Although these areas of medicine are diverse, key similarities exist across almost all areas. Many medical areas have, within their standards of care, at least one GDPM test for a genetic variant of strong effect that aids the identification/diagnosis of a more homogeneous subset within a larger disease group or identifies a subset with different therapeutic requirements. However, for almost all complex diseases, the majority of patients do not carry established single-gene mutations with large effects. Thus, research is underway that seeks to determine the polygenic basis of many complex diseases. Nevertheless, most complex diseases are caused by the interplay of genetic, behavioural and environmental risk factors, which will likely necessitate models for prediction and diagnosis that incorporate genetic and non-genetic data.

[New recommendation on pharmacogenetic screening prior to 5-fluorouracil based cancer treatment - Swedish experience indicates less adverse effects and healthcare cost savings]

5-fluorouracil (5-FU) is still a cornerstone in drug treatment for cancer. Some patients starting standard dosed 5-FU will experience severe adverse events (SAEs). One mechanism behind SAEs is impaired dihydropyrimidine dehydrogenase (DPD) activity, resulting in an accumulation of cytotoxic metabolites. Pre-emptive testing of DPD enzyme activity or genetic variation in its gene, DPYD,  is recommended since 2020 in Sweden. We report experience from DPYD testing in 368 patients planned for 5-FU treatment. DPYD variants associated with reduced DPD activity were observed in 28 patients (8%), which is close to the expected frequency. These patients tolerated 5-FU treatment when doses were reduced according to guidelines. However, 4 out of 5 variant allele carriers starting 5-FU at standard dose due to late arrival of test results experienced SAEs. Pre-emptive testing was calculated to be cost saving and thus beneficial from a healthcare economy perspective.

PharmVar GeneFocus: CYP2C9

The Pharmacogene Variation Consortium (PharmVar) catalogues star (*) allele nomenclature for the polymorphic human CYP2C9 gene. Genetic variation within the CYP2C9 gene locus impacts the metabolism or bioactivation of many clinically important drugs, including nonsteroidal anti-inflammatory drugs, phenytoin, antidiabetic agents, and angiotensin receptor blockers. Variable CYP2C9 activity is of particular importance regarding efficacy and safety of warfarin and siponimod as indicated in their package inserts. This GeneFocus provides a comprehensive overview and summary of CYP2C9 and describes how haplotype information catalogued by PharmVar is utilized by the Pharmacogenomics Knowledgebase and the Clinical Pharmacogenetics Implementation Consortium.

[Pharmacogenomics - a cornerstone of Precision Medicine. Genomic Medicine Sweden analyses genotypes associated with serious drug toxicity or therapeutic failure]

Serious adverse drug reactions, drug intolerance, and lack of effect are major problems in healthcare. Pharmacogenomics is the part of precision medicine that aims to develop predictive risk markers in this respect and establish such testing in clinical practice. The nation-wide project Genomic Medicine Sweden (GMS) is undertaking large-scale sequencing to predict risk of drug toxicity and lack of efficacy in malignant diseases. The aim is to facilitate an improved, individualized treatment with increased patient safety. In addition to accurate genotyping, other technical or infrastructure-related aspects need to be considered for a successful implementation in healthcare, for example electronic accessibility and visibility of pharmacogenomic data of long-standing relevance for an individual's ongoing and future drug treatment.

Genetic Risk Factors in Drug-Induced Liver Injury Due to Isoniazid-Containing Antituberculosis Drug Regimens

Drug-induced liver injury (DILI) is a complication of treatment with antituberculosis (TB) drugs, especially in isoniazid (INH)-containing regimens. To investigate genetic risk factors, we performed a genomewide association study (GWAS) involving anti-TB DILI cases (55 Indian and 70 European) and controls (1,199 Indian and 10,397 European). Most cases were treated with a standard anti-TB drug regimen; all received INH. We imputed single nucleotide polymorphism and HLA genotypes and performed trans-ethnic meta-analysis on GWAS and candidate gene genotypes. GWAS found one significant association (rs117491755) in Europeans only. For HLA, HLA-B*52:01 was significant (meta-analysis odds ratio (OR) 2.67, 95% confidence interval (CI) 1.63-4.37, P = 9.4 × 10^-5 ). For N-acetyltransferase 2 (NAT2), NAT2*5 frequency was lower in cases (OR 0.69, 95% CI 0.57-0.83, P = 0.01). NAT2*6 and NAT2*7 were more common, with homozygotes for NAT2*6 and/or NAT2*7 enriched among cases (OR 1.89, 95% CI 0.84-4.22, P = 0.004). We conclude HLA genotype makes a small contribution to TB drug-related DILI and that the NAT2 contribution is complex, but consistent with previous reports when differences in the metabolic effect of NAT2*5 compared with those of NAT2*6 and NAT2*7 are considered.

Exome Sequencing Reveals Common and Rare Variants in F5 Associated With ACE Inhibitor and Angiotensin Receptor Blocker-Induced Angioedema

Angioedema occurring in the head and neck region is a rare and sometimes life-threatening adverse reaction to angiotensin-converting enzyme inhibitors (ACEIs) and angiotensin receptor blockers (ARBs). Few studies have investigated the association of common variants with this extreme reaction, but none have explored the combined influence of rare variants yet. Adjudicated cases of ACEI-induced angioedema (ACEI-AE) or ARB-induced angioedema (ARB-AE) and controls were recruited at five different centers. Sequencing of 1,066 samples (408 ACEI-AE, ARB-AE, and 658 controls) was performed using exome-enriched sequence data. A common variant of the F5 gene that causes an increase in blood clotting (rs6025, p.Arg506Gln, also called factor V Leiden), was significantly associated with both ACEI-AE and ARB-AE (odds ratio: 2.85, 95% confidence interval (CI), 1.89-4.25). A burden test analysis of five rare missense variants in F5 was also found to be associated with ACEI-AE or ARB-AE, P = 2.09 × 10-3 . A combined gene risk score of these variants, and the common variants rs6025 and rs6020, showed that individuals carrying at least one variant had 2.21 (95% CI, 1.49-3.27, P = 6.30 × 10-9 ) times the odds of having ACEI-AE or ARB-AE. The increased risk due to the common Leiden allele was confirmed in a genome-wide association study from the United States. A high risk of angioedema was also observed for the rs6020 variant that is the main coagulation defect-causing variant in black African and Asian populations. We found that deleterious missense variants in F5 are associated with an increased risk of ACEI-AE or ARB-AE.

Genome-Wide Association Study of Metamizole-Induced Agranulocytosis in European Populations

Agranulocytosis is a rare yet severe idiosyncratic adverse drug reaction to metamizole, an analgesic widely used in countries such as Switzerland and Germany. Notably, an underlying mechanism has not yet been fully elucidated and no predictive factors are known to identify at-risk patients. With the aim to identify genetic susceptibility variants to metamizole-induced agranulocytosis (MIA) and neutropenia (MIN), we conducted a retrospective multi-center collaboration including cases and controls from three European populations. Association analyses were performed using genome-wide genotyping data from a Swiss cohort (45 cases, 191 controls) followed by replication in two independent European cohorts (41 cases, 273 controls) and a joint discovery meta-analysis. No genome-wide significant associations (p < 1 × 10⁻⁷) were observed in the Swiss cohort or in the joint meta-analysis, and no candidate genes suggesting an immune-mediated mechanism were identified. In the joint meta-analysis of MIA cases across all cohorts, two candidate loci on chromosome 9 were identified, rs55898176 (OR = 4.01, 95%CI: 2.41–6.68, p = 1.01 × 10⁻⁷) and rs4427239 (OR = 5.47, 95%CI: 2.81–10.65, p = 5.75 × 10⁻⁷), of which the latter is located in the SVEP1 gene previously implicated in hematopoiesis. This first genome-wide association study for MIA identified suggestive associations with biological plausibility that may be used as a stepping-stone for post-GWAS analyses to gain further insight into the mechanism underlying MIA.

High-Throughput Sequencing to Investigate Associations Between HLA Genes and Metamizole-Induced Agranulocytosis

Background and Objective: Agranulocytosis is a rare and potentially life-threatening complication of metamizole (dipyrone) intake that is characterized by a loss of circulating neutrophil granulocytes. While the mechanism underlying this adverse drug reaction is not well understood, involvement of the immune system has been suggested. In addition, associations between genetic variants in the Human Leukocyte Antigen (HLA) region and agranulocytosis induced by other drugs have been reported. The aim of the present study was to assess whether genetic variants in classical HLA genes are associated with the susceptibility to metamizole-induced agranulocytosis (MIA) in a European population by targeted resequencing of eight HLA genes.

Design: A case-control cohort of Swiss patients with a history of neutropenia or agranulocytosis associated with metamizole exposure (n = 53), metamizole-tolerant (n = 39) and unexposed controls (n = 161) was recruited for this study. A high-throughput resequencing (HTS) and high-resolution typing method was used to sequence and analyze eight HLA loci in a discovery subset of this cohort (n = 31 cases, n = 38 controls). Identified candidate alleles were investigated in the full Swiss cohort as well as in two independent cohorts from Germany and Spain using HLA imputation from genome-wide SNP array data. In addition, variant calling based on HTS data was performed in the discovery subset for the class I genes HLA-A, -B, and -C using the HLA-specific mapper hla-mapper.

Results: Eight candidate alleles (p < 0.05) were identified in the discovery subset, of which HLA-C^∗04:01 was associated with MIA in the full Swiss cohort (p < 0.01) restricted to agranulocytosis (ANC < 0.5 × 10⁹/L) cases. However, no candidate allele showed a consistent association in the Swiss, German and Spanish cohorts. Analysis of individual sequence variants in class I genes produced consistent results with HLA typing but did not reveal additional small nucleotide variants associated with MIA.

Conclusion: Our results do not support an HLA-restricted T cell-mediated immune mechanism for MIA. However, we established an efficient high-resolution (three-field) eight-locus HTS HLA resequencing method to interrogate the HLA region and demonstrated the feasibility of its application to pharmacogenetic studies.

MTHFR, TYMS and SLCO1B1 polymorphisms and adverse liver effects of methotrexate in rheumatoid arthritis

Aims: To investigate whether variants of MTHFR, TYMS and SLCO1B1 are associated with ALT elevation in rheumatoid arthritis patients starting methotrexate (MTX). Patients & methods: Clinical and laboratory data were collected from the start of MTX treatment. Genotyping of MTHFR, TYMS and SLCO1B1 was performed. Univariate and multiple logistic regression were used for statistical analysis. Results: 34 out of 369 patients experienced ALT >1.5 × ULN less than 6 months from start. MTHFR A1298C (rs1801131) was nominally associated with an ALT >1.5 × ULN within 6 months after the start of MTX (OR = 1.7 [95% CI: 1.04-2.9]; p = 0.03), but did not pass correction for multiple testing. A multiple model containing MTHFR 1298C and clinical factors predicted the outcome (C-statistic 0.735). TYMS and SLCO1B1 were not associated with the outcome. Conclusions: A model containing MTHFR 1298C and clinical factors might predict risk of early ALT elevation.

SWEDEGENE-a Swedish nation-wide DNA sample collection for pharmacogenomic studies of serious adverse drug reactions

SWEDEGENE is a Swedish nation-wide sample collection established to facilitate studies of clinical and genetic risk factors for adverse drug reactions (ADRs). Most cases are recruited among patients reported to the ADR registry at the Swedish Medical Products Agency by health-care professionals. Clinical data are collected both from medical and laboratory records and through interviews using standardized questionnaires. Genome-wide scans and whole-genome sequencing are done, and association studies are conducted using mainly controls from the Swedish TwinGene biobank with data on diagnoses and prescribed drugs. SWEDEGENE was established in 2008 and currently contains DNA and information from about 2550 adults who have experienced specific ADRs, and from 580 drug exposed controls. Results from genome-wide association studies have now been published, and data from whole-genome sequencing are being analyzed. SWEDEGENE has the potential to offer a new means of developing individualized and safe drug therapy through patient pre-treatment screening.

Shared Genetic Risk Factors Across Carbamazepine-Induced Hypersensitivity Reactions

Carbamazepine (CBZ) causes life‐threating T‐cell‐mediated hypersensitivity reactions, including serious cutaneous adverse reactions (SCARs) and drug‐induced liver injury (CBZ‐DILI). In order to evaluate shared or phenotype‐specific genetic predisposing factors for CBZ hypersensitivity reactions, we performed a meta‐analysis of two genomewide association studies (GWAS) on a total of 43 well‐phenotyped Northern and Southern European CBZ‐SCAR cases and 10,701 population controls and a GWAS on 12 CBZ‐DILI cases and 8,438 ethnically matched population controls. HLA‐A*31:01 was identified as the strongest genetic predisposing factor for both CBZ‐SCAR (odds ratio (OR) = 8.0; 95% CI 4.10–15.80; P = 1.2 × 10⁻⁹) and CBZ‐DILI (OR = 7.3; 95% CI 2.47–23.67; P = 0.0004) in European populations. The association with HLA‐A*31:01 in patients with SCAR was mainly driven by hypersensitivity syndrome (OR = 12.9; P = 2.1 × 10⁻⁹) rather than by Stevens‐Johnson syndrome/toxic epidermal necrolysis cases, which showed an association with HLA‐B*57:01. We also identified a novel risk locus mapping to ALK only for CBZ‐SCAR cases, which needs replication in additional cohorts and functional evaluation.

A Genome-Wide Association Study of Bisphosphonate-Associated Atypical Femoral Fracture

Atypical femoral fracture is a well-documented adverse reaction to bisphosphonates. It is strongly related to duration of bisphosphonate use, and the risk declines rapidly after drug withdrawal. The mechanism behind bisphosphonate-associated atypical femoral fracture is unclear, but a genetic predisposition has been suggested. With the aim to identify common genetic variants that could be used for preemptive genetic testing, we performed a genome-wide association study. Cases were recruited mainly through reports of adverse drug reactions sent to the Swedish Medical Products Agency on a nation-wide basis. We compared atypical femoral fracture cases (n = 51) with population-based controls (n = 4891), and to reduce the possibility of confounding by indication, we also compared with bisphosphonate-treated controls without a current diagnosis of cancer (n = 324). The total number of single-nucleotide polymorphisms after imputation was 7,585,874. A genome-wide significance threshold of p < 5 × 10^-8 was used to correct for multiple testing. In addition, we performed candidate gene analyses for a panel of 29 genes previously implicated in atypical femoral fractures (significance threshold of p < 5.7 × 10^-6). Compared with population controls, bisphosphonate-associated atypical femoral fracture was associated with four isolated, uncommon single-nucleotide polymorphisms. When cases were compared with bisphosphonate-treated controls, no statistically significant genome-wide association remained. We conclude that the detected associations were either false positives or related to the underlying disease, i.e., treatment indication. Furthermore, there was no significant association with single-nucleotide polymorphisms in the 29 candidate genes. In conclusion, this study found no evidence of a common genetic predisposition for bisphosphonate-associated atypical femoral fracture. Further studies of larger sample size to identify possible weakly associated genetic traits, as well as whole exome or whole-genome sequencing studies to identify possible rare genetic variation conferring a risk are warranted.

Pharmacogenomics of statin-related myopathy: Meta-analysis of rare variants from whole-exome sequencing

AIMS

Statin-related myopathy (SRM), which includes rhabdomyolysis, is an uncommon but important adverse drug reaction because the number of people prescribed statins world-wide is large. Previous association studies of common genetic variants have had limited success in identifying a genetic basis for this adverse drug reaction. We conducted a multi-site whole-exome sequencing study to investigate whether rare coding variants confer an increased risk of SRM.

METHODS AND RESULTS

SRM 3-5 cases (N = 505) and statin treatment-tolerant controls (N = 2047) were recruited from multiple sites in North America and Europe. SRM 3-5 was defined as symptoms consistent with muscle injury and an elevated creatine phosphokinase level >4 times upper limit of normal without another likely cause of muscle injury. Whole-exome sequencing and variant calling was coordinated from two analysis centres, and results of single-variant and gene-based burden tests were meta-analysed. No genome-wide significant associations were identified. Given the large number of cases, we had 80% power to identify a variant with minor allele frequency of 0.01 that increases the risk of SRM 6-fold at genome-wide significance.

CONCLUSIONS

In this large whole-exome sequencing study of severe statin-related muscle injury conducted to date, we did not find evidence that rare coding variants are responsible for this adverse drug reaction. Larger sample sizes would be required to identify rare variants with small effects, but it is unclear whether such findings would be clinically actionable.

Methotrexate treatment in rheumatoid arthritis and elevated liver enzymes: A long-term follow-up of predictors, surveillance, and outcome in clinical practice

Aim

To assess predictors of alanine aminotransferase (ALT) elevation in methotrexate (MTX) treated rheumatoid arthritis (RA) patients, and to describe the monitoring of liver enzymes, including handling and outcome of elevated ALT.

Methods

All RA patients starting MTX in January, 2005 to April, 2013 at a rheumatology clinic, (Uppsala University Hospital, Sweden) were identified from electronic medical records. Clinical and laboratory data were obtained from medical records, supplemented by telephone interviews. Predictors for ALT >1.5× over the upper limit of normal (ULN) were identified by multiple regression analysis.

Results

The study comprised 213 RA patients starting MTX. During a mean follow‐up of 4.3 years, 6288 ALT tests were performed; 7% of tests with ALT were >ULN. ALT >1.5× ULN was observed in 44 (21%) patients and the strongest predictor was a pre‐treatment elevation of ALT (adjusted odds ratio = 6.8, 95% CI 2.2‐20.5). Recurrent elevations occurred in 70% of patients who continued treatment, and the proportion was similar in those with and without interventions, for example MTX dose reduction (67% vs 73%, P = 0.43). Seven patients (3%) permanently stopped MTX due to ALT elevation, and two were eventually diagnosed with non‐alcoholic fatty liver disease. No patient developed hepatic failure.

Conclusion

Only a small number of ALT tests performed during MTX therapy in RA capture an elevation. A pre‐treatment elevation of ALT was the strongest predictor for early and recurrent ALT elevations during therapy. This study supports a more individualized approach to monitoring and handling of ALT elevations during MTX therapy in RA than recommended in current guidelines.

Fusidic Acid: A Neglected Risk Factor for Statin-Associated Myopathy

Background:

Statins are widely used lipid-lowering drugs used for the prevention of cardiovascular disease. Statins are known to cause myopathy, an adverse drug reaction with various clinical features rhabdomyolysis.

Objective:

To describe clinical characteristics of statin-treated individuals who experienced myopathy and identify risk factors of statin-associated myopathy.

Methods:

A retrospective study was conducted on cases of statin-associated myopathy reported to the Swedish Medical Products Agency. Clinical factors were compared between cases and statin-treated controls not diagnosed with myopathy. Statistical methods were univariate and multivariate logistic regression and results were presented as odds ratio (OR) with 95% confidence interval (CI). To correct for multiple comparisons, the cutoff for statistical significance was set to P < .0017.

Results:

In total, 47 cases of statin-associated myopathy were compared with 3871 treated controls. Rhabdomyolysis was diagnosed in 51% of the cases. Markers for cardiovascular disease were more common in cases than controls. Statistical analysis revealed the following independent risk factors for myopathy: high statin dose (OR = 1.54, calculated using the standard deviation 19.82, 95% CI = 1.32-1.80, P < .0001), and concomitant treatment with fusidic acid (OR = 1002, 95% CI = 54.55-18 410, P < .0001), cyclosporine (OR = 34.10, 95% CI = 4.43-262.45, P = .0007), and gemfibrozil (OR = 12.35, 95% CI = 2.38-64.10, P = .0028).

Conclusions:

The risk of myopathy increases with statin dose and cotreatment with cyclosporine and gemfibrozil. Concomitant fusidic acid has previously only been noted in a few case reports. Considering that use of fusidic acid may become more frequent, it is important to remind of this risk factor for statin-associated myopathy.

Clinical factors predicting drug-induced liver injury due to flucloxacillin

Objectives

Drug-induced liver injury (DILI) is a serious adverse reaction due to flucloxacillin. The pathogenesis is not fully understood. Female sex, age over 60 years, and a longer treatment duration have been suggested to be predisposing factors. Carriers of HLA-B*57:01 have an 80-fold increased risk, but due to the rarity of the reaction, testing of all patients is not cost-effective. We aimed to validate and detect clinical risk factors for flucloxacillin DILI.

Methods

Clinical characteristics of flucloxacillin-treated patients with (n=50) and without DILI (n=2,330) were compared in a retrospective case control study. Cases were recruited from the Swedish database of spontaneously reported adverse drug reactions. Treated controls were selected from the Swedish Twin Registry. Statistical comparisons were made using chi-squared test and logistic regression. The significance threshold was set to P<0.00357 to correct for multiple comparisons. Reliable variables were tested in a multiple regression model.

Results

DILI was associated with female sex, OR 2.79, 95% CI 1.50-5.17, P=0.0011, and with a history of kidney stones, OR 5.51, 95% CI 2.21-13.72, P=0.0003. Cases were younger than controls, OR per increase in years 0.91, 95% CI 0.88-0.94, P<0.0001, probably due to selection bias. No difference in treatment duration was detected, OR 1.03, 95% CI 0.98-1.08, P=0.1790.

Conclusion

We established female sex as a risk factor for flucloxacillin-induced DILI, and a history of kidney stones was identified as a potential risk factor. Clinical risk factors for flucloxacillin-induced DILI could be used to indicate whom to test for HLA-B*57:01 before treatment.

CKM Glu83Gly Is Associated With Blunted Creatine Kinase Variation, but Not With Myalgia

BACKGROUND

To test the association of a recently reported variant in the creatine kinase (CK) muscle gene, CKM Glu83Gly (rs11559024) with constitutive creatine phosphokinase (CK) levels, CK variation, and inducibility. Given the diagnostic importance of CK in determining muscle damage, we tested the association of the variant with myalgia.

METHODS AND RESULTS

Meta-analysis between longitudinal cohort GoDARTS (Genetics of Diabetes Audit and Research, Tayside Scotland), minor allele frequency (=0.02), and randomized clinical trial (JUPITER [Justification for the Use of Statins in Primary Prevention: An Intervention Trial Evaluating Rosuvastatin], minor allele frequency=0.018) was used to replicate the association with baseline CK measures. GoDARTS was used to study the relationship with CK variability. Myalgia was studied in JUPITER trial participants. Baseline and SDs of CK were on average 18% (P value=6×10^-63) and 24% (P value=2×10^-5) lower for carriers of the variant, respectively. The variant was not associated with myalgia (odds ratio, 0.84; 95% confidence interval, 0.52-1.38).

CONCLUSIONS

This study highlights that a genetic factor known to be associated with constitutive CK levels is also associated with CK variability and inducibility. This is discussed in the context of evidence to suggest that the variant has an impact on inducibility of CK by trauma through a previously reported case of a homozygous carrier. However, the lack of association between the variant and myalgia suggests that it cannot reliably be used as a biomarker for muscle symptoms.

A common missense variant of LILRB5 is associated with statin intolerance and myalgia

Aims

A genetic variant in LILRB5 (leukocyte immunoglobulin-like receptor subfamily-B) (rs12975366: T > C: Asp247Gly) has been reported to be associated with lower creatine phosphokinase (CK) and lactate dehydrogenase (LDH) levels. Both biomarkers are released from injured muscle tissue, making this variant a potential candidate for susceptibility to muscle-related symptoms. We examined the association of this variant with statin intolerance ascertained from electronic medical records in the GoDARTS study.

Methods and results

In the GoDARTS cohort, the LILRB5 Asp247 variant was associated with statin intolerance (SI) phenotypes; one defined as having raised CK and being non-adherent to therapy [odds ratio (OR) 1.81; 95% confidence interval (CI): 1.34-2.45] and the other as being intolerant to the lowest approved dose of a statin before being switched to two or more other statins (OR 1.36; 95% CI: 1.07-1.73). Those homozygous for Asp247 had increased odds of developing both definitions of intolerance. Importantly the second definition did not rely on CK elevations. These results were replicated in adjudicated cases of statin-induced myopathy in the PREDICTION-ADR consortium (OR1.48; 95% CI: 1.05-2.10) and for the development of myalgia in the JUPITER randomized clinical trial of rosuvastatin (OR1.35, 95% CI: 1.10-1.68). A meta-analysis across the studies showed a consistent association between Asp247Gly and outcomes associated with SI (OR1.34; 95% CI: 1.16-1.54).

Conclusion

This study presents a novel immunogenetic factor associated with statin intolerance, an important risk factor for cardiovascular outcomes. The results suggest that true statin-induced myalgia and non-specific myalgia are distinct, with a potential role for the immune system in their development. We identify a genetic group that is more likely to be intolerant to their statins.

Cohort-specific imputation of gene expression improves prediction of warfarin dose for African Americans

Background

Genome-wide association studies are useful for discovering genotype–phenotype associations but are limited because they require large cohorts to identify a signal, which can be population-specific. Mapping genetic variation to genes improves power and allows the effects of both protein-coding variation as well as variation in expression to be combined into “gene level” effects.

Methods

Previous work has shown that warfarin dose can be predicted using information from genetic variation that affects protein-coding regions.

Here, we introduce a method that improves dose prediction by integrating tissue-specific gene expression. In particular, we use drug pathways and expression quantitative trait loci knowledge to impute gene expression—on the assumption that differential expression of key pathway genes may impact dose requirement. We focus on 116 genes from the pharmacokinetic and pharmacodynamic pathways of warfarin within training and validation sets comprising both European and African-descent individuals.

Results

We build gene-tissue signatures associated with warfarin dose in a cohort-specific manner and identify a signature of 11 gene-tissue pairs that significantly augments the International Warfarin Pharmacogenetics Consortium dosage-prediction algorithm in both populations.

Conclusions

Our results demonstrate that imputed expression can improve dose prediction and bridge population-specific compositions.

MATLAB code is available at https://github.com/assafgo/warfarin-cohort

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-017-0495-0) contains supplementary material, which is available to authorized users.

Sulfasalazine-Induced Agranulocytosis Is Associated With the Human Leukocyte Antigen Locus

Agranulocytosis is a serious, although rare, adverse reaction to sulfasalazine, which is used to treat inflammatory joint and bowel disease. We performed a genome-wide association study comprising 9,380,034 polymorphisms and 180 HLA alleles in 36 cases of sulfasalazine-induced agranulocytosis and 5,170 population controls. Sulfasalazine-induced agranulocytosis was significantly associated with the HLA region on chromosome 6. The top hit (rs9266634) was located close to HLA-B, odds ratio (OR) 5.36 (95% confidence interval (CI) (2.97, 9.69) P = 2.55 × 10^-8 ). We HLA-sequenced a second cohort consisting of 40 cases and 142 treated controls, and confirmed significant associations with HLA-B*08:01, OR = 2.25 (95% CI (1.02, 4.97) P = 0.0439), in particular the HLA-B*08:01 haplotype HLA-DQB1*02:01-DRB1*03:01-B*08:01-C*07:01, OR = 3.79 (95% CI (1.63, 8.80) P = 0.0019), and with HLA-A*31:01, OR = 4.81 (95% CI (1.52, 15.26) P = 0.0077). The number needed to test for HLA-B*08:01 and HLA-A*31:01 to avoid one case was estimated to be 1,500. We suggest that intensified monitoring or alternative treatment should be considered for known carriers of HLA-B*08:01 or HLA-A*31:01.

Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for Pharmacogenetics-Guided Warfarin Dosing: 2017 Update

This document is an update to the 2011 Clinical Pharmacogenetics Implementation Consortium (CPIC) guideline for CYP2C9 and VKORC1 genotypes and warfarin dosing. Evidence from the published literature is presented for CYP2C9, VKORC1, CYP4F2, and rs12777823 genotype-guided warfarin dosing to achieve a target international normalized ratio of 2-3 when clinical genotype results are available. In addition, this updated guideline incorporates recommendations for adult and pediatric patients that are specific to continental ancestry.

Comparison of Clinical Factors Between Patients With Angiotensin-Converting Enzyme Inhibitor-Induced Angioedema and Cough

BACKGROUND

Angioedema is a rare and serious adverse drug reaction (ADR) to angiotensin-converting enzyme (ACE) inhibitor treatment. Dry cough is a common side effect of ACE inhibitors and has been identified as a possible risk factor for angioedema.

OBJECTIVE

We compared characteristics between patients with ACE inhibitor-induced angioedema and cough with the aim of identifying risk factors that differ between these adverse events.

METHODS

Data on patients with angioedema or cough induced by ACE inhibitors were collected from the Swedish database of spontaneously reported ADRs or from collaborating clinicians. Wilcoxon rank sum test, Fisher's exact test, and odds ratios (ORs) with 95% CIs were used to test for between-group differences. The significance threshold was set to P <0.00128 to correct for multiple comparisons.

RESULTS

Clinical characteristics were compared between 168 patients with angioedema and 121 with cough only. Smoking and concomitant selective calcium channel blocker treatment were more frequent among patients with angioedema than cough: OR = 4.3, 95% CI = 2.1-8.9, P = 2.2 × 10^-5, and OR = 3.7, 95% CI = 2.0-7.0, P = 1.7 × 10^-5. Angioedema cases were seen more often in male patients (OR = 2.2, 95% CI = 1.4-3.6, P = 1.3 × 10^-4) and had longer time to onset and higher doses than those with cough ( P = 3.2 × 10^-10 and P = 2.6 × 10^-4). A multiple model containing the variables smoking, concurrent calcium channel blocker treatment, male sex, and time to onset accounted for 26% of the variance between the groups.

CONCLUSION

Smoking, comedication with selective calcium channel blockers, male sex, and longer treatment time were associated with ACE inhibitor-induced angioedema rather than cough.

Association of Liver Injury From Specific Drugs, or Groups of Drugs, With Polymorphisms in HLA and Other Genes in a Genome-Wide Association Study

BACKGROUND & AIMS

We performed a genome-wide association study (GWAS) to identify genetic risk factors for drug-induced liver injury (DILI) from licensed drugs without previously reported genetic risk factors.

METHODS

We performed a GWAS of 862 persons with DILI and 10,588 population-matched controls. The first set of cases was recruited before May 2009 in Europe (n = 137) and the United States (n = 274). The second set of cases were identified from May 2009 through May 2013 from international collaborative studies performed in Europe, the United States, and South America. For the GWAS, we included only cases with patients of European ancestry associated with a particular drug (but not flucloxacillin or amoxicillin-clavulanate). We used DNA samples from all subjects to analyze HLA genes and single nucleotide polymorphisms. After the discovery analysis was concluded, we validated our findings using data from 283 European patients with diagnosis of DILI associated with various drugs.

RESULTS

We associated DILI with rs114577328 (a proxy for A*33:01 a HLA class I allele; odds ratio [OR], 2.7; 95% confidence interval [CI], 1.9-3.8; P = 2.4 × 10^-8) and with rs72631567 on chromosome 2 (OR, 2.0; 95% CI, 1.6-2.5; P = 9.7 × 10^-9). The association with A*33:01 was mediated by large effects for terbinafine-, fenofibrate-, and ticlopidine-related DILI. The variant on chromosome 2 was associated with DILI from a variety of drugs. Further phenotypic analysis indicated that the association between DILI and A*33:01 was significant genome wide for cholestatic and mixed DILI, but not for hepatocellular DILI; the polymorphism on chromosome 2 was associated with cholestatic and mixed DILI as well as hepatocellular DILI. We identified an association between rs28521457 (within the lipopolysaccharide-responsive vesicle trafficking, beach and anchor containing gene) and only hepatocellular DILI (OR, 2.1; 95% CI, 1.6-2.7; P = 4.8 × 10^-9). We did not associate any specific drug classes with genetic polymorphisms, except for statin-associated DILI, which was associated with rs116561224 on chromosome 18 (OR, 5.4; 95% CI, 3.0-9.5; P = 7.1 × 10^-9). We validated the association between A*33:01 terbinafine- and sertraline-induced DILI. We could not validate the association between DILI and rs72631567, rs28521457, or rs116561224.

CONCLUSIONS

In a GWAS of persons of European descent with DILI, we associated HLA-A*33:01 with DILI due to terbinafine and possibly fenofibrate and ticlopidine. We identified polymorphisms that appear to be associated with DILI from statins, as well as 2 non-drug-specific risk factors.

Time to onset of bisphosphonate-related osteonecrosis of the jaws: a multicentre retrospective cohort study

OBJECTIVES

Osteonecrosis of the jaw (ONJ) is a potentially severe adverse effect of bisphosphonates (BP). Although the risk of ONJ increases with increasing duration of BP treatment, there are currently no reliable estimates of the ONJ time to onset (TTO). The objective of this study was to estimate the TTO and associated risk factors in BP-treated patients.

SUBJECTS AND METHODS

Retrospective analysis of data from 22 secondary care centres in seven countries relevant to 349 patients who developed BP-related ONJ between 2004 and 2012.

RESULTS

The median (95%CI) TTO was 6.0 years in patients treated with alendronate (n = 88) and 2.2 years in those treated with zoledronate (n = 218). Multivariable Cox regression showed that dentoalveolar surgery was inversely associated, and the use of antiangiogenics directly associated, with the TTO in patients with cancer treated with zoledronate.

CONCLUSIONS

The incidence of ONJ increases with the duration of BP therapy, with notable differences observed with respect to BP type and potency, route of administration and underlying disease. When data are stratified by BP type, a time of 6.0 and 2.2 years of oral alendronate and intravenous zoledronate therapy, respectively, is required for 50% of patients to develop ONJ. After stratification by disease, a time of 5.3 and 2.2 years of BP therapy is required for 50% of patients with osteoporosis and cancer, respectively, to develop ONJ. These findings have significant implications for the design of future clinical studies and the development of risk-reduction strategies aimed at either assessing or modulating the risk of ONJ associated with BP.

Influence of APOE Gene Polymorphism on Interindividual and Interethnic Warfarin Dosage Requirement: A Systematic Review and Meta-Analysis

BACKGROUND

Warfarin is the most extensively used coumarin anticoagulant. It has been shown that the anticoagulant effect of warfarin is associated with genetic variation. Apolipoprotein E (ApoE) is a possible candidate to influence the maintenance dose of warfarin. ApoE affects the vitamin K cycle by mediating the uptake of vitamin K into the liver. The vitamin K cycle is the drug target of warfarin. However, the association between genetic variants of the APOE gene and warfarin dose requirement is still controversial.

METHODS

Revman 5.3 software was used to analyze the relationship between APOE genotypes and warfarin dose requirements.

RESULTS

In our meta-analysis, the E2/E2 genotype was significantly associated with warfarin dose. E2/E2 patients required 12% (P = 0.0002) lower mean daily warfarin dose than E3/E3 carriers. In addition, subgroup analysis showed that Asians with the E4/E4 genotype tended to need lower warfarin maintenance doses, while the African American E4/E4 carriers needed slightly higher doses than E3/E3 carriers; however, these subgroups were very small.

CONCLUSION

This is the first meta-analysis of the association between APOE genotypes and warfarin dose. APOE E2/E2 might be one of the factors affecting warfarin dose requirements. The effect of APOE may vary between ethnicities.

Genetic variants associated with angiotensin-converting enzyme inhibitor-induced cough: a genome-wide association study in a Swedish population

Aim: We conducted a genome-wide association study on angiotensin-converting enzyme inhibitor-induced cough and used our dataset to replicate candidate genes identified in previous studies. Patients & methods: A total of 124 patients and 1345 treated controls were genotyped using Illumina arrays. The genome-wide significance level was set to p < 5 × 10-8. Results: We identified nearly genome-wide significant associations in CLASP1, PDE11A, KCNMB2, TGFA, SLC38A6 and MMP16. The strongest association was with rs62151109 in CLASP1 (odds ratio: 3.97; p = 9.44 × 10-8). All top hits except two were located in intronic or noncoding DNA regions. None of the candidate genes were significantly associated in our study. Conclusion: Angiotensin-converting enzyme inhibitor-induced cough is potentially associated with genes that are independent of bradykinin pathways.

Genetic determinants of warfarin maintenance dose and time in therapeutic treatment range: a RE-LY genomics substudy

AIMS

We investigated associations between genetic variation in candidate genes and on a genome-wide scale with warfarin maintenance dose, time in therapeutic range (TTR), and risk of major bleeding.

MATERIALS & METHODS

In total, 982 warfarin-treated patients from the RE-LY trial were studied.

RESULTS

After adjusting for SNPs in VKORC1 and CYP2C9, SNPs in DDHD1 (rs17126068) and NEDD4 (rs2288344) were associated with dose. Adding these SNPs and CYP4F2 (rs2108622) to a base model increased R(2) by 2.9%. An SNP in ASPH (rs4379440) was associated with TTR (-6.8% per minor allele). VKORC1 was associated with time less than INR 2.0. VKORC1 and CYP2C9 were associated with time more than INR 3.0, but not with major bleeding.

CONCLUSIONS

We identified two novel genes associated with warfarin maintenance dose and one gene associated with TTR. These genes need to be replicated in an independent cohort.

Genetic variants associated with antithyroid drug-induced agranulocytosis: a genome-wide association study in a European population

BACKGROUND

Drug-induced agranulocytosis is a potentially life-threatening adverse reaction. Genome-wide association studies (GWASs) in ethnic Chinese people in Taiwan and Hong Kong have shown an association between agranulocytosis induced by antithyroid drugs and the HLA alleles HLA-B*38:02 and HLA-DRB1*08:03. We aimed to identify genetic variants associated with antithyroid drug-induced agranulocytosis in a white European population.

METHODS

We did a GWAS in 234 European adults with any non-chemotherapy drug-induced agranulocytosis (absolute neutrophil count ≤0·5 × 10(9)/L [≤500/μL]) and 5170 population controls. 39 of the 234 patients had agranulocytosis that was induced by antithyroid drugs (thiamazole [methimazole], carbimazole, or propylthiouracil). After imputation and HLA allele prediction, 9 380 034 single nucleotide polymorphisms (SNPs) and 180 HLA alleles were tested for association. The genome-wide significance threshold was p<5 × 10(-8).

FINDINGS

Agranulocytosis induced by non-chemotherapy drugs in general was significantly associated with the HLA region on chromosome 6, with odds ratios (ORs) of 3·24 (95% CI 2·31-4·55, p=1·20 × 10(-11)) for HLA-B*27:05 and 3·57 (2·61-4·90, p=2·32 × 10(-15)) for the top SNP (rs114291795). Drug-specific analysis showed that the association with HLA-B*27:05 was largely driven by cases induced by antithyroid drugs. In a multiple logistic regression model, the OR for HLA-B*27:05 was 7·30 (3·81-13·96) when antithyroid drug-induced agranulocytosis was compared with population controls (p=1·91 × 10(-9)) and 16·91 (3·44-83·17) when compared with a small group of hyperthyroid controls (p=5·04 × 10(-4)). Three SNPs were strongly associated with antithyroid drug-induced agranulocytosis: rs652888 (OR 4·73, 95% CI 3·00-7·44, p=1·92 × 10(-11)) and rs199564443 (17·42, 7·38-41·12, p=7·04 × 10(-11)), which were independent of HLA-B*27:05, and rs1071816 (5·27, 3·06-9·10, p=2·35 × 10(-9)) which was in moderate linkage disequilibrium with HLA-B*27:05. In heterozygous carriers of all three SNPs, the predicted probability of antithyroid drug-induced agranulocytosis was about 30% (OR 753, 95% CI 105-6812). To avoid one case of agranulocytosis, based on the possible risk reduction if all three SNPs are genotyped and carriers are treated or monitored differently from non-carriers, roughly 238 patients would need to be genotyped.

INTERPRETATION

In white European people, antithyroid drug-induced agranulocytosis was associated with HLA-B*27:05 and with other SNPs on chromosome 6. In the future, carriers of these variants could be placed under intensified monitoring or offered alternative treatment for hyperthyroidism.

FUNDING

Swedish Research Council, Swedish Heart and Lung Foundation, Clinical Research Support at Uppsala University, German Federal Institute for Drugs and Medical Devices, Carlos III Spanish Health Institute, European Regional Development Fund, UK National Institute for Health Research, The Selander's Foundation, Thuréus Foundation, European Commission, and Science for Life Laboratory.

Novel regulatory variant detected on the VKORC1 haplotype that is associated with warfarin dose

AIM

Warfarin dose requirement is associated with VKORC1 rs9923231, and we studied whether it is a functional variant.

MATERIALS & METHODS

We selected variants in linkage disequilibrium with rs9923231 that bind transcription factors in an allele-specific way. Representative haplotypes were cloned or constructed, nuclear protein binding and transcriptional activity were evaluated.

RESULTS

rs56314408C>T and rs2032915C>T were detected in a liver enhancer in linkage disequilibrium with rs9923231. The rs56314408-rs2032915 C-C haplotype preferentially bound nuclear proteins and had higher transcriptional activity than T-T and the African-specific T-C. A motif for TFAP2A/C was disrupted by rs56314408T. No difference in transcriptional activity was detected for rs9923231G>A.

CONCLUSION

Our results supported an activating role for rs56314408C, while rs9923231G>A had no evidence of being functional.

Characterizing variability in warfarin dose requirements in children using modelling and simulation

Aims

Although genetic, clinical and demographic factors have been shown to explain approximately half of the inter-individual variability in warfarin dose requirement in adults, less is known about causes of dose variability in children. This study aimed to identify and quantify major genetic, clinical and demographic sources of warfarin dose variability in children using modelling and simulation.

Methods

Clinical, demographic and genetic data from 163 children with a median age of 6.3 years (range 0.06–18.9 years), covering over 183 years of warfarin therapy and 6445 INR observations were used to update and optimize a published adult pharmacometric warfarin model for use in children.

Results

Genotype effects in children were found to be comparable with what has been reported for adults, with CYP2C9 explaining up to a four-fold difference in dose (CYP2C9 *1/*1 vs. *3/*3) and VKORC1 explaining up to a two-fold difference in dose (VKORC1 G/G vs. A/A), respectively. The relationship between bodyweight and warfarin dose was non-linear, with a three-fold difference in dose for a four-fold difference in bodyweight. In addition, age, baseline and target INR, and time since initiation of therapy, but not CYP4F2 genotype, had a significant impact on typical warfarin dose requirements in children.

Conclusions

The updated model provides quantitative estimates of major clinical, demographic and genetic factors impacting on warfarin dose variability in children. With this new knowledge more individualized dosing regimens can be developed and prospectively evaluated in the pursuit of improving both efficacy and safety of warfarin therapy in children.

Genome-wide association study of warfarin maintenance dose in a Brazilian sample

AIM

Extreme discordant phenotype and genome-wide association (GWA) approaches were combined to explore the role of genetic variants on warfarin dose requirement in Brazilians.

METHODS

Patients receiving low (≤ 20 mg/week; n = 180) or high stable warfarin doses (≥ 42.5 mg/week; n = 187) were genotyped with Affymetrix Axiom(®) Biobank arrays. Imputation was carried out using data from the combined 1000 Genomes project.

RESULTS

Genome-wide signals (p ≤ 5 × 10(-8)) were identified in the well-known VKORC1 (lead SNP, rs749671; OR: 20.4; p = 1.08 × 10(-33)) and CYP2C9 (lead SNP, rs9332238, OR: 6.8 and p = 4.4 × 10(-13)) regions. The rs9332238 polymorphism is in virtually perfect LD with CYP2C9*2 (rs1799853) and CYP2C9*3 (rs1057910). No other genome-wide significant regions were identified in the study.

CONCLUSION

We confirmed the important role of VKORC1 and CYP2C9 polymorphisms in warfarin dose. Original submitted 14 January 2015; Revision submitted 26 May 2015.

Oral anticoagulation: a critique of recent advances and controversies

There have recently been significant advances in the field of oral anticoagulation, but these have also led to many controversies. Warfarin is still the commonest drug used for clotting disorders but its use is complicated owing to wide inter-individual variability in dose requirement and its narrow therapeutic index. Warfarin dose requirement can be influenced by both genetic and environmental factors. Two recent randomized controlled trials (RCTs) came to different conclusion regarding the utility of genotype-guided dosing; we critically explore the reasons for the differences. The new generation of oral anticoagulants have been demonstrated to be as efficacious as warfarin, but further work is needed to evaluate their safety in real clinical settings.

A Bayesian decision support tool for efficient dose individualization of warfarin in adults and children

BACKGROUND

Warfarin is the most widely prescribed anticoagulant for the prevention and treatment of thromboembolic events. Although highly effective, the use of warfarin is limited by a narrow therapeutic range combined with a more than ten-fold difference in the dose required for adequate anticoagulation in adults. An optimal dose that leads to a favourable balance between the wanted antithrombotic effect and the risk of bleeding as measured by the prothrombin time International Normalised Ratio (INR) must be found for each patient. A model describing the time-course of the INR response can be used to aid dose selection before starting therapy (a priori dose prediction) and after therapy has been initiated (a posteriori dose revision).

RESULTS

In this paper we describe a warfarin decision support tool. It was transferred from a population PKPD-model for warfarin developed in NONMEM to a platform independent tool written in Java. The tool proved capable of solving a system of differential equations that represent the pharmacokinetics and pharmacodynamics of warfarin with a performance comparable to NONMEM. To estimate an a priori dose the user enters information on body weight, age, baseline and target INR, and optionally CYP2C9 and VKORC1 genotype. By adding information about previous doses and INR observations, the tool will suggest a new dose a posteriori through Bayesian forecasting. Results are displayed as the predicted dose per day and per week, and graphically as the predicted INR curve. The tool can also be used to predict INR following any given dose regimen, e.g. a fixed or an individualized loading-dose regimen.

CONCLUSIONS

We believe that this type of mechanism-based decision support tool could be useful for initiating and maintaining warfarin therapy in the clinic. It will ensure more consistent dose adjustment practices between prescribers, and provide efficient and truly individualized warfarin dosing in both children and adults.

Incorporation of pharmacogenomics into routine clinical practice: the Clinical Pharmacogenetics Implementation Consortium (CPIC) guideline development process

The Clinical Pharmacogenetics Implementation Consortium (CPIC) publishes genotype-based drug guidelines to help clinicians understand how available genetic test results could be used to optimize drug therapy. CPIC has focused initially on well-known examples of pharmacogenomic associations that have been implemented in selected clinical settings, publishing nine to date. Each CPIC guideline adheres to a standardized format and includes a standard system for grading levels of evidence linking genotypes to phenotypes and assigning a level of strength to each prescribing recommendation. CPIC guidelines contain the necessary information to help clinicians translate patient-specific diplotypes for each gene into clinical phenotypes or drug dosing groups. This paper reviews the development process of the CPIC guidelines and compares this process to the Institute of Medicine’s Standards for Developing Trustworthy Clinical Practice Guidelines.