CYP2D6 Protein Level Is the Major Contributor to Interindividual Variability in CYP2D6-Mediated Drug Metabolism in Healthy Human Liver Tissue

A frequent CYP2D6 variant promotes skipping of exon 3 and reduces CYP2D6 protein expression in human liver samples

CYP2D6 is one of the most polymorphic drug-metabolizing enzymes in the liver. While genetic CYP2D6 variants serve as clinical biomarkers to predict CYP2D6 activity, large inter-person variability in CYP2D6 expression remains unaccounted for. Previous results suggest that there is variable expression of a CYP2D6 splice isoform with an in-frame deletion of exon 3 (CYP2D6ΔE3) encoding a protein lacking numerous active site residues. Here, using fragment analysis and RT-qPCR, we revealed that rs1058164 G (MAF = 27%-43%) is associated with increased formation of CYP2D6∆E3 in human liver samples (1.4-2.5-fold) and transfected cells. Furthermore, western blots showed that rs1058164 G was associated with a 50% decrease in full-length hepatic CYP2D6 protein expression. In addition, by studying a larger liver cohort, we confirmed our previous results that rs16947 (CYP2D6*2) reduces full-length CYP2D6 mRNA by increasing the production of an unstable splice isoform lacking exon 6 (CYP2D6ΔE6) and that the impact of CYP2D6ΔE6 is offset in carriers of the downstream enhancer variant rs5758550. The three frequent SNPs (rs1058164, rs16947, and rs5758550) form various 3-SNP-haplotypes, each with distinct CYP2D6 expression characteristics. Using an expression score (ES) system, we tested the impact of the 3-SNP-haplotype on improving the standard model to predict hepatic CYP2D6 protein expression based on genotype. A model that incorporates the 3-SNP-haplotype provided the best fit for CYP2D6 expression and also accounted for more variability in CYP2D6 protein levels (59%) than a model based on the accepted standard (36%) or one that only adds rs16947 and rs5758550 (42%). Clinical studies are needed to determine whether including the 3-SNP-haplotype alongside current standard CYP2D6 models improves the predictive value of CYP2D6 panels.

Categorization of Cytochrome P4502D6 Activity Score by Urinary Amphetamine/Methamphetamine Ratios

Methamphetamine (MA) level in urine has been used for judgment in MA consumption. Metabolism and intoxication of MA are correlated with the activity of cytochrome P450 2D6 (CYP2D6). The activity score (AS) is a potential tool for predicting exposure and personalized dose of drugs metabolized by CYP2D6. Prediction of the CYP2D6 activity score might be described as MA intoxication. The objective of this study was to categorize the CYP2D6 activity score using the urinary amphetamine (AM)/MA ratio. Urine samples (n = 23,258) were collected. The levels of MA and AM were determined by a gas chromatography-nitrogen-phosphorus detector. The log AS was calculated by an AM/MA ratio and classified into four groups following the percentile position: lower than the 2.5th, the 2.5th-the 50th, the 50th-97.5th, and greater than the 97.5th percentile, respectively. The AS value for males presented was less than 0.024, 0.024-0.141, 0.141-0.836, and greater than 0.836. Meanwhile, the AS values were revealed to be lower than 0.023, 0.023-0.148, 0.148-0.850, and higher than 0.850 for females. The AS value of CYP2D6 can be applied to describe the toxicity of MA in forensic crime scenes and relapse behavior.

From Croatian Roma to 1000 Genomes: The Story of the CYP2D6 Gene Promoter and Enhancer SNPs

The CYP2D6 gene encodes an enzyme responsible for the metabolism of ~20% of clinically prescribed drugs. In this study, 18 SNPs from the enhancer and promoter regions of CYP2D6 in 323 Roma from Croatia were genotyped, to find out whether the demographic history of Roma affected the distribution of the studied SNPs and their linkage disequilibrium (LD) values, with the major SNPs defining the CYP2D6 star alleles. No differences were found between the three Roma groups in allele and genotype frequencies. The distribution of LD values of Roma was compared with LD values of European and Asian populations. Regulatory CYP2D6 SNPs (rs5758550, rs28624811, rs1080985 and rs1080983) showed similar distribution and the highest LDs with rs16947 from the gene-coding region in all populations. In the promoter region, a complete LD between rs1080989 and rs28588594, and between rs1080983 and rs28624811, was found in Croatian Roma and investigated populations from 1000 genomes. A high LD was also found between rs1080985 from the promoter and rs5758550 from the enhancer region. SNP rs28735595 from the gene promoter region had the highest LD, with two gene region SNPs, rs1058164 and rs1135840. To conclude, the Croatian Roma population shows an LD pattern of the CYP2D6 gene region similar to the 1000 Genomes European and Asian populations.

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.

