Gene-based association analysis of a large patient cohort provides insights into genetics of atypical femur fractures.

Several small genetic association studies have been conducted for atypical femur fracture (AFF) without replication of results. We assessed previously implicated and novel genes associated with AFFs in a larger set of unrelated AFF cases using whole exome sequencing (WES). We performed gene-based association analysis on 139 European AFF cases and 196 controls matched for bisphosphonate use. We tested all rare, protein-altering variants using both candidate gene and hypothesis-free approaches. In the latter, genes suggestively associated with AFFs (uncorrected p-values <.01) were investigated in a Swedish whole-genome sequencing replication study and assessed in 46 non-European cases. In the candidate gene analysis, PLOD2 showed a suggestive signal. The hypothesis-free approach revealed 10 tentative associations, with XRN2, SORD, and PLOD2 being the most likely candidates for AFF. XRN2 and PLOD2 showed consistent direction of effect estimates in the replication analysis, albeit not statistically significant. Three SNPs associated with SORD expression according to the GTEx portal were in linkage disequilibrium (R2 ≥ 0.2) with an SNP previously reported in a genome-wide association study of AFF. The prevalence of carriers of variants for both PLOD2 and SORD was higher in Asian versus European cases. While we did not identify genes enriched for damaging variants, we found suggestive evidence of a role for XRN2, PLOD2, and SORD, which requires further investigation. Our findings indicate that genetic factors responsible for AFFs are not widely shared among AFF cases. The study provides a stepping-stone for future larger genetic studies of AFF.

We investigated the genetic factors contributing to atypical femur fractures (AFF), which are rare and unusual fractures in the thigh bone. These fractures are related to the use of bisphosphonates (BP), which are prescribed to prevent fractures caused by osteoporosis. Previous studies suggested potential genetic links, but their findings were not confirmed in larger groups. To address this, we analyzed genetic data from 139 European individuals with AFF and 196 individuals without AFF, all of whom used BP, using a genetic technique called whole exome sequencing. Our results suggested three genes­XRN2, SORD, and PLOD2­might be linked to AFF, although the evidence was not conclusive. Importantly, our findings suggest that AFF may be caused by different genes in different individuals. A much larger sample size is now needed to fully understand the genetic architecture of AFF. These findings may guide future research into the genetic causes of AFF.

Precision Diagnostics in Myeloid Malignancies: Development and Validation of a National Capture-Based Gene Panel.

Gene panel sequencing has become a common diagnostic tool for detecting somatically acquired mutations in myeloid neoplasms. However, many panels have restricted content, provide insufficient sensitivity levels, or lack clinically validated workflows. We here describe the development and validation of the Genomic Medicine Sweden myeloid gene panel (GMS-MGP), a capture-based 191 gene panel including mandatory genes in contemporary guidelines as well as emerging candidates. The GMS-MGP displayed uniform coverage across all targets, including recognized difficult GC-rich areas. The validation of 117 previously described somatic variants showed a 100% concordance with a limit-of-detection of a 0.5% variant allele frequency (VAF), achieved by utilizing error correction and filtering against a panel-of-normals. A national interlaboratory comparison investigating 56 somatic variants demonstrated highly concordant results in both detection rate and reported VAFs. In addition, prospective analysis of 323 patients analyzed with the GMS-MGP as part of standard-of-care identified clinically significant genes as well as recurrent mutations in less well-studied genes. In conclusion, the GMS-MGP workflow supports sensitive detection of all clinically relevant genes, facilitates novel findings, and is, based on the capture-based design, easy to update once new guidelines become available. The GMS-MGP provides an important step toward nationally harmonized precision diagnostics of myeloid malignancies.

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.

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.

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.

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.

The Clinical Pharmacogenetics Implementation Consortium Guideline for SLCO1B1, ABCG2, and CYP2C9 genotypes and Statin-Associated Musculoskeletal Symptoms.

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

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.

[New recommendation on pharmacogenetic screening prior to 5-fluorouracil based cancer treatment - Swedish experience indicates less adverse effects and healthcare cost savings]. / Fördelaktigt med genetisk analys före behandling med 5-fluorouracil.

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

Association of Variants Near the Bradykinin Receptor B2 Gene With Angioedema in Patients Taking ACE Inhibitors.

BACKGROUND: Angioedema is a rare but potentially life-threatening adverse reaction associated with angiotensin-converting enzyme (ACE) inhibitors. Identification of potential genetic factors related to this adverse event may help identify at-risk patients. OBJECTIVES: The aim of this study was to identify genetic factors associated with ACE inhibitor-associated angioedema. METHODS: A genomewide association study involving patients of European descent, all taking ACE inhibitors, was conducted in a discovery cohort (Copenhagen Hospital Biobank), and associations were confirmed in a replication cohort (Swedegene). Cases were defined as subjects with angioedema events and filled prescriptions for ACE inhibitors ≤180 days before the events. Control subjects were defined as those with continuous treatment with ACE inhibitors without any history of angioedema. Odds ratios (ORs) and 95% confidence intervals (CIs) were computed for angioedema risk using logistic mixed model regression analysis. Summary statistics from the discovery and replication cohorts were analyzed using a fixed-effects meta-analysis model. RESULTS: The discovery cohort consisted of 462 cases and 53,391 ACE inhibitor-treated control subjects. The replication cohort consisted of 142 cases and 1,345 ACE inhibitor-treated control subjects. In the discovery cohort, 1 locus, residing at chromosome 14q32.2, was identified that associated with angioedema at the genomewide significance level of P <5 × 10-8. The lead variant at this locus, rs34485356, is an intergenic variant located 60 kb upstream of BDKRB2 (OR: 1.62; 95% CI: 1.38 to 1.90; P = 4.3 × 10-9). This variant was validated in our replication cohort with a similar direction and effect size (OR: 1.60; 95% CI: 1.13 to 2.25; P = 7.2 × 10-3). We found that carriers of the risk allele had significantly lower systolic (-0.46 mm Hg per T allele; 95% CI: -0.83 to -0.10; P = 0.013) and diastolic (-0.26 mm Hg per T allele; 95% CI: -0.46 to -0.05; P = 0.013) blood pressure. CONCLUSIONS: In this genomewide association study involving individuals treated with ACE inhibitors, we found that common variants located in close proximity to the bradykinin receptor B2 gene were associated with increased risk for ACE inhibitor-related angioedema.

PharmVar GeneFocus: CYP2C9.

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

[Pharmacogenomics - a cornerstone of Precision Medicine. Genomic Medicine Sweden analyses genotypes associated with serious drug toxicity or therapeutic failure]. / Farmakogenomik ­ individuell anpassning av läkemedel och dos.

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

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

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

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

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