Skip to main content
  • Research article
  • Open access
  • Published:

Functional phenotyping of the CYP2D6 probe drug codeine in the horse



In humans, the drug metabolizing enzyme CYP2D6 is highly polymorphic resulting in substantial differences in the metabolism of drugs including anti-arrhythmics, neuroleptics, and opioids. The objective of this study was to phenotype a population of 100 horses from five different breeds and assess differences in the metabolic activity of the equine CYP2D6 homolog using codeine as a probe drug. Administration of a probe drug is a common method used for patient phenotyping in human medicine, whereby the ratio of parent drug to metabolite (metabolic ratio, MR) can be used to compare relative enzyme function between individuals. A single oral dose of codeine (0.6 mg/kg) was administered and plasma concentrations of codeine and its metabolites were determined using liquid chromatography mass spectrometry. The MR of codeine O-demethylation [(codeine)/(morphine + morphine-3-glucuronide + morphine-6-glucuronide)] was determined using the area under the plasma concentration-time curve extrapolated from time zero to infinity (AUC0-∞) for each analyte and used to group horses into predicted phenotypes (high-, moderate-, and low-MR).


The MR of codeine O-demethylation ranged from 0.002 to 0.147 (median 0.018) among all horses. No significant difference in MR was observed between breeds, age, or sex. Of the 100 horses, 11 were classified as high-MR, 72 moderate-MR, and 17 low-MR. Codeine AUC0-∞ and O-demethylation MR were significantly different (p < 0.05) between all three groups. The mean ± SD MR was 0.089 ± 0.027, 0.022 ± 0.011, and 0.0095 ± 0.001 for high-, moderate-, and low-MR groups, respectively. The AUC for the morphine metabolites morphine-3-glucuronide and morphine-6-glucuronide were significantly different between high-and low-MR groups (p < 0.004 and p < 0.006).


The MR calculated from plasma following codeine administration allowed for classification of horses into metabolic phenotypes within a large population. The range of codeine metabolism observed among horses suggests the presence of genetic polymorphisms in CYP2D82 of which codeine is a known substrate. Additional studies including CYP2D82 genotyping of high- and low-MR individuals are necessary to determine the presence of CYP2D polymorphisms and their functional implications with respect to the metabolism of therapeutics.


Cytochrome 2D6 (CYP2D6) belongs to the cytochrome P450 (CYP450) superfamily of enzymes, which is involved in the metabolism of drugs and toxicants, as well as the metabolism and synthesis of endogenous substrates such as eicosanoids and cholesterol [1]. While constituting only 2% of total CYP450 content in the liver, CYP2D6 is responsible for metabolizing 25–30% of drugs including antiarrhythmics, antidepressants, anti-cancer drugs, neuroleptics, and opioids [2, 3]. In addition, human CYP2D6 displays high inter-individual variability, with a 1000-fold difference in enzyme activity, largely due to the prevalence of polymorphisms within the CYP2D6 gene [4]. Clinical implications of these polymorphisms are well documented and can result in lack of therapeutic efficacy due to rapid drug metabolism or toxicity and even death in slow metabolizers due to drug bioaccumulation [5,6,7]. This large variation in CYP2D6 activity has led to the classification of three metabolizer phenotypes; poor (PM), extensive (EM), and ultra-rapid (UM) [8]. Administration of a CYP2D6-selective probe drug and subsequent calculation of the metabolic ratio (MR) of parent drug to metabolite from plasma or urine samples is a common method for the characterization of CYP2D6 phenotypes. Correspondingly, poor, extensive, and ultra-rapid metabolizers are described as having high-, moderate-, and low-MR, respectively. Individuals classified as poor or ultra-rapid metabolizers, based on the MR, are at risk for adverse drug reactions or lack of therapeutic response, respectively when administered CYP2D6 substrates [6, 9, 10].

In the horse, reports of polymorphisms in drug metabolizing enzymes such as CYP450s, and their effect on therapeutic outcomes are sparse. However, previous reports have identified mutations in an equine CYP2D ortholog, CYP2D50, as well as phenotypic differences in the metabolism of CYP2D50 substrates, suggesting a range of enzyme function [11, 12]. While humans possess one active CYP2D gene, other species, such as the marmoset, rat and mouse, have multiple CYP2D genes [13,14,15]. In the horse, six CYP2D genes have been identified, however, the expression and activity for all has not yet been established [16]. Recently, another CYP2D6 ortholog, CYP2D82 was found to catalyze the O-demethylation of codeine to morphine in vitro [17]. The investigators concluded that codeine was a selective substrate for CYP2D82, as no metabolites were detected in incubations with other expressed equine CYP2D enzymes [17]. Additionally, in vivo studies have demonstrated that codeine is metabolized to the same five main metabolites found in humans, including the CYP2D-generated metabolite morphine (Fig. 1) [18].

Fig. 1
figure 1

Major codeine metabolites in the human and enzymes responsible for metabolite formation

