The incidence of narcolepsy rose sharply after the swine influenza A (H1N1) vaccination campaign with Pandemrix. Narcolepsy is an immune-related disorder with excessive daytime sleepiness. The most frequent form is strongly associated with HLA-DQB1*06:02, but only a minority of carriers develop narcolepsy. We aimed to identify genetic markers that predispose to Pandemrix-induced narcolepsy.

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

A report on the 9th Annual Cold Spring Harbor/Wellcome Trust meeting 'Pharmacogenomics and Personalized Medicine', Hinxton, Cambridge, UK, 29 September to 2 October 2011.

We report an extensive study of variability in genes encoding proteins that are believed to be involved in the action and biotransformation of warfarin. Warfarin is a commonly prescribed anticoagulant that is difficult to use because of the wide interindividual variation in dose requirements, the narrow therapeutic range and the risk of serious bleeding. We genotyped 201 patients for polymorphisms in 29 genes in the warfarin interactive pathways and tested them for association with dose requirement. In our study, polymorphisms in or flanking the genes VKORC1, CYP2C9, CYP2C18, CYP2C19, PROC, APOE, EPHX1, CALU, GGCX and ORM1-ORM2 and haplotypes of VKORC1, CYP2C9, CYP2C8, CYP2C19, PROC, F7, GGCX, PROZ, F9, NR1I2 and ORM1-ORM2 were associated with dose (P < 0.05). VKORC1, CYP2C9, CYP2C18 and CYP2C19 were significant after experiment-wise correction for multiple testing (P < 0.000175), however, the association of CYP2C18 and CYP2C19 was fully explained by linkage disequilibrium with CYP2C9*2 and/or *3. PROC and APOE were both significantly associated with dose after correction within each gene. A multiple regression model with VKORC1, CYP2C9, PROC and the non-genetic predictors age, bodyweight, drug interactions and indication for treatment jointly accounted for 62% of variance in warfarin dose. Weaker associations observed for other genes could explain up to approximately 10% additional dose variance, but require testing and validation in an independent and larger data set. Translation of this knowledge into clinical guidelines for warfarin prescription will be likely to have a major impact on the safety and efficacy of warfarin.

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

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.

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

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.

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

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

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

We report the first genome-wide association study (GWAS) whose sample size (1,053 Swedish subjects) is sufficiently powered to detect genome-wide significance (p<1.5 x 10(-7)) for polymorphisms that modestly alter therapeutic warfarin dose. The anticoagulant drug warfarin is widely prescribed for reducing the risk of stroke, thrombosis, pulmonary embolism, and coronary malfunction. However, Caucasians vary widely (20-fold) in the dose needed for therapeutic anticoagulation, and hence prescribed doses may be too low (risking serious illness) or too high (risking severe bleeding). Prior work established that approximately 30% of the dose variance is explained by single nucleotide polymorphisms (SNPs) in the warfarin drug target VKORC1 and another approximately 12% by two non-synonymous SNPs (*2, *3) in the cytochrome P450 warfarin-metabolizing gene CYP2C9. We initially tested each of 325,997 GWAS SNPs for association with warfarin dose by univariate regression and found the strongest statistical signals (p<10(-78)) at SNPs clustering near VKORC1 and the second lowest p-values (p<10(-31)) emanating from CYP2C9. No other SNPs approached genome-wide significance. To enhance detection of weaker effects, we conducted multiple regression adjusting for known influences on warfarin dose (VKORC1, CYP2C9, age, gender) and identified a single SNP (rs2108622) with genome-wide significance (p = 8.3 x 10(-10)) that alters protein coding of the CYP4F2 gene. We confirmed this result in 588 additional Swedish patients (p<0.0029) and, during our investigation, a second group provided independent confirmation from a scan of warfarin-metabolizing genes. We also thoroughly investigated copy number variations, haplotypes, and imputed SNPs, but found no additional highly significant warfarin associations. We present power analysis of our GWAS that is generalizable to other studies, and conclude we had 80% power to detect genome-wide significance for common causative variants or markers explaining at least 1.5% of dose variance. These GWAS results provide further impetus for conducting large-scale trials assessing patient benefit from genotype-based forecasting of warfarin dose.

