Next Article in Journal
The Role of Salmonella Genomic Island 4 in Metal Tolerance of Salmonella enterica Serovar I 4,[5],12:i:- Pork Outbreak Isolate USDA15WA-1
Next Article in Special Issue
A Review of the Important Role of CYP2D6 in Pharmacogenomics
Previous Article in Journal
Platinum-Quality Mitogenome Haplotypes from United States Populations
Previous Article in Special Issue
Digital Health Applications for Pharmacogenetic Clinical Trials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by Anca Liliana Cismaru 1,2, Deborah Rudin 3,4, Luisa Ibañez 5, Evangelia Liakoni 6,7, Nicolas Bonadies 8, Reinhold Kreutz 9, Alfonso Carvajal 10, Maria Isabel Lucena 11, Javier Martin 12, Esther Sancho Ponce 13, Mariam Molokhia 14, Niclas Eriksson 15, EuDAC collaborators, Stephan Krähenbühl 3, Carlo R. Largiadèr 1, Manuel Haschke 6,7, Pär Hallberg 16, Mia Wadelius 16 and Ursula Amstutz 1,*
1
Department of Clinical Chemistry, Inselspital Bern University Hospital, University of Bern, 3010 Bern, Switzerland
2
Graduate School for Cellular and Biomedical Sciences, University of Bern, 3012 Bern, Switzerland
3
Department of Clinical Pharmacology & Toxicology, University Hospital Basel, University of Basel, 4031 Basel, Switzerland
4
Department of Biomedicine, University of Basel, 4051 Basel, Switzerland
5
Clinical Pharmacology Service, Hospital Universitari Vall d’Hebron, Department of Pharmacology, Therapeutics and Toxicology, Autonomous University of Barcelona, Fundació Institut Català de Farmacología, 08035 Barcelona, Spain
6
Department of Clinical Pharmacology & Toxicology, Inselspital Bern University Hospital, University of Bern, 3010 Bern, Switzerland
7
Institute of Pharmacology, University of Bern, 3012 Bern, Switzerland
8
Department of Hematology and Central Hematology Laboratory, Inselspital Bern University Hospital, University of Bern, 3010 Bern, Switzerland
9
Charité–Universitätsmedizin Berlin, Corporate Member of Freie Universität Berlin, Humboldt-Universität zu Berlin, and Berlin Institute of Health, Institut für Klinische Pharmakologie und Toxikologie, 10117 Berlin, Germany
10
Centro de Estudios sobre la Seguridad de los Medicamentos, Universidad de Valladolid, 47005 Valladolid, Spain
11
Servicio Farmacologia Clinica, Instituto de Investigación Biomedica de Málaga, Hospital Universitario Virgen de la Victoria, Universidad de Málaga, 29010 Málaga, Spain
12
Instituto de Parasitología y Biomedicina Lopez-Neyra, Consejo Superior de Investigaciones Cientiíficas, 18016 Granada, Spain
13
Servei d’Hematologia i Banc de Sang, Hospital General de Catalunya, 08190 Sant Cugat del Vallès, Spain
14
Department of Population Health Sciences, King’s College London, London WC2R 2LS, UK
15
Uppsala Clinical Research Center and Department of Medical Sciences, Uppsala University, 751 85 Uppsala, Sweden
16
Department of Medical Sciences, Clinical Pharmacology and Science for Life Laboratory, Uppsala University, 751 85 Uppsala, Sweden
*
Author to whom correspondence should be addressed.
EuDAC collaborative authors are listed in the Acknowledgements.
Submission received: 18 September 2020 / Revised: 23 October 2020 / Accepted: 27 October 2020 / Published: 29 October 2020

Abstract

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

Graphical Abstract

1. Introduction

Metamizole, or dipyrone, is an analgesic and antipyretic drug used to manage different types of pain as well as fever and often serves as an alternative to therapy with conventional non-steroidal anti-inflammatory drugs (NSAIDs). Despite a well-known favorable gastrointestinal safety profile [1,2,3], spontaneous adverse drug reaction (ADR) reports with risk estimates varying between different countries [4,5,6] provide accounts of metamizole-induced blood dyscrasias. Specifically, these include metamizole-induced neutropenia (MIN) and a more severe form, agranulocytosis (MIA) characterized by a decrease in circulating neutrophil granulocytes below 1.5 × 109 cells/L, and 0.5 × 109 cells/L, respectively [7,8]. Consequently, the use of metamizole is currently subjected to contrasting regulations put in place by different government authorities, ranging from market withdrawal to non-prescription use [4,5,9,10,11,12]. To date, the pathogenesis and biologic markers of risk for metamizole-induced agranulocytosis (MIA) remain poorly understood and there are no effective strategies for prediction or prevention.
Genome-wide association studies (GWAS) have served as a cost-effective approach to identify associations of single nucleotide polymorphisms (SNPs) with a variety of human phenotypes ranging from complex diseases to drug-related outcomes. In particular, previous GWAS have identified associations of genetic variants with rather large effect sizes for the susceptibility to various rare and serious adverse drug reactions (ADRs) such as flucloxacillin-induced liver injury, statin-induced myopathy and carbamazepine-induced skin injury [13,14]. In more recent years, studies of idiosyncratic drug-induced agranulocytosis identified genes associated with increased risk either in the human leukocyte antigen (HLA) region or in other regions involved in immune responses for culprit drugs such as clozapine, carbimazole and sulfasalazine [15,16,17,18,19].
Although the biological significance of many of the reported associations has only partially been unmasked, these findings would tend to favor the involvement of the immune system in these adverse drug reactions. Conversely, for metamizole, recent evidence from mechanistic in vitro studies is not in agreement with an immune-system driven mechanistic hypothesis as it has been shown that the main metamizole metabolite MAA (N-methyl-4-aminoantipyrine) reacts with the hemoglobin break down product hemin to produce a reactive intermediate that is cytotoxic for promyelocytic HL60 cells, but not mature neutrophil granulocytes [20,21,22]. Given this conflicting evidence, exploring genetic associations with MIA could be helpful to shed more light on whether to support or reject either hypothesis.
To identify loci associated with metamizole-induced neutropenia/agranulocytosis, we performed single-marker tests of association by means of a GWAS using genotype data from the largest MIA/MIN cohort to date (86 cases, 464 controls). Genetic variants in 84 candidate genes were analyzed in a Swiss discovery cohort (45 cases, 191 controls), with replication of findings and a genome-wide meta-analysis in two independent cohorts genotyped by the European Drug-induced Agranulocytosis Consortium (EuDAC).

2. Materials and Methods

2.1. Ethical Statement

This study was approved by the local ethics committees “Ethikkommission Nordwest- und Zentralschweiz” and “Kanton Bern” in Switzerland (EKNZ BASEC 2015–00231, KEK-Nr 2015-00231). Research was conducted according to the latest update of the Declaration of Helsinki. Collection of biological materials and clinical information was undertaken with a written informed consent from all participants. Samples and data of German (EuDAC-DE) and Spanish (EuDAC-ES) patient cohorts had previously been collected through the European Drug-induced Agranulocytosis Consortium (EuDAC) with approval by the local ethics committees (22 December 2014, Málaga, Spain; RTF011, Barcelona, Spain; Charité–Universitätsmedizin Berlin, Germany).

2.2. Study Design and Participants

This retrospective observational case-control study included genotype data from a Swiss discovery cohort (MIA/MIN-CH) that consisted of a total of 53 cases, 39 tolerant controls and 161 unexposed controls recruited at the University Hospitals Basel and Bern. MIN and MIA were defined by an absolute neutrophil count (ANC) below 1.5 × 109 cells/L, and 0.5 × 109 cells/L respectively. Tolerant controls included in the study had received at least 500 mg metamizole per day for a minimum of twenty-eight consecutive days, a treatment duration encompassing the latency time observed for a majority of cases for the occurrence of MIA based on a previous report [6]. Clinical drug tolerability during metamizole therapy was assessed by the absence of symptoms including fever, sore throat, or mucositis [23]. EDTA blood samples for genetic analyses were collected and coded at the time of recruitment. Clinical data, including patient characteristics and concomitant drug therapy (such as antibiotics, analgesics and β-lactam antibiotics) were retrieved from medical charts. A more detailed description of participants recruited for this study has been previously published [23]. In addition, participants genotyped as part of a project of the EuDAC were assigned to two independent replication cohorts based on their recruitment site: Germany (EuDAC-DE) and Spain (EuDAC-ES). Additional description of these cohorts is available in previous reports [16]. The majority of subjects were of self-reported European ancestry, which was consistent with multidimensional scaling (MDS) analyses of population structure using genome-wide genotype data (see further details below).

2.3. Genotype Data and Quality Control