While codeine is not commonly used as a probe drug for CYP2D6 in humans, there are numerous reports demonstrating that genetic polymorphisms in CYP2D6 can lead to differing concentrations of the active metabolite morphine between individuals [7, 19]. Differences in morphine concentrations can result in a large degree of individual variability with respect to analgesia, following codeine administration [7, 20]. Furthermore, studies have shown that the MR for codeine O-demethylation in humans is highly correlated with CYP2D6 genotype and that codeine is as effective as the established probe drug, desbrisoquine, in CYP2D6 phenotyping [21, 22]. Given previous evidence of variable CYP2D metabolism in the horse, and the specificity of codeine for the novel equine CYP2D enzyme, CYP2D82, we hypothesized that codeine will be an effective probe drug for the phenotyping of horses. Due to the large number of drugs metabolized by CYP2D6 in the human, if members of the CYP2D family in the horse prove to be as important in the metabolism of therapeutic substances, significant variations in CYP2D enzyme metabolism may result in drug toxicity and lack of efficacy in the horse as is observed in humans [6, 9, 10]. In addition to the clinical implications, decreases in clearance as a result of altered metabolism of CYP2D substrates could lead to positive regulatory findings, despite adherence to established withdrawal times, in performance horses. The goal of the currently reported study is to add to previously published reports describing individual differences in CYP2D metabolic activity between horses.


Quality control samples prepared in replicates (n = 6) for codeine, C6G, norcodeine, morphine, M3G, and M6G were analyzed to determine the intra- and inter-day precision and accuracy of the assay (Table 1). The concentration-response relationship (relationship between calibrators and the LC-MS/MS instrument response) for all analytes was linear (r = 0.99). The limit of quantitation for the assay was 0.1 ng/ mL and the limit of detection approximately 0.03 ng/mL for all analytes.

