Citations

Despite advances in pharmacogenomics that have led to innovative therapeutic strategies, substantial interindividual variability in drug efficacy and toxicity persists. Both endogenous and exogenous factors contribute to this unexplained variation, resulting in suboptimal treatment effectiveness, unnecessary healthcare costs, preventable mortality, and significant human suffering. The emerging field of exposomics is aimed at addressing the complex exposures that impact human health and offers an opportunity to interrogate how these exposures influence drug actions and outcomes. Here, we introduce pharmacoexposomics as a complementary framework that systematically quantifies the body's metabolic response to drug and environmental exposures using high-resolution mass spectrometry. We then present IndiPHARM (Individualized Pharmacology), a multi-institution, Advanced Research Projects Agency for Health–supported platform designed to integrate population-scale pharmacoexposomics with pharmacogenomics and clinical data to elucidate interindividual variability in drug response, identify adverse reactions, and inform tailored, effective treatment strategies that advance precision medicine.

Introduction Semaglutide, a glucagon-like peptide-1 receptor agonist (GLP-1RA) used to treat type 2 diabetes mellitus (T2D), has potential associations with higher rates of diabetic retinopathy (DR) complications including proliferative DR (PDR) and diabetic macular edema (DME). The purpose of this study was to determine whether an association exists between semaglutide and PDR and treatment-requiring DR/DME. Research design and methods This was a retrospective cohort study of 14 databases (six administrative claims and eight electronic health records) in the Observational Health Data Sciences and Informatics Evidence Network. Adults with T2D on semaglutide, other GLP-1RA (dulaglutide, exenatide), or non-GLP-1RA medications (empagliflozin, sitagliptin, glipizide) from 1 December 2017 to 31 December 2023 were included. The association between semaglutide and PDR or treatment-requiring DR/DME was assessed using an active-comparator cohort design comparing new users of semaglutide as second-line T2D treatment to those on other GLP-1RAs and non-GLP-1RAs. Propensity score-adjusted Cox proportional hazards models were used to estimate hazard ratios (HRs). Network-wide HR estimates were generated using a random-effects meta-analysis. Results The study included 810 390 new semaglutide users for T2D. PDR risk for semaglutide was similar to dulaglutide (HR 0.81, 95% CI 0.42 to 1.54, p=0.51), empagliflozin (HR 0.83, 95% CI 0.53 to 1.30, p=0.41) and sitagliptin (HR 0.83, 95% CI 0.45 to 1.55, p=0.57) but was lower than glipizide (HR 0.59, 95% CI 0.39 to 0.88, p=0.01). The risk for treatment-requiring DR/DME for semaglutide was similar to empagliflozin (HR 0.66, 95% CI 0.43 to 1.02, p=0.06) but lower than dulaglutide (HR 0.53, 95% CI 0.31 to 0.91, p=0.02), sitagliptin (HR 0.46, 95% CI 0.26 to 0.81, p=0.008) and glipizide (HR 0.55, 95% CI 0.33 to 0.91, p=0.02). Conclusions and relevance We did not identify increased risk for either PDR or treatment-requiring DR/DME comparing semaglutide with other GLP-1RAs or non-GLP-1RAs. Patients with T2D should still undergo close eye care follow-up, particularly when initiating new antihyperglycemic medications.

Importance Semaglutide, a glucagonlike peptide-1 receptor agonist (GLP-1RA), has recently been implicated in cases of nonarteritic anterior ischemic optic neuropathy (NAION), raising safety concerns in the treatment of type 2 diabetes (T2D). Objective To investigate the potential association between semaglutide and NAION in the Observational Health Data Sciences and Informatics (OHDSI) network. Design, Setting, and Participants This was a retrospective study across 14 databases (6 administrative claims and 8 electronic health records). Included were adults with T2D taking semaglutide, other GLP-1RA (dulaglutide, exenatide), or non–GLP-1RA medications (empagliflozin, sitagliptin, glipizide) from December 1, 2017, to December 31, 2023. The incidence proportion and rate of NAION were calculated. Association between semaglutide and NAION was assessed using 2 approaches: an active-comparator cohort design comparing new users of semaglutide with those taking other GLP-1RAs and non–GLP-1RA drugs, and a self-controlled case-series (SCCS) analysis to compare individuals’ risks during exposure and nonexposure periods for each drug. The cohort design used propensity score–adjusted Cox proportional hazards models to estimate hazard ratios (HRs). The SCCS used conditional Poisson regression models to estimate incidence rate ratios (IRRs). Network-wide HR and IRR estimates were generated using a random-effects meta-analysis model. Exposures GLP-1RA and non–GLP-1RAs. Main Outcomes and Measures NAION under 2 alternative definitions based on diagnosis codes: one more inclusive and sensitive, the other more restrictive and specific. Results The study included 37.1 million individuals with T2D, including 810 390 new semaglutide users. Of the 43 620 new users of semaglutide in the Optum’s deidentified Clinformatics Data Mart Database, 24 473 (56%) were aged 50 to 69 years, and 26 699 (61%) were female. The incidence rate of NAION was 14.5 per 100 000 person-years among semaglutide users. The HR for NAION among new users of semaglutide was not different compared with that of the non–GLP-1RAs using the sensitive NAION definition—empagliflozin (HR, 1.44; 95% CI, 0.78-2.68; P = .12), sitagliptin (HR, 1.30; 95% CI, 0.56-3.01; P = .27), and glipizide (HR, 1.23; 95% CI, 0.66-2.28; P = .25). The risk was higher only compared with patients taking empagliflozin (HR, 2.27; 95% CI, 1.16-4.46; P = .02) using the specific definition. SCCS analysis of semaglutide exposure showed an increased risk of NAION (meta-analysis IRR, 1.32; 95% CI, 1.14-1.54; P < .001). Conclusions and Relevance Results of this study suggest a modest increase in the risk of NAION among individuals with T2D associated with semaglutide use, smaller than that previously reported, and warranting further investigation into the clinical implications of this association.

