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−7) 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−7) and rs4427239 (OR = 5.47, 95%CI: 2.81–10.65, p = 5.75 × 10−7), 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.This research was funded by Swiss National Science Foundation, grant number 31003A_160206. The EuDAC study was funded by Carlos III Spanish Health Institute [FIS10/02632], the European Regional Development Fund FEDER, the Swedish Research Council [Medicine 521-2011-2440,521-2014-3370 and 2018-03307], the Swedish Heart and Lung Foundation [20120557,20140291 and 20170711] and a grant from the Federal Institute for Drugs and Medical Devices (Bonn, Germany). The EUDRAGENE collaboration has received support from the EC 5th Framework program [QLRI-CR-2002-02757] and the Serious Adverse Events Consortium, SAEC, a collaboration for academia and industry. M.M. is supported by the National Institute for Health Research Biomedical Research Center at Guy’s and St Thomas’ National Health Service Foundation Trust and King’s College London

Drug-Induced Liver Injury due to Flucloxacillin:Relevance of Multiple Human Leukocyte Antigen Alleles

© 2019 The Authors Clinical Pharmacology & Therapeutics © 2019 American Society for Clinical Pharmacology and Therapeutics Some patients prescribed flucloxacillin (~0.01%) develop drug-induced liver injury (DILI). HLA-B*57:01 is an established genetic risk factor for flucloxacillin DILI. To consolidate this finding, identify additional genetic factors, and assess relevance of risk factors for flucloxacillin DILI in relation to DILI due to other penicillins, we performed a genomewide association study involving 197 flucloxacillin DILI cases and 6,835 controls. We imputed single-nucleotide polymorphism and human leukocyte antigen (HLA) genotypes. HLA-B*57:01 was the major risk factor (allelic odds ratio (OR)=36.62; P=2.67×10−97). HLA-B*57:03 also showed an association (OR=79.21; P=1.2×10−6). Within the HLA-B protein sequence, imputation showed valine97, common to HLA-B*57:01 and HLA-B*57:03, had the largest effect (OR=38.1; P=9.7×10−97). We found no HLA-B*57 association with DILI due to other isoxazolyl penicillins (n=6) or amoxicillin (n=15) and no significant non-HLA signals for any penicillin-related DILI

Common variation near IRF6 is associated with IFN-β-induced liver injury in multiple sclerosis

Multiple sclerosis (MS) is a disease of the central nervous system treated with disease-modifying therapies, including the biologic, interferon-β (IFN-β). Up to 60% of IFN-β-exposed MS patients develop abnormal biochemical liver test results1,2, and 1 in 50 experiences drug-induced liver injury3. Since genomic variation contributes to other forms of drug-induced liver injury4,5, we aimed to identify biomarkers of IFN-β-induced liver injury using a two-stage genome-wide association study. The rs2205986 variant, previously linked to differential expression of IRF6, surpassed genome-wide significance in the combined two-stage analysis (P = 2.3 × 10-8, odds ratio = 8.3, 95% confidence interval = 3.6-19.2). Analysis of an independent cohort of IFN-β-treated MS patients identified via electronic medical records showed that rs2205986 was also associated with increased peak levels of aspartate aminotransferase (P = 7.6 × 10-5) and alkaline phosphatase (P = 4.9 × 10-4). We show that these findings may be applicable to predicting IFN-β-induced liver injury, offering insight into its safer use

Genome-wide association study of angioedema induced by angiotensin-converting enzyme inhibitor and angiotensin receptor blocker treatment

Angioedema in the mouth or upper airways is a feared adverse reaction to angiotensin-converting enzyme inhibitor (ACEi) and angiotensin receptor blocker (ARB) treatment, which is used for hypertension, heart failure and diabetes complications. This candidate gene and genome-wide association study aimed to identify genetic variants predisposing to angioedema induced by these drugs. The discovery cohort consisted of 173 cases and 4890 controls recruited in Sweden. In the candidate gene analysis, ETV6, BDKRB2, MME, and PRKCQ were nominally associated with angioedema (p < 0.05), but did not pass Bonferroni correction for multiple testing (p < 2.89 × 10−5). In the genome-wide analysis, intronic variants in the calcium-activated potassium channel subunit alpha-1 (KCNMA1) gene on chromosome 10 were significantly associated with angioedema (p < 5 × 10−8). Whilst the top KCNMA1 hit was not significant in the replication cohort (413 cases and 599 ACEi-exposed controls from the US and Northern Europe), a meta-analysis of the replication and discovery cohorts (in total 586 cases and 1944 ACEi-exposed controls) revealed that each variant allele increased the odds of experiencing angioedema 1.62 times (95% confidence interval 1.05–2.50, p = 0.030). Associated KCNMA1 variants are not known to be functional, but are in linkage disequilibrium with variants in transcription factor binding sites active in relevant tissues. In summary, our data suggest that common variation in KCNMA1 is associated with risk of angioedema induced by ACEi or ARB treatment. Future whole exome or genome sequencing studies will show whether rare variants in KCNMA1 or other genes contribute to the risk of ACEi- and ARB-induced angioedema

