Clinical pharmacogenomic testing and reporting: A technical standard of the American College of Medical Genetics and Genomics (ACMG).

Pharmacogenomic testing interrogates germline sequence variants implicated in interindividual drug response variability to infer a drug response phenotype and to guide medication management for certain drugs. Specifically, discrete aspects of pharmacokinetics, such as drug metabolism, and pharmacodynamics, as well as drug sensitivity, can be predicted by genes that code for proteins involved in these pathways. Pharmacogenomics is unique and differs from inherited disease genetics because the drug response phenotype can be drug-dependent and is often unrecognized until an unexpected drug reaction occurs or a patient fails to respond to a medication. Genes and variants with sufficiently high levels of evidence and consensus may be included in a clinical pharmacogenomic test; however, result interpretation and phenotype prediction can be challenging for some genes and medications. This document provides a resource for laboratories to develop and implement clinical pharmacogenomic testing by summarizing publicly available resources and detailing best practices for pharmacogenomic nomenclature, testing, result interpretation, and reporting.

Ontogeny of Scaling Factors for Pediatric Physiology-Based Pharmacokinetic Modeling and Simulation: Microsomal Protein Per Gram of Liver.

Microsomal protein per gram of liver (MPPGL) is an important scaling factor for bottom-up physiology-based pharmacokinetic modeling and simulation, but data in pediatrics are limited. Therefore, MPPGL was determined in 160 liver samples from pediatric (n = 129) and adult (n = 31) donors obtained from four sources: the University of Maryland Brain and Tissue Bank (UMBTB), tissue retrieval services at the University of Minnesota and University of Pittsburgh, and Sekisui-Xenotech. Tissues were homogenized and subjected to differential centrifugation to prepare microsomes, and cytochrome c reductase activities in tissue homogenates and microsomes were used to estimate cytochrome P450 reductase (POR) activity as a marker of microsomal recovery; microsomal POR content was also assessed by quantitative proteomics. MPPGL values varied 5- to 10-fold within various age groups/developmental stages, and tissue source was identified as a contributing factor. Using a "trimmed" dataset comprised of samples ranging from 3 to 18 years of age common to the four sources, POR protein abundance and activity in microsomes and POR activity in homogenates was lower in UMBTB samples (autopsy) compared with other sources (perfused/flash-frozen). Regression analyses revealed that the UMBTB samples were driving an apparent age effect as no effect of age on log-transformed MPPGL values was observed when the UMBTB samples were excluded. We conclude that a mean±SD MPPGL value of 30.4±1.7 mg/g is representative between one month postnatal age and early adulthood. Potential source effects should be considered for studies involving tissue samples from multiple sources with different procurement and processing procedures. SIGNIFICANCE STATEMENT: Microsomal protein per gram of liver (MPPGL) is an important scaling factor for bottom up PBPK modeling and simulation, but data in pediatrics are limited. Although MPPGL varies 5- to 10-fold at a given developmental stage, a value of 30.4 ± 1.7 mg/g (mean ± SD) is representative between one month postnatal age and early adulthood. However, when tissue samples are obtained from multiple sources, different procurement and processing procedures may influence the results and should be taken into consideration.

A cross-sectional study of the relationship between CYP2D6 and CYP2C19 variations and depression symptoms, for women taking SSRIs during pregnancy.

Depression during pregnancy affects 10-15% of women, and 5% of women take antidepressants during pregnancy. Clinical guidelines provide recommendations for selective serotonin reuptake inhibitor (SSRI) drug choice and dose based on CYP2D6 and CYP2C19 genotype; however, they are based on evidence from non-pregnant cohorts. This study aimed to test the hypothesis that women with function-altering variants (increased, decreased, or no function) in these pharmacogenes, taking SSRIs prenatally, would have more depression symptoms than women whose pharmacogenetic variants are associated with normal SSRI metabolism. Comprehensive CYP2D6 and CYP2C19 genotyping using a range of methods, including gene copy number analysis, was performed as secondary analyses on two longitudinal cohorts of pregnant women (N = 83) taking the SSRIs paroxetine, citalopram, escitalopram, or sertraline. The Kruskal-Wallis test compared mean depression scores across four predicted metabolizer groups: poor (n = 5), intermediate (n = 10), normal (n = 53), and ultrarapid (n = 15). There were no significant differences between mean depression scores across the four metabolizer groups (H(3) = .73, p = .87, eta-squared = .029, epsilon-squared = .0089). This is the first study of the relationship in pregnancy between CYP2C19 pharmacogenetic variations and depression symptoms in the context of SSRI use. Findings from this initial study do not support the clinical use of pharmacogenetic testing for SSRI use during the second or third trimesters of pregnancy, but these findings should be confirmed in larger cohorts. There is an urgent need for further research to clarify the utility of pharmacogenetic testing for pregnant women, especially as companies offering direct-to-consumer genetic testing expand their marketing efforts.