Abstract Importance Chronic kidney disease (CKD) is a genetically complex disease determined by an interplay of monogenic, polygenic, and environmental risks. The most common forms of monogenic kidney disorders include autosomal dominant polycystic kidney disease (ADPKD), caused by mutations in the PKD1 or PKD2 genes, and COL4A-associated nephropathy (COL4A-AN), caused by mutations in COL4A3, COL4A4, or COL4A5 genes. These disorders have incomplete penetrance and variable expressivity. It is presently unknown if some of the variability in penetrance can be attributed to polygenic factors. Objective Our objectives were to investigate the interplay of polygenic and monogenic risk of kidney disease and to test the utility of polygenic risk scores for risk stratification in ADPKD and COL4-AN. Design We combined exome/genome sequencing, SNP micro-array, and electronic health record data from the UK Biobank (N = 469,835 participants) and the All of Us (N = 98,622 participants) datasets to test the effects of the genome-wide polygenic score (GPS) for CKD in ADPKD and COL4-AN variant carriers defined by strict variant classification criteria. We used the eMERGE-III electronic CKD phenotype to define cases (estimated glomerular filtration rate (eGFR) < 60 mL/min/1.73m2 or kidney failure) and controls (eGFR > 90 mL/min/1.73m2 in the absence of kidney disease diagnoses). The GPS was tested as a predictor of CKD in qualifying variant carriers and non-carriers using logistic regression after adjustment for age, sex, diabetes, and genetic ancestry. Results The GPS was predictive of CKD in ADPKD variant carriers (ORmeta=2.28 per SD, 95%CI: 1.55–3.37, P = 2.6E-05) as well as non-carriers (ORmeta=1.72 per SD, 95% CI = 1.69–1.76, P < E-300) independent of age, sex, diabetes, and genetic ancestry. Compared to the middle tertile of the GPS distribution for non-carriers, ADPKD variant carriers in the top tertile had a 54-fold increased risk of CKD, while ADPKD variant carriers in the bottom tertile had only a 3-fold increased risk of CKD. Similarly, the GPS was predictive of CKD in both COL4-AN carriers (ORmeta=1.78, 95% CI = 1.22–2.58, P = 2.38E-03) and non-carriers (ORmeta=1.70, 95%CI: 1.68–1.73 P < E-300). The carriers in the top GPS tertile had a 2.5-fold higher risk of CKD while the risk for carriers in the bottom GPS tertile was not different from the average population risk. Conclusions and Relevance Variable penetrance of kidney disease in ADPKD and COL4-AN is partially explained by the differences in polygenic risk profiles. Accounting for polygenic factors improves risk stratification in monogenic kidney disease and may have clinical implications for genetic counseling.

AbstractObjectiveThe aim of this study was to analyze a publicly available sample of rule-based phenotype definitions to characterize and evaluate the variability of logical constructs used.Materials and MethodsA sample of 33 preexisting phenotype definitions used in research that are represented using Fast Healthcare Interoperability Resources and Clinical Quality Language (CQL) was analyzed using automated analysis of the computable representation of the CQL libraries.ResultsMost of the phenotype definitions include narrative descriptions and flowcharts, while few provide pseudocode or executable artifacts. Most use 4 or fewer medical terminologies. The number of codes used ranges from 5 to 6865, and value sets from 1 to 19. We found that the most common expressions used were literal, data, and logical expressions. Aggregate and arithmetic expressions are the least common. Expression depth ranges from 4 to 27.DiscussionDespite the range of conditions, we found that all of the phenotype definitions consisted of logical criteria, representing both clinical and operational logic, and tabular data, consisting of codes from standard terminologies and keywords for natural language processing. The total number and variety of expressions are low, which may be to simplify implementation, or authors may limit complexity due to data availability constraints.ConclusionsThe phenotype definitions analyzed show significant variation in specific logical, arithmetic, and other operators but are all composed of the same high-level components, namely tabular data and logical expressions. A standard representation for phenotype definitions should support these formats and be modular to support localization and shared logic.

