CYP2D6, SULT1A1 and UGT2B17 copy number variation: quantitative detection by multiplex PCR

PharmVar Tutorial on CYP2D6 Structural Variation Testing and Recommendations on Reporting

The Pharmacogene Variation Consortium (PharmVar) provides nomenclature for the highly polymorphic human CYP2D6 gene locus and a comprehensive summary of structural variation. CYP2D6 contributes to the metabolism of numerous drugs and, thus, genetic variation in its gene impacts drug efficacy and safety. To accurately predict a patient's CYP2D6 phenotype, testing must include structural variants including gene deletions, duplications, hybrid genes, and combinations thereof. This tutorial offers a comprehensive overview of CYP2D6 structural variation, terms, and definitions, a review of methods suitable for their detection and characterization, and practical examples to address the lack of standards to describe CYP2D6 structural variants or any other pharmacogene. This PharmVar tutorial offers practical guidance on how to detect the many, often complex, structural variants, as well as recommends terms and definitions for clinical and research reporting. Uniform reporting is not only essential for electronic health record-keeping but also for accurate translation of a patient's genotype into phenotype which is typically utilized to guide drug therapy.

Copy number variations of cytochrome P450 genes in Kinh Vietnamese

Background

The cytochrome P450 (CYP450) family is well known as a major group of drug metabolizing enzymes. The polymorphism of CYP450 genes is the main factor having an impact on the interindividual difference in drug response, including drug efficacy and drug safety. The single nucleotide polymorphism (SNPs) of Vietnamese Kinh has been widely studied, but information about the copy number variations (CNVs) of other CYP450 genes is still unknown.

Objective

To identify the CNV variability of CYP450 in 154 healthy unrelated Kinh Vietnamese, except eCYP2D6, which was previously reported.

Methods

Multiplex Ligation-Dependent Probe Amplification (MLPA) was applied for determination of copy number of 10 CYP450 genes. Later, PCR or quantitative PCR (qPCR) was used to confirm the detected CNVs in randomly chosen subjects.

Results

Of the 154 subjects, along with CYP2D6, 4 other CYP450 genes showed CNVs including duplications (CYP1B1), deletions (CYP2A6 and CYP2C9), and both duplications and deletions (CYP2E1). Among these, CYP2A6 exhibited the greatest frequency of CNVs compared with other CYP450, in which CYP2A6Del accounted for 11%. Meanwhile, allele CYP2E1Del showed the lowest frequency with only 0.3%.

Conclusions

The present study provides new insight into CYP450 CNVs in the Kinh Vietnamese cohort. Our data have contributed to genetic profiling of CYP450 CNVs in Vietnam, which would be helpful for facilitating implementation of pharmacogenetics in drug dosing adjustment in Vietnam.

A longitudinal study of cytochrome P450 2D6 (CYP2D6) activity during adolescence

CYP2D6 substrates are among the most highly prescribed medications in teenagers and also commonly associated with serious adverse events. To investigate the relative contributions of genetic variation, growth, and development on CYP2D6 activity during puberty, healthy children and adolescents 7-15 years of age at enrollment participated in a longitudinal phenotyping study involving administration of 0.3 mg/kg dextromethorphan (DM) and 4-h urine collection every 6 months for 3 years (7 total visits). At each visit, height, weight, and sexual maturity were recorded, and CYP2D6 activity was determined as the urinary molar ratio of DM to its metabolite dextrorphan (DX). A total of 188 participants completed at least one visit, and 102 completed all seven study visits. Following univariate analysis, only CYP2D6 activity score (p < 0.001), urinary pH (p < 0.001), weight (p = 0.018), and attention-deficit/hyperactivity disorder (ADHD) diagnosis (p < 0.001) were significantly correlated with log(DM/DX). Results of linear mixed model analysis with random intercept, random slope covariance structure revealed that CYP2D6 activity score had the strongest effect on log(DM/DX), with model-estimated average log(DM/DX) being 3.8 SDs higher for poor metabolizers than for patients with activity score 3. A moderate effect on log(DM/DX) was observed for sex, and smaller effects were observed for ADHD diagnosis and urinary pH. The log(DM/DX) did not change meaningfully with age or pubertal development. CYP2D6 genotype remains the single, largest determinant of variability in CYP2D6 activity during puberty. Incorporation of genotype-based dosing guidelines should be considered for CYP2D6 substrates given the prevalent use of these agents in this pediatric age group.

Analytical Validation of a Computational Method for Pharmacogenetic Genotyping from Clinical Whole Exome Sequencing

Germline whole exome sequencing from molecular tumor boards has the potential to be repurposed to support clinical pharmacogenomics. However, accurately calling pharmacogenomics-relevant genotypes from exome sequencing data remains challenging. Accordingly, this study assessed the analytical validity of the computational tool, Aldy, in calling pharmacogenomics-relevant genotypes from exome sequencing data for 13 major pharmacogenes. Germline DNA from whole blood was obtained for 164 subjects seen at an institutional molecular solid tumor board. All subjects had whole exome sequencing from Ashion Analytics and panel-based genotyping from an institutional pharmacogenomics laboratory. Aldy version 3.3 was operationalized on the LifeOmic Precision Health Cloud with copy number fixed to two copies per gene. Aldy results were compared with those from genotyping for 56 star allele-defining variants within CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, CYP3A5, CYP4F2, DPYD, G6PD, NUDT15, SLCO1B1, and TPMT. Read depth was >100× for all variants except CYP3A4∗22. For 75 subjects in the validation cohort, all 3393 Aldy variant calls were concordant with genotyping. Aldy calls for 736 diplotypes containing alleles assessed by both platforms were also concordant. Aldy identified additional star alleles not covered by targeted genotyping for 139 diplotypes. Aldy accurately called variants and diplotypes for 13 major pharmacogenes, except for CYP2D6 variants involving copy number variations, thus allowing repurposing of whole exome sequencing to support clinical pharmacogenomics.

The Impact of the CYP2D6 "Enhancer" Single Nucleotide Polymorphism on CYP2D6 Activity

rs5758550 has been associated with enhanced transcription and suggested to be a useful marker of CYP2D6 activity. As there are limited and inconsistent data regarding the utility of this distant "enhancer" single nucleotide polymorphism (SNP), our goal was to further assess the impact of rs5758550 on CYP2D6 activity toward two probe substrates, atomoxetine (ATX) and dextromethorphan (DM), using in vivo urinary metabolite (DM; n = 188) and pharmacokinetic (ATX; n = 70) and in vitro metabolite formation (ATX and DM; n = 166) data. All subjects and tissues were extensively genotyped, the "enhancer" SNP phased with established CYP2D6 haplotypes either computationally or experimentally, and the impact on CYP2D6 activity investigated using several linear models of varying complexity to determine the proportion of variability in CYP2D6 activity captured by each model. For all datasets and models, the "enhancer" SNP had no or only a modest impact on CYP2D6 activity prediction. An increased effect, when present, was more pronounced for ATX than DM suggesting potential substate-dependency. In addition, CYP2D6*2 alleles with the "enhancer" SNP were associated with modestly higher metabolite formation rates in vitro, but not in vivo; no effect was detected for CYP2D6*1 alleles with "enhancer" SNP. In summary, it remains inconclusive whether the small effects detected in this investigation are indeed caused by the "enhancer" SNP or are rather due to its incomplete linkage with other variants within the gene. Taken together, there does not appear to be sufficient evidence to warrant the "enhancer" SNP be included in clinical CYP2D6 pharmacogenetic testing.