CYP2D6 genotype and reduced codeine analgesic effect in real-world clinical practice.

Cytochrome P450 2D6 (CYP2D6) O-demethylates codeine to the active drug, morphine. However, the utility of testing for CYP2D6 metabolizer status in patients receiving codeine in real-world clinical practice is poorly defined. Using data from a DNA bank linked to de-identified electronic health records, we studied 157 patients with a baseline pain score higher than 4 (0-10 scale) who received codeine. Based on CYP2D6 genotyping, 69 were classified as poor/intermediate and 88 as normal/ultrarapid CYP2D6 metabolizers. Pain response was defined as a score of 4 or lower while receiving codeine. In a propensity-score adjusted model, poor/intermediate metabolizers had lower odds (OR = 0.35, p = 0.02) of achieving a pain response than normal/ultrarapid metabolizers. To discriminate between codeine responders and nonresponders, a score including CYP2D6 phenotype and clinical variables was built. The response rate was 38.5% among patients in the high, 17.3% in the intermediate, and 9.4% in the low-score groups, respectively (p = 0.001).

Transfer learning enables prediction of CYP2D6 haplotype function.

Cytochrome P450 2D6 (CYP2D6) is a highly polymorphic gene whose protein product metabolizes more than 20% of clinically used drugs. Genetic variations in CYP2D6 are responsible for interindividual heterogeneity in drug response that can lead to drug toxicity and ineffective treatment, making CYP2D6 one of the most important pharmacogenes. Prediction of CYP2D6 phenotype relies on curation of literature-derived functional studies to assign a functional status to CYP2D6 haplotypes. As the number of large-scale sequencing efforts grows, new haplotypes continue to be discovered, and assignment of function is challenging to maintain. To address this challenge, we have trained a convolutional neural network to predict functional status of CYP2D6 haplotypes, called Hubble.2D6. Hubble.2D6 predicts haplotype function from sequence data and was trained using two pre-training steps with a combination of real and simulated data. We find that Hubble.2D6 predicts CYP2D6 haplotype functional status with 88% accuracy in a held-out test set and explains 47.5% of the variance in in vitro functional data among star alleles with unknown function. Hubble.2D6 may be a useful tool for assigning function to haplotypes with uncurated function, and used for screening individuals who are at risk of being poor metabolizers.

Review and Consensus on Pharmacogenomic Testing in Psychiatry.

The implementation of pharmacogenomic (PGx) testing in psychiatry remains modest, in part due to divergent perceptions of the quality and completeness of the evidence base and diverse perspectives on the clinical utility of PGx testing among psychiatrists and other healthcare providers. Recognizing the current lack of consensus within the field, the International Society of Psychiatric Genetics assembled a group of experts to conduct a narrative synthesis of the PGx literature, prescribing guidelines, and product labels related to psychotropic medications as well as the key considerations and limitations related to the use of PGx testing in psychiatry. The group concluded that to inform medication selection and dosing of several commonly-used antidepressant and antipsychotic medications, current published evidence, prescribing guidelines, and product labels support the use of PGx testing for 2 cytochrome P450 genes (CYP2D6, CYP2C19). In addition, the evidence supports testing for human leukocyte antigen genes when using the mood stabilizers carbamazepine (HLA-A and HLA-B), oxcarbazepine (HLA-B), and phenytoin (CYP2C9, HLA-B). For valproate, screening for variants in certain genes (POLG, OTC, CSP1) is recommended when a mitochondrial disorder or a urea cycle disorder is suspected. Although barriers to implementing PGx testing remain to be fully resolved, the current trajectory of discovery and innovation in the field suggests these barriers will be overcome and testing will become an important tool in psychiatry.

