Clinical associations with a polygenic predisposition to benign lower white blood cell counts.

Polygenic variation unrelated to disease contributes to interindividual variation in baseline white blood cell (WBC) counts, but its clinical significance is uncharacterized. We investigated the clinical consequences of a genetic predisposition toward lower WBC counts among 89,559 biobank participants from tertiary care centers using a polygenic score for WBC count (PGSWBC) comprising single nucleotide polymorphisms not associated with disease. A predisposition to lower WBC counts was associated with a decreased risk of identifying pathology on a bone marrow biopsy performed for a low WBC count (odds-ratio = 0.55 per standard deviation increase in PGSWBC [95%CI, 0.30-0.94], p = 0.04), an increased risk of leukopenia (a low WBC count) when treated with a chemotherapeutic (n = 1724, hazard ratio [HR] = 0.78 [0.69-0.88], p = 4.0 × 10-5) or immunosuppressant (n = 354, HR = 0.61 [0.38-0.99], p = 0.04). A predisposition to benign lower WBC counts was associated with an increased risk of discontinuing azathioprine treatment (n = 1,466, HR = 0.62 [0.44-0.87], p = 0.006). Collectively, these findings suggest that there are genetically predisposed individuals who are susceptible to escalations or alterations in clinical care that may be harmful or of little benefit.

Composite CYP3A (CYP3A4 and CYP3A5) phenotypes and influence on tacrolimus dose adjusted concentrations in adult heart transplant recipients.

