|Home | About | Journals | Submit | Contact Us | Français|
The aim of this study was to examine the association between the CYP2D6 and CYP2C19 genotype-predicted combined phenotypes and short-term measures of psychotropic efficacy and toxicity.
A rater-blinded, retrospective genotype association design examined a cohort of hospitalized pediatric psychiatric patients genotyped for CYP2D6 and CYP2C19 as part of clinical care. These combined genotypes were used to predict a combined phenotype. The primary efficacy outcome measure was the behavior intervention score (BIS), a function of the number of recorded timeouts/seclusions, therapeutic holds, and physical restraints. Drug tolerability was defined as the total number of recorded adverse drug reactions.
Primary analysis was performed on 279 pediatric patients taking CYP2D6- or CYP2C19- dependent psychotropics. Combined phenotype was associated with BIS (p=0.01) and number of adverse drug reactions (p=0.03). Combined poor metabolizers treated with psychotropics had the lowest BIS (highest efficacy) and the highest number of adverse drug reactions. Combined ultrarapid metabolizers had the highest BIS (lowest efficacy) and the lowest number of adverse drug reactions.
Common variants in CYP2D6 and CYP2C19 are associated with the short-term efficacy and tolerability of psychotropic medications in hospitalized pediatric patients.
Mental health problems among American children have been referred to as a national public health crisis by the President's New Freedom Commission on Mental Health, subcommittee on Children and Family (New Freedom Commission on Mental Health 2003). Mental illness, significant enough to impair normal development and functioning, affects approximately 10% of children ages 9–17 years (see Goldman et al. 1999). Many children with mental illness require more than outpatient treatment; as a result, psychiatric disorders are the leading cause for hospitalizations for children between 5 and 19 years old (Geller and Biebel 2006). Aggression, a common symptom of many psychiatric diagnoses and often the reason for hospitalization, poses a significant risk and cost to patient and caretaker (Jensen et al. 2007).
Patient response to commonly used psychotropic medications demonstrates significant variability; only 30–75% of patients experience efficacy, whereas 65–75% encounter adverse events (Kirchheiner et al. 2004; Emslie et al. 2006; Kratochvil et al. 2006; Hetrick et al. 2007). Many psychotropic medications are metabolized by cytochrome P450 (CYP) enzymes coded for by the polymorphic genes CYP2D6 (+124030) and CYP2C19 (*124020). The relationships between the distinct CYP2D6 or CYP2C19 genotypes and the pharmacokinetics of psychotropic medications are the basis for genotype-based dosing recommendations for some medications (Kirchheiner et al. 2001).
Despite the prevalence of the problem, the contribution of drug-metabolizing enzyme genotypes to the variability in aggressive behavior treatment response in hospitalized pediatric psychiatric patients is unknown. The purpose of this study was to examine the association between CYP2D6 and CYP2C19 genotype-predicted metabolizing phenotypes and short-term measures of drug efficacy (incidence of behavioral interventions for aggressive behavior) and toxicity (adverse drug reactions) in an inpatient pediatric psychiatric population. We hypothesized common variants in these two genes would be important factors in psychotropic medication clinical response variability in these patients.
The study used a rater-blinded, retrospective genotype association design that involved a cohort of hospitalized psychiatric patients genotyped for CYP2D6 and CYP2C19 as part of clinical care. The study was approved by the first author's Institutional Review Board.
Figure Figure11 details the sampling procedure. Patients were eligible for inclusion in the study if they had genotyping for both CYP2D6 and CYP2C19 between January 1, 2005, and September 30, 2005; had their first inpatient admission to the psychiatric service at the authors' institution between January 1, 2000, and September 30, 2005; and had a combination of CYP2D6/CYP2C19-predicted metabolizing phenotypes that matched one of four (out of a possible six) predetermined combinations (Table 1, stage 1). Patients were excluded if: (1) Complete records from their first inpatient psychiatric hospitalization were not available; (2) the patient was enrolled in a blinded placebo-controlled clinical drug trial; or (3) the patient had a primary medical condition that could contribute to psychiatric symptoms (e.g., substance abuse, brain injury, congenital brain abnormality).
CYP2D6 and CYP2C19 genotyping and phenotype prediction were completed prior to this study as part of clinical care. Genotyping was performed in a Clinical Laboratory Improvement Amendments (CLIA)-approved, College of American Pathologists (CAP)-certified laboratory using DNA extracted from peripheral blood or buccal swab samples. The TaqMan® allelic discrimination system (Applied Biosystems, Forest City, CA) and long polymerase chain reaction (PCR) (Lovlie et al. 1996) were used to analyze DNA for the alleles *1, *3,*4, *5 (deletion) and duplication in CYP2D6 and the *1, *2 alleles in CYP2C19. Genotype results were used to classify a patient's predicted metabolizing phenotype (Table 1, stage 1).
Data were collected from the medical records of the patient's first inpatient psychiatric hospitalization. All chart reviewers were trained by the principal investigator (PI) and were blinded to the patient's CYP2D6 and CYP2C19 genotype and corresponding predicted metabolizing phenotype. The PI and study psychiatric pharmacist were blinded when they re-reviewed >50% of the charts to assure integrity and consistency of extracted data.
Each patient's discharge diagnoses, using the Diagnostic and Statistical Manual of Mental Disorders, 4th edition, Text Revision (DSM-IV-TR) (American Psychiatric Association 2000) multiaxial classification system were taken from the medical record. After all patient data were entered but before the blind was broken, the study's child and adolescent psychiatrist categorized all abstracted axis I diagnoses into one of nine major diagnostic categories (Table 2).
A psychotropic medication was classified as a CYP2D6 and/or CYP2C19 substrate if two or more of the following references listed it as a substrate: Drug-Interaction Table (Flockhart 2007), Lexi-Comp Drug Information Handbook (Lacy et al. 2005), and Micromedex (Micromedex Healthcare Series). The number of psychotropics was obtained by totaling all scheduled psychotropics a patient was taking during the admission.
An algorithm to classify patients by their genotype-predicted metabolizing phenotype groups (called “combined phenotype” groups) was developed to identify patients considered at highest risk for poor efficacy or poor tolerability to routinely prescribed psychotropics (Table 1, stage 2, steps 1–4).
The study's primary efficacy outcome was a measure of aggression severity. Because the clinicians had not administered any objective aggression rating scale prospectively, a transformation was used based on the sum of the recorded behavioral interventions (BI) for disruptive/aggressive behavior and included the number of timeouts or seclusions (BI1), therapeutic holds (BI2), and physical restraints (BI3). A weight-free approach combined the individual recorded behavioral interventions and was called the Behavioral Intervention Score (BIS),
This score is considered weight-free because it is unaffected by different weights applied to the individual BI items (Elston 1963).
Three secondary efficacy outcome measures included the patient's total number of as-needed (PRN) psychotropic medication doses, length of stay (LOS), and change in global assessment of functioning (GAF) score (American Psychiatric Association 2000) from admission to discharge.
The study's primary tolerability outcome measure was the number of recorded adverse drug reactions (ADRs). Symptoms were counted as ADRs if documented during the time frame when scheduled psychotropics were administered and there was: (1) Documentation by a physician/nurse that a symptom was an ADR; or (2) documentation that parent or patient reported the symptom was drug related and there was no physician documentation of an alternative etiology for the symptom; or (3) no documented history of a symptom prior to hospital admission and no documented explanation for the symptom; or (4) a history of a symptom, but during the medication course the symptom increased in frequency and/or severity. Each ADR was classified as mild or severe by the study's child and adolescent psychiatrist. The operational definition for a severe ADR was the need for physician evaluation and intervention.
Descriptive statistics were calculated for patients' characteristics. The following values were used to reduce the influence of outliers: Number of PRN doses >30=30, number of time outs/seclusions >20=20, number of holds >5=5, number of restraints >3=3, number of ADRs >10=10, and LOS >31 days=31.
Statistical models were developed for all outcomes and modeled as a function of age, sex (Dean et al. 2007), measure of severity (multiple diagnoses), number of psychotropic medications (Jensen et al. 2007), admit GAF (American Psychiatric Association 2000), and the primary predictor of interest, combined phenotype. Battery reduction was performed to reduce the dimensionality of the nine diagnosis categories to produce valid statistical models (D'Agostino et al.1995). In all regression models, predictors were transformed using either logs (admit GAF score) or natural splines (age, number of psychotropic medications, and predicted phenotypes) to relax the linearity assumption. Nonlinear terms were removed if the p value was greater than 0.20.
For the primary statistical analysis using combined phenotype as well as the subset analyses using individual CYP2D6- and CYP2C19-predicted phenotypes, multivariable regression models were fit to examine the effects of the predictors (described above) on the primary efficacy (BIS) and primary tolerability (number of ADRs) outcome variables. With the advice of the study psychiatrist, race was not included as a model predictor because race was not considered a contributing factor to clinical course variability. Multivariable models were also used to examine the same model predictors on the secondary outcomes (total number of PRN psychotropic doses during admission, change in GAF score, and LOS).
A standard formula using a one-way analysis of variance (ANOVA) with a planned contrast was used to determine sample size and suggested 376 subjects would give investigators 81% power to detect a linear trend among the combined phenotype groups with moderate effect sizes between the combined poor metabolizer (C-PM) and combined extensive metabolizer (C-EM) groups. To maximize the potential of detecting a linear trend among the groups, all patients in the C-PM, combined intermediate metabolizer (C-IM), and combined ultrarapid metabolizer (C-UM) groups (n=268) and 279 patients from the C-EM group (a number roughly equal to the total of the other three groups) were included for eligibility review (total n=547; Fig. 1). The subset of C-EM group patients' charts was selected by the study statisticians using a list generated by a simple randomization and presented along with the other charts to the raters with the genotype information concealed.
The primary analysis was performed on patients meeting inclusion/exclusion criteria who received psychotropics metabolized by either CYP2D6 or CYP2C19. All statistical analyses were conducted using S-Plus version 8.0 (TIBCO Software INC, Palo Alto, CA). The Hmisc and Design Libraries (Harrell 2009) were used for statistical modeling.
As part of routine clinical care, 1322 (78%) of the 1701 unique patients admitted to the inpatient psychiatric service between January 1, 2005, and September 30, 2005, were genotyped for both CYP2D6 and CYP2C19. A total of 369 children met study inclusion/exclusion criteria; the primary analyses were performed on the subset of 279 patients (ages 3–18 years, median age 13 years) who received psychotropic medications metabolized by CYP2D6 or CYP2C19 or both (Fig. 1). The clinical characteristics are shown in Table 2. The majority of patients were on either one (36%) or two (29%) psychotropics. The ten most commonly used psychotropic medications metabolized by CYP2D6 and/or CYP2C19 are listed in Table 3. Battery reduction, used to reduce dimensionality, indicated that four diagnostic categories (psychotic disorders, anxiety disorders, impulse control disorders, pervasive developmental disorder) accounted for 77% of the original variance. These were selected as predictors in the subsequent analyses.
At least one ADR was experienced by 50% of patients. The most common were sleep disturbances (28%), gastrointestinal symptoms (15%), headache (13%), and difficulty concentrating (8%); all were rated as mild. Severe ADRs were less common and included mood change (8%), dizziness (4%), extrapyramidal symptoms (4%), aggressive behavior (2%), rash (1%), and shortness of breath (1%). There were single episodes of neuroleptic malignant syndrome and convulsions.
One or more PRN psychotropics were used by 60%. The most common PRN medication was diphenhydramine (92%) followed by (in descending frequency ranging between 11% and 3%) olanzapine, risperidone, quetiapine, ziprasidone, melatonin, and haloperidol.
Individuals from the subset of 279 patients were placed into one of four possible phenotype groups using the algorithm described in Table 1, stage 2. Slightly over half (52%) of the patients had a predicted C-EM phenotype, 22% had a C-IM phenotype, 20% had a C-PM phenotype, and 6% had a C-UM phenotype . There were no differences among combined phenotype groups in sex, diagnostic categories, admit GAF score, or the distribution of the total number of scheduled psychotropic drugs used. There was greater use of scheduled psychotropic medications classified as CYP2D6 substrates in the C-UM and C-PM phenotype groups (both 100%) than the C-IM phenotype group (87%) and the C-EM phenotype group (89%) (p=0.03). There was no difference among the four combined phenotype groups in the use of scheduled psychotropic medications classified as CYP2C19 substrates.
For each patient, the maximal dose of each scheduled psychotropic medication was identified. Among the four combined phenotype groups, there was no difference in the maximal dose for the four most commonly used CYP2D6 substrate psychotropics (risperidone, fluoxetine, aripiprazole, atomoxetine) and the two most commonly used CYP2C19 substrate psychotropics (fluoxetine, escitalopram).
After adjustment for model predictors, there was a significant relationship between the model and the study's primary efficacy outcome variable, BIS (p<0.0001, r2=0.30). Five of the factors were independently significant—age, sex, number of psychotropics, impulse control disorder diagnosis, and the combined phenotype.
There was a statistically significant relationship between combined phenotype and BIS (p=0.01, Table 4) when adjusted for other predictors in the model. Patients in the C-PM group had lower BIS (higher efficacy), whereas patients in the C-UM group had the highest BIS (lowest efficacy). Although none of the secondary outcome variables reached significance, the total number of PRN medication doses (p=0.14) showed a steady increase across groups (from C-PM to C-UM) that paralleled the effect detected by the BIS. There was no difference among groups in change in GAF scores (p=0.90).
As a secondary analysis, the effect of each gene's predicted metabolizing phenotype was examined (Table 5). A relationship between CYP2D6-predicted metabolizing phenotype and BIS was noted (p=0.01). In contrast, no relationship was detected between CYP2C19-predicted metabolizing phenotype and BIS (p=0.57).
A total of 90 out of 369 patients met all inclusion/exclusion criteria but did not receive psychotropics dependent on CYP2D6 or CYP2C19 (Fig. 1). These 90 patients were analyzed using statistical methods identical to the study cohort to determine if a similar linear relationship between BIS and genotype existed in patients taking psychotropics not metabolized by either CYP2D6 or CYP2C19. No difference among the four combined phenotype groups was detected (p=0.42).
There was a significant relationship between combined predicted phenotype and the number of ADRs (primary tolerability variable, p=0.03; Table 4) when adjusted for other predictors in the model. As a secondary analysis, a relationship between CYP2C19-predicted metabolizing phenotype and number of ADRs was noted (p=0.01; Table 5). An association between CYP2C19-predicted metabolizing phenotype and the type of ADRs (severe vs. mild vs. none) was also seen using the nonparametric Kruskal–Wallis test (p=0.04).
This is the first study to detect a pharmacogenetic association between drug-metabolizing enzyme genotypes and short-term psychotropic medication efficacy and tolerability in hospitalized children. As metabolizing capability increased across the four combined metabolizing phenotype groups, the BIS increased steadily, indicating decreased drug efficacy (p=0.01), and the number of adverse drug reactions (ADRs) decreased, indicating less drug toxicity (p=0.03). When only the CYP2C19-predicted phenotype was considered, there was a significant increase in the total number (p=0.01) and severe ADRs (p=0.04) as metabolizing ability decreased, indicating increased drug toxicity.
Aggressive behavior is a common nonspecific symptom in many psychiatric disorders (Masters et al. 2002). In one study, 33% of youths in an inpatient psychiatric facility were involved in an assault on hospital staff (Ryan et al. 2004). Because patients' aggressive behaviors put them, their families, and their inpatient caretakers at great risk, limiting and decreasing aggressive behavior is a key objective of patients' hospitalization3 (Foster et al. 2007; Jensen et al. 2007) This study used BIS as the measure of short-term drug efficacy because as medication efficacy improves, the BIS should decrease.
Clinicians at the study site did not rate aggression in the medical records using a standard scale, which created the need to use a proxy, BIS, as the measure of short-term drug efficacy. The BIS is based on the objective, reliable, and verifiable interventions of time outs, seclusions, therapeutic holds, and physical restraints. By hospital policy, these four interventions are accepted clinical interventions for aggression against others or self, and each use of any one of these must be documented in the patient's hospital chart. As such, retrospective identification of any of these four interventions is a reliable, identifiable proxy for patient aggression. In general, a child requiring a time out was considered to have less severe aggression than a child requiring a therapeutic hold or physical restraint making it unreasonable to simply sum the total number of behavioral interventions for each child. A search of the literature and consultation with the study's child adolescent psychiatrist did not provide evidence for assigning differential weights to BI1, BI2, and BI3. Therefore, a composite BIS was calculated using a weight-free index approach,
that does not change under different weights applied to the individual BI items (Elston 1963). Internal consistency of the individual BI items demonstrated a Cronbach alpha coefficient of 0.72, which is considered acceptable (Nunnaly 1978) when combining items. As scheduled psychotropics effectively treat the underlying diagnoses and associated aggression, the need for behavioral interventions should decrease as clinical condition improves. Therefore, BIS is a logical measure of short-term drug efficacy: As medication efficacy improves, BIS should decrease.
The study found that patients with C-UM phenotypes had the highest BIS and fewest ADRs; as metabolizing ability decreased across combined phenotype groups, the BIS decreased (improved efficacy) and ADRs increased (greater toxicity). Two potential explanations include: (1) An underlying relationship between CYP2D6/CYP2C19 genotype and psychiatric disease severity or (2) a genotype–pharmacokinetic association where differential drug exposure occurs across combined phenotype groups. Two findings from this study indicate a genotype–pharmacokinetic relationship is more likely than a genotype–disease severity phenotype: (1) The 90 patients who did not receive CYP2D6- or CYP2C19-dependent psychotropic medications showed no significant differences in BIS across the four combined phenotype groups and (2) the maximal dose for six of the most commonly used CYP2D6/CYP2C19-dependent psychotropics was the same across combined phenotype groups.
Specifically, inadequate drug efficacy, as measured by the BIS, coupled with low ADR rates particularly in the C-UM patients strongly suggests these patients are being underdosed. The reverse is true for C-PM patients. These results are similar to analyses of pooled data from children and adolescents who were treated with atomoxetine in any of four acute, double-blind, and placebo-controlled studies (Michelson et al. 2007). The investigators found that patients who had genotypes consistent with CYP2D6-predicted poor metabolism had greater symptom reduction when compared to CYP2D6 EM patients at the cost of a modestly greater increase in some adverse events. The investigators noted that more EM patients discontinued because of lack of efficacy than PMs who discontinued due to adverse events.
This is the first study to identify a significant relationship between the total number of ADRs (p=0.01) and severe ADRs (p=0.04) across CYP2C19 genotypes for the subgroups of patients taking psychotropics metabolized through the CYP2C19 pathway. While not statistically significant, a trend of increasing ADRs from CYP2D6 UMs to CYP2D6 PMs was demonstrated. Similar to this study, negative CYP2D6 results have been noted by others (Binder and Holsboer 2006). Yet, some studies (using smaller sample sizes, different methodologies, and patient populations) detected a relationship between the CYP2D6 genotype and the incidence of ADRs to medications metabolized through the CYP2D6 pathway (Chen et al. 1996; de Leon et al. 1998; Chou et al. 2000; Kirchheiner et al. 2001; Mrazek 2006).
In this study, CYP2D6 substrates were mostly atypical antipsychotics and selective serotonin reuptake inhibitor (SSRI) antidepressants rather than typical antipsychotics or tricyclic antidepressants. The most commonly occurring ADRs with these psychotropics are weight gain and metabolic effects (Correll and Carlson 2006) and were likely underreported in this study due to insufficient LOS. Common, yet highly subjective, ADRs such as akathisia were not detected frequently, possibly due to underreporting or inadequate documentation in the medical record.
This study is limited by its retrospective design, the number of different medications patients received, the lack of long-term follow up, the lack of standardized aggression and ADR scale use, the exclusion of certain “mixed” genotype combinations (Fig.1), and the limited number of CYP2D6 and CYP2C19 alleles tested. Retrospective studies have the potential for more bias than prospective studies. We attempted to minimize bias by using the rater-blinded methodology. The lack of differences across combined phenotype groups in both the medications used and their maximal drug doses implies that medication selection and short-term usage patterns were similar across the groups. This reduces, but does not eliminate, the potential bias introduced by lack of a prospective protocol-driven medication selection and titration schedule approach. Clearly, prospective, longitudinal studies in pediatric patients would contribute to further understanding of the combined phenotype–clinical response relationships. Although ideal, standardized ADR and aggression scales are not commonly used in the routine inpatient psychiatric clinical care setting and were not available for this study.
The CYP2D6 and CYP2C19 genotype data used in this study presented the results of a genetic analysis performed for clinical care not research purposes. As such, the panel of alleles tested in 2005 was not comprehensive but represented a balance between established clinical relevance of allelic variation (Bradford 2002), need for results within a 2-day window, and cost. Future prospective studies should include a larger panel of CYP2D6 alleles (e.g., the functional allele CYP2D6*2 (Gaedigk et al. 2003) and the reduced-function alleles that occur more commonly in non-Caucasian populations, for example CYP2D6*10 (Luo et al. 2005; Sistonen et al. 2007), CYP2D6*17, CYP2D6*41 (Cai et al. 2006; Gaedigk et al. 2007), along with CYP2C19*17, a recently discovered polymorphism associated with CYP2C19 ultrarapid metabolism (Sim et al. 2006; Rudberg et al. 2008)).
Identifying factors underlying the variability in drug efficacy or tolerability is a key component for optimizing the patient's response to therapy. This is the first study to demonstrate that, for children hospitalized for psychiatric conditions, CYP2D6 and CYP2C19 genetic variation contributes to psychotropic medication's short-term efficacy against aggression and the occurrence of ADRs. Consideration of a patient's genotype at the onset of a psychiatric hospitalization could play a significant role in personalizing and improving subsequent therapy. Prospective longitudinal studies are necessary to better inform how to optimally incorporate this genetic information into the medical management of patients with aggression. In summary, this study indicates that CYP2D6 and CYP2C19 genotypes are important in the clinical care of children with psychiatric diagnoses requiring medications that are metabolized through these two enzyme pathways.
The manuscript's authors disclose the following corporate/commercial relations that might pose a conflict of interest: Cynthia A. Prows is a consultant for PTC Therapeutics; Sanjeev Pathak currently works at AstraZeneca Pharmaceuticals, Wilmington, DE, and during study planning, implementation, and results interpretation was a psychiatrist at Cincinnati Children's Hospital Medical Center; Alexander A. Vinks is a consultant for NPS pharmaceuticals, has grant support from Roche Laboratories, and has a laboratory contract with Isotechnika, Inc.; Tracy A. Glauser is a consultant for UCB Pharma, Ortho-McNeil Neurologics, Eisai Pharmaceuticals, and Jazz Pharmaceutical; Todd G. Nick, Shannon N. Saldana, Kejian Zhang, Chunyan Liu, and Zachary S. Daniels have no conflicts of interest or financial ties to disclose.
Supported in part by NIH grant 5U10H0037249=10(SNS, AAV, TA6).
The authors would like to thank: Richard Wenstrup, M.D., for his advice during the development and implementation of the study; Marie Malgaz, Kasia Bryc, B.S., and Steve Fordyce for their assistance with identifying eligible patients; Debbie Baker, Luke Botting, B.A., and Carol Hetteberg, M.S.N., for their tireless and careful chart reviews; and Thomas Boat, M.D., Michael Sorter, M.D., Robert Kowatch, M.D., Ph.D., and Gregory Grabowski, M.D., for their critical and insightful reviews.