Next-generation phenotyping: introducing phecodeX for enhanced discovery research in medical phenomics.

MOTIVATION: Phecodes are widely used and easily adapted phenotypes based on International Classification of Diseases codes. The current version of phecodes (v1.2) was designed primarily to study common/complex diseases diagnosed in adults; however, there are numerous limitations in the codes and their structure. RESULTS: Here, we present phecodeX, an expanded version of phecodes with a revised structure and 1,761 new codes. PhecodeX adds granularity to phenotypes in key disease domains that are under-represented in the current phecode structure-including infectious disease, pregnancy, congenital anomalies, and neonatology-and is a more robust representation of the medical phenome for global use in discovery research. AVAILABILITY AND IMPLEMENTATION: phecodeX is available at https://github.com/PheWAS/phecodeX.

Association between nitric oxide synthase 3 genetic variant and acute kidney injury following pediatric cardiac surgery.

BACKGROUND: Acute kidney injury (AKI) complicates 30% to 50% of cardiac surgeries in pediatric patients. Genetic variants that affect renal blood flow and inflammation have been associated with AKI after cardiac surgery in diverse populations of adults but have not been studied in children. The objective of this study is to test the hypothesis that common candidate genetic variants are associated with AKI following pediatric cardiac surgery. METHODS: This is a retrospective cohort study at a single tertiary referral children's hospital of 2,062 individual patients undergoing surgery for congenital heart disease from September 2007 to July 2020. Pre-specified variants in candidate genes (AGTR1, APOE, IL6, NOS3, and TNF) were chosen. AKI was defined using Kidney Disease: Improving Global Outcomes serum creatinine criteria in the first week following surgery. Outcomes were analyzed by univariate and multivariable analysis of demographic, clinical, and genetic factors. RESULTS: The study population had median age of 6 (interquartile range [IQR], 1-53) months, 759 (37%) of whom met criteria for postoperative AKI. In unadjusted analyses of each genetic variant, only NOS3 (rs2070744) was associated with lower risk for AKI (OR 0.75, 95% CI 0.62-0.9, P = .002). In logistic regression analyses adjusting for body surface area, previously identified genetic syndrome, Society of Thoracic Surgeons-European Association for Cardio-Thoracic Surgery (STAT) score, cardiopulmonary bypass time, and nephrotoxic medication exposure, the NOS3 variant remained protective against AKI (OR 0.7, 95% CI 0.58-0.85, P<.001). CONCLUSIONS: A common variant in NOS3 is associated with decreased incidence of AKI in children undergoing cardiac surgery. Further analysis of the genetic contributions to postoperative AKI may help identify individual risk in the pediatric population.

Genome-Wide Approach to Measure Variant-Based Heritability of Drug Outcome Phenotypes.

Pharmacogenomic studies have successfully identified variants-typically with large effect sizes in drug target and metabolism enzymes-that predict drug outcome phenotypes. However, these variants may account for a limited proportion of phenotype variability attributable to the genome. Using genome-wide common variation, we measured the narrow-sense heritability ( h SNP 2 ) of seven pharmacodynamic and five pharmacokinetic phenotypes across three cardiovascular drugs, two antibiotics, and three immunosuppressants. We used a Bayesian hierarchical mixed model, BayesR, to model the distribution of genome-wide variant effect sizes for each drug phenotype as a mixture of four normal distributions of fixed variance (0, 0.01%, 0.1%, and 1% of the total additive genetic variance). This model allowed us to parse h SNP 2 into bins representing contributions of no-effect, small-effect, moderate-effect, and large-effect variants, respectively. For the 12 phenotypes, a median of 969 (range 235-6,304) unique individuals of European ancestry and a median of 1,201,626 (range 777,427-1,514,275) variants were included in our analyses. The number of variants contributing to h SNP 2 ranged from 2,791 to 5,356 (median 3,347). Estimates for h SNP 2 ranged from 0.05 (angiotensin-converting enzyme inhibitor-induced cough) to 0.59 (gentamicin concentration). Small-effect and moderate-effect variants contributed a majority to h SNP 2 for every phenotype (range 61-95%). We conclude that drug outcome phenotypes are highly polygenic. Thus, larger genome-wide association studies of drug phenotypes are needed both to discover novel variants and to determine how genome-wide approaches may improve clinical prediction of drug outcomes.

Evidence for Pharmacogenomic Effects on Risperidone Outcomes in Pediatrics.

