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

Identification of the effects of pathogenic genetic variations of human CYP2C9 and CYP2D6: an in silico approach

Genetic variations in the human cytochrome P450 family 2 subfamily C member 9 (CYP2C9) and cytochrome P450 family 2 subfamily D member 6 (CYP2D6) genes may affect drug metabolism and lead to alterations in phenotypes. Genetic variations are associated with toxicity, adverse drug reactions, inefficient treatment. Various in silico tools were combined to investigate the deleterious effects of missense non-synonymous single nucleotide polymorphisms (nsSNPs) of the human CYP2C9 and CYP2D6. The structural and functional effects of the high-risk non-synonymous SNPs in the human CYP2C9 and CYP2D6 were predicted by numerous computational mutation analysis methods. Out of 24 pathogenic missense SNPs in the CYP2C9, 22 nsSNPs had a decreasing effect on protein stability and 13 SNPs were showed to be located at conserved positions. Out of 27 high-risk deleterious non-synonymous SNPs in the human CYP2D6, 21 SNPs decreased protein stability and 16 nsSNPs were predicted to be positioned at conserved regions. Our present study suggests that the identified functional SNPs may affect drug metabolism associated with CYP2C9 and CYP2D6 enzymes.

Molecular identification of vivax malaria relapse patients in the Yunnan Province based on homology analysis of the Plasmodium vivax circumsporozoite protein gene

More than 85% of the malaria burden in the Yunnan Province is caused by imported vivax malaria, and Yunnan is also where the majority of vivax malaria patients are diagnosed in China. Timely removal of the infection sources of Plasmodium vivax and its breeding environment remains the key to eliminating the secondary transmission of imported malaria. To that end, blood samples were collected from cases diagnosed and revalidated as single species infection with P. vivax in the Yunnan Province from 2013 to 2020. Specifically, samples from vivax malaria patients with suspected relapses episodes were subjected to PCR amplification, product sequencing, and analysis of the P. vivax circumsporozoite protein (pvcsp) gene. In total, 77 suspected relapse patients were identified out of 2484 cases infected with P. vivax, with a total of 81 recurrent episodes. A total of 156 CDS (coding DNA sequence) chains were obtained through PCR amplification and sequencing of the pvcsp gene from 159 blood samples, 121 of which can be matched to the paired sequences of 59 vivax malaria patients with both primary attack and recurrent experience. Of the 59 pairs of pvcsp gene sequences, every one of 31 pairs showed only one haplotype and no variant sites (VS), meaning every two paired sequence was completely homologous. Every one of the remaining 28 paired sequences had two haplotypes but no length polymorphism, indicating that the paired sequences was "weakly heterologous" with no fragment insertions (or deletions). All 59 vivax malaria patients with recurrences were caused by the activation of P. vivax hypnozoites originated from the same population as the primary infection. The paired analysis of the similarity between high variant genes allowed the identification of relapse episodes caused by P. vivax homologous hypnozoites and also demonstrated pvcsp gene as one of the candidate molecular markers for tracing infection origin.

Effects of CYP2D6 genotypes on Plasmodium vivax recurrence after primaquine treatment: A meta-analysis

OBJECTIVES

To elucidate the relationship between CYP2D6 polymorphisms and Plasmodium vivax recurrence in patients receiving primaquine-based treatment through systematic review and meta-analysis.

METHODS

We searched the PubMed, EMBASE, Cochrane Library, and Web of Science databases for eligible studies published up to August of 2021. We included studies investigating the associations between CYP2D6 polymorphisms and P. vivax recurrence. We evaluated the pooled odds ratio (OR) and 95% confidence interval (CI).

RESULTS

Data from nine studies, including 970 patients, were analyzed. We found that CYP2D6 poor metabolizers (PMs), intermediate metabolizers (IMs), or normal metabolizers slow (NM-Ss) were associated with a 1.8-fold (95% CI, 1.34-2.45; P = 0.0001) higher recurrence of P. vivax than normal metabolizers fast (NM-Fs), extensive metabolizers (EMs), or ultrarapid metabolizer (UMs). Subgroup analysis showed that studies on both Brazilian and Southeast or East Asian individuals had similar results to the main results. Sensitivity analysis by sequentially excluding individual studies also showed robust results (OR range: 1.63-2.01).

CONCLUSIONS

This meta-analysis confirmed that CYP2D6 PMs, IMs, or NM-Ss increased the risk of P. vivax recurrence compared to NM-Fs, EMs, or UMs. The results of this study could be used to predict P. vivax recurrence and suggest CYP2D6 genotype-based primaquine dosing.

Genetic Variation of G6PD and CYP2D6: Clinical Implications on the Use of Primaquine for Elimination of Plasmodium vivax

Primaquine, an 8-aminoquinoline, is the only medication approved by the World Health Organization to treat the hypnozoite stage of Plasmodium vivax and P. ovale malaria. Relapse, triggered by activation of dormant hypnozoites in the liver, can occur weeks to years after primary infection, and provides the predominant source of transmission in endemic settings. Hence, primaquine is essential for individual treatment and P. vivax elimination efforts. However, primaquine use is limited by the risk of life-threatening acute hemolytic anemia in glucose-6-phosphate dehydrogenase (G6PD) deficient individuals. More recently, studies have demonstrated decreased efficacy of primaquine due to cytochrome P450 2D6 (CYP2D6) polymorphisms conferring an impaired metabolizer phenotype. Failure of standard primaquine therapy has occurred in individuals with decreased or absent CYP2D6 activity. Both G6PD and CYP2D6 are highly polymorphic genes, with considerable geographic and interethnic variability, adding complexity to primaquine use. Innovative strategies are required to overcome the dual challenge of G6PD deficiency and impaired primaquine metabolism. Further understanding of the pharmacogenetics of primaquine is key to utilizing its full potential. Accurate CYP2D6 genotype-phenotype translation may optimize primaquine dosing strategies for impaired metabolizers and expand its use in a safe, efficacious manner. At an individual level the current challenges with G6PD diagnostics and CYP2D6 testing limit clinical implementation of pharmacogenetics. However, further characterisation of the overlap and spectrum of G6PD and CYP2D6 activity may optimize primaquine use at a population level and facilitate region-specific dosing strategies for mass drug administration. This precision public health approach merits further investigation for P. vivax elimination.

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.