Association between the number of coadministered P-glycoprotein inhibitors and serum digoxin levels in patients on therapeutic drug monitoring

BACKGROUND: The ABC transporter P-glycoprotein (P-gp) is recognized as a site for drug-drug interactions and provides a mechanistic explanation for clinically relevant pharmacokinetic interactions with digoxin. The question of whether several P-gp inhibitors may have additive effects has not yet been addressed. METHODS: We evaluated the effects on serum concentrations of digoxin (S-digoxin) in 618 patients undergoing therapeutic drug monitoring. P-gp inhibitors were classified as Class I, with a known effect on digoxin kinetics, or Class II, showing inhibition in vitro but no documented effect on digoxin kinetics in humans. Mean S-digoxin values were compared between groups of patients with different numbers of coadministered P-gp inhibitors by a univariate and a multivariate model, including the potential covariates age, sex, digoxin dose and total number of prescribed drugs. RESULTS: A large proportion (47%) of the digoxin patients undergoing therapeutic drug monitoring had one or more P-gp inhibitor prescribed. In both univariate and multivariate analysis, S-digoxin increased in a stepwise fashion according to the number of coadministered P-gp inhibitors (all P values < 0.01 compared with no P-gp inhibitor). In multivariate analysis, S-digoxin levels were 1.26 ± 0.04, 1.51 ± 0.05, 1.59 ± 0.08 and 2.00 ± 0.25 nmol/L for zero, one, two and three P-gp inhibitors, respectively. The results were even more pronounced when we analyzed only Class I P-gp inhibitors (1.65 ± 0.07 for one and 1.83 ± 0.07 nmol/L for two). CONCLUSIONS: Polypharmacy may lead to multiple drug-drug interactions at the same site, in this case P-gp. The S-digoxin levels increased in a stepwise fashion with an increasing number of coadministered P-gp inhibitors in patients taking P-gp inhibitors and digoxin concomitantly. As coadministration of digoxin and P-gp inhibitors is common, it is important to increase awareness about P-gp interactions among prescribing clinicians

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

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

Pharmacogenomics of Antihypertensive Treatment &amp; Clinical Pharmacological Studies of Digoxin Treatment

In Part I we found that the CYP2C9 genotype appears to influence the diastolic blood pressure response to the angiotensin II-receptor antagonist irbesartan in patients with hypertension and left ventricular hypertrophy. Those with the *1/*2 genotype (slower metabolism) responded better than those with the *1/*1 genotype (normal metabolism), likely due to a slower elimination of the drug. We further found that a +9/-9 exon 1 polymorphism of the B2 bradykinin receptor gene – shown to affect mRNA expression - appears to influence the regression of left ventricular mass during therapy with irbesartan or the beta-blocker atenolol in the same patients. Subjects with the -9/-9 genotype (higher mRNA expression) had a greater regression than carriers of the +9 allele. In Part II we found that women on digoxin therapeutic drug monitoring have higher serum digoxin concentrations (SDCs) as compared to men (1.54±0.04 [nmol/L±SE] vs 1.20±0.05 [nmol/L±SE], p&lt;0.001), which could be of importance since an SDC &gt;1.4 nmol/L has been associated with increased mortality. We further found that coadministration of P-glycoprotein inhibitors with digoxin was common (47%) among the same patients, and that the SDC increased in a stepwise fashion with the number of P-glycoprotein inhibitors (20-60%). Lastly, we found that patients admitted to Swedish coronary care units with atrial fibrillation without heart failure and who had been given digoxin had a higher 1-year mortality than those not given digoxin (RR 1.44 [95% CI 1.29-1.60], adjustment made for potential confounders). In conclusion, Part I represents a further step in the pharmacogenomic prospect of tailoring antihypertensive therapy. Part II indicates that heightened attention to the digoxin-dose is warranted in women, that there is a need for awareness about P-glycoprotein interactions with digoxin, and that long-term therapy with digoxin is an independent risk factor for death among patients with atrial fibrillation without heart failure