Background: Sequencing Mendelian arrhythmia genes in individuals without an indication for arrhythmia genetic testing can identify carriers of pathogenic or likely pathogenic (P/LP) variants. However, the extent to which these variants are associated with clinically meaningful phenotypes before or after return of variant results is unclear. In addition, the majority of discovered variants are currently classified as variants of uncertain significance, limiting clinical actionability. Methods: The eMERGE-III study (Electronic Medical Records and Genomics Phase III) is a multicenter prospective cohort that included 21 846 participants without previous indication for cardiac genetic testing. Participants were sequenced for 109 Mendelian disease genes, including 10 linked to arrhythmia syndromes. Variant carriers were assessed with electronic health record–derived phenotypes and follow-up clinical examination. Selected variants of uncertain significance (n=50) were characterized in vitro with automated electrophysiology experiments in HEK293 cells. Results: As previously reported, 3.0% of participants had P/LP variants in the 109 genes. Herein, we report 120 participants (0.6%) with P/LP arrhythmia variants. Compared with noncarriers, arrhythmia P/LP carriers had a significantly higher burden of arrhythmia phenotypes in their electronic health records. Fifty-four participants had variant results returned. Nineteen of these 54 participants had inherited arrhythmia syndrome diagnoses (primarily long-QT syndrome), and 12 of these 19 diagnoses were made only after variant results were returned (0.05%). After in vitro functional evaluation of 50 variants of uncertain significance, we reclassified 11 variants: 3 to likely benign and 8 to P/LP. Conclusions: Genome sequencing in a large population without indication for arrhythmia genetic testing identified phenotype-positive carriers of variants in congenital arrhythmia syndrome disease genes. As the genomes of large numbers of people are sequenced, the disease risk from rare variants in arrhythmia genes can be assessed by integrating genomic screening, electronic health record phenotypes, and in vitro functional studies. Registration: URL: https://www.clinicaltrials.gov ; Unique identifier; NCT03394859.

Abstract IntroductionCurrently, one of the commonly used methods for disseminating electronic health record (EHR)-based phenotype algorithms is providing a narrative description of the algorithm logic, often accompanied by flowcharts. A challenge with this mode of dissemination is the potential for under-specification in the algorithm definition, which leads to ambiguity and vagueness. MethodsThis study examines incidents of under-specification that occurred during the implementation of 34 narrative phenotyping algorithms in the electronic Medical Record and Genomics (eMERGE) network. We reviewed the online communication history between algorithm developers and implementers within the Phenotype Knowledge Base (PheKB) platform, where questions could be raised and answered regarding the intended implementation of a phenotype algorithm. ResultsWe developed a taxonomy of under-specification categories via an iterative review process between two groups of annotators. Under-specifications that lead to ambiguity and vagueness were consistently found across narrative phenotype algorithms developed by all involved eMERGE sites. Discussion & ConclusionOur findings highlight that under-specification is an impediment to the accuracy and efficiency of the implementation of current narrative phenotyping algorithms, and we propose approaches for mitigating these issues and improved methods for disseminating EHR phenotyping algorithms.

Background: Coronary heart disease (CHD) is a leading cause of death globally. Although therapy with statins decreases circulating levels of low-density lipoprotein cholesterol and the incidence of CHD, additional events occur despite statin therapy in some individuals. The genetic determinants of this residual cardiovascular risk remain unknown. Methods: We performed a 2-stage genome-wide association study of CHD events during statin therapy. We first identified 3099 cases who experienced CHD events (defined as acute myocardial infarction or the need for coronary revascularization) during statin therapy and 7681 controls without CHD events during comparable intensity and duration of statin therapy from 4 sites in the Electronic Medical Records and Genomics Network. We then sought replication of candidate variants in another 160 cases and 1112 controls from a fifth Electronic Medical Records and Genomics site, which joined the network after the initial genome-wide association study. Finally, we performed a phenome-wide association study for other traits linked to the most significant locus. Results: The meta-analysis identified 7 single nucleotide polymorphisms at a genome-wide level of significance within the LPA/PLG locus associated with CHD events on statin treatment. The most significant association was for an intronic single nucleotide polymorphism within LPA/PLG (rs10455872; minor allele frequency, 0.069; odds ratio, 1.58; 95% confidence interval, 1.35–1.86; P =2.6×10 − 10 ). In the replication cohort, rs10455872 was also associated with CHD events (odds ratio, 1.71; 95% confidence interval, 1.14–2.57; P =0.009). The association of this single nucleotide polymorphism with CHD events was independent of statin-induced change in low-density lipoprotein cholesterol (odds ratio, 1.62; 95% confidence interval, 1.17–2.24; P =0.004) and persisted in individuals with low-density lipoprotein cholesterol ≤70 mg/dL (odds ratio, 2.43; 95% confidence interval, 1.18–4.75; P =0.015). A phenome-wide association study supported the effect of this region on coronary heart disease and did not identify noncardiovascular phenotypes. Conclusions: Genetic variations at the LPA locus are associated with CHD events during statin therapy independently of the extent of low-density lipoprotein cholesterol lowering. This finding provides support for exploring strategies targeting circulating concentrations of lipoprotein(a) to reduce CHD events in patients receiving statins.