Patients positive for P. vivax infection treated at the Tropical Medicine Foundation of Amazonas hospital in Manaus, Brazil with a confirmed diagnosis of a monoinfection were enrolled in this study of Duffy/G6PD genotypes and their implications in disease pathogenesis. Patients with comorbidities, hemoglobinopathies, mixed Plasmodium infections, and viral infections were excluded. The World Health Organization (WHO) criteria: uncomplicated and severe were used to classify malaria patients, as described previously (33). There were 225 patients diagnosed with P. vivax malaria included in the study, 52.4% were male.
Genotyping of c.202G > A resulted in 29 (12.9%) positive patients, with 4.4% heterozygous and 0.44% homozygous females, while 8% of males were hemizygous. The c.376A > G was identified in 43 (19.1%) patients, with 18 (8%) heterozygous and 2 (0.88%) homozygous females, while 23 (10.22%) males were hemizygous.
Analysis of the Duffy system in the malaria patients showed a phenotypic distribution of 70 (31.1% - Fy(a + b-), 98 (43.55%) - Fy(a + b+), 56 (24.9%) - Fy(a-b+) and 1 (0.44%) - Fy(a-b-). Table 1 shows the genotypic distribution for the Duffy system, with predominance of the genotype FY*01/FY*02 in the population. These results indicate that the Manaus-Amazonas region has individuals who express three Duffy phenotypes Fy(a + b+), Fy(a + b-) and Fy(a-b+) more frequently with expression of FY*01 and FY*02 antigens.
Table 1
Allele frequency of the Duffy blood group from vivax malaria patients.
Phenotypes | Genotypes | Percentage (%) |
Fy(a + b+) | FY * 01 / FY * 02 | 42,66 |
Fy(a + b-) | FY * 01 / FY * 02N.01 | 16,00 |
FY * 01 / FY * 01 | 15,11 |
Fy(a-b+) | FY * 02 / FY * 02 | 16,00 |
FY * 02 / FY * 02N.01 | 8,88 |
Fy(a + w) | FY*01 / FY * 01W.02 | 0,88 |
Fy(a-b-) | FY * 02N.01 / FY * 02N.01 | 0,44 |
The FY*01/FY*02 genotype was present at higher frequencies in uncomplicated and severe malaria, 45.3% and 39.2%, respectively (Fig. 1). The presentation was more expressive and had a highly severe clinic course when the c.376A > G variant was present. In females, the FY*01/FY*02 allele in combination with the c.376A > G variant was present at a higher frequency in cases of uncomplicated malaria, and in males this combination of variants occurred at a higher frequency in severe malaria cases. Furthermore, when the clinical severity of malaria was compared with G6PD genotypes c.202/c.376A > G which demonstrated a higher frequency of severe malaria when compared with c.202G > A (RR = 5.57, 95% CI: 2.25–17.79, p < 0.001); c.376A > G (RR: 4.49, 95% CI: 1.99–8.43, p < .001) strongly associated with previous malaria episodes (p < 0.001) (Table 2).
Table 2
Genotypic frequency of G6PD variants among patients with uncomplicated and severe vivax malaria.
G6PD Genotypes | MALARIA N (%) | RR (CI) | p-value | PREVIOUS MALARIA EPISODES N (%) | RR (CI) | p-value |
Severe | Uncomplicated | Yes | No |
202GA/AA | 22 (26.2) | 07 (6.0) | 5.57 (2.25–17.79) | < .001* | 17 (25.8) | 08 (6.3) | 5.11 (2.02–12.63) | < .001* |
202GG | 62 (73.8) | 110 (94,0) | 49 (74.2) | 118 (93.7) |
376AG/GG | 30 (32.6) | 13 (10,6) | 4.09 (1.99–8.43) | < .001 | 25 (33.8) | 17 (12.6) | 3.54 (1.75–7.13) | < .001* |
376AA | 62 (67.4) | 110 (89.4) | 49 (66.2) | 118 (87.4) |
202GA/AA / 376AG/GG | 17 (21,5) | 03 (2,7) | 10.54 (2.83–35.66) | < .001** | 13 (21.0) | 07 (5.6) | 4.47 (1.68–11.88) | < .001* |
202GG / 376AA | 62 (78,5) | 110 (97,3) | 49 (79.0) | 118 (94.4) |
RR: Relative Risk (CI): Interval of Confidence * χ2 test (Yates’s corrected) * * Fisher’s exact test N: cases |
For c.202G > A, the highest frequency among females was the Fy(a + b+) and Fy(a + b-) phenotype in severe malaria cases and it was not found at the same frequency in the Fy(a-b+) phenotype, which was frequent among cases of uncomplicated malaria. In males, all Duffy phenotypes were more frequent in severe malaria, with Fy(a + b-) representing the largest number of cases. For c.376A > G, we observed a higher number of cases with the polymorphism for the Fy(a + b) phenotype in patients with severe malaria, with 47.36% of cases reported in males (data not showed).
The only allele FY*02N.01/FY*02N.01 responsible for the Fy(a-b-) phenotype, was found in a single female patient with uncomplicated malaria, without the c.202G > A and c.376A > G variants. For uncomplicated malaria, the frequency of c.202G > A and c.376A > G variants decreased, particularly when these polymorphisms were concomitant. Severe malaria showed no variations in frequency for the c.202G > A variant, however, in the presence of A376G there was a slight decrease in frequency. When combined with the GATA variant, the c.202G > A and c.376A > G variants were observed at a lower frequency in uncomplicated malaria, but this finding did not reach statistical significance, despite a slight decrease in frequency. For severe malaria, the presence of the combined polymorphisms showed a moderate variation in frequency (Fig. 2).
Results showed no statistically significant association for women with P. vivax infection, nor with the frequency of mutated Duffy alleles. The authors do not speculate about the race of the participants due to the strong regional ancestral mix found in the Amazonian Caboclos, which originate with the arrival of Caucasians and Blacks to indigenous lands (Data not showed).
There was no significant association of parasite count and density with the Duffy antigen genotypes and phenotypes. However, However, Duffy phenotype Fy(a-b+) (p = 0.003) and genotype FY*02/FY*02 (p = 0.007) presented the highest values parasitemia density of the vivax malaria (Figs. 3, 4).
In this study, the highest frequency of the FY*01/FY*02 genotype and the Fy(a + b+) and Fy(a + b-) phenotypes was found in malaria cases. These findings corroborate studies reported by Cavasini et al. (2007), which correlated the high frequency of the FY*01 and FY*02 alleles among P. vivax malaria patients. They concluded that the FY*01/FY*02 and FY*01/FY*01 genotypes were associated with a high frequency of Plasmodium vivax infection, and suggested that these individuals have a higher risk of disease (34).
The FY*01/FY*02 and FY*01/FY*01 genotypes are associated with increased frequency of P. vivax infection, while FY*01/ FY*02W.01 and FY*02/ FY*02W.01 were associated with low parasitic density levels (21).
The presence of a single case with null Duffy variants in this study supports a possible advantageous selection, driven by defense mechanisms against P. vivax. in endemic areas, which is corroborated by literature that has reported the Duffy negative blood group as a protective factor against P. vivax malaria infection (35).
In a study conducted in São Paulo, the phenotypic frequencies found in blood donors for Duffy blood system antigens were 19.8% for the phenotype Fy(a + b-) in Caucasians and 14.0% in Blacks, Fy(a + b+) in 41.4% of Caucasians and 1.6% of Blacks, Fy(a-b+) in 37.8% of Caucasians and 17.5% of Blacks and Fy(a-b-) in 1.1% of Caucasians and 66.9% of Blacks (36).
The results and clinical presentation of patients in this study support the hypothesis that G6PD deficiency does not confer a decreased risk of severe malaria infection (Table 3). Many controversies remain as hemizygous males and heterozygous females failed to have changes in frequency of uncomplicated malaria. However, the results showed that both male or female sex is independent and provides no protection against severe malaria.
Table 3
Genotypic frequency distribution of G6PD variants among uncomplicated and severe vivax malaria patients by the presence of Duffy GATA normal or mutated variants.
Malaria | Duffy GATA | c.202 G > A N (%) | p-value | c.376 A > G N (%) | p-value | c.202 G > A - c.376 A > G N (%) | p-value |
GG | GA/AA | AA | AG/GG | GG/AA | GA/AA - AG/GG |
Uncomplicated | Normal | 86 (94.5) | 05 (5.5) | .487** | 86 (88.7) | 11 (11.3) | .452** | 86 (97.7) | 02 (2.3) | .531** |
Mutated | 24 (92.3) | 02 (7.7) | 24 (92.3) | 02 (7.7) | 24 (96.0) | 01 (4.0) |
Severe | Normal | 45 (77.6) | 13 (22.4) | .181* | 45 (70.3) | 19 (29.7) | .252* | 45 (83.3) | 09 (16.7) | .107* |
Mutated | 17 (65.4) | 09 (34.6) | 17 (60.7) | 11 (39.3) | 17 (68.0) | 08 (32.0) |
Total | Normal | 131 (87.9) | 18 (12.1) | .087* | 131 (81.4) | 30 (18.6) | .249* | 131 (92.3) | 11 (7.7) | .043* |
Mutated | 41 (78.8) | 11 (21.2) | 41 (75.9) | 13 (24.1) | 41 (82.0) | 09 (18.0) |
* χ2 test (Yates’s corrected) * * Fisher’s exact test N: cases |
For the process of parasite invasion into red blood cells, Duffy protein is functionally important, and it has been shown that in regions where the malaria transmission rates are average, as in the inhabitants of the Amazon, this protein is naturally immunogenic (37).
In a study with P. vivax malaria patients the FY*01 and FY*02 alleles were found to have low, medium, and high parasitic density, but in the presence of the GATA variant, genotypes with alleles FY*02N.01 and FY*02W.01, were found only in patients with low parasitic density and low symptomatology (21).
A study performed in the state of Pará with a population of African descent demonstrated the presence of the c.202G > A variant was 0.060, the Duffy blood group included 24.3% Duffy negative and 41.3% individuals heterozygous for FY*02w. The frequency of the FY*02w allele was 41.0%. These findings support the monitoring of individuals with G6PD deficiency for use of primaquine during the routine care of Afro-descendant communities of the Trombetas, Erepecuru, and Cumná rivers, to assess risks of hemolytic crisis in recurrent cases of malaria in the region (38).
In recent decades, studies have revealed the complexity of Duffy system phenotypic and genotypic variation and variants of G6PD (A-) have a significant impact on the distribution of human populations in areas where malaria is endemic. The Duffy system and G6PD are, therefore, polymorphic systems that offer great challenges to researchers not only due to their academic importance, but also given their potential applications in the treatment of vivax malaria (39). That is, where natural selection occurs within a population endemic for malaria, natural adaptations may result from genetic variation that provides a partial defense mechanism against P. vivax infections (40).