CYP2D6 gene resequencing in the Malagasy, a population at the crossroads between Asia and Africa: a pilot study

Background:

Plasmodium vivax malaria is endemic in Madagascar, where populations have genetic inheritance from Southeast Asia and East Africa. Primaquine, a drug of choice for vivax malaria, is metabolized principally via CYP2D6. CYP2D6 variation was characterized by locus-specific gene sequencing and was compared with TaqMan™ genotype data.

Materials & methods:

Long-range PCR amplicons were generated from 96 Malagasy samples and subjected to next-generation sequencing.

Results:

The authors observed high concordance between TaqMan™-based CYP2D6 genotype calls and the base calls from sequencing. In addition, there are new variants and haplotypes present in the Malagasy.

Conclusion:

Sequencing unique admixed populations provides more detailed and accurate insights regarding CYP2D6 variability, which may help optimize primaquine treatment across human genetic diversity.

Simple and Robust Detection of CYP2D6 Gene Deletions and Duplications Using CYP2D8P as Reference

Genotyping of the CYP2D6 gene is the most commonly applied pharmacogenetic test globally. Significant economic interests have led to the development of a plurality of assays, available for almost any genotyping platform or DNA detection chemistry. Of all the genetic variants, copy number variations are particular difficult to detect by polymerase chain reaction. Here, we present two simple novel approaches for the identification of samples carrying either deletions or duplications of the CYP2D6 gene; by relative quantification using a singleplex 5′nuclease real-time PCR assay, and by high-resolution melting of PCR products. These methods make use of universal primers, targeting both the CYP2D6 and the reference gene CYP2D8P, which is necessary for the analysis. The assays were validated against a reference method using a large set of samples. The singleplex nature of the 5′nuclease real-time PCR ensures that the primers anneal with equal affinity to both the sequence of the CYP2D6 and the reference gene. This facilitates robust identification of gene deletions and duplications based on the cycle threshold value. In contrast, the high-resolution melting assay is an end-point PCR, where the identification relies on variations between the amount of product generated from each of the two genes.

Methodology for clinical genotyping of CYP2D6 and CYP2C19

Many antidepressants, atomoxetine, and several antipsychotics are metabolized by the cytochrome P450 enzymes CYP2D6 and CYP2C19, and guidelines for prescribers based on genetic variants exist. Although some laboratories offer such testing, there is no consensus regarding validated methodology for clinical genotyping of CYP2D6 and CYP2C19. The aim of this paper was to cross-validate multiple technologies for genotyping CYP2D6 and CYP2C19 against each other, and to contribute to feasibility for clinical implementation by providing an enhanced range of assay options, customizable automated translation of data into haplotypes, and a workflow algorithm. AmpliChip CYP450 and some TaqMan single nucleotide variant (SNV) and copy number variant (CNV) data in the Genome-based therapeutic drugs for depression (GENDEP) study were used to select 95 samples (out of 853) to represent as broad a range of CYP2D6 and CYP2C19 genotypes as possible. These 95 included a larger range of CYP2D6 hybrid configurations than have previously been reported using inter-technology data. Genotyping techniques employed were: further TaqMan CNV and SNV assays, xTAGv3 Luminex CYP2D6 and CYP2C19, PharmacoScan, the Ion AmpliSeq Pharmacogenomics Panel, and, for samples with CYP2D6 hybrid configurations, long-range polymerase chain reactions (L-PCRs) with Sanger sequencing and Luminex. Agena MassARRAY was also used for CYP2C19. This study has led to the development of a broader range of TaqMan SNV assays, haplotype phasing methodology with TaqMan adaptable for other technologies, a multiplex genotyping method for efficient identification of some hybrid haplotypes, a customizable automated translation of SNV and CNV data into haplotypes, and a clinical workflow algorithm.

Influence of SULT1A1*2 Polymorphism on Plasma Efavirenz Concentration in Thai HIV-1 Patients

Purpose

Plasma efavirenz (EFV) concentrations within therapeutic levels are essential to successfully treat patients suffering from human immunodeficiency virus (HIV) type 1. In addition to the drug-metabolizing enzyme CYP2B6, other phase II drug-metabolizing enzymes and transporters may have an important role in the pharmacokinetics of EFV. Thus, the influence of phase II drug-metabolizing enzymes and drug transporters on plasma EFV levels was investigated in Thai HIV patients receiving EFV.

Patients and Methods

Genotyping was performed by TaqMan^® real-time PCR in 149 HIV-infected Thai adults, and plasma efavirenz concentration was measured by a validated high-performance liquid chromatography in 12 hours after dosing steady-state plasma samples at week 12 and 24.

Results

Patients with three or more copies of SULT1A1 had significantly lower median plasma EFV concentrations than those carrying two copies at week 12 (p=0.046) and SULT1A1*2 (c.638G>A) carriers had significantly lower median plasma EFV concentrations compared to those not carrying the variant at week 24 (p=0.048). However, no significant association was found after adjusting for CYP2B6 genotype.

Conclusion

Genetic variation in a combination of SULT1A1*2 and SULT1A1 copy number may contribute to variability in EFV metabolism and thereby may impact drug response. The influence of a combination between the SULT1A1 and CYP2B6 genotype on EFV pharmacokinetics should be further investigated in a larger study population.

Resolving discordant CYP2D6 genotyping results in Thai subjects: platform limitations and novel haplotypes

Aim: Several CYP2D6 Luminex xTAG genotype calls were identified as inconsistent or suspicious among Thai subjects and further characterized to identify the root causes. Material & methods: Forty-eight subjects were followed-up with long-range-PCR, quantitative copy number assays and/or Sanger sequencing. Results: Most of the Luminex-duplication calls were either negative or had hybrid structures involving CYP2D6*36 in various configurations. Ten samples were inaccurately called as CYP2D6*2, *29 or *35 alleles. Sequencing revealed three novel haplotypes, CYP2D6*142, *143 and *144 of which two are nonfunctional. Conclusion: The Luminex platform produced a relatively high number of false genotype calls for Thai subjects. Our findings underscore the need for the systematic characterization of the CYP2D6 locus in diverse populations and rigorous platform validation.

