Temporal changes in genetic diversity of msp-1, msp-2, and msp-3 in Plasmodium falciparum isolates from Grande Comore Island after introduction of ACT
Malaria Journal volume 17, Article number: 83 (2018)
Malaria is still one of the serious public health problems in Grande Comore Island, although the number of annual cases has been greatly reduced in recent years. A better understanding of malaria parasite population diversity and transmission dynamics is critical for assessing the effectiveness of malaria control measures. The objective of this study is to investigate temporal changes in genetic diversity of Plasmodium falciparum populations and multiplicity of infection (MOI) in Grande Comore 10 years after introduction of ACT.
A total of 232 P. falciparum clinical isolates were collected from the Grande Comore Island during two sampling periods (118 for 2006‒2007 group, and 114 for 2013‒2016 group). Parasite isolates were characterized for genetic diversity and complexity of infection by genotyping polymorphic regions in merozoite surface protein gene 1 (msp-1), msp-2, and msp-3 using nested PCR and DNA sequencing.
Three msp-1 alleles (K1, MAD20, and RO33), two msp-2 alleles (FC27 and 3D7), and two msp-3 alleles (K1 and 3D7) were detected in parasites of both sampling periods. The RO33 allele of msp-1 (84.8%), 3D7 allele of msp-2 (90.8%), and K1 allele of msp-3 (66.7%) were the predominant allelic types in isolates from 2006–2007 group. In contrast, the RO33 allele of msp-1 (63.4%), FC27 allele of msp-2 (91.1%), and 3D7 allele of msp-3 (53.5%) were the most prevalent among isolates from the 2013–2016 group. Compared with the 2006‒2007 group, polyclonal infection rates of msp-1 (from 76.7 to 29.1%, P < 0.01) and msp-2 (from 62.4 to 28.3%, P < 0.01) allelic types were significantly decreased in those from 2013‒2016 group. Similarly, the MOIs for both msp-1 and msp-2 were higher in P. falciparum isolates in the 2006–2007 group than those in 2013–2016 group (MOI = 3.11 vs 1.63 for msp-1; MOI = 2.75 vs 1.35 for msp-2). DNA sequencing analyses also revealed reduced numbers of distinct sequence variants in the three genes from 2006‒2007 to 2013‒2016: msp-1, from 32 to 23 (about 28% decline); msp-2 from 29 to 21 (about 28% decline), and msp-3 from 11 to 3 (about 72% decline).
The present data showed dramatic reduction in genetic diversity and MOI among Grande Comore P. falciparum populations over the course of the study, suggesting a trend of decreasing malaria transmission intensity and genetic diversity in Grande Comore Island. These data provide valuable information for surveillance of P. falciparum infection and for assessing the appropriateness of the current malarial control strategies in the endemic area.
Malaria is a major infectious disease that led to ~ 212 million clinical cases and about 429,000 deaths worldwide in 2016 . Plasmodium falciparum malaria had been widely distributed throughout the Union of Comoros (Grande Comore, Moheli, and Anjouan Islands) and posed a serious impediment to socioeconomic development historically . To effectively control malaria in Comoros, many malaria control measures have been deployed since 2000s, including indoor residual sprayings (IRS), long-lasting insecticide nets (LLINs), artemisinin-based combination therapy (ACT), intermittent presumptive treatment (IPT) for all pregnant women, and, particularly, mass drug administration (MDA) of ACT. These malaria control measures have resulted in substantial decrease malaria infection, from 108,260 cases in 2006 to 1072 in 2015 (about 99.0% decline) in Comoros, with no malaria-related deaths. However, despite the great efforts in malaria control, the annual malaria cases increased from 2015 (1072 cases) to 2016 (1372 cases) in Comoros, and the threat of future malaria outbreak remains. Furthermore, malaria transmission intensity differs among the three islands of Comoros (Grande Comore, Moheli, and Anjouan Islands). In Anjouan and Moheli, there was a limited numbers of malaria annual cases during 2014 to 2016 (7 and 5 in 2014; 3 and 8 in 2015; 4 and 6 in 2016, respectively) without local malaria infection; in contrast, the Grande Comore accounted for about 99% of the total of malaria annual cases reported in Comoros during 2013–2016 (e.g. 53,979 in 2013; 2130 in 2014; 1061 in 2015; 1362 in 2016) due to low coverage level of ACT-based MDA. To achieve an ambitious goal of completely eliminating malaria by 2020 in Comoros, there is an urgent need to develop effective and affordable malaria control and treatment strategies.
To date, several malaria pre-erythrocytic (RTS/S and PfSPZ Vaccine) or erythrocytic (MSP-1, MSP-2, and MSP-3) stage vaccines have been designed to induce immunity against the pre-erythrocytic or erythrocytic stage of the malaria parasites, respectively [3, 4]. Although several vaccines are now being tested in clinical Phase I and II trials (MSP-1, MSP-2, and MSP-3) or even have completed the pivotal Phase III clinical testing (RTS/S), the efficacies of these vaccines have been low, with limited impact against clinical malaria [5, 6]. One of the difficulties in developing an effective vaccine against P. falciparum parasite is the extensive genetic diversity of vaccine targets allowing parasites with mutated genes to escape from the host’s immune response [7, 8]. Thus, studying genetic diversity of malaria parasites in endemic areas may provide important information to improve vaccine design. Additionally, the genetic diversity of P. falciparum parasites has been widely used as an indicator of level of malaria transmission intensity in endemic regions, thus serving as a tool to evaluate the effectiveness of malaria control and intervention.
Polymorphic genetic marks, such as microsatellites and genes encoding merozoite surface proteins (msp-1, msp-2, and msp-3) have been widely used for characterization of parasite genetic diversity [9, 10]. Currently, only one study described the genetic structure of P. falciparum parasites collected from Comoros Archipelago (Grande Comore, Moheli, Anjouan, and Mayotte) using microsatellite loci , showing that microsatellite genotypes of the P. falciparum populations in Grande Comore were substantially different from those in other two islands (Moheli and Anjouan). Currently, no data on temporal changes in genetic diversity of P. falciparum isolates from Grande Comoros after introduction of ACT are available. Herein, the objective of this study is to investigate the dynamics of genetic diversity and multiplicity of infection (MOI) in clinical P. falciparum isolates from Grande Comore during two different periods (2006‒2007 and 2013‒2016) using polymorphic markers of msp-1, msp-2, and msp-3. The data in this study provide insights on parasite diversity and MOI after various malaria control measures.
This study was approved by the Ethics Committees of Comoros Ministry of Health (No. 07-123/VP-MSSPG/DNS) and Guangzhou University of Chinese Medicine (No. 2012L0816). Blood samples were collected from children after obtaining written informed consent from their parents or legal guardians.
Study sites and sample collection
This study was conducted on the Grande Comore Island, Union of Comoros (Fig. 1), that is located in the Indian Ocean off the south-east coast of Africa, to the east of Mozambique and north-west of Madagascar (11°00′–12°00′S, 43°10′–43°35′E). This island has an area of 1147 km2 with about 420,000 inhabitants (2012 estimate). A tropical hot and rainy season occurs from November to April, and a cooler dry season runs from May to October. Annual temperature ranges from 11 to 35 °C and rainfall ranges from 1000 to 3000 mm per year. Malaria transmission on this island is perennial with most of infections occurring during the rainy season. P. falciparum is the dominant malaria species, with occasional Plasmodium malariae and Plasmodium vivax infections .
A total of 232 blood samples from microscopically P. falciparum positive patients visiting local healthcare facilities were collected in two different periods (March 2006–October 2007 and March 2013–December 2016). Among these samples, 118 samples from Mitsoudje Centre Hospital (50), National Malaria Centre (30), and Mitsamiouli Centre Hospital (38) were collected during March 2006–October 2007; the other 114 samples from Mitsoudje Centre Hospital (41), National Malaria Centre (24), Mitsamiouli Centre Hospital (29), and Mbeni District Health Centre (20) were collected during March 2013–December 2016. Venous blood sample of 1.0 ml was collected from each patient in an EDTA tube and stored at −20 °C until DNA extraction.
PCR amplification and allelic analysis of the pfmsp-1, pfmsp-2, and pfmsp-3 genes
Genomic DNA of each blood sample was extracted using Takara DNA Blood Mini Kit according to the manufacturer’s instructions (Takara, Kyoto, Japan). Extracted parasite DNA was dissolved in TE buffer (10 mM Tris–HCl, 0.1 M EDTA, pH 8.0) and stored in microfuge tubes at −20 °C. Segments of the pfmsp-1 (block 2), pfmsp-2 (block 3), and pfmsp-3 were amplified using nest PCR, as described previously [10, 13]. An initial amplification of the outer regions of the three genes was followed by a nested PCR with sequence specific primer pairs. All the reactions were carried out in final volume of 25 μl containing 10.0 μl of dH2O, 0.5 μl of each primer (0.4 µM), 12.5 μl of Taq PCR Mast Mix (2.5 U) following the manufacturer’s instructions (Sangon Bio Inc., Shanghai, China) on a S1000 Thermal cycler (Bio-Rad, Hercules, USA). In the primary amplification reactions, 2.0 μl of template genomic DNA were added as a temple. In the nested reaction, 0.5 μl of primary PCR product was added as a temple. The nested PCR products were separated on 2.0% agarose gel (Sangon Bio Inc., Shanghai, China) and visualized under ultraviolet (UV) trans-illumination. A 100-bp DNA ladder was used to determine the size of PCR products (Sangon Bio Inc., Shanghai, China). MOI of the msp-1 or msp-2 genes was calculated by averaging number of amplified bands per positive P. falciparum isolate as described previously . For DNA sequencing, amplified DNA fragments representing different alleles were purified using a PCR purification kit (Takara, Kyoto, Japan). Purified PCR products from selected isolates representing different alleles of msp-1, msp-2, and msp-3 were directly sequenced in both directions with the primers in the secondary PCR using an ABI PRISM3730 DNA sequencer (Sangon Bio Inc., Shanghai, China). The sequences were also used to correct the estimated molecular weight and to confirm the nature of the amplified product.
Statistical analysis was determined using SPSS (version 13.0) software. Comparisons of MOI of msp-1 or msp-2 in isolates collected between 2006–2007 and 2013–2016 were made using t test. Mann–Whitney U test was used to compare in the frequencies of the mutations and alleles of the msp-1, msp-2 and msp-3 in isolates collected between 2006–2007 and 2013–2016. P < 0.05 was considered indicative of a statistically significant difference.
Allelic polymorphism of msp-1 and msp-2
Approximately 98% of P. falciparum isolates (222/232) were successfully amplified from msp-1 block 2 region (112 for 2006–2007 group; and 110 for 2013–2016 group), and three allelic types (K1, MAD20, and RO33) were identified from the samples (Table 1). The band sizes were 150–300 bp for the K1 type, 140–280 bp for the MAD20 type, and only one band size (~ 150 bp) for RO33 allelic type. For 2006–2007 group, 58 isolates (51.8%) were K1 type; 48 (42.9%) were MAD20 type; and 95 (84.8%) were RO33 type. For the 2013–2016 group, 46 (41.8%), 26 (23.6%), and 70 (63.4%) isolates had K1, MAD20, and RO33 allelic types, respectively. For individual infections among the 2006–2007 group, 6.3% carried K1 type, 2.7% carried MAD20 type, 14.3% carried RO33 type, and 76.7% had multiple allelic types (i.e. K1/MAD20, K1/RO33, MAD20/RO33, and K1/MAD20/RO33). In contrast, 23.6% of the 2013–2016 infections had K1 type, 12.7% had MAD20 type, 34.5% had RO33 type, and 29.1% contained multiple allelic types (i.e. K1/MAD20, K1/RO33, and MAD20/RO33). Compared with the 2006–2007 group, the frequencies of isolates with multiple allelic types was significantly decreased (from 76.7 to 29.1%, P < 0.01) in those from 2013–2016 group, whereas the frequencies of isolates with only K1, MAD20, or RO33 allelic types was significantly increased (P < 0.01).
A total of 222 P. falciparum isolates had a positive PCR outcome for msp-2 block 3, including 109 samples collected during 2006–2007 and 113 samples during 2013–2016. Two allelic types (FC27 and 3D7) were identified in both two periods (Table 1), with PCR product sizes varying from 250 to 550 bp for FC27 allelic type and from 400 to 660 bp for 3D7 allelic type, respectively. For samples collected during 2006–2007, 71.6 and 90.8% carried FC27 and 3D7 allelic types, while 91.1 and 37.1% of the 2013–2016 samples had FC27 and 3D7 allelic types, respectively. For individual infections, 9.2% contained only FC27 type, 28.4% carried only 3D7 type, and the remaining 62.4% had two allelic types in 2006–2007 group (i.e. FC27/3D7). For 2013–2016 group, 62.8% of the infections contained only FC27 type, 8.8% carried only 3D7 type, and the remaining 28.3% had multiple allelic types (i.e. FC27/3D7), again showing a decreasing trend in mixed genotype infection (62.4% in 2006–2007, and 28.3% in 2013–2016; P < 0.01). MOIs for both msp-1 (3.11) and msp-2 (2.75) were higher in P. falciparum isolates in 2006–2007 group than those in 2013–2016 group (MOI = 1.63 and 1.35, respectively).
Sequence analysis of the msp-1 gene
Sequence analysis of msp-1 block 2 region revealed that the Grande Comore isolates contained a limited number of different tripeptide repeat units (SAG, SGT, SGA, and SGP for K1 allelic type; and SGG, SVA, SKG, and SVT for MAD20 allelic type). For the K1 allelic type, the tripeptide repeat region always started with SAQ and terminated with SGT (Fig. 2); while the tripeptide repeat region of the MAD20 type alleles usually started with one of two tripeptide repeat units (SGG or SKG) and always ended with two tripeptide repeat units (SVASGG) (Fig. 3). The allelic diversity in K1 and MAD20 allelic types could be mainly caused by duplications or deletions of these repeat motifs. For the RO33 allelic type, the polymorphisms could be mainly due to amino acid substitutions in the KDGANTQVVAKPA/DAVSTQSAKNPPGATVPSGTASTKGAIRSPGAANPS sequences (Fig. 4).
Sequence analysis of msp-1 block 2 region showed that the isolates could be divided into three allelic types (K1, MAD20, and RO33) (Table 2). In total, 32 distinct sequence variants were found in 2006–2007 group, with 14 variant subtypes for the K1 type, 9 variant subtypes for MAD20, and 9 variant subtypes for RO33. Of these variant subtypes, K1-18 (20.0%), MAD20-1 (25.0%), and RO33-3 (24.2%) subtypes were the most prevalent in the 2006–2007 group. For the 2013–2016 group, a total of 23 distinct variants in msp-1 gene were detected: 11 variant subtypes for the K1 type, 3 for MAD20 type, and 9 for RO33 type. Of the 23 sequence variants, the K1-13, MAD20-1/MAD20-3, and RO33-1 subtypes were the most prevalent, accounting for 20.5, 41.6, and 59.3%, respectively. Compared with the 2006–2007 group, the frequencies of msp-1 K1-1 (P < 0.05), K1-4 (P < 0.01), K1-6 (P < 0.01), K1-10 (P < 0.01), K1-15 (P < 0.05), K1-17 (P < 0.01), and K1-18 (P < 0.05) subtypes were significantly decreased in the 2013–2017 group, whereas the K1-2 (P < 0.05), K1-5 (P < 0.01), K1-8 (P < 0.01), K1-11 (P < 0.01), K1-13 (P < 0.01), and K1-14 (P < 0.05) were significantly increased. For msp-1 MAD20 allelic types, significant decrease in MAD20-4 to MAD20-9 allelic types (P < 0.05 or P < 0.01) and significant increase in MAD20-1 to MAD20-3 types (P < 0.01) were observed in 2013–2016 group compared with those of the 2006–2007 group. Similarly, in comparison to 2006–2007 group, significantly reduced frequencies of RO33-3, RO33-5, RO33-6, RO33-8, and RO33-10 were observed in the 2013–2016 group (P < 0.01), whereas significantly increased of RO33-1 (P < 0.01), RO33-7 (P < 0.05), and RO33-11 (P < 0.01) types were detected in the 2013–2016 group.
Sequence analysis of the msp-2 gene
Sequence analysis of msp-2 block 3 region showed that the FC27 alleles contained a family-specific region (ADTIASGSQSSTNSASTSTTNNGESQTTTPTA), a conserved 8 amino acid region (ADTPTATE), a 12 amino acid repeating unit (SNSRSPPITTE or SNSPSPPITITE, n = 0–4), and followed by a conserved region (SSGNAPNK) at the end (Fig. 5). The polymorphisms of FC27 allelic type were mainly due to the non-synonymous amino acid changes in these regions, whereas the polymorphic character of 3D7 allelic type were mainly variable repeating units of 4–10 amino acids as well as different numbers (n = 5–8) of threonine residues (T) (Fig. 6).
Fragments of the P. falciparum msp-2 block 3 were successfully amplified and sequenced from 118 samples of the 2006–2007 group and 111 samples of the 2013–2016 group. Amino acid sequence analysis revealed that 52.5% (62/118) and 47.5% (56/118) of the samples had FC27 and 3D7 allelic types in the 2006–2007 group (Table 3), with a total of 26 distinct haplotypes (16 for FC27 allelic type; and 13 for 3D7 allelic type). For the 2006–2007 group, the FC27-5/10 (12.9%) and 3D7-2 (34.0%) haplotype was the most abundant. However, in the 2013–2016 group, 80 (72.1%) and 33 (27.9%) samples had the FC27 and 3D7 allelic types of msp-2, respectively. There were 13 different FC27 haplotypes and 8 3D7 haplotypes in the 2013–2016 group (Table 3). Of these haplotypes, the FC27-14 and 3D7-9 were the most abundant, with an overall frequency of 58.8 and 48.5%, respectively. Compared with the 2006–2007 group, the frequencies of FC27-5 (P < 0.01), FC27-6 (P < 0.01), FC27-7 (P < 0.01), FC27-10 (P < 0.01), FC27-11 (P < 0.01), FC27-16 (P < 0.01), and FC27-19 (P < 0.01), and 3D7-2 (P < 0.01), 3D7-4 (P < 0.01), 3D7-5 (P < 0.05), 3D7-13 (P < 0.01), and 3D7-15 (P < 0.05) were significantly decreased in the 2013–2016 group, whereas the frequencies of FC27-14, 3D7-7, 3D7-9, and 3D7-12 haplotypes were significantly increased (P < 0.01). Some of the P. falciparum Grande Comore isolates with msp-2 FC27 and 3D7 haplotypes in this study showed 100% identity with other strains from Vietnam (AAG47596 with FC27-20), while other haplotypes (FC27-1 to FC27-19, FC27-19, 3D7-1 to 3D7-15) were new alleles identified in this study.
Sequence analysis of msp-3
Of the 232 confirmed P. falciparum samples obtained from Grande Comore, 204 were successfully amplified and sequenced for msp-3 (105 for the 2006–2007 group and 99 for the 2013–2016 group). Sequence analysis of msp-3 showed that the allelic diversity of msp-3 allelic types (3D7 and K1) was attributed to variation in size of the AHR region, with K1 allelic type being 366 or 369 bp long and 3D7 type having 339 bp (Fig. 7). For the 2006–2007 group, 66.7% (70/105) and 33.3% (35/105) of the samples contained K1 and 3D7 allelic types, with 8 (K1-1 to K1-8) and 3 (3D7-1 to-3D7-3) distinct haplotypes, respectively (Table 4). Of these haplotypes, the K1-5 and 3D7-1 were the most abundant haplotypes, with an overall frequency of 20.0 and 15.2%, respectively. In contrast, 46 (46.5%) and 53 (53.5%) samples in the 2013–2016 group had K1 and 3D7 allelic types, respectively (Table 4). There were only two different K1 haplotypes (K1 and K2) and one 3D7 haplotype in the 2013–2016 group. The K1-2 and 3D7-1 were the most abundant haplotypes, accounting for 37.4.0 and 53.5% of the isolates examined, respectively. Compared with the 2006–2007 group, the frequencies of K1-3, K1-5, K1-6, 3D7-2 and 3D7-3 was significantly decreased in 2013–2014 group (P < 0.01), whereas the frequencies of K1-2 and 3D7-1 haplotypes was significantly increased (P < 0.01).
Dramatic reduction in annual malaria cases has been achieved in Grande Comore through the use of ACT for the treatment of uncomplicated P. falciparum patients, ACT-based MDA, and other malaria control interventions. However, malaria continues to be one of most important public health problems on this island, which calls for monitoring changes in drug resistance status and parasite population dynamics. Determining P. falciparum genetic diversity and MOI from field samples is important for understanding the impacts of malaria control measures on parasite populations and for developing strategies to better control malaria infection. The present study investigates the temporal change of genetic diversity and MOI of Grande Comore P. falciparum populations based on msp-1, msp-2, and msp-3 genes that have been used to monitor parasite population widely [10, 15]. The present data showed that the frequencies of allelic diversity and MOI of msp-1, msp-2, and msp-3 significantly decreased in the 2013–2016 group when compared with those the 2006–2007 group, which may reflect a trend associated with decreasing malaria transmission intensity on this island.
Some msp-1 K1, MAD20, or RO33 haplotypes found in the present study have been reported in other regions of the world, including Brazil (AFS44739 with K1-5), Senegal (ABS84524 with K1-5, ABS84428 with K1-6, and ABS84432 with K1-11), Malawi (ADQ74224 with K1-12 and ADQ74227 with RO33-1), India (AFR61077 with K1-12 and AFR61093 with MAD20-3), Tanzania (BAM84405 with K1-12 and AAC69748 with RO33-1), Colombia (ACS26173 with MAD20-6), and Myanmar (ACB69813 with MAD20-6); whereas the other allele variants of msp-1 were new alleles identified in this study. The present data reveal RO33 being the dominant allelic type of msp-1 gene among Grande Comore P. falciparum isolates in both groups. RO33 has also been reported to the dominant allele in parasites collected from Malaysia , Brazil , and Gabon , whereas MAD20 allele was the most prevalent in Myanmar [18, 19], Thailand , Iran , Pakistan , and Colombia , Senegal [23, 24]. For parasites from French Guiana , Kenya , and Peru , K1 was the dominant msp-1 allelic type. For msp-2, the 3D7 type was the predominant among isolates in 2006–2007 group, similar to reports from several African countries, including Gambia , Cameroon , Congo , Ghana [31, 32], Senegal [14, 23, 24], Burkina Faso , Malawi , Uganda , and Tanzania . The 3D7 allele was also major type in some Southeast Asian countries such as Cambodia , Iran , Malaysia [13, 36], Myanmar , Pakistan , Papua New Guinea , Thailand , as well as Thai-Myanmar borders . In contrast, the FC27 was the dominant allele for parasites from Gabon , Cameroon , Nigeria . Several previous reports indicate that the FC27 allele is associated with disease severity , and the 3D7 type may be the common genotype circulating in high disease transmission areas . Here the data in this study show that over 10 years after the introduction of ACT, the prevalence of 3D7 allelic type in msp-2 gene is dramatically decreased, from 90.8 to 37.1%, whereas FC27 allelic type increased from 71.6 to 91.1%.
Previous reports show that polyclonal infection is more common in areas with high endemicity, and 50–100% of infections are polyclonal infections in mesoendemic and holoendemic areas [43,44,45]. Furthermore, a significant association between the complexity of infection and polyclonal infections with the asymptomatic malaria was observed in malaria endemic area of Congo . In the present study, more than 76 and 62% of the samples examined harboured polyclonal infections (two or three allelic types) of the msp-1 and msp-2 gene, respectively, in 2006–2007 group. The frequencies of polyclonal infections were reduced to about 29 and 28%, respectively, in the 2013–2016 group, which again suggests decreasing population diversity and/or transmission intensity. MOI is conventional index to measure of complexity of infection and intensity of transmission. A high MOI value is often observed in a hyperendemic region with high malaria transmission [21, 31, 47, 48]. In the present study, the MOI values decreased from 3.11 to 1.63 for msp-1 and from 2.75 to 1.35 for msp-2, respectively. The findings in this study are similar to those reported in southeastern Senegal  and Congo . The present data suggest a progressive decrease of P. falciparum transmission on this island. In fact, according to a report from the Comoros Ministry of Health, the numbers of annual malaria cases in Grande Comore dramatically decreased from 92,480 (in 2006) with high level incidence (about 35%) to 1362 (in 2016) with low level incidence (about 0.3%). Therefore, the present data confirm that MOI can be used as a useful indicator for monitoring malaria transmission level in the endemic areas.
Sequence analysis revealed that the decline in the number of msp-1 haplotypes (32 for the 2013–2016 and 23 for the 2006–2007 group) among Grande Comore isolates. Similarly, the total number of haplotypes in msp-2 dramatically decreased from 29 in the 2006–2007 group (16 for FC27 and 13 for 3D7 allelic types) to 21 (13 for FC27 and 8 for 3D7 allelic types) in 2013–2016 group (about 28% decline). This is in agreement with previous reports in other countries with declining endemicity [49, 50]. However, studies from Senegal, Mozambique, and Iran indicated that the introduction of ACT in Congo has reduced the MOI but not the genetic diversity of msp-2 gene among P. falciparum isolates from children living in Southern districts of Brazzaville . Again, haplotype analysis supports reduced genetic diversity and transmission on the Grande Comore island.
The polymorphism of pfmsp-3 is predominantly confined to sequence diversity in the N-terminal domain within the heptad-repeats (insertion/deletion and nucleotide substitutions) . The present study detected both msp-3 K1 and 3D7, but not recombinant type, similar to those reported from Thailand, Papua New Guinea, India, Keyan . Recombinant msp-3 alleles were detected in Iran and African countries at a very low frequency [51, 52]. Some of the Grande Comore parasites collected in this study had new msp-3 alleles (K1-3 to K1-8, 3D7-2 and 3D7-3) that have not been reported previously. However, many parasites showed 100% identity with those from Asia, Africa, and South America reported previously, such as Thailand (AOT86948 with K1-1, AOT86951 with K1-2, and AOT86944 with 3D7-1), India (AEI28718 with K1-1, AEI28725 with K1-2, and AEI28765 with 3D7-1), Kenya (AMM75906 with K1-1, AMM75893 with K1-2, and AMM75927 with 3D7-1), Nigeria (CAJ44166 with K1-1, CAJ44194 with K1-2, and CAJ44184 with 3D7-1), China (AAF04099 with K1-1), Indonesia (AAF59914 with K1-2), Papua New Guinea (AAC47670 with K1-2, and AAC47662 with 3D7-1), Vietnam (AAK94780 with K1-2), and Brazil (AFP75269 with K1-2). In the present study, the msp-3 3D7-1 haplotype was the most prevalent in both 2006–2007 and 2013–2016 groups. The present data were in line with the findings from Thailand, India, and Nigeria, where 3D7-1 haplotypes was the most abundant types . In the present study, K1 allelic type was the predominant (66.7%) in 2006–2007 group. The data in the present study are in some degree consistent with the reports of K1 being the most prevalent type in the Thailand , Thailand–Myanmar border , and Thailand–Cambodia border , but is contrast to previous reports from in Thailand–Laos border  and Peru , with the 3D7 being the most prevalent type. Over the course of 10 years (from 2006 to 2016), the frequencies of K1 type dramatically decreased from 66.7 to 46.5% (P < 0.01), while the 3D7 type dramatically increased from 33.3 to 53.5% (P < 0.01), suggesting that parasites with the msp-3 3D7 type may survive better after introduction of ACTs in Grande Comore. The total number of haplotypes in msp-3 gene changed from 11 in 2006–2007 to 3 in 2013–2016 (a 60% decline), suggesting a decreasing tend in genetic diversity of msp-3 in Grande Comore after 10 years of use of ACT.
The observation of increased frequencies of msp-2 FC27 and msp-3 3D7 allelic types when general population genetic diversity and other allelic types have decreased are interesting, although we do not know the reason for the shift of the alleles. One remote possibility is that the msp-2 FC27 and/or msp-3 3D7 alleles or some unknown genes nearby (linked to msp-2 and/or msp-3) play a role in parasite response to ACT. Parasites carrying these specific alleles/genes can survive better under drug pressure and increase frequency. Another possibility is that the 2013–2016 parasite populations might consist of some parasites carrying these alleles imported from nearby endemic regions after reduction in the original parasite populations. These issues require further investigations.
This study investigated the temporal change in genetic diversity and MOI of P. falciparum populations in Grande Comore Island after the introduction of ACT using the polymorphic genetic markers (MSP-1, MSP-2, and MSP-3). Results from the current study showed that the prevalence of genetic diversity and MOI in msp-1, msp-2, or msp-3 decreased over the course of the study (July 2006 to July 2016). The data in this study suggest a progressive decrease in genetic diversity likely due to lower malaria transmission intensity. The data presented here provide a valuable information for assessing the appropriateness of the current malarial control strategies in this endemic area.
multiplicity of infection
- msp-1 :
P. falciparum merozoite surface protein gene 1
- msp-2 :
P. falciparum merozoite surface protein gene 2
- msp-3 :
P. falciparum merozoite surface protein gene 3
artemisinin-based combination therapy
WHO. World malaria report 2016. Geneva: World Health Organization; 2016.
Ouledi A. Epidemiology and control of malaria in the Federal Islamic Republic of Comoros (in French). Sante. 1995;5:368–71.
Hoffman S, Vekemans J, Richie T, Duffy P. The march toward malaria vaccines. Vaccine. 2015;33(Suppl 4):D13–23.
Matuschewski K. Vaccines against malaria-still a long way to go. FEBS J. 2017;284:2560–8.
Gosling R, von Seidlein L. The future of the RTS, S/AS01 malaria vaccine: an alternative development plan. PLoS Med. 2016;13:e1001994.
Beeson J, Drew D, Boyle M, Feng G, Fowkes F, Richards J. Merozoite surface proteins in red blood cell invasion, immunity and vaccines against malaria. FEMS Microbiol Rev. 2016;40:343–72.
Genton B, Betuela I, Felger I, Al-Yaman F, Anders R, Saul A, et al. A recombinant blood-stage malaria vaccine reduces Plasmodium falciparum density and exerts selective pressure on parasite populations in a phase 1-2b trial in Papua New Guinea. J Infect Dis. 2002;185:820–7.
Healer J, Murphy V, Hodder A, Masciantonio R, Gemmill A, Anders R, et al. Allelic polymorphisms in apical membrane antigen-1 are responsible for evasion of antibody-mediated inhibition in Plasmodium falciparum. Mol Microbiol. 2004;52:159–68.
Kiwanuka G. Genetic diversity in Plasmodium falciparum merozoite surface protein 1 and 2 coding genes and its implications in malaria epidemiology: a review of published studies from 1997–2007. J Vector Borne Dis. 2009;46:1–12.
Pattaradilokrat S, Sawaswong V, Simpalipan P, Kaewthamasorn M, Siripoon N, Harnyuttanakorn P. Genetic diversity of the merozoite surface protein-3 gene in Plasmodium falciparum populations in Thailand. Malar J. 2016;15:517.
Rebaudet S, Bogreau H, Silaï R, Lepere J, Bertaux L, Pradines B, et al. Genetic structure of Plasmodium falciparum and elimination of malaria, Comoros archipelago. Emerg Infect Dis. 2010;16:1686–94.
Papa Mze N, Ahouidi A, Diedhiou C, Silai R, Diallo M, Ndiaye D, et al. Distribution of Plasmodium species on the island of Grande Comore on the basis of DNA extracted from rapid diagnostic tests. Parasite. 2016;23:34.
Atroosh W, Al-Mekhlafi H, Mahdy M, Saif-Ali R, Al-Mekhlafi A, Surin J. Genetic diversity of Plasmodium falciparum isolates from Pahang, Malaysia based on MSP-1 and MSP-2 genes. Parasit Vectors. 2011;4:233.
Zwetyenga J, Rogier C, Tall A, Fontenille D, Snounou G, Trape J, et al. No influence of age on infection complexity and allelic distribution in Plasmodium falciparum infections in Ndiop, a Senegalese village with seasonal, mesoendemic malaria. Am J Trop Med Hyg. 1998;59:726–35.
Rich S, Ayala F. Population structure and recent evolution of Plasmodium falciparum. Proc Natl Acad Sci USA. 2000;97:6994–7001.
Kimura E, Mattei D, Mana di Santa S, Scherf A. Genetic diversity in the major merozoite surface antigen of Plasmodium falciparum high prevalence of a third polymorphic form detected in strains derived from malaria patients. Gene. 1990;91:57–62.
Kun F, Schmidt-Ott R, Lehman L, Lell B, Luckner D, Greve B, et al. Merozoite surface antigen 1 and 2 genotypes and rosetting of Plasmodium falciparum in severe and mild malaria in Lambarene, Gabon. Trans R Soc Trop Med Hyg. 1998;92:110–4.
Snounou G, Zhu X, Siripoon N, Jarra W, Thaithong S, Brown KN, et al. Biased distribution of msp1 and msp2 allelic variants in Plasmodium falciparum populations in Thailand. Trans R Soc Trop Med Hyg. 1999;93:369–74.
Kang J, Moon S, Kim J, Cho S, Lin K, Sohn W, et al. Genetic polymorphism of merozoite surface protein-1 and merozoite surface protein-2 in Plasmodium falciparum field isolates from Myanmar. Malar J. 2010;9:131.
Zakeri S, Bereczky S, Naimi P, Pedro G, Djadid N, Färnert A, et al. Multiple genotypes of the merozoite surface proteins 1 and 2 in Plasmodium falciparum infections in a hypoendemic area in Iran. Trop Med Int Health. 2005;10:1060–4.
Ghanchi N, Matensson A, Ursing J, Jafri S, Bereczky S, Hussain R, et al. Genetic diversity among Plasmodium falciparum field isolates in Pakistan measured with PCR genotyping of the merozoite surface protein 1 and 2. Malar J. 2010;9:1.
Gomez D, Chaparro J, Rubiano C, Rojas M, Wasserman M. Genetic diversity of Plasmodium falciparum field samples from an isolated Colombian village. Am J Trop Med Hyg. 2002;67:611–6.
Niang M, Loucoubar C, Sow A, Diagne M, Faye O, Faye O, et al. Genetic diversity of Plasmodium falciparum isolates from concurrent malaria and arbovirus co-infections in Kedougou, southeastern Senegal. Malar J. 2016;15:155.
Niang M, Thiam L, Loucoubar C, Sow A, Sadio B, Diallo M, et al. Spatio-temporal analysis of the genetic diversity and complexity of Plasmodium falciparum infections in Kedougou, southeastern Senegal. Parasit Vectors. 2017;10:33.
Ariey F, Chalvet W, Hommel D, Peneau C, Hulin A, Mercereau-Puijalon O, et al. Plasmodium falciparum parasites in French Guiana: limited genetic diversity and high selfing rate. Am J Trop Med Hyg. 1999;61:978–85.
Takala S, Branch O, Escalante A, Kariuki S, Wootton J, Lal A. Evidence for intragenic recombination in Plasmodium falciparum: identification of a novel allele family in block 2 of merozoite surface protein-1: asembo Bay Area Cohort Project XIV. Mol Biochem Parasitol. 2002;125:163–71.
Chenet S, Branch O, Escalante A, Lucas C, Bacon D. Genetic diversity of vaccine candidate antigens in Plasmodium falciparum isolates from the Amazon basin of Peru. Malar J. 2008;7:93.
Conway D, McBride J. Population genetics of Plasmodium falciparum within a malaria hyperendemic area. Parasitology. 1991;103(Pt 1):7–16.
Basco L, Tahar R, Escalante A. Molecular epidemiology of malaria in Cameroon. XVIII. Polymorphisms of the Plasmodium falciparum merozoite surface antigen-2 gene in isolates from symptomatic patients. Am J Trop Med Hyg. 2004;70:238–44.
Ibara-Okabande R, Koukouikila-Koussounda F, Ndounga M, Vouvoungui J, Malonga V, Casimiro PN, et al. Reduction of multiplicity of infections but no change in msp2 genetic diversity in Plasmodium falciparum isolates from Congolese children after introduction of artemisinin-combination therapy. Malar J. 2012;11:410.
Agyeman-Budu A, Brown C, Adjei G, Adams M, Dosoo D, Dery D, et al. Trends in multiplicity of Plasmodium falciparum infections among asymptomatic residents in the middle belt of Ghana. Malar J. 2013;12:22.
Duah N, Matrevi S, Quashie N, Abuaku B, Koram K. Genetic diversity of Plasmodium falciparum isolates from uncomplicated malaria cases in Ghana over a decade. Parasit Vectors. 2016;9:416.
Mwingira K, Nkwengulila G, Schoepflin S, Sumari D, Beck H, Snounou G, et al. Plasmodium falciparum msp1, msp2 and glurp allele frequency and diversity in sub-Saharan Africa. Malar J. 2011;10:79.
Kidima W, Nkwengulila G. Plasmodium falciparum msp2 genotypes and multiplicity of infections among children under 5 years with uncomplicated malaria in Kibaha, Tanzania. J Parasitol Res. 2015;2015:1–6.
Gosi P, Lanteri C, Tyner S, Se Y, Lon C, Spring M, et al. Evaluation pf parasite subpopulations and genetic diversity of msp1, msp2 and glurp genes during and following artesunate monotherapy treatment of Plasmodium falciparum malaria in Western Cambodia. Malar J. 2013;12:403.
Razak M, Sastu U, Norahmad N, Abdul-Karim A, Muhammad A, Munlandy P, et al. Genetic diversity of Plasmodium falciparum populations in malaria declining areas of Sabah, East Malaysia. PLoS One. 2016;11:e0152415.
Barry A, Schultz L, Senn N, Nale J, Kiniboro B, Siba P, et al. High levels of genetic diversity of Plasmodium falciparum populations in Papua New Guinea despite variable infection prevalence. Am J Trop Med Hyg. 2013;88:718–25.
Congpuong K, Sukaram R, Prompan Y, Dornae A. Genetic diversity of the msp-1, msp-2, and glurp genes of Plasmodium falciparum isolates along the Thai–Myanmar borders. Asian Pac J Trop Biomed. 2014;4:598–602.
Issifou S, Rogier C, Adjagba-Olakpo M, Chabi-Worou N, Ntoumi F. Complexity and genetic diversity of Plasmodium falciparum infections in young children living in urban areas of Central and West Africa. Parasitol Res. 2003;90:423–8.
Njama-Meya D, Kamya M, Dorsey G. Asymptomatic parasitemia as a risk factor for symptomatic malaria in a cohort of Ugandan children. Trop Med Int Health. 2004;9:862–8.
Oyedeji S, Awobode H, Anumudu C, Kun J. Genetic diversity of Plasmodium falciparum isolates from naturally infected children in north-central Nigeria using the merozoite surface protein-2 as molecular marker. Asian Pac J Trop Med. 2013;6(8):589–94.
Soulama I, Nébié I, Ouédraogo A, Gansane A, Diarra A, Tiono A, et al. Plasmodium falciparum genotypes diversity in symptomatic malaria of children living in an urban and a rural setting in Burkina Faso. Malar J. 2009;8:135.
Babiker H, Abdel-Muhsin A, Ranford-Cartwright L, Satti G, Walliker D. Characteristics of Plasmodium falciparum parasites that survive the lengthy dry season in eastern Sudan where malaria transmission is markedly seasonal. Am J Trop Med Hyg. 1998;59:582–90.
Paul R, Day K. Mating patterns of Plasmodium falciparum. Parasitol Today. 1998;14:197–202.
Legrand E, Volney B, Lavergne A, Tournegros C, Florent L, Accrombessi D, et al. Molecular analysis of two local falciparum malaria outbreaks on the French Guiana coast confirms the msp1 B-K1/varD genotype association with severe malaria. Malar J. 2005;4:26.
Ekala M, Jouin H, Lekoulou F, Issifou S, Mercereau-Puijalon O, Ntoumi F. Plasmodium falciparum merozoite surface protein 1 (MSP1): genotyping and humoral responses to allele-specific variants. Acta Trop. 2002;81:33–46.
Babiker H, Walliker D. Current views on the population structure of Plasmodium falciparum: implications for control. Parasitol Today. 1997;13:262–7.
Branch O, Takala S, Kariuki S, Nahlen B, Kolczak M, Hawley W, et al. Plasmodium falciparum genotypes, low complexity of infection, and resistance to subsequent malaria in participants in the Asembo bay cohort project. Infect Immun. 2001;69:7783–92.
Anderson T, Williams J, Estrada-Franco J, Richardson L, Mollinedo R, Bockarie M, et al. Microsatellites reveal a spectrum of population structure in the malaria parasite Plasmodium falciparum. Mol Biol Evol. 2000;17:1467–82.
Anthony T, Cox-Singh J, Matusop A, Ratnam S, Shamsul S, Singh B. Fragmented population structure of Plasmodium falciparum in a region of declining endemicity. J Infect Dis. 2005;191:1558–64.
Ebrahimzadeh A, Mohammadi S, Jamshidi A. Allelic forms of merozoite surface protein-3 in Plasmodium falciparum isolates from southeast of Iran, Jundishapur. J Microbiol. 2014;7:e9829.
Soulama I, Bigoga J, Ndiaye M, Bougouma E, Quagraine J, Casimiro P, et al. Genetic diversity of polymorphic vaccine candidate antigens (apical membrane antigen-1, merozoite surface protein-3, and erythrocyte binding antigen-175) in Plasmodium falciparum isolates from western and central Africa. Am J Trop Med Hyg. 2011;84:276–84.
Gardner M, Hall N, Fung E, White O, Berriman M, Hyman R, et al. Genome sequence of the human malaria parasite Plasmodium falciparum. Nature. 2002;419:498–511.
Jordan S, Oliveira A, Hernandez J, Oster R, Chattopadhyay D, Branch O, et al. Malaria immunoepidemiology in low transmission: correlation of infecting genotype and immune response to domains of Plasmodium falciparum merozoite surface protein 3. Infect Immun. 2011;79:2070–8.
BH, JS, and CD designed, organized, supervised the study and analysed data. BH, FT, WW, YL, GW, SH, QZ, HZ, ML, AB, KSA, AMM, QW, ZY, SZ, and QX carried out the field work and preliminary data analysis. BH wrote and drafted manuscript. X-zS analysed data and wrote manuscript. All authors read and approved the final manuscript.
The authors gratefully acknowledge all the doctors, nurses, and patients in Mitsoudje Centre Hospital, National Malaria Centre, Mitsamiouli Centre Hospital, and Mbeni District Health Centre for their participation in this study. The authors also thank Embassy of the People’s Republic of China in the Union of Comoros for the help in this study.
The authors declare that they have no competing interests.
Availability of data and materials
The data supporting the findings of this article are included in the main text.
Consent for publication
All participants have consented to the publication of data in this study.
Ethics approval and consent to participate
This study was approved by the Ethics Committees of Comoros Ministry of Health (No. 07-123/VP-MSSPG/DNS) and Guangzhou University of Chinese Medicine (No. 2012L0816).
This work was supported in part by grants from National Natural Science Foundation of China [No. 81702020], Guangdong Provincial Science Foundation [No. 2015A030310107], and Guangdong Provincial Medicine Science Foundation [No. A2016315] to BH, Guangdong Provincial Science Technology Programme [Nos. 2014B050502013 and 2015B020234003] and State Administration of Traditional Chinese Medicine Project [No. GJZX2016003] to JS, National Natural Science Foundation of China [No. 81403295] and YangFan Innovative and Entrepreneurial Research Team Project [No. 2014YT02S008] to CD, and by the Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health (X-zS). The funds had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Huang, B., Tuo, F., Liang, Y. et al. Temporal changes in genetic diversity of msp-1, msp-2, and msp-3 in Plasmodium falciparum isolates from Grande Comore Island after introduction of ACT. Malar J 17, 83 (2018). https://doi.org/10.1186/s12936-018-2227-3
- Malaria population
- Grande Comore
- Genetic diversity
- Multiplicity of infection
- DNA sequencing