Skip to main content

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

Abstract

Background

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.

Methods

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.

Results

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).

Conclusions

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.

Background

Malaria is a major infectious disease that led to ~ 212 million clinical cases and about 429,000 deaths worldwide in 2016 [1]. 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 [2]. 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 [11], 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.

Methods

Ethics clearance

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 [12].

Fig. 1
figure 1

Maps of Grande Comore Island, Union of Comoros, showing the studied locations: Mitsoudje Centre Hospital (triangle), Mbeni District Health Centre (star), National Malaria Centre (square), and Mitsamiouli Centre Hospital (circle)

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 [14]. 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

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.

Results

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).

Table 1 Prevalence and multiplicity of infection (MOI) of msp-1 and msp-2 allelic types in Plasmodium falciparum isolates from Grande Comore Island in two different periods

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).

Fig. 2
figure 2

Schematic representation of msp-1 K1 distinct allelic variants of Grande Comore Plasmodium falciparum isolates. There were 18 distinct variants based on the number and arrangement of SAQ, SGT, SGA and SGP motifs, with 14 from the 2006‒2007 group (K1-1, 3, 4, 6, 7, 9, 10–13, and 15–18) and 11 from the 2013‒2016 group (K1-2, 5–9, 11–14, and 18)

Fig. 3
figure 3

Schematic representation of msp-1 MAD20 allelic types of Grande Comore Plasmodium falciparum isolates. A total of 9 distinct variants were detected based on the number and arrangement of SGG, SVT, SKG, and SVA motifs, with 9 from the 2006‒2007 group (MAD20-1 to 9) and 3 from the 2013‒2016 group (MAD20-1 to 3)

Fig. 4
figure 4

Sequence alignments of amino acid sequences of msp-1 RO33 distinct allelic variants among Plasmodium falciparum isolates from Grande Comore. A total of 12 distinct variants were found by sequence analysis of msp-1 block 2, including 9 from the 2006‒2007 group (RO33-1 to 6, 8, 10, and 12) and 9 form the 2013‒2016 group (RO33-1 to 5, 7, 9, 11, and 12). Dots and dashes represent identical residues and deletions, respectively

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.

Table 2 Prevalence of msp-1 K1, MAD20, and RO33 allelic variants collected from Plasmodium falciparum isolates along Grande Comore Island in 2006–2007 and 2013–2016 periods

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).

Fig. 5
figure 5

Sequence alignments of amino acid sequences of msp-2 FC27 distinct allelic variants of Plasmodium falciparum isolates in Grande Comore. A total of 21 distinct variants were identified by sequence analysis of msp-2 block 3 in two groups, including 16 for the 2006‒2007 group (FC27-1 to 3, 5–7, 9–14, 16, and 19–21) and 13 for the 2013‒2016 group (FC27-1 to 4, 8–10, 12–15, 17, and 18). Dots and dashes represent identical residues and deletions, respectively

Fig. 6
figure 6

Sequence alignments of amino acid sequences of msp-2 3D7 distinct allelic variants of Plasmodium falciparum isolates in Grande Comore. A total of 15 distinct variants were identified by sequence analysis of msp-2 block 3, including 13 for the 2006‒2007 group (3D7-1 to 6, 8–11, and 13–15) and 8 for the 2013‒2016 group (3D7-1 to 3, 7–9, 11, and 12). Dots and dashes represent identical residues and deletions, respectively

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.

Table 3 Prevalence of msp-2 FC27 and 3D7 allelic variants collected from Plasmodium falciparum isolates along Grande Comore Island in 2006–2007 and 2013–2016 periods

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).

Fig. 7
figure 7

Sequence alignments of amino acid sequences of msp-3 distinct allelic variants of Plasmodium falciparum isolates in Grande Comore. A total of 11 distinct variants were found, including 11 for the 2006‒2007 group (K1-1 to 8; and 3D7-1 to 3) and 3 for the 2013‒2016 group (K1-1 and 2; and 3D7-1). Dots and dashes represent identical residues and deletions, respectively

Table 4 Prevalence of msp-3 K1 and 3D7 allelic variants collected from Plasmodium falciparum isolates along Grande Comore Island in 2006–2007 and 2013–2016 periods

