- Open Access
High prevalence of very-low Plasmodium falciparum and Plasmodium vivax parasitaemia carriers in the Peruvian Amazon: insights into local and occupational mobility-related transmission
Malaria Journal volume 16, Article number: 415 (2017)
The incidence of malaria due both to Plasmodium falciparum and Plasmodium vivax in the Peruvian Amazon has risen in the past 5 years. This study tested the hypothesis that the maintenance and emergence of malaria in hypoendemic regions such as Amazonia is determined by submicroscopic and asymptomatic Plasmodium parasitaemia carriers. The present study aimed to precisely quantify the rate of very-low parasitaemia carriers in two sites of the Peruvian Amazon in relation to transmission patterns of P. vivax and P. falciparum in this area.
This study was carried out within the Amazonian-ICEMR longitudinal cohort. Blood samples were collected for light microscopy diagnosis and packed red blood cell (PRBC) samples were analysed by qPCR. Plasma samples were tested for total IgG reactivity against recombinant PvMSP-10 and PfMSP-10 antigens by ELISA. Occupation and age 10 years and greater were considered surrogates of occupation-related mobility. Risk factors for P. falciparum and P. vivax infections detected by PRBC-qPCR were assessed by multilevel logistic regression models.
Among 450 subjects, the prevalence of P. vivax by PRBC-PCR (25.1%) was sixfold higher than that determined by microscopy (3.6%). The prevalence of P. falciparum infection was 4.9% by PRBC-PCR and 0.2% by microscopy. More than 40% of infections had parasitaemia under 5 parasites/μL. Multivariate analysis for infections detected by PRBC-PCR showed that participants with recent settlement in the study area (AOR 2.1; 95% CI 1.03:4.2), age ≥ 30 years (AOR 3.3; 95% CI 1.6:6.9) and seropositivity to P. vivax (AOR 1.8; 95% CI 1.0:3.2) had significantly higher likelihood of P. vivax infection, while the odds of P. falciparum infection was higher for participants between 10 and 29 years (AOR 10.7; 95% CI 1.3:91.1) and with a previous P. falciparum infection (AOR 10.4; 95% CI 1.5:71.1).
This study confirms the contrasting transmission patterns of P. vivax and P. falciparum in the Peruvian Amazon, with stable local transmission for P. vivax and the source of P. falciparum to the study villages dominated by very low parasitaemia carriers, age 10 years and older, who had travelled away from home for work and brought P. falciparum infection with them.
Over the past decade international funding, ongoing political commitment, and improved diagnostic, prevention and treatment strategies have facilitated worldwide reduction of malaria incidence and mortality . Unfortunately in contrast to this trend, the Loreto Region of Peru (which accounts for 97% of the Peruvian countrywide malaria burden) experienced a fivefold increase in malaria cases, as reported through passive surveillance to the Peruvian Ministry of Health (MoH) between 2011 and 2015, reaching a peak of 54,823 Plasmodium vivax cases in 2014 and 12,646 Plasmodium falciparum cases in 2015 [2, 3].
While the Peruvian Amazon basin has traditionally been classified as a low transmission setting [3,4,5,6,7,8,9], foci of high malaria transmission have been identified in some villages across the Loreto Region [10, 11], where the modulation of parasite density by acquired immunity in the host must occur due to repeated infections as is common in high transmission settings  or due to clonal infections for long periods within these areas [13,14,15].
In contrast to patients with acute malaria disease, who generally seek treatment at MoH facilities (thus treatment immediately interrupts transmission, particularly with P. vivax), individuals with premunition or “clinical immunity” (i.e. lack of typical symptoms of malaria such as fever, chills, sweats and headache in the presence of low parasite density) do not experience overt systemic inflammation or acute malarial disease. Because such individuals do not present for medical attention and anti-malarial treatment, they remain as potential parasite reservoirs that maintain local transmission with micro-geographic or regional movement of parasites [9, 16, 17], undermining regional control and elimination efforts.
In this scenario, where malaria transmission is maintained by a high prevalence of submicroscopic infections, the development and application of tools to identify these reservoirs is an important priority. This is especially important given that standard malaria control systems in the Americas and elsewhere, including that used by the Peruvian MoH, rely on passive case detection (PCD) of symptomatic infections confirmed only by light microscopy at health facilities [18, 19]. Several studies have suggested that multiple metrics are needed to reflect regional and national trends and current malaria transmission [11, 20,21,22,23,24]. Implementation of these must be guided by the precision, accuracy and cost of each metric .
The present study tested the hypothesis that very low parasitaemia of both P. vivax and P. falciparum is common in contrasting epidemiological contexts of the Peruvian Amazon. Further asymptomatic parasite carriage and human occupation-related mobility across areas with different levels of endemicity leads to local reintroduction of Plasmodium infections was assessed. Detailed data to support these hypotheses would be of high generalizable importance in explaining how hypoendemic malaria is maintained on a regional basis. Further such data would support the notion that occupationally mobile, asymptomatically parasitaemic individuals contribute to local movement and regional reintroductions of malaria parasites, making elimination challenging in the absence of targeting anti-malarial interventions at such populations. To address these questions, this study used whole blood samples to estimate an accurate burden of P. vivax and P. falciparum infections in two villages of the Peruvian Amazon through a parasitological and serological survey.
Study area and population
This study was conducted in two sites of Loreto, near the capital city of Iquitos in the Peruvian Amazon: Cahuide and Lupuna (Fig. 1). Cahuide (04°13.785′ S, 73°276′ W) is located 60 km from Iquitos city on the Iquitos-Nauta road in the district of San Juan Bautista (south of Iquitos city). This site is divided into three villages: 12 de April, Cahuide and La Habana, distributed along the Iquitos-Nauta road. Lupuna (03° 44.591′ S, 73° 19.615′ W) is a network of three villages: San José de Lupuna, San Pedro, and Santa Rita, located on the banks of the Nanay River in the district of Punchana (north of Iquitos city); access to this site is only by river. The population of both sites consists mainly of mestizos (individuals that cannot be clearly identified as belonging to a specific ethnic minority but generally being genetic admixture of Hispanic ancestry and indigenous peoples). Living standards are generally based on a subsistence economy including agriculture, fishing and occasional hunting [25, 26].
Malaria transmission in Loreto near to Iquitos has been generally thought to be seasonal. In Lupuna, malaria peaks from November to May, the rainier season [27, 28]. In Cahuide, malaria generally has the same pattern, but, after an outbreak in 2011, has decreased overall . Anopheles darlingi is the principal mosquito vector in both villages with diverse breeding sites, including slow moving parts of larger rivers, smaller streams, pools, and swamps .
This study was carried out within the International Centers of Excellence for Malaria Research (ICEMR) Amazonia project based on a longitudinal, population-based cohort. For this study, a cross-sectional survey was conducted in January 2013, in which subjects with complete surveillance data from the Amazonian-ICEMR census (August 2012) were invited to participate. Of these, participants who consented to donate a blood sample by venipuncture were included in the study. The Amazonian-ICEMR cohort collected blood samples for light microscopy and on filter papers for PCR diagnosis. Participants with incomplete blood samples in December 2012 or January 2013 were excluded (Fig. 2).
A full census of the study population was conducted in August 2012. Individual and household data on socio-demographics (age, gender, education, occupation), self-reported previous history of clinical malaria, and structural characteristics of the household were collected. All households and inhabitants were encoded and geo-referenced using a Global Positioning System (GPS) handheld device (Garmin’s GPSMAP 60CSx, Garmin International Inc., USA).
During the cross-sectional survey conducted in January 2013, data regarding clinical examinations for fever and other malaria symptoms were recorded. Blood samples were taken by finger prick on slides for immediate microscopic examination (thick and thin blood smears). If a participant agreed, an additional 6 mL for adults or 3 mL for children of whole blood was collected by venipuncture in tubes with EDTA (BD Vacutainer, BD Franklin Lakes, USA) as preservative. Venipuncture blood samples were separated by centrifugation (3500 rpm) into plasma for serological analysis and packed red blood cells (PRBC) for molecular diagnosis. Both samples were, respectively, stored at 4 and − 70 °C until processed at the Institute of Tropical Medicine “Alexander von Humboldt”, Lima (ITM-AvH) for molecular (PCR) and serological (ELISA) analyses.
Thick and thin smears were stained for 10 min with a 10% Giemsa solution (Giemsa’s azur eosin methylene blue, Merck KGaA, Germany) using a standard procedure . Microscopic examination was immediately performed in the field by an expert microscopist after the sample collection. Parasite density for each species (parasites/μL blood) was determined by the number of parasites after counting a total of 200 white blood cells (WBC) and assuming an average of 6000 WBC/μL according to the national guidelines. Microscopy fields were read to count at least 500 WBCs before an individual was diagnosed as negative. Quality control was done blindly on all positive slides and 10% of randomly chosen negative slides by a senior technician at ITM-AvH.
Molecular testing for malaria parasitaemia by PCR
PRBC samples were processed using the QIAamp DNA Mini Kit of QIAGEN (PRBC-PCR). For both sample types, subsequent amplification was performed by a real-time quantitative PCR (qPCR) method targeting the 18s rRNA gene region. Oligonucleotides 5-TAACGAACGAGATCTTAA-3 and 5-GTTCCTCTAAGAAGCTTT-3 were used as primers as reported by Mangold et al.  and the PCR conditions consisted of an initial denaturation at 95 °C for 2 min, followed by amplification for 45 cycles of 20 s at 95 °C, 20 s at 52 °C, and 30 s at 68 °C. Amplification was immediately followed by a melt programme consisting of 5 s at 65 °C and a stepwise temperature increase of 0.5 °C/s until 95 °C for species discrimination. Ambiguous melting results were confirmed by using a nested ssPCR method described elsewhere . Parasite density was determined by using a standard curve from sevenfold serial dilutions to 1:10 of a culture sample at concentrations of 2 × 106 parasites/µL down to 2 parasites/µL in blood from an uninfected donor.
Recombinant MSP-10 proteins for P. vivax (PvMSP-10) and P. falciparum (PfMSP-10) (SalI and PF3D7 strains, respectively) were produced in HEK-293 mammal cells (Aragen Biosciences, California, USA). Plasma samples were analysed for total IgG reactivity to recombinant PvMSP-10 and Pf-MSP10 antigens by enzyme-linked immunosorbent assay (ELISA). 96-well ELISA plates (Clear Flat-Bottom Immuno, Thermo Scientific, USA) were coated with 0.5 μg/mL of the recombinant protein in 0.05 M carbonate buffer pH 9.6 overnight at 4 °C. Plates were washed 5 times with phosphate-buffered saline (PBS) containing 0.05% (0.3% for PfMSP-10) Tween 20 (PBST) and blocked for 1 h with 5% skimmed milk in PBST at room temperature. Plates were washed 5 times with PBST and both test and control plasma samples (1:500 diluted in blocking buffer) were added to each well in duplicate for 1 h at room temperature. A 1:10,000 dilution of peroxidase-conjugated goat anti-human IgG (Kirkegaard and Perry Laboratories, Inc., Gaithersburg, MD) was added to each well as a second antibody and incubated for 1 h. After three washings with PBST, bound antibodies were detected by adding the SureBlue tetramethyl benzidine substrate (Kirkegaard and Perry Laboratories, Inc., Gaithersburg, MD). The reaction was then stopped by the addition of 2 M H2SO4. Optical density (OD) at 450 nm was measured using an iMark™ microplate absorbance Reader (Bio-Rad Laboratories).
The corrected OD of each sample was calculated by subtracting the background OD from the corresponding non-coated wells. Serum samples from healthy non-exposed Peruvians were used as negative controls, and positive controls included samples from six different P. vivax and two P. falciparum infected individuals. A mixture model was used to determine OD cutoffs for the seropositivity of each Plasmodium spp. Briefly, mixture models allows to split the OD distribution into two Gaussian-distributed populations using a maximum likelihood approach (a narrow distribution of seronegatives and a broader distribution of seropositives). The cut-off to define seropositivity was the mean OD corresponding to the seronegative population plus 3 standard deviations [33,34,35].
Statistical analyses were conducted in STATA 14 (StataCorp, 2015. Stata Statistical Software: Release 14. College Station, TX). The significance level was defined at 5 and 95% confidence intervals (CI) were estimated whenever appropriate. Fisher’s exact test was used for significance testing of categorical factors for each Plasmodium spp. and a negative binomial regression was used for significance testing of continuous skewed data (i.e. parasite density) . All factors were obtained from the structured questionnaires or laboratory tests, and the report of a previous Plasmodium spp. infection detected by PCR for each participant was obtained from the Amazonian-ICEMR cohort data. Special attention was paid to time spent in the community; it was computed for permanent inhabitants (≥ 6 months living in the community) and a dummy variable was assigned according to recent (≤ 2 years) or long-term (> 2 years) settlement. This cut-off was set in order to evaluate the effect of malaria rebound since 2011 in the Peruvian Amazon Region [11, 29].
To handle the nested structure of sampled data—450 individuals nested within 235 households in 2 communities—a generalized linear mixed effects model (GLMM) was used. Two species-specific models were constructed with Plasmodium infection as the outcome. The dependent variable was defined as the PRBC-PCR (binary) result, and subjects with a different Plasmodium spp. Infection, detected by PRBC-PCR, were excluded in each model.
The univariate and multivariate analyses were fitted with a mixed-effects logistic regression, calculating odds ratios (OR). The suitability of the multilevel structure (two levels: individuals within households; or three levels: individuals within households within communities) was evaluated based on the variance components, intra-class correlation coefficient (ICC), and median odds ratio (MOR) of the null model . The final multilevel structure was applied to both, univariate and multivariate regressions. Associated variables at p < 0.2 in the univariate mixed-effects logistic models for each Plasmodium species were included in its multivariate model construction. Final model variables were retained if p < 0.2 in a backward stepwise process.
Socio-demographic and household characteristics
A total of 450 individuals (Fig. 2) from 252 households were enrolled in Cahuide (62%) and Lupuna (38%). Of these, 46% were male, and participants under 10, between 10 and 30, and over 30 years old represented 32, 39 and 30%, respectively. Just over half of adults (≥ 18 years) had primary school (53%) and the rest secondary school (47%) education. Most adults did not report a single specific occupation (51%), while 24% were farmers and 7% worked in the forest environment (loggers or charcoal workers). Finally, 16% of participants had recently moved and settled in the study areas (≤ 2 years) (Table 1).
Malaria parasite prevalence and exposure
Of the participants with complete surveillance data, a significant number of participants had previous Plasmodium spp. infection diagnosed by PCR on December 2012 (16% P. vivax and 2% P. falciparum). As of January 2013, 28 and 28% of subjects had evidence of recent exposure to P. vivax and P. falciparum, as determined by seropositivity to recombinant MSP-10.
The prevalence of P. vivax by PRBC-PCR (25%) was sevenfold higher than that determined by microscopy (3.6%). For its part, prevalence of P. falciparum by PRBC-PCR (5%) was 24-fold higher than that determined by microscopy (0.2%). The vast majority of subjects with any form of parasitaemia were asymptomatic; few participants (4%) had fever at the time of sampling (Table 2). No mixed-species infections were identified during the study.
Parasite density analysis
Parasite densities by PRBC-PCR were low (estimated geometric mean 57 parasites/μL, 95% CI 29–112). More than 40% of Plasmodium infections detected by PRBC PCR had fewer than 5 parasites/μL and 52% less than 10 parasites/μL. Important differences were observed between P. vivax and P. falciparum infections. The estimated geometric mean was higher in P. vivax (90 parasites/μL, 95% CI 42–190) compared to P. falciparum (8 parasites/μL, 95% CI 2–25). Parasite densities below 10 parasites/μL represented 50% of P. vivax and 68% of P. falciparum infections. The stratified distribution of infections with fewer than 10- and 5-parasites/μL according to some epidemiological features is presented in Fig. 3.
Recent immigrants to the area (≤ 2 years) had significantly higher logs of expected P. vivax parasite densities (coefficient 3.1; 95% CI 1.5:4.7; p < 0.01) compared to permanent settlers after adjusting for community, age, fever and seropositivity while for P. falciparum there were not statistically differences between immigrants and permanent settlers. For P. falciparum, the most important factor associated with a higher logs of expected P. falciparum parasite densities were fever (coefficient 6.4; 95% CI 4.3:8.4; p < 0.01) and seropositivity to PfMSP-10 (coefficient 3.1; 95% CI 1.7:4.6; p < 0.01), after adjusting for community, age and time in community (Fig. 4, Table 3).
Multilevel analysis for malaria infection
The hierarchical structure of data best fit with a multilevel structure with only two levels (individuals within households). Table 4 presents the results of the univariate and multivariate models for malaria infection with this multilevel structure for each Plasmodium species. Univariate analysis showed that participants over 30 years old (OR 3.7; 95% CI 1.9–7.4), seropositives to P. vivax (OR 2.3; 95% CI 1.3–3.8) and those with a P. vivax infection in the previous month (OR 2.1; 95% CI 1.1–3.8) were independently associated with P. vivax infection. On the other hand, only age (≥ 30 years, OR 9.7; 95% CI 1.8–51.5) was independently associated with P. falciparum infection.
The multivariate analysis shows that participants with recent (≤ 2 years) arrival and settlement in the study area [Adjusted odds ratio (AOR) 2.1; 95% CI 1.03–4.2] had significantly higher odds of P. vivax infection, after adjusting for age (≥ 30 years, AOR 3.3; 95% CI 1.6–6.9) and seropositivity to P. vivax (AOR 1.8; 95% CI 1.03–3.2). Regarding P. falciparum infections, participants between 10 and 29 years (AOR 10.7; 95% CI 1.3–91.1) and ≥ 30 years (AOR 16.3; 95% CI 1.9–138.2) had significantly higher odds after adjusting for previous infection (AOR 10.4; 95% CI 1.5–71.1), time in community and fever symptoms (Table 4).
The random effects for each species-specific model are presented in Table 5. The ICC and the MOR of the adjusted models show an important clustering of P. vivax infections in households, an effect that was not observed in P. falciparum.
This study shows, consistent with previous studies, that malaria transmission in the Amazon Region is distributed in highly heterogeneous micro-geographic patterns, with the data suggesting that P. vivax transmitted more locally within villages and sources of P. falciparum more often being distantly acquired and transported on a regional basis. This conclusion is based on the observation that age 10 years and old, a surrogate for occupation-related riverine travel in our study population, is a strong, albeit non-exclusion risk factor for P. vivax vs. P. falciparum infection. While it is not possible to exclude entirely that some P. vivax transmission is acquired elsewhere and conversely that P. falciparum transmission occurs within home villages, the predominant pattern of transmission seems to follow the pattern that P. falciparum is mostly acquired away from and is reintroduced. Overall, the most important finding in this study was the demonstration that a large burden of sub-microscopic infections were detected in a traditionally low-transmission season. As a consequence, these results highlight the high proportion of P. vivax infections found in recent (non-immune) migrants, but that a higher proportion of P. falciparum infection was primarily found in subjects 10 and older, consistent with acquiring infection away from home, related to work. This study reported that relatively less-immune mobile people harbour higher P. vivax parasite densities while non-mobile village inhabitants harbour lower P. vivax parasite densities because they are more likely to develop clinical immunity due to local P. vivax transmission. In contrast, higher P. falciparum parasite densities were observed in inhabitants recently exposed to P. falciparum (PfMSP-10) and also with clinical symptoms, suggesting a lack of clinical immunity. These findings have direct relevance to malaria control and elimination strategies at the population level because they indicate that new public health strategies are required both to prevent local malaria transmission and the regional movement of parasites, particularly P. falciparum.
The sample collection method used in this study allowed for confident identification of associated factors to very low-level parasitaemia infections, which included, for P. vivax, age, time living in community and PvMSP10 seropositivity, and for P. falciparum, age and previous infections. The higher odds for P. vivax infections observed in migrant population (less than 2 years in the community) in comparison to permanent settlers suggest a certain level of naturally-acquired immunity to P. vivax due to local transmission exposure [38, 39]. These estimates were adjusted by a marker of recent exposure (PvMSP-10), hence the increased proportion of P. vivax infections observed in migrants participants was not an effect of recent exposure to parasite, but a lack of development of strain-specific acquired immunity. By contrast, the increased odds of P. falciparum infections according to age, instead of predominantly young people suggest that the local host population lacks immunity against the acquired P. falciparum parasites, for example as observed in imported malaria infections detected in the Solomon Islands [40, 41]. It is important to note that this increase in the P. falciparum infections was observed in populations over 10 years, the typically age where inhabitants in this study site start occupational related activities, and as consequence a high mobility outside community. In addition, the clustering of P. vivax infections at household level observed with the multilevel regression model, also support the hypothesis that P. vivax transmission occurs at local level. This household clustering pattern was observed in other studies with ongoing local malaria transmission [42,43,44]. Importantly, this effect was not observed for P. falciparum infections, suggesting a different transmission pattern for each Plasmodium species.
In areas with local malaria transmission, repeated exposure to infective mosquitoes’ bites carrying malaria parasites, leads to the fitness of the host’ immune system, and thus became both capable to control parasite density and clinical symptoms. Conversely, high parasite densities and clinical manifestations were observed as consequence to the exposition to a new parasite strain. Data reported here suggest that the differences in the parasite densities between P. vivax and P. falciparum presumably arose from the underlying transmission dynamic of each Plasmodium species. The higher P. vivax parasite densities in recent migrants compared to permanent settlers observed in this study is most consistent with a scenario of predominantly local, home-village-based P. vivax transmission. Inhabitants with recent exposure to P. falciparum (PfMSP-10) had higher parasite densities, suggesting a lack of immunity to P. falciparum strains possibly acquired outside communities due to occupational-related activities, as explained above. This important characteristic was not observed in P. vivax, because recent exposure to P. vivax (PvMSP-10) was not associated with high parasite densities. Relapses from hypnozoites could not be excluded as an explanation for the lack of association between parasite density and recent exposure. However, relapses in these areas presumably arose from infections in the recent past (4–10 weeks) [45,46,47,48]. Of note, albeit less marked than in this study, a high burden of sub-microscopic and very-low-densities infections has been reported in other settings in the Peruvian Amazon Region [6, 7, 11]. However, the results of this study show a parasite prevalence of 30% in a region characterized as low-transmission intensity. These findings challenge the traditional classification of the Peruvian Amazon as a hypoendemic malaria transmission setting  and indicate the existence of a sub-microscopic reservoir of infections between peaks of clinical malaria cases.
Taken together, the evidence presented here confirms previous data that malaria is distributed among scattered micro-geographic foci of infection in the Amazon Region . This notion of spatially-constrained transmission is also supported by population genetics studies that report a strong population structure and high diversity and differentiation among communities [14, 50–52]. Studies on the A. darlingi biting behaviours also demonstrated that this vector, the most important in the Peruvian Amazon, has high behavioural plasticity [26, 53, 54], in particular animal biting preferences in diverse habitats [55–57], with biting taking place as far away as 400–500 m from breeding sites [58, 59]. In addition, remarkable differences in the human biting rate were reported at small- to moderate- spatial scales .
This study provides important insights into the contrasting transmission patterns of P. vivax and P. falciparum in the Peruvian Amazon context, particularly points of exposure, with relevance and potential generalizability to other low-transmission contexts. Therefore, future analytical approaches ought to focus on locations and populations at high-risk of Plasmodium species-specific infections where such infections co-exist. The data presented here suggest that optimal control of P. vivax transmission must address the interruption of transmission at the local level on subpatent parasitemics as well as hypnozoite carriers. Indoor residual spraying (IRS) and/or delivery of insecticide-treated bed nets (ITNs) are suitable adjunct control strategies due to highly-efficient sporogonic development of P. vivax within Anophelines . Nonetheless, the hypnozoite reservoir remains a barrier to malaria control due to lack of suitable compounds for mass drug administration capable of eliminating hypnozoites. The most widely used drug to target hypnozoites is primaquine (PQ), but individuals with clinically significant glucose-6-phosphate dehydrogenase deficiency (G6PDd) have a high risk of haemolysis. In Peru, the Ministry of Health mandated treatment with PQ is not typically supervised , but the available evidence suggests low levels of clinically important G6PD mutations among Peruvians . For its part, P. falciparum transmission needs a comprehensive framework to identify ‘source’ and ‘sink’ areas [63–65] to better allocate resources, and an importation vigilance in free-malaria areas to prevent reintroduction.
This study has some limitations. First, novel antigens (PvMSP-10 and PfMSP-10) with proven sensitivity for P. falciparum and P. vivax [66, 67] was used; however, local and regional variability in protective host immunes responses cannot be excluded as a contribution to the findings. Future studies with the use of more antigens as surrogate markers of clinical immunity [35, 68, 69] could address this limitation. Second, because this study was cross-sectional, associations with malaria status should be interpreted with caution as they do not imply causality. Still, the findings presented here are largely in agreement with previous work and local understanding of malaria transmission dynamics in this region. Third, the selection of participants was not based on a random sampling that would assure that missing data were completely at random (MCAR). However, since the exclusion of participants was not be due to having malaria or not, it is reasonable to assume that missing data was at random (MAR). While ideally MCAR is preferable, MAR re-assures us of a low likelihood of bias, at a level that would not mislead the data interpretation or conclusions. Furthermore, no statistically significant differences were observed between census demographics (sex, age, education, and occupation) of participants included and excluded in the present study. Finally, despite the shared vector, A. darlingi, in the Peruvian Amazon, no mixed infections were detected. This is likely due to several factors, including overall low transmission intensity, the micro-geographic localization of P. falciparum, and perhaps even intraspecific competition of both species within the human host. Finding no mixed infections is also possible because of fluctuating parasitaemias and the single point sampling of this study as well as technical limitations of the qPCR method at low-level parasitaemias. Thus, only one dominant Plasmodium species was observed in each test.
This study highlights the varying transmission patterns of P. vivax and P. falciparum in the Peruvian Amazon Region. These findings suggest strong and stable local P. vivax transmission and imported P. falciparum transmission due to occupational-related activities. The detection of a high burden of very low parasitaemia carriers in this study opposes the traditional classification of the Peruvian Amazon Region as a hypo-endemic area and urges a reformulation of current malaria control policies.
Ministry of Health
passive case detection
packed red blood cells
International Centers of Excellence for Malaria Research
intra-class correlation coefficient
median odds ratio
generalized linear mixed effects model
Cohen’s Kappa coefficient
adjusted odds ratio
indoor residual spraying
insecticide-treated bed net
glucose-6-phosphate dehydrogenase deficiency
WHO. World Malaria Report 2016. Geneva: World Health Organization; 2016.
Dirección General de Epidemiología. Ministerio de Salud. Sala Situacional para el Análisis de Situación de Salud 2016-SE03 [Internet]. 2016. http://www.dge.gob.pe/portal/index.php?option=com_content&view=article&id=14&Itemid=121. Cited Jan 8 2016.
Rosas-Aguirre A, Gamboa D, Manrique P, Conn JE, Moreno M, Lescano AG, et al. Epidemiology of Plasmodium vivax Malaria in Peru. Am J Trop Med Hyg. 2016;95:133–44.
Alves FP, Durlacher RR, Menezes MJ, Krieger H, Silva LHP, Camargo EP. High prevalence of asymptomatic Plasmodium vivax and Plasmodium falciparum infections in native Amazonian populations. Am J Trop Med Hyg. 2002;66:641–8.
Vinetz JM, Gilman RH. Asymptomatic Plasmodium parasitemia and the ecology of malaria transmission. Am J Trop Med Hyg. 2002;66:639–40.
Roshanravan B, Kari E, Gilman RH, Cabrera L, Lee E, Metcalfe J, et al. Endemic malaria in the Peruvian Amazon region of Iquitos. Am J Trop Med Hyg. 2003;69:45–52.
Branch O, Casapia WM, Gamboa DV, Hernandez JN, Alava FF, Roncal N, et al. Clustered local transmission and asymptomatic Plasmodium falciparum and Plasmodium vivax malaria infections in a recently emerged, hypoendemic Peruvian Amazon community. Malar J. 2005;4:27.
Chuquiyauri R, Paredes M, Peñataro P, Torres S, Marin S, Tenorio A, et al. Socio-demographics and the development of malaria elimination strategies in the low transmission setting. Acta Trop. 2012;121:292–302.
da Silva-Nunes M, Moreno M, Conn JE, Gamboa D, Abeles S, Vinetz JM, et al. Amazonian malaria: asymptomatic human reservoirs, diagnostic challenges, environmentally driven changes in mosquito vector populations, and the mandate for sustainable control strategies. Acta Trop. 2012;121:281–91.
Parker BS, Paredes Olortegui M, Peñataro Yori P, Escobedo K, Florin D, Rengifo Pinedo S, et al. Hyperendemic malaria transmission in areas of occupation-related travel in the Peruvian Amazon. Malar J. 2013;12:178.
Rosas-Aguirre A, Speybroeck N, Llanos-Cuentas A, Rosanas-Urgell A, Carrasco-Escobar G, Rodriguez H, et al. Hotspots of malaria transmission in the Peruvian Amazon: rapid assessment through a parasitological and serological survey. PLoS ONE. 2015;10:e0137458.
Bousema T, Okell L, Felger I, Drakeley C. Asymptomatic malaria infections: detectability, transmissibility and public health relevance. Nat Rev Microbiol. 2014;12:833–40.
Clark EH, Silva CJ, Weiss GE, Li S, Padilla C, Crompton PD, et al. Plasmodium falciparum malaria in the Peruvian Amazon, a region of low transmission, is associated with immunologic memory. Infect Immun. 2012;80:1583–92.
Delgado-Ratto C, Soto-Calle VE, Van den Eede P, Gamboa D, Rosas A, Abatih EN, et al. Population structure and spatio-temporal transmission dynamics of Plasmodium vivax after radical cure treatment in a rural village of the Peruvian Amazon. Malar J. 2014;13:8.
Lin JT, Saunders DL, Meshnick SR. The role of submicroscopic parasitemia in malaria transmission: what is the evidence? Trends Parasitol. 2014;30:183–90.
Sattabongkot J, Maneechai N, Phunkitchar V, Eikarat N, Khuntirat B, Sirichaisinthop J, et al. Comparison of artificial membrane feeding with direct skin feeding to estimate the infectiousness of Plasmodium vivax gametocyte carriers to mosquitoes. Am J Trop Med Hyg. 2003;69:529–35.
Collins WE, Jeffery GM, Roberts JM. A retrospective examination of the effect of fever and microgametocyte count on mosquito infection on humans infected with Plasmodium vivax. Am J Trop Med Hyg. 2004;70:638–41.
Ministerio de Salud. Norma Técnica de Salud para la Atención de Malaria y Malaria Grave en el Perú. 2007.
WHO. A framework for malaria elimination [Internet]. Geneva: World Health Organization; 2017 [cited 2017 Apr 11]. http://www.who.int/malaria/publications/atoz/9789241511988/en/.
Wanjala CL, Waitumbi J, Zhou G, Githeko AK. Identification of malaria transmission and epidemic hotspots in the western Kenya highlands: its application to malaria epidemic prediction. Parasit Vectors. 2011;4:81.
Cunha MG, Silva ES, Sepúlveda N, Costa SPT, Saboia TC, Guerreiro JF, et al. Serologically defined variations in malaria endemicity in Pará state, Brazil. PloS ONE. 2014;9:e113357.
Tusting LS, Bousema T, Smith DL, Drakeley C. Measuring changes in Plasmodium falciparum transmission: precision, accuracy and costs of metrics. Adv Parasitol. 2014;84:151–208.
Tietje K, Hawkins K, Clerk C, Ebels K, McGray S, Crudder C, et al. The essential role of infection-detection technologies for malaria elimination and eradication. Trends Parasitol. 2014;30:259–66.
Sturrock HJW, Bennett AF, Midekisa A, Gosling RD, Gething PW, Greenhouse B. Mapping malaria risk in low transmission settings: challenges and opportunities. Trends Parasitol. 2016;32:635–45.
Grietens KP, Muela Ribera J, Soto V, Tenorio A, Hoibak S, Aguirre AR, et al. Traditional nets interfere with the uptake of long-lasting insecticidal nets in the Peruvian Amazon: the relevance of net preference for achieving high coverage and use. PLoS ONE. 2013;8:e50294.
Moreno M, Saavedra MP, Bickersmith SA, Lainhart W, Tong C, Alava F, et al. Implications for changes in Anopheles darlingi biting behaviour in three communities in the peri-Iquitos region of Amazonian Peru. Malar J. 2015;14:290.
Ministerio del Ambiente. Servicio Nacional de Metereología e Hidrología del Perú [Internet]. http://www.senamhi.gob.pe/.
Aramburú Guarda J, Ramal Asayag C, Witzig R. Malaria reemergence in the Peruvian Amazon region. Emerg Infect Dis. 1999;5:209–15.
Dirección General de Epidemiología. Ministerio de Salud. Sala Situacional para el Análisis de Situación de Salud 2013-SE03 [Internet]. http://www.dge.gob.pe/portal/index.php?option=com_content&view=article&id=14&Itemid=121. Cited Jan 8 2016.
Ministerio de Salud. Instituto Nacional de Salud. Manual de procedimientos de laboratorio para el diagnóstico de malaria. Lima: MINSA.
Mangold KA, Manson RU, Koay ESC, Stephens L, Regner M, Thomson RB, et al. Real-time PCR for detection and identification of Plasmodium spp. J Clin Microbiol. 2005;43:2435–40.
Singh B, Bobogare A, Cox-Singh J, Snounou G, Abdullah MS, Rahman HA. A genus- and species-specific nested polymerase chain reaction malaria detection assay for epidemiologic studies. Am J Trop Med Hyg. 1999;60:687–92.
Drakeley CJ, Corran PH, Coleman PG, Tongren JE, McDonald SLR, Carneiro I, et al. Estimating medium- and long-term trends in malaria transmission by using serological markers of malaria exposure. Proc Natl Acad Sci USA. 2005;102:5108–13.
Corran PH, Cook J, Lynch C, Leendertse H, Manjurano A, Griffin J, et al. Dried blood spots as a source of anti-malarial antibodies for epidemiological studies. Malar J. 2008;7:195.
Cook J, Speybroeck N, Sochanta T, Somony H, Sokny M, Claes F, et al. Sero-epidemiological evaluation of changes in Plasmodium falciparum and Plasmodium vivax transmission patterns over the rainy season in Cambodia. Malar J. 2012;11:86.
Alexander N. Analysis of parasite and other skewed counts. Trop Med Int Health. 2012;17:684–93.
Merlo J, Chaix B, Ohlsson H, Beckman A, Johnell K, Hjerpe P, et al. A brief conceptual tutorial of multilevel analysis in social epidemiology: using measures of clustering in multilevel logistic regression to investigate contextual phenomena. J Epidemiol Community Health. 2006;60:290–7.
Mueller I, Galinski MR, Tsuboi T, Arevalo-Herrera M, Collins WE, King CL. Natural acquisition of immunity to Plasmodium vivax: epidemiological observations and potential targets. Adv Parasitol. 2013;81:77–131.
Longley RJ, Sattabongkot J, Mueller I. Insights into the naturally acquired immune response to Plasmodium vivax malaria. Parasitology. 2016;143:154–70.
Gray K-A, Dowd S, Bain L, Bobogare A, Wini L, Shanks GD, et al. Population genetics of Plasmodium falciparum and Plasmodium vivax and asymptomatic malaria in Temotu Province, Solomon Islands. Malar J. 2013;12:429.
Waltmann A, Darcy AW, Harris I, Koepfli C, Lodo J, Vahi V, et al. High rates of asymptomatic, sub-microscopic Plasmodium vivax infection and disappearing Plasmodium falciparum malaria in an area of low transmission in Solomon Islands. PLoS Negl Trop Dis. 2015;9:e0003758.
Pongvongsa T, Nonaka D, Iwagami M, Nakatsu M, Phongmany P, Nishimoto F, et al. Household clustering of asymptomatic malaria infections in Xepon district, Savannakhet province, Lao PDR. Malar J. 2016;15:508.
Rulisa S, Kateera F, Bizimana JP, Agaba S, Dukuzumuremyi J, Baas L, et al. Malaria prevalence, spatial clustering and risk factors in a low endemic area of Eastern Rwanda: a cross sectional study. PLoS ONE. 2013;8:e69443.
Carter R, Mendis KN, Roberts D. Spatial targeting of interventions against malaria. Bull World Health Organ. 2000;78:1401–11.
White MT, Shirreff G, Karl S, Ghani AC, Mueller I. Variation in relapse frequency and the transmission potential of Plasmodium vivax malaria. Proc R Soc B. 2016;283:20160048.
Lover AA, Coker RJ. Quantifying effect of geographic location on epidemiology of Plasmodium vivax malaria. Emerg Infect Dis. 2013;19:1058–65.
Battle KE, Karhunen MS, Bhatt S, Gething PW, Howes RE, Golding N, et al. Geographical variation in Plasmodium vivax relapse. Malar. J. 2014;13:144.
White NJ. Determinants of relapse periodicity in Plasmodium vivax malaria. Malar J. 2011;10:297.
Carrasco-Escobar G, Gamboa D, Castro MC, Bangdiwala SI, Rodriguez H, Contreras-Mancilla J, Alava F, Speybroeck N, Lescano AG, Vinetz JM, Rosas-Aguirre A, Llanos-Cuentas A. Micro-epidemiology and spatial heterogeneity of P. vivax parasitaemia in riverine communities of the Peruvian Amazon: A multilevel analysis. Sci Rep. 2017;7(1):8082. doi:10.1038/s41598-017-07818-0.
Kosek M, Yori PP, Gilman RH, Calderon M, Zimic M, Chuquiyauri R, et al. High degree of Plasmodium vivax diversity in the Peruvian Amazon demonstrated by tandem repeat polymorphism analysis. Am J Trop Med Hyg. 2012;86:580–6.
den Eede PV, der Auwera GV, Delgado C, Huyse T, Soto-Calle VE, Gamboa D, et al. Multilocus genotyping reveals high heterogeneity and strong local population structure of the Plasmodium vivax population in the Peruvian Amazon. Malar J. 2010;9:151.
Delgado-Ratto C, Gamboa D, Soto-Calle VE, den Eede PV, Torres E, Sánchez-Martínez L, et al. Population genetics of Plasmodium vivax in the Peruvian Amazon. PLoS Negl Trop Dis. 2016;10:e0004376.
Lainhart W, Bickersmith SA, Nadler KJ, Moreno M, Saavedra MP, Chu VM, et al. Evidence for temporal population replacement and the signature of ecological adaptation in a major Neotropical malaria vector in Amazonian Peru. Malar J. 2015;14:375.
Naranjo-Díaz N, Conn JE, Correa MM. Behavior and population structure of Anopheles darlingi in Colombia. Infect Genet Evol. 2016;39:64–73.
Villarreal-Treviño C, Penilla-Navarro RP, Vázquez-Martínez MG, Moo-Llanes DA, Ríos-Delgado JC, Fernández-Salas I, et al. Larval habitat characterization of Anopheles darlingi from its northernmost geographical distribution in Chiapas, Mexico. Malar J. 2015;14:517.
Hiwat H, Bretas G. Ecology of Anopheles darlingi root with respect to vector importance: a review. Parasit Vectors. 2011;4:177.
Moreno M, Saavedra MP, Bickersmith SA, Prussing C, Michalski A, Tong Rios C, et al. Intensive trapping of blood-fed Anopheles darlingi in Amazonian Peru reveals unexpectedly high proportions of avian blood-meals. PLoS Negl Trop Dis. 2017;11:e0005337.
de Barros FSM, Honório NA. Man biting rate seasonal variation of malaria vectors in Roraima, Brazil. Mem Inst Oswaldo Cruz. 2007;102:299–302.
Achee NL, Grieco JP, Andre RG, Rejmankova E, Roberts DR. A mark-release-recapture study using a novel portable hut design to define the flight behavior of Anopheles darlingi in Belize, Central America. J Am Mosq Control Assoc. 2005;21:366–79.
Sattabongkot J, Tsuboi T, Zollner GE, Sirichaisinthop J, Cui L. Plasmodium vivax transmission: chances for control? Trends Parasitol. 2004;20:192–8.
Grietens KP, Soto V, Erhart A, Ribera JM, Toomer E, Tenorio A, et al. Adherence to 7-Day primaquine treatment for the radical cure of P. vivax in the Peruvian Amazon. Am J Trop Med Hyg. 2010;82:1017–23.
Beck H-P, Wampfler R, Carter N, Koh G, Osorio L, Rueangweerayut R, et al. Estimation of the antirelapse efficacy of tafenoquine, using Plasmodium vivax genotyping. J Infect Dis. 2016;213:794–9.
Wesolowski A, Eagle N, Tatem AJ, Smith DL, Noor AM, Snow RW, et al. Quantifying the impact of human mobility on malaria. Science. 2012;338:267–70.
Wesolowski A, Stresman G, Eagle N, Stevenson J, Owaga C, Marube E, et al. Quantifying travel behavior for infectious disease research: a comparison of data from surveys and mobile phones. Sci Rep. 2014;4:5678.
Tatem AJ, Huang Z, Narib C, Kumar U, Kandula D, Pindolia DK, et al. Integrating rapid risk mapping and mobile phone call record data for strategic malaria elimination planning. Malar J. 2014;13:52.
Torres KJ, Castrillon CE, Moss EL, Saito M, Tenorio R, Molina DM, et al. Genome-level determination of Plasmodium falciparum blood-stage targets of malarial clinical immunity in the Peruvian Amazon. J Infect Dis. 2015;211:1342–51.
Chuquiyauri R, Molina DM, Moss EL, Wang R, Gardner MJ, Brouwer KC, et al. Genome-scale protein microarray comparison of human antibody responses in Plasmodium vivax relapse and reinfection. Am J Trop Med Hyg. 2015;93:801–9.
Corran P, Coleman P, Riley E, Drakeley C. Serology: a robust indicator of malaria transmission intensity? Trends Parasitol. 2007;23:575–82.
Smith DL, Drakeley CJ, Chiyaka C, Hay SI. A quantitative analysis of transmission efficiency versus intensity for malaria. Nat Commun. 2010;1:108.
GCE and JMA conceived and designed the study. JMA and CFM conducted the qPCR assays. DG, KT and MC supervised the serology assays. GCE analysed the data. GCE, JMV, JMA, CFM, KCB and ALC wrote the manuscript. All authors read and approved the final manuscript.
We thank all study participants and the local authorities for their enthusiastic support provided to our study as well as all field workers for their dedication during the surveys. We also thank all members of the Amazonian-ICEMR project in Lima and Iquitos. We particularly thank Dr. Kailash Patra for providing the HEK-produced MSP10 proteins, which he designed and carried out.
The authors declare that they have no competing interests.
Availability of data and materials
All datasets are presented in the main paper and the supplementary materials.
Ethics approval and consent to participate
This study is part of the Amazonian-ICEMR and was approved by the Ethics Review Board of the Regional Health Direction of Loreto, Universidad Peruana Cayetano Heredia in Lima and the Human Subjects Protection Program of the University of California, San Diego, USA. Every adult provided written informed consent and participants under 18 years old provided informed assent in addition to their parents informed consent, prior to subject enrollment and blood sampling.
This work was supported by the International Centers of Excellence for Malaria Research (ICEMR) grant U19AI089681, Fondo Nacional de Desarrollo Científico Tecnológico y de Innovación Tecnológica (FONDECYT) Contract N° 031-2015-FONDECYT and BAYER AG Peru social responsibility grant.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Carrasco-Escobar, G., Miranda-Alban, J., Fernandez-Miñope, C. et al. High prevalence of very-low Plasmodium falciparum and Plasmodium vivax parasitaemia carriers in the Peruvian Amazon: insights into local and occupational mobility-related transmission. Malar J 16, 415 (2017). https://doi.org/10.1186/s12936-017-2063-x
- Plasmodium vivax
- Plasmodium falciparum
- Molecular epidemiology
- Human mobility