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

Association of Polymorphic Cytochrome P450 Enzyme Pathways With Falls in Multimedicated Older Adults

OBJECTIVES

Dose exposure is considered relevant for drug-associated falls in older adults, pointing to an importance of drug metabolism. Aim was to analyze individual factors altering drug metabolism such as enzyme saturation by drug exposure and pharmacogenetics in the context of drug-associated falls.

DESIGN

Prospective population-based study (ActiFE-Ulm study).

SETTING AND PARTICIPANTS

Community-dwelling older adults.

METHODS

Focus was laid on the metabolism by polymorphic cytochrome P450 (CYP) enzymes CYP2C19, 2C9, and 2D6. Relevant variants of pharmacogenes were analyzed. Logistic binary regression analysis was used to calculate odds ratios (ORs) and 95% CIs for falls observed prospectively over a 1-year period with drug metabolism characteristics.

RESULTS

In total, 1377 participants were included in the analysis. Although the phenotype predicted by the genotype was not, the use of drugs metabolized by CYP2C19 was associated with falls. Drugs not known as fall risk-increasing drugs (FRIDs; ie, non-FRIDs), but metabolized by CYP2C19, showed an OR of 1.46 (1.11-1.93) in adjusted analysis. Significant effect modification was observed for a reduced CYP2C19 activity phenotype with non-FRIDs metabolized by CYP2C19.

CONCLUSIONS AND IMPLICATIONS

This study suggests an association between the occurrence of falls in older adults and the metabolic capacity of CYP2C19. Thus, an important step toward prevention of falls might be to personalize dosage and treatment length of the main drug classes known to be CYP2C19 substrates, such as many antidepressants, opioids, and sedatives, but also proton pump inhibitors in particular in poor and intermediate metabolizers.

Prevalence of CYP2D6 structural variation in large retrospective study

CYP2D6 is a highly polymorphic gene with clinically important structural variations. Commonly, only exon 9 is assayed on clinical pharmacogenomics panels, as it allows for accurate functional characterization even in the presence of a CYP2D6::CYP2D7 conversion. However, this method does not capture CYP2D7::CYP2D6 (CYP2D6*13) conversions, possibly leading to inaccurate phenotype assignment. The study's purpose was to determine the frequency of structural variations in CYP2D6 utilizing multiple copy number variation (CNV) assay locations to quantify the potential impact on clinical phenotype classification. A retrospective analysis was conducted of de-identified pharmacogenomics data submitted through the Translational Software, Inc. platform. Samples with CYP2D6 CNV data for exon 9 and at least one additional CNV location (5'UTR, exon 1, intron 2, exon 5 or intron 6) were included. CYP2D7::CYP2D6 and CYP2D6::CYP2D7 conversions were classified according to PharmVar nomenclature. The CYP2D6 copies were capped at four total copies to account for assay limitations in detecting more than four copies. A total of 106,474 samples were included for analysis. CYP2D7::CYP2D6 conversions were present in approximately 2.44% of samples, and 5.84% of samples had CYP2D6::CYP2D7 conversions. Many samples did not have a CYP2D7 conversion detected (91.5%; 97,462/106,474). A full gene deletion was detected in 0.15%, and 5.98% had a duplication or multiplication present. This retrospective study underscores the importance of testing more than one CNV site for CYP2D6 . Over 2% of patients were found to have a CYP2D7::CYP2D6 conversion. This translates into potentially misclassified phenotype classification and incongruent clinical recommendations.

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.

Characterization of complex structural variation in the CYP2D6-CYP2D7-CYP2D8 gene loci using single-molecule long-read sequencing

Complex regions in the human genome such as repeat motifs, pseudogenes and structural (SVs) and copy number variations (CNVs) present ongoing challenges to accurate genetic analysis, particularly for short-read Next-Generation-Sequencing (NGS) technologies. One such region is the highly polymorphic CYP2D loci, containing CYP2D6, a clinically relevant pharmacogene contributing to the metabolism of >20% of common drugs, and two highly similar pseudogenes, CYP2D7 and CYP2D8. Multiple complex SVs, including CYP2D6/CYP2D7-derived hybrid genes are known to occur in different configurations and frequencies across populations and are difficult to detect and characterize accurately. This can lead to incorrect enzyme activity assignment and impact drug dosing recommendations, often disproportionally affecting underrepresented populations. To improve CYP2D6 genotyping accuracy, we developed a PCR-free CRISPR-Cas9 based enrichment method for targeted long-read sequencing that fully characterizes the entire CYP2D6-CYP2D7-CYP2D8 loci. Clinically relevant sample types, including blood, saliva, and liver tissue were sequenced, generating high coverage sets of continuous single molecule reads spanning the entire targeted region of up to 52 kb, regardless of SV present (n = 9). This allowed for fully phased dissection of the entire loci structure, including breakpoints, to accurately resolve complex CYP2D6 diplotypes with a single assay. Additionally, we identified three novel CYP2D6 suballeles, and fully characterized 17 CYP2D7 and 18 CYP2D8 unique haplotypes. This method for CYP2D6 genotyping has the potential to significantly improve accurate clinical phenotyping to inform drug therapy and can be adapted to overcome testing limitations of other clinically challenging genomic regions.

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 CYP2D6 Pharmacogenetic Variation in Sub-Saharan African Populations

Cytochrome P450 2D6 (CYP2D6) is a key enzyme in drug response owing to its involvement in the metabolism of ~ 25% of clinically prescribed medications. The encoding CYP2D6 gene is highly polymorphic, and many pharmacogenetics studies have been performed worldwide to investigate the distribution of CYP2D6 star alleles (haplotypes); however, African populations have been relatively understudied to date. In this study, the distributions of CYP2D6 star alleles and predicted drug metabolizer phenotypes-derived from activity scores-were examined across multiple sub-Saharan African populations based on bioinformatics analysis of 961 high-depth whole genome sequences. This was followed by characterization of novel star alleles and suballeles in a subset of the participants via targeted high-fidelity Single-Molecule Real-Time resequencing (Pacific Biosciences). This study revealed varying frequencies of known CYP2D6 alleles and predicted phenotypes across different African ethnolinguistic groups. Twenty-seven novel CYP2D6 star alleles were predicted computationally and two of them were further validated. This study highlights the importance of studying variation in key pharmacogenes such as CYP2D6 in the African context to better understand population-specific allele frequencies. This will aid in the development of better genotyping panels and star allele detection approaches with a view toward supporting effective implementation of precision medicine strategies in Africa and across the African diaspora.

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.

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.