Single nucleotide and structural variants of CYP2D6 gene in Kinh Vietnamese population

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.

CYP2D6 and CYP2C19 Genes Associated with Tricontinental and Latin American Ancestry of Pe-ruvians

Precision medicine seeks to individualize the dose from the beginning of phar-macological therapy based on the characteristics of each patient, genes involved in the metabolic phenotype, ethnicity or miscegenation, with the purpose to minimize adverse effects and optimize drug efficacy. The objective was to re-view studies that describe the association of the CYP2D6 and CYP2C19 genes with the tricontinental and Latin American ancestry of Peruvians. A biblio-graphic search was carried out in PubMed/Medline and SciELO, with various descriptors in Spanish and English. The results of this review confirm that the ethnic origin of Peruvians is triconti-nental due to European (mainly Spanish), African and Asian migration, in addi-tion to Latin American migration, being 60.2% mixed, 25.8% Amerindian, 5.9% white, 3.6% African descent, 1.2% Chinese and Japanese descent, and 3.3% unspecified. Studies on CYP2C19*3, CYP2D6*2, *3 and *6 have been reported in Peruvians, and the frequency is similar to that studied in Ecuadori-ans and Colombians. The CYP2C19*3, CYP2D6*3, and CYP2D6*6 alleles found in Peruvians are common in Europeans, Africans, and Asians; while CYP2D6*4 in Africans and CYP2D6*2 related to Asians. In some studies, the ethnic/gene association has not been demonstrated; while others have shown a significant association, which is why further investigation is warranted. It is concluded that the studies on CYP2D6 and CYP2C19 genes associated with the tricontinental and Latin American ancestry of Peruvians are little, and ac-cording to what has been investigated, the CYP2C19*3, CYP2D6*2, *3, *4 and *6 alleles have more related to their ancestry.

Global perspectives on CYP2D6 associations with primaquine metabolism and Plasmodium vivax radical cure

Clinical trial and individual patient treatment outcomes have produced accumulating evidence that effective primaquine (PQ) treatment of Plasmodium vivax and P. ovale liver stage hypnozoites is associated with genetic variation in the human cytochrome P450 gene, CYP2D6. Successful PQ treatment of individual and population-wide infections by the Plasmodium species that generate these dormant liver stage forms is likely to be necessary to reach elimination of malaria caused by these parasites globally. Optimizing safe and effective PQ treatment will require coordination of efforts between the malaria and pharmacogenomics research communities.

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.

Exploring the Kinh Vietnamese genomic database for the polymorphisms of the P450 genes toward precision public health

BACKGROUND

Human cytochrome P450 (CYPs) genes are essential in metabolizing drugs. Due to their high polymorphism, population-specific studies are of great interest.

AIM

This research examined the six CYP genes, including CYP2B6, CYP2C9, CYP2C19, CYP2D6, CYP3A5, and CYP4F2 in the Kinh Vietnamese (KHV) for population-scale precision medicine.

SUBJECTS AND METHODS

We processed data from a genomics database of 206 healthy and unrelated KHV individuals to calculate CYP allele frequencies. First, we compared the CYP genes of the KHV to six other populations retrieved from the 1000 Genomes Project. Second, we searched the PharmGBK database for drug-CYP interaction data to compile a drug dosage recommendation for KHV.

RESULTS

We observed diverging trends in the genetic variations of CYP2B6, CYP2D6, and CYP3A5 in KHV. In terms of the phenotypic drug responses in KHV, CYP2C19 exhibited all of the metabolic phenotypes at a non-trivial frequency. CYP3A5 metabolized drugs at a lower rate than the other five CYPs.

CONCLUSION

This is the first large-scale study to investigate multiple CYP genes in the KHV for precision medicine from a public health perspective. Differences found in the distributions of metabolizers for the KHV suggest careful prescriptions for CYP2C19 and CYP3A5-metabolized drugs.

Genophenotypic Factors and Pharmacogenomics in Adverse Drug Reactions

