Molecular surveillance over 14 years confirms reduction of Plasmodium vivax and falciparum transmission after implementation of Artemisinin-based combination therapy in Papua, Indonesia

Genetic epidemiology can provide important insights into parasite transmission that can inform public health interventions. The current study compared long-term changes in the genetic diversity and structure of co-endemic Plasmodium falciparum and P. vivax populations. The study was conducted in Papua Indonesia, where high-grade chloroquine resistance in P. falciparum and P. vivax led to a universal policy of Artemisinin-based Combination Therapy (ACT) in 2006. Microsatellite typing and population genetic analyses were undertaken on available isolates collected between 2004 and 2017 from patients with uncomplicated malaria (n = 666 P. falciparum and n = 615 P. vivax). The proportion of polyclonal P. falciparum infections fell from 28% (38/135) before policy change (2004-2006) to 18% (22/125) at the end of the study (2015-2017); p<0.001. Over the same period, polyclonal P. vivax infections fell from 67% (80/119) to 35% (33/93); p<0.001. P. falciparum strains persisted for up to 9 years compared to 3 months for P. vivax, reflecting higher rates of outbreeding in the latter. Sub-structure was observed in the P. falciparum population, but not in P. vivax, confirming different patterns of outbreeding. The P. falciparum population exhibited 4 subpopulations that changed in frequency over time. Notably, a sharp rise was observed in the frequency of a minor subpopulation (K2) in the late post-ACT period, accounting for 100% of infections in late 2016-2017. The results confirm epidemiological evidence of reduced P. falciparum and P. vivax transmission over time. The smaller change in P. vivax population structure is consistent with greater outbreeding associated with relapsing infections and highlights the need for radical cure to reduce recurrent infections. The study emphasizes the challenge in disrupting P. vivax transmission and demonstrates the potential of molecular data to inform on the impact of public health interventions.

IgM in human immunity to Plasmodium falciparum malaria

Most studies on human immunity to malaria have focused on the roles of immunoglobulin G (IgG), whereas the roles of IgM remain undefined. Analyzing multiple human cohorts to assess the dynamics of malaria-specific IgM during experimentally induced and naturally acquired malaria, we identified IgM activity against blood-stage parasites. We found that merozoite-specific IgM appears rapidly in Plasmodium falciparum infection and is prominent during malaria in children and adults with lifetime exposure, together with IgG. Unexpectedly, IgM persisted for extended periods of time; we found no difference in decay of merozoite-specific IgM over time compared to that of IgG. IgM blocked merozoite invasion of red blood cells in a complement-dependent manner. IgM was also associated with significantly reduced risk of clinical malaria in a longitudinal cohort of children. These findings suggest that merozoite-specific IgM is an important functional and long-lived antibody response targeting blood-stage malaria parasites that contributes to malaria immunity.

Genomic analysis of a pre-elimination Malaysian Plasmodium vivax population reveals selective pressures and changing transmission dynamics

The incidence of Plasmodium vivax infection has declined markedly in Malaysia over the past decade despite evidence of high-grade chloroquine resistance. Here we investigate the genetic changes in a P. vivax population approaching elimination in 51 isolates from Sabah, Malaysia and compare these with data from 104 isolates from Thailand and 104 isolates from Indonesia. Sabah displays extensive population structure, mirroring that previously seen with the emergence of artemisinin-resistant P. falciparum founder populations in Cambodia. Fifty-four percent of the Sabah isolates have identical genomes, consistent with a rapid clonal expansion. Across Sabah, there is a high prevalence of loci known to be associated with antimalarial drug resistance. Measures of differentiation between the three countries reveal several gene regions under putative selection in Sabah. Our findings highlight important factors pertinent to parasite resurgence and molecular cues that can be used to monitor low-endemic populations at the end stages of P. vivax elimination.

Passively versus Actively Detected Malaria: Similar Genetic Diversity but Different Complexity of Infection