Monitoring Plasmodium vivax resistance to antimalarials: Persisting challenges and future directions

Emerging antimalarial drug resistance may undermine current efforts to control and eliminate Plasmodium vivax, the most geographically widespread yet neglected human malaria parasite. Endemic countries are expected to assess regularly the therapeutic efficacy of antimalarial drugs in use in order to adjust their malaria treatment policies, but proper funding and trained human resources are often lacking to execute relatively complex and expensive clinical studies, ideally complemented by ex vivo assays of drug resistance. Here we review the challenges for assessing in vivo P. vivax responses to commonly used antimalarials, especially chloroquine and primaquine, in the presence of confounding factors such as variable drug absorption, metabolism and interaction, and the risk of new infections following successful radical cure. We introduce a simple modeling approach to quantify the relative contribution of relapses and new infections to recurring parasitemias in clinical studies of hypnozoitocides. Finally, we examine recent methodological advances that may render ex vivo assays more practical and widely used to confirm P. vivax drug resistance phenotypes in endemic settings and review current approaches to the development of robust genetic markers for monitoring chloroquine resistance in P. vivax populations.

•

Plasmodium vivax resistance to chloroquine may undermine malaria elimination efforts.

•

Plasmodium vivax resistance to schizontocides has been mostly monitored in therapeutic efficacy studies.

•

In vivo studies to determine the anti-relapse efficacy of primaquine are challenging to design and execute.

•

Ex vivo assays to determine Plasmodium vivax resistance to schizontocides remain limited to research settings.

•

Robust molecular markers to monitor Plasmodium vivax drug resistance are currently lacking.

Tri-Allelic Haplotypes Determine and Differentiate Functionally Normal Allele CYP2D6*2 and Impaired Allele CYP2D6*41

CYP2D6 metabolizes 20-25% of all clinically used drugs and its complex genetic polymorphism is a major determinant of drug safety and efficacy. We investigated the basis for the functional difference between the two common alleles *2 (g.2851C>T + g.4181G>C, normal function) and *41 (additional intronic g.2989G>A, reduced function). A recently reported far-distant enhancer polymorphism rs5758550A/G linked to *2 has been suggested to play a decisive role. Genotyping of two white cohorts confirmed strong linkage of rs5758550G to *2, whereas no influence was found on metabolic ratio of sparteine or hepatic expression. Genomic plasmid constructs carrying individual variants or combinations thereof were expressed in COS1 and Huh7 cells. Both g.2851C>T(R296C) and g.2989G>A reduced enzyme activity and protein levels similarly by ~ 50-65% compared to reference (*1), whereas the double variant had only ~ 20% activity. Although the unexpected loss of function caused by g.2851C>T was compensated by g.4181G>C (mimicking the EM-phenotype of *2), the additional loss of function due to intronic g.2989G>A in the triple variant was not compensated (mimicking the IM-phenotype of *41). We also confirmed increased erroneous splicing in carriers of g.2989G>A but not of g.2851C>T as a likely explanation for the impaired function of *41. In conclusion, our data demonstrate g.2989G>A as causal variant of impaired allele CYP2D6*41 whereas triple-haplotypes have to be considered to explain the functional difference between *2 and *41. These data are important for genotyping strategies and clinical implementation of CYP2D6 pharmacogenetics.

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.

Determinants of Cytochrome P450 2D6 mRNA Levels in Healthy Human Liver Tissue

Cytochrome P450 2D6 (CYP2D6) is a major drug‐metabolizing enzyme that exhibits large interindividual variability. Recent studies suggest that differential transcriptional regulation of CYP2D6 in part may be responsible for the variability. In this study, we characterized potential determinants of CYP 2D6  transcript levels in healthy human liver tissue samples (n = 115), including genetic polymorphisms in CYP2D6 and the genes encoding transcription regulators for CYP2D6 expression; mRNA expression of the transcription factors and their known target genes; and hepatic levels of bile acids and retinoids, agents that modulate the expression/activity of the transcription factors. Their associations with CYP2D6 mRNA levels in the tissues were examined. Results from multivariable linear regression analysis revealed CYP8B1 mRNA level and rs3892097, the single‐ nucleotide polymorphism defining the nonfunctional CYP2D6*4 allele, as the two most significant predictors of CYP2D6 mRNA levels in the liver tissue samples, explaining 30% of the variability.

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.