The Respective Roles of CYP3A4 and CYP2D6 in the Metabolism of Pimozide to Established and Novel Metabolites.

Pimozide is a dopamine receptor antagonist indicated for the treatment of Tourette syndrome. Prior in vitro studies characterized N-dealkylation of pimozide to 1,3-dihydro-1-(4-piperidinyl)-2H-benzimidazol-2-one (DHPBI) via CYP3A4 and, to a lesser extent, CYP1A2 as the only notable routes of pimozide biotransformation. However, drug-drug interactions between pimozide and CYP2D6 inhibitors and CYP2D6 genotype-dependent effects have since been observed. To reconcile these incongruities between the prior in vitro and in vivo studies, we characterized two novel pimozide metabolites: 5-hydroxypimozide and 6-hydroxypimozide. Notably, 5-hydroxypimozide was the major metabolite produced by recombinant CYP2D6 (Km ∼82 nM, V max ∼0.78 pmol/min per picomoles), and DHPBI was the major metabolite produced by recombinant CYP3A4 (apparent Km ∼1300 nM, V max ∼2.6 pmol/min per picomoles). Kinetics in pooled human liver microsomes (HLMs) for the 5-hydroxylation (Km ∼2200 nM, V max ∼59 pmol/min per milligram) and N-dealkylation (Km ∼3900 nM, V max ∼600 pmol/min per milligram) reactions were also determined. Collectively, formation of DHPBI, 5-hydroxypimozide, and 6-hydroxypimozide accounted for 90% of pimozide depleted in incubations of NADPH-supplemented pooled HLMs. Studies conducted in HLMs isolated from individual donors with specific cytochrome P450 isoform protein abundances determined via mass spectrometry revealed that 5-hydroxypimozide (r 2 = 0.94) and 6-hydroxypimozide (r 2 = 0.86) formation rates were correlated with CYP2D6 abundance, whereas the DHPBI formation rate (r 2 = 0.98) was correlated with CYP3A4 abundance. Furthermore, the HLMs differed with respect to their capacity to form 5-hydroxypimozide relative to DHPBI. Collectively, these data confirm a role for CYP2D6 in pimozide clearance via 5-hydroxylation and provide an explanation for a lack of involvement when only DHPBI formation was monitored in prior in vitro studies. SIGNIFICANCE STATEMENT: Current CYP2D6 genotype-guided dosing information in the pimozide label is discordant with available knowledge regarding the primary biotransformation pathways. Herein, we characterize the CYP2D6-dependent biotransformation of pimozide to previously unidentified metabolites. In human liver microsomes, formation rates for the novel metabolites and a previously identified metabolite were determined to be a function of CYP2D6 and CYP3A4 content, respectively. These findings provide a mechanistic basis for observations of CYP2D6 genotype-dependent pimozide clearance in vivo.

Structural variation at the CYP2C locus: Characterization of deletion and duplication alleles.

The human CYP2C locus harbors the polymorphic CYP2C18, CYP2C19, CYP2C9, and CYP2C8 genes, and of these, CYP2C19 and CYP2C9 are directly involved in the metabolism of ~15% of all medications. All variant CYP2C19 and CYP2C9 star (*) allele haplotypes currently cataloged by the Pharmacogene Variation (PharmVar) Consortium are defined by sequence variants. To determine if structural variation also occurs at the CYP2C locus, the 10q23.33 region was interrogated across deidentified clinical chromosomal microarray (CMA) data from 20,642 patients tested at two academic medical centers. Fourteen copy number variants that affected the coding region of CYP2C genes were detected in the clinical CMA cohorts, which ranged in size from 39.2 to 1,043.3 kb. Selected deletions and duplications were confirmed by MLPA or ddPCR. Analysis of the clinical CMA and an additional 78,839 cases from the Database of Genomic Variants (DGV) and ClinGen (total n = 99,481) indicated that the carrier frequency of a CYP2C structural variant is ~1 in 1,000, with ~1 in 2,000 being a CYP2C19 full gene or partial-gene deletion carrier, designated by PharmVar as CYP2C19*36 and *37, respectively. Although these structural variants are rare in the general population, their detection will likely improve metabolizer phenotype prediction when interrogated for research and/or clinical testing.