The surveillance of malaria is generally undertaken on the assumption that samples passively collected at health facilities are comparable to or representative of the broader Plasmodium reservoir circulating in the community. Further characterization and comparability of the hidden asymptomatic parasite reservoir are needed to inform on the potential impact of sampling bias. This study explores the impact of sampling strategy on molecular surveillance by comparing the genetic make-up of Plasmodium falciparum and Plasmodium vivax isolates collected by passive versus active case detection. Sympatric isolates of P. falciparum and P. vivax were collected from a large community survey and ongoing clinical surveillance studies undertaken in the hypomesoendemic setting of Mimika District (Papua, Indonesia). Plasmodium falciparum isolates were genotyped at nine microsatellite loci and P. vivax at eight loci. Measures of diversity and differentiation were used to compare different patient and parasitological sample groups. The results demonstrated that passively detected cases (symptomatic) had comparable population diversity to those circulating in the community (asymptomatic) in both species. In addition, asymptomatic patent infections were as diverse as subpatent infections. However, a significant difference in multiplicity of infection (MOI) and percentage of polyclonal infections was observed between actively and passively detected P. vivax cases (mean MOI: 1.7 ± 0.7 versus 1.4 ± 1.4, respectively; P = 0.001). The study findings infer that, in hypomesoendemic settings, passive sampling is appropriate for molecular parasite surveillance strategies using the predominant clone in any given infection; however, the findings suggest caution when analyzing complexity of infection. Further evaluation is required in other endemic settings.

Chloroquine efficacy for Plasmodium vivax in Myanmar in populations with high genetic diversity and moderate parasite gene flow

BACKGROUND

Plasmodium vivax malaria remains a major public health burden in Myanmar. Resistance to chloroquine (CQ), the first-line treatment for P. vivax, has been reported in the country and has potential to undermine local control efforts.

METHODS

Patients over 6 years of age with uncomplicated P. vivax mono-infection were enrolled into clinical efficacy studies in Myawaddy in 2014 and Kawthoung in 2012. Study participants received a standard dose of CQ (25 mg/kg over 3 days) followed by weekly review until day 28. Pvmdr1 copy number (CN) and microsatellite diversity were assessed on samples from the patients enrolled in the clinical study and additional cross-sectional surveys undertaken in Myawaddy and Shwegyin in 2012.

RESULTS

A total of 85 patients were enrolled in the CQ clinical studies, 25 in Myawaddy and 60 in Kawthoung. One patient in Myawaddy (1.2%) had an early treatment failure and two patients (2.3%) in Kawthoung presented with late treatment failures on day 28. The day 28 efficacy was 92.0% (95% CI 71.6-97.9) in Myawaddy and 98.3% (95% CI 88.7-99.8) in Kawthoung. By day 2, 92.2% (23/25) in Myawaddy and 85.0% (51/60) in Kawthoung were aparasitaemic. Genotyping and pvmdr1 CN assessment was undertaken on 43, 52 and 46 clinical isolates from Myawaddy, Kawthoung and Shwegyin respectively. Pvmdr1 amplification was observed in 3.2% (1/31) of isolates in Myawaddy, 0% (0/49) in Kawthoung and 2.5% (1/40) in Shwegyin. Diversity was high in all sites (H _E 0.855-0.876), with low inter-population differentiation (F _ST 0.016-0.026, P < 0.05).

CONCLUSIONS

Treatment failures after chloroquine were observed following chloroquine monotherapy, with pvmdr1 amplification present in both Myawaddy and Shwegyin. The results emphasize the importance of ongoing P. vivax drug resistance surveillance in Myanmar, particularly given the potential connectivity between parasite population at different sites.

Genetic micro-epidemiology of malaria in Papua Indonesia: Extensive P. vivax diversity and a distinct subpopulation of asymptomatic P. falciparum infections

BACKGROUND

Genetic analyses of Plasmodium have potential to inform on transmission dynamics, but few studies have evaluated this on a local spatial scale. We used microsatellite genotyping to characterise the micro-epidemiology of P. vivax and P. falciparum diversity to inform malaria control strategies in Timika, Papua Indonesia.

METHODS

Genotyping was undertaken on 713 sympatric P. falciparum and P. vivax isolates from a cross-sectional household survey and clinical studies conducted in Timika. Standard population genetic measures were applied, and the data was compared to published data from Kalimantan, Bangka, Sumba and West Timor.

RESULTS

Higher diversity (HE = 0.847 vs 0.625; p = 0.017) and polyclonality (46.2% vs 16.5%, p<0.001) were observed in P. vivax versus P. falciparum. Distinct P. falciparum substructure was observed, with two subpopulations, K1 and K2. K1 was comprised solely of asymptomatic infections and displayed greater relatedness to isolates from Sumba than to K2, possibly reflecting imported infections.

CONCLUSIONS