Adverse drug reactions (ADRs) rank as one of the top 10 leading causes of death and illness in developed countries. ADRs show differential features depending upon genotype, age, sex, race, pathology, drug category, route of administration, and drug–drug interactions. Pharmacogenomics (PGx) provides the physician effective clues for optimizing drug efficacy and safety in major problems of health such as cardiovascular disease and associated disorders, cancer and brain disorders. Important aspects to be considered are also the impact of immunopharmacogenomics in cutaneous ADRs as well as the influence of genomic factors associated with COVID-19 and vaccination strategies. Major limitations for the routine use of PGx procedures for ADRs prevention are the lack of education and training in physicians and pharmacists, poor characterization of drug-related PGx, unspecific biomarkers of drug efficacy and toxicity, cost-effectiveness, administrative problems in health organizations, and insufficient regulation for the generalized use of PGx in the clinical setting. The implementation of PGx requires: (i) education of physicians and all other parties involved in the use and benefits of PGx; (ii) prospective studies to demonstrate the benefits of PGx genotyping; (iii) standardization of PGx procedures and development of clinical guidelines; (iv) NGS and microarrays to cover genes with high PGx potential; and (v) new regulations for PGx-related drug development and PGx drug labelling.

Prediction of the CYP2D6 enzymatic activity based on investigating of the CYP2D6 genotypes around the vivax malaria patients in Yunnan Province, China

Background

In recent years, the incidence rate of vivax malaria recurrence still had 3.1% in Yunnan Province population after eradication therapy using primaquine (PQ). In order to understand the specific failure reasons for preventing vivax malaria relapses, a preliminary exploration on the CYP2D6 enzyme activity was carried out in the vivax malaria patients in Yunnan Province population by analysing mutational polymorphism in the coding region of CYP2D6 gene.

Methods

Blood samples were collected from vivax malaria patients with suspected relapse (SR) and non-relapsed (NR) malaria in Yunnan Province. The DNA fragments containing 9 exons regions of human CYP2D6 gene were amplified by performing PCR and sequenced. The sequencing results were aligned by using DNAStar 11.0 to obtain the coding DNA sequence (CDS) of CYP2D6 gene. DnaSP 6.11.01 software was used to identify mutant polymorphisms and haplotypes of the CDS chain. The waterfall function of GenVisR package in R was utilized to visualize the mutational landscape. The alleles of CYP2D6 gene were identified according to the criteria prescribed by Human Cytochrome P450 (CYP) Allele Nomenclature Committee Database and the CYP2D6 enzyme activity was predicted based on diploid genotype.

Results

