The coupling of electronic medical records (EMR) with genetic data has created the potential for implementing reverse genetic approaches in humans, whereby the function of a gene is inferred from the shared pattern of morbidity among homozygotes of a genetic variant. We explored the feasibility of this approach to identify phenotypes associated with low frequency variants using Vanderbilt's EMR-based BioVU resource. We analyzed 1,658 low frequency non-synonymous SNPs (nsSNPs) with a minor allele frequency (MAF)<10% collected on 8,546 subjects. For each nsSNP, we identified diagnoses shared by at least 2 minor allele homozygotes and with an association p<0.05. The diagnoses were reviewed by a clinician to ascertain whether they may share a common mechanistic basis. While a number of biologically compelling clinical patterns of association were observed, the frequency of these associations was identical to that observed using genotype-permuted data sets, indicating that the associations were likely due to chance. To refine our analysis associations, we then restricted the analysis to 711 nsSNPs in genes with phenotypes in the On-line Mendelian Inheritance in Man (OMIM) or knock-out mouse phenotype databases. An initial comparison of the EMR diagnoses to the known in vivo functions of the gene identified 25 candidate nsSNPs, 19 of which had significant genotype-phenotype associations when tested using matched controls. Twleve of the 19 nsSNPs associations were confirmed by a detailed record review. Four of 12 nsSNP-phenotype associations were successfully replicated in an independent data set: thrombosis (F5,rs6031), seizures/convulsions (GPR98,rs13157270), macular degeneration (CNGB3,rs3735972), and GI bleeding (HGFAC,rs16844401). These analyses demonstrate the feasibility and challenges of using reverse genetics approaches to identify novel gene-phenotype associations in human subjects using low frequency variants. As increasing amounts of rare variant data are generated from modern genotyping and sequence platforms, model organism data may be an important tool to enable discovery.

Although risperidone is increasingly used for behavioral indications in children, the associated adverse events (AEs) are not well defined in this population.

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.

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.

Medication-related osteonecrosis of the jaw (MRONJ) is a rare but serious drug-related adverse event. To identify pharmacogenomic markers of MRONJ associated with bisphosphonate therapy, we conducted a genomewide association study (GWAS) meta-analysis followed by functional analysis of 5,008 individuals of European ancestry treated with bisphosphonates, which includes the largest number of MRONJ cases to date (444 cases and 4,564 controls). Discovery GWAS was performed in randomly selected 70% of the patients with cancer and replication GWAS was performed in the remaining 30% of the patients with cancer treated with intravenous bisphosphonates followed by meta-analysis of all 3,639 patients with cancer. GWAS was also performed in 1,369 patients with osteoporosis treated with oral bisphosphonates. The lead single-nucleotide polymorphism (SNP), rs2736308 on chromosome 8, was associated with an increased risk of MRONJ with an odds ratio (OR) of 2.71 and 95% confidence interval (CI) of 1.90-3.86 (P = 3.57*10-8 ) in the meta-analysis of patients with cancer. This SNP was validated in the MRONJ GWAS in patients with osteoporosis (OR: 2.82, 95% CI: 1.55-4.09, P = 6.84*10-4 ). The meta-analysis combining patients with cancer and patients with osteoporosis yielded the same lead SNP rs2736308 on chromosome 8 as the top SNP (OR: 2.74, 95% CI: 2.09-3.39, P = 9.65*10-11 ). This locus is associated with regulation of the BLK, CTSB, and FDFT1 genes, which had been associated with bone mineral density. FDFT1 encodes a membrane-associated enzyme, which is implicated in the bisphosphonate pathway. This study provides insights into the potential mechanism of MRONJ.

In a large cohort of 373 pediatric patients with Marfan syndrome (MFS) with a severe cardiovascular phenotype, we explored the proportion of patients with MFS with a pathogenic FBN1 variant and analyzed whether the type/location of FBN1 variants was associated with specific clinical characteristics and response to treatment. Patients were recruited on the basis of the following criteria: aortic root z-score > 3, age 6 months to 25 years, no prior or planned surgery, and aortic root diameter < 5 cm.