OBJECTIVE: To determine the association between genetic variants reported to affect risperidone and adverse events (AEs) in children and adolescents. METHODS: Individuals aged 18 years or younger with ≥4 weeks of risperidone exposure in a deidentified DNA biobank were included. The primary outcome was AE frequency as a function of genotype. Individuals were classified according to metabolizer status for CYP2D6, CYP3A4, and CYP3A5; wild type, heterozygote, or homozygote for specific single nucleotide variants for DRD2, DRD3, HTR2A, and HTR2C; and wild type versus nonwild type for multiple uncommon variants in ABCG2, ABCB1, and HTR2C. Tests of association of each classification to AEs were performed using a Fisher exact test and logistic regression, and statistically significant classifications were included in a final logistic regression. RESULTS: The final cohort included 257 individuals. AEs were more common in CYP2D6 poor/intermediate metabolizers (PMs/IMs) than normal/rapid/ultrarapid metabolizers (NMs/RMs/UMs) in univariate and multivariate analysis. HTR2A-rs6311 heterozygotes and homozygotes had fewer AEs than wild types in logistic regression but not in univariate analysis. In the final multivariable model adjusting for age, race, sex, and risperidone dose, AEs were associated with CYP2D6 (adjusted odds ratio [AOR] 2.6, 95% CI 1.1-5.5, for PMs/IMs vs. NMs/RMs/UMs) and HTR2A-rs6311 (AOR 0.6, 95% CI 0.4-0.9, for each variant allele), both consistent with previous studies. CONCLUSION: Children and adolescents who are CYP2D6 PMs/IMs may have an increased risk for risperidone AEs. Of the genes and variants studied, only CYP2D6 has consistent association and sufficient data for clinical use, whereas HTR2A-rs6311 has limited data and requires further study.

Characterizing Pharmacogenetic Testing Among Children's Hospitals.

Although pharmacogenetic testing is becoming increasingly common across medical subspecialties, a broad range of utilization and implementation exists across pediatric centers. Large pediatric institutions that routinely use pharmacogenetics in their patient care have published their practices and experiences; however, minimal data exist regarding the full spectrum of pharmacogenetic implementation among children's hospitals. The primary objective of this nationwide survey was to characterize the availability, concerns, and barriers to pharmacogenetic testing in children's hospitals in the Children's Hospital Association. Initial responses identifying a contact person were received from 18 institutions. Of those 18 institutions, 14 responses (11 complete and 3 partial) to a more detailed survey regarding pharmacogenetic practices were received. The majority of respondents were from urban institutions (72%) and held a Doctor of Pharmacy degree (67%). Among all respondents, the three primary barriers to implementing pharmacogenetic testing identified were test reimbursement, test cost, and money. Conversely, the three least concerning barriers were potential for genetic discrimination, sharing results with family members, and availability of tests in certified laboratories. Low-use sites rated several barriers significantly higher than the high-use sites, including knowledge of pharmacogenetics (P = 0.03), pharmacogenetic interpretations (P = 0.04), and pharmacogenetic-based changes to therapy (P = 0.03). In spite of decreasing costs of pharmacogenetic testing, financial barriers are one of the main barriers perceived by pediatric institutions attempting clinical implementation. Low-use sites may also benefit from education/outreach in order to reduce perceived barriers to implementation.

Pharmacogenetics to Predict Adverse Events Associated With Antidepressants.

OBJECTIVES: To determine the association between cytochrome P450 2C19 (CYP2C19) metabolizer status and risk for escitalopram and citalopram, collectively termed (es)citalopram, and sertraline adverse events (AEs) in children. METHODS: In this retrospective cohort study, we used deidentified electronic health records linked to DNA. The cohort included children ≤18 years with ≥2 days of (es)citalopram or ≥7 days of sertraline exposure. The primary outcome was AEs assessed by manual chart review. CYP2C19 was genotyped for functional variants (*2, *3, *4, *6, *8, and *17), and individuals were assigned metabolizer status. Association between AEs and metabolizer status was determined by using Cox regression adjusting for age, race, ethnicity, dose, and concomitant CYP2C19-inhibiting medications. RESULTS: The cohort included 249 sertraline-exposed and 458 (es)citalopram-exposed children, with a median age of 14.2 years (interquartile range 11.2-16.2) and 13.4 years (interquartile range 10.1-15.9), respectively. Sertraline AEs were more common in normal metabolizers (NMs) compared to poor metabolizers (PMs) or intermediate metabolizers (IMs) (hazard ratio [HR] 1.8; 95% confidence interval [CI] 1.01-3.2; P = .047) in unadjusted analysis and after adjustment (HR 1.9; CI 1.04-3.4; P = .04). For (es)citalopram, more AEs were observed in NMs than PMs and IMs without statistically significant differences (unadjusted HR 1.6; CI 0.95-2.6; P = .08; adjusted HR 1.6; CI 0.95-2.6; P = .08). CONCLUSIONS: In contrast to adults, in our pediatric cohort, CYP2C19 NMs experienced increased sertraline AEs than PMs and IMs. (Es)citalopram AEs were not associated with CYP2C19 status in the primary analysis. The mechanism underlying this pediatric-specific finding is unknown but may be related to physiologic differences of adolescence. Further research is required to inform genotype-guided prescribing for these drugs in children.