The results demonstrate the greater refractoriness of P. vivax versus P. falciparum to control measures, and risk of distinct parasite subpopulations persisting in the community undetected by passive surveillance. These findings highlight the need for complimentary new surveillance strategies to identify transmission patterns that cannot be detected with traditional malariometric methods.

Molecular Epidemiology of P. vivax in Iran: High Diversity and Complex Sub-Structure Using Neutral Markers, but No Evidence of Y976F Mutation at pvmdr1

BACKGROUND

Malaria remains endemic at low levels in the south-eastern provinces of Iran bordering Afghanistan and Pakistan, with the majority of cases attributable to P. vivax. The national guidelines recommend chloroquine (CQ) as blood-stage treatment for uncomplicated P. vivax, but the large influx of imported cases enhances the risk of introducing CQ resistance (CQR).

METHODOLOGY AND PRINCIPAL FINDINGS

The genetic diversity at pvmdr1, a putative modulator of CQR, and across nine putatively neutral short tandem repeat (STR) markers were assessed in P. vivax clinical isolates collected between April 2007 and January 2013 in Hormozgan Province, south-eastern Iran. One hundred blood samples were collected from patients with microscopy-confirmed P. vivax enrolled at one of five district clinics. In total 73 (73%) were autochthonous cases, 23 (23%) imported cases from Afghanistan or Pakistan, and 4 (4%) with unknown origin. 97% (97/100) isolates carried the F1076L mutation, but none carried the Y976F mutation. STR genotyping was successful in 71 (71%) isolates, including 57(57%) autochthonous and 11 (11%) imported cases. Analysis of population structure revealed 2 major sub-populations, K1 and K2, with further sub-structure within K2. The K1 sub-population had markedly lower diversity than K2 (HE = 0.06 vs HE = 0.82) suggesting that the sub-populations were sustained by distinct reservoirs with differing transmission dynamics, possibly reflecting local versus imported/introduced populations. No notable separation was observed between the local and imported cases although the sample size was limited.

CONCLUSIONS

The contrasting low versus high diversity in the two sub-populations (K1 and K2) infers that a combination of local transmission and cross-border malaria from higher transmission regions shape the genetic make-up of the P. vivax population in south-eastern Iran. There was no molecular evidence of CQR amongst the local or imported cases, but ongoing clinical surveillance is warranted.

Submicroscopic and Asymptomatic Plasmodium Parasitaemia Associated with Significant Risk of Anaemia in Papua, Indonesia

Submicroscopic Plasmodium infections are an important parasite reservoir, but their clinical relevance is poorly defined. A cross-sectional household survey was conducted in southern Papua, Indonesia, using cluster random sampling. Data were recorded using a standardized questionnaire. Blood samples were collected for haemoglobin measurement. Plasmodium parasitaemia was determined by blood film microscopy and PCR. Between April and July 2013, 800 households and 2,830 individuals were surveyed. Peripheral parasitaemia was detected in 37.7% (968/2,567) of individuals, 36.8% (357) of whom were identified by blood film examination. Overall the prevalence of P. falciparum parasitaemia was 15.4% (396/2567) and that of P. vivax 18.3% (471/2567). In parasitaemic individuals, submicroscopic infection was significantly more likely in adults (adjusted odds ratio (AOR): 3.82 [95%CI: 2.49-5.86], p<0.001) compared to children, females (AOR = 1.41 [1.07-1.86], p = 0.013), individuals not sleeping under a bednet (AOR = 1.4 [1.0-1.8], p = 0.035), and being afebrile (AOR = 3.2 [1.49-6.93], p = 0.003). The risk of anaemia (according to WHO guidelines) was 32.8% and significantly increased in those with asymptomatic parasitaemia (AOR 2.9 [95% 2.1-4.0], p = 0.007), and submicroscopic P. falciparum infections (AOR 2.5 [95% 1.7-3.6], p = 0.002). Asymptomatic and submicroscopic infections in this area co-endemic for P. falciparum and P. vivax constitute two thirds of detectable parasitaemia and are associated with a high risk of anaemia. Novel public health strategies are needed to detect and eliminate these parasite reservoirs, for the benefit both of the patient and the community.

Genomic Analysis Reveals a Common Breakpoint in Amplifications of the Plasmodium vivax Multidrug Resistance 1 Locus in Thailand