Stargazer: a software tool for calling star alleles from next-generation sequencing data using CYP2D6 as a model.

PURPOSE: Genotyping CYP2D6 is important for precision drug therapy because the enzyme it encodes metabolizes approximately 25% of drugs, and its activity varies considerably among individuals. Genotype analysis of CYP2D6 is challenging due to its highly polymorphic nature. Over 100 haplotypes (star alleles) have been defined for CYP2D6, some involving a gene conversion with its nearby nonfunctional but highly homologous paralog CYP2D7. We present Stargazer, a new bioinformatics tool that uses next-generation sequencing (NGS) data to call star alleles for CYP2D6 ( https://stargazer.gs.washington.edu/stargazerweb/ ). Stargazer is currently being extended for other pharmacogenes. METHODS: Stargazer identifies star alleles from NGS data by detecting single nucleotide variants, insertion-deletion variants, and structural variants. Stargazer detects structural variation, including gene deletions, duplications, and conversions, by calculating paralog-specific copy numbers from read depths. RESULTS: We applied Stargazer to the NGS data of 32 ethnically diverse HapMap trios that were genotyped by TaqMan assays, long-range polymerase chain reaction, quantitative multiplex polymerase chain reaction, high-resolution melting analysis, and/or Sanger sequencing. CYP2D6 genotyping by Stargazer was 99.0% concordant with the data obtained by these methods, and showed that 28.1% of the samples had structural variation including CYP2D6/CYP2D7 hybrids. CONCLUSION: Accurate genotyping of pharmacogenes with NGS and subsequent allele calling with Stargazer will aid the implementation of precision drug therapy.

Ontogeny of Hepatic Sulfotransferases and Prediction of Age-Dependent Fractional Contribution of Sulfation in Acetaminophen Metabolism.

Cytosolic sulfotransferases (SULTs), including SULT1A, SULT1B, SULT1E, and SULT2A isoforms, play noteworthy roles in xenobiotic and endobiotic metabolism. We quantified the protein abundances of SULT1A1, SULT1A3, SULT1B1, and SULT2A1 in human liver cytosol samples (n = 194) by liquid chromatography-tandem mass spectrometry proteomics. The data were analyzed for their associations by age, sex, genotype, and ethnicity of the donors. SULT1A1, SULT1B1, and SULT2A1 showed significant age-dependent protein abundance, whereas SULT1A3 was invariable across 0-70 years. The respective mean abundances of SULT1A1, SULT1B1, and SULT2A1 in neonatal samples was 24%, 19%, and 38% of the adult levels. Interestingly, unlike UDP-glucuronosyltransferases and cytochrome P450 enzymes, SULT1A1 and SULT2A1 showed the highest abundance during early childhood (1 to <6 years), which gradually decreased by approx. 40% in adolescents and adults. SULT1A3 and SULT1B1 abundances were significantly lower in African Americans compared with Caucasians. Multiple linear regression analysis further confirmed the association of SULT abundances by age, ethnicity, and genotype. To demonstrate clinical application of the characteristic SULT ontogeny profiles, we developed and validated a proteomics-informed physiologically based pharmacokinetic model of acetaminophen. The latter confirmed the higher fractional contribution of sulfation over glucuronidation in the metabolism of acetaminophen in children. The study thus highlights that the ontogeny-based age-dependent fractional contribution (fm) of individual drug-metabolizing enzymes has better potential in prediction of drug-drug interactions and the effect of genetic polymorphisms in the pediatric population.

Hepatic Abundance and Activity of Androgen- and Drug-Metabolizing Enzyme UGT2B17 Are Associated with Genotype, Age, and Sex.

The major objective of this study was to investigate the association of genetic and nongenetic factors with variability in protein abundance and in vitro activity of the androgen-metabolizing enzyme UGT2B17 in human liver microsomes (n = 455). UGT2B17 abundance was quantified by liquid chromatography-tandem mass spectrometry proteomics, and enzyme activity was determined by using testosterone and dihydrotestosterone as in vitro probe substrates. Genotyping or gene resequencing and mRNA expression were also evaluated. Multivariate analysis was used to test the association of UGT2B17 copy number variation, single nucleotide polymorphisms (SNPs), age, and sex with its mRNA expression, abundance, and activity. UGT2B17 gene copy number and SNPs (rs7436962, rs9996186, rs28374627, and rs4860305) were associated with gene expression, protein levels, and androgen glucuronidation rates in a gene dose-dependent manner. UGT2B17 protein (mean ± S.D. picomoles per milligram of microsomal protein) is sparsely expressed in children younger than 9 years (0.12 ± 0.24 years) but profoundly increases from age 9 years to adults (∼10-fold) with ∼2.6-fold greater abundance in males than in females (1.2 vs. 0.47). Association of androgen glucuronidation with UGT2B15 abundance was observed only in the low UGT2B17 expressers. These data can be used to predict variability in the metabolism of UGT2B17 substrates. Drug companies should include UGT2B17 in early phenotyping assays during drug discovery to avoid late clinical failures.

Genetic and Nongenetic Factors Associated with Protein Abundance of Flavin-Containing Monooxygenase 3 in Human Liver.

Hepatic flavin-containing mono-oxygenase 3 (FMO3) metabolizes a broad array of nucleophilic heteroatom (e.g., N or S)-containing xenobiotics (e.g., amphetamine, sulindac, benzydamine, ranitidine, tamoxifen, nicotine, and ethionamide), as well as endogenous compounds (e.g., catecholamine and trimethylamine). To predict the effect of genetic and nongenetic factors on the hepatic metabolism of FMO3 substrates, we quantified FMO3 protein abundance in human liver microsomes (HLMs; n = 445) by liquid chromatography-tandem mass chromatography proteomics. Genotyping/gene resequencing, mRNA expression, and functional activity (with benzydamine as probe substrate) of FMO3 were also evaluated. FMO3 abundance increased 2.2-fold (13.0 ± 11.4 pmol/mg protein vs. 28.0 ± 11.8 pmol/mg protein) from neonates to adults. After 6 years of age, no significant difference in FMO3 abundance was found between children and adults. Female donors exhibited modestly higher mRNA fragments per kilobase per million reads values (139.9 ± 76.9 vs. 105.1 ± 73.1; P < 0.001) and protein FMO3 abundance (26.7 ± 12.0 pmol/mg protein vs. 24.1 ± 12.1 pmol/mg protein; P < 0.05) compared with males. Six single nucleotide polymorphisms (SNPs), including rs2064074, rs28363536, rs2266782 (E158K), rs909530 (N285N), rs2266780 (E308G), and rs909531, were associated with significantly decreased protein abundance. FMO3 abundance in individuals homozygous and heterozygous for haplotype 3 (H3), representing variant alleles for all these SNPs (except rs2066534), were 50.8% (P < 0.001) and 79.5% (P < 0.01), respectively, of those with the reference homozygous haplotype (H1, representing wild-type). In summary, FMO3 protein abundance is significantly associated with age, gender, and genotype. These data are important in predicting FMO3-mediated heteroatom-oxidation of xenobiotics and endogenous biomolecules in the human liver.

Prediction of CYP2D6 phenotype from genotype across world populations.

PURPOSE: Owing to its highly polymorphic nature and major contribution to the metabolism and bioactivation of numerous clinically used drugs, CYP2D6 is one of the most extensively studied drug-metabolizing enzymes and pharmacogenes. CYP2D6 alleles confer no, decreased, normal, or increased activity and cause a wide range of activity among individuals and between populations. However, there is no standard approach to translate diplotypes into predicted phenotype. METHODS: We exploited CYP2D6 allele-frequency data that have been compiled for Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines (>60,000 subjects, 173 reports) in order to estimate genotype-predicted phenotype status across major world populations based on activity score (AS) assignments. RESULTS: Allele frequencies vary considerably across the major ethnic groups predicting poor metabolizer status (AS = 0) between 0.4 and 5.4% across world populations. The prevalence of genotypic intermediate (AS = 0.5) and normal (AS = 1, 1.5, or 2) metabolizers ranges between 0.4 and 11% and between 67 and 90%, respectively. Finally, 1 to 21% of subjects (AS >2) are predicted to have ultrarapid metabolizer status. CONCLUSIONS: This comprehensive study summarizes allele frequencies, diplotypes, and predicted phenotype across major populations, providing a rich data resource for clinicians and researchers. Challenges of phenotype prediction from genotype data are highlighted and discussed.Genet Med 19 1, 69-76.

Age-dependent Protein Abundance of Cytosolic Alcohol and Aldehyde Dehydrogenases in Human Liver.

Hepatic cytosolic alcohol and aldehyde dehydrogenases (ADHs and ALDHs) catalyze the biotransformation of xenobiotics (e.g., cyclophosphamide and ethanol) and vitamin A. Because age-dependent hepatic abundance of these proteins is unknown, we quantified protein expression of ADHs and ALDH1A1 in a large cohort of pediatric and adult human livers by liquid chromatography coupled with tandem mass spectrometry proteomics. Purified proteins were used as calibrators. Two to three surrogate peptides per protein were quantified in trypsin digests of liver cytosolic samples and calibrator proteins under optimal conditions of reproducibility. Neonatal levels of ADH1A, ADH1B, ADH1C, and ALDH1A1 were 3-, 8-, 146-, and 3-fold lower than the adult levels, respectively. For all proteins, the abundance steeply increased during the first year of life, which mostly reached adult levels during early childhood (age between 1 and 6 years). Only for ADH1A protein abundance in adults (age > 18 year) was ∼40% lower relative to the early childhood group. Abundances of ADHs and ALDH1A1 were not associated with sex in samples with age > 1 year compared with males. Known single nucleotide polymorphisms had no effect on the protein levels of these proteins. Quantification of ADHs and ALDH1A1 protein levels could be useful in predicting disposition and response of substrates of these enzymes in younger children.

Age-Dependent Absolute Abundance of Hepatic Carboxylesterases (CES1 and CES2) by LC-MS/MS Proteomics: Application to PBPK Modeling of Oseltamivir In Vivo Pharmacokinetics in Infants.

The age-dependent absolute protein abundance of carboxylesterase (CES) 1 and CES2 in human liver was investigated and applied to predict infant pharmacokinetics (PK) of oseltamivir. The CES absolute protein abundance was determined by liquid chromatography-tandem mass spectrometry proteomics in human liver microsomal and cytosolic fractions prepared from tissue samples obtained from 136 pediatric donors and 35 adult donors. Two surrogate peptides per protein were selected for the quantification of CES1 and CES2 protein abundance. Purified CES1 and CES2 protein standards were used as calibrators, and the heavy labeled peptides were used as the internal standards. In hepatic microsomes, CES1 and CES2 abundance (in picomoles per milligram total protein) increased approximately 5-fold (315.2 vs. 1664.4) and approximately 3-fold (59.8 vs. 174.1) from neonates to adults, respectively. CES1 protein abundance in liver cytosol also showed age-dependent maturation. Oseltamivir carboxylase activity was correlated with protein abundance in pediatric and adult liver microsomes. The protein abundance data were then used to model in vivo PK of oseltamivir in infants using pediatric physiologically based PK modeling and incorporating the protein abundance-based ontogeny function into the existing pediatric Simcyp model. The predicted pediatric area under the curve, maximal plasma concentration, and time for maximal plasma concentration values were below 2.1-fold of the clinically observed values, respectively.

Long-Read Single Molecule Real-Time Full Gene Sequencing of Cytochrome P450-2D6.

The cytochrome P450-2D6 (CYP2D6) enzyme metabolizes ∼25% of common medications, yet homologous pseudogenes and copy number variants (CNVs) make interrogating the polymorphic CYP2D6 gene with short-read sequencing challenging. Therefore, we developed a novel long-read, full gene CYP2D6 single molecule real-time (SMRT) sequencing method using the Pacific Biosciences platform. Long-range PCR and CYP2D6 SMRT sequencing of 10 previously genotyped controls identified expected star (*) alleles, but also enabled suballele resolution, diplotype refinement, and discovery of novel alleles. Coupled with an optimized variant-calling pipeline, CYP2D6 SMRT sequencing was highly reproducible as triplicate intra- and inter-run nonreference genotype results were completely concordant. Importantly, targeted SMRT sequencing of upstream and downstream CYP2D6 gene copies characterized the duplicated allele in 15 control samples with CYP2D6 CNVs. The utility of CYP2D6 SMRT sequencing was further underscored by identifying the diplotypes of 14 samples with discordant or unclear CYP2D6 configurations from previous targeted genotyping, which again included suballele resolution, duplicated allele characterization, and discovery of a novel allele and tandem arrangement. Taken together, long-read CYP2D6 SMRT sequencing is an innovative, reproducible, and validated method for full-gene characterization, duplication allele-specific analysis, and novel allele discovery, which will likely improve CYP2D6 metabolizer phenotype prediction for both research and clinical testing applications.

Common CYP2D6 polymorphisms affecting alternative splicing and transcription: long-range haplotypes with two regulatory variants modulate CYP2D6 activity.

Cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of 25% of clinically used drugs. Genetic polymorphisms cause substantial variation in CYP2D6 activity and serve as biomarkers guiding drug therapy. However, genotype-phenotype relationships remain ambiguous except for poor metabolizers carrying null alleles, suggesting the presence of yet unknown genetic variants. Searching for regulatory CYP2D6 polymorphisms, we find that a SNP defining the CYP2D6*2 allele, rs16947 [R296C, 17-60% minor allele frequency (MAF)], previously thought to convey normal activity, alters exon 6 splicing, thereby reducing CYP2D6 expression at least 2-fold. In addition, two completely linked SNPs (rs5758550/rs133333, MAF 13-42%) increase CYP2D6 transcription more than 2-fold, located in a distant downstream enhancer region (>100 kb) that interacts with the CYP2D6 promoter. In high linkage disequilibrium (LD) with each other, rs16947 and the enhancer SNPs form haplotypes that affect CYP2D6 enzyme activity in vivo. In a pediatric cohort of 164 individuals, rs16947 alone (minor haplotype frequency 28%) was associated with reduced CYP2D6 metabolic activity (measured as dextromethorphan/metabolite ratios), whereas rs5758550/rs133333 alone (frequency 3%) resulted in increased CYP2D6 activity, while haplotypes containing both rs16947 and rs5758550/rs133333 were similar to the wild-type. Other alleles used in biomarker panels carrying these variants such as CYP2D6*41 require re-evaluation of independent effects on CYP2D6 activity. The occurrence of two regulatory variants of high frequency and in high LD, residing on a long haplotype, highlights the importance of gene architecture, likely shaped by evolutionary selection pressures, in determining activity of encoded proteins.

Characterization of Atomoxetine Biotransformation and Implications for Development of PBPK Models for Dose Individualization in Children.

Atomoxetine (ATX) is a second-line nonstimulant medication used to control symptoms of attention deficit hyperactivity disorder (ADHD). Inconsistent therapeutic efficacy has been reported with ATX, which may be related to variable CYP2D6-mediated drug clearance. We characterized ATX metabolism in a panel of human liver samples as a basis for a bottom-up PBPK model to aid in ATX exposure prediction and control. Km, Vmax, and Clint values in pooled human liver microsomes (HLMs) were 2.4 µM, 479 pmol/min/mg protein, and 202 µl/min/mg protein, respectively. Mean population values of kinetic parameters are not adequate to describe variability in a population, given that Km, Vmax, and Clint values from single-donor HLMs ranged from 0.93 to 79.2 µM, 20.0 to 1600 pmol/min/mg protein, and 0.3 to 936 µl/min/mg protein. All kinetic parameters were calculated from 4-hydroxyatomoxetine (4-OH-ATX) formation. CYP2E1 and CYP3A contributed to 4-OH-ATX formation in livers with CYP2D6 intermediate and poor metabolizer status. In HLMs with lower CYP2D6 activity levels, 2-hydroxymethylatomoxetine (2-CH2OH-ATX) formation became a more predominant pathway of metabolism, which appeared to be catalyzed by CYP2B6. ATX biotransformation at clinically relevant plasma concentrations was characterized in a panel of pediatric HLM (n = 116) samples by evaluating primary metabolites. Competing pathways of ATX metabolism [N-desmethylatomoxetine (NDM-ATX) and 2-CH2OH-ATX formation] had increasing importance in livers with lesser CYP2D6 activity, but, overall ATX clearance was still compromised. Modeling ATX exposure to individualize therapy would require comprehensive knowledge of factors that affect CYP2D6 activity as well as an understanding of competing pathways, particularly for individuals with lower CYP2D6 activity.

Developmental Expression of CYP2B6: A Comprehensive Analysis of mRNA Expression, Protein Content and Bupropion Hydroxylase Activity and the Impact of Genetic Variation.

Although CYP2B6 catalyzes the biotransformation of many drugs used clinically for children and adults, information regarding the effects of development on CYP2B6 expression and activity are scarce. Utilizing a large panel of human liver samples (201 donors: 24 fetal, 141 pediatric, and 36 adult), we quantified CYP2B6 mRNA and protein expression levels, characterized CYP2B6 (bupropion hydroxylase) activity in human liver microsomes (HLMs), and performed an extensive genotype analysis to differentiate CYP2B6 haplotypes such that the impact of genetic variation on these parameters could be assessed. Fetal livers contained extremely low levels of CYP2B6 mRNA relative to postnatal samples and fetal HLMs did not appear to catalyze bupropion hydroxylation; however, fetal CYP2B6 protein levels were not significantly different from postnatal levels. Considerable interindividual variation in CYP2B6 mRNA expression, protein levels, and activity was observed in postnatal HLMs (mRNA, ∼40,000-fold; protein, ∼300-fold; activity, ∼600-fold). The extremely wide range of interindividual variability in CYP2B6 expression and activity was significantly associated with age (P < 0.01) following log transformation of the data. Our data suggest that CYP2B6 activity appears as early as the first day of life, increases through infancy, and by 1 year of age, CYP2B6 levels and activity may approach those of adults. Surprisingly, CYP2B6 interindividual variability was not significantly associated with genetic variation in CYP2B6, nor was it associated with differences in gender or ethnicity, suggesting that factors other than these are largely responsible for the wide range of variability in CYP2B6 expression and activity observed among a large group of individuals/samples.

Population pharmacokinetics of oral baclofen in pediatric patients with cerebral palsy.

OBJECTIVE: To characterize the population pharmacokinetics (PK) of oral baclofen and assess impact of patient-specific covariates in children with cerebral palsy (CP) in order to support its clinical use. SUBJECTS DESIGN: Children (2-17 years of age) with CP received a dose of titrated oral baclofen from 2.5 mg 3 times a day to a maximum tolerated dose of up to 20 mg 4 times a day. PK sampling followed titration of 10-12 weeks. Serial R- and S-baclofen plasma concentrations were measured for up to 16 hours in 49 subjects. Population PK modeling was performed using NONMEM 7.1 (ICON PLC; Ellicott City, Maryland). RESULTS: R- and S-baclofen showed identical concentration-time profiles. Both baclofen enantiomers exhibited linear and dose/kg-proportional PK, and no sex differences were observed. Average baclofen terminal half-life was 4.5 hours. A 2-compartment PK model with linear elimination and transit absorption steps adequately described concentration-time profiles of both baclofen enantiomers. The mean population estimate of apparent clearance/F was 0.273 L/h/kg with 33.4% inter-individual variability (IIV), and the apparent volume of distribution (Vss/F) was 1.16 L/kg with 43.9% IIV. Delayed absorption was expressed by a mean transit time of 0.389 hours with 83.7% IIV. Body weight, a possible genetic factor, and age were determinants of apparent clearance in these children. CONCLUSION: The PK of oral baclofen exhibited dose-proportionality and were adequately described by a 2-compartment model. Our population PK findings suggest that baclofen dosage can be based on body weight (2 mg/kg per day) and the current baclofen dose escalation strategy is appropriate in the treatment of children with CP older than 2 years of age.