The mutation polymorphism of G6PD gene coding region in vivax malaria patients in Yunnan Province, China

doi:10.12140/j.issn.1000-7423.2020.02.003

Abstract

Abstract:

Objective To analyze the polymorphism of G6PD (glucose-6-phosphate dehydrogenase) gene mutation in coding regions and its correlations with the proneness to primaquine hemolysis in vivax malaria patients in Yunnan Province, China.Methods Blood samples were collected from the vivax malaria patients receiving an treatment regime of “eight-day therapy of chloroquine/primaquine” in Yunnan Province in 2018. Patients with “tea-colored” urine and reduced G6PD enzyme activity during the primaquine therapy were defined as primaquine-induced acute hemolysis. Human genomic DNA was extracted from the blood samples, and of the G6PD gene fragments containing 12 exon were amplified by PCR and sequenced. The sequences were compared with the nonmutation-type sequence, edited and spliced by DNAStar 11.0 and BioEdit 7.2.5 softwares to generate the cDNA sequence of G6PD gene. MEGA 5.04 and DnaSP 5.10 softwares were used to analyze the mutation polymorphism the selection effect of the cDNAs. The correlation coefficient (r) and odds ratio (OR) for primaquine-induced hemolysis was estimated.Results A total of 184 blood samples were collected, of them 44 generated complete G6PD cDNA strand sequence (length, 1 545 bp), including one sample from a hemolysis case. Alignment of the cDNA chain with the nonmutation-type sequence revealed the mutation sites including c.461 T > A, c.574 C > T, c.786 C > T, c.1059 C > T, c.1311 T > C, and c.1376 G > T, with a frequency of 1.5/100 000, 1.5/100 000, 1.5/100 000, 1.5/100 000, 45.6/100 000 and 1.5/100 000, respectively. The two-site mutation linkage between c.1311 and c.1376 was only present in the genome of the hemolysis cases, and the c.1376 site mutation was positively correlated with the primaquine-induced hemolysis (r = 1.000, P > 0.05). The 44 cDNA strands were defined into 6 haplotypes (Hap_1-Hap_6), with an expected heterozygosity, nucleic acid diversity index π, and Ka/Ks ratio were 0.493, 0.001, and 0.062, respectively. The Hap_2 type containing c.1311 single-site mutation accounted for the highest proportion of (68.2%, 30/44). The site c.1311 with high frequency of mutation had the most complex allele composition, comprising 5 types of zygote including nonmutation hemizygote (9/44, 20.4%), mutant hemizygote (16/44, 36.4%), nonmutation homozygote (5/44, 11.5%), mutant heterozygote (10/44, 22.7%), and mutant homozygote (4/44, 9.1%). The OR value of c.1311 mutation to the primaquine-induced hemolysis was 0.6879 (P > 0.05), showing no correlation in between. The D value of Tjima’s neutral test was -1.414 (P > 0.05), and the D value of the whole fragment was < 1, indicating that the mutations detected in the sequence were not a non-neutral mutation under directional selection pressure; meanwhile the Ka/Ks ratio was 0.062, indicating that the sequence was highly conservative and the frequency of missense mutation was lower than that of synonymous mutation.Conclusion In malaria patients in Yunnan Province, the T > C base substitution at c.1311 is a common mutation in G6PD coding region, but only the G > T mutation at c.1376 is positively correlated with the occurrence of oxidant induced hemolysis.

Key words: Glucose-6-phosphate dehydrogenase gene, Coding region, Mutation, Primaquine, Influence, Correlation

CLC Number:

R531.31

Shu-ping LIU, Ying DONG, Yan-chun XU, Yan LIU, Yan DENG, Cang-lin ZHANG, Meng-ni CHEN. The mutation polymorphism of G6PD gene coding region in vivax malaria patients in Yunnan Province, China[J]. CHINESE JOURNAL OF PARASITOLOGY AND PARASITIC DISEASES, 2020, 38(2): 146-151.

Figures/Tables 2

References 30