VKORC1 and CYP2C9 are important contributors to warfarin dose variability, but explain less variability for individuals of African descent than for those of European or Asian descent. We aimed to identify additional variants contributing to warfarin dose requirements in African Americans.

BACKGROUND: Drug-induced agranulocytosis, a severe side effect marked by a deficit or absolute lack of granulocytic white blood cells, is a rare side-effect of the anti-inflammatory drug sulphasalazine. Mutations in the human neutrophil elastase gene (ELA2), causing increased intracellular concentration of this serine protease, inhibits neutrophil differentiation in severe congenital neutropenia (SCN). Since the clinical symptoms of agranulocytosis and SCN are similar, we hypothesized that it may origin from a common genetic variation in ELA2 or that sulphasalazine may affect human neutrophil elastase activity and protein expression. METHODS: We screened for genetic differences in ELA2 in DNA from 36 patients who had suffered from sulphasalazine-induced agranulocytosis, and compared them with 72 patients treated with sulphasalazine without blood reactions. We also performed in vitro studies of the blood cell lines HL60 and U937 after sulphasalazine exposure with respect to cell survival index, neutrophil elastase protein expression and activity. RESULTS: None of the mutations in ELA2, which previously have been reported to be associated with SCN, was found in this material. Protein expression of human neutrophil elastase in lymphoma U937 cells was not affected by treatment with concentrations equivalent to therapeutic doses. Cell survival of lymphoma U937 and promyelocytic leukemia HL-60 cells was not affected in this concentration range, but exhibited a decreased proliferative capacity with higher sulphasalazine concentrations. Interestingly the promyelocytic cells were more sensitive to sulphasalazine than the lymphoma cell line. CONCLUSION: Neutrophil elastase expression and ELA2 mutations do, however, not seem to be involved in the etilogy of sulphasalazine-induced agranulocytosis. Why sulphasalazine is more toxic to promyelocytes than to lymphocytes remains to be explained.

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

Background and Objective: Agranulocytosis is a rare and potentially life-threatening complication of metamizole (dipyrone) intake that is characterized by a loss of circulating neutrophil granulocytes. While the mechanism underlying this adverse drug reaction is not well understood, involvement of the immune system has been suggested. In addition, associations between genetic variants in the Human Leukocyte Antigen (HLA) region and agranulocytosis induced by other drugs have been reported. The aim of the present study was to assess whether genetic variants in classical HLA genes are associated with the susceptibility to metamizole-induced agranulocytosis (MIA) in a European population by targeted resequencing of eight HLA genes. Design: A case-control cohort of Swiss patients with a history of neutropenia or agranulocytosis associated with metamizole exposure (n = 53), metamizole-tolerant (n = 39) and unexposed controls (n = 161) was recruited for this study. A high-throughput resequencing (HTS) and high-resolution typing method was used to sequence and analyze eight HLA loci in a discovery subset of this cohort (n = 31 cases, n = 38 controls). Identified candidate alleles were investigated in the full Swiss cohort as well as in two independent cohorts from Germany and Spain using HLA imputation from genome-wide SNP array data. In addition, variant calling based on HTS data was performed in the discovery subset for the class I genes HLA-A, -B, and -C using the HLA-specific mapper hla-mapper. Results: Eight candidate alleles (p < 0.05) were identified in the discovery subset, of which HLA-C∗04:01 was associated with MIA in the full Swiss cohort (p < 0.01) restricted to agranulocytosis (ANC < 0.5 × 109/L) cases. However, no candidate allele showed a consistent association in the Swiss, German and Spanish cohorts. Analysis of individual sequence variants in class I genes produced consistent results with HLA typing but did not reveal additional small nucleotide variants associated with MIA. Conclusion: Our results do not support an HLA-restricted T cell-mediated immune mechanism for MIA. However, we established an efficient high-resolution (three-field) eight-locus HTS HLA resequencing method to interrogate the HLA region and demonstrated the feasibility of its application to pharmacogenetic studies.

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

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

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.