In regions of coendemicity for Plasmodium falciparum and Plasmodium vivax where mefloquine is used to treat P. falciparum infection, drug pressure mediated by increased copy numbers of the multidrug resistance 1 gene (pvmdr1) may select for mefloquine-resistant P. vivax Surveillance is not undertaken routinely owing in part to methodological challenges in detection of gene amplification. Using genomic data on 88 P. vivax samples from western Thailand, we identified pvmdr1 amplification in 17 isolates, all exhibiting tandem copies of a 37.6-kilobase pair region with identical breakpoints. A novel breakpoint-specific polymerase chain reaction assay was designed to detect the amplification. The assay demonstrated high sensitivity, identifying amplifications in 13 additional, polyclonal infections. Application to 132 further samples identified the common breakpoint in all years tested (2003-2015), with a decline in prevalence after 2012 corresponding to local discontinuation of mefloquine regimens. Assessment of the structure of pvmdr1 amplification in other geographic regions will yield information about the population-specificity of the breakpoints and underlying amplification mechanisms.

Asymptomatic Vivax and Falciparum Parasitaemia with Helminth Co-Infection: Major Risk Factors for Anaemia in Early Life

BACKGROUND

Anaemia in children under five years old is associated with poor health, growth and developmental outcomes. In Papua, Indonesia, where the burden of anaemia in infants is high, we conducted a community survey to assess the association between Plasmodium infection, helminth carriage and the risk of anaemia.

METHODS

A cross sectional household survey was carried out between April and July 2013 in 16 villages in the District of Mimika using a multistage sampling procedure. A total of 629 children aged 1-59 months from 800 households were included in the study. Demographic, symptom and anthropometry data were recorded using a standardized questionnaire. Blood and stool samples were collected for examination.

RESULTS

Of the 533 children with blood film examination, 8.8% (47) had P. vivax parasitaemia and 3.9% (21) had P. falciparum; the majority of children with malaria were asymptomatic (94.4%, 68/72). Soil transmitted helminth (STH) infection was present in 43% (105/269) of children assessed; those with STH were at significantly greater risk of P. vivax parasitaemia compared to those without STH (OR = 3.7 [95%CI 1.5-9.2], p = 0.004). Anaemia (Hb<10 g/dl) was present in 24.5% (122/497) of children and associated with P. vivax parasitaemia (OR = 2.9 [95%CI, 1.7-4.9], p = 0.001), P. falciparum parasitaemia (OR = 4.3 [95%CI, 2.0-9.4], p<0.001), hookworm carriage (OR = 2.6 [95%CI, 1.2-5.8], p = 0.026), Plasmodium-helminth coinfection (OR 4.0 [95%CI, 1.4-11.3], p = 0.008) and severe stunting (OR = 1.9 ([95%CI, 1.1-3.3], p = 0.012).

CONCLUSIONS

Asymptomatic P. vivax and P. falciparum infections and hookworm all contribute to risk of paediatric anaemia in coendemic areas and should be targeted with prevention and treatment programs. The relationship between helminth infections and the increased risk of P. vivax parasitaemia should be explored prospectively.

Characterization of blood dendritic and regulatory T cells in asymptomatic adults with sub-microscopic Plasmodium falciparum or Plasmodium vivax infection

BACKGROUND

Plasmodium falciparum and Plasmodium vivax infections compromise dendritic cell (DC) function and expand regulatory T (Treg) cells in both clinical disease (malaria) and experimental human sub-microscopic infection. Conversely, in asymptomatic microscopy-positive (patent) P. falciparum or P. vivax infection in endemic areas, blood DC increase or retain HLA-DR expression and Treg cells exhibit reduced activation, suggesting that DC and Treg cells contribute to the control of patent asymptomatic infection. The effect of sub-microscopic (sub-patent) asymptomatic Plasmodium infection on DC and Treg cells in malaria-endemic area residents remains unclear.

METHODS

In a cross-sectional household survey conducted in Papua, Indonesia, 162 asymptomatic adults were prospectively evaluated for DC and Treg cells using field-based flow cytometry. Of these, 161 individuals (99 %) were assessed retrospectively by polymerase chain reaction (PCR), 19 of whom had sub-microscopic infection with P. falciparum and 15 with sub-microscopic P. vivax infection. Flow cytometric data were re-analysed after re-grouping asymptomatic individuals according to PCR results into negative controls, sub-microscopic and microscopic parasitaemia to examine DC and Treg cell phenotype in sub-microscopic infection.

RESULTS

Asymptomatic adults with sub-microscopic P. falciparum or P. vivax infection had DC HLA-DR expression and Treg cell activation comparable to PCR-negative controls. Sub-microscopic P. falciparum infection was associated with lower peripheral CD4(+) T cells and lymphocytes, however sub-microscopic Plasmodium infection had no apparent effect on DC sub-set number or Treg cell frequency.

CONCLUSIONS

In contrast to the impairment of DC maturation/function and the activation of Treg cells seen with sub-microscopic parasitaemia in primary experimental human Plasmodium infection, no phenotypic evidence of dysregulation of DC and Treg cells was observed in asymptomatic sub-microscopic Plasmodium infection in Indonesian adults. This is consistent with DC and Treg cells retaining their functional capacity in sub-microscopic asymptomatic infection with P. falciparum or P. vivax in malaria-endemic areas.

Analysis of ex vivo drug response data of Plasmodium clinical isolates: the pros and cons of different computer programs and online platforms

Background

In vitro drug susceptibility testing of malaria parasites remains an important component of surveillance for anti-malarial drug resistance. The half-maximal inhibition of growth (IC₅₀) is the most commonly reported parameter expressing drug susceptibility, derived by a variety of statistical approaches, each with its own advantages and disadvantages.

Methods

In this study, licensed computer programs WinNonlin and GraphPad Prism 6.0, and the open access programs HN-NonLin, Antimalarial ICEstimator (ICE), and In Vitro Analysis and Reporting Tool (IVART) were tested for their ease of use and ability to estimate reliable IC₅₀ values from raw drug response data from 31 Plasmodium falciparum and 29 P. vivax clinical isolates tested with five anti-malarial agents: chloroquine, amodiaquine, piperaquine, mefloquine, and artesunate.

Results

The IC₅₀ and slope estimates were similar across all statistical packages for all drugs tested in both species. There was good correlation of results derived from alternative statistical programs and non-linear mixed-effects modelling (NONMEM) which models all isolate data simultaneously. The user-friendliness varied between packages. While HN-NonLin and IVART allow users to enter the data in 96-well format, IVART and GraphPad Prism 6.0 are capable to analyse multiple isolates and drugs in parallel. WinNonlin, GraphPad Prism 6.0, IVART, and ICE provide alerts for non-fitting data and incorrect data entry, facilitating data interpretation. Data analysis using WinNonlin or ICE took the longest computationally, whilst the offline ability of GraphPad Prism 6.0 to analyse multiple isolates and drugs simultaneously made it the fastest among the programs tested.

Conclusion

IC₅₀ estimates obtained from the programs tested were comparable. In view of processing time and ease of analysis, GraphPad Prism 6.0 or IVART are best suited for routine and large-scale drug susceptibility testing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1173-1) contains supplementary material, which is available to authorized users.