A total of 320 maternal CDS chains, including 63 from SR group and 257 from NR group, were obtained. Twelve mutant loci, including c.31 (rs769259), c.100 (rs1065852), c.271 (rs28371703), c.281 (rs28371704), c.294 (rs28371705), c.297 (rs200269944), c.336 (rs1081003), c.408 (rs1058164), c.505 (rs5030865), c.801 (rs28371718), c.886 (rs16947), and c.1,457 (rs1135840) were observed on the 640 CDS chains (including 320 maternal and 320 paternal chains). The high-frequency mutation at rs1135840 (0.703) and low-frequency mutation, such as rs28371703, were detected only in the SR group. The frequency of mutant rs1058164 and rs1135840 were significantly increased in the SR group (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${x}^{2}$$\end{document}x2= 4.468, 5.889, P < 0.05), as opposed to the NR group. Of the 23 haplotypes (from Hap_1 to Hap_23), the nomenclatures of 11 allelic forms could be found: Hap_3 was non-mutant, Hap_2 accounted for the highest frequency (36.9%, 236/640), and Hap_9 had the most complex sequence structure, containing 7 loci mutations. Allele *10 was the most frequent among these genotypes (0.423). Among the allele *10 standard named genotypes, *1/*10, *1/*1 and *2/*10 were significantly more frequent in the NR group (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${x}^{2}$$\end{document}x2= 3.911, P < 0.05) and all showed uncompromised enzyme activity; the impaired genotype *10/*39 was more frequent in the SR group (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${x}^{2}$$\end{document}x2= 10.050, P < 0.05), and genotype *4/*4was detected only in the SR group.

Conclusion

In the patients receiving PQ dosage in Yunnan Province population, both rs1135840 single nucleotide polymorphism and *10 allele form was common in the CYP2D6 gene. Low-frequency mutation sites, such as rs28371703, were only presented in patients with vivax malaria relapse.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03988-5.

The association of CYP2D6 gene polymorphisms in the full-length coding region with higher recurrence rate of vivax malaria in Yunnan Province, China

Background

Accumulating evidence suggest that compromised CYP2D6 enzyme activity caused by gene mutation could contribute to primaquine failure for the radical cure of vivax malaria. The current study aims to preliminarily reveal the association between the recurrence of vivax malaria in Yunnan Province and CYP2D6 gene mutation by analysing polymorphisms in the entire coding region of human CYP2D6 gene.

Methods

Blood samples were collected from patients with vivax malaria, who received "chloroquine and 8-day course of primaquine therapy" in Yunnan Province. The suspected relapsed cases were determined by epidemiological approaches and gene sequence alignment. PCR was conducted to amplify the CYP2D6 gene in the human genome, and the amplified products were then sequenced to compare with the non-mutation “reference” sequence, so as to ensure correct sequencing results and to determine 9 exon regions. Subsequently, the DNA sequences of 9 exons were spliced into the coding DNA sequence (CDS), which, by default, is known as maternal CDS. The paternal CDS was obtained by adjusting the bases according to the sequencing peaks. The mutation loci, haplotypes (star alleles), genotypes and odds ratios (OR) of all the CDSs were analysed.

Results

Of the119 maternal CDS chains in total with 1491 bp in length, 12 mutation sites in the 238 maternal and paternal CDS chains were detected. The c.408G > C mutation was most frequently detected in both suspected relapsed group (SR) and non-relapsed group (NR), reaching 85.2% (75/88) and 76.0% (114/150), respectively. The c.886C > T mutation was most closely related to the recurrence of vivax malaria (OR = 2.167, 95% CI 1.104–4.252, P < 0.05). Among the 23 haplotypes (Hap_1 ~ Hap_23), Hap_3 was non-mutant, and the sequence structure of Hap_9 was the most complicated one. Five star alleles, including *1, *2, *4, *10 and *39, were confirmed by comparison, and CYP2D6*10 allele accounted for the largest percentage (45.4%, 108/238). The frequency of CYP2D6*2 allele in the SR group was significantly higher than that in the NR group (Χ² = 16.177, P < 0.05). Of the defined 24 genotypes, 8 genotypes, including *4/*4, *4/*o, *2/*39, *39/*m, *39/*x, *1/*r, *1/*n, and *v/*10, were detected only in the SR group.

Conclusion

Mutation of CYP2D6*10 allele accounts for the highest proportion of vivax malaria cases in Yunnan Province. The mutations of c. 886C > T and CYP2D6*2 allele, which correspond to impaired PQ metabolizer phenotype, are most closely related to the relapse of vivax malaria. In addition, the genotype *4/*4 with null CYP2D6 enzyme function was only detected in the SR group. These results reveal the risk of defected CYP2D6 enzyme activity that diminishes the therapeutic effect of primaquine on vivax malaria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03685-3.

Relationship between CYP2D6 genotype, activity score and phenotype in a pediatric Thai population treated with risperidone

Recently, the Clinical Pharmacogenetics Implementation Consortium (CPIC) have revised recommendations for the translation of CYP2D6 genotype to phenotype. Changes affect phenotype grouping, as well as the value used to calculate activity score for the CYP2D6*10 allele to better reflect the substantially decreased activity of this allele which is the most frequent allele found in Asian populations. This study aimed to evaluate whether the lower value for CYP2D6*10 as recommended, and the revised phenotype groupings improve the relationship between CYP2D6 genotype and risperidone measures. One hundred and ninety-nine children and adolescents with autism treated with a risperidone-based regimen for at least four weeks were included. CYP2D6 genotype was determined using the Luminex xTAG CYP2D6 Kit assay and translated into phenotype using different translation methods. Plasma concentrations of risperidone and 9-hydroxyrisperidone were measured using LC/MS/MS. Plasma levels of risperidone, risperidone concentration/dose ratio, and risperidone/9-hydroxyrisperidone ratio in patients with an activity score < 1 were significantly higher than those ≥ 1 (P value < 0.001 for all three parameters). Plasma risperidone levels and risperidone concentration/dose ratios were significantly higher in intermediate metabolizers (defined as AS = 0.25–0.75) than normal metabolizer (defined as AS = 1–2) patients (1.44 vs. 0.23 ng/ml, P < 0.001 and 1.63 vs. 0.29 ng/ml/ng, P < 0.001, respectively) as well as risperidone/9-hydroxyrisperidone ratio (0.20 vs. 0.04, P < 0.001). This is the first study in an Asian population utilizing the revised CPIC-recommended method for translating the CYP2D6 genotype to phenotype. In addition to validating that CYP2D6 genetic variation significantly impacts risperidone metabolism, we demonstrated that revised value for the CYP2D6*10 was superior for genotype to phenotype translation. However, at least for risperidone, subjects with an activity score of 1 presented as phenotypic normal, and not intermediate metabolizers, suggesting that phenotype classification is substrate dependent.

Pharmacogenetics of tamoxifen therapy in Asian populations: from genetic polymorphism to clinical outcomes

BACKGROUND

Compared with western countries, Asian breast cancer patients have unique pathological and biological characteristics. Most of them are premenopausal women with HR positive. Tamoxifen as the first-line drug for premenopausal women with HR+ is involved in multiple enzymes and transporters during metabolizing and transporting process. Variants that cause decreased or inactive gene products leading to abnormal responses in tamoxifen therapy have well been studied in western countries, whereas such information is much less reported in Asian populations.

OBJECTIVE

In order to elucidate the relationship between genetic variants and tamoxifen-induced individual drug reactions in different Asian populations and further identify genotypes/phenotypes with potential therapeutic significance.

METHODS

We reviewed the frequencies of genetic variants in major enzymes and transporter genes involved in the metabolism and transport of tamoxifen across Asian populations as well as significant correlations between genotypes/metabolic phenotypes and metabolites concentrations or BC clinical outcomes.

RESULTS

Significant inter-ethnic differences in allele frequencies was found among Asian populations, such as CYP2D6*4, *10, *41, CYP2C9*2, ABCB1 C3435T and SLCO1B1*5, and CYP2D6*10/*10 is the most common genotype correlated with adverse clinical outcomes. Moreover, we summarized the barriers and controversies of implementing pharmacogenetics in tamoxifen therapy and concluded that more population-specific pharmacogenetic studies are needed in the future.

CONCLUSION

This review revealed more systematic pharmacogenomics of genes involved in the metabolism and transport besides CYP2D6, are required to optimize the genotyping strategies and guide the personalized tamoxifen therapy in Asian populations.

Prevalence of CYP2D6 Genotypes and Predicted Phenotypes in a Cohort of Cambodians at High Risk for Infections with Plasmodium vivax

Clinical failure of primaquine (PQ) has been demonstrated in people with CYP450 2D6 genetic polymorphisms that result in reduced or no enzyme activity. The distribution of CYP2D6 genotypes and predicted phenotypes in the Cambodian population is not well described. Surveys in other Asian countries have shown an approximate 50% prevalence of the reduced activity CYP2D6 allele *10, which could translate into increased risk of PQ radical cure failure and repeated relapses, making interruption of transmission and malaria elimination difficult to achieve. We determined CYP2D6 genotypes from 96 volunteers from Oddor Meanchey Province, Cambodia, an area endemic for Plasmodium vivax. We found a 54.2% frequency of the *10 allele, but in approximately half of our subjects, it was paired with a normal activity allele, either *1 or *2. The prevalence of *5, a null allele, was 9.4%. Overall predicted phenotype percentages were normal metabolizers, 46%; intermediate metabolizers, 52%; and poor metabolizers, 1%.