CYP3A5 genetic variants are associated with tacrolimus metabolism. Controversy remains on whether CYP3A4 increased [*1B (rs2740574), *1 G (rs2242480)] and decreased function [*22 (rs35599367)] genetic variants provide additional information. This retrospective cohort study aims to address whether tacrolimus dose-adjusted trough concentrations differ between combined CYP3A (CYP3A5 and CYP3A4) phenotype groups. Heart transplanted patients (n = 177, between 2008 and 2020) were included and median age was 54 years old. Significant differences between CYP3A phenotype groups in tacrolimus dose-adjusted trough concentrations were found in the early postoperative period and continued to 6 months post-transplant. In CYP3A5 nonexpressers, carriers of CYP3A4*1B or *1 G variants (Group 3) compared to CYP3A4*1/*1 (Group 2) patients were found to have lower tacrolimus dose-adjusted trough concentrations at 2 months. In addition, significant differences were found among CYP3A phenotype groups in the dose at discharge and time to therapeutic range while time in therapeutic range was not significantly different. A combined CYP3A phenotype interpretation may provide more nuanced genotype-guided TAC dosing in heart transplant recipients.

Association Between Neighborhood Socioeconomic Factors and Length of Stay After Surgery for Congenital Heart Disease.

OBJECTIVES: Neighborhood socioeconomic status, as measured by area deprivation index (ADI) is associated with longer length of stay (LOS) after surgery for hypoplastic left heart syndrome. We tested the hypothesis that LOS is associated with ADI in a large cohort of congenital heart disease (CHD) surgical cases of varying severity and sought to determine which other components of the ADI accounted for any associations. DESIGN: Retrospective analysis of a curated dataset. The Brokamp ADI was determined using residential addresses. Overall, ADI and each of its six individual components were dichotomized, and LOS compared between groups above versus below the median for the entire cohort and after stratifying by surgical The Society of Thoracic Surgeons-European Association for Cardio-Thoracic Surgery (STAT) severity category. SETTING: Single-center academic pediatric teaching hospital. PATIENTS: CHD patients who underwent surgical repair/palliation between September 2007 and August 2022. INTERVENTIONS: None. MEASUREMENTS AND MAIN RESULTS: A total of 2462 patients (52.7% male) were included. Median age was 254 (interquartile range [IQR] 95-1628) days and median LOS in the hospital was 8 (IQR 5-18) days. We failed to identify an association between Brokamp ADI, above versus below the median for the entire cohort, and LOS; nor in STAT categories 1-4. However, in STAT category 5 (n = 129) those with ADI above the median (more deprived) had a significantly longer LOS (48 [20-88] vs. 36 [18-49] d, p = 0.034). Of the individual components of the ADI, only percent below poverty level and percent vacant houses were associated with LOS in STAT category 5. CONCLUSIONS: LOS after CHD surgery is associated with Brokamp ADI in STAT category 5 cases, we failed to identify an association in lower-risk cardiac operations.

Developing and evaluating pediatric phecodes (Peds-Phecodes) for high-throughput phenotyping using electronic health records.

OBJECTIVE: Pediatric patients have different diseases and outcomes than adults; however, existing phecodes do not capture the distinctive pediatric spectrum of disease. We aim to develop specialized pediatric phecodes (Peds-Phecodes) to enable efficient, large-scale phenotypic analyses of pediatric patients. MATERIALS AND METHODS: We adopted a hybrid data- and knowledge-driven approach leveraging electronic health records (EHRs) and genetic data from Vanderbilt University Medical Center to modify the most recent version of phecodes to better capture pediatric phenotypes. First, we compared the prevalence of patient diagnoses in pediatric and adult populations to identify disease phenotypes differentially affecting children and adults. We then used clinical domain knowledge to remove phecodes representing phenotypes unlikely to affect pediatric patients and create new phecodes for phenotypes relevant to the pediatric population. We further compared phenome-wide association study (PheWAS) outcomes replicating known pediatric genotype-phenotype associations between Peds-Phecodes and phecodes. RESULTS: The Peds-Phecodes aggregate 15 533 ICD-9-CM codes and 82 949 ICD-10-CM codes into 2051 distinct phecodes. Peds-Phecodes replicated more known pediatric genotype-phenotype associations than phecodes (248 vs 192 out of 687 SNPs, P < .001). DISCUSSION: We introduce Peds-Phecodes, a high-throughput EHR phenotyping tool tailored for use in pediatric populations. We successfully validated the Peds-Phecodes using genetic replication studies. Our findings also reveal the potential use of Peds-Phecodes in detecting novel genotype-phenotype associations for pediatric conditions. We expect that Peds-Phecodes will facilitate large-scale phenomic and genomic analyses in pediatric populations. CONCLUSION: Peds-Phecodes capture higher-quality pediatric phenotypes and deliver superior PheWAS outcomes compared to phecodes.

Developing and Evaluating Pediatric Phecodes (Peds-Phecodes) for High-Throughput Phenotyping Using Electronic Health Records.

Objective: Pediatric patients have different diseases and outcomes than adults; however, existing phecodes do not capture the distinctive pediatric spectrum of disease. We aim to develop specialized pediatric phecodes (Peds-Phecodes) to enable efficient, large-scale phenotypic analyses of pediatric patients. Materials and Methods: We adopted a hybrid data- and knowledge-driven approach leveraging electronic health records (EHRs) and genetic data from Vanderbilt University Medical Center to modify the most recent version of phecodes to better capture pediatric phenotypes. First, we compared the prevalence of patient diagnoses in pediatric and adult populations to identify disease phenotypes differentially affecting children and adults. We then used clinical domain knowledge to remove phecodes representing phenotypes unlikely to affect pediatric patients and create new phecodes for phenotypes relevant to the pediatric population. We further compared phenome-wide association study (PheWAS) outcomes replicating known pediatric genotype-phenotype associations between Peds-Phecodes and phecodes. Results: The Peds-Phecodes aggregate 15,533 ICD-9-CM codes and 82,949 ICD-10-CM codes into 2,051 distinct phecodes. Peds-Phecodes replicated more known pediatric genotype-phenotype associations than phecodes (248 versus 192 out of 687 SNPs, p<0.001). Discussion: We introduce Peds-Phecodes, a high-throughput EHR phenotyping tool tailored for use in pediatric populations. We successfully validated the Peds-Phecodes using genetic replication studies. Our findings also reveal the potential use of Peds-Phecodes in detecting novel genotype-phenotype associations for pediatric conditions. We expect that Peds-Phecodes will facilitate large-scale phenomic and genomic analyses in pediatric populations. Conclusion: Peds-Phecodes capture higher-quality pediatric phenotypes and deliver superior PheWAS outcomes compared to phecodes.

Clinical consequences of a polygenic predisposition to benign lower white blood cell counts: Consequences of benign WBC count genetics.

Polygenic variation unrelated to disease contributes to interindividual variation in baseline white blood cell (WBC) counts, but its clinical significance is undefined. We investigated the clinical consequences of a genetic predisposition toward lower WBC counts among 89,559 biobank participants from tertiary care centers using a polygenic score for WBC count (PGSWBC) comprising single nucleotide polymorphisms not associated with disease. A predisposition to lower WBC counts was associated with a decreased risk of identifying pathology on a bone marrow biopsy performed for a low WBC count (odds-ratio=0.55 per standard deviation increase in PGSWBC [95%CI, 0.30 - 0.94], p=0.04), an increased risk of leukopenia (a low WBC count) when treated with a chemotherapeutic (n=1,724, hazard ratio [HR]=0.78 [0.69 - 0.88], p=4.0×10-5) or immunosuppressant (n=354, HR=0.61 [0.38 - 0.99], p=0.04). A predisposition to benign lower WBC counts was associated with an increased risk of discontinuing azathioprine treatment (n=1,466, HR=0.62 [0.44 - 0.87], p=0.006). Collectively, these findings suggest that a WBC count polygenic score identifies individuals who are susceptible to escalations or alterations in clinical care that may be harmful or of little benefit.

Population genetic testing and SERPINA1 sequencing identifies unidentified alpha-1 antitrypsin deficiency alleles and gene-environment interaction with hepatitis C infection.

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.

Pediatric Persistent Inflammation, Immunosuppression, and Catabolism Syndrome Prevalence in Sepsis-Related Mortalities: A 23-Year Institutional History.

BACKGROUND: Delayed mortality in sepsis often is linked to a lack of resolution in the inflammatory cascade termed persistent inflammation, immunosuppression, and catabolism syndrome (PICS). Limited research exists on PICS in pediatric patients with sepsis. RESEARCH QUESTION: What is the prevalence of pediatric PICS (pPICS) in patients who died of sepsis-related causes and what associated pathogen profiles and comorbidities did they have compared with those patients without pPICS who died from sepsis? STUDY DESIGN AND METHODS: A retrospective study of a single institution using a de-identified database from 1997 through 2020 for all patients aged 21 years or younger who died of culture-positive sepsis from a known source and who had laboratory data available were evaluated for the presence of pPICS. RESULTS: Among records extracted from the institutional database, 557 patients had culture-positive sepsis, with 262 patients having pPICS (47%). Patients with pPICS were more likely to have underlying hematologic or oncologic disease or cardiac disease. In addition, patients who had pPICS showed increased odds of associated fungal infection compared with those patients who did not (OR, 2.69; 95% CI, 1.59-4.61; P < .001). When assessing laboratory criteria, having a sustained absolute lymphocyte count of < 1.0 × 103/µL was most closely associated with having pPICS compared with other laboratory parameters. Finally, the results of multivariate logistic regression analysis indicated that patients with pPICS were more common in the cardiac ICU, as opposed to the PICU (OR, 3.43; CI, 1.57-7.64; P = .002). INTERPRETATION: Pediatric patients who died of a sepsis-related cause have a pPICS phenotype nearly one-half of the time. These patients are more likely to be in the cardiac ICU than the pediatric ICU and have associated fungal infections. Special attention should be directed toward this population in future research.

Pediatric considerations for pharmacogenetic selective serotonin reuptake inhibitors clinical decision support.

Pharmacogenetic testing for psychiatry is growing at a rapid pace, with multiple sites utilizing results to help clinical decision-making. Genotype-guided dosing and drug selection have been implemented at several sites, including Vanderbilt University Medical Center, where clinical decision support (CDS) based on pharmacogenetic results went live for selective serotonin reuptake inhibitors in 2020 for both adult and pediatric patients. Effective and appropriate implementation of CYP2D6- and CYP2C19-guided CDS for the pediatric population requires consideration of the evidence for the pharmacogenetic associations, medication indications, and appropriate alternative therapies to be used when a pharmacogenetic contraindication is identified. In this article, we review these pediatric pharmacogenetic considerations for selective serotonin reuptake inhibitor CDS. We include a case study, the current literature supporting clinical recommendations, considerations when designing pediatric CDS, future implications, and examples of sertraline, (es)citalopram, paroxetine, and fluvoxamine alerts.

Tutorial: Using Community Engagement Studios to Enhance Pharmacogenetic Study Design for Maximizing Enrollment of Diverse Children and Pregnant People.

Most pharmacogenetic research is conducted in adult, non-pregnant populations of European ancestry. Study of more diverse and special populations is necessary to validate findings and improve health equity. However, there are significant barriers to recruitment of diverse populations for genetic studies, such as mistrust of researchers due to a history of unethical research and ongoing social inequities. Engaging communities and understanding community members' perspectives may help to overcome these barriers and improve research quality. Here, we highlight one method for engaging communities, the Community Engagement Studio (CES), a consultative session that allows researchers to obtain guidance and feedback based on community members' lived experiences. We also provide an example of its use in pharmacogenetic studies. In designing a survey study of knowledge and attitudes around pharmacogenetic testing among children with chronic conditions and pregnant individuals, we sought input from diverse community stakeholders through CESs at Vanderbilt University Medical Center. We participated in two CESs with community stakeholders representing study target populations. Our goals were to learn specific concerns about pharmacogenetic testing and preferred recruitment strategies for these communities. Concerns were expressed about how genetic information would be used beyond the immediate study. Participants emphasized the importance of clarity and transparency in communication to overcome participation hesitancy and mistrust of the study team. Recruitment strategy recommendations ranged from informal notices posted in healthcare settings to provider referrals. The CES enabled us to modify our recruitment methods and research materials to better communicate with populations currently under-represented in pharmacogenetics research.

Pediatric CYP2D6 metabolizer status and post-tonsillectomy nausea and vomiting after ondansetron.

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.

Bronchopulmonary dysplasia is associated with polyhydramnios in a scan for novel perinatal risk factors.

BACKGROUND: The pathogenesis of bronchopulmonary dysplasia (BPD) is multifactorial, and there are limited data about prenatal exposures and risk of BPD. STUDY DESIGN: Our study performed parallel analyses using a logistic regression model in a cohort of 4527 infants with data from a curated registry and using a phenome wide association study (PheWAS) based on ICD9/10-based phecodes. We examined 20 prenatal exposures from a neonatal intensive care unit (NICU) curated registry database related to pregnancy and maternal health as well as 94 maternal diagnosis phecodes with a PheWAS analysis. RESULT: In both the curated registry and PheWAS analyses, polyhydramnios was associated with an increased risk of BPD (OR 5.70, 95% CI 2.78-11.44, p = 1.37 × 10-6). CONCLUSION: Our data suggest that polyhydramnios may be a clinical indicator of premature infants at increased risk for bronchopulmonary dysplasia. Combining curated registry data with PheWAS analysis creates a valuable tool to generate hypotheses. IMPACT: Polyhydramnios was significantly associated with bronchopulmonary dysplasia in both a curated registry and by ICD coding analysis with a phenome wide association study (PheWAS). Preterm polyhydramnios may be a clinical indicator of infants at increased risk for developing bronchopulmonary dysplasia after preterm birth. Combining curated registry with PheWAS analysis creates a valuable tool to generate hypotheses about perinatal risk factors and morbidities associated with preterm birth.

Medication history-wide association studies for pharmacovigilance of pregnant patients.

Background: Systematic exclusion of pregnant people from interventional clinical trials has created a public health emergency for millions of patients through a dearth of robust safety data for common drugs. Methods: We harnessed an enterprise collection of 2.8 M electronic health records (EHRs) from routine care, leveraging data linkages between mothers and their babies to detect drug safety signals in this population at full scale. Our mixed-methods signal detection approach stimulates new hypotheses for post-marketing surveillance agnostically of both drugs and diseases-by identifying 1,054 drugs historically prescribed to pregnant patients; developing a quantitative, medication history-wide association study; and integrating a qualitative evidence synthesis platform using expert clinician review for integration of biomedical specificity-to test the effects of maternal exposure to diverse drugs on the incidence of neurodevelopmental defects in their children. Results: We replicated known teratogenic risks and existing knowledge on drug structure-related teratogenicity; we also highlight 5 common drug classes for which we believe this work warrants updated assessment of their safety. Conclusion: Here, we present roots of an agile framework to guide enhanced medication regulations, as well as the ontological and analytical limitations that currently restrict the integration of real-world data into drug safety management during pregnancy. This research is not a replacement for inclusion of pregnant people in prospective clinical studies, but it presents a tractable team science approach to evaluating the utility of EHRs for new regulatory review programs-towards improving the delicate equipoise of accuracy and ethics in assessing drug safety in pregnancy.

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.

Implementing a pragmatic clinical trial to tailor opioids for acute pain on behalf of the IGNITE ADOPT PGx investigators.

Opioid prescribing for postoperative pain management is challenging because of inter-patient variability in opioid response and concern about opioid addiction. Tramadol, hydrocodone, and codeine depend on the cytochrome P450 2D6 (CYP2D6) enzyme for formation of highly potent metabolites. Individuals with reduced or absent CYP2D6 activity (i.e., intermediate metabolizers [IMs] or poor metabolizers [PMs], respectively) have lower concentrations of potent opioid metabolites and potentially inadequate pain control. The primary objective of this prospective, multicenter, randomized pragmatic trial is to determine the effect of postoperative CYP2D6-guided opioid prescribing on pain control and opioid usage. Up to 2020 participants, age ≥8 years, scheduled to undergo a surgical procedure will be enrolled and randomized to immediate pharmacogenetic testing with clinical decision support (CDS) for CYP2D6 phenotype-guided postoperative pain management (intervention arm) or delayed testing without CDS (control arm). CDS is provided through medical record alerts and/or a pharmacist consult note. For IMs and PM in the intervention arm, CDS includes recommendations to avoid hydrocodone, tramadol, and codeine. Patient-reported pain-related outcomes are collected 10 days and 1, 3, and 6 months after surgery. The primary outcome, a composite of pain intensity and opioid usage at 10 days postsurgery, will be compared in the subgroup of IMs and PMs in the intervention (n = 152) versus the control (n = 152) arm. Secondary end points include prescription pain medication misuse scores and opioid persistence at 6 months. This trial will provide data on the clinical utility of CYP2D6 phenotype-guided opioid selection for improving postoperative pain control and reducing opioid-related risks.

Frequency of benign neutropenia among Black versus White individuals undergoing a bone marrow assessment.

Healthy individuals in the United States identified as having Black race have lower neutrophil counts, on average, than individuals identified as having White race, which could result in more negative diagnostic evaluations for neutropenia. To test this hypothesis, the proportion of evaluations where the final diagnosis was clinically insignificant neutropenia for Black and White individuals who underwent an evaluation by a haematologist that included a bone marrow (BM) biopsy to investigate neutropenia was assessed. 172 individuals without prior haematological diagnoses who underwent a haematological evaluation to investigate neutropenia. Individuals diagnosed with clinically insignificant neutropenia between Black and White individuals were compared using a propensity-score-adjusted logistic regression. Of 172 individuals, 42 (24%) were classified as Black race, 86 (50%) were males, and the 79 (46%) were over 18 years old. A BM biopsy did not identify pathology in 95% (40 of 42) of Black individuals and 68% (89 of 130) of White Individuals. Black individuals (25 of 42 [60%]) received a final diagnosis of clinically insignificant neutropenia, compared to White individuals (12 of 130 [9%]) (adjusted odds ratio =7.9, 95% CI: 3.1 - 21.1). We conclude that black individuals were more likely to receive a diagnosis of clinically insignificant neutropenia after haematological assessment.

Race, Genotype, and Azathioprine Discontinuation : A Cohort Study.

BACKGROUND: Thiopurines are an important class of immunosuppressants despite their risk for hematopoietic toxicity and narrow therapeutic indices. Benign neutropenia related to an ACKR1 variant (rs2814778-CC) is common among persons of African ancestries. OBJECTIVE: To test whether rs2814778-CC was associated with azathioprine discontinuation attributed to hematopoietic toxicity and lower thiopurine dosing. DESIGN: Retrospective cohort study. SETTING: Two tertiary care centers. PATIENTS: Thiopurine users with White or Black race. MEASUREMENTS: Azathioprine discontinuation attributed to hematopoietic toxicity. Secondary outcomes included weight-adjusted final dose, leukocyte count, and change in leukocyte count. RESULTS: The rate of azathioprine discontinuation attributed to hematopoietic toxicity was 3.92 per 100 person-years among patients with the CC genotype (n = 101) and 1.34 per 100 person-years among those with the TT or TC genotype (n = 1365) (hazard ratio [HR] from competing-risk model, 2.92 [95% CI, 1.57 to 5.41]). The risk remained significant after adjustment for race (HR, 2.61 [CI, 1.01 to 6.71]). The risk associated with race alone (HR, 2.13 [CI, 1.21 to 3.75]) was abrogated by adjustment for genotype (HR, 1.13 [CI, 0.48 to 2.69]). Lower last leukocyte count and lower dosing were significant among patients with the CC genotype. Lower dosing was validated in an external cohort of 94 children of African ancestries prescribed the thiopurine 6-mercaptopurine (6-MP) for acute lymphoblastic leukemia. The CC genotype was independently associated with lower 6-MP dose intensity relative to the target daily dose of 75 mg/m2 (median, 0.83 [IQR, 0.70 to 0.94] for the CC genotype vs. 0.94 [IQR, 0.72 to 1.13] for the TT or TC genotype; P = 0.013). LIMITATIONS: Unmeasured confounding; data limited to tertiary centers. CONCLUSION: Patients with the CC genotype had higher risk for azathioprine discontinuation attributed to hematopoietic toxicity and lower thiopurine doses. Genotype was associated with those risks, even after adjustment for race. PRIMARY FUNDING SOURCE: National Institutes of Health.