Review | Open | Published:
Genetic diversity of Plasmodium vivax isolates from Azerbaijan
Malaria Journalvolume 3, Article number: 40 (2004)
Plasmodium vivax, although causing a less serious disease than Plasmodium falciparum, is the most widespread of the four human malarial species. Further to the recent recrudescence of P. vivax cases in the Newly Independent States (NIS) of central Asia, a survey on the genetic diversity and dissemination in Azerbaijan was undertaken. Azerbaijan is at the crossroads of Asia and, as such, could see a rise in the number of cases, although an effective malaria control programme has been established in the country.
Thirty-six P. vivax isolates from Central Azerbaijan were characterized by analysing the genetic polymorphism of the circumsporozoite protein (CSP) and the merozoite surface protein 1 (MSP-1) genes, using PCR amplifications and amplicons sequencing.
Analysis of CSP sequences showed that all the processed isolates belong to the VK 210 type, with variations in the alternation of alanine residue (A) or aspartic acid residue (D) in the repeat motif GDRA(A/D)GQPA along the sequence. As far as MSP-1 genotyping is concerned, it was found that the majority of isolates analysed belong to Belem and Sal I types. Five recombinant isolates were also identified. Combined analysis with the two genetic markers allowed the identification of 19 plasmodial sub-types.
The results obtained in the present study indicate that there are several P. vivax clones circulating in Azerbaijan and, consequently, a careful malaria surveillance could be of paramount importance to identify, at early stage, the occurrence of possible P. vivax malaria outbreaks.
Plasmodium vivax is the most widely distributed human parasite, with an estimated burden of 70–80 million cases annually . In some parts of the world (Asia, South America), it is the most prevalent form of the four human malarial parasites. Although it causes a less severe disease than Plasmodium falciparum, being rarely lethal, P. vivax affects the working capacity of the population and the lack of efficient drug distribution favors the onset of drug resistant strains [2, 3]. Imported malaria is an increasing health problem in Western Europe, where about 6,500 cases are reported annually in Germany, France, Italy and the United Kingdom . Although P. falciparum infections account for the majority of cases (64%), P. vivax is responsible for an additional 23% . Presence in this area of residual anopheline populations susceptible to P. vivax infection represents a permanent risk for the occurrence of P. vivax indigenous malaria cases, as recently occurred in central Italy [5, 6]. Since 1970, malaria had been eradicated in central Asia, except for some residual foci in two countries belonging to the Newly Independent States (NIS), i. e. Azerbaijan and Tajikistan (WHO, Regional Office for Europe, unpublished document). At the beginning of the 1990s, the situation changed dramatically due to the re-emergence of malaria in the NIS area and especially in Tajikistan, where at present an epidemic is still in progress [7, 8]. In these countries, the existing state of the primary health care system is extremely precarious, especially in rural areas and in small villages. Malaria is a common disease, which can easily re-establish itself when a lack of control occurs.
In comparison with P. falciparum, molecular studies of the genetic diversity and dissemination of P. vivax are scanty. Recently, 33 polymorphic tandem repeats (TRs) of P. vivax and a P. vivax polymorphic microsatellite have been identified and shown to be useful in population studies [9, 10]. The merozoite surface protein 3α (MSP3-α) gene also seems to be a good candidate for studying the genetic diversity of P. vivax populations, since PCR-RFLP products indicate the presence of up to 13 alleles [11, 12]. However, the circumsporozoite protein (CSP) and merozoite surface protein 1 (MSP-1) genes still remain the most studied molecular markers in genetic epidemiological surveys carried out in P. vivax endemic areas.
In the frame of a malaria research project funded by the European Commission, a molecular study was undertaken in Azerbaijan, aimed at collecting information on the genetic make-up of P. vivax natural populations present in this endemic country. For this purpose the extent of polymorphism of CSP and MSP-1 genes were analysed in parasite isolates from five localities of central Azerbaijan by using PCR amplification and sequencing.
Materials and Methods
Study area and samples collection
Azerbaijan covers an area of 29.540 Km2, with a populations of approximatively 2,5 millions. The climate is typically continental with an average temperature between 12 and 15°C and a rainfall between 200 and 600 mm per year. Climatic and agro-ecological conditions of this area make the environment favourable to mosquito vectors breeding. The major malaria vector is Anopheles sacharovi that breeds preferably in lakes, swamps, irrigation canals and pools. Although it prefers well-oxygenated water, it is known to tolerate moderate salty water. Other anopheline species found in this area are A. maculipennis, A. subalpinus, A. superpictus and A. hyrcanus, all of which are considered secondary vectors of malaria transmission . Malaria transmission occurs in Azerbaijan mainly from June to October. In the last years, number of malaria cases showed a negative trend, accounting for 610 cases in 2000, 418 cases in 2001, 203 cases in 2002.
During summer 2002, a malaria epidemiological survey was performed in central Azerbaijan, in the frame of a Malaria Surveillance Programme launched in year 2001 by the Ministry of Health in collaboration with WHO. Active case detection was carried out in five districts included in previously identified sentinel sites, namely Mingaçevir, Beylagan, Imisli, Saatli and Sabirabad. A map of Azerbaijan with the study area is shown in Figure 1.
All individuals who visited the district health centres or were found in villages with a history of recent fever and no history of travel in the past few months were considered. In this context, a total of 36 infected individuals with positives blood smears at the microscopic examination, collected between August and September 2002, were selected for the genetic study. The age of patients ranged from 6 to 78 years and parasitaemia varied from 288 to 12,800 parasites/μl. For the molecular analysis, a blood sample of about 1 ml was taken from by venipuncture before drug treatment was given. Patients or the guardians of children were informed about the study. According to the international rules for research involving human subjects, any information which would identify a participant was removed in order to keep each sample processed anonymous. Number of samples for each district is shown in Table 1.
Plasmodial DNA was extracted from 200 μl of each infected blood sample using QIAamp DNA blood kit following the manufacturer's instructions (Qiagen, CA).
Circumsporozoite (CSP) marker analysis
The CSP gene was amplified for the most part of samples using PV5 and PV6 primers . Samples that did not provide good PCR products with this set of primers were processed a second time by using CSP-A2  as forward primer and PV6bis (5'-CACAGGTTACACTGCATGGAGT-3') as original reverse primer. PCR amplification was performed in a reaction mixture of 50 μl containing the parasite DNA, 1x reaction buffer, 2.5 mM MgC12, 80 μM of each deoxynucleotide triphosphate, 6 pmol of each primer and 1.3 U of Taq polymerase (Promega, Madison, USA). The PCR programme was: denaturation at 94°C for five minutes; 34 cycles of one minute at 94°C, one minute at 54°C and two minutes at 72°C. The PCR products were separated using electrophoresis on a 1.5 % NuSieve gel and the band of interest was cut out and purified using the QIAquick PCR purification Kit (Qiagen). The purified product was sequenced in both directions using an ABI-PRISM 373 sequencer. Nucleotide or amino acid sequences were aligned first using the CLUSTAL X programme  with manual editing and adjustments made using the MUST package . The ExPASy Molecular Biology Server http://us.expasy.org/ was used to convert nucleotide sequences into amino acid sequences. The GenBank accession numbers of the eight sub-types of VK210 type are from AY792359 to AY 792366.
Merozoite surface protein 1 (MSP-1) marker analyses
A portion of the MSP-1 gene (the region encompassing the interspecies conserved blocks ICB5 and ICB6) was amplified using a nested PCR with, respectively, the two outer primers A5 and A6  and the two inner primers MSP1N1 forward and MSP1N2 reverse . The first round of amplification was performed in a reaction mixture of 50 μl containing parasite DNA, 1x reaction buffer, 2.5 mM MgC12, 200 μM of each deoxynucleotide triphosphate, 30 pmol of each primer and 2.5 U of Taq polymerase (Promega, Madison, USA). For the second round, 10 μl of the first amplification product was added to a fresh PCR mixture with 30 pmoles of each inner primer. The thermal profile was: denaturation at 94°C for five minutes; 35 cycles of 94°C for one minute, 60°C for one minute and 75°C for three minutes. All nested-PCR products were purified by Microcon-PCR (Millipore), following the manufacturer's instructions and sequenced in both directions at the MWG Biotech. The results were analysed by means of Omiga 2.0 (ACCELRYS, Cambrige) and Mega 2 (S. Kumar, K. Tamura, M. Nei and Pennsylvania State University) computer programmes. The GenBank accession numbers of the 36 nucleotide sequences from P. vivax isolates are from AY789657 to AY789692.
The aligned nucleotide sequences of CSP were converted to a distance matrix (% of differences) using the Net algorithm of the MUST package . The dendrogram was generated using the neighbour-joining method . Bootstrap proportions were used to assess the robustness of the tree with 1,000 bootstrap replications .
CSP sequences obtained from 36 Azerbaijan P. vivax isolates were found to belong to the VK210 type . The isolates tested displayed variations in the peptide repeat motifs GDRA(A/D)GQPA with different alternations of non-synonymous codons GCT or GAT, respectively, coding for alanine (A) and aspartic acid (D) (Figure 2). All our sequence types had the same three repeat units (GDRAAGQPA) at the 3' end, identical to that of the VK210 type. Furthermore four non-synonymous mutations were found, one being the RDRADGQPA variant (sequence named in the present study as sub-type 1), already described in North Korean and Chinese isolates . In summary, eight different sub-types of VK 210 were observed (Figure 2 and Figure 4). Among all 36 azeri isolates analysed, 24 isolates were found to have identical sequence (sub-type 4, Table 1 and Figure 4). In particular, the Beylagan (n = 5) and Mingaçevir (n = 7) isolates appeared the most diversified since they displayed four and five different sub-types respectively. The Imishli (n = 9), Saatli (n = 5) and Sabirabad (n = 10) isolates only showed two different sub-types each-one. Figure 4 clearly shows that the genetic diversity of CSP is relatively small inside the Azerbaijan isolates when compared to the South Korean and Chinese isolates.
The majority of Azerbaijan isolates (Table 1) belong to either the Belem (22 isolates, all with the same poly-Q region of 21 repeats) or the Sal I (9 isolates) types already described . Only 5 P. vivax isolates were identified as recombinant types. Isolate Satl1 (sequence named in the present study as sub-type T) could be ascribed to the type 3a (accession no. D85252, ), while isolates Bey 14 and Min7 (sub-type S) and isolates Sat5 and Sat7 (sub-type U) seem to be the result of recombinant events between the recombinant type 3a and Sal I. All the three recombinant types showed a different number of poly-Q repeats (Table 1). In addition to these sources of diversity, nucleotide substitutions could be observed, allowing the identification of 17 sub-types (Table 1 and Figure 3). The Imisli and Sabirabad districts appeared to be less diversified, accounting for five different genotypes for 9 isolates and six genotypes for 10 isolates, respectively. Finally, Saatli district was found to have the greatest variability, with four different genotypes for 5 isolates.
Combined analysis between the two markers
By combining the results of genotyping obtained by CSP and MSP-1, 19 P. vivax sub-types (Figure 5) were identified as circulating in the central region of Azerbaijan. The sub-type named G/4 with the greatest representation (n = 14 isolates), was detected in all districts investigated. Genotypes identified as M/4 and U/5 were observed twice in the districts of Sabirabad and Saatli, respectively. Genotype F/4 was detected once in Imisli and Saatli districts, as was for genotype S/2, detected once in Beylagan and Mingacevir districts.
Discussion and conclusions
For CSP, the main variations already reported in the literature consist of two variant sequences, VK210 and VK247 that show a variable number of repeat units, GDRA(A/D)GQPAA and ANGAGNQPG, respectively, with some variant positions within the sequence [22, 23]. These two variants have a worldwide distribution, and locally their distribution have been also correlated with climatic gradients or with the Anophelines vector specificity . Studies carried out in some Asian endemic countries, i.e. South and North Korea, China, the Philippines, the Solomon Islands and Thailand [27, 15, 28, 12] suggest that CSP has a limited value as a molecular marker for genetic variability when used alone. In the present study, all the analysed P. vivax isolates from Azerbaijan were found to belong to the type VK210 and they constitute a group of eight CSP sub-types that closely linked in the dendrogramme shown in Figure 4. Differently from what reported by other Authors  who showed that the CSP sequence analysis allows detecting the geographic origin of plasmodial isolates, our results did not support its use for tracking the geographic origin of Azeri isolates since we dealt with a limited number of samples studied. No genotype association with particular sampling districts was observed since, for example, the most common sub-type identified (G/4) was present in all five districts.
Our results confirmed that MSP-1 is a good polymorphic marker. In particular, the region of the gene known to be highly polymorphic and discriminative between the Belem and Sal I types was analysed . A variable poly-glutamine (poly-Q) region is characteristic of the Belem type and represents the principal source of genetic diversity of this marker. Moreover, a poli-Q region is the recombination site between the two types Belem and Sal I and, as shown in the literature, interallelic recombinations between the two types are frequent . The total genetic diversity observed when including the nucleotide substitutions is relatively important taking into account the low endemicity of studied area. Similar results were observed in Southeastern Iran and in Thailand [29, 12], low endemic countries for P. vivax malaria as well, where the authors detected the two types Belem and Sal I, together with several recombinant types. In particular, in the study carried out in Iran by Zakeri et al., the analysis of MSP1 genetic diversity on a total of 16 plasmodial isolates leaded to the identification of 14 genetic sub-types. It is worth noting that such a high level of diversity is probably due to the small sample size.
Our results show 17 genetic sub-types on a total of 36 isolates analysed and a quite high MSP-1 polymorphism also in Azerbaijan. As suggested in other studies  and also reported by Zakeri et al., it is possible to speculate that the observed genetic diversity could be also explained considering the studied area, i.e. central Azerbaijan, as transit road of the country and also of neighboring Asian country, where the circulating P. vivax populations show considerable MSP-1 genetic diversity. However, further studies aimed at collecting more information about people moving within the whole country and to closer countries are needed to verify the above hypothesis. The combined analysis of CSP and MSP1 sequence polymorphism has led to the identification of a total of 19 P. vivax sub-types, confirming that the simultaneous use of more than one genetic marker in this kind of study enhances the knowledge of genetic diversity existing in the parasite populations. The results of the current study show the circulation of multiple plasmodial clones in the studied area thus leading to the conclusion that malaria surveillance activities must be maintained in Azerbaijan in order to avoid serious disease outbreaks in the future.
The understanding of the polymorphism extent in surface antigens as CSP and MSP-1 and the resulting genetic diversity in P. vivax field populations could help in implementing malaria control activities being a crucial step for the development of a malaria vaccine.
Sina B: Focus on Plasmodium vivax. Trends Parasitol. 2002, 18: 287-289. 10.1016/S1471-4922(02)02329-2.
Wilairatana P, Silachamroon U, Krudsood S, Singhasivanon P, Treeprasertsuk S, Bussaratid S, Phumratanaprapin W, Srivilirit S, Looareesuwan S: Efficacy of primaquine regimens for primaquine-resistant Plasmodium vivax in Thailand. Am J Trop Med Hyg. 1999, 61: 973-977.
Imwong M, Pukrittakayamee S, Looareesuwan S, Pasvol G, Poirreiz J, White NJ, Snounou G: Association of genetic mutations in Plasmodium vivax dhfr with resistance to sulfadoxine-pyrimethamine: geographical and clinical correlates. Antimicrob Agents Chemother. 2001, 45: 3122-3127. 10.1128/AAC.45.11.3122-3127.2001.
Anonymous: Euro Surveillance. 1998, 3: 35-42.
Baldari M, Tamburro A, Sabatinelli G, Romi R, Severini C, Cuccagna G, Fiorilli G, Allegri MP, Buriani C, Toti M: Malaria in Maremma, Italy. Lancet. 1998, 351: 1246-1247. 10.1016/S0140-6736(97)10312-9.
Severini C, Menegon M, Gradoni L, Majori G: Use of Plasmodium vivax merozoite surface protein 1 gene sequence analysis in the investigation of an introduced malaria case in Italy. Acta Trop. 2002, 84: 151-157. 10.1016/S0001-706X(02)00186-9.
Sabatinelli G, Ejov M, Joergensen P: Malaria in the WHO European Region (1971–1999). Euro Surveillance. 2001, 6: 61-65.
Sabatinelli G, Aliev S, Ejov M, Verbitsky V: The malaria epidemiological situation in Tajikistan. Parassitologia. 2000, 42 (suppl 1): 148-
Feng X, Carlton JM, Joy DA, Mu J, Furuya T, Suh BB, Wang Y, Barnwell JW, Su XZ: Single-nucleotide polymorphisms and genome diversity in Plasmodium vivax. Proc Natl Acad Sci USA. 2003, 100: 8502-8507. 10.1073/pnas.1232502100.
Gomez JC, McNamara DT, Bockarie MJ, Baird JK, Carlton JM, Zimmerman PA: Identification of a polymorphic Plasmodium vivax microsatellite marker. Am J Trop Med Hyg. 2003, 69: 377-379.
Bruce MC, Galinski MR, Barnwell JW, Snounou G, Day KP: Polymorphism at the merozoite protein -3α locus of Plasmodium vivax: global and local diversity. Am J Trop Med Hyg. 1999, 61: 518-525.
Cui L, Mascorro CN, Fan Q, Rzomp KA, Khuntirat B, Zhou G., Chen H, Yan G, Sattabongkot J: Genetic diversity and multiple infections of Plasmodium vivax malaria in Western Thailand. Am J Trop Med Hyg. 2003, 68: 613-619.
Temel T: Malaria from the gap: need for cross-sector co-operation in Azerbaijan. Acta Trop. 2004, 89: 249-259. 10.1016/j.actatropica.2003.08.005.
Kain KC, Brown AE, Webster HK, Wirtz RA, Keystone JS, Rodriguez MH, Kinahan J, Rowland M, Lanar DE: Circumsporozoite genotyping of global isolates of Plasmodium vivax from dried blood specimens. J Clin Microbiol. 1992, 30: 1863-1866.
Kim T, Kim YJ, Song KJ, Song JW, Cha SH, Kim YK, Shin YK, Suh IB, Lim CS: The molecular characteristics of circumsporozoite protein gene subtypes from Plasmodium vivax isolates in Republic of Korea. Parasitol Res. 2002, 88: 1051-1054. 10.1007/s00436-002-0699-z.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The clustal X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucl Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.
Philippe H: MUST: a computer package of management utilities for sequences and trees. Nucl Acids Res. 1993, 21: 5264-5272.
Premawansa S, Snewin VA, Khouri E, Mendis KN, David PH: Plasmodium vivax: recombination between potential allelic types of the merozoite surface protein MSP1 in parasites isolated from patients. Exp Parasitol. 1993, 76: 192-199. 10.1006/expr.1993.1022.
Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.
Felsenstein J: Confidence limits on phylogenetics: an approach using the bootstrap. Evol. 1985, 39: 783-791.
Cavalli-Sforza LL, Edwards AWF: Phylogenetic analysis: model and estimation procedures. Am J Hum Genet. 1967, 19: 233-257.
Rosenberg R, Wirtz RA, Lanar DE, Sattabongkot J, Hall T, Waters AP, Prasittisuk C: Circumsporozoite protein heterogeneity in the human malaria parasite Plasmodium vivax. Science. 1989, 245: 973-976.
Mann VH, Good MF, Saul A: Diversity in the circumsporozoite protein of Plasmodium vivax: Does it matter?. Parasitol Today. 1995, 1: 33-36. 10.1016/0169-4758(95)80107-3.
Del Portillo HA, longacre S, Khouri E, David P: Primary structure of the merozoite surface antigen 1 of Plasmodium vivax reveals sequences conserved between different Plasmodium species. Proc Natl Acad Sci USA. 1991, 88: 4030-4034.
Putaporntip C, Jongwutiwes S, Tanabe K, Thaithong S: Interallelic recombination in the merozoite surface protein 1 (MSP1) gene of Plasmodium vivax from Thai isolates. Mol Bioch Parasitol. 1997, 84: 49-56. 10.1016/S0166-6851(96)02786-7.
Rodriguez MH, Gonzalez-Ceron L, Hernandez JE, Nettel JA, Villarreal C, Kain KC, Wirtz RA: Different prevalences of Plasmodium vivax phenotypes VK210 and VK247 associated with the distribution of Anopheles albimanus and Anopheles pseudopunctipennis in Mexico. Am J Trop Med Hyg. 2000, 62: 122-127.
Lim CS, Kim YK, Lee KN, Kim SH, Hoffman KJ, Song KJ, Song JW: The analysis of Circumsporozoite protein gene sequences from South Korean isolates of Plasmodium vivax. Ann Trop Med Parasitol. 2001, 95: 229-235. 10.1080/00034980120053997.
Mann VH, Huang T, Cheng Q, Saul A: Sequence variation in the circumsporozoite- protein gene of Plasmodium vivax appears to be regionally biased. Mol Bioch Parasitol. 1994, 68: 45-52. 10.1016/0166-6851(94)00148-0.
Zakeri S, Djadid ND, Zeinali S: Sequence heterogeneity of the merozoite surface protein-1 gene (MSP-1) of Plasmodium vivax wild isolates in Southeastern Iran. Acta Trop. 2003, 88: 91-97. 10.1016/S0001-706X(03)00192-X.
We are grateful to N. G. Eyvazov (MOH/RCHE, Baku) and F. Abdullayev (WHO/EURO, Baku) for their valuable support to our research activities in Azerbaijan. We thank F. Severini (ISS, Rome) for the figures editing and N. Billote (CIRAD, Montpellier) for advises improving the manuscript.
The work was supported by a grant from the European Commission, INCO-Copernicus 2 project contract ICA2-CT-2000-10046 (acronym VIVAXNIS).
S. Mammadov, N. Aliyev, E. Gasimov were involved in field collection of blood samples and microscopy examinations. M.C. Leclerc and A. Cligny did the CSP sequence analysis and M. Menegon did the MSP-1 sequence analysis. M.C. Leclerc and J.L. Noyer did the distance analyses. M.C. Leclerc wrote the report with major contributions of C Severini, M Menegon and G. Majori. G. Majori coordinated the field activities carried out in Azerbaijan. C. Severini, as scientific coordinator of the VIVAXNIS project mentioned below, got the financial support.