Discussion

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 [13], Brazil [16], and Gabon [17], whereas MAD20 allele was the most prevalent in Myanmar [18, 19], Thailand [18], Iran [20], Pakistan [21], and Colombia [22], Senegal [23, 24]. For parasites from French Guiana [25], Kenya [26], and Peru [27], 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 [28], Cameroon [29], Congo [30], Ghana [31, 32], Senegal [14, 23, 24], Burkina Faso [33], Malawi [33], Uganda [33], and Tanzania [34]. The 3D7 allele was also major type in some Southeast Asian countries such as Cambodia [35], Iran [20], Malaysia [13, 36], Myanmar [19], Pakistan [21], Papua New Guinea [37], Thailand [18], as well as Thai-Myanmar borders [38]. In contrast, the FC27 was the dominant allele for parasites from Gabon [39], Cameroon [40], Nigeria [41]. Several previous reports indicate that the FC27 allele is associated with disease severity [42], and the 3D7 type may be the common genotype circulating in high disease transmission areas [32]. 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 [46]. 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 [24] and Congo [30]. 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 [30]. 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) [10]. 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 [10]. 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 [10]. 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 [10], Thailand–Myanmar border [53], and Thailand–Cambodia border [53], but is contrast to previous reports from in Thailand–Laos border [53] and Peru [54], 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.

Conclusion

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.

Abbreviations

MOI:

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

ACT:

artemisinin-based combination therapy

References

  1. WHO. World malaria report 2016. Geneva: World Health Organization; 2016.

    Google Scholar 

  2. Ouledi A. Epidemiology and control of malaria in the Federal Islamic Republic of Comoros (in French). Sante. 1995;5:368–71.

    CAS  PubMed  Google Scholar 

  3. Hoffman S, Vekemans J, Richie T, Duffy P. The march toward malaria vaccines. Vaccine. 2015;33(Suppl 4):D13–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Matuschewski K. Vaccines against malaria-still a long way to go. FEBS J. 2017;284:2560–8.

    Article  CAS  PubMed  Google Scholar 

  5. Gosling R, von Seidlein L. The future of the RTS, S/AS01 malaria vaccine: an alternative development plan. PLoS Med. 2016;13:e1001994.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. 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.

    Article  PubMed  Google Scholar 

  8. 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.

    Article  CAS  PubMed  Google Scholar 

  9. 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.

    CAS  PubMed  Google Scholar 

  10. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  13. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. 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.

    Article  CAS  PubMed  Google Scholar 

  15. Rich S, Ayala F. Population structure and recent evolution of Plasmodium falciparum. Proc Natl Acad Sci USA. 2000;97:6994–7001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. 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.

    Article  CAS  PubMed  Google Scholar 

  17. 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.

    Article  CAS  PubMed  Google Scholar 

  18. 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.

    Article  CAS  PubMed  Google Scholar 

  19. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  20. 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.

    Article  CAS  PubMed  Google Scholar 

  21. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 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.

    Article  CAS  PubMed  Google Scholar 

  23. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  24. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  25. 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.

    Article  CAS  PubMed  Google Scholar 

  26. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Conway D, McBride J. Population genetics of Plasmodium falciparum within a malaria hyperendemic area. Parasitology. 1991;103(Pt 1):7–16.

    Article  PubMed  Google Scholar 

  29. 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.

    CAS  PubMed  Google Scholar 

  30. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  32. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 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.

    Article  Google Scholar 

  35. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  36. 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.

    Article  Google Scholar 

  37. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  38. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. 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.

    Article  PubMed  Google Scholar 

  40. 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.

    Article  PubMed  Google Scholar 

  41. 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.

    Article  CAS  PubMed  Google Scholar 

  42. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  43. 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.

    Article  CAS  PubMed  Google Scholar 

  44. Paul R, Day K. Mating patterns of Plasmodium falciparum. Parasitol Today. 1998;14:197–202.

    Article  CAS  PubMed  Google Scholar 

  45. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  46. 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.

    Article  CAS  PubMed  Google Scholar 

  47. Babiker H, Walliker D. Current views on the population structure of Plasmodium falciparum: implications for control. Parasitol Today. 1997;13:262–7.

    Article  CAS  PubMed  Google Scholar 

  48. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. 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.

    Article  CAS  PubMed  Google Scholar 

  50. 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.

    Article  PubMed  Google Scholar 

  51. 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.

    Google Scholar 

  52. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  53. 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.

    Article  CAS  PubMed  Google Scholar 

  54. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

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.

Acknowledgements

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.

Competing interests

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).

Funding

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.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Jianping Song or Changsheng Deng.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-018-2227-3

Keywords