CYP2D6 Genetic Variation and Its Implication for Vivax Malaria Treatment in Madagascar

Plasmodium vivax is one of the five human malaria parasite species, which has a wide geographical distribution and can cause severe disease and fatal outcomes. It has the ability to relapse from dormant liver stages (hypnozoites), weeks to months after clearance of the acute blood-stage infection. An 8-aminoquinoline drug primaquine (PQ) can clear the hypnozoites, and thus can be used as an anti-relapse therapeutic agent. Recently, a number of studies have found that its efficacy is compromised by polymorphisms in the cytochrome P450 2D6 (CYP2D6) gene; decreased or absence of CYP2D6 activity contributes to PQ therapeutic failure. The present study sought to characterize CYP2D6 genetic variation in Madagascar, where populations originated from admixture between Asian and African populations, vivax malaria is endemic, and PQ can be deployed soon to achieve national malaria elimination. In a total of 211 samples collected from two health districts, CYP2D6 decreased function alleles CYP2D6*10, *17, *29, *36+*10, and *41 were observed at frequencies of 3.55-17.06%. In addition, nonfunctional alleles were observed, the most common of which were CYP2D6*4 (2.13%), *5 (1.66%), and the *4x2 gene duplication (1.42%). Given these frequencies, 34.6% of the individuals were predicted to be intermediate metabolizers (IM) with an enzyme activity score (AS) ≤ 1.0; both the IM phenotype and AS ≤ 1.0 have been found to be associated with PQ therapeutic failure. Furthermore, the allele and genotype frequency distributions add to the archaeological and genomic evidence of Malagasy populations constituting a unique, Asian-African admixed origin. The results from this exploratory study provide fresh insights about genomic characteristics that could affect the metabolism of PQ into its active state, and may enable optimization of PQ treatment across human genetic diversity, which is critical for achieving P. vivax elimination.

Identification of CYP2D6 Haplotypes that Interfere with Commonly Used Assays for Copy Number Variation Characterization

Cytochrome P450 2D6 (CYP2D6) copy number (CN) variation affects the metabolism of numerous prescribed drugs. Sequence variation within primer or probe target regions of hydrolysis probe CN assays can generate false-positive calls for CN loss. Furthermore, CYP2D6-CYP2D7 hybrids and gene conversions can cause difficult to interpret discordant CN calls. The identification of haplotypes with CN variations and structural arrangements is important to predict phenotype accurately. During clinical testing with hydrolysis probe assays targeting three CYP2D6 regions (intron 2, intron 6, and exon 9), samples with haplotypes causing inconsistent CN calls were identified. To resolve these cases, next-generation sequencing and allele-specific Sanger sequencing was performed. Sequence analysis of 16 samples, all but one from subjects of African descent, identified six novel suballeles containing single-nucleotide polymorphisms, which cause false-positive calls for CN loss in introns 2 and 6. Five samples with an exon 9 CN loss contained CYP2D6/CYP2D7 hybrids (∗13 or ∗36) and one sample was found to have a novel haplotype, CYP2D6∗141. Interestingly, CYP2D6∗141 contains a CYP2D7-derived exon 9 conversion and core single-nucleotide polymorphisms that are otherwise found in CYP2D6∗17 and ∗27. Although these variants are rare, they can cause inconsistent CN calls that typically are reported as no calls or indeterminant, and thus may deprive patients, particularly those of African descent, from taking full benefit of pharmacogenetic testing.

Cyrius: accurate CYP2D6 genotyping using whole-genome sequencing data

Responsible for the metabolism of ~21% of clinically used drugs, CYP2D6 is a critical component of personalized medicine initiatives. Genotyping CYP2D6 is challenging due to sequence similarity with its pseudogene paralog CYP2D7 and a high number and variety of common structural variants (SVs). Here we describe a novel bioinformatics method, Cyrius, that accurately genotypes CYP2D6 using whole-genome sequencing (WGS) data. We show that Cyrius has superior performance (96.5% concordance with truth genotypes) compared to existing methods (84–86.8%). After implementing the improvements identified from the comparison against the truth data, Cyrius’s accuracy has since been improved to 99.3%. Using Cyrius, we built a haplotype frequency database from 2504 ethnically diverse samples and estimate that SV-containing star alleles are more frequent than previously reported. Cyrius will be an important tool to incorporate pharmacogenomics in WGS-based precision medicine initiatives.

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.

Long-Distance Phasing of a Tentative "Enhancer" Single-Nucleotide Polymorphism With CYP2D6 Star Allele Definitions

BACKGROUND

The CYP2D6 gene locus has been extensively studied over decades, yet a portion of variability in CYP2D6 activity cannot be explained by known sequence variations within the gene, copy number variation, or structural rearrangements. It was proposed that rs5758550, located 116 kb downstream of the CYP2D6 gene locus, increases gene expression and thus contributes to variability in CYP2D6 activity. This finding has, however, not been validated. The purpose of the study was to address a major technological barrier, i.e., experimentally linking rs5758550, also referred to as the "enhancer" single-nucleotide polymorphism (SNP), to CYP2D6 haplotypes >100 kb away. To overcome this challenge is essential to ultimately determine the contribution of the "enhancer" SNP to interindividual variability in CYP2D6 activity.

METHODS

A large ethnically mixed population sample (n=3,162) was computationally phased to determine linkage between the "enhancer" SNP and CYP2D6 haplotypes (or star alleles). To experimentally validate predicted linkages, DropPhase2D6, a digital droplet PCR (ddPCR)-based method was developed. 10X Genomics Linked-Reads were utilized as a proof of concept.

RESULTS