[1]	Lander ES, Linton LM, Birren B , et al. Initial sequencing and analysis of the human genome[J]. Nature, 2001,409(6822):860-921.
[2]	Vulliamy T, Mason P, Luzzatto L . The molecular basis of glucose-6-phosphate dehydrogenase deficiency[J]. Trends Genet, 1992,8(4):138-143.
[3]	Alving AS, Carson PE, Flanagan CL . Ickes CE: enzymatic deficiency in primaquine-sensitive erythrocytes[J]. Science, 1956,124(3220):484-485.
[4]	Howes E, Dewi M, Piel FB , et al. Spatial distribution of G6PD deficiency variants across malaria-endemic regions[J]. Malar J, 2013,12:418.
[5]	Lee J, Kim TI, Kang JM , et al. Prevalence of glucose-6-phosphate dehydrogenase (G6PD) deficiency among malaria patients in Upper Myanmar[J]. BMC Infect Dis, 2018,18(1):131.
[6]	Khim N, Benedet C, Kim S , et al. G6PD deficiency in Plasmodium falciparum and Plasmodium vivax malaria-infected Cambodian patients[J]. Malar J, 2013,12:171.
[7]	Zhong Z, Wu H, Li B , et al. Analysis of glucose-6-phosphate dehydrogenase genetic polymorphism in the Hakka population in southern China[J]. Med Sci Monit, 2018,24:7316-7321.
[8]	WHO Malaria Policy Advisory Committee and Secretaria. Malaria policy advisory committee to the WHO: conclusions and recommendations of september 2012 meeting[J]. Malar J, 2012,1:424.
[9]	White NJ, Qiao LG, Qi G , et al. Rationale for recommending a lower dose of primaquine as a Plasmodium falciparum gametocytocide in populations where G6PD deficiency is common[J]. Malar J, 2012,11(1):418.
[10]	Kim S, Nguon C, Guillard B , et al. Performance of the CareStart ^TM G6PD deficiency screening test, a point-of-care diagnostic for pimaquine therapy screening [J]. PLoS One, 2011,6(12):e28357.
[11]	Dong Y, Liu SP, Xu YC , et al. Mutations and predicted structure change of G6PD isolated from a patient with primaquine-induced hemolysis in Yunnan Province[J]. Chin J Parasitol Parasit Dis, 2019,37(4):399-405. (in Chinese)
	( 董莹, 刘淑萍, 徐艳 , 等. 云南省1例伯氨喹诱发溶血患者的G6PD基因编码区突变分析及酶蛋白空间结构预测[J]. 中国寄生虫学与寄生虫病杂志, 2019,37(4):399-405.)
[12]	Rueangweerayut R, Bancone G, Harrell EJ , et al. Hemolytic potential of tafenoquine in female volunteers heterozygous for glucose-6-phosphate dehydrogenase (G6PD) deficiency (G6PD mahidol variant) versus G6PD-normal volunteers[J]. Am J Trop Med Hyg, 2017,97(3):702-711.
[13]	Chu CS, Bancone G, Nosten F , et al. Primaquine-induced haemolysis in females heterozygous for G6PD deficiency[J]. Malar J, 2018,17(1):101.
[14]	Howes RE, Chan ER, Rakotomanga TA , et al. Prevalence and genetic variants of G6PD deficiency among two Malagasy populations living in Plasmodium vivax-endemic areas[J]. Malar J, 2017,16(1):139.
[15]	Minucci A, Moradkhani K, Hwang MJ , et al. Glucose-6-phosphate dehydrogenase (G6PD) mutations database: review of the “old” and update of the new mutations[J]. Blood Cells Mol Dis, 2012,48(3):154-165.
[16]	Howes RE, Piel FB, Patil AP , et al. G6PD deficiency prevalence and estimates of affected populations in malaria endemic countries: a geostatistical model-based map[J]. PLoS medicine, 2012,9(11):e1001339.
[17]	Pamba A, Richardson ND, Carter N , et al. Clinical spectrum and severity of hemolytic anemia in glucose 6-phosphate dehydrogenase-defcient children receiving dapsone[J]. Blood, 2012,120:4123-4133.
[18]	Carter N, Pamba A, Duparc S , et al. Frequency of glucose-6-phosphate dehydrogenase deficiency in malaria patients from six African countries enrolled in two randomized anti-malarial clinical trials[J]. Malar J, 2011,10:241.
[19]	Rigano P, Fabiano C, Pojero F , et al. Glucose 6-phosphate dehydrogenase Palermo R257M: a novel variant associated with chronic non-spherocytic haemolytic anaemia[J]. Br J Haematol, 2010,149(2):296-297.
[20]	Jalloh A, Jalloh M, Gamanga I , et al. G6PD deficiency assessment in Freetown, Sierra Leone, reveals further insight into the molecular heterogeneity of G6PDA[J]. J Hum Genet, 2008,53(7):675-679.
[21]	Nantakomol D, Paul R, Palasuwan A , et al. Evaluation of the phenotypic test and genetic analysis in the detection of glucose-6-phosphate dehydrogenase deficiency Duangdao[J]. Malar J, 2013,12:289.
[22]	He YS, Shi J, Huang YK , et al. The two common G6PD gene mutations in Chinese and their clinical manifestations[J]. Chin J Pediatr, 1999,37(12):733-725. (in Chinese)
	( 何永蜀, 施进, 黄永坤 , 等. 国人常见的两种葡萄糖-6-磷酸脱氢酶基因突变及其临床表现[J]. 中华儿科杂志, 1999,37(12):733-725.)
