Utilizing Large Electronic Medical Record Data Sets to Identify Novel Drug-Gene Interactions for Commonly Used Drugs.

Real-world prescribing of drugs differs from the experimental systems, physiological-pharmacokinetic models, and clinical trials used in drug development and licensing, with drugs often used in patients with multiple comorbidities with resultant polypharmacy. The increasing availability of large biobanks linked to electronic healthcare records enables the potential to identify novel drug-gene interactions in large populations of patients. In this study we used three Scottish cohorts and UK Biobank to identify drug-gene interactions for the 50 most commonly used drugs and 162 variants in genes involved in drug pharmacokinetics. We defined two phenotypes based upon prescribing behavior-drug-stop or dose-decrease. Using this approach, we replicate 11 known drug-gene interactions including, for example, CYP2C9/CYP2C8 variants and sulfonylurea/thiazolidinedione prescribing and ABCB1/ABCG2 variants and statin prescribing. We identify eight novel associations after Bonferroni correction, three of which are replicated or validated in the UK Biobank or have other supporting results: The C-allele at rs4918758 in CYP2C9 was associated with a 25% (15-44%) lower odds of dose reduction of quinine, P = 1.6 × 10-5 ; the A-allele at rs9895420 in ABCC3 was associated with a 46% (24-62%) reduction in odds of dose reduction with doxazosin, P = 1.2 × 10-4 , and altered blood pressure response in the UK Biobank; the CYP2D6*2 variant was associated with a 30% (18-40%) reduction in odds of stopping ramipril treatment, P = 1.01 × 10-5 , with similar results seen for enalapril and lisinopril and with other CYP2D6 variants. This study highlights the scope of using large population bioresources linked to medical record data to explore drug-gene interactions at scale.

Association Between Schizophrenia-Related Polygenic Liability and the Occurrence and Level of Mood-Incongruent Psychotic Symptoms in Bipolar Disorder.

Importance: Bipolar disorder (BD) overlaps schizophrenia in its clinical presentation and genetic liability. Alternative approaches to patient stratification beyond current diagnostic categories are needed to understand the underlying disease processes and mechanisms. Objective: To investigate the association between common-variant liability for schizophrenia, indexed by polygenic risk scores (PRSs), and psychotic presentations of BD. Design, Setting, and Participants: This case-control study in the United Kingdom used multinomial logistic regression to estimate differential PRS associations across categories of cases and controls. Participants included in the final analyses were 4436 cases of BD from the Bipolar Disorder Research Network. These cases were compared with the genotypic data for 4976 cases of schizophrenia and 9012 controls from the Type 1 Diabetes Genetics Consortium study and the Generation Scotland study. Data were collected between January 1, 2000, and December 31, 2013. Data analysis was conducted from March 1, 2016, to February 28, 2017. Exposures: Standardized PRSs, calculated using alleles with an association threshold of P < .05 in the second Psychiatric Genomics Consortium genome-wide association study of schizophrenia, were adjusted for the first 10 population principal components and genotyping platforms. Main Outcomes and Measures: Multinomial logit models estimated PRS associations with BD stratified by Research Diagnostic Criteria subtypes of BD, by lifetime occurrence of psychosis, and by lifetime mood-incongruent psychotic features. Ordinal logistic regression examined PRS associations across levels of mood incongruence. Ratings were derived from the Schedules for Clinical Assessment in Neuropsychiatry interview and the Bipolar Affective Disorder Dimension Scale. Results: Of the 4436 cases of BD, 2966 (67%) were female patients, and the mean (SD) age at interview was 46 [12] years. Across clinical phenotypes, there was an exposure-response gradient, with the strongest PRS association for schizophrenia (risk ratio [RR] = 1.94; 95% CI, 1.86-2.01), followed by schizoaffective BD (RR = 1.37; 95% CI, 1.22-1.54), bipolar I disorder subtype (RR = 1.30; 95% CI, 1.24-1.36), and bipolar II disorder subtype (RR = 1.04; 95% CI, 0.97-1.11). Within BD cases, there was an effect gradient, indexed by the nature of psychosis. Prominent mood-incongruent psychotic features had the strongest association (RR = 1.46; 95% CI, 1.36-1.57), followed by mood-congruent psychosis (RR = 1.24; 95% CI, 1.17-1.33) and BD with no history of psychosis (RR = 1.09; 95% CI, 1.04-1.15). Conclusions and Relevance: For the first time to date, a study shows a polygenic-risk gradient across schizophrenia and BD, indexed by the occurrence and level of mood-incongruent psychotic symptoms.