Phasing predicted that the "enhancer" SNP can occur on numerous CYP2D6 haplotypes including CYP2D6*1, *2, *5, and *41 and suggested that linkage is incomplete, i.e., a portion of these alleles do not have the "enhancer" SNP. Phasing also revealed differences among the European and African ancestry data sets regarding the proportion of alleles with and without the "enhancer" SNP. DropPhase2D6 was utilized to confirm or refute the predicted "enhancer" SNP location for individual samples, e.g., of n=3 samples genotyped as *1/*41, rs5758550 was on the *41 allele of two samples and on the *1 allele of one sample. Our findings highlight that the location of the "enhancer" SNP must not be assigned by "default." Furthermore, linkage between the "enhancer" SNP and CYP2D6 star allele haplotypes was confirmed with 10X Genomics technology.

CONCLUSIONS

Since the "enhancer" SNP can be present on a portion of normal, decreased, or no function alleles, the phase of the "enhancer" SNP must be considered when investigating the impact of the "enhancer" SNP on CYP2D6 activity.

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.

Normalized Testosterone Glucuronide as a Potential Urinary Biomarker for Highly Variable UGT2B17 in Children 7-18 Years

UDP-glucuronosyltransferase 2B17 (UGT2B17) is a highly variable androgen-metabolizing and drug-metabolizing enzyme. UGT2B17 exhibits a unique ontogeny profile characterized by a dramatic increase in hepatic protein expression from prepubertal age to adulthood. Age, sex, copy number variation (CNV), and single nucleotide polymorphisms only explain 26% of variability in protein expression, highlighting the need for a phenotypic biomarker for predicting interindividual variability in glucuronidation of UGT2B17 substrates. Here, we propose testosterone glucuronide (TG) normalized by androsterone glucuronide (TG/AG) as a urinary UGT2B17 biomarker, and examine the associations among urinary TG/AG and age, sex, and CNV. We performed targeted metabolomics of 12 androgen conjugates with liquid-chromatography tandem mass spectrometry in 63 pediatric subjects ages 7-18 years followed over 7 visits in 3 years. Consistent with the reported developmental trajectory of UGT2B17 protein expression, urinary TG/AG is significantly associated with age, sex, and CNV. In conclusion, TG/AG shows promise as a phenotypic urinary UGT2B17 biomarker.

PharmVar GeneFocus: CYP2D6

The Pharmacogene Variation Consortium (PharmVar) provides nomenclature for the highly polymorphic human CYP2D6 gene locus. CYP2D6 genetic variation impacts the metabolism of numerous drugs and, thus, can impact drug efficacy and safety. This GeneFocus provides a comprehensive overview and summary of CYP2D6 genetic variation and describes how the information provided by PharmVar is utilized by the Pharmacogenomics Knowledgebase (PharmGKB) and the Clinical Pharmacogenetics Implementation Consortium (CPIC).

Copy number variation profiling in pharmacogenes using panel-based exome resequencing and correlation to human liver expression

Structural variants including copy number variations (CNV) have gained widespread attention, especially in pharmacogenomics but for several genes functional relevance and clinical evidence are still lacking. Detection of CNVs in next-generation sequencing data is challenging but offers widespread applications. We developed a cohort-based CNV detection workflow to extract CNVs from read counts of targeted NGS of 340 genes involved in absorption, distribution, metabolism and excretion (ADME) of drugs. We applied our method to 150 human liver tissue samples and correlated identified CNVs to mRNA expression levels. In total, we identified 445 deletions (73%) and 167 duplications (27%) in 36 pharmacogenes including all well-known CNVs of CYPs, GSTs, SULTs, UGTs, numerous described rare CNVs of CYP2E1, SLC16A3 or UGT2B15 as well as novel observations, e.g., for SLC22A12, SLC22A17 and GPS2 (G Protein Pathway Suppressor 2). We were able to fine-map complex CNVs of CYP2A6 and CYP2D6 with exon resolution. Correlation analysis confirmed known expression patterns for common CNVs and suggested an influence on expression variability for some rare CNVs. Our straightforward CNV detection workflow can be easily applied to any NGS coverage data and helped to analyze CNVs in an ADME-NGS panel of 340 pharmacogenes to improve genotype-phenotype correlations.

Polymorphic Expression of UGT1A9 is Associated with Variable Acetaminophen Glucuronidation in Neonates: A Population Pharmacokinetic and Pharmacogenetic Study

INTRODUCTION

Acetaminophen (paracetamol, APAP) is widely used as an analgesic and antipyretic drug in children and neonates. A number of enzymes contribute to the metabolism of acetaminophen, and genetic factors might be important to explain variability in acetaminophen metabolism among individuals.

METHODS

The current investigation utilized a previously published parent-metabolite population pharmacokinetic model describing acetaminophen glucuronidation, sulfation, and oxidation to examine the potential role of genetic variability on the relevant metabolic pathways. Neonates were administered 30-min intravenous infusions of acetaminophen 15 mg/kg every 12 h (< 28 weeks' gestational age [GA]) or every 8 h (≥ 28 weeks GA) for 48 h. A total of 18 sequence variations (SVs) in UDP-glucuronosyltransferase (UGT), sulfotransferase (SULT), and cytochrome P450 (CYP) genes from 33 neonates (aged 1-26 days) were examined in a stepwise manner for an effect on the metabolic formation clearance of acetaminophen by glucuronidation (UGT), sulfation (SULT), and oxidation (CYP). The stepwise covariate modeling procedure was performed using NONMEM^® version 7.3.

RESULTS

Incorporation of genotype as a covariate for one SV located in the UGT1A9 gene promoter region (rs3832043, - 118 > insT, T₉ > T₁₀) significantly improved model fit (likelihood ratio test, p < 0.001) and reduced between-subject variability in glucuronide formation clearance. Individuals with the UGT1A9 T₁₀ polymorphism, indicating insertion of an additional thymidine nucleotide, had a 42% reduction in clearance to APAP-glucuronide as compared to their wild-type counterparts.

CONCLUSION

This study shows a pharmacogenetic effect of an SV in the UGT1A9 promoter region on the metabolism of acetaminophen in neonates.

Multiple relapses of Plasmodium vivax malaria acquired from West Africa and association with poor metabolizer CYP2D6 variant: a case report

Background

Plasmodium vivax transmission in West Africa, dominant for the Duffy-negative blood group, has been increasingly recognized from both local residents as well as international travelers who contracted P. vivax malaria there. However, the relapsing pattern and sensitivity to antimalarial treatment of P. vivax strains originated from this region are largely unknown. There is evidence that the efficacy of primaquine for radical cure of relapsing malaria depends on host factors such as the hepatic enzyme cytochrome P450 (CYP) 2D6.

Case presentation

A 49-year-old Chinese man was admitted to the Shanglin County Hospital in Guangxi Province, China, on December 19, 2016, 39 days after he returned from Ghana, where he stayed for one and a half years. He was diagnosed by microscopy as having uncomplicated P. vivax malaria. Treatment included 3 days of intravenous artesunate (420 mg total), and 3 days of chloroquine (1550 mg total), and 8 days of primaquine (180 mg total). Although parasites and symptoms were cleared rapidly and he was malaria-negative for almost two months, he suffered four relapses with relapse intervals ranging from 58 to 232 days. The last relapse occurred at 491 days from his first vivax attack. For the first three relapses, he was treated similarly with chloroquine and primaquine, sometimes supplemented with additional artemisinin combination therapies (ACTs). For the last relapse, he was treated with intravenous artesunate, 3 days of an ACT, and 7 days of azithromycin, and had remained healthy for 330 days. Molecular studies confirmed P. vivax infections for all the episodes. Although this patient was diagnosed to have normal glucose-6-phosphate dehydrogenase (G6PD) activity, his CYP2D6 genotype corresponded to a *2A/*36 allele variant suggesting of an impaired primaquine metabolizer phenotype.

Conclusions

This clinical case suggests that P. vivax malaria originating from West Africa may produce multiple relapses extending beyond one year. The failures of primaquine as an anti-relapse therapy may be attributed to the patient’s impaired metabolizer phenotype of the CYP2D6. This highlights the importance of knowing the host G6PD and CYP2D6 activities for effective radical cure of relapsing malaria by primaquine.

Characterization of Reference Materials for Genetic Testing of CYP2D6 Alleles: A GeT-RM Collaborative Project

Pharmacogenetic testing increasingly is available from clinical and research laboratories. However, only a limited number of quality control and other reference materials currently are available for the complex rearrangements and rare variants that occur in the CYP2D6 gene. To address this need, the Division of Laboratory Systems, CDC-based Genetic Testing Reference Material Coordination Program, in collaboration with members of the pharmacogenetic testing and research communities and the Coriell Cell Repositories (Camden, NJ), has characterized 179 DNA samples derived from Coriell cell lines. Testing included the recharacterization of 137 genomic DNAs that were genotyped in previous Genetic Testing Reference Material Coordination Program studies and 42 additional samples that had not been characterized previously. DNA samples were distributed to volunteer testing laboratories for genotyping using a variety of commercially available and laboratory-developed tests. These publicly available samples will support the quality-assurance and quality-control programs of clinical laboratories performing CYP2D6 testing.

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.

CYP2D6 genotype analysis of a Thai population: platform comparison

The highly polymorphic CYP2D6 gene locus leads to a wide range of enzyme activity. Since there are limited data for Thai, the major aim was to investigate CYP2D6 genetic variation in a large Thai population (n = 920). CYP2D6 genotyping was performed using four different platforms. Genotype call rates of the Luminex xTAG® and AmpliChip CYP450 test were 96.5% and 87.4%, respectively. Based on Luminex xTAG® data, the most common alleles and genotypes were *1 0 (49.6%), *1 (24.6%), *2 (10.8%), *5 (6.7%), *41 (6.5%) and *1/*10 (23.9%), *10/*10 (21.5%), *2/*10 (9.4%), *5/*10 (6.9%), *10/*41 (5.7%), respectively. Challenges and limitations of the platforms evaluated are discussed. These data add to our knowledge regarding interethnic variability in CYP2D6 activity and contribute to improving drug therapy in the Thai.

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 it 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 non-functional 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 number from read depth.

Results

We applied Stargazer to NGS data of 32 ethnically diverse HapMap trios that were genotyped by TaqMan assays, long-range PCR, quantitative multiplex PCR, High Resolution Melt analysis, and/or Sanger sequencing. CYP2D6 genotyping by Stargazer was 99.0% concordant with data obtained by these methods and showed 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.

Points-to-consider documents: Scientific information on the evaluation of genetic polymorphisms during non-clinical studies and phase I clinical trials in the Japanese population

Pharmacotherapy shows striking individual differences in pharmacokinetics and pharmacodynamics, involving drug efficacy and adverse reactions. Recent genetic research has revealed that genetic polymorphisms are important intrinsic factors for these inter-individual differences. This pharmacogenomic information could help develop safer and more effective precision pharmacotherapies and thus, regulatory guidance/guidelines were developed in this area, especially in the EU and US. The Project for the Promotion of Progressive Medicine, Medical Devices, and Regenerative Medicine by the Ministry of Health, Labour and Welfare, performed by Tohoku University, reported scientific information on the evaluation of genetic polymorphisms, mainly on drug metabolizing enzymes and transporters, during non-clinical studies and phase I clinical trials in Japanese subjects/patients. We anticipate that this paper will be helpful in drug development for the regulatory usage of pharmacogenomic information, most notably pharmacokinetics.

Frequency of CYP2D6 Alleles Including Structural Variants in the United States

The CYP2D6 gene encodes an enzyme important in the metabolism of many commonly used medications. Variation in CYP2D6 is associated with inter-individual differences in medication response, and genetic testing is used to optimize medication therapy. This report describes a retrospective study of CYP2D6 allele frequencies in a large population of 104,509 de-identified patient samples across all regions of the United States (US). Thirty-seven unique CYP2D6 alleles including structural variants were identified. A majority of these alleles had frequencies which matched published frequency data from smaller studies, while eight had no previously published frequencies. Importantly, CYP2D6 structural variants were observed in 13.1% of individuals and accounted for 7% of the total variants observed. The majority of structural variants detected (73%) were decreased-function or no-function alleles. As such, structural variants were found in approximately one-third (30%) of CYP2D6 poor metabolizers in this study. This is the first CYP2D6 study to evaluate, with a consistent methodology, both structural variants and single copy alleles in a large US population, and the results suggest that structural variants have a substantial impact on CYP2D6 function.

Incorporating population-level genetic variability within laboratory models in toxicology: From the individual to the population

Humans respond to chemical exposures differently due to many factors, such as previous and concurrent stressors, age, sex, and genetic background. The vast majority of laboratory-based toxicology studies, however, have not considered the impact of population-level variability within dose-response relationships. The lack of data dealing with the influence of genetic diversity on the response to chemical exposure provides a difficult challenge for risk assessment as individuals within the population will display a wide-range of responses following toxicant challenge. Notably, the genetic background of individuals plays a major role in the variability seen in a population-level response to a drug or chemical and, thus, there is growing interest in including genetic diversity into laboratory-models. Here we outline several laboratory-based models that can be used to assay the influence of genetic variability on an individual's response to chemicals: 1) genetically-diverse cell lines, 2) human primary cells, 3) and genetically-diverse mouse panels. We also provide a succinct review for several seminal studies to highlight the capability, feasibility, and power of each of these models. This article is intended to highlight the need to include population-level genetic diversity into toxicological study designs via laboratory-based models with the goal to provide and supplement evidence in assessing the risk posed by chemicals to the human population. As such, incorporation of genetic variability will positively impact human-based risk assessment and provide empirical data to aid and influence decision-making processes in relation to chemical exposures.

Characterization of ADME gene variation in 21 populations by exome sequencing

OBJECTIVE

Proteins involving absorption, distribution, metabolism, and excretion (ADME) play a critical role in drug pharmacokinetics. The type and frequency of genetic variation in the ADME genes differ among populations. The aim of this study was to systematically investigate common and rare ADME coding variation in diverse ethnic populations by exome sequencing.

MATERIALS AND METHODS

Data derived from commercial exome capture arrays and next-generation sequencing were used to characterize coding variation in 298 ADME genes in 251 Northeast Asians and 1181 individuals from the 1000 Genomes Project.

RESULTS

Approximately 75% of the ADME coding sequence was captured at high quality across the joint samples harboring more than 8000 variants, with 49% of individuals carrying at least one 'knockout' allele. ADME genes carried 50% more nonsynonymous variation than non-ADME genes (P=8.2×10) and showed significantly greater levels of population differentiation (P=7.6×10). Out of the 2135 variants identified that were predicted to be deleterious, 633 were not on commercially available ADME or general-purpose genotyping arrays. Forty deleterious variants within important ADME genes, with frequencies of at least 2% in at least one population, were identified as candidates for future pharmacogenetic studies.

CONCLUSION

Exome sequencing was effective in accurately genotyping most ADME variants important for pharmacogenetic research, in addition to identifying rare or potentially de novo coding variants that may be clinically meaningful. Furthermore, as a class, ADME genes are more variable and less sensitive to purifying selection than non-ADME genes.

Association between CYP2D6 Genotypes and the Risk of Antidepressant Discontinuation, Dosage Modification and the Occurrence of Maternal Depression during Pregnancy

Importance: Polymorphic expression of drug metabolizing enzymes affects the metabolism of antidepressants, and thus can contribute to drug response and/or adverse events. Pregnancy itself can affect CYP2D6 activity with profound variations determined by CYP2D6 genotype. Objective: To investigate the association between CYP2D6 genotype and the risk of antidepressant discontinuation, dosage modification, and the occurrence of maternal CYP2D6, Antidepressants, Depression during pregnancy. Setting: Data from the Organization of Teratology Information Specialists (OTIS) Antidepressants in Pregnancy Cohort, 2006-2010, were used. Women were eligible if they were within 14 completed weeks of pregnancy at recruitment and exposed to an antidepressant or having any exposures considered non-teratogenic. Main Outcomes and Measures: Gestational antidepressant usage was self-reported and defined as continuous/discontinued use, and non-use; dosage modification was further documented. Maternal depression and anxiety were measured every trimester using the telephone interviewer-administered Edinburgh Postnatal Depression Scale and the Beck Anxiety Inventory, respectively. Saliva samples were collected and used for CYP2D6 genotype analyses. Logistic regression models were used to calculate crude and adjusted odds ratios (OR) with 95% confidence intervals. Results: A total of 246 pregnant women were included in the study. The majority were normal metabolizers (NM, n = 204, 83%); 3.3% (n = 8) were ultrarapid metabolizers (UM), 5.7% (n = 14) poor metabolizers (PM), and 8.1% (n = 20) intermediate metabolizers (IM). Among study subjects, 139 women were treated with antidepressants at the beginning of pregnancy, and 21 antidepressant users (15%) discontinued therapy during pregnancy. Adjusting for depressive symptoms, and other potential confounders, the risk of discontinuing antidepressants during pregnancy was nearly four times higher in slow metabolizers (poor or intermediate metabolizers) compared to those with a faster metabolism rate (normal or ultrarapid metabolizers), aOR = 3.57 (95% CI: 1.15-11.11). Predicted CYP2D6 metabolizer status did not impact dosage modifications. Compared with slow metabolizers, significantly higher proportion of women in the fast metabolizer group had depressive symptom in the first trimester (19.81 vs. 5.88%, P = 0.049). Almost 21% of treated women remained depressed during pregnancy (14.4% NM-UM; 6.1% PM-IM). Conclusions and Relevance: Prior knowledge of CYP2D6 genotype may help to identify pregnant women at greater risk of antidepressant discontinuation. Twenty percent of women exposed to antidepressants during pregnancy remained depressed, indicating an urgent need for personalized treatment of depression during pregnancy.

Methyleugenol DNA adducts in human liver are associated with SULT1A1 copy number variations and expression levels

Methyleugenol is a rodent hepatocarcinogen occurring in many herbs and spices as well as essential oils used for flavoring. Following metabolic activation by cytochromes P450 (CYPs) and sulfotransferases (SULTs), methyleugenol can form DNA adducts. Previously, we showed that DNA adduct formation by methyleugenol in mouse liver is dependent on SULT1A1 expression and that methyleugenol DNA adducts are abundant in human liver specimens. In humans, SULT1A1 activity is affected by genetic polymorphisms, including single-nucleotide polymorphisms (SNPs) and copy number variations (CNVs). Here we investigated the relationship between individual methyleugenol DNA adduct levels and SULT1A1 in human liver samples. Using isotope-dilution ultraperformance liquid chromatography coupled with tandem mass spectrometry, we quantified methyleugenol DNA adducts in 121 human surgical liver samples. Frequent CNVs, including deletions (f = 3.3%) and duplications (f = 36.4%) of SULT1A1, were identified using qPCR and TaqMan assays in the donors' genomic DNA. SULT1A1 mRNA and protein levels were quantified using microarray data and Western blot analysis, respectively. Methyleugenol DNA adducts were detected in all 121 liver samples studied. Their levels varied 122-fold between individuals and were significantly correlated to both mRNA and protein levels of SULT1A1 (r _s = 0.43, and r _s = 0.44, respectively). Univariate and multivariate statistical analysis identified significant associations of SULT1A1 CNVs with mRNA (p = 1.7 × 10^-06) and protein (p = 4.4 × 10^- 10) levels as well as methyleugenol DNA adduct levels (p = 0.003). These data establish the importance of SULT1A1 genotype for hepatic methyleugenol DNA adducts in humans, and they confirm a strong impact of SULT1A1 CNVs on SULT1A1 hepatic phenotype.

QTL Mapping and Identification of Candidate Genes in DO Mice: A Use Case Model Derived from a Benzene Toxicity Experiment

Diversity Outbred (DO) mice are a multiparental advanced generation intercross population derived from eight inbred strains which are genetically very diverse. They are maintained as an outbred population using a randomized mating design. Thus DO mice represent an ideal population to map phenotypic traits. Here, we provide a case study in which male DO mice were exposed to benzene and phenotyped for the number of micronucleated reticulocytes. We provide step-by-step R scripts for the analysis of phenotypes, genotypes, mapping of resistance gene loci and identification of candidate genes.

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.

Constellation: a tool for rapid, automated phenotype assignment of a highly polymorphic pharmacogene, CYP2D6, from whole-genome sequences

An important component of precision medicine-the use of whole-genome sequencing (WGS) to guide lifelong healthcare-is electronic decision support to inform drug choice and dosing. To achieve this, automated identification of genetic variation in genes involved in drug absorption, distribution, metabolism, excretion and response (ADMER) is required. CYP2D6 is a major enzyme for drug bioactivation and elimination. CYP2D6 activity is predominantly governed by genetic variation; however, it is technically arduous to haplotype. Not only is the nucleotide sequence of CYP2D6 highly polymorphic, but the locus also features diverse structural variations, including gene deletion, duplication, multiplication events and rearrangements with the nonfunctional, neighbouring CYP2D7 and CYP2D8 genes. We developed Constellation, a probabilistic scoring system, enabling automated ascertainment of CYP2D6 activity scores from 2×100 paired-end WGS. The consensus reference method included TaqMan genotyping assays, quantitative copy-number variation determination and Sanger sequencing. When compared with the consensus reference Constellation had an analytic sensitivity of 97% (59 of 61 diplotypes) and analytic specificity of 95% (116 of 122 haplotypes). All extreme phenotypes, i.e., poor and ultrarapid metabolisers were accurately identified by Constellation. Constellation is anticipated to be extensible to functional variation in all ADMER genes, and to be performed at marginal incremental financial and computational costs in the setting of diagnostic WGS.

Single dose, CYP2D6 genotype-stratified pharmacokinetic study of atomoxetine in children with ADHD

The effect of CYP2D6 genotype on the dose-exposure relationship for atomoxetine has not been well characterized in children. Children 6-17 years of age diagnosed with attention-deficit hyperactivity disorder (ADHD) were stratified by CYP2D6 genotype into groups with 0 (poor metabolizers [PMs], n = 4), 0.5 (intermediate metabolizers [IMs], n = 3), one (extensive metabolizer [EM]1, n = 8) or two (EM2, n = 8) functional alleles and administered a single 0.5 mg/kg oral dose of atomoxetine (ATX). Plasma and urine samples were collected for 24 (IM, EM1, and EM2) or 72 hours (PMs). Dose-corrected ATX systemic exposure (area under the curve [AUC]0-∞ ) varied 29.6-fold across the study cohort, ranging from 4.4 ± 2.7 μM*h in EM2s to 5.8 ± 1.7 μM*h, 16.3 ± 2.9 μM*h, and 50.2 ± 7.3 μM*h in EM1s, IMs, and PMs, respectively (P < 0.0001). Simulated steady state profiles at the maximum US Food and Drug Administration (FDA)-recommended dose suggest that most patients are unlikely to attain adequate ATX exposures. These data support the need for individualized dosing strategies for more effective use of the medication.

CYP2A6 genotyping methods and strategies using real-time and end point PCR platforms

CYP2A6 genotyping is of clinical importance--CYP2A6 gene variants influence nicotine metabolism and are associated with nicotine dependence, cigarettes per day, smoking cessation and the risk for tobacco-associated cancers. CYP2A6 gene variants also influence the metabolism of therapeutic drugs, such as the anticancer agents, tegafur and letrozole. Over the years, CYP2A6 genotyping methods have evolved to incorporate novel gene variants and to circumvent genotyping errors resulting from the high degree of homology between CYP2A6 and neighboring CYP2A genes. Herein, CYP2A6 genotyping strategies are described for commonly genotyped functionally significant alleles including SNPs, small insertions/deletions and more complex structural variants. The methods presented utilize higher throughput SYBR green real-time PCR technology in addition to standard thermocycling.

CYP2D6 Haplotype Determination Using Long Range Allele-Specific Amplification: Resolution of a Complex Genotype and a Discordant Genotype Involving the CYP2D6*59 Allele

Cytochrome P450 (CYP) 2D6, a major contributor to the metabolism and bioactivation of many clinically used drugs, is encoded by a complex, highly polymorphic gene locus. To aid in the characterization of CYP2D6 allelic variation, we developed allele-specific long-range PCR (ASXL-PCR) to amplify only the allele of interest for further characterization by PCR. This development was achieved utilizing single-nucleotide polymorphisms in the upstream region of CYP2D6 and a universal CYP2D6-specific reverse primer. This approach was assessed and optimized on samples with known genotypes. The application of ASXL-PCR clarified a case with a complex genotype (CYP2D6*2x2/*4N+*4) by amplifying the duplicated gene units separately for subsequent analysis. Furthermore, ASXL-PCR and subsequent sequence analysis also resolved genotype discord in a mother/daughter relationship by revealing the presence of the CYP2D6*59 allelic variant in both individuals. Finally, we demonstrated that the 2939G>A single-nucleotide polymorphism present on CYP2D6*59 interfered with the TaqMan genotype assay that detected 2850C>T, causing false genotype assignments. Assay interference was resolved using an alternative TaqMan genotype assay currently available as a custom-made assay. These examples demonstrate the utility of ASXL-PCR for improved CYP2D6 allele/haplotype characterization. This fast, easy-to-perform method is not limited to CYP2D6 but can be adapted to any gene locus for which polymorphic sites are known.

Resolution of a clinical AmpliChip CYP450 Test™ no call: discovery and characterization of novel CYP2D6*1 haplotypes

A Han Chinese patient failed CYP2D6 genotype analysis with the AmpliChip CYP450 Test™. The CYP2D6 gene locus of the patient and her son were extensively genotyped including copy number variation and gene resequencing. Two SNPs were discovered on the patient's CYP2D6*1 allele, -498C>A and 1661G>C, while the son's CYP2D6*1 allele had -498C>A only. AmpliChip failure was attributed to the presence of a CYP2D6*1 allele carrying the 1661G>C SNP. Functional analyses of -498C>A did not reveal altered activity in vitro or in vivo suggesting that both novel CYP2D6*1 subvariants are functional. The implementation of pharmacogenetics-guided drug therapy relies on accurate clinical-grade genotype analysis. Although the AmpliChip is a reliable platform, numerous allelic (sub)variants and gene arrangements are not detected or may trigger no calls. While such cases may be rare, the clinical/genetic testing community must be aware of the challenges of CYP2D6 testing on the AmpliChip platform and implications regarding accuracy of test results.

Diversity Outbred Mice Identify Population-Based Exposure Thresholds and Genetic Factors that Influence Benzene-Induced Genotoxicity

BACKGROUND

Inhalation of benzene at levels below the current exposure limit values leads to hematotoxicity in occupationally exposed workers.

OBJECTIVE

We sought to evaluate Diversity Outbred (DO) mice as a tool for exposure threshold assessment and to identify genetic factors that influence benzene-induced genotoxicity.

METHODS

We exposed male DO mice to benzene (0, 1, 10, or 100 ppm; 75 mice/exposure group) via inhalation for 28 days (6 hr/day for 5 days/week). The study was repeated using two independent cohorts of 300 animals each. We measured micronuclei frequency in reticulocytes from peripheral blood and bone marrow and applied benchmark concentration modeling to estimate exposure thresholds. We genotyped the mice and performed linkage analysis.

RESULTS

We observed a dose-dependent increase in benzene-induced chromosomal damage and estimated a benchmark concentration limit of 0.205 ppm benzene using DO mice. This estimate is an order of magnitude below the value estimated using B6C3F1 mice. We identified a locus on Chr 10 (31.87 Mb) that contained a pair of overexpressed sulfotransferases that were inversely correlated with genotoxicity.

CONCLUSIONS

The genetically diverse DO mice provided a reproducible response to benzene exposure. The DO mice display interindividual variation in toxicity response and, as such, may more accurately reflect the range of response that is observed in human populations. Studies using DO mice can localize genetic associations with high precision. The identification of sulfotransferases as candidate genes suggests that DO mice may provide additional insight into benzene-induced genotoxicity.

A novel simple method for determining CYP2D6 gene copy number and identifying allele(s) with duplication/multiplication

Background

Cytochrome P450 2D6 (CYP2D6) gene duplication and multiplication can result in ultrarapid drug metabolism and therapeutic failure or excessive response in patients. Long range polymerase chain reaction (PCR), restriction fragment length polymorphism (RFLP) and sequencing are usually used for genotyping CYP2D6 duplication/multiplications and identification, but are labor intensive, time consuming, and costly.

Methods

We developed a simple allele quantification-based Pyrosequencing genotyping method that facilitates CYP2D6 copy number variation (CNV) genotyping while also identifying allele-specific CYP2D6 CNV in heterozygous samples. Most routine assays do not identify the allele containing a CNV. A total of 237 clinical and Coriell DNA samples with different known CYP2D6 gene copy numbers were genotyped for CYP2D6 *2, *3, *4, *6, *10, *17, *41 polymorphisms and CNV determination.

Results

The CYP2D6 gene allele quantification/identification were determined simultaneously with CYP2D6*2, *3, *4, *6, *10, *17, *41 genotyping. We determined the exact CYP2D6 gene copy number, identified which allele had the duplication or multiplication, and assigned the correct phenotype and activity score for all samples.

Conclusions

Our method can efficiently identify the duplicated CYP2D6 allele in heterozygous samples, determine its copy number in a fraction of time compared to conventional methods and prevent incorrect ultrarapid phenotype calls. It also greatly reduces the cost, effort and time associated with CYP2D6 CNV genotyping.

Characterization of a complex CYP2D6 genotype that caused an AmpliChip CYP450 Test no-call in the clinical setting

BACKGROUND

CYP2D6, a major drug-metabolizing enzyme, is encoded by a highly polymorphic and complex gene locus. We have identified a patient who failed to produce a CYP2D6 genotype with the AmpliChip P450 Test (AmpliChip), whereas his CYP2C19 genotype was readily determined. The aim of this investigation was to fully characterize the patient's CYP2D6 gene locus to resolve the AmpliChip no-call.

METHODS

The case, a brother, and son were genotyped with the AmpliChip and subsequently characterized using long-range (XL)-PCR coupled with TaqMan assay technology. Copy number variation was assessed by XL-PCR and quantitative PCR. Selected XL-PCR products were sequenced.

RESULTS

The AmpliChip also produced a no-call for the son; the brother produced a result. The two alleles of the case were subsequently found to carry additional gene units that likely caused the AmpliChip no-calls. One was characterized as a CYP2D6*68+*4 tandem (CYP2D6*68 is a hybrid gene composed of 2D6 and 2D7), the other as a rare CYP2D6*13+*2 tandem (CYP2D6*13 is a 2D7/2D6 hybrid formerly known as CYP2D6*77). A novel CYP2D6*2 subvariant was identified in the son; the brother also carried the CYP2D6*68+*4 tandem.

CONCLUSIONS

The implementation of pharmacogenetics-guided drug therapy relies on accurate clinical-grade genotype analysis. Although the AmpliChip is deemed to be a reliable platform, numerous more recently discovered allelic variants and gene arrangements are not detected or trigger no-calls. Although such cases may be rare, the clinical/genetic testing community must be aware of the challenges of CYP2D6 testing on the AmpliChip platform and implications regarding accuracy of test results.

Detection of an endogenous urinary biomarker associated with CYP2D6 activity using global metabolomics

AIM

We sought to discover endogenous urinary biomarkers of human CYP2D6 activity.

PATIENTS & METHODS

Healthy pediatric subjects (n = 189) were phenotyped using dextromethorphan and randomized for candidate biomarker selection and validation. Global urinary metabolomics was performed using liquid chromatography quadrupole time-of-flight mass spectrometry. Candidate biomarkers were tested in adults receiving fluoxetine, a CYP2D6 inhibitor.

RESULTS

A biomarker, M1 (m/z 444.3102) was correlated with CYP2D6 activity in both the pediatric training and validation sets. Poor metabolizers had undetectable levels of M1, whereas it was present in subjects with other phenotypes. In adult subjects, a 9.56-fold decrease in M1 abundance was observed during CYP2D6 inhibition.

CONCLUSION

Identification and validation of M1 may provide a noninvasive means of CYP2D6 phenotyping.

Distribution of CYP2D6 alleles and phenotypes in the Brazilian population

The CYP2D6 enzyme is one of the most important members of the cytochrome P450 superfamily. This enzyme metabolizes approximately 25% of currently prescribed medications. The CYP2D6 gene presents a high allele heterogeneity that determines great inter-individual variation. The aim of this study was to evaluate the variability of CYP2D6 alleles, genotypes and predicted phenotypes in Brazilians. Eleven single nucleotide polymorphisms and CYP2D6 duplications/multiplications were genotyped by TaqMan assays in 1020 individuals from North, Northeast, South, and Southeast Brazil. Eighteen CYP2D6 alleles were identified in the Brazilian population. The CYP2D6*1 and CYP2D6*2 alleles were the most frequent and widely distributed in different geographical regions of Brazil. The highest number of CYPD6 alleles observed was six and the frequency of individuals with more than two copies ranged from 6.3% (in Southern Brazil) to 10.2% (Northern Brazil). The analysis of molecular variance showed that CYP2D6 is homogeneously distributed across different Brazilian regions and most of the differences can be attributed to inter-individual differences. The most frequent predicted metabolic status was EM (83.5%). Overall 2.5% and 3.7% of Brazilians were PMs and UMs respectively. Genomic ancestry proportions differ only in the prevalence of intermediate metabolizers. The IM predicted phenotype is associated with a higher proportion of African ancestry and a lower proportion of European ancestry in Brazilians. PM and UM classes did not vary among regions and/or ancestry proportions therefore unique CYP2D6 testing guidelines for Brazilians are possible and could potentially avoid ineffective or adverse events outcomes due to drug prescriptions.

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.

Impact of metabolizing enzymes on drug response of endocrine therapy in breast cancer

Estrogen-receptor positive breast cancer accounts for 75% of diagnosed breast cancers worldwide. There are currently two major options for adjuvant treatment: tamoxifen and aromatase inhibitors. Variability in metabolizing enzymes determines their pharmacokinetic profile, possibly affecting treatment response. Therefore, prediction of therapy outcome based on genotypes would enable a more personalized medicine approach, providing optimal therapy for each patient. In this review, the authors will discuss the current evidence on the most important metabolizing enzymes in endocrine therapy, with a special focus on CYP2D6 and its role in tamoxifen metabolism.