[23]	Lv RJ, Chen XW, Zhang M , et al. Detection of gene mutation in glucose-6-phosphate dehydrogenase deficiency by RT-PCR sequencing[J]. Chin J Contemp Pediatr, 2016,18(7):630-634. (in Chinese)
	( 吕荣钰, 陈小文, 张民 , 等. RT-PCR测序法检测G6PD缺乏症基因突变[J]. 中国当代儿科杂志, 2016,18(7):630-634.)
[24]	CHEN YS, LI CG, CHEN XW , et al. Study on the enzymatic activity of G6PD G1388A and G1376T common gene in China[J]. China Trop Med, 2007,7(5):666-668. (in Chinese)
	( 陈运生, 李长钢, 陈小文 , 等. 中国人两种常见G6PD基因突变型的酶活性研究[J]. 中国热带医学, 2007,7(5):666-668.)
[25]	Louicharoen C, Patin E, Paul R , et al. Positively selected G6PD-Mahidol mutation reduces Plasmodium vivax density in southeast Asians[J]. Science, 2009,326(5959):1546-1549.
[26]	Guindo A, Fairhurst RM, Doumbo OK , et al. X-linked G6PD deficiency protects hemizygous males but not heterozygous females against severe malaria[J]. PLoS Med, 2007,4(3):e66.
[27]	Goo YK, Ji SY, Shin HI , et al. First evaluation of glucose-6-phosphate dehydrogenase (G6PD) deficiency in vivax malaria endemic regions in the republic of Korea[J]. PLoS One, 2014,9(5):e97390.
[28]	Assefa A, Ali A, Deressa W , et al. Glucose-6-phosphate dehydrogenase (G6PD) deficiency in Ethiopia: absence of common African and Mediterranean allelic variants in a nationwide study[J]. Malar J, 2018,17(1):388.
[29]	Kießling N, Brintrup J, Zeynudin A , et al. Glucose-6-phosphate dehydrogenase activity measured by spectrophotometry and associated genetic variants from the Oromiya zone, Ethiopia[J]. Malar J, 2018,17:358.
[30]	Gómez-Manzo S, Marcial-Quino J, Vanoye-Carlo A , et al. Glucose-6-phosphate dehydrogenase: update and analysis of new mutations around the world[J]. Int J Mol Sci, 2016,17(12):E2069.

位点Loci	单位点突变Singleton variable sites				位点等位Alleles of locus
	碱基替换 Base substitution	密码改变Codon change	变异Variation	频率/10^-5 Frequency/10^-5	半合子Hemizygote				合子 Zygote
					非突变（NM）Nonmutation (NM)		突变（M） Mutation (M)		非突变纯合（NM/NM） NMH^b (NM/NM)		突变杂合（NM/M）NMH^c (NM/M)		突变纯合（M/M）MH^e (M/M)
					数量（构成比/%） No. (CR/%)^a	G6PDd^d	数量（构成比/%） No. (CR/%)^a	G6PDd^d	数量（构成比/%）No. (CR/%)^a	G6PDd^d	数量（构成比/%）No. (CR/%)^a	G6PDd^d	数量（构成比/%）No. (CR/%)^a	G6PDd^d
c.461	T > A	ATT > AAT	I154N	1.5	24 (54.5)	0	1 (2.3)	0	19 (43.2)	0	-	-	-	-
c.574	C > T	CGT > TGT	R192C	1.5	24 (54.5)	0	1 (2.3)	0	19 (43.2)	0	-	-	-	-
c.786	C > T	AAC > ATG	N262N	1.5	24 (54.5)	0	1 (2.3)	0	19 (43.2)	0	-	-	-	-
c.1059	C > T	CAC > AAT	P353P	1.5	25 (56.8)	0	-		18 (40.9)	0	-	-	1 (2.3)	0
c.1311	T > C	TAT > TAC	Y437Y	45.6	9 (20.4)	0	16 (36.4)	1	5 (11.4)	0	10 (22.7)	0	4 (9.1)	0
c.1376	G > T	CGT > CTT	R459L	1.5	24 (54.5)	0	1 (2.3)	1	19 (43.2)	0	-	-	-	-