Acute Kidney Injury During Treatment with Intravenous Acyclovir for Suspected or Confirmed Neonatal Herpes Simplex Virus Infection.

OBJECTIVE: To describe the epidemiology of and risk factors associated with acute kidney injury (AKI) during acyclovir treatment in neonates and infants. STUDY DESIGN: We conducted a multicenter (n = 4), retrospective cohort study of all hospitalized infants age <60 days treated with intravenous acyclovir (≥1 dose) for suspected or confirmed neonatal herpes simplex virus disease from January 2011 to December 2015. Infants with serum creatinine measured both before acyclovir (baseline) and during treatment were included. We classified AKI based on changes in creatinine according to published neonatal AKI criteria and performed Cox regression analysis to evaluate risk factors for AKI during acyclovir treatment. RESULTS: We included 1017 infants. The majority received short courses of acyclovir (median, 5 doses). Fifty-seven infants (5.6%) developed AKI during acyclovir treatment, with an incidence rate of AKI at 11.6 per 1000 acyclovir days. Cox regression analysis identified having confirmed herpes simplex virus disease (OR, 4.35; P = .002), receipt of ≥2 concomitant nephrotoxic medications (OR, 3.07; P = .004), receipt of mechanical ventilation (OR, 5.97; P = .001), and admission to an intensive care unit (OR, 6.02; P = .006) as risk factors for AKI during acyclovir treatment. CONCLUSIONS: Among our cohort of infants exposed to acyclovir, the rate of AKI was low. Sicker infants and those exposed to additional nephrotoxic medications seem to be at greater risk for acyclovir-induced toxicity and warrant closer monitoring.

CYP2C19 Phenotype and Risk of Proton Pump Inhibitor-Associated Infections.

OBJECTIVES: Proton pump inhibitors (PPIs) are often used in pediatrics to treat common gastrointestinal disorders, and there are growing concerns for infectious adverse events. Because CYP2C19 inactivates PPIs, genetic variants that increase CYP2C19 function may decrease PPI exposure and infections. We tested the hypothesis that CYP2C19 metabolizer phenotypes are associated with infection event rates in children exposed to PPIs. METHODS: This retrospective biorepository cohort study included individuals aged 0 to 36 months at the time of PPI exposure. Respiratory tract and gastrointestinal tract infection events were identified by using International Classification of Diseases codes in the year after the first PPI mention. Variants defining CYP2C19 *2, *3, *4, *8, *9, and *17 were genotyped, and all individuals were classified as CYP2C19 poor or intermediate, normal metabolizers (NMs), or rapid or ultrarapid metabolizers (RM/UMs). Infection rates were compared by using univariate and multivariate analyses. RESULTS: In all, 670 individuals were included (median age 7 months; 44% girls). CYP2C19 NMs (n = 267; 40%) had a higher infection rate than RM/UMs (n = 220; 33%; median 2 vs 1 infections per person per year; P = .03). There was no difference between poor or intermediate (n = 183; 27%) and NMs. In multivariable analysis of NMs and RM/UMs adjusting for age, sex, PPI dose, and comorbidities, CYP2C19 metabolizer status remained a significant risk factor for infection events (odds ratio 0.70 [95% confidence interval 0.50-0.97] for RM/UMs versus NMs). CONCLUSIONS: PPI therapy is associated with higher infection rates in children with normal CYP2C19 function than in those with increased CYP2C19 function, highlighting this adverse effect of PPI therapy and the relevance of CYP2C19 genotypes to PPI therapeutic decision-making.

CYP2C9*2 is associated with indomethacin treatment failure for patent ductus arteriosus.

Aims: To identify clinical andgenetic factors associated with indomethacin treatment failure in preterm neonates with patent ductus arteriosus (PDA). Patients & Methods: This is a multicenter cohort study of 144 preterm infants (22-32 weeks gestational age) at three centers who received at least one treatment course of indomethacin for PDA. Indomethacin failure was defined as requiring subsequent surgical intervention. Results: In multivariate analysis, gestational age (AOR 0.76, 95% CI 0.60-0.96), surfactant use (AOR 9.77, 95% CI 1.15-83.26), and CYP2C9*2 (AOR 3.74; 95% CI 1.34-10.44) were each associated with indomethacin failure. Conclusion: Age, surfactant use, and CYP2C9*2 influence indomethacin treatment outcome in preterm infants with PDA. This combination of clinical and genetic factors may facilitate targeted indomethacin use for PDA.