Propafenone is an antiarrhythmic drug metabolized primarily by cytochrome P450 2D6 (CYP2D6). In adults, propafenone adverse events (AEs) are associated with CYP2D6 poor metabolizer status; however, pediatric data are lacking. Subjects were tested for 10 CYP2D6 allelic variants and copy number status, and activity scores assigned to each genotype. Seventy-six individuals (median 0.3 [range 0-26] years old) were included. Propafenone AEs occurred in 29 (38%); 14 (18%) required drug discontinuation due to AE. The most common AEs were QRS (n = 10) and QTc (n = 6) prolongation. Those with AEs were older at the time of propafenone initiation (1.58 [0.13-9.92] vs. 0.20 [0.08-2.01] years old; p = 0.042). CYP2D6 activity scores were not associated with presence of an AE (odds ratio [OR] 0.48 [0.22-1.03]; p = 0.055) but with the total number of AE (β1  = -0.31 [-0.60, -0.03]; p = 0.029), systemic AEs (OR 0.33 [0.13-0.88]; p = 0.022), and drug discontinuation for systemic AEs (OR 0.28 [0.09-0.83]; p = 0.017). Awareness of CYP2D6 activity score and patient age may aid in determining an individual's risk for an AE with propafenone administration.

This review evaluates the pediatric evidence for pharmacogenetic associations for drugs that are commonly prescribed by or encountered by pediatric clinicians across multiple subspecialties, organized from most to least pediatric evidence. We begin with the pharmacogenetic research that led to the warning of increased risk of death in certain pediatric populations ("ultrarapid metabolizers") who are prescribed codeine after tonsillectomy or adenoidectomy. We review the evidence for genetic testing for thiopurine metabolism, which has become routine in multiple pediatric subspecialties. We discuss the pharmacogenetic research in proton pump inhibitors, for which clinical guidelines have recently been made available. With an increase in the prevalence of behavioral health disorders including attention deficit hyperactivity disorder (ADHD), we review the pharmacogenetic literature on selective serotonin reuptake inhibitors, selective norepinephrine reuptake inhibitors, and ADHD medications. We will conclude this section on the current pharmacogenetic data on ondansetron. We also provide our perspective on how to integrate the current research on pharmacogenetics into clinical care and what further research is needed. We discuss how institutions are managing pharmacogenetic test results and implementing them clinically, and how the electronic health record can be leveraged to ensure testing results are available and taken into consideration when prescribing medications. IMPACT: While many reviews of pharmacogenetics literature are available, there are few focused on pediatrics. Pediatricians across subspecialties will become more comfortable with pharmacogenetics terminology, know resources they can use to help inform their prescribing habits for drugs with known pharmacogenetic associations, and understand the limitations of testing and where further research is needed.

There is growing interest in utilizing pharmacogenetic (PGx) testing to guide antidepressant use, but there is lack of clarity on how to implement testing into clinical practice. We administered two surveys at 17 sites that had implemented or were in the process of implementing PGx testing for antidepressants. Survey 1 collected data on the process and logistics of testing. Survey 2 asked sites to rank the importance of Consolidated Framework for Implementation Research (CFIR) constructs using best-worst scaling choice experiments. Of the 17 sites, 13 had implemented testing and four were in the planning stage. Thirteen offered testing in the outpatient setting, and nine in both outpatient/inpatient settings. PGx tests were mainly ordered by psychiatry (92%) and primary care (69%) providers. CYP2C19 and CYP2D6 were the most commonly tested genes. The justification for antidepressants selected for PGx guidance was based on Clinical Pharmacogenetics Implementation Consortium guidelines (94%) and US Food and Drug Administration (FDA; 75.6%) guidance. Both institutional (53%) and commercial laboratories (53%) were used for testing. Sites varied on the methods for returning results to providers and patients. Sites were consistent in ranking CFIR constructs and identified patient needs/resources, leadership engagement, intervention knowledge/beliefs, evidence strength and quality, and the identification of champions as most important for implementation. Sites deployed similar implementation strategies and measured similar outcomes. The process of implementing PGx testing to guide antidepressant therapy varied across sites, but key drivers for successful implementation were similar and may help guide other institutions interested in providing PGx-guided pharmacotherapy for antidepressant management.

The goal of this study was to determine whether CYP2D6 metabolizer status within the ondansetron-treated pediatric tonsillectomy population is associated with risk of postoperative nausea and vomiting (PONV) in the post-anesthesia care unit. We conducted a retrospective cohort study of pediatric patients (<18 years) who underwent tonsillectomy and received ondansetron on the day of the procedure. Data were obtained from BioVU, an institutional biobank that links DNA to de-identified electronic health record data. Subjects were tested for 10 CYP2D6 allelic variants and copy number variation, and genotype data translated into CYP2D6 metabolizer status. The cohort included 652 individuals, 105 (16.1%) of whom had PONV. Rates of PONV were similar across groups: ultrarapid metabolizers (UMs), 1 of 9 (11.1%); normal metabolizers (NMs), 64 of 354 (18.1%); intermediate metabolizers (IMs), 33 of 234 (14.1%); poor metabolizers (PMs), 6 of 39 (15.4%); and ambiguous phenotypes, 1 of 16 (6.3%). In multivariable analysis adjusted for age, sex, and time under anesthesia, CYP2D6 metabolizer status was not associated with PONV, with an odds ratio of 1.37 (95% confidence interval 0.9, 2.1) when comparing PM/IM versus NM/UM. In this large pediatric population, no significant differences were detected for PONV based on CYP2D6 metabolizer status. Further investigation is needed to determine mechanisms for ondansetron inefficacy in children.