For the samples of the MIA/MIN-CH discovery cohort, genomic DNA from whole blood or buffy coat was extracted and purified using the QIAamp DNA Blood Maxi or Midi Kits (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Eluted DNA was quantified with the Qubit 4 Fluorometer (ThermoFisher Scientific, Waltham, MA, USA) and quality was assessed by measuring absorbance at A230, A260 and A280 on a NanoDrop 1000 Spectrophotometer (ThermoFisher Scientific).
Genotyping in each of the contributing cohorts was conducted independently at different time points using various high-density SNP arrays (Figure 1). For the MIA/MIN-CH cohort, genotyping was performed using the Infinium Human CoreExome-24 BeadChip and processed by the iScan System together with a customized Tecan liquid-handling robot (both Illumina, San Diego, CA, USA) at the iGE3 Facility (University of Geneva, Geneva, Switzerland). Genotype calls were generated with the GenomeStudio software (Illumina). For replication and meta-analyses, existing genotype data was obtained through the EuDAC.
In all datasets, SNP positions were based on National Center for Biotechnology Information (NCBI, Bethesda, MD, USA) build 37 (hg19) and alleles were labeled on the positive strand of the reference genome. Similar genotyping quality control (QC) procedures [24,25,26,27,28] were used for each cohort using PLINK v1.9 [29]. Specifically, individual samples were removed if more than 2% of SNPs failed genotyping (call rate < 98%), and individual SNPs were removed if they showed a Hardy-Weinberg Equilibrium (HWE) p-value < 0.001 [30,31,32,33,34], a minor allele frequency (MAF) < 1%, if more than 2% of samples were flagged as having failed genotype calling for any given SNP (call rate < 98%) and if any duplicate or triplicate discordance was detected.
The rs111876221 polymorphism in the SERINC5 gene was additionally genotyped in the MIA-CH cohort using Polymerase Chain Reaction (PCR) TaqMan qPCR Master Mix 2× and a custom TaqMan SNP Genotyping assay kit (cat. 4331349, Assay ID-ANZTJD4, ThermoFisher Scientific), according to the manufacturer’s recommendations. Real-time PCR reactions were performed using 10 μL final reaction volumes, consisting of 5 μL 2× master mix, 4.75 μL water, 0.25 μL 40× probe assay and 1.0 μL gDNA. The amplification was performed with the QuantStudio™ 6 Flex System (ThermoFisher Scientific). The PCR program used was as follows: 10 min at 95 °C with at least 40× (15 s at 95 °C and 1 min at 60 °C). TaqMan assay results were validated (seven samples homozygous for the major allele, all samples with heterozygous genotype and other samples with an undetermined genotype) using Sanger sequencing. Details pertaining to the SNP primer sequences utilized for genotype validation can be found in Supplementary Table S1.

2.4. Multidimensional Scaling and Identification of Genetic Outliers

Multi-dimensional scaling (MDS) analysis was performed independently for each cohort using PLINK 1.9 [29]. SNPs that passed quality control were pruned using the linkage disequilibrium (LD) pruning parameters of r2 < 0.2 over a window size of 50. Genome-wide identity-by-descent (IBD) given identity-by-state (IBS) information was calculated using all SNPs that remained after pruning. Five nearest neighbors were identified for each individual based upon the pairwise IBS distance and the distance to the nth nearest neighbor was converted to a Z score. Fourteen individuals (nine in MIA/MIN-CH, two in EuDAC-ES and three in EuDAC-DE) with a minimum Z score less than −2 among the five nearest neighbors were excluded from analysis as population outliers. Furthermore, three pairs of EuDAC individuals (one in EuDAC-DE and two in EuDAC-ES) showed high levels of IBD sharing and one individual from each pair was randomly selected to be excluded from subsequent analyses due to this cryptic relatedness.

2.5. Imputation

Imputation of SNPs was performed by minimac3 [35,36,37] on the Michigan Imputation Server [38] using the Haplotype Reference Consortium (HRC 1.1 2016) reference panel [39] on the basis of pre-imputation SNP-level validated data (https://www.well.ox.ac.uk/~wrayner/tools/) from each cohort (Supplementary Table S2). HRC-based imputation was performed separately for each cohort, which increased the genome-wide SNP densities to 39′004′932, 38′994′811 and 39′008′189 SNPs for MIA/MIN-CH, EuDAC-ES and EuDAC-DE respectively. After performing quality control as described above, a total of 7′608′978, 7′675′043 and 7′422′709 SNPs were retained for the meta-analysis performed using a fixed-effects model in the MIA-CH, EuDAC-ES and EuDAC-DE cohorts respectively.

2.6. Candidate Gene and Genome-Wide Association Analyses

All analyses were carried out with PLINK (version 2.0) using logistic regression with the first four MDS dimensions and sex (only for EuDAC-ES) chosen as covariates. An additive genetic model was assumed for all variants and the genomic inflation factor was calculated for each analysis to assess potential confounding effects of population stratification [40,41].
Association analyses were carried out in multiple steps: Initial discovery analyses were carried out both in the entire MIA/MIN-CH cohort (N = 45 MIA/MIN cases, N = 191 tolerant/unexposed controls) and including only agranulocytosis cases of this cohort (N = 30 MIA cases, N = 191 tolerant/unexposed controls) focusing first on specific candidate gene sets followed by genome-wide association analyses. Candidate variants identified in the MIA/MIN-CH cohort were subsequently investigated for replication in the EuDAC cohorts. In a second discovery analysis, a genome-wide meta-analysis was carried out including association summary statistics from all three cohorts (Figure 1). In all analyses, tolerant and unexposed controls were combined as a control group in the MIA/MIN-CH cohort.
To perform candidate gene analyses, we analyzed different genes following a two-stage approach. First, eight genes were previously identified by association studies [15,42,43,44,45,46,47] of agranulocytosis induced by other drugs (NOX3, SERINC5, NQO2, GSTM1, SLCO1B3, SLCO1B7, TAP2, AGER; Phase I set) were selected and then 50 additional candidate genes (Phase II set) were included for their hypothesized relationship to drug-induced agranulocytosis based on their functional annotation potentially related to the oxidative metabolism of metamizole [48], the anti-oxidant defense [20,22] and the immune system defense response specific to genes expressed in granulocytes [49]. For each gene, we examined SNPs in a window of 10kb upstream and downstream of the gene. We used a statistical significance threshold p-value based on the method of Li et al. [34], which calculates the effective number of independent tests (Me) and applies a Bonferroni correction based on that number. A list of all genes part of the Phase II gene set is reported in a separate Supplementary Table S3.
For all genome-wide analyses, a threshold p-value of 1 × 10−7 was used to declare whether an association signal was statistically significant at a genome-wide level. This threshold was used based on an evaluation of empirical replication success of discovery GWAS, demonstrating good replication success for borderline associations with p < 1 × 10−7 [50]. The SNP association p-values from the three cohorts were subjected to meta-analysis with METAL using an inverse variance weighted scheme [51]. Identified candidate loci were further investigated for other traits associated with variation in these regions using publicly available databases of GWAS results (https://genetics.opentargets.org) and for effects on regulatory motifs or on expression from eQTL studies [52].

3. Results

3.1. Cohort Characteristics

Of the 96 cases of metamizole-induced agranulocytosis identified, ten were excluded after adjudication: eight cases in MIA/MIN-CH (two because they were identified as suggested genetic outlier i.e., nearest neighbor analyses suggesting non-European ancestry, two because of ongoing infectious diseases, four because of immunosuppressive therapy with cytotoxic drugs) and two cases in EuDAC-ES (because of suggested genetic outlier). No cases were excluded in the EuDAC-DE cohort.
Of the 479 controls identified, data from fifteen samples were excluded after adjudication: nine controls in MIA/MIN-CH (one individual was recruited twice as determined by IBD analysis, seven individuals because they were identified as suggested genetic outlier, and one because of sample genotype call rate < 50%), two controls in EuDAC-ES (because of cryptic relatedness) and four controls in EuDAC-DE (three because of suggested genetic outlier and one because of cryptic relatedness). Demographic characteristics of the tolerant and unexposed controls in the discovery cohort (MIA/MIN-CH) are presented in Supplementary Table S4. Comparison of demographic and clinical factors between the three independent cohorts are shown in Table 1. Cases were balanced by sex and age. Data on sex and age was unavailable for controls in the EuDAC-DE and EuDAC-ES cohorts.

3.2. Association Analyses in the Discovery Cohort (MIA/MIN-CH and MIA-CH)

3.2.1. Candidate Gene Analyses

In the Swiss discovery cohort of 236 individuals, the genomic inflation factor (λ) was λ = 1.00 and no systemic bias was detected in any of the analyses that were conducted. The standard per-SNP significance threshold for a family-wise error rate (FWER) of α = 0.05 was estimated at 5.4 × 10−4 for the phase I candidate gene analysis and 1.08 × 10−4 for phase II. No statistically significant SNPs associated with MIA/MIN and MIA, respectively were identified in either phase I or phase II when including all 45 MIA and MIN cases (Supplementary Tables S5 and S6) or only the 30 MIA cases (Supplementary Tables S7 and S8). The top finding in both phase I analyses was a common intronic variant in the serine incorporator 5 (SERINC5) gene (rs10041917, p = 3.4 × 10−3, OR = 0.44, 95%-CI 0.25–0.76). This signal remained the top finding also in the analysis of the extended number of candidate genes in phase II. The rs10041917 variant in SERINC5 was in linkage disequilibrium with a variant (rs111876221) previously reported to be associated with sulfasalazine-induced agranulocytosis (D’ = 1.0, p = 0.0068) as determined using the LDpair tool in LDlink using 1000 Genomes project data from European populations (https://ldlink.nci.nih.gov/?tab=ldpair). To evaluate a potential association with MIN/MIA, rs111876221 was separately genotyped in the entire MIA/MIN-CH cohort as this SNP was not included on the Infinium array and imputation results were only obtained in a later stage of the data analysis. Direct genotyping of rs111876221 in the MIA/MIN-CH cohort revealed that the A allele at that locus, which is linked to the A allele of rs10041917, was found more frequently in the case group (MAF = 0.055) compared to the controls (MAF = 0.016) indicating the same direction of effect. This association was shown to be significant with a one-sided Fisher’s exact test p-value of 0.03 (OR = 3.85, 95% 1.12–13.25).
Among the top SNPs from the phase I analysis in agranulocytosis cases only (Supplementary Table S7), several additional variants were identified as of potential interest. Specifically, the same missense SNP in the ATP binding cassette subfamily B member (TAP2) gene, rs2228391 (p = 5.1 × 10−3, OR = 377.7, 95% CI 5.93–24028), as observed here at a higher frequency in MIA cases was previously associated with susceptibility of a northern Chinese Han population to antithyroid drug-induced agranulocytosis [53]. In addition, another missense rs4148876 SNP in TAP2 (p = 6.1 × 10−3, OR = 3.06, 95% CI 1.37–6.83) has been reported to be correlated with gene expression of the HLA-DOB and TAP2 genes in several tissues including whole blood [54,55,56,57], Similarly, the intergenic variant rs6588432 SNP near the glutathione peroxidase 7 (GPX7) gene, has been reported to be associated with gene expression of GPX7 and origin recognition complex subunit 1 (ORC1L) in several tissues including blood [54,57]. Due to their potential functional relevance, these additional variants were selected to assess replication in the EuDAC cohorts.

3.2.2. Genome-Wide Association Analyses

We performed two genome-wide association analyses of 304,704 genotyped variants, first in the entire MIA/MIN-CH cohort (N = 45 cases, N = 191 tolerant/unexposed controls) and then including only MIA cases (N = 30 cases) versus the same 191 controls. None of the SNPs showed an association at the genome-wide significance level (p < 1 × 10−7) in these analyses (Figure 2). In the analysis including both MIA and MIN cases (Figure 2a), the leading SNP with the best evidence for a suggestive association was in an intergenic region on chromosome 6 (rs191786, p = 8.6 × 10−6, OR = 0.24, 95% CI 0.13–0.45) near the TENT5A gene. In the analysis with only MIA cases (Figure 2b), the leading SNP was also in an intergenic region on chromosome 6 (rs9366076, p = 2.89 × 10−6, OR = 4.81, 95% CI 2.49–9.30).
Interestingly, this lead SNP was previously reported to be associated with gene expression of the nearby ribonuclease T2 (RNASET2) gene in multiple tissues [55,57,58] and in neutrophils [58]. Supplementary Tables S9 and S10 summarizes the leading associations of the genome-wide analyses.

3.2.3. Replication in Independent Cohorts

To evaluate potential replication of these candidate associations in independent cohorts, directly typed and imputed genotype data from the EuDAC-ES and EuDAC-DE cohorts was used (Supplementary Table S11). Based on genome-wide association analyses for 1,786,726 genotyped SNPs in the EuDAC-ES cohort and 650,136 SNPs in the EuDAC-DE cohort, no evidence of systemic bias was detected in these datasets with genomic inflation factors (λ) of λ = 1.016 and λ = 1.047 in the EuDAC-ES and EuDAC-DE cohorts, respectively.
For the candidate SNPs selected from the phase I and phase II candidate gene analyses in the MIA/MIN-CH and MIA-CH cohort, we did not observe any significant associations with MIA in the two independent EuDAC cohorts (Supplementary Table S11). Due to reduced imputation accuracy in the EuDAC-ES (RSQ = 0.74) and EuDAC-DE (RSQ = 0.26) cohorts, replication of the rs111876221 variant previously associated with sulfasalazine-induced agranulocytosis could not be reliably assessed in those cohorts. Specifically, the reduced genotype imputation accuracy for this relatively uncommon SNP was illustrated in the MIA/MIN-CH cohort where the genotype was not correctly imputed for two out of eleven heterozygous minor allele carriers when compared to direct genotyping, in spite of an imputation RSQ of 0.91.
Associations of the leading SNPs from both genome-wide analyses in the MIA/MIN-CH and MIA-CH discovery cohort (Supplementary Tables S9 and S10), rs191786 and rs9366076 near RNASET2, previously reported as an eQTL of this gene, could not be replicated in the two independent European cohorts (Supplementary Table S11).

3.3. GWAS Meta-Analysis across Independent Cohorts

For the meta-analyses, genome-wide association results from all three cohorts based on imputed genotype data were used. Also, in the meta-analysis, two analyses were conducted, first including MIN cases in the MIA/MIN-CH cohort and a second analysis only considering MIA cases. Manhattan plots of the results of both meta-analyses are shown in Figure 3. The first meta-analysis performed in all 86 cases and 464 controls across all three cohorts revealed that a leading SNP with the best evidence for association in an intronic region of the transforming growth factor β receptor 3 (TGFBR3) gene on chromosome 1 (rs11583606, p = 1.72 × 10−7, OR = 7.95% CI 3.37–14.5) (Supplementary Figures S1 and S2A). The second strongest independent signal was positioned in an intronic region in the protein tyrosine phosphatase receptor type O (PTPRO) gene on chromosome 12 (rs112917452, p = 9.92 × 10−7, OR = 7.3, 95% CI 3.29–16.20) (Supplementary Figures S2B and S3). No association reached statistical significance at the genome-wide level (Table 2).
The second meta-analysis (Table 3) performed in the 71 MIA cases and 464 controls revealed a leading SNP with the best evidence for association in an intergenic region on chromosome 9 (rs55898176, p = 1.01 × 10−7, OR = 4.01, 95% CI 2.41–6.68) (Supplementary Figure S4A), only marginally not meeting the genome-wide significance threshold. This SNP was identified to be located near a long non-coding RNA with CAAP1 as the closest coding gene at a distance of 177 kilobases to the canonical Transcript Start Site (TSS) of the gene (Figure 4).
The second strongest independent signal was positioned in an intronic region in the Sushi von Willebrand Factor Type A, EGF and Pentraxin Domain Containing 1 (SVEP1) gene, also on chromosome 9 (rs4427239, p = 5.75 × 10−7, OR = 5.47, 95% CI 2.81–10.65) (Supplementary Figure S4B). Several other genes located in proximity of SVEP1 include the thioredoxin (TXN) gene, involved in the cellular antioxidant defense (Figure 4). An evaluation of the frequencies of the associated alleles for the lead signals from both meta-analyses in the MIN vs. MIA cases of the MIA/MIN-CH cohort is shown in Supplementary Tables S5 and S6, respectively. Interestingly, the observed allele frequencies of the leading SNPs show that the risk allele distributions of rs11583606 and rs112917452 were higher in the MIN cases than in the MIA cases (Supplementary Figure S5). In contrast, the risk allele distribution of rs55898176 and rs4427239 were higher in the MIA cases than in the MIN cases (Supplementary Figure S6). Post stage GWAS analysis of these candidate loci using HaploReg [52] and GTEx [59] did not indicate any potential regulatory effects associated with the top SNPs across the available gene expression quantitative trait loci (eQTL) datasets. Furthermore, no GWAS lead variants were found to be linked to any of the candidate loci. However, several hematological measurement traits such as hemoglobin concentration, autoimmune hemolytic anemias and monocyte count reported by previous studies were found to be associated (p < 0.05) with the rs112917452, rs55898176 and rs4427239 variants in the UK Biobank and GWAS catalog summary statistics repository respectively. Also, several association studies [60] reported associations between lead variants in the SVEP1 gene and hematological measurements obtained from blood plasma including neutrophil percentage of granulocytes (p = 4.5 × 10−14), white blood cell count (p = 2.0 × 10−15), platelet count (p = 2.2 × 10−11) and eosinophil percentage of granulocytes (p = 9.1 × 10−17).

4. Discussion

We performed candidate gene-based and genome-wide association analyses for metamizole-induced agranulocytosis in three datasets totaling 86 MIN/MIA cases and 464 controls. To our knowledge, this is the first study to investigate genetic associations with MIA at a genome-wide level in the largest patient cohort available to date. In contrast to previous GWAS for agranulocytosis related to other pharmaceutical agents, we did not identify any genome-wide significant genetic associations that would implicate an immune-mediated mechanism for MIA. However, we report suggestive evidence for an association of two loci on chromosome 9, one of which implicates anti-oxidant components as well as hematopoiesis, that could serve as plausible candidates for studies to come.
Overall, the genetic susceptibility to rare yet potentially lethal adverse drug reactions, ranging from skin and liver injury to agranulocytosis, has been repeatedly shown to be specific to each offending drug [61,62,63,64,65]. Nevertheless, common underlying T cell-mediated mechanisms have been identified for some drug-induced serious skin reactions and drug-induced liver injury. Similarly, the strongest pharmacogenetic associations pertaining to drug-induced agranulocytosis up to date mostly involved immune-related genes such as specific HLA alleles [15,16,17,44,45,46,66,67], even though these findings have yet to be corroborated by mechanistic studies. Here, on the other hand, in line with a previous association study focused on classical HLA genes in the same cohort [68], none of the candidate loci identified as leading markers in our genome-wide meta-analysis involved genes that would support an immune-mediated mechanism underlying MIA. These findings thus suggest that the underlying mechanism for MIA may differ from other agranulocytosis-inducing drugs. Such a lack of evidence for an immune-mediated mechanism underlying MIA is concordant with recent in vitro studies suggesting that the cytotoxicity in immature myeloid cells stems from a reactive electrophilic intermediate formed from the main metabolite of metamizole, MAA and hemin [20,21,22]. Genetic variation in enzymes involved in the heme degradation and antioxidant defense pathways could thus impact an individual’s susceptibility to develop MAA-associated myelotoxicity.
In this context, the suggestive association of the SNP rs4427239 (chr9p13) in the analysis of only MIA cases vs. tolerant/population controls is notable due to the reported associations of variants at this locus with hematopoietic traits and the role of its nearby gene TXN, encoding for thioredoxin (TXN), in redox control [69,70,71,72,73,74,75]. As part of the thioredoxin system, TXN is an antioxidant enzyme found in a variety of cells including granulocytes [76,77], and is involved in the defense against oxidative stress by controlling cellular free radicals and reactive oxygen species [78,79,80,81]. Functional investigations of TXN in vitro cell culture models, in particular involving granulocyte precursor cells may thus be of interest [20,22]. While mRNA expression of SVEP1, the gene in which the lead SNP rs4427239 is located, was not detected in any immune cells, it was reported to be overexpressed in a hematopoiesis-supporting splenic stromal cell line [82]. Furthermore, its encoding protein SVEP1 was identified as a ligand of integrin α 9 β 1 (ITGA9) [83], which is involved in signal transduction in hematopoietic stem cells [84], and for which an important role has been suggested for granulopoiesis [85]. Given these findings, SVEP1 also represents an interesting starting point for experimental investigations aiming to clarify the impact of metamizole and its metabolites on cell adhesion between granulocyte precursors and stromal cells during maturation.
While the signal related to rs4427239 thus includes two nearby genes with biological plausibility, the strongest association signal we observed in the meta-analysis including only MIA cases was the intergenic variant rs55898176 (chr9p13) located near an uncharacterized long non-coding RNA (LOC105376000). The nearest coding gene to this variant is CAAP1, which encodes conserved anti-apoptotic proteins that modulate the mitochondrial apoptosis pathway [86]. The possible biological implications of this association with respect to MIA thus still remains to be elucidated.
Interestingly, the frequencies of both of the above-mentioned lead SNPs, rs4427239 and rs55898176, were substantially lower in the Swiss MIN cases compared to MIA cases and similar to the frequencies observed in controls (Supplementary Table S6). The allele frequencies of these SNPs in the MIN cases thus did not follow an intermediate distribution between the frequencies in the MIA cases and the controls, as has been observed for genetic associations with other adverse drug reactions for cases with intermediate severity [87]. If these signals indeed represent true causal associations, a potential underlying mechanism involving the intergenic locus on chromosome 9, or the SVEP1 locus in agranulocytosis may not apply to metamizole-induced neutropenia.
Conversely, we observed that the frequencies of the risk alleles for the rs11583606 and rs112917452 SNPs identified as leading markers the meta-analysis of all cases versus controls were substantially higher in the Swiss MIN cases compared to the Swiss MIA cases (Supplementary Figure S5). These signals thus appear to be driven primarily by the MIN cases in the Swiss cohort. The leading loci from this meta-analysis, TGFBR3 and PTPRO thus showed less consistent effects among the MIA cases across the three cohorts and as a consequence may not represent the most promising signals for further investigation. Similarly, both candidate gene analyses (phase I and phase II) revealed leading markers in the Swiss cohort positioned in intronic regions of the same gene, SERINC5, a gene coding for transmembrane proteins reducing Human Immunodeficiency Virus (HIV-1) infectivity [88]. In this gene, a genome-wide significant intronic SNP rs111876221 was previously reported to be associated with sulfasalazine-induced agranulocytosis [15]. However, given that the associations in this gene were not replicated in the Spanish or German cohorts and the lack of obvious functional implications of the SERINC5 gene as it relates to MIA, we do not consider this signal as a strong candidate for further investigation.
None of the identified lead SNPs were located in protein-coding regions or were found to be eQTLs in previous studies including the 48 human tissues analyzed in the eQTL analysis of the Genotype-Tissue Expression (GTEx) project [59]. Importantly, however, it should be considered that the biological interpretation of the candidate loci potentially associated with MIA is hampered by the scarcity of functional data on the impact of genetic variation on gene expression specifically in bone marrow or in granulocyte precursors, where cytotoxic effects of metamizole may take place [20,22]. As a consequence, potential regulatory effects of these candidate SNPs at the level of the bone marrow and specifically during granulocyte maturation cannot be excluded. With the increasing availability of single-cell RNA sequencing studies providing data at the level of individual cell types or cell differentiation stages [89,90,91], additional functional annotation of the identified candidate loci may be an option in the future.
Main strengths of this study are the fairly homogenous discovery cohort and the possibility of replication in a second and third cohort. Although the case-control cohorts originated from three different European countries, by limiting our analyses to individuals of European ancestry using genetic information and adjusting for residual population structure, the likelihood of false positive associations due to population stratification was reduced. Also, by adhering to the case definition criteria of the EuDAC studies for the selection of individuals of the MIA-CH cohort, the three cohorts were similar in their definition of the ADR phenotype, which provided the advantage of increased power for the genome-wide meta-analysis.
Although this study has several strengths, it should also be interpreted in the context of several limitations. First, with the limited sample size, our study was not powered to detect associations with smaller effect size for uncommon variants. Despite the assembly of large patient cohorts by large international consortia, studies on rare idiosyncratic adverse drug reactions with incidences of 1 in 10,000 patients or less, finding it difficult to recruit large enough numbers of cases and tolerant comparison controls, have been faced with the reality of possibly being underpowered and to date have commonly included less than hundred cases [92]. Nevertheless, several of these studies have successfully identified associations due to large effect sizes being commonly observed for adverse drug reactions compared to genetic associations with common diseases [87,92,93,94,95]. To try to attenuate this limitation, recruitment of more cases and tolerant controls would have been a wishful but very challenging endeavor given the rarity of MIA. Second, another limitation lies in the use of population controls instead of drug-tolerant controls. However, the use of population controls instead of tolerant controls in genetic association studies is suitable for rare ADRs (incidence < 1%) as is the case for MIA [96]. Third, any residual selection bias stemming from the different recruitment channels for cases and controls cannot be fully ruled out, even though no systematic genetic differences between cases and controls were identified by MDS. Finally, considering the lack of statistically significant associations and the lack of an additional independent replication cohort to further investigate the lead signals from the genome-wide meta-analysis, it is possible that some of the identified candidate loci, or possibly even all, represent false positive associations, despite their consistent effects observed in the three investigated cohorts. Similarly, genetic variation associated with MIA should be explored in other populations as well, as the frequency of the identified variants and their clinical impact could vary between different ethnicities.

5. Conclusions

In this study, we associated genotyped and imputed markers with MIA in three independent cohorts of European individuals. Taken together, although our findings do not have immediate clinical utility, they point towards novel candidate genes with a possible role in the mechanism of MIA in individuals of European ancestry. Interestingly, our findings do not implicate biological pathways that have been previously identified in association studies for other causal drugs of agranulocytosis studied in populations of the same ethnicity, for which current evidence is consistent with an immune-mediated mechanism. Although further studies are warranted to follow up on these suggestive candidate genes and pathways, our findings are a stepping-stone in the ongoing endeavor of better understanding the genetic background behind MIA and its underlying causal mechanisms. Ultimately, this could pave the way to a better screening or treatment for a safer use of metamizole.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/11/11/1275/s1, Figure S1: Regional association plot for regions around rs11583606, Figure S2: Forest plot of associations with metamizole-induced agranulocytosis and neutropenia in the three independent cohorts (all cases and controls), Figure S3: Regional association plot for regions around rs112917452, Figure S4: Forest plot of associations with metamizole-induced agranulocytosis in the three independent cohorts (agranulocytosis cases and controls), Figure S5: Allele frequency distribution of the leading SNPs associated with metamizole-induced agranulocytosis in the three independent cohorts (all cases versus controls), Figure S6: Allele frequency distribution of the leading SNPs associated with metamizole-induced agranulocytosis in the three independent cohorts (agranulocytosis cases versus controls), Table S1: Primers used for confirmation of rs111876221 SNP in MIA-CH cohort, Table S2: Pre-imputation data checking, Table S3 (provided as separate file): Phase II gene set, Table S4: Characteristics of cases, tolerant and unexposed controls in the discovery cohort (MIA/MIN-CH), Table S5: Top associations of the Phase I association analysis in the MIA/MIN-CH cohort, Table S6: Top associations of the Phase II association analysis in the MIA/MIN-CH cohort, Table S7: Top associations of the Phase I association analysis in the MIA-CH cohort, Table S8: Top associations of the Phase II association analysis in the MIA-CH cohort, Table S9: Top associations of the genome-wide association analysis in the MIA/MIN-CH cohort, Table S10: Top associations of the genome-wide association analysis in the MIA-CH cohort, Table S11: Replication of association signals from MIA/MIN-CH and MIA-CH analyses in the EuDAC cohorts.

Author Contributions

Conceptualization, S.K., M.H., P.H., M.W. and U.A.; Data curation, A.L.C. and N.E.; Formal analysis, A.L.C.; Funding acquisition, L.I., R.K., M.M., S.K., M.H., P.H., M.W. and U.A.; Investigation, A.L.C. and U.A.; Methodology, A.L.C., C.R.L. and U.A.; Project administration, A.L.C., M.H. and U.A.; Resources, D.R., L.I., E.L., N.B., R.K., A.C., M.I.L., J.M., E.S.P., M.M., S.K., C.R.L., M.H., P.H., M.W. and U.A.; Software, A.L.C.; Supervision, C.R.L., M.H. and U.A.; Validation, A.L.C. and U.A.; Visualization, A.L.C.; Writing—original draft, A.L.C.; Writing—review & editing, A.L.C., D.R., L.I., E.L., N.B., R.K., A.C., M.I.L., J.M., E.S.P., M.M., N.E., S.K., C.R.L., M.H., P.H., M.W. and U.A. All authors have read and agreed to the published version of the manuscript.

Funding

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.

Acknowledgments

The authors would like to cordially thank the laboratory technicians of the University Institute of Clinical Chemistry Gisela Andrey, Daniel Schärer, and Daniela Sommer as well as our master student, Livia Grimm, for their support with the DNA extractions. We further thank Zoltan Kutalik and Eleonora Porcu for their helpful suggestions concerning imputation and meta-analysis. Last but not least, we would like to extend our thanks to William Rayner for his kind and valuable troubleshooting support on pre-imputation data handling. Finally, we thank all study participants without whom this effort would not have been possible. The graphical abstract was created with Biorender.com. EuDAC collaborative authors: Lourdes Vendrell (Fundació Institut Català de Farmacologia, Hospital Universitari Vall d’Hebron, Barcelona, Spain); Ramon Puig Treserra (Fundació Institut Català de Farmacologia, Barcelona, Spain); Jose Luis Caro (Banc de Sang i Teixits, Barcelona, Spain); Eduard Palou (Banc de Sang i Teixits, Barcelona, Spain); María José Herrero (Banc de Sang i Teixits, Barcelona, Spain); Esmeralda de la Banda Ledrado (Hospital Universitario Bellvitge, Barcelona, Spain); Eva Montané (Hospital Universitari Germans Trias i Pujol, Universitat Autònoma de Barcelona); Elisa Orna Montero (Institut Català d’Oncologia-Hospital Germans Trias I Pujol, Badalona, Spain); José Tomás Navarro Ferrando (Catalan Institute of Oncology; Francesc Salvador Rudilla, Banc de Sang i Teixits, Barcelona, Spain); Consuelo Pedrós Cholvi (Hospital Universitario Bellvitge, Barcelona, Spain).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Arellano, F.; Sacristán, J.A. Metamizole: Reassessment of its therapeutic role. Eur. J. Clin. Pharmacol. 1990, 38, 617–619. [Google Scholar] [CrossRef]
  2. Andrade, S.E.; Martinez, C.; Walker, A.M. Comparative safety evaluation of non-narcotic analgesics. J. Clin. Epidemiol. 1998, 51, 1357–1365. [Google Scholar] [CrossRef]
  3. Sánchez, S.; De la Lastra, C.A.; Ortiz, P.; Motilva, V.; Martín, M.J. Gastrointestinal tolerability of metamizol, acetaminophen, and diclofenac in subchronic treatment in rats. Dig. Dis. Sci. 2002, 47, 2791–2798. [Google Scholar] [CrossRef]
  4. Hedenmalm, K.; Spigset, O. Agranulocytosis and other blood dyscrasias associated with dipyrone (metamizole). Eur. J. Clin. Pharmacol. 2002, 58, 265–274. [Google Scholar] [CrossRef]
  5. Ibáñez, L.; Vidal, X.; Ballarín, E.; Laporte, J.R. Agranulocytosis associated with dipyrone (metamizol). Eur. J. Clin. Pharmacol. 2005, 60, 821–829. [Google Scholar] [CrossRef]
  6. Blaser, L.S.; Tramonti, A.; Egger, P.; Haschke, M.; Krähenbühl, S.; Rätz Bravo, A.E. Hematological safety of metamizole: Retrospective analysis of WHO and Swiss spontaneous safety reports. Eur. J. Clin. Pharmacol. 2015, 71, 209–217. [Google Scholar] [CrossRef] [PubMed]
  7. Blaser, L.; Hassna, H.; Hofmann, S.; Holbro, A.; Haschke, M.; Rätz Bravo, A.; Zeller, A.; Krähenbühl, S.; Taegtmeyer, A. Leucopenia associated with metamizole: A case-control study. Swiss Med. Wkly. 2017, 147, w14438. [Google Scholar] [CrossRef] [PubMed]
  8. Andrès, E.; Maloisel, F. Idiosyncratic drug-induced agranulocytosis or acute neutropenia. Curr. Opin. Hematol. 2008, 15, 15–21. [Google Scholar] [CrossRef] [PubMed]
  9. Van Der Klauw, M.M.; Goudsmit, R.; Halie, M.R.; Van’t Veer, M.B.; Herings, R.M.C.; Wilson, J.H.P.; Stricker, B.H.C. A population-based case-cohort study of drug-associated agranulocytosis. Arch. Intern. Med. 1999, 159, 369–374. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Basak, G.W.; Drozd-Sokołowska, J.; Wiktor-Jedrzejczak, W. Update on the incidence of metamizole sodium-induced blood dyscrasias in Poland. J. Int. Med. Res. 2010, 38, 1374–1380. [Google Scholar] [CrossRef]
  11. Chan, T.Y.K.; Chan, A.W.K. Aminopyrine-induced blood dyscrasias-still a problem in many parts of the world. Pharmacoepidemiol. Drug Saf. 1996, 5, 215–219. [Google Scholar] [CrossRef]
  12. Stammschulte, T.; Ludwig, W.D.; Mühlbauer, B.; Bronder, E.; Gundert-Remy, U. Metamizole (dipyrone)-associated agranulocytosis. An analysis of German spontaneous reports 1990–2012. Eur. J. Clin. Pharmacol. 2015, 71, 1129–1138. [Google Scholar] [CrossRef] [PubMed]
  13. Zhou, K.; Pearson, E.R. Insights from Genome-Wide Association Studies of Drug Response. Annu. Rev. Pharmacol. Toxicol. 2013, 53, 299–310. [Google Scholar] [CrossRef]
  14. Daly, A.K. Genome-wide association studies in pharmacogenomics. Nat. Rev. Genet. 2010, 11, 241–246. [Google Scholar] [CrossRef] [PubMed]
  15. Wadelius, M.; Eriksson, N.; Kreutz, R.; Bondon-Guitton, E.; Ibañez, L.; Carvajal, A.; Lucena, M.I.; Sancho Ponce, E.; Molokhia, M.; Martin, J.; et al. Sulfasalazine-Induced Agranulocytosis Is Associated With the Human Leukocyte Antigen Locus. Clin. Pharmacol. Ther. 2017, 103, 843–853. [Google Scholar] [CrossRef] [Green Version]
  16. Hallberg, P.; Eriksson, N.; Ibañez, L.; Bondon-Guitton, E.; Kreutz, R.; Carvajal, A.; Lucena, M.I.; Ponce, E.S.; Molokhia, M.; Martin, J.; et al. Genetic variants associated with antithyroid drug-induced agranulocytosis: A genome-wide association study in a European population. Lancet Diabetes Endocrinol. 2016, 4, 507–516. [Google Scholar] [CrossRef] [Green Version]
  17. Goldstein, J.I.; Fredrik Jarskog, L.; Hilliard, C.; Alfirevic, A.; Duncan, L.; Fourches, D.; Huang, H.; Lek, M.; Neale, B.M.; Ripke, S.; et al. Clozapine-induced agranulocytosis is associated with rare HLA-DQB1 and HLA-B alleles. Nat. Commun. 2014, 5, 1–9. [Google Scholar] [CrossRef]
  18. Athanasiou, M.C.; Dettling, M.; Cascorbi, I.; Mosyagin, I.; Salisbury, B.A.; Pierz, K.A.; Zou, W.; Whalen, H.; Malhotra, A.K.; Lencz, T.; et al. Candidate gene analysis identifies a polymorphism in HLA-DQB1 associated with clozapine-induced agranulocytosis. J. Clin. Psychiatry 2011, 72, 458–463. [Google Scholar] [CrossRef]
  19. Pavlos, R.; Mallal, S.; Ostrov, D.; Buus, S.; Metushi, I.; Peters, B.; Phillips, E. T Cell–Mediated Hypersensitivity Reactions to Drugs. Annu. Rev. Med. 2015, 66, 439–454. [Google Scholar] [CrossRef] [Green Version]
  20. Rudin, D.; Lanzilotto, A.; Bachmann, F.; Housecroft, C.E.; Constable, E.C.; Drewe, J.; Haschke, M.; Krähenbühl, S. Non-immunological toxicological mechanisms of metamizole-associated neutropenia in HL60 cells. Biochem. Pharmacol. 2019, 163, 345–356. [Google Scholar] [CrossRef]
  21. Rudin, D.; Haschke, M.; Krahenbuhl, S. Mechanisms of cytotoxicity involved in metamizole-induced neutropenia. Clin. Toxicol. 2017, 55, 371–544. [Google Scholar] [CrossRef]
  22. Rudin, D.; Roos, N.J.; Duthaler, U.; Krähenbühl, S. Toxicity of metamizole on differentiating HL60 cells and human neutrophil granulocytes. Toxicology 2019, 426, 152254. [Google Scholar] [CrossRef] [PubMed]
  23. Rudin, D.; Spoendlin, J.; Cismaru, A.L.; Liakoni, E.; Bonadies, N.; Amstutz, U.; Meier, C.R.; Krähenbühl, S.; Haschke, M. Metamizole-associated neutropenia: Comparison of patients with neutropenia and metamizole-tolerant patients. Eur. J. Intern. Med. 2019, 68, 36–43. [Google Scholar] [CrossRef] [PubMed]
  24. Igo, R.P.; Cooke Bailey, J.N.; Romm, J.; Haines, J.L.; Wiggs, J.L. Quality control for the illumina humanexome beadchip. Curr. Protoc. Hum. Genet. 2016, 90, 2–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Coleman, J.R.I.; Euesden, J.; Patel, H.; Folarin, A.A.; Newhouse, S.; Breen, G. Quality control, imputation and analysis of genome-wide genotyping data from the Illumina HumanCoreExome microarray. Brief. Funct. Genom. 2016, 15, 298–304. [Google Scholar] [CrossRef] [Green Version]
  26. Turner, S.; Armstrong, L.L.; Bradford, Y.; Carlsony, C.S.; Crawford, D.C.; Crenshaw, A.T.; de Andrade, M.; Doheny, K.F.; Haines, J.L.; Hayes, G.; et al. Quality control procedures for genome-wide association studies. Curr. Protoc. Hum. Genet. 2011, 68, 1.19.1–1.19.18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Alonso, N.; Lucas, G.; Hysi, P. Big data challenges in bone research: Genome-wide association studies and next-generation sequencing. Bonekey Rep. 2015, 4, 635. [Google Scholar] [CrossRef]
  28. Marees, A.T.; de Kluiver, H.; Stringer, S.; Vorspan, F.; Curis, E.; Marie-Claire, C.; Derks, E.M. A tutorial on conducting genome-wide association studies: Quality control and statistical analysis. Int. J. Methods Psychiatr. Res. 2018, 27, e1608. [Google Scholar] [CrossRef]
  29. Purcell, S.; Neale, B.; Todd-Brown, K.; Thomas, L.; Ferreira, M.A.R.; Bender, D.; Maller, J.; Sklar, P.; De Bakker, P.I.W.; Daly, M.J.; et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 2007, 81, 559–575. [Google Scholar] [CrossRef] [Green Version]
  30. Reddy, S.; Madireddy, S.; Nair, P. Association of Genetic Markers Contributing to Dyslexia Susceptibility in Indian Population. J. Neurol. Neurosci. 2016. [Google Scholar] [CrossRef]
  31. Graffelman, J.; Nelson, S.; Gogarten, S.M.; Weir, B.S. Exact inference for Hardy-Weinberg proportions with missing genotypes: Single and multiple imputation. G3 Genes Genomes Genet. 2015, 5, 2365–2373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Miyagawa, T.; Nishida, N.; Ohashi, J.; Kimura, R.; Fujimoto, A.; Kawashima, M.; Koike, A.; Sasaki, T.; Tanii, H.; Otowa, T.; et al. Appropriate data cleaning methods for genome-wide association study. J. Hum. Genet. 2008, 53, 886–893. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Graffelman, J.; Jain, D.; Weir, B. A genome-wide study of Hardy–Weinberg equilibrium with next generation sequence data. Hum. Genet. 2017, 136, 727–741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Li, M.X.; Yeung, J.M.Y.; Cherny, S.S.; Sham, P.C. Evaluating the effective numbers of independent tests and significant p-value thresholds in commercial genotyping arrays and public imputation reference datasets. Hum. Genet. 2012, 131, 747–756. [Google Scholar] [CrossRef] [Green Version]
  35. Zhang, B.; Zhi, D.; Zhang, K.; Gao, G.; Limdi, N.A.; Liu, N. Practical consideration of genotype imputation: Sample size, window size, reference choice, and untyped rate. Stat. Interface 2011, 4, 339–351. [Google Scholar] [PubMed]
  36. Li, Y.; Willer, C.J.; Ding, J.; Scheet, P.; Abecasis, G.R. MaCH: Using sequence and genotype data to estimate haplotypes and unobserved genotypes. Genet. Epidemiol. 2010, 34, 816–834. [Google Scholar] [CrossRef] [Green Version]
  37. Marchini, J.; Howie, B. Genotype imputation for genome-wide association studies. Nat. Rev. Genet. 2010, 11, 499–511. [Google Scholar] [CrossRef]
  38. Das, S.; Forer, L.; Schönherr, S.; Sidore, C.; Locke, A.E.; Kwong, A.; Vrieze, S.I.; Chew, E.Y.; Levy, S.; McGue, M.; et al. Next-generation genotype imputation service and methods. Nat. Genet. 2016, 48, 1284–1287. [Google Scholar] [CrossRef] [Green Version]
  39. McCarthy, S.; Das, S.; Kretzschmar, W.; Delaneau, O.; Wood, A.; Teumer, A.; Kang, H.M.; Fuchsberger, C.; Danecek, P.; Sharp, K.; et al. A reference panel of 64,976 haplotypes for genotype imputation. bioRxiv 2015. [Google Scholar] [CrossRef]
  40. Freedman, M.L.; Reich, D.; Penney, K.L.; McDonald, G.J.; Mignault, A.A.; Patterson, N.; Gabriel, S.B.; Topol, E.J.; Smoller, J.W.; Pato, C.N.; et al. Assessing the impact of population stratification on genetic association studies. Nat. Genet. 2004, 36, 388–393. [Google Scholar] [CrossRef]
  41. Devlin, B.; Roeder, K. Genomic control for association studies. Biometrics 1999, 55, 997–1004. [Google Scholar] [CrossRef] [PubMed]
  42. Legge, S.E.; Hamshere, M.L.; Ripke, S.; Pardinas, A.F.; Goldstein, J.I.; Rees, E.; Richards, A.L.; Leonenko, G.; Jorskog, L.F.; Chambert, K.D.; et al. Genome-wide common and rare variant analysis provides novel insights into clozapine-associated neutropenia. Mol. Psychiatry 2017, 22, 1502–1508. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Saito, T.; Ikeda, M.; Mushiroda, T.; Ozeki, T.; Kondo, K.; Shimasaki, A.; Kawase, K.; Hashimoto, S.; Yamamori, H.; Yasuda, Y.; et al. Pharmacogenomic Study of Clozapine-Induced Agranulocytosis/Granulocytopenia in a Japanese Population. Biol. Psychiatry 2016, 80, 636–642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Chen, P.L.; Shih, S.R.; Wang, P.W.; Lin, Y.C.; Chu, C.C.; Lin, J.H.; Chen, S.C.; Chang, C.C.; Huang, T.S.; Tsai, K.S.; et al. Genetic determinants of antithyroid drug-induced agranulocytosis by human leukocyte antigen genotyping and genome-wide association study. Nat. Commun. 2015, 6, 7633. [Google Scholar] [CrossRef] [PubMed]
  45. Plantinga, T.S.; Arts, P.; Knarren, G.H.; Mulder, A.H.; Wakelkamp, I.M.; Hermus, A.R.; Joosten, L.A.; Netea, M.G.; Bisschop, P.H.; de Herder, W.W.; et al. Rare NOX3 Variants Confer Susceptibility to Agranulocytosis During Thyrostatic Treatment of Graves’ Disease. Clin. Pharmacol. Ther. 2017, 102, 1017–1024. [Google Scholar] [CrossRef]
  46. Ostrousky, O.; Meged, S.; Loewenthal, R.; Valevski, A.; Weizman, A.; Carp, H.; Gazit, E. NQO2 gene is associated with clozapine-induced agranulocytosis. Tissue Antigens 2003, 62, 483–491. [Google Scholar] [CrossRef]
  47. Van Der Weide, K.; Loovers, H.; Pondman, K.; Bogers, J.; Van Der Straaten, T.; Langemeijer, E.; Cohen, D.; Commandeur, J.; Van Der Weide, J. Genetic risk factors for clozapine-induced neutropenia and agranulocytosis in a Dutch psychiatric population. Pharm. J. 2017, 17, 471–478. [Google Scholar] [CrossRef]
  48. García-Martín, E.; Esguevillas, G.; Blanca-López, N.; García-Menaya, J.; Blanca, M.; Amo, G.; Canto, G.; Martínez, C.; Cordobés, C.; Agúndez, J.A.G. Genetic determinants of metamizole metabolism modify the risk of developing anaphylaxis. Pharmacogenet. Genom. 2015, 25, 462–464. [Google Scholar] [CrossRef]
  49. Skelly, D.A.; Squiers, G.T.; McLellan, M.A.; Bolisetty, M.T.; Robson, P.; Rosenthal, N.A.; Pinto, A.R. Single-Cell Transcriptional Profiling Reveals Cellular Diversity and Intercommunication in the Mouse Heart. Cell Rep. 2018, 22, 600–610. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Panagiotou, O.A.; Ioannidis, J.P.A.; Hirschhorn, J.N.; Abecasis, G.R.; Frayling, T.M.; McCarthy, M.I.; Lindgren, C.M.; Beaty, T.H.; Eriksson, N.; Polychronakos, C.; et al. What should the genome-wide significance threshold be? Empirical replication of borderline genetic associations. Int. J. Epidemiol. 2012, 41, 273–286. [Google Scholar] [CrossRef]
  51. Willer, C.J.; Li, Y.; Abecasis, G.R. METAL: Fast and efficient meta-analysis of genomewide association scans. Bioinformatics 2010, 26, 2190–2191. [Google Scholar] [CrossRef] [PubMed]
  52. Ward, L.D.; Kellis, M. HaploReg: A resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 2012, 40, D930–D934. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. He, Y.; Zheng, J.; Zhang, Q.; Hou, P.; Zhu, F.; Yang, J.; Li, W.; Chen, P.; Liu, S.; Zhang, B.; et al. Association of HLA-B and HLA-DRB1 polymorphisms with antithyroid drug-induced agranulocytosis in a Han population from northern China. Sci. Rep. 2017, 7, 1–8. [Google Scholar] [CrossRef] [Green Version]
  54. Westra, H.J.; Peters, M.J.; Esko, T.; Yaghootkar, H.; Schurmann, C.; Kettunen, J.; Christiansen, M.W.; Fairfax, B.P.; Schramm, K.; Powell, J.E.; et al. Systematic identification of trans eQTLs as putative drivers of known disease associations. Nat. Genet. 2013, 45, 1238–1243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Traynor, K. Suvorexant approved for insomnia. Am. J. Health Pharm. 2014, 71, 1524. [Google Scholar] [CrossRef] [PubMed]
  56. Lappalainen, T.; Sammeth, M.; Friedländer, M.R.; ’T Hoen, P.A.C.; Monlong, J.; Rivas, M.A.; Gonzàlez-Porta, M.; Kurbatova, N.; Griebel, T.; Ferreira, P.G.; et al. Transcriptome and genome sequencing uncovers functional variation in humans. Nature 2013, 501, 506–511. [Google Scholar] [CrossRef]
  57. Ardlie, K.G.; DeLuca, D.S.; Segrè, A.V.; Sullivan, T.J.; Young, T.R.; Gelfand, E.T.; Trowbridge, C.A.; Maller, J.B.; Tukiainen, T.; Lek, M.; et al. The Genotype-Tissue Expression (GTEx) pilot analysis: Multitissue gene regulation in humans. Science 2015, 348, 648–660. [Google Scholar] [CrossRef] [Green Version]
  58. Andiappan, A.K.; Melchiotti, R.; Poh, T.Y.; Nah, M.; Puan, K.J.; Vigano, E.; Haase, D.; Yusof, N.; San Luis, B.; Lum, J.; et al. Genome-wide analysis of the genetic regulation of gene expression in human neutrophils. Nat. Commun. 2015, 6, 7971. [Google Scholar] [CrossRef] [Green Version]
  59. Carithers, L.J.; Moore, H.M. The Genotype-Tissue Expression (GTEx) Project. Biopreserv. Biobank. 2015, 13, 307–308. [Google Scholar] [CrossRef]
  60. Carvalho-Silva, D.; Pierleoni, A.; Pignatelli, M.; Ong, C.K.; Fumis, L.; Karamanis, N.; Carmona, M.; Faulconbridge, A.; Hercules, A.; McAuley, E.; et al. Open Targets Platform: New developments and updates two years on. Nucleic Acids Res. 2019, 47, D1056–D1065. [Google Scholar] [CrossRef]
  61. Daly, A. Human Leukocyte Antigen (HLA) Pharmacogenomic Tests: Potential and Pitfalls. Curr. Drug Metab. 2014, 15, 196–201. [Google Scholar] [CrossRef] [PubMed]
  62. Johnston, A.; Uetrecht, J. Current understanding of the mechanisms of idiosyncratic drug-induced agranulocytosis. Expert Opin. Drug Metab. Toxicol. 2015, 11, 243–257. [Google Scholar] [CrossRef] [PubMed]
  63. Tesfa, D.; Keisu, M.; Palmblad, J. Idiosyncratic drug-induced agranulocytosis: Possible mechanisms and management. Am. J. Hematol. 2009, 84, 428–434. [Google Scholar] [CrossRef] [PubMed]
  64. Mallal, S.; Phillips, E.; Carosi, G.; Molina, J.-M.; Workman, C.; Tomažič, J.; Jägel-Guedes, E.; Rugina, S.; Kozyrev, O.; Cid, J.F.; et al. HLA-B*5701 Screening for Hypersensitivity to Abacavir. N. Engl. J. Med. 2008, 358, 568–579. [Google Scholar] [CrossRef] [Green Version]
  65. Hetherington, S.; Hughes, A.R.; Mosteller, M.; Shortino, D.; Baker, K.L.; Spreen, W.; Lai, E.; Davies, K.; Handley, A.; Dow, D.J.; et al. Genetic variations in HLA-B region and hypersensitivity reactions to abacavir. Lancet 2002, 359, 1121–1122. [Google Scholar] [CrossRef]
  66. Cheung, C.L.; Sing, C.W.; Tang, C.S.M.; Cheng, V.K.F.; Pirmohamed, M.; Choi, C.H.; Hung, C.S.; Lau, E.Y.F.; Lee, K.F.; Mak, M.W.H.; et al. HLA-B∗38:02:01 predicts carbimazole/methimazole-induced agranulocytosis. Clin. Pharmacol. Ther. 2016, 99, 555–561. [Google Scholar] [CrossRef]
  67. Spraggs, C.F.; Budde, L.R.; Briley, L.P.; Bing, N.; Cox, C.J.; King, K.S.; Whittaker, J.C.; Mooser, V.E.; Preston, A.J.; Stein, S.H.; et al. HLA-DQA1*02:01 is a major risk factor for lapatinib-induced hepatotoxicity in women with advanced breast cancer. J. Clin. Oncol. 2011, 29, 667–673. [Google Scholar] [CrossRef]
  68. Cismaru, A.L.; Grimm, L.; Rudin, D.; Ibañez, L.; Liakoni, E.; Bonadies, N.; Kreutz, R.; Hallberg, P.; Wadelius, M.; Haschke, M.; et al. High-Throughput Sequencing to Investigate Associations Between HLA Genes and Metamizole-Induced Agranulocytosis. Front. Genet. 2020, 11. [Google Scholar] [CrossRef]
  69. Astle, W.J.; Elding, H.; Jiang, T.; Allen, D.; Ruklisa, D.; Mann, A.L.; Mead, D.; Bouman, H.; Riveros-Mckay, F.; Kostadima, M.A.; et al. The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease. Cell 2016, 167, 1415–1429.e19. [Google Scholar] [CrossRef] [Green Version]
  70. Kichaev, G.; Bhatia, G.; Loh, P.R.; Gazal, S.; Burch, K.; Freund, M.K.; Schoech, A.; Pasaniuc, B.; Price, A.L. Leveraging Polygenic Functional Enrichment to Improve GWAS Power. Am. J. Hum. Genet. 2019, 104, 65–75. [Google Scholar] [CrossRef] [Green Version]
  71. Canela-Xandri, O.; Rawlik, K.; Tenesa, A. An atlas of genetic associations in UK Biobank. Nat. Genet. 2018, 50, 1593–1599. [Google Scholar] [CrossRef] [PubMed]
  72. Liang, C.; Wang, X.; Hu, J.; Lian, X.; Zhu, T.; Zhang, H.; Gu, N.; Li, J. PTPRO Promotes Oxidized Low-Density Lipoprotein Induced Oxidative Stress and Cell Apoptosis through Toll-Like Receptor 4/Nuclear Factor κb Pathway. Cell. Physiol. Biochem. 2017, 42, 495–505. [Google Scholar] [CrossRef] [PubMed]
  73. Collet, J.F.; Messens, J. Structure, function, and mechanism of thioredoxin proteins. Antioxid. Redox Signal. 2010, 13, 1205–1216. [Google Scholar] [CrossRef] [PubMed]
  74. Mitsui, A.; Hamuro, J.; Nakamura, H.; Kondo, N.; Hirabayashi, Y.; Ishizaki-Koizumi, S.; Hirakawa, T.; Inoue, T.; Yodoi, J. Overexpression of human thioredoxin in transgenic mice controls oxidative stress and life span. Antioxid. Redox Signal. 2002, 4, 693–696. [Google Scholar] [CrossRef] [PubMed]
  75. Lillig, C.H.; Holmgren, A. Thioredoxin and related molecules-From biology to health and disease. Antioxid. Redox Signal. 2007, 9, 25–47. [Google Scholar] [CrossRef]
  76. Uhlén, M.; Fagerberg, L.; Hallström, B.M.; Lindskog, C.; Oksvold, P.; Mardinoglu, A.; Sivertsson, Å.; Kampf, C.; Sjöstedt, E.; Asplund, A.; et al. Tissue-based map of the human proteome. Science 2015, 347, 1260419. [Google Scholar] [CrossRef]
  77. Uhlen, M.; Karlsson, M.J.; Zhong, W.; Tebani, A.; Pou, C.; Mikes, J.; Lakshmikanth, T.; Forsström, B.; Edfors, F.; Odeberg, J.; et al. A genome-wide transcriptomic analysis of protein-coding genes in human blood cells. Science 2019, 366, eaax9198. [Google Scholar] [CrossRef]
  78. Matsuzawa, A. Thioredoxin and redox signaling: Roles of the thioredoxin system in control of cell fate. Arch. Biochem. Biophys. 2017, 617, 101–105. [Google Scholar] [CrossRef]
  79. Lu, J.; Holmgren, A. The thioredoxin antioxidant system. Free Radic. Biol. Med. 2014, 66, 75–87. [Google Scholar] [CrossRef]
  80. Lu, J.; Holmgren, A. Thioredoxin system in cell death progression. Antioxid. Redox Signal. 2012, 17, 1738–1747. [Google Scholar] [CrossRef]
  81. Ravi, D.; Muniyappa, H.; Das, K.C. Endogenous thioredoxin is required for redox cycling of anthracyclines and p53-dependent apoptosis in cancer cells. J. Biol. Chem. 2005, 280, 40084–40096. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Periasamy, P.; Tran, V.; O’Neill, H.C. Identification of genes which regulate stromadependent in vitro hematopoiesis. PLoS ONE 2018, 13, e0205583. [Google Scholar] [CrossRef] [PubMed]
  83. Sato-Nishiuchi, R.; Nakano, I.; Ozawa, A.; Sato, Y.; Takeichi, M.; Kiyozumi, D.; Yamazaki, K.; Yasunaga, T.; Futaki, S.; Sekiguchi, K. Polydom/SVEP1 is a ligand for integrin α9β1. J. Biol. Chem. 2012, 287, 25615–25630. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Schreiber, T.D.; Steinl, C.; Essl, M.; Abele, H.; Geiger, K.; Müller, C.A.; Aicher, W.K.; Klein, G. The integrin α9β1 on hematopoietic stem and progenitor cells: Involvement in cell adhesion, proliferation and differentiation. Haematologica 2009, 94, 1493–1501. [Google Scholar] [CrossRef]
  85. Chen, C.; Huang, X.; Atakilit, A.; Zhu, Q.S.; Corey, S.J.; Sheppard, D. The Integrin α9β1 Contributes to Granulopoiesis by Enhancing Granulocyte Colony-Stimulating Factor Receptor Signaling. Immunity 2006, 25, 895–906. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Zhang, Y.; Johansson, E.; Miller, M.L.; Jänicke, R.U.; Ferguson, D.J.; Plas, D.; Meller, J.; Anderson, M.W. Identification of a conserved anti-apoptotic protein that modulates the mitochondrial apoptosis pathway. PLoS ONE 2011, 6, e25284. [Google Scholar] [CrossRef]
  87. Wright, G.E.B.; Amstutz, U.; Drögemöller, B.I.; Shih, J.; Rassekh, S.R.; Hayden, M.R.; Carleton, B.C.; Ross, C.J.D.; Visscher, H.; Aminkeng, F.; et al. Pharmacogenomics of Vincristine-Induced Peripheral Neuropathy Implicates Pharmacokinetic and Inherited Neuropathy Genes. Clin. Pharmacol. Ther. 2019, 105, 402–410. [Google Scholar] [CrossRef] [Green Version]
  88. Usami, Y.; Wu, Y.; Göttlinger, H.G. SERINC3 and SERINC5 restrict HIV-1 infectivity and are counteracted by Nef. Nature 2015, 526, 218–223. [Google Scholar] [CrossRef] [Green Version]
  89. Regev, A.; Teichmann, S.A.; Lander, E.S.; Amit, I.; Benoist, C.; Birney, E.; Bodenmiller, B.; Campbell, P.; Carninci, P.; Clatworthy, M.; et al. The human cell atlas. Elife 2017, 6. [Google Scholar] [CrossRef]
  90. Popescu, D.M.; Botting, R.A.; Stephenson, E.; Green, K.; Webb, S.; Jardine, L.; Calderbank, E.F.; Polanski, K.; Goh, I.; Efremova, M.; et al. Decoding human fetal liver haematopoiesis. Nature 2019, 574, 365–371. [Google Scholar] [CrossRef]
  91. Setty, M.; Kiseliovas, V.; Levine, J.; Gayoso, A.; Mazutis, L.; Pe’er, D. Palantir characterizes cell fate continuities in human hematopoiesis. bioRxiv 2018, 385328. [Google Scholar] [CrossRef]
  92. Daly, A.K. Using genome-wide association studies to identify genes important in serious adverse drug reactions. Annu. Rev. Pharmacol. Toxicol. 2012, 52, 21–35. [Google Scholar] [CrossRef] [PubMed]
  93. Diouf, B.; Crews, K.R.; Lew, G.; Pei, D.; Cheng, C.; Bao, J.; Zheng, J.J.; Yang, W.; Fan, Y.; Wheeler, H.E.; et al. Association of an inherited genetic variant with vincristine-related peripheral neuropathy in children with acute lymphoblastic leukemia. JAMA J. Am. Med. Assoc. 2015, 313, 815–823. [Google Scholar] [CrossRef]
  94. Aminkeng, F.; Bhavsar, A.P.; Visscher, H.; Rassekh, S.R.; Li, Y.; Lee, J.W.; Brunham, L.R.; Caron, H.N.; Van Dalen, E.C.; Kremer, L.C.; et al. A coding variant in RARG confers susceptibility to anthracycline-induced cardiotoxicity in childhood cancer. Nat. Genet. 2015, 47, 1079–1084. [Google Scholar] [CrossRef] [PubMed]
  95. Kowalec, K.; Wright, G.E.B.; Drögemöller, B.I.; Aminkeng, F.; Bhavsar, A.P.; Kingwell, E.; Yoshida, E.M.; Traboulsee, A.; Marrie, R.A.; Kremenchutzky, M.; et al. Common variation near IRF6 is associated with IFN-β-induced liver injury in multiple sclerosis. Nat. Genet. 2018, 50, 1081–1085. [Google Scholar] [CrossRef] [PubMed]
  96. Nelson, M.R.; Bacanu, S.A.; Mosteller, M.; Li, L.; Bowman, C.E.; Roses, A.D.; Lai, E.H.; Ehm, M.G. Genome-wide approaches to identify pharmacogenetic contributions to adverse drug reactions. Pharm. J. 2009, 9, 23–33. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Study design. Cases and controls were genotyped separately using the indictaed Illumina arrays. MIA/MIN-CH = Swiss cohort, EuDAC-ES = Spanish cohort, EuDAC-DE = German cohort, SNP = single nucleotide polymorphisms, K = 1000, M = 100,000.
Figure 1. Study design. Cases and controls were genotyped separately using the indictaed Illumina arrays. MIA/MIN-CH = Swiss cohort, EuDAC-ES = Spanish cohort, EuDAC-DE = German cohort, SNP = single nucleotide polymorphisms, K = 1000, M = 100,000.
Genes 11 01275 g001
Figure 2. Manhattan plot of genome-wide association analyses in the MIA/MIN-CH cohort. Results from 304,704 genotyped SNPs after quality control are shown, adjusted by multidimensional scaling (MDS) dimensions 1–4: (a) Analysis of all 45 cases (MIA and MIN) versus 191 controls (tolerant and unexposed combined). The top SNP was rs191786, located in an intergenic region nearest to the TENT5A gene on chromosome 6; (b) Analysis of 30 MIA cases versus 191 tolerant/unexposed controls. The top SNP was rs9366076, located in an intergenic region nearest to the RNASET2 gene on chromosome 6. The red line shows the threshold for genome-wide significance of 1 × 10−7.
Figure 2. Manhattan plot of genome-wide association analyses in the MIA/MIN-CH cohort. Results from 304,704 genotyped SNPs after quality control are shown, adjusted by multidimensional scaling (MDS) dimensions 1–4: (a) Analysis of all 45 cases (MIA and MIN) versus 191 controls (tolerant and unexposed combined). The top SNP was rs191786, located in an intergenic region nearest to the TENT5A gene on chromosome 6; (b) Analysis of 30 MIA cases versus 191 tolerant/unexposed controls. The top SNP was rs9366076, located in an intergenic region nearest to the RNASET2 gene on chromosome 6. The red line shows the threshold for genome-wide significance of 1 × 10−7.
Genes 11 01275 g002
Figure 3. Manhattan plot of GWAS meta-analyses in the three independent cohorts (MIA/MIN-CH, EuDAC-ES and EuDAC-DE). Results from approximately 7 million SNPs after imputation and quality control are shown, adjusted by multidimensional scaling (MDS) dimensions 1–4 and sex (only for EuDAC-ES). The red line shows the threshold for genome-wide significance of 1 × 10−7. (a) Analysis of all 86 cases (MIA and MIN) versus 464 controls (tolerant and unexposed). The top SNP was rs11583606, located in an intron of the TGFBR3 gene on chromosome 1. (b) Analysis of 71 agranulocytosis cases (MIA) versus 464 controls (tolerant and unexposed). The top SNP was rs55898176, located in an intergenic region nearest to the CAAP1 gene on chromosome 9.
Figure 3. Manhattan plot of GWAS meta-analyses in the three independent cohorts (MIA/MIN-CH, EuDAC-ES and EuDAC-DE). Results from approximately 7 million SNPs after imputation and quality control are shown, adjusted by multidimensional scaling (MDS) dimensions 1–4 and sex (only for EuDAC-ES). The red line shows the threshold for genome-wide significance of 1 × 10−7. (a) Analysis of all 86 cases (MIA and MIN) versus 464 controls (tolerant and unexposed). The top SNP was rs11583606, located in an intron of the TGFBR3 gene on chromosome 1. (b) Analysis of 71 agranulocytosis cases (MIA) versus 464 controls (tolerant and unexposed). The top SNP was rs55898176, located in an intergenic region nearest to the CAAP1 gene on chromosome 9.
Genes 11 01275 g003aGenes 11 01275 g003b
Figure 4. Regional association plots for regions around rs55898176 and rs4427239. Plots were produced in Locus Zoom and show the most strongly associated independent SNPs: rs55898176 (a) and rs4427239 (b) from the GWAS meta-analyses of the 71 MIA cases versus 464 controls (tolerant and unexposed). Different colors represent the strength of the LD of each SNP with the most significant SNP represented by a diamond.
Figure 4. Regional association plots for regions around rs55898176 and rs4427239. Plots were produced in Locus Zoom and show the most strongly associated independent SNPs: rs55898176 (a) and rs4427239 (b) from the GWAS meta-analyses of the 71 MIA cases versus 464 controls (tolerant and unexposed). Different colors represent the strength of the LD of each SNP with the most significant SNP represented by a diamond.
Genes 11 01275 g004
Table 1. Characteristics of cases and controls in the three independent cohorts. Data are N(%) or median (range). a cases, b controls, * latency missing for 3 MIA cases and 1 MIN case. ANC = lowest absolute neutrophil count, BMI = body mass index, NA = not applicable/not available.
Table 1. Characteristics of cases and controls in the three independent cohorts. Data are N(%) or median (range). a cases, b controls, * latency missing for 3 MIA cases and 1 MIN case. ANC = lowest absolute neutrophil count, BMI = body mass index, NA = not applicable/not available.
CohortMIA/MIN-CHEuDAC-DEEuDAC-ES
CasesControlsCasesControlsCasesControls
N = 45N = 191N = 12N = 92N = 29N = 181
ANC < 500/uL30 (67)-12-29-
Sex, male (%)13 (42)17 (45)4 (33)41 (44.5)6 (21)87 (48)
Age, years (%)
<257 (23)1 (3)2 (16.6)NA3 (10.3)NA
25–4411 (35)6 (16)6 (50)NA6 (21)NA
45–6410 (32)15 (39)2 (16.6)NA12 (41.4)NA
65–743 (10)9 (24)2 (16.6)NA5 (17)NA
>74-7 (18)-NA3 (10.3)NA
BMI, median (range)24 (19–47)28 (16–39)NANANANA
Latency time * a/treatment duration b, days17 (1–204)25 (1–5297)33.5 (4–9855)NA11.5 (1–235)NA
Table 2. Top associations of the GWAS meta-analysis with MIA/MIN in all cases versus controls in the three independent European cohorts. Top GWAS meta-analysis results based on approximately 7 million SNPs after imputation in 83 MIA/MIN cases vs. all 464 controls (tolerant and unexposed). All results were adjusted for four genetic multidimensional scaling (MDS) components and sex (only for EuDAC-ES). Chromosomal location is according to the Genome Reference Consortium human assembly GRCh37. CHR = chromosome, SNP = single nucleotide polymorphism, BP = base pair, MAF = minor allele frequency, OR [95%] = odds ratio with 95% confidence interval, HetISq = I2 statistic which measures heterogeneity on a scale of 0–100%.
Table 2. Top associations of the GWAS meta-analysis with MIA/MIN in all cases versus controls in the three independent European cohorts. Top GWAS meta-analysis results based on approximately 7 million SNPs after imputation in 83 MIA/MIN cases vs. all 464 controls (tolerant and unexposed). All results were adjusted for four genetic multidimensional scaling (MDS) components and sex (only for EuDAC-ES). Chromosomal location is according to the Genome Reference Consortium human assembly GRCh37. CHR = chromosome, SNP = single nucleotide polymorphism, BP = base pair, MAF = minor allele frequency, OR [95%] = odds ratio with 95% confidence interval, HetISq = I2 statistic which measures heterogeneity on a scale of 0–100%.
CHRSNPAlleles (Minor/Major)BPMAF Cases
MIA/MIN-CH| EuDAC-ES| EuDAC-DE
MAF Controls
MIA/MIN-CH| EuDAC-ES| EuDAC-DE
OR [95%]p-ValueGene RegionHetISq
1rs11583606T/C923492470.10| 014| 0.0830.023| 0.025| 0.0277.0 [3.37–14.5]1.72 × 10−7TGFBR30
1rs149072800C/T924457200.089| 0.12| 0.0830.020| 0.019| 0.0227.81 [3.57–17.1]2.66 × 10−7BRDT0
1rs146378328G/A925280470.089| 0.10| 0.0830.018| 0.019| 0.0218.31 [3.69–18.7]2.93 × 10−7EPHX40
1rs75499485G/A924862740.089| 0.10| 0.0830.020| 0.019| 0.0217.96 [3.56–17.8]4.23 × 10−7EPHX40
12rs112917452C/A156388580.067| 0.15| 0.0420.016| 0.027| 0.0167.30 [3.29–16.20]9.92 × 10−7PTPRO0
12rs118135416A/G156389140.067| 0.15| 0.0420.016| 0.027| 0.0167.30 [3.29–16.20]9.92 × 10−7PTPRO0
12rs7135120T/C156269200.067| 0.15| 0.0420.016| 0.027| 0.0167.30 [3.29–16.20]9.92 × 10−7PTPRO0
12rs143843248T/C156338120.067| 0.15| 0.0410.016| 0.027| 0.0167.30 [3.29–16.20]9.92 × 10−7PTPRO0
13rs73163933A/G339680200.12| 0.086| 0.210.036| 0.027| 0.0545.36 [2.73–10.5]1.01 × 10−6STARD130
1rs78201766G/A923790780.10| 0.14| 0.0830.026| 0.030| 0.0275.65 [2.81–11.34]1.12 × 10−6TGFBR30
Table 3. Top associations of the GWAS meta-analysis in MIA cases versus controls in the three independent European cohorts. Top GWAS meta-analysis results based on approximately 7 million SNPs after imputation in 71 MIA cases vs. all 464 controls (tolerant and unexposed). All results were adjusted for four genetic multidimensional scaling (MDS) components and sex (only for EuDAC-ES). Chromosomal location is according to the Genome Reference Consortium human assembly GRCh37. CHR = chromosome, SNP = single nucleotide polymorphism, BP = base pair, MAF = minor allele frequency, OR [95%] = odds ratio with 95% confidence interval, HetISq = I2 statistic which measures heterogeneity on scale of 0–100%.
Table 3. Top associations of the GWAS meta-analysis in MIA cases versus controls in the three independent European cohorts. Top GWAS meta-analysis results based on approximately 7 million SNPs after imputation in 71 MIA cases vs. all 464 controls (tolerant and unexposed). All results were adjusted for four genetic multidimensional scaling (MDS) components and sex (only for EuDAC-ES). Chromosomal location is according to the Genome Reference Consortium human assembly GRCh37. CHR = chromosome, SNP = single nucleotide polymorphism, BP = base pair, MAF = minor allele frequency, OR [95%] = odds ratio with 95% confidence interval, HetISq = I2 statistic which measures heterogeneity on scale of 0–100%.
CHRSNPAlleles (Minor/Major)BPMAF Cases
MIA/MIN-CH| EuDAC-ES| EuDAC-DE
MAF Controls
MIA/MIN-CH| EuDAC-ES| EuDAC-DE
OR [95%]p-ValueGene RegionHetISq
9rs55898176T/C267152940.27| 0.24| 0.250.065| 0.12| 0.0814.01 [2.41–6.68]1.01 × 10−7-21.7
9rs112223975C/G267158280.27| 0.24| 0.250.065| 0.12| 0.0813.89 [2.34–6.48]1.50 × 10−7-28.6
9rs11790418G/A267130120.27| 0.24| 0.250.065| 0.12| 0.0983.81 [2.29–6.35]2.54 × 10−7-29.3
9rs1434481G/C267111340.27| 0.24| 0.250.065| 0.11| 0.0983.81 [2.29–6.35]2.54 × 10−7-29.3
9rs28475568G/C267099330.27| 0.24| 0.250.065| 0.11| 0.0983.81 [2.29–6.35]2.54 × 10−7-29.3
9rs28649995A/G267099120.27| 0.24| 0.250.065| 0.11| 0.0983.81 [2.29–6.35]2.54 × 10−7-14.6
9rs56285046A/G267149500.28| 0.24| 0.250.073| 0.12| 0.0923.70 [2.27–6.05]3.25 × 10−7-0
9rs77949268A/G267383660.27| 0.26| 0.250.078| 0.12| 0.103.59 [2.20–5.87]3.74 × 10−7-0
9rs4427239A/G1132706010.16| 0.15| 0.0410.029| 0.041| 0.0165.47 [2.81–10.65]5.75 × 10−7SVEP10
9rs10759436C/T1132686500.15| 0.15| 0.0410.029| 0.041| 0.0165.81 [2.92–11.54]8.13 × 10−7SVEP10
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cismaru, A.L.; Rudin, D.; Ibañez, L.; Liakoni, E.; Bonadies, N.; Kreutz, R.; Carvajal, A.; Lucena, M.I.; Martin, J.; Sancho Ponce, E.; et al. Genome-Wide Association Study of Metamizole-Induced Agranulocytosis in European Populations. Genes 2020, 11, 1275. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11111275

AMA Style

Cismaru AL, Rudin D, Ibañez L, Liakoni E, Bonadies N, Kreutz R, Carvajal A, Lucena MI, Martin J, Sancho Ponce E, et al. Genome-Wide Association Study of Metamizole-Induced Agranulocytosis in European Populations. Genes. 2020; 11(11):1275. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11111275

Chicago/Turabian Style

Cismaru, Anca Liliana, Deborah Rudin, Luisa Ibañez, Evangelia Liakoni, Nicolas Bonadies, Reinhold Kreutz, Alfonso Carvajal, Maria Isabel Lucena, Javier Martin, Esther Sancho Ponce, and et al. 2020. "Genome-Wide Association Study of Metamizole-Induced Agranulocytosis in European Populations" Genes 11, no. 11: 1275. https://0-doi-org.brum.beds.ac.uk/10.3390/genes11111275

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop