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

Investigation of the Differences in the Pharmacokinetics of CYP2D6 Substrates, Desipramine, and Dextromethorphan in Healthy African Subjects Carrying the Allelic Variants CYP2D6*17 and CYP2D6*29, When Compared with Normal Metabolizers

This study investigated the differences in the pharmacokinetics (PK) of dextromethorphan and desipramine in healthy African volunteers to understand the effect of allelic variants of the human cytochrome P450 2D6 (CYP2D6) enzyme, namely the diplotypes of CYP2D6*1/*2 (*1*1, *1*2, *2*2) and the genotypes of CYP2D6*17*17 and CYP2D6*29*29. Overall, 28 adults were included and split into 3 cohorts after genotype screening: CYP2D6*1/*2 (n = 12), CYP2D6*17*17 (n = 12), and CYP2D6*29*29 (n = 4). Each subject received a single oral dose of dextromethorphan 30 mg syrup on day 1 and desipramine 50 mg tablet on day 8. The PK parameters of area under the plasma concentration-time curve from time of dosing to time of last quantifiable concentration (AUClast), and extrapolated to infinity (AUCinf), and the maximum plasma concentration (Cmax) were determined. For both dextromethorphan and desipramine, AUCinf and Cmax were higher in subjects of the CYP2D6*29*29 and CYP2D6*17*17 cohorts, as compared with subjects in the CYP2D6*1/*2 diplotype cohort and with normal metabolizers from the literature. All PK parameters, including AUCinf, Cmax, and the elimination half-life, followed a similar trend: CYP2D6*17*17 > CYP2D6*29*29 > CYP2D6*1/*2. The plasma and urinary drug/metabolite exposure ratios of both drugs were higher in subjects of the CYP2D6*17*17 and CYP2D6*29*29 cohorts, when compared with subjects in the CYP2D6*1/*2 diplotype cohort. All adverse events were mild, except in 1 subject with CYP2D6*17*17 who had moderately severe headache with desipramine. These results indicate that subjects with CYP2D6*17*17 and CYP2D6*29*29 genotypes were 5-10 times slower metabolizers than those with CYP2D6*1/*2 diplotypes. These findings suggest that dose optimization may be required when administering CYP2D6 substrate drugs in African patients. Larger studies can further validate these findings.

Determination of the Duplicated CYP2D6 Allele Using Real-Time PCR Signal: An Alternative Approach

CYP2D6 duplication has important pharmacogenomic implications. Reflex testing with long-range PCR (LR-PCR) can resolve the genotype when a duplication and alleles with differing activity scores are detected. We evaluated whether visual inspection of plots from real-time-PCR-based targeted genotyping with copy number variation (CNV) detection could reliably determine the duplicated CYP2D6 allele. Six reviewers evaluated QuantStudio OpenArray CYP2D6 genotyping results and the TaqMan Genotyper plots for seventy-three well-characterized cases with three copies of CYP2D6 and two different alleles. Reviewers blinded to the final genotype visually assessed the plots to determine the duplicated allele or opt for reflex sequencing. Reviewers achieved 100% accuracy for cases with three CYP2D6 copies that they opted to report. Reviewers did not request reflex sequencing in 49-67 (67-92%) cases (and correctly identified the duplicated allele in each case); all remaining cases (6-24) were marked by at least one reviewer for reflex sequencing. In most cases with three copies of CYP2D6, the duplicated allele can be determined using a combination of targeted genotyping using real-time PCR with CNV detection without need for reflex sequencing. In ambiguous cases and those with >3 copies, LR-PCR and Sanger sequencing may still be necessary for determination of the duplicated allele.

Characterization of Novel CYP2D6 Alleles across Sub-Saharan African Populations

The CYP2D6 gene has been widely studied to characterize variants and/or star alleles, which account for a significant portion of variability in drug responses observed within and between populations. However, African populations remain under-represented in these studies. The increasing availability of high coverage genomes from African populations has provided the opportunity to fill this knowledge gap. In this study, we characterized computationally predicted novel CYP2D6 star alleles in 30 African subjects for whom DNA samples were available from the Coriell Institute. CYP2D6 genotyping and resequencing was performed using a variety of commercially available and laboratory-developed tests in a collaborative effort involving three laboratories. Fourteen novel CYP2D6 alleles and multiple novel suballeles were identified. This work adds to the growing catalogue of validated African ancestry CYP2D6 allelic variation in pharmacogenomic databases, thus laying the foundation for future functional studies and improving the accuracy of CYP2D6 genotyping, phenotype prediction, and the refinement of clinical pharmacogenomic implementation guidelines in African and global settings.

Cas9 targeted nanopore sequencing with enhanced variant calling improves CYP2D6-CYP2D7 hybrid allele genotyping

CYP2D6 is a very important pharmacogene as it is responsible for the metabolization or bioactivation of 20 to 30% of the clinically used drugs. However, despite its relatively small length of only 4.4 kb, it is one of the most challenging pharmacogenes to genotype due to the high similarity with its neighboring pseudogenes and the frequent occurrence of CYP2D6-CYP2D7 hybrids. Unfortunately, most current genotyping methods are therefore not able to correctly determine the complete CYP2D6-CYP2D7 sequence. Therefore, we developed a genotyping assay to generate complete allele-specific consensus sequences of complex regions by optimizing the PCR-free nanopore Cas9-targeted sequencing (nCATS) method combined with adaptive sequencing, and developing a new comprehensive long read genotyping (CoLoRGen) pipeline. The CoLoRGen pipeline first generates consensus sequences of both alleles and subsequently determines both large structural and small variants to ultimately assign the correct star-alleles. In reference samples, our genotyping assay confirms the presence of CYP2D6-CYP2D7 large structural variants, single nucleotide variants (SNVs), and small insertions and deletions (INDELs) that go undetected by most current assays. Moreover, our results provide direct evidence that the CYP2D6 genotype of the NA12878 DNA should be updated to include the CYP2D6-CYP2D7 *68 hybrid and several additional single nucleotide variants compared to existing references. Ultimately, the nCATS-CoLoRGen genotyping assay additionally allows for more accurate gene function predictions by enabling the possibility to detect and phase de novo mutations in addition to known large structural and small variants.

During the last decades, the usefulness of personalized medicine has become increasingly apparent. Directly linked to that is the need for accurate genotyping assays to determine the pharmacogenetic profile of patients. Continuing research has led to the development of genotyping assays that perform quite robustly. However, complex genes remain an issue when it comes to determining the complete sequence correctly. An example of such a complex but very important pharmacogene is CYP2D6. Therefore, we developed a genotyping assay in an attempt to generate complete allele-specific consensus sequences of CYP2D6, by optimizing a targeted amplification-free long-read sequencing method and developing a new analysis pipeline. In reference samples, we showed that our genotyping assay performed accurately and confirmed the presence of variants that go undetected by most current assays. However, the implementation of this assay in practice is still hampered as the selected enrichment strategies inherently lead to a low percentage of on-target reads, resulting in low on-target sequencing depths. Further optimization and validation of the assay is thus needed, but definitely worth considering for follow-up research as we already demonstrated the added value for generating more complete genotypes, which on its turn will result in more accurate gene function predictions.

Genetic Variation among Pharmacogenes in the Sardinian Population

Pharmacogenetics (PGx) aims to identify the genetic factors that determine inter-individual differences in response to drug treatment maximizing efficacy while decreasing the risk of adverse events. Estimating the prevalence of PGx variants involved in drug response, is a critical preparatory step for large-scale implementation of a personalized medicine program in a target population. Here, we profiled pharmacogenetic variation in fourteen clinically relevant genes in a representative sample set of 1577 unrelated sequenced Sardinians, an ancient island population that accounts for genetic variation in Europe as a whole, and, at the same time is enriched in genetic variants that are very rare elsewhere. To this end, we used PGxPOP, a PGx allele caller based on the guidelines created by the Clinical Pharmacogenetics Implementation Consortium (CPIC), to identify the main phenotypes associated with the PGx alleles most represented in Sardinians. We estimated that 99.43% of Sardinian individuals might potentially respond atypically to at least one drug, that on average each individual is expected to have an abnormal response to about 17 drugs, and that for 27 drugs the fraction of the population at risk of atypical responses to therapy is more than 40%. Finally, we identified 174 pharmacogenetic variants for which the minor allele frequency was at least 10% higher among Sardinians as compared to other European populations, a fact that may contribute to substantial interpopulation variability in drug response phenotypes. This study provides baseline information for further large-scale pharmacogenomic investigations in the Sardinian population and underlines the importance of PGx characterization of diverse European populations, such as Sardinians.

Influence of CYP2D6, CYP2C19 and CYP3A5 polymorphisms on plasma levels of tamoxifen metabolites in Algerian women with ER+ breast cancer

Background

Tamoxifen, a selective estrogen receptor modulator, is indicated for breast cancer developed in response to estrogen.

Findings

In the current study we explored the relationship between the different variants of CYP2D6, CYP2C19, CYP3A5 and plasma Endoxifen levels in Algerian patients with ER + breast cancer. We further conducted the relationship between the candidate genes and the recurrences rate. Endoxifen levels differed significantly (p < .005) between carriers of two functional alleles and patients genotyped as CYP2D6*10, CYP2D6*17, CYP2D6*41 or CYP2D6*5/*5. Patients with elevated Endoxifen concentrations were significantly more likely to not report recurrences than patients with reduced or nul alleles. Such nul/nul, red/red, and red/nul diplotypes have been associated with a higher rate of recurrences than other genotypes during treatment.

Conclusion

Our findings suggest that the CYP2D6 genotype should be considered in tamoxifen-treated women. While quantitatively, CYP2D6 represents only a minor fraction of the total drug metabolizing capacity of the liver, it is polymorphic and, therefore, may alter the balance of metabolism of tamoxifen toward the activation pathways. Breast cancer patients with the CYP2D6 nul/nul or red/nul diplotype may benefit less from Tamoxifen treatment and are more likely to develop recurrences. Comprehensive CYP2D6 genotyping has a good predictive value for CYP2D6 activity. Common variants in CYP2C19 and CYP3A5 did not have a significant impact on the recurrences in this cohort of patients with ER + breast cancer.

The Identification of Novel CYP2D6 Variants in US Hmong: Results From Genome Sequencing and Clinical Genotyping

Objective: Hmong individuals represent a unique East Asian subpopulation in whom limited information concerning pharmacogenetic variation exists. The objectives of this study were to comprehensively characterize the highly polymorphic CYP2D6 gene in Hmong, estimate allele and phenotype frequencies and to compare results between two testing platforms. Methods: DNA from 48 self-identified Hmong participants were sequenced using a targeted next-generation sequencing (NGS) panel. Star allele calls were made using Astrolabe, manual inspection of NGS variant calls and confirmatory Sanger sequencing. Structural variation was determined by long-range (XL)-PCR and digital droplet PCR (ddPCR). The consensus diplotypes were subsequently translated into phenotype utilizing the activity score system. Clinical grade pharmacogenetic testing was obtained for 12 of the 48 samples enabling an assessment of concordance between the consensus calls and those determined by clinical testing platforms. Results: A total of 13 CYP2D6 alleles were identified. The most common alleles were CYP2D6*10 and its structural arrangements (37.5%, 36/96) and the *5 gene deletion (13.5%, 13/96). Three novel suballeles (*10.007, *36.004, and *75.002) were also identified. Phenotype frequencies were as follows: ultrarapid metabolizers (4.2%, 2/48), normal metabolizers (41.7%, 20/48) and intermediate metabolizers (52.1%, 25/48); none of the 48 participants were predicted to be poor metabolizers. Concordance of diplotype and phenotype calls between the consensus and clinical testing were 66.7 and 50%, respectively. Conclusion: Our study to explore CYP2D6 genotypes in the Hmong population suggests that this subpopulation is unique regarding CYP2D6 allelic variants; also, a higher portion of Hmong participants (50%) are predicted to have an intermediate metabolizer phenotype for CYP2D6 compared to other East Asians which range between 27 and 44%. Results from different testing methods varied considerably. These preliminary findings underscore the importance of thoroughly interrogating unique subpopulations to accurately predict a patient's CYP2D6 metabolizer status.

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.

Rare Variants in RPPH1 Real-Time Quantitative PCR Control Assay Binding Sites Result in Incorrect Copy Number Calls

Real-time quantitative PCR (qPCR) using RPPH1 as a reference gene is a standard method for assessment and validation of genomic copy number variations. However, variants in the reference amplicon may cause errors, which was investigated herein. While conducting copy number variation validations for birth defects research studies, 13 of 1634 specimens with multiple loci that appeared to be present as three copies were unexpectedly detected. This apparent trisomy was hypothesized to be an amplification artifact caused by a variant in the RPPH1 amplicon. Sequencing revealed all 13 individuals carried one of the four different variants within the RPPH1 amplicon. These variants could produce allelic dropout or altered reaction efficiency, causing an inaccurate measurement of copy number. Additional genotyping predicted a low frequency of the most common variant (rs3093876; 14/3562 alleles; minor allele frequency, 0.39%). Laboratories should recognize the potential for inaccurate results when using a single qPCR control assay. Overestimated CFTR and SMN2 copy numbers identified during newborn screening that otherwise would have been incorrectly called were also detected. Variants in reference loci may produce false-negative normal results for test loci when real deletions are present. For clinical laboratories screening for heterozygous deletions for diagnostic testing or prenatal/carrier screening via qPCR, the most cost-effective solution to maximize sensitivity is to run triplex reactions targeting the region of interest with two control 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.

Recommendations for Clinical CYP2D6 Genotyping Allele Selection: A Joint Consensus Recommendation of the Association for Molecular Pathology, College of American Pathologists, Dutch Pharmacogenetics Working Group of the Royal Dutch Pharmacists Association, and the European Society for Pharmacogenomics and Personalized Therapy

The goals of the Association for Molecular Pathology Clinical Practice Committee's Pharmacogenomics (PGx) Working Group are to define the key attributes of pharmacogenetic alleles recommended for clinical testing, and to determine a minimal set of variants that should be included in clinical PGx genotyping assays. This document series provides recommendations on a minimal panel of variant alleles (Tier 1) and an extended panel of variant alleles (Tier 2) that will aid clinical laboratories in designing assays for PGx testing. When developing these recommendations, the Association for Molecular Pathology PGx Working Group considered the functional impact of the variant alleles, allele frequencies in multiethnic populations, the availability of reference materials, as well as other technical considerations with regard to PGx testing. The ultimate goal of this Working Group is to promote standardization of PGx gene/allele testing across clinical laboratories. This document is focused on clinical CYP2D6 PGx testing that may be applied to all cytochrome P450 2D6-metabolized medications. These recommendations are not meant to be interpreted as prescriptive but to provide a reference guide for clinical laboratories that may be either implementing PGx testing or reviewing and updating their existing platform.

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.

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.

Determination of novel CYP2D6 haplotype using the targeted sequencing followed by the long-read sequencing and the functional characterization in the Japanese population

Next-generation sequencing (NGS) has identified variations in cytochrome P450 (CYP) 2D6 associated with drug responses. However, determination of novel haplotypes is difficult because of the short reads generated by NGS. We aimed to identify novel CYP2D6 variants in the Japanese population and predict the CYP2D6 phenotype based on in vitro metabolic studies. Using a targeted NGS panel (PKSeq), 990 Japanese genomes were sequenced, and then novel CYP2D6 haplotypes were determined. K_m, V_max, and intrinsic clearance (V_max/K_m) of N-desmethyl-tamoxifen 4-hydroxylation were calculated by in vitro metabolic studies using cDNA-expressed CYP2D6 proteins. After determination of the CYP2D6 diplotypes, phenotypes of the individuals were predicted based on the in vitro metabolic activities. Targeted NGS identified 14 CYP2D6 variants not registered in the Pharmacogene Variation Consortium (PharmVar) database. Ten novel haplotypes were registered as CYP2D6*128 to *137 alleles in the PharmVar database. Based on the V_max/K_m value of each allele, *128, *129, *130, *131, *132, and *133 were predicted to be nonfunctional alleles. According to the results of the present study, six normal metabolizers (NM) and one intermediate (IM) metabolizers were designated as IM and poor metabolizers (PM), respectively. Our findings provide important insights into novel haplotypes and haplotypes of CYP2D6 and the effects on in vitro metabolic activities.

Nanopore sequencing of the pharmacogene CYP2D6 allows simultaneous haplotyping and detection of duplications

Aim: Long read sequencing offers the promise of overcoming some of the challenges in accurate genotyping of complex genes, along with the advantage of straightforward variant phasing. We have established methods for sequencing and haplotyping of the whole CYP2D6 gene using nanopore sequencing. Materials and methods: 32 samples covering various haplotypes including gene duplication were sequenced on the GridION platform. Results: Haplotypes of 52 alleles matched accurately to known star (*) allele subvariants, with the remaining 12 being assigned as new alleles, or new subvariants of known alleles. Duplicated alleles could be detected by analyzing the allelic balance. Conclusion: Nanopore sequencing of CYP2D6 offers a high throughput method for accurate haplotyping, detection of new variants and determination of duplicated alleles.

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.

CYP2D6 haplotypes with enhancer single-nucleotide polymorphism rs5758550 and rs16947 (*2 allele): implications for CYP2D6 genotyping panels

INTRODUCTION

CYP2D6 metabolizes ∼25% of all clinically used drugs, with numerous genetic polymorphisms affecting enzyme activity and drug response. Clinical utility of current CYP2D6 genotyping is partially compromised the unresolved complex haplotype structure of the CYP2D6 locus. We have identified a distal enhancer single-nucleotide polymorphism rs5758550 that robustly increases CYP2D6 expression, whereas rs16947 (CYP2D6*2), previously considered inert, reduces correct mRNA splicing and expression, thereby affecting presumed activity of other alleles on the *2 haplotype.

OBJECTIVE

This study aims to determine the structure and frequency of haplotypes containing either rs5758550 or rs16947, or both, together with other relevant CYP2D6 alleles, assigning predictive enzyme activity scores to each, and addressing ambiguities in estimating diplotypes in different populations.

METHODS

The structure and frequency of haplotypes containing rs5758550 and/or rs16947 in different populations were determined by using phased genotype data from 'The 1000 Genomes Project'. The assigned haplotype-phenotype relationship was tested by associating assigned CYP2D6 activity score with CYP2D6 enzyme activity in a cohort of 122 human liver microsomes.

RESULTS

Addition of enhancer single-nucleotide polymorphism rs5758550 and *2 to a CYP2D6 panel improves prediction of CYP2D6 activity. Moreover, the haplotype containing rs5758550 and rs16947 predict extensive CYP2D6 activity more accurately than CYP2D6*2A, a surrogate marker for extensive activity.

CONCLUSION

With further studies, the results support possible incorporation of rs5758550 and rs16947 into CYP2D6 biomarker panels for more accurate prediction of CYP2D6 metabolizer status.

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).

Interrogation of CYP2D6 Structural Variant Alleles Improves the Correlation Between CYP2D6 Genotype and CYP2D6-Mediated Metabolic Activity

The cytochrome P450 2D6 (CYP2D6) gene locus is challenging to accurately genotype due to numerous single nucleotide variants and complex structural variation. Our goal was to determine whether the CYP2D6 genotype‐phenotype correlation is improved when diplotype assignments incorporate structural variation, identified by the bioinformatics tool Stargazer, with next‐generation sequencing data. Using CYP2D6 activity measured with substrates dextromethorphan and metoprolol, activity score explained 40% and 34% of variability in metabolite formation rates, respectively, when diplotype calls incorporated structural variation, increasing from 36% and 31%, respectively, when diplotypes did not incorporate structural variation. We also investigated whether the revised Clinical Pharmacogenetics Implementation Consortium (CPIC) recommendations for translating genotype to phenotype improve CYP2D6 activity predictions over the current system. Although the revised recommendations do not improve the correlation between activity score and CYP2D6 activity, perhaps because of low frequency of the CYP2D6*10 allele, the correlation with metabolizer phenotype group was significantly improved for both substrates. We also measured the function of seven rare coding variants: one (A449D) exhibited decreased (44%) and another (R474Q) increased (127%) activity compared with reference CYP2D6.1 protein. Allele‐specific analysis found that A449D is part of a novel CYP2D6*4 suballele, CYP2D6*4.028. The novel haplotype containing R474Q was designated CYP2D6*138 by PharmVar; another novel haplotype containing R365H was designated CYP2D6*139. Accuracy of CYP2D6 phenotype prediction is improved when the CYP2D6 gene locus is interrogated using next‐generation sequencing coupled with structural variation analysis. Additionally, revised CPIC genotype to phenotype translation recommendations provides an improvement in assigning CYP2D6 activity.

CYP2D6 Polymorphisms and the Safety and Gametocytocidal Activity of Single-Dose Primaquine for Plasmodium falciparum

Single-dose primaquine (PQ) clears mature gametocytes and reduces the transmission of Plasmodium falciparum after artemisinin combination therapy. Genetic variation in CYP2D6, the gene producing the drug-metabolizing enzyme cytochrome P450 2D6 (CYP2D6), influences plasma concentrations of PQ and its metabolites and is associated with PQ treatment failure in Plasmodium vivax malaria.

Single-dose primaquine (PQ) clears mature gametocytes and reduces the transmission of Plasmodium falciparum after artemisinin combination therapy. Genetic variation in CYP2D6, the gene producing the drug-metabolizing enzyme cytochrome P450 2D6 (CYP2D6), influences plasma concentrations of PQ and its metabolites and is associated with PQ treatment failure in Plasmodium vivax malaria. Using blood and saliva samples of varying quantity and quality from 8 clinical trials across Africa (n = 1,076), we were able to genotype CYP2D6 for 774 samples (72%). We determined whether genetic variation in CYP2D6 has implications for PQ efficacy in individuals with gametocytes at the time of PQ administration (n = 554) and for safety in glucose-6-phosphate dehydrogenase (G6PD)-deficient individuals treated with PQ (n = 110). Individuals with a genetically inferred CYP2D6 poor/intermediate metabolizer status had a higher gametocyte prevalence on day 7 or 10 after PQ than those with an extensive/ultrarapid CYP2D6 metabolizer status (odds ratio [OR] = 1.79 [95% confidence interval {CI}, 1.10, 2.90]; P = 0.018). The mean minimum hemoglobin concentrations during follow-up for G6PD-deficient individuals were 11.8 g/dl for CYP2D6 extensive/ultrarapid metabolizers and 12.1 g/dl for CYP2D6 poor/intermediate metabolizers (P = 0. 803). CYP2D6 genetically inferred metabolizer status was also not associated with anemia following PQ treatment (P = 0.331). We conclude that CYP2D6 poor/intermediate metabolizer status may be associated with prolonged gametocyte carriage after treatment with single-low-dose PQ but not with treatment safety.

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.

Integrated CYP2D6 interrogation for multiethnic copy number and tandem allele detection

AIM

To comprehensively interrogate CYP2D6 by integrating genotyping, copy number analysis and novel strategies to identify CYP2D6*36 and characterize CYP2D6 duplications.

METHODS

Genotyping of 16 CYP2D6 alleles, multiplex ligation-dependent probe amplification (MLPA) and CYP2D6*36 and duplication allele-specific genotyping were performed on 427 African-American, Asian, Caucasian, Hispanic, and Ashkenazi Jewish individuals.

RESULTS

A novel PCR strategy determined that almost half of all CYP2D6*10 (100C>T) alleles are actually *36 (isolated or in tandem with *10) and all identified duplication alleles were characterized. Integrated results from all testing platforms enabled the refinement of genotype frequencies across all studied populations.

CONCLUSION

The polymorphic CYP2D6 gene requires comprehensive interrogation to characterize allelic variation across ethnicities, which was enabled in this study by integrating multiplexed genotyping, MLPA copy number analysis, novel PCR strategies and duplication allele-specific genotyping.

User considerations in assessing pharmacogenomic tests and their clinical support tools

Pharmacogenomic (PGx) testing is gaining recognition from physicians, pharmacists and patients as a tool for evidence-based medication management. However, seemingly similar PGx testing panels (and PGx-based decision support tools) can diverge in their technological specifications, as well as the genetic factors that determine test specificity and sensitivity, and hence offer different values for users. Reluctance to embrace PGx testing is often the result of unfamiliarity with PGx technology, a lack of knowledge about the availability of curated guidelines/evidence for drug dosing recommendations, and an absence of wide-spread institutional implementation efforts and educational support. Demystifying an often confusing and variable PGx marketplace can lead to greater acceptance of PGx as a standard-of-care practice that improves drug outcomes and provides a lifetime value for patients. Here, we highlight the key underlying factors of a PGx test that should be considered, and discuss the current progress of PGx implementation.

Ten Years' Experience with the CYP2D6 Activity Score: A Perspective on Future Investigations to Improve Clinical Predictions for Precision Therapeutics

The seminal paper on the CYP2D6 Activity Score (AS) was first published ten years ago and, since its introduction in 2008, it has been widely accepted in the field of pharmacogenetics. This scoring system facilitates the translation of highly complex CYP2D6 diplotype data into a patient’s phenotype to guide drug therapy and is at the core of all CYP2D6 gene/drug pair guidelines issued by the Clinical Pharmacogenetics Implementation Consortium (CPIC). The AS, however, only explains a portion of the variability observed among individuals and ethnicities. In this review, we provide an overview of sources in addition to CYP2D6 genotype that contribute to the variability in CYP2D6-mediated drug metabolism and discuss other factors, genetic and non-genetic, that likely contribute to the observed variability in CYP2D6 enzymatic activity.

Sequencing the CYP2D6 gene: from variant allele discovery to clinical pharmacogenetic testing

CYP2D6 is one of the most studied enzymes in the field of pharmacogenetics. The CYP2D6 gene is highly polymorphic with over 100 catalogued star (*) alleles, and clinical CYP2D6 testing is increasingly accessible and supported by practice guidelines. However, the degree of variation at the CYP2D6 locus and homology with its pseudogenes make interrogating CYP2D6 by short-read sequencing challenging. Moreover, accurate prediction of CYP2D6 metabolizer status necessitates analysis of duplicated alleles when an increased copy number is detected. These challenges have recently been overcome by long-read CYP2D6 sequencing; however, such platforms are not widely available. This review highlights the genomic complexities of CYP2D6, current sequencing methods and the evolution of CYP2D6 from allele discovery to clinical pharmacogenetic testing.

In vivo characterization of CYP2D6*12, *29 and *84 using dextromethorphan as a probe drug: a case report

CYP2D6*84 was first described in a Black South African subject, however, its function remains unknown. Astrolabe, a probabilistic scoring tool developed in our laboratory to call genotypes from whole genome sequence, identified CYP2D6*84 in a trio. The father presented with intermediate metabolism when challenged with the CYP2D6 probe drug dextromethorphan (DM/dextrorphan [DX] = 0.0839). Since his second allele, CYP2D6*12, is nonfunctional, the observed activity is derived by CYP2D6*84. This finding suggests that the allele's hallmark P267H causes decreased activity toward DM and that this allele should receive a value of 0.5 for Activity Score calculations. The mother's DM/DX of 0.0543 was consistent with the decreased activity classification of CYP2D6*29. The child, a critically ill neonate, was not phenotyped, but predicted to be a normal metabolizer.

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.

CYP2D7 Sequence Variation Interferes with TaqMan CYP2D6 (*) 15 and (*) 35 Genotyping

TaqMan™ genotyping assays are widely used to genotype CYP2D6, which encodes a major drug metabolizing enzyme. Assay design for CYP2D6 can be challenging owing to the presence of two pseudogenes, CYP2D7 and CYP2D8, structural and copy number variation and numerous single nucleotide polymorphisms (SNPs) some of which reflect the wild-type sequence of the CYP2D7 pseudogene. The aim of this study was to identify the mechanism causing false-positive CYP2D6^*15 calls and remediate those by redesigning and validating alternative TaqMan genotype assays. Among 13,866 DNA samples genotyped by the CompanionDx® lab on the OpenArray platform, 70 samples were identified as heterozygotes for 137Tins, the key SNP of CYP2D6^*15. However, only 15 samples were confirmed when tested with the Luminex xTAG CYP2D6 Kit and sequencing of CYP2D6-specific long range (XL)-PCR products. Genotype and gene resequencing of CYP2D6 and CYP2D7-specific XL-PCR products revealed a CC>GT dinucleotide SNP in exon 1 of CYP2D7 that reverts the sequence to CYP2D6 and allows a TaqMan assay PCR primer to bind. Because CYP2D7 also carries a Tins, a false-positive mutation signal is generated. This CYP2D7 SNP was also responsible for generating false-positive signals for rs769258 (CYP2D6^*35) which is also located in exon 1. Although alternative CYP2D6^*15 and ^*35 assays resolved the issue, we discovered a novel CYP2D6^*15 subvariant in one sample that carries additional SNPs preventing detection with the alternate assay. The frequency of CYP2D6^*15 was 0.1% in this ethnically diverse U.S. population sample. In addition, we also discovered linkage between the CYP2D7 CC>GT dinucleotide SNP and the 77G>A (rs28371696) SNP of CYP2D6^*43. The frequency of this tentatively functional allele was 0.2%. Taken together, these findings emphasize that regardless of how careful genotyping assays are designed and evaluated before being commercially marketed, rare or unknown SNPs underneath primer and/or probe regions can impact the performance of PCR-based genotype assays, including TaqMan. Regardless of the test platform used, it is prudent to confirm rare allele calls by an independent method.

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.