There are challenges in achieving and maintaining therapeutic tacrolimus levels after solid organ transplantation (SOT). The purpose of this genome-wide association study was to generate an integrated clinical and genetic prediction model for tacrolimus levels in pediatric SOT.

Defining the full spectrum of human disease associated with a biomarker is necessary to advance the biomarker into clinical practice. We hypothesize that associating biomarker measurements with electronic health record (EHR) populations based on shared genetic architectures would establish the clinical epidemiology of the biomarker. We use Bayesian sparse linear mixed modeling to calculate SNP weightings for 53 biomarkers from the Atherosclerosis Risk in Communities study. We use the SNP weightings to computed predicted biomarker values in an EHR population and test associations with 1139 diagnoses. Here we report 116 associations meeting a Bonferroni level of significance. A false discovery rate (FDR)-based significance threshold reveals more known and undescribed associations across a broad range of biomarkers, including biometric measures, plasma proteins and metabolites, functional assays, and behaviors. We confirm an inverse association between LDL-cholesterol level and septicemia risk in an independent epidemiological cohort. This approach efficiently discovers biomarker-disease associations.

High-dose (HD) methotrexate (MTX) is a critical component of treatment for hematologic malignancies in children and young adults. Therapeutic drug monitoring is necessary due to substantial interindividual variation in MTX clearance. Common function-altering polymorphisms in SLCO1B1 (encodes OATP1B1, which transports MTX) may contribute to clearance variability. We performed pharmacokinetic modeling using data for 106 children and young adults treated with HD MTX for hematologic malignancies; of 396 total courses of HD MTX, 360 consisted of 5 g/m2 over 24 hours. We evaluated the contribution of clinical covariates and SLCO1B1 genotype (388A>G and 521T>C) to MTX clearance variability. Of the clinical covariates studied, patient weight improved the pharmacokinetic model most significantly (P < 0.001). The addition of the SLCO1B1 variants individually further improved the model (P < 0.05 for each). An interaction between these variants was suggested when both were included (P = 0.017). SLCO1B1 genotype should be considered in efforts to personalize HD MTX dosing.

Pediatric acute kidney injury (AKI) is common and associated with increased morbidity, mortality, and length of stay. We performed a pragmatic randomized trial testing the hypothesis that AKI risk alerts increase AKI screening.

Finding the optimal treatment for a chronic condition can be a complex and lengthy endeavor for both the patient and the clinician. To address this challenge, we developed an "N-of-1" quality improvement infrastructure to aid providers and patients in personalized treatment decision-making using systematic assessment of patient-reported outcomes during routine care.

Alpha-1 antitrypsin deficiency (AATD), a relatively common autosomal recessive genetic disorder, is underdiagnosed in symptomatic individuals. We sought to compare the risk of liver transplantation associated with hepatitis C infection with AATD heterozygotes and homozygotes and determine if SERPINA1 sequencing would identify undiagnosed AATD. We performed a retrospective cohort study in a deidentified Electronic Health Record (EHR)-linked DNA biobank with 72,027 individuals genotyped for the M, Z, and S alleles in SERPINA1. We investigated liver transplantation frequency by genotype group and compared with hepatitis C infection. We performed SERPINA1 sequencing in carriers of pathogenic AATD alleles who underwent liver transplantation. Liver transplantation was associated with the Z allele (ZZ: odds ratio [OR] = 1.31, p<2e-16; MZ: OR = 1.02, p = 1.2e-13) and with hepatitis C (OR = 1.20, p<2e-16). For liver transplantation, there was a significant interaction between genotype and hepatitis C (ZZ: interaction OR = 1.23, p = 4.7e-4; MZ: interaction OR = 1.11, p = 6.9e-13). Sequencing uncovered a second, rare, pathogenic SERPINA1 variant in six of 133 individuals with liver transplants and without hepatitis C. Liver transplantation was more common in individuals with AATD risk alleles (including heterozygotes), and AATD and hepatitis C demonstrated evidence of a gene-environment interaction in relation to liver transplantation. The current AATD screening strategy may miss diagnoses whereas SERPINA1 sequencing may increase diagnostic yield for AATD, stratify risk for liver disease, and inform clinical management for individuals with AATD risk alleles and liver disease risk factors.