Adverse Events Associated with Risperidone Use in Pediatric Patients: A Retrospective Biobank Study.

BACKGROUND: Although risperidone is increasingly used for behavioral indications in children, the associated adverse events (AEs) are not well defined in this population. OBJECTIVE: We determined the incidence of and risk factors for AEs among children treated with risperidone at our institution, an academic medical center with inpatient, outpatient, generalist, and specialist pediatric care. METHODS: The study included children aged ≤ 18 years with ≥ 4 weeks of risperidone exposure. Data were obtained using de-identified electronic health records. AEs were defined as any untoward event attributed to risperidone reported by the patient, parent/guardian, or physician or detected following a laboratory investigation. Associations between AEs and clinical variables were determined using univariate and multivariate analyses. RESULTS: The study cohort included 371 individuals (median age 7.8 years [interquartile range 5.9-10.2]; 271 [73.0%] male). The two most common primary diagnoses were attention-deficit/hyperactivity disorder (160 [43.1%]) and autism (102 [27.5%]). The most frequent indications for risperidone were aggression (166 [44.7%]) and behavioral problems (114 [30.7%]). Altogether, 110 (29.6%) individuals had 156 AEs. Weight gain (32 [20.5%]) and extrapyramidal symptoms (23 [14.7%]) were the most common AEs. Aggression, irritability, and self-injurious behavior were positively associated with AEs, and concomitant analgesics and antibiotics were negatively associated. In multivariate analysis, associations remained significant for self-injurious behavior (adjusted odds ratio [aOR] 3.1; 95% confidence interval [CI] 1.7-5.4) and concomitant antibiotics (aOR 0.2; 95% CI 0.1-0.9). CONCLUSIONS: Nearly one in three children treated with risperidone for ≥ 1 month experienced one or more AEs. Particular vigilance is warranted for children with self-injurious behavior.

CYP2D6 genotype and adverse events to risperidone in children and adolescents.

BACKGROUND: There are few and conflicting data on the role of cytochrome P450 2D6 (CYP2D6) polymorphisms in relation to risperidone adverse events (AEs) in children. This study assessed the association between CYP2D6 metabolizer status and risk for risperidone AEs in children. METHODS: Children ≤18 years with at least 4 weeks of risperidone exposure were identified using BioVU, a de-identified DNA biobank linked to electronic health record data. The primary outcome of this study was AEs. After DNA sequencing, individuals were classified as CYP2D6 poor, intermediate, normal, or ultrarapid CYP2D6 metabolizers. RESULTS: For analysis, the 257 individuals were grouped as poor/intermediate metabolizers (n = 33, 13%) and normal/ultrarapid metabolizers (n = 224, 87%). AEs were more common in poor/intermediate vs. normal/ultrarapid metabolizers (15/33, 46% vs. 61/224, 27%, P = 0.04). In multivariate analysis adjusting for age, sex, race, and initial dose, poor/intermediate metabolizers had increased AE risk (adjusted odds ratio 2.4, 95% confidence interval 1.1-5.1, P = 0.03). CONCLUSION: Children with CYP2D6 poor or intermediate metabolizer phenotypes are at greater risk for risperidone AEs. Pre-prescription genotyping could identify this high-risk subset for an alternate therapy, risperidone dose reduction, and/or increased monitoring for AEs.

Clinical Pharmacogenetics of Cytochrome P450-Associated Drugs in Children.

Cytochrome P450 (CYP) enzymes are commonly involved in drug metabolism, and genetic variation in the genes encoding CYPs are associated with variable drug response. While genotype-guided therapy has been clinically implemented in adults, these associations are less well established for pediatric patients. In order to understand the frequency of pediatric exposures to drugs with known CYP interactions, we compiled all actionable drug-CYP interactions with a high level of evidence using Clinical Pharmacogenomic Implementation Consortium (CPIC) data and surveyed 10 years of electronic health records (EHR) data for the number of children exposed to CYP-associated drugs. Subsequently, we performed a focused literature review for drugs commonly used in pediatrics, defined as more than 5000 pediatric patients exposed in the decade-long EHR cohort. There were 48 drug-CYP interactions with a high level of evidence in the CPIC database. Of those, only 10 drugs were commonly used in children (ondansetron, oxycodone, codeine, omeprazole, lansoprazole, sertraline, amitriptyline, citalopram, escitalopram, and risperidone). For these drugs, reports of the drug-CYP interaction in cohorts including children were sparse. There are adequate data for implementation of genotype-guided therapy for children for three of the 10 commonly used drugs (codeine, omeprazole and lansoprazole). For the majority of commonly used drugs with known CYP interactions, more data are required to support pharmacogenomic implementation in children.