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.

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.

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.

[Bisphosphonate-associated atypical fractures of the femur: an update of the current evidence]

Atypical fracture of the femur is a well-documented adverse reaction to antiresorptive treatment with bisphosphonates. Although there has been significant gain of knowledge during the past decade, the pathogenesis of this type of fracture is still poorly understood. We present an update on the evidence in regard to epidemiology, pathophysiology, and management of atypical fractures.

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.

A Missense Variant in PTPN22 is a Risk Factor for Drug-induced Liver Injury

BACKGROUND & AIMS

We performed genetic analyses of a multiethnic cohort of patients with idiosyncratic drug-induced liver injury (DILI) to identify variants associated with susceptibility.

METHODS

We performed a genome-wide association study of 2048 individuals with DILI (cases) and 12,429 individuals without (controls). Our analysis included subjects of European (1806 cases and 10,397 controls), African American (133 cases and 1,314 controls), and Hispanic (109 cases and 718 controls) ancestry. We analyzed DNA from 113 Icelandic cases and 239,304 controls to validate our findings.

RESULTS

We associated idiosyncratic DILI with rs2476601, a nonsynonymous polymorphism that encodes a substitution of tryptophan with arginine in the protein tyrosine phosphatase, nonreceptor type 22 gene (PTPN22) (odds ratio [OR] 1.44; 95% confidence interval [CI] 1.28-1.62; P = 1.2 × 10^-9 and replicated the finding in the validation set (OR 1.48; 95% CI 1.09-1.99; P = .01). The minor allele frequency showed the same effect size (OR > 1) among ethnic groups. The strongest association was with amoxicillin and clavulanate-associated DILI in persons of European ancestry (OR 1.62; 95% CI 1.32-1.98; P = 4.0 × 10^-6; allele frequency = 13.3%), but the polymorphism was associated with DILI of other causes (OR 1.37; 95% CI 1.21-1.56; P = 1.5 × 10^-6; allele frequency = 11.5%). Among amoxicillin- and clavulanate-associated cases of European ancestry, rs2476601 doubled the risk for DILI among those with the HLA risk alleles A*02:01 and DRB1*15:01.

CONCLUSIONS

In a genome-wide association study, we identified rs2476601 in PTPN22 as a non-HLA variant that associates with risk of liver injury caused by multiple drugs and validated our finding in a separate cohort. This variant has been associated with increased risk of autoimmune diseases, providing support for the concept that alterations in immune regulation contribute to idiosyncratic DILI.

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.

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.

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.

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

Identification of Risk Factors for Bisphosphonate-Associated Atypical Femoral Fractures and Osteonecrosis of the Jaw in a Pharmacovigilance Database

BACKGROUND

Atypical femoral fractures (AFs) and osteonecrosis of the jaw (ONJ) are well-known adverse drug reactions (ADRs) associated with bisphosphonates. To prevent these ADRs and to aid in the search for pathogenic mechanisms, knowledge of risk factors can be helpful.

OBJECTIVE

To identify risk factors for bisphosphonate-related ONJ and AF.

METHODS

In this case-control study of reports of bisphosphonate-related ADRs from February 16, 1984, to October 16, 2013, in the Swedish national database of ADRs, we compared characteristics for cases of ONJ (n = 167) and AF (n = 55) with all other bisphosphonate-related ADRs (n = 565) with regard to demographic variables, clinical characteristics, and concomitant drug treatments. We adjusted for multiple comparisons with Bonferroni correction.

RESULTS

Time to onset of ADRs differed statistically significantly between cases of AF and controls (2156 vs 111 days). For ONJ versus controls, differences were statistically significant for time to onset (1240 vs 111 days), intravenous administration (40% vs 20%), dental procedures (49% vs 0.2%) and prostheses (5% vs 0%), cancer disease (44% vs 12%), multiple myeloma (21% vs 1%), rheumatoid arthritis (14% vs 5%), and treatment with antineoplastic agents and oxycodone.

CONCLUSION

These results lend further evidence to previously identified risk factors for ONJ-that is, intravenous bisphosphonate administration; invasive dental procedures and dental prostheses; cancer disease, in particular multiple myeloma; and possibly, long-term bisphosphonate treatment. A putative further risk factor is rheumatoid arthritis. Only long-term bisphosphonate treatment was more common among AF cases. The lack of overlap of risk factors between ONJ and AF suggests different pathogenic mechanisms.

Mortality After Atypical Femoral Fractures: A Cohort Study

Although osteoporotic fracture rates can be reduced by bisphosphonates, prolonged therapy is associated with higher risk of atypical femoral fractures. Ordinary fragility fractures are linked to high mortality rates. We aimed to determine whether atypical femoral fractures also confer excess mortality. Radiographs were reviewed for all patients aged ≥55 years who had experienced a subtrochanteric or femoral shaft fracture in Sweden in 2008 to 2010. The fractures were classified as either atypical or ordinary. Data on medication use, coexisting conditions, and date of death were obtained from national registers. We estimated multivariable-adjusted relative risks of death after atypical femoral fractures compared with ordinary subtrochanteric or femoral shaft fractures and calculated age- and sex-standardized mortality ratios (SMRs) for atypical and ordinary fractures compared with the population average. During a mean of 4 years of follow-up, 39 of 172 (23%) patients with an atypical fracture had died compared with 588 of 952 (62%) with an ordinary fracture, corresponding to a relative risk of 0.51 (95% confidence interval [CI] 0.38-0.68). The lower risk was evident in both users and nonusers of bisphosphonates. No patient with atypical fracture died in the first year after fracture. Individuals with an ordinary fracture had a higher mortality risk than the general population (SMR = 1.82; 95% CI 1.69-1.99), but no excess risk was found in patients with atypical fracture (SMR = 0.92; 95% CI 0.65-1.26). We conclude that in contrast to ordinary subtrochanteric and femoral shaft fractures, atypical femoral fractures are not associated with excess mortality.

Prodromal Symptoms in Patients with Bisphosphonate-Associated Atypical Fractures of the Femur

Symptoms have been reported to precede bisphosphonate-associated atypical fractures (AFs) of the femoral shaft. We aimed to determine the frequency and clinical characteristics of such prodromal symptoms. We searched the Swedish national database of spontaneously reported adverse drug reactions for all cases of AF associated with bisphosphonates from January 2006 to March 2013. To confirm diagnostic accuracy and to characterize and determine the frequency of any prodromal symptoms we retrieved copies of medical journals and radiographs for patients who consented to participate in the study. The frequency of prodromal symptoms was compared with that of patients where information was based only on narratives from the adverse drug reaction case reports. A total of 45 reports of AF were identified. We were able to obtain medical records and x-rays for 21 cases and diagnostic accuracy was confirmed for all. Medical records revealed prodromal symptoms in 86 % (n = 18), most commonly pain in the ipsilateral thigh (14 out of 18 patients) preceding the fracture for weeks or longer. Awareness of such symptoms may facilitate early diagnosis and possible prevention of the AF.

Phenotype standardization of angioedema in the head and neck region caused by agents acting on the angiotensin system

Angioedema is a potentially life-threatening adverse reaction to angiotensin-converting enzyme inhibitors and angiotensin receptor blockers. To study the genetic etiology of this rare adverse event, international consortia and multicenter recruitment of patients are needed. To reduce patient heterogeneity, we have standardized the phenotype. In brief, it comprises swelling in the head and neck region that first occurs during treatment. It should not coincide with urticaria or have another likely cause such as hereditary angioedema.

Identification of risk factors for carbamazepine-induced serious mucocutaneous adverse reactions: A case-control study using data from spontaneous adverse drug reaction reports

OBJECTIVES

To identify risk factors other than genetic for severe carbamazepine-induced mucocutaneous reactions, that is, SJS, TEN, and exfoliative dermatitis (ED).

MATERIALS AND METHODS

We did a case-control study using data from the Swedish national database of spontaneously reported adverse drug reactions (ADRs). We selected all patients who had been reported from January 1, 1965 to March 31, 2010 as having experienced SJS (n = 78), TEN (n = 6), or ED (n = 8), and assessed as at least possibly related to carbamazepine. We also included diagnoses possibly representative of early signs of these serious conditions, that is, erythema multiforme (EM, n = 34) and scaly rash (n = 13). We compared data on demographics, drug treatment, and clinical features for these patients (cases, n = 139) with those from patients who had experienced any other type of ADR from carbamazepine during the same time period (controls, n = 887).

RESULTS

After adjustment for multiple comparisons, alcohol abuse was statistically significantly more common among cases than controls (34.5% vs 8.7%, odds ratio 5.5 [95% confidence interval 3.6-8.4], P = 3.14 × 10(-14) ). The same was seen for SJS and EM individually.

CONCLUSION

Alcohol abuse is a possible risk factor for serious carbamazepine-induced mucocutaneous reactions.

Surgical approach to snus-induced injury of the oral mucosa

Snus (Swedish moist snuff) causes lesions in the oral mucosa at the location where pinches are regularly placed. In addition, some patients develop irreversible local gingival recession and sometimes ulcers with perforations to the roots. Such injuries lead to denuded roots that are at risk for caries and periodontal disease, with subsequent esthetic consequences. Therapy for irreversible local gingival recession is currently lacking. In the present report, we describe two cases of successful surgical treatment for irreversible lesions caused by snus.