BackgroundAcute kidney injury (AKI) is common in pediatric inpatients and is associated with increased morbidity, mortality, and length of stay. Its early identification can reduce severity.MethodsTo create and validate an electronic health record (EHR)-based AKI screening tool, we generated temporally distinct development and validation cohorts using retrospective data from our tertiary care children's hospital, including children aged 28 days through 21 years with sufficient serum creatinine measurements to determine AKI status. AKI was defined as 1.5-fold or 0.3 mg/dl increase in serum creatinine. Age, medication exposures, platelet count, red blood cell distribution width, serum phosphorus, serum transaminases, hypotension (ICU only), and pH (ICU only) were included in AKI risk prediction models.ResultsFor ICU patients, 791/1,332 (59%) of the development cohort and 470/866 (54%) of the validation cohort had AKI. In external validation, the ICU prediction model had a c-statistic=0.74 (95% confidence interval 0.71-0.77). For non-ICU patients, 722/2,337 (31%) of the development cohort and 469/1,474 (32%) of the validation cohort had AKI, and the prediction model had a c-statistic=0.69 (95% confidence interval 0.66-0.72).ConclusionsAKI screening can be performed using EHR data. The AKI screening tool can be incorporated into EHR systems to identify high-risk patients without serum creatinine data, enabling targeted laboratory testing, early AKI identification, and modification of care.

The increased availability of clinical pharmacogenetic (PGx) guidelines and decreasing costs for genetic testing have slowly led to increased utilization of PGx testing in clinical practice. Pre-emptive PGx testing, where testing is performed in advance of drug prescribing, is one means to ensure results are available at the time of prescribing decisions. However, the most efficient and effective methods to clinically implement this strategy remain unclear.

Vancomycin, a commonly used antibiotic, can be nephrotoxic. Known risk factors such as age, creatinine clearance, vancomycin dose / dosing interval, and concurrent nephrotoxic medications fail to accurately predict nephrotoxicity. To identify potential genomic risk factors, we performed a genome-wide association study (GWAS) of serum creatinine levels while on vancomycin in 489 European American individuals and validated findings in three independent cohorts totaling 439 European American individuals. In primary analyses, the chromosome 6q22.31 locus was associated with increased serum creatinine levels while on vancomycin therapy (most significant variant rs2789047, risk allele A, β = -0.06, p = 1.1 x 10(-7)). SNPs in this region had consistent directions of effect in the validation cohorts, with a meta-p of 1.1 x 10(-7). Variation in this region on chromosome 6, which includes the genes TBC1D32/C6orf170 and GJA1 (encoding connexin43), may modulate risk of vancomycin-induced kidney injury.

A single mutation can alter cellular and global homeostatic mechanisms and give rise to multiple clinical diseases. We hypothesized that these disease mechanisms could be identified using low minor allele frequency (MAF<0.1) non-synonymous SNPs (nsSNPs) associated with "mechanistic phenotypes", comprised of collections of related diagnoses. We studied two mechanistic phenotypes: (1) thrombosis, evaluated in a population of 1,655 African Americans; and (2) four groupings of cancer diagnoses, evaluated in 3,009 white European Americans. We tested associations between nsSNPs represented on GWAS platforms and mechanistic phenotypes ascertained from electronic medical records (EMRs), and sought enrichment in functional ontologies across the top-ranked associations. We used a two-step analytic approach whereby nsSNPs were first sorted by the strength of their association with a phenotype. We tested associations using two reverse genetic models and standard additive and recessive models. In the second step, we employed a hypothesis-free ontological enrichment analysis using the sorted nsSNPs to identify functional mechanisms underlying the diagnoses comprising the mechanistic phenotypes. The thrombosis phenotype was solely associated with ontologies related to blood coagulation (Fisher's p = 0.0001, FDR p = 0.03), driven by the F5, P2RY12 and F2RL2 genes. For the cancer phenotypes, the reverse genetics models were enriched in DNA repair functions (p = 2×10-5, FDR p = 0.03) (POLG/FANCI, SLX4/FANCP, XRCC1, BRCA1, FANCA, CHD1L) while the additive model showed enrichment related to chromatid segregation (p = 4×10-6, FDR p = 0.005) (KIF25, PINX1). We were able to replicate nsSNP associations for POLG/FANCI, BRCA1, FANCA and CHD1L in independent data sets. Mechanism-oriented phenotyping using collections of EMR-derived diagnoses can elucidate fundamental disease mechanisms.