Quantification of Plasmodium ex vivo drug susceptibility by flow cytometry

BACKGROUND

The emergence and spread of multidrug-resistant Plasmodium falciparum and Plasmodium vivax highlights the need for objective measures of ex vivo drug susceptibility. Flow cytometry (FC) has potential to provide a robust and rapid quantification of ex vivo parasite growth.

METHODS

Field isolates from Papua, Indonesia, underwent ex vivo drug susceptibility testing against chloroquine, amodiaquine, piperaquine, mefloquine, and artesunate. A single nucleic acid stain (i.e., hydroethidine (HE) for P. falciparum and SYBR Green I (SG) for P. vivax) was used to quantify infected red blood cells by FC-based signal detection. Data derived by FC were compared to standard quantification by light microscopy (LM). A subset of isolates was used to compare single and double staining techniques.

RESULTS

In total, 57 P. falciparum and 23 P. vivax field isolates were collected for ex vivo drug susceptibility testing. Reliable paired data between LM and FC was obtained for 88 % (295/334) of these assays. The median difference of derived IC50 values varied from -5.4 to 6.1 nM, associated with 0.83-1.23 fold change in IC50 values between LM and FC. In 15 assays (5.1 %), the derived difference of IC50 estimates was beyond the 95 % limits of agreement; in eleven assays (3.7 %), this was attributable to low parasite growth (final schizont count < 40 %), and in four assays (1.4 %) due to low initial parasitaemia at the start of assay (<2000 µl(-1)). In a subset of seven samples, LM, single and double staining FC techniques generated similar IC50 values.

CONCLUSIONS

A single staining FC-based assay using a portable cytometer provides a simple, fast and versatile platform for field surveillance of ex vivo drug susceptibility in clinical P. falciparum and P. vivax isolates.