Table 1 Intra- and interday accuracy and precision (values for liquid chromatography–tandem mass spectrometry (LC-MS/MS) analysis of codeine and its metabolites, norcodeine, codeine-6-glucuronide (C6G), morphine, morphine-6-glucuronide (M3G), and morphine-6-glucuronide (M6G) in equine plasma

The normal probability plot of log-transformed MRs showed a nonlinear plot indicative of a non-normal distribution [23] (Fig. 2). Visual analysis of the normal probability plot indicated steep break points at both ends corresponding to anti-modes of − 2.1 and − 1.2 using methods previously described [11, 24]. The frequency distribution histogram of log-transformed metabolic ratios confirmed the presence of anti-modes at − 2.1 and − 1.2 by visual analysis (Fig. 3). Anti-modes determined from these plots were used to group horses into predicted phenotypes with 17 classified as low-MR, 11 as high-MR, and the remaining 72 as moderate-MR. The low-MR group consisted of 12 mares and 5 geldings including 11 Thoroughbreds, two Quarter horses, and four warmbloods aged 2 to 20 years. In the high-MR group there were four mares and seven geldings, ten of which were Thoroughbreds and one Standardbred, aged 3–11 years.

Fig. 2
figure 2

Normal probability plot of logarithmic metabolic ratio [log (AUC0-∞ codeine) / (AUC0-∞ morphine + AUC0-∞ M3G + AUC0-∞ M6G)] in 100 horses after a single oral administration of codeine (0.6 mg/kg)

Fig. 3
figure 3

Frequency distribution histogram of log metabolic ratio [log (AUC0-∞ codeine) / (AUC0-∞ morphine + AUC0-∞ M3G + AUC0-∞ M6G)] in 100 horses after a single oral administration of codeine (0.6 mg/kg) grouped by predicted phenotype. Dashed vertical lines denote anti-modes

The mean (±SD), median, and range of MR for O-demethylation for each breed and age group are presented in Table 2. No significant differences were found for the metabolic ratios between breed, age, or sex. The AUC (mean ± SD and range) for codeine, C6G, morphine, M3G, and M6G, the metabolic ratio, and codeine elimination half-life for the predicted phenotype groups are listed in Table 3 The Cmax for codeine, C6G, morphine, M3G, and M6G ranged from 4.3–443.1, 18.1–320.9, 0.7–8.3, 113.2–1024.5, and 7.0–87.2 ng/mL, respectively. Metabolic ratios for codeine O-demethylation ranged from 0.002 to 0.147 with a median of 0.018. The AUC for codeine and the MR for O-demethylation were found to be significantly different (p < 0.05) between all three predicted phenotypes. Additionally, the AUC for M3G and M6G and the MR for glucuronidation were significantly different (p < 0.05) between high- and low-MR groups.

Table 2 Metabolic ratio for codeine O-demethylation following administration of a single oral dose to a population of horses of differing breeds and ages. The metabolic ratio calculated as (codeine AUC0 − ∞)/(morphine AUC0 − ∞ + M3G AUC0 − ∞ + M6G AUC0 − ∞)
Table 3 Pharmacokinetic values (mean ± SD) for codeine and codeine 6-glucuronide (C6G), morphine, morphine 3-glucuronide (M3G), and morphine 6-glucuronide (M6G) after a single oral administration of codeine (0.6 mg/kg) to horses. Parameters are sorted by predicted phenotypic groups

Average plasma concentrations of codeine and its metabolites were plotted for each phenotypic group (Fig. 4a-f). The average plasma concentrations for codeine and C6G were highest for the high-MR group, followed by the moderate-MR, and lastly the low-MR group. Morphine plasma concentrations were similar for moderate- and low-MR groups, but lower for the high-MR group. Morphine-3-glucuronide plasma concentrations were highest in the low-MR individuals and lowest in the high-MR. The average Cmax for M6G was 14.84, 25.28, and 25.96 ng/mL for high-, moderate-, and low-MR groups.

Fig. 4
figure 4

Mean plasma-concentration time curves for codeine (a), morphine (b), C6G (c), M3G (d), norcodeine (e), and M6G (f) based on predicted phenotype in horses following a single oral dose of codeine (0.6 mg/kg)


Similar to what has been described in humans, there are reports suggestive of varying metabolic capabilities between individual horses [12]. The goal of the current study was to add to the limited knowledge regarding differences in drug metabolism in this species, specifically that mediated by CYP2D82, the equine homolog to the highly polymorphic, CYP2D6. To that end, codeine, which is metabolized to morphine by CYP2D82 [17], was used as a probe drug to screen for differing CYP2D82 metabolic phenotypes in a population of 100 horses. The pharmacokinetics and metabolic profile of codeine has been described previously [18]. In agreement with the previous report, the primary metabolite produced in the current study was the morphine metabolite M3G, followed by C6G, with morphine and M6G present at lower concentrations.

As described previously, the AUC for codeine and its metabolites were determined and subsequently used to calculate a MR (codeine: metabolite) to classify horses into metabolic phenotypes [23]. The codeine t1/2 was not significantly different between the predicted metabolic phenotypes and was in accordance with previous reports. Interestingly, the codeine AUC0-∞ for moderate-MR individuals was lower than previously reported [18]. While the reason for this is not immediately evident, possible explanations include inter-individual variability within study animals coupled with the smaller sample size in the previous report. Additionally, differences in the pharmacokinetic modeling approach between the studies might account for the observed differences in codeine AUC0-∞. The AUC0-∞ for codeine was significantly different between predicted metabolic groups (high-, moderate-, and low-MR) with high-MR > moderate-MR > low-MR. Although the AUC0-∞ for both C6G and morphine followed the expected pattern based on phenotype groupings (low-MR > moderate-MR > high-MR), this difference was not statistically significant. When evaluating the AUCs from the different metabolic groups, it is important to consider that the AUC can be affected by other factors, notably the extent of absorption and distribution. For this reason, the MR is often considered a better indicator of metabolic capabilities and is used routinely in phenotyping studies. The MR for determination of CYP2D82 activity (O-demethylation) was calculated as the AUC0-∞ codeine / (AUC0-∞ morphine + AUC0-∞ M3G + AUC0-∞ M6G) as described previously [6, 25]. Production of M3G and M6G from codeine first requires biotransformation to morphine, therefore the glucuronidated metabolites were included as products of the O-demethylation metabolic pathway in the determination of the MR [25].

The MR frequency distribution in the current study is best described as unimodal, however, the absence of an obvious multi-modal distribution upon inspection of the histogram is not definitive evidence that such a distribution does not exist. Although they can be valuable tools for determination of a non-normal or multi-modal distribution and are commonly utilized in phenotyping studies, the shape of a histogram is dependent on sample size and can therefore be misleading. A normal probability plot is another, often times, more definitive way of assessing normality. In the current study, visual analysis of the normal probability plot of the log MR representing O-demethylation, revealed breakpoints suggestive of three phenotypic groups (high-, moderate-, and low-MR) within the population of horses tested. Furthermore, apparent anti-modes on both sides of the frequency distribution histogram suggest the existence of high- and low-MR groups. In the current analysis of 100 horses, the incidence of predicted CYP2D82 high- and low-MR individuals as classified using the probe drug, codeine, was determined to be 11 and 17%, respectively. Comparably, within the Caucasian population for CYP2D6, 5–10% are PMs (high-MR) and 3–5% are UMs (low-MR) [26]. While, notably, this is the first study of its kind in horses and sample size is limited compared to human studies, results of the current study suggest the existence of differing metabolic phenotypes in horses. Similar to what is observed in humans, there is a high likelihood that differing metabolic capabilities may influence therapeutic responses to CYP2D substrates [27, 28].

In humans, the prevalence of CYP2D6 PMs and UMs can vary significantly among different ethnicities. For example, while 5–10% of the Caucasian population are PMs, only 1% of Chinese are PMs [26, 29]. Similar to the interethnic variation observed in CYP2D6 phenotypes, it is plausible that the prevalence of metabolic phenotypes may vary between horse breeds. The present study evaluated a total of one hundred horses from five different breeds to assess the overall prevalence of polymorphisms within a mixed population. Inter-individual variation in MR was observed within the breeds, while most horses fell within the category of moderate-MR, there were high- and low-MR individuals from several different breeds. While no significant differences were observed in the MR between the breeds, it is important to note that the majority of horses studied were of one breed (Thoroughbreds). This, combined with the presumed relatively low frequency of high- and low-MR individuals, may confound the identification of differences between breeds. Additional studies including larger numbers of horses of other breeds would be needed to determine whether phenotype distribution varies between breeds.

In the current study, the primary focus was on identifying differences in the biotransformation of codeine to morphine that could be attributed to altered CYP2D82 activity. Although it has been well established, based on screening of recombinant equine CYPs, that the biotransformation of codeine to morphine is catalyzed by CYP2D82 [17], it is important to note that glucuronidation of the parent drug (codeine) and N-demethylation to norcodeine also contributes to the overall clearance of this drug in horses [18]. These alternate metabolic pathways may explain the similar codeine terminal half-life between metabolic phenotypes observed in the present study. A decrease in O-demethylation capabilities may be compensated for by the codeine glucuronidation and N-demethylation pathways leading to comparable terminal half-lives between the high- and low-MR groups. Inspection of plasma concentration-time curves for the metabolites supports this theory with slightly higher concentrations of C6G and norcodeine noted in the high-MR group compared to low-MR group.

While the current study theorizes that the higher conversion of codeine to C6G may be a possible compensatory clearance pathway, it should also be noted that studies have identified polymorphisms in genes that code for UGT enzymes in humans [30, 31]. Although the UGT enzyme/s responsible for codeine glucuronidation in the horse have not yet been identified, in humans this reaction is catalyzed by UGT2B7 [32]. The gene coding for UGT2B7 in humans has been shown to be polymorphic and, while there are several studies describing the implications of these polymorphisms on morphine glucuronidation and elimination, there are only a few describing effects on overall codeine clearance (both generation of C6G and morphine glucuronides) and those results appear mixed [33, 34]. Nonetheless, the potential for UGT2B7 polymorphisms complicates classification of horses into CYP2D82 metabolic phenotypes using codeine as a probe drug. Studies conducted in our laboratory, have shown that similar to humans, morphine is glucuronidated by equine UGT2B7 to M3G (unpublished data). If, as in humans, generation of C6G from codeine is also catalyzed by UGT2B7, and this enzyme is polymorphic in horses, it is plausible that a presumed CYP2D82 high-MR individual may actually be a UGT2B7 UM. In the current study, there does not appear to be a difference in the morphine: M3G MR (calculated as morphine AUC0-∞ /M3G AUC0-∞), suggesting that the rate-limiting step is production of morphine, which is catalyzed by CYP2D82 and not glucuronidation. Furthermore, the MRs for the codeine glucuronidation pathway (codeine: C6G) are comparable between proposed phenotype groups suggesting that differences between groups are not related to genetic mutations in genes coding for UGT enzymes. However, additional metabolism studies as well as genotyping studies focused on CYP2D82 and UGTs are warranted to definitively prove this theory.

While screening of recombinant equine CYPs have demonstrated that CYP2D82 catalyzes the biotransformation of codeine to morphine, it is important to note that the reported in vitro Km of CYP2D82 was greater than the maximum observed codeine plasma concentrations in the presented report [17]. However, plasma concentrations of codeine may differ significantly from hepatic concentrations, where the majority of CYPs reside. In humans, hepatic concentrations of codeine have been found to be up to 100-times greater than in plasma [35]. Furthermore, differences in in vitro and in vivo conditions may affect enzyme kinetics, such as the less efficient coupling of CYPs with the coenzyme P450 reductase in in vitro experiments. Given this limitation, potential involvement of other phase 1 drug metabolizing enzymes such as flavin-containing monooxygenases (FMOs), or additional unidentified CYPs in the metabolism of codeine in the horse warrants further exploration.

With the reduced cost, and development of rapid DNA sequencing methods, modern phenotyping studies frequently include genotyping to characterize polymorphic alleles within the population and assign activity scores based on in vivo activity [36]. In order to determine a causal link between allele function and metabolic phenotype, genotyping of suspected polymorphic alleles is necessary. However, prior to widespread genetic sequencing and characterization of polymorphic CYP alleles in humans, phenotyping studies used methods of visual analysis of normal probability plots and frequency histograms, as applied in the presented study, to describe proposed metabolic phenotypes [24, 37]. Ongoing studies involving whole genome sequencing are currently underway to determine whether genetic polymorphisms are present in the CYP2D82 gene of horses identified as high- or low-MR in this report. Previous sequencing of the equine CYP2D50 gene identified several exonic SNPs that may affect enzyme activity and metabolism of CYP2D substrates [11, 12]. It is important to note that other factors may influence the metabolism of drugs in vivo, including epigenetic changes, polymorphisms within other genes such as transporters involved in drug absorption [38], age, and environment, such as diet and toxicant exposure [39]. While phenotyping is a valuable tool to assess in vivo metabolism of CYP2D substrates, genotyping horses identified as high- or low-MRs would determine whether observed differences in metabolism arise from altered CYP2D82 activity as opposed to nongenetic or epigenetic differences.


The administration of codeine, a known substrate of equine CYP2D82, as a probe drug suggests three metabolic phenotypes (high-, moderate-, and low-MR) within a mixed-breed population of horses, as evidenced by significant differences between the MR representing O-demethylation. While the study of polymorphisms in genes coding for metabolic enzymes is very much in its infancy in horses, the potential therapeutic implications, as have been described in human medicine, warrant further study including additional drug metabolism, phenotyping and genotyping studies.



One hundred horses (2 American Paint Horses, 14 Quarter Horses, 4 Standardbreds, 66 Thoroughbreds, and 14 Warmbloods) were used in this study. Of the 100 horses, 52 were geldings and 48 were mares. Horses ranged in age from 2 to 28 years (11 ± 7) with a mean ± SD weight of 530.6 ± 51.0 kg. The horses used in this study included University-owned research horses and horses from a private Northern California Thoroughbred training facility. The horses were housed in a mixture of pasture, indoor stalls, and stalls with outside runs attached, depending on the facility. The horses were fed twice daily, consisting of AM and PM feedings of grass or alfalfa hay, according to their regular schedule. Prior to participation in the study, all horses were determined healthy by physical examination, complete blood count, and a serum biochemistry panel. The horses did not receive any medications for at least 2 weeks prior to the start of the study. Additionally, no anthelminthics were administered for at least 8 weeks prior to codeine administration. Following the conclusion of the study, all animals were returned to the University research herd or to training. This study was approved by the Institutional Animal Care and Use Committee of the University of California, Davis (#20319).

Drug administration and instrumentation

Codeine sulfate tablets (30 mg tablets; West-Ward Pharmaceuticals Corp., NJ) were suspended in 6 mL of water in a 20 mL dosing syringe and administered as a single dose of 0.6 mg/kg into the oral cavity. Manual shaking of the dosing syringe was performed to ensure adequate suspension prior to administration. The dosing syringe was rinsed with an additional 6 mL of water and delivered to the oral cavity immediately after dosing. Horses were not fasted prior to drug administration and were weighed immediately prior to drug administration.

Sample collection

Blood samples were collected by direct jugular venipuncture into ethylenediaminetetraacetic acid (EDTA) containing blood tubes (Covidien, MA, USA) immediately before drug administration (time 0) and at 0.5, 2, 4, 8, 24, 48, and 72 h post administration. Samples were placed on ice until centrifugation at 3000 g for 10 min. Plasma was transferred to cryovials (Phenix Research Products, Chandler, NC, USA) and stored at − 20 °C until analysis.

Codeine plasma concentration determination

Quantitative analysis of plasma samples was performed using liquid chromatography tandem mass spectrometry (LC-MS/MS). Plasma samples, calibrators, and negative controls were all prepared as previously described [18]. The previously published method was adapted for the TSQ Altis triple quadrupole mass spectrometer coupled with a Vanquish liquid chromatography system (Thermo Scientific, San Jose, CA).

The analyte concentrations were measured in equine plasma by LC-MS/MS using positive heated electrospray ionization (HESI(+)). Quantitative analysis was done on a TSQ Altis triple quadrupole mass spectrometer with a Vanquish liquid chromatography system (Thermo Scientific, San Jose, CA). The parameters were as follows; spray voltage 3500 V, vaporizer temperature 350 °C, sheath gas 50 (arbitrary units), and auxillary gas 10 (arbitrary units). Product masses and collision energies of each analyte were optimized for the TSQ Altis by infusion of the standards into the machine. Chromatography was performed using a Zorbax Eclipse-XDB-Phenyl 3x100mm, 3 μm column (Agilent Technologies) and a linear gradient of ACN in water with 0.2% formic acid and a 0.45 ml/min flow rate. The ACN concentration was held at 5% for 0.5 min initially, then gradually increased to 30% over 4.0 min, then to 90% over 0.01 min, kept at 90% for 0.2 min, and returned to 5% for 2.3 min.

Selective Reaction Monitoring (SRM) was used to detect and quantify the initial precursor ion for codeine (mass to charge ratio 300.1 (m/z)), C6G (476.1 (m/z)), norcodeine (286.1 (m/z)), morphine (286.1 (m/z)), M3G (462.2 (m/z)), M6G (462.1 (m/z)), and the internal standards d6-codeine (306.1 (m/z)), d6-morphine (292.1 (m/z)), and d3-M6G (465.2 (m/z)). Responses for the following product ions were plotted and Quanbrowser software (Thermo Scientific) was used to integrate the peaks at the corresponding retention time: codeine (m/z 165, 215, 199), C6G (m/z 300, 215, 165), norcodeine (m/z 165, 201), morphine (m/z 165, 201), M3G (m/z 286), M6G (m/z 286) and the internal standards d6-codeine (m/z 165, 218), d6-morphine (m/z 201), and d3-M6G (m/z 289). Linear regression using Quanbrowser software was employed to create the calibration curves and quantitate the analytes using a 1/X weighting factor.

Pharmacokinetic analysis

Noncompartmental analysis of plasma codeine, C6G, morphine, M3G and M6G concentrations were performed using Phoenix WinNonlin Version 8.2 (Certara, Princeton, NJ, USA). The analysis included plasma concentrations above the limit of quantitation (LOQ) for the respective analytes. Elimination half-life (t1/2) was calculated with the formula t1/2 = (ln2)/λz, where λz is the elimination rate constant determined as the terminal slope of the plasma concentration-time curve. Area under the curve from time 0 to infinity (AUC0-∞) was determined using the log up-linear down trapezoidal method extrapolated to infinity by dividing the last measured plasma concentration by the terminal slope.

Metabolic ratios (MR) were calculated using the following formulas:

$$ \mathrm{MR}\ \mathrm{codeine}\ \mathrm{O}-\mathrm{demethylation}=\left({\mathrm{AUC}}_{0-\infty }\ \mathrm{codeine}\right)/\left({\mathrm{AUC}}_{0-\infty }\ \mathrm{morphine}+{\mathrm{AUC}}_{0-\infty }\ \mathrm{M}3\mathrm{G}+{\mathrm{AUC}}_{0-\infty }\ \mathrm{M}6\mathrm{G}\right) $$
$$ \mathrm{MR}\ \mathrm{codeine}\ \mathrm{glucuronidation}=\left({\mathrm{AUC}}_{0-\infty }\ \mathrm{codeine}\right)/\left({\mathrm{AUC}}_{0-\infty }\ \mathrm{codeine}-6-\mathrm{glucuronide}\right) $$

The Kolmogorov-Smirnov test was used to assess the normality of the MR data. Metabolic ratios were logarithmically transformed and plotted as a frequency histogram. Probit transformations of the data were made by plotting log MR against the corresponding percent area under the probability curve on probability paper. Percent area under the probability curve was calculated as 100 x i / (n + 1) where i is the rank of the MR for each horse from lowest to highest value and n is the number of samples. Visual inspection of the graphs for breakpoints indicated anti-modes separating phenotypes [24].

Statistical analyses were conducted using commercially available software (Stata/IC 15.1, StataCorp LP, TX). Statistical comparisons of mean differences between the AUC0-∞, MRs, and t1/2 of the phenotypic groups were performed using a one-way analysis of variance (ANOVA). If the overall group effect was significant (p < 0.05), then pairwise comparisons were performed using Scheffe’s method for multiple comparison adjustment.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.



Area under the curve


Cytochrome P450




Extensive metabolizer


Liquid chromatography mass spectrometry/mass spectrometry


Limit of quantitation






Metabolic ratio


Poor metabolizer


Ultra-rapid metabolizer


  1. Nebert DW, Russell DW. Clinical importance of the cytochromes P450. Lancet. 2002;360:1155–62.

    Article  CAS  Google Scholar 

  2. Ingelman-Sundberg M. Genetic polymorphisms of cytochrome P450 2D6 (CYP2D6): clinical consequences, evolutionary aspects and functional diversity. Pharmacogenomics J. 2005;5:6–13.

    Article  CAS  Google Scholar 

  3. Ingelman-Sundberg M, Sim SC, Gomez A, Rodriguez-Antona C. Influence of cytochrome P450 polymorphisms on drug therapies: Pharmacogenetic, pharmacoepigenetic and clinical aspects. Pharmacol Ther. 2007;116:496–526.

    Article  CAS  Google Scholar 

  4. McGraw J, Gerhardt A, Morris TC. Opportunities and obstacles in genotypic prediction of cytochrome P450 phenotypes. Expert Opin Drug Metab Toxicol. 2018;14:659–61.

    Article  CAS  Google Scholar 

  5. Desmeules J, Gascon MP, Dayer P, Magistris M. Impact of environmental and genetic factors on codeine analgesia. Eur J Clin Pharmacol. 1991;41:23–6.

    Article  CAS  Google Scholar 

  6. Kirchheiner J, Nickchen K, Bauer M, Wong ML, Licino J, Roots I, et al. Pharmacogenetics of antidepressants and antipsychotics: the contribution of allelic variations to the phenotype of drug response. Mol Psychiatry. 2004;9:442–73.

    Article  CAS  Google Scholar 

  7. Dean L. Codeine therapy and CYP2D6 genotype. In: Pratt VM, McLeod HL, Rubinstein WS, Scott SA, Dean LC, Kattman BL, Malheiro AJ, editors. Medical Genetics Summaries. Bethesda (MD): National Center for Biotechnology Information (US); 2012. p. 1–9.

    Google Scholar 

  8. Neafsey P, Ginsberg G, Hattis D, Sonawane B. Genetic polymorphism in cytochrome P450 2D6 (CYP2D6): population distribution of CYP2D6 activity. J Toxicol Environ Health B Crit Rev. 2009;12:334–61.

    Article  CAS  Google Scholar 

  9. Gasche Y, Daali Y, Fathi M, Chiappe A, Cottini S, Dayer P, et al. Codeine intoxication associated with ultrarapid CYP2D6 metabolism. N Engl J Med. 2004;351:2827–31.

    Article  CAS  Google Scholar 

  10. Zhou SF. Polymorphism of human cytochrome p450 2D6 and its clinical significance: part I. Clin Pharmacokinet. 2009;48:689–723.

    Article  CAS  Google Scholar 

  11. Corado CR, McKemie DS. Knych HK (2016a) dextromethorphan and debrisoquine metabolism and polymorphism of the gene for cytochrome P450 isozyme 2D50 in thoroughbreds. Am J Vet Res. 2016a;77:1029–35.

    Article  CAS  Google Scholar 

  12. Corado CR, McKemie DS, Young A, Knych HK. Evidence for polymorphism in the cytochrome P450 2D50 gene in horses. J Vet Pharmacol Ther. 2016b;39:245–54.

    Article  CAS  Google Scholar 

  13. Uehara S, Uno Y, Hagihira Y, Murayama M, Shimizu M, Inoue T, et al. Marmoset cytochrome P450 2D8 in livers and small intestines metabolizes typical human P450 2D6 substrates, metoprolol, bufuralol and dextromethorphan. Xenobiotica. 2015;45:766–72.

    Article  CAS  Google Scholar 

  14. Matsunaga E, Zanger UM, Hardwick JP, Gelboin HV, Meyer UA, Gonzalez FJ. The CYP2D gene subfamily: analysis of the molecular basis of the desbrisoquine 4-hydroxylase deficiency in DA rats. Biochemistry. 1989;28:7349–55.

    Article  CAS  Google Scholar 

  15. Nelson DR, Zeldin DC, Hoffman SM, Maltais LJ, Wain HM, Nebert DW. Comparison of cytochrome P450 (CYP) genes from the mouse and human genomes, including nomenclature recommendations for genes, pseudogenes and alternative-splice variants. Pharmacogenetics. 2004;14:1–18.

    Article  CAS  Google Scholar 

  16. Yasukochi Y, Satta Y. Evolution of the CYP2D gene cluster in humans and four non-human primates. Genes Genet Syst. 2011;86:109–16.

    Article  CAS  Google Scholar 

  17. Knych HK, Baden RW, Gretler SR, McKemie DS. Characterization of the in vitro CYP450 mediated metabolism of the polymorphic CYP2D6 probe drug codeine in horses. Biochem Pharmacol. 2019;168:184–92.

    Article  CAS  Google Scholar 

  18. Gretler SR, Finno CJ, McKemie DS, Kass PH, Knych HK. Metabolism, pharmacokinetics and selected pharmacodynamic effects of codeine following a single oral administration to horses. Vet Anaesth Analg. 2020;

  19. Poulsen L, Brosen K, Arendt-Nielsen L, Gram LF, Elbaek K, Sindrup SH. Codeine and morphine in extensive and poor metabolizers of sparteine: pharmacokinetics, analgesic effect and side effects. Eur J Clin Pharmacol. 1996;51:289–95.

    Article  CAS  Google Scholar 

  20. Sindrup SH, Arendt-Nielsen L, Brosen K, Bjerring P, Angelo HR, Eriksen B, et al. The effect of quinidine on the analgesic effect of codeine. Eur J Clin Pharmacol. 1992;42:587–91.

    Article  CAS  Google Scholar 

  21. Kirchheiner J, Schmidt H, Tzvetkov M, Keulen JT, Lotsch J, Roots I, et al. Pharmacokinetics of codeine and its metabolite morphine in ultra-rapid metabolizers due to CYP2D6 duplication. Pharmacogenomics J. 2007;7:257–65.

    Article  CAS  Google Scholar 

  22. Yue QY, Svensson JO, Alm C, Sjoqvist F, Sawe J. Codeine O-demethylation co-segregates with polymorphic debrisoquine hydroxylation. Br J Clin Pharmacol. 1989;28:639–45.

    Article  CAS  Google Scholar 

  23. Jackson PR, Tucker GT, Woods HF. Testing for bimodality in frequency distributions of data suggesting polymorphisms of drug metabolism - histograms and probit plots. Br J Clin Pharmacol. 1989;28:647–53.

    Article  CAS  Google Scholar 

  24. García-Quetglas E, Azanza JR, Sádaba B, Munoz MJ, Gil I, Campanero MA. Pharmacokinetics of tramadol enantiomers and their respective phase I metabolites in relation to CYP2D6 phenotype. Pharmacol Res. 2007;55:122–30.

    Article  Google Scholar 

  25. Yue QY, Svensson JO, Sjoqvist F, Säwe J. A comparison of the pharmacokinetics of codeine and its metabolites in healthy Chinese and Caucasian extensive hydroxylators of desbrisoquine. Br J Clin Pharmacol. 1991;31:635–42.

    Article  CAS  Google Scholar 

  26. Sachse C, Brockmoller J, Bauer S, Roots I. Cytochrome P450 2D6 variants in a Caucasian population: allele frequencies and phenotypic consequences. Am J Hum Genet. 1997;60:284–95.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Crews KR, Gaedigk A, Dunnenberger HM, Leeder JS, Klein TE, Caudle KE, et al. Clinical pharmacogenetics consortium guidelines for cytochrome P450 2D6 genotype and codeine therapy: 2014 update. Clin Pharmacol Ther. 2014;95:376–82.

    Article  CAS  Google Scholar 

  28. Haufroid V, Hantson P. CYP2D6 genetic polymorphisms and their relevance for poisoning due to amfetamines, opioid analgesics, and antidepressants. Clin Toxicol. 2015;53:501–10.

    Article  CAS  Google Scholar 

  29. Zhou SF. Cytochrome P450 2D6 structure, function, regulation, and polymorphism. 1st ed. Boca Raton: CRC Press; 2016. p. 415–31.

    Book  Google Scholar 

  30. Lampe JW, Bigler J, Bush AC, Potter JD. Prevalence of polymorphisms in the human UDP-glucuronosyltransferase 2B family: UGT2B4(D458E), UGT2B7(H268Y), and UGT2B15(D85Y). Cancer Epidemiol Biomark Prev. 2000;9:329–33.

    CAS  Google Scholar 

  31. Coffman BL, King CD, Rios GR, Tephly TR. The glucuoronidation of opioids, other xenobiotics, and androgens by human UGT2B7Y(268) and UGT2B7H(268). Drug Metab Dispos. 1998;26:73–7.

    CAS  PubMed  Google Scholar 

  32. Thorn CF, Klein TE, Altman RB. Codeine and morphine pathway. Pharmacogenet Genomics. 2009;19:556–8.

    Article  CAS  Google Scholar 

  33. He YJ, Brockmoller J, Schmidt H, Roots I, Kirchheiner J. CYP2D6 ultrarapid metabolism and morphine/codeine ratios in blood: was it codeine or heroin? J Anal Toxicol. 2008;32:178–82.

    Article  CAS  Google Scholar 

  34. Madadi P, Koren G, Cairns J, Chitayat D, Gaedigk A, Leeder JS, et al. Safety of codeine during breastfeeding: fatal morphine poisoning in the breastfed neonate of a mother prescribed codeine. Can Fam Physician. 2007;53:33–5.

    PubMed  PubMed Central  Google Scholar 

  35. Nakamura GR, Griesemer EC, Noguchi TT. Antemortem conversion of codeine to morphine in man. J Forensic Sci. 1976;21:518–24.

    Article  CAS  Google Scholar 

  36. Gaedigk A, Simon S, Pearce R, Bradford LD, Kennedy MJ, Leeder JS. The CYP2D6 activity score: translating genotype information into a qualitative measure of phenotype. Clin Pharmacol Ther. 2008;83:234–42.

    Article  CAS  Google Scholar 

  37. Frank D, Jaehde U, Fuhr U. Evaluation of probe drugs and pharmacokinetic metrics for CYP2D6 phenotyping. Eur J Clin Pharmacol. 2007;63:321–33.

    Article  CAS  Google Scholar 

  38. Afsar NA, Bruckmueller H, Werk AN, Nisar MK, Ahmad HR, Cascorbi I. Implications of genetic variation of common drug metabolizing enzymes and ABC transporters among the Pakistani population. Sci Rep. 2019;9:1–12.

    Article  CAS  Google Scholar 

  39. O’Mahoney MS, Woodhouse KW. Age, environmental factors and drug metabolism. Pharmacol Ther. 1994;61:279–87.

    Article  Google Scholar 

Download references


The authors would like to thank Kelsey Seminoff, Kirsten Kanarr, Stacy Steinmetz, Russell Baden, Briana Hamamoto-Hardman and Sandy Yim for technical assistance.

Competing interests and declarations

Philip Kass is a member of the BMC Veterinary Research Editorial Board (Associate Editor). No other authors have any competing interests or declarations.


Financial support for this project was provided by the Center for Equine Health with funds provided by the State of California pari-mutuel fund. Additional funding was provided by the California Horse Racing board. Support for CJF was provided by the National.

Institutes of Health (L40 TR001136). The funding organizations had no role in study design, execution, analysis or interpretation of the data.

Author information

Authors and Affiliations



SRG- study execution, data analysis, manuscript preparation, reviewed and approved final manuscript CJF-study design, reviewed and approved final manuscript; PHK-data analysis, manuscript preparation, reviewed and approved final manuscript; HKK-study design, study execution, data analysis, manuscript preparation, reviewed and approved final manuscript.

Corresponding author

Correspondence to H. K. Knych.

Ethics declarations

Ethics approval and consent to participate

This study has been reviewed and approved by the Institutional Animal Care and Use Committee of the University of California, Davis. Informed consent was provided by owners or the owner’s representative for all privately owned animals utilized in this study.

Consent for publication

Not applicable.

Competing interests

The authors have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gretler, S.R., Finno, C.J., Kass, P.H. et al. Functional phenotyping of the CYP2D6 probe drug codeine in the horse. BMC Vet Res 17, 77 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: