Next Article in Journal
The Effects of Short-Term Visual Feedback Training on the Stability of the Roundhouse Kicking Technique in Young Karatekas
Previous Article in Journal
Effects of the Health-Awareness-Strengthening Lifestyle Program in a Randomized Trial of Young Adults with an At-Risk Mental State
Previous Article in Special Issue
Spatio-Temporal Dynamics of Plasmodium falciparum and Plasmodium vivax in French Guiana: 2005–2019
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Malaria in Cambodia: A Retrospective Analysis of a Changing Epidemiology 2006–2019

1
National Center for Epidemiology and Population Health, Research School of Population Health, College of Health and Medicine, Australian National University, Canberra 0200, Australia
2
Institut Pasteur du Cambodge, Phnom Penh 12201, Cambodia
3
Espace-Dev, IRD, University Antilles, University Guyane, University Montpellier, University Réunion, Phnom Penh 12201, Cambodia
4
School of Public Health, The National Institute of Public Health, Phnom Penh 12152, Cambodia
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(4), 1960; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18041960
Submission received: 10 December 2020 / Revised: 22 January 2021 / Accepted: 12 February 2021 / Published: 18 February 2021
(This article belongs to the Special Issue Geo-Epidemiology of Malaria)

Abstract

:
Background: In Cambodia, malaria persists with changing epidemiology and resistance to antimalarials. This study aimed to describe how malaria has evolved spatially from 2006 to 2019 in Cambodia. Methods: We undertook a secondary analysis of existing malaria data from all government healthcare facilities in Cambodia. The epidemiology of malaria was described by sex, age, seasonality, and species. Spatial clusters at the district level were identified with a Poisson model. Results: Overall, incidence decreased from 7.4 cases/1000 population in 2006 to 1.9 in 2019. The decrease has been drastic for females, from 6.7 to 0.6/1000. Adults aged 15–49 years had the highest malaria incidence among all age groups. The proportion of Plasmodium (P.) falciparum + Mixed among confirmed cases declined from 87.9% (n = 67,489) in 2006 to 16.6% (n = 5290) in 2019. Clusters of P. falciparum + Mixed and P. vivax + Mixed were detected in forested provinces along all national borders. Conclusions: There has been a noted decrease in P. falciparum cases in 2019, suggesting that an intensification plan should be maintained. A decline in P. vivax cases was also noted, although less pronounced. Interventions aimed at preventing new infections of P. vivax and relapses should be prioritized. All detected malaria cases should be captured by the national surveillance system to avoid misleading trends.

1. Introduction

Malaria, an Anopheles mosquito-borne disease, remains a serious public health issue globally, although it is preventable and curable [1]. According to the World Health Organization (WHO), an estimated 228 million malaria cases occurred worldwide in 2018, compared with 251 million cases in 2010 and 231 million cases in 2017 [2]. Malaria was responsible for the loss of 405,000 lives in 2018, 416,000 lives in 2017, and 585,000 lives in 2010 worldwide. Morbidity and mortality rates due to malaria vary from one global region to another. The overall malaria incidence rate was 57.4 per 1000 at-risk population, ranging from 2.6 per 1000 at-risk population in the WHO Western Pacific Region covering 37 countries of Oceania, East Asia, and Southeast Asia to 229.3 per 1000 at-risk population in the WHO African Region covering 47 countries in Africa in 2018 [2,3,4]. The deaths due to malaria were lowest at 0.4 per 100,000 population in the WHO Region of the Americas and highest at 41.0 per 100,000 in the WHO Africa Region. The global mortality rate was 10.2 per 100,000 population in 2018. [2,3,5].
Historically, Cambodia has had a high disease burden due to malaria since the 1950s [6,7]. Cambodia benefits from a tropical climate and is home to a variety of mosquito species, including the Anopheles species [8]. Malaria vectors mainly live in forests close to the Vietnamese, Laotian, and Thai borders [9,10]. People living or working in or near forested areas are at higher risk of contracting malaria [11,12].
Cambodia has been widely recognized as successful in the fight against malaria because of the dramatic drop in malaria cases from a peak of over 100,000 cases or 7.4 per 1000 population in 2006 to over 62,000 cases or 3.9 per 1000 population in 2018. Deaths due to malaria also dropped from around 400 reported deaths in 2006 to zero reported malaria deaths in 2017 and 2018 [13,14]. However, Cambodia still has the highest malaria burden compared to its neighboring countries—Vietnam, Thailand, and Laos—and the WHO Western Pacific Region overall [2]. The decline in malaria morbidity and mortality led the Ministry of Health to acknowledge the impact of malaria interventions combined with prevention, economic growth, improved infrastructure, and strengthening of the health system [13], all of which have enabled universal access to timely malaria diagnosis and treatment [13].
Built on this successful experience, since the year 2011, the government has worked with local and international partners on an elimination strategy [6,7]. Implementations were guided by the “National Strategic Plan for Elimination of Malaria in the Kingdom of Cambodia, 2011–2025” and the revised “Cambodia Malaria Elimination Action Framework, 2016–2020”. Cambodia has divided the country’s malaria-endemic areas into four zones based on the annual parasite index and malarial multi-drug resistance status [6,7]. All zones rely on quality diagnosis and treatment at health facilities, village malaria workers (VMWs), and vector control by promoting malaria prevention education and distribution of long-lasting insecticide nets [6,7]. However, to be more efficient and cost-effective, specific interventions were added for each zone [6,7]. First, the elimination zone focused on “transmission interruption” through household testing, vector control, and the expansion of VMWs [6,7]. Second, the pre-elimination zone focused on universal access to diagnosis and treatment by expanding VMWs [6,7]. Third, the reduction zone focused on the intensification plan—re-activating and scaling up VMWs, switching to new antimalarial drugs (dihydroartemisinin–piperaquine (DHA-PIP) to artesunate–mefloquine (AS-MQ)), and creating mobile malaria workers (MMWs) dedicated to the diagnosis and treatment of malaria among forest goers [7]. Finally, the non-endemic zone, defined as areas where no local transmission had been reported, uses government-run healthcare facilities’ basic services for diagnosis and treatment [6,7].
As part of the elimination plan, the Ministry of Health licensed the private sector to test and treat malaria patients between 2011 and 2018 to increase early diagnosis and treatment accessibility. This initiative was known as “Private Public Mix (PPM)” [7]. PPM provided private point-of-care, including private health facilities, pharmacies, and other retailers with malaria test kits and antimalarial drugs free of charge or with a subsidy [7,15,16]. Private point-of-care sites charged patients for their service provision [6,7,13]. These were rapidly scaled up in all endemic areas [15]. According to available data through malaria outlet surveys between 2011 and 2015, private point-of-care sites detected over 50% of malaria patients in the country [15]. PPM was later banned from implementing malaria testing and treatment in April 2018 [7,15]. The reason for the ban was due to malarial drug resistance concerns and the inflexibility experienced in changing testing and treatment policy in the private sector compared to government-run healthcare facilities [7].
The road to success has been fraught with challenges. One main challenge has been antimalarial drug resistance [17,18,19,20,21,22]. Like other countries in the Greater Mekong sub-region—China, Myanmar, Thailand, Laos, and Vietnam—Cambodia is facing the continual threat of antimalarial drug resistance [17,18,19,20,21,22,23. Cambodia first detected parasite resistance to artemisinin in 2006, but retrospective analysis of molecular markers suggested the resistance may have occurred as early as 2001 [6,7,23,24,25]. After nine years of usage, the first-line treatment—AS-MQ—was replaced by DHA-PPQ in Pailin province, in 2008, and nationally in 2010 [6,7,26]. However, the efficacy of DHA-PPQ dropped after a short period of usage [19,27,28,29]. In 2014, the cure rate of DHA-PPQ ranged from 37.5% in Siem Reap province to 89.9% in Pursat province [6,7,28]. In consultation with its partners, the Ministry of Health decided to re-introduce AS-MQ as the first-line treatment nationally in 2016 [7]. However, full coverage was delayed until 2018 due to procurement issues and funding gaps [7,30]. Cambodia relies heavily on external funding. About 70% of malaria control funding was from the Global Fund, while the remaining 30% was from a combination of government funds and other partners [6,7]. Core activities were severely affected during an external funding interruption in 2015–2016, with all activities supported by the Global Fund being suspended for more than a year. This likely influenced the increase in malaria notifications that followed in 2017 and 2018 [2,14].
In the global fight against malaria, the WHO, in the global technical strategy for malaria 2016–2030, has called for the use of surveillance systems as a core intervention in all malaria-endemic settings [31]. Data from the malaria surveillance systems are crucial to guide interventions, planning, and resource allocation [32,33,34,35,36,37,38,39,40,41]. In 2019, malaria cases in Cambodia were reported through two surveillance systems under the Ministry of Health’s umbrella. Those two systems are the Health Management Information System (HMIS) [42] and the Malaria Information System [43]. Created in 1993, with multiple enhancements, HMIS captures data from all health services, including out-patient and in-patient services at government-run healthcare facilities [42]. Specifically for malaria, a more reliable HMIS was developed in 2004. Based on these data, this study aimed to describe how malaria incidence has evolved spatially and how the sequence of untoward events and malaria control interventions may have affected geospatial trends.

2. Materials and Methods

We undertook a secondary data analysis of malaria data from all government healthcare facilities in Cambodia between 2006 and 2019. Data from private providers are not captured in this dataset.

2.1. Data Sources

2.1.1. Malaria Data

The Epidemiology Unit of the National Centre for Parasitology, Entomology and Malaria Control provided HMIS malaria data for the cases recorded between 1 January 2006 and 31 December 2019. These aggregated data included monthly malaria cases categorized by age group, gender, type of diagnosis (rapid diagnostic test (RDT), microscope), species (Plasmodium (P.) falciparum, P. vivax, and Mixed meaning having both P. falciparum and P. vivax), points of care for diagnosis (public hospitals, health centers, and VMWs), case classification (uncomplicated, severe, death [44]), and hospitalization (yes, no).
HMIS is under the supervision of the Ministry of Health’s Department of Planning and Health Information. It is the largest health information system in Cambodia, collecting information from all government-run services, including out-patient, in-patient, and national program data. At the hospital level, where computers are available, data are collected and entered by hospital staff. However, at the health center or health post level, data are recorded in patients’ logbooks. Health center data are merged with data from health posts and VMWs by health center staff and aggregated using a standardized line-listing form before monthly submission to the operational district. Operational districts are responsible for entering data in the HMIS system, which is automatically synchronized to the national database of the Department of Planning and Health Information. The operational districts were created by the Ministry of Health to provide a local level to Cambodia’s health system administration (central/national level, provincial level, and operational district level) [42]. One operational district manages health facilities in one or more administrative districts and covers populations between 100,000 and 200,000 [42]. Of note, HMIS did not record malaria cases diagnosed by the private sector.

2.1.2. Population Data and Maps

The population data, disaggregated by gender and age group at the district level, were obtained from an open-access dataset provided by the United Nations for Coordination of Humanitarian Affairs (OCHA) [45]. Population data were only available for 2016 [45]. To obtain population data for other years, we projected backward and forward using the National Institute of Statistics’ national population growth rate [46]. We also used the district map of Cambodia from OCHA [45].

2.1.3. Untoward Events

We gathered untoward events to contextualize the malaria surveillance data. Information such as treatment failures, key interventions, and events that may have impacted the case detection were obtained from the gray literature review. We mostly relied on governmental documents, including strategic plans and strategic plan reviews for additional untoward events presented in this paper.

2.2. Inclusion Criteria and Setting

All government-run healthcare facilities (hospitals and health centers) reporting to the HMIS were included in this study. The number of healthcare facilities increased from 1183 in 2006 to 1350 in 2019 (Table 1). During this period, some health centers were upgraded to referral hospitals.
These health facilities are located in 24 provinces and the capital city of Phnom Penh. This includes 197 administrative districts and a total of 94 operational districts. A health center covers a catchment area of between 10,000 to 20,000 people [42].

2.3. Case Definition

According to the National Treatment Guideline in Cambodia, confirmed malaria cases were persons who tested positive for P. falciparum, P. vivax, or both (“Mixed”) by RDT or microscopy [47].
Before 2014, unconfirmed malaria cases were also counted in the HMIS system. The definition of an unconfirmed case was a person who had not been tested for malaria with a diagnosis based only on signs and symptoms. Unconfirmed cases were counted if they:
had fever, chills, or sweats or two of the following: headache, nausea, vomiting, diarrhea, AND
any of the following: traveled to the forest in the previous month, had confirmed malaria in the past 28 days, traveled to a malaria-endemic area from a non-endemic area, or lived or worked around others with a recently confirmed malaria diagnosis [47].
Our analysis included all cases (confirmed and unconfirmed) for overall incidence rates and rates by sex, age, district, and seasonality. We, however, used only confirmed cases to disaggregate by Plasmodium species as it was not possible to identify the species among unconfirmed cases.

2.4. Data Analysis

We calculated the annual malaria incidence (disaggregated by species) at the national and district levels by dividing each year’s confirmed cases by their respective populations and obtaining incidences by gender and age group.
We investigated the existence of clusters [48,49] by testing two hypotheses: (1) H0 = malaria cases in each district throughout the country are proportional to the population, and (2) HA = malaria cases in one or more districts are statistically higher than expected proportional to the population. We ran the cluster analysis using a spatial Poisson model with SaTScan™ version 9.6 (Information Management Services Inc., Calverton, MD, USA) at the administrative district scale, using annual malaria cases and population data [48,49]. We introduced the assumption that the cluster’s radius had to be smaller than 100 km or cover less than 50% of the country’s total malaria cases to mitigate the effect of large cluster sizes on the analysis. Such clusters were split into smaller clusters. Another assumption we introduced was that the relative risk had to be greater than three. This restriction allowed us to focus the analysis on high-burden clusters. The default of Monte Carlo replication of 999 was used with a cut-off p-value < 0.001 [48,49]. In this publication, we use the term “cluster” in reference to SaTScan’s “most likely cluster”, which are clusters detected during the first run. In the Poisson model, after the first run, SaTScan removes all data inside the detected clusters and treats them as “no location”, “no case”, “no population”. It then runs again to detect clusters with the remaining data. These newly detected clusters are called “secondary clusters”. Since we are only interested in high-burden areas or “most likely clusters”, we did not allow the option to display “secondary clusters” in our analysis result.
We used Stata version 15.1 (StataCorp, College Station, TX, USA) for descriptive data analysis and R version 4.0 (R Foundation for Statistical Computing, Vienna, Austria) to produce graphics and maps.

3. Results

3.1. Overall Malaria Cases and Incidence Rate

A total of 737,210 malaria cases (confirmed or unconfirmed) was reported to HMIS between 2006 and 2019 (3.6/1000 population) (Figure 1). The highest incidence was reported in 2006, with 7.4 cases/1000 population (100,322 cases), and the lowest incidence was reported in 2016, with 1.5 cases/1000 (23,367 cases). However, the incidence did not decrease regularly over the time period, with a sharp increase to 3.9 cases/1000 in 2018 (62,582 cases) followed by a drop to 1.9 in 2019 (32,597 cases). Between 2006 and 2013, unconfirmed malaria cases were included in the HMIS. Unconfirmed cases made up a large proportion of all reported cases between 2006 and 2010, and subsequently dropped between 2011 and 2013. Nevertheless, the number of confirmed cases was particularly high in 2006, and stable between 2007 and 2011. Between 2007 and 2011, the total number of malaria cases was quite stable, with a drop in unconfirmed cases being compensated by an increase in the confirmed cases. Since 2014, all reported malaria cases have met the confirmed case definition.
After the Private Public Mix (PPM) program was initiated in 2011, the malaria incidence trends (in which data from the private sector were not included) dropped dramatically from 4.3 cases/1000 in 2011 to 1.6 cases/1000 in 2013. As a significant proportion of malaria cases were detected by private sector during this period, the decreasing trend in incidence is unlikely to represent the actual incidence of malaria in Cambodia for this period of time.
Therefore, data represented in this study must be considered in light of the changes in the surveillance case definition over time and changes in reporting sources.
The proportion of P. falciparum + Mixed peaked at 87.9% (n = 67,489) of all confirmed malaria cases recorded in 2006 to reach the lowest rate of 16.6% (n = 5290) of all recorded cases in 2019. The proportion of P. falciparum + Mixed was highest between 2006 and 2009, when the country was using P. falciparum-only RDTs. The proportion of P. vivax + Mixed exceeded the proportion of P. falciparum + Mixed between 2011 and 2014, and then again from 2017 onwards. The change occurred a year after Cambodia started using dual RDT in late 2009, which detected both P. falciparum and P. vivax. More detailed information can be found in the Supplementary Materials (Figure S1).

3.2. Incidence by Sex

The incidence of malaria (confirmed + unconfirmed cases per 1000 population) among women decreased steadily from 6.7/1000 in 2006 to 0.5/1000 in 2019 (average decrease of 4.6% per year since 2006). This was not observed for men, whose incidence fluctuated widely between a maximum of 8.2/1000 in 2006 and a minimum of 2.5/1000 in 2013 (Figure 2). Malaria incidence among males was consistently higher throughout the analysis period.

3.3. Incidence by Age Group

Adults aged 15–49 years had the highest malaria incidence (confirmed + unconfirmed) among all age groups. The incidence in this age group decreased from 9.7 cases/1000 in 2006 to 3.1/1000 in 2019 and was lowest at 2.5 in 2013 and 2016 (Figure 3). In other age groups (under five years, 5–14 years, 49 years and above), malaria incidence rates presented similar trends, albeit much lower rates.

3.4. Seasonality

Both P. falciparum + Mixed and P. vivax + Mixed presented a similar seasonal pattern over time, except in 2016 and 2019 (Figure 4). Fewer P. falciparum + Mixed and P. vivax + Mixed cases were recorded between February and May, while the high malaria season was observed between June and January. A more detailed visualization of cases by month is available in the Supplementary Materials (Figure S2).

3.5. Geographical Distribution

3.5.1. Clusters of P. falciparum + Mixed

The number of confirmed P. falciparum + Mixed cases ranged from zero to 4002 by district per year (Figure 5). The districts with a high burden of malaria, over 1000 cases per year, were seen along the national borders, in the western provinces (Pursat, Kampong Speu, Koh Kong, and Pailin), northern provinces (Preah Vihear and Udor Mean Chey), northeastern provinces (Kratie, Stung Treng), and eastern provinces (Mondul Kiri and Ratanakiri) (Figure 5). The majority of districts with a low malaria burden, lower than 250 cases per year, were observed in central and southern provinces (Figure 5).
The cluster analysis revealed a decrease in the size and number of the clusters of P. falciparum + Mixed cases during the last three years, 2017, 2018, and 2019 (Figure 6). Between 2017 and 2019, these clusters were detected in districts within the northeastern and eastern provinces (Kratie, Stung Treng, Mondul Kiri, and Ratanak Kiri) and districts within the western provinces (Kampong Speu and Pursat).

3.5.2. Clusters of P. vivax + Mixed

The number of confirmed P. vivax + Mixed cases ranged from zero to 5800 cases by district per year. As for P. falciparum + Mixed, the districts with a high burden of malaria, over 1000 cases per year, were seen along the national borders, in the western provinces (Pursat, Kampong Speu, Koh Kong, and Pailin), the northern provinces (Preah Vihear and Udor Mean Chey), northeastern provinces (Kratie, and Stung Treng), and the eastern provinces (Mondul Kiri and Ratanak Kiri) (Figure 7).
As for P. falciparum, fewer and smaller clusters of P. vivax + Mixed cases were detected in the last three years (2017–2019) (Figure 8). Between 2017 and 2019, clusters were also seen along the national borders, in the northern provinces (Preah Vihear and Udor Mean Chey), northeastern provinces (Kratie, and Stung Treng), eastern provinces (Mondul Kiri and Ratanak Kiri), and the western provinces (Pursat, Kampong Speu, and Koh Kong).

4. Discussion

This study describes how malaria has evolved spatially from 2006 to 2019 in Cambodia. The peak malaria seasons were found to be between June and January. A similar finding was described by Maude et al. in 2014, based on Cambodia’s national malaria data between 2004 and 2013 [50]. The rainy season starts one month before the malaria season commences (May vs. June), while it ends two months before the malaria season ends (November vs. January). Previous studies in China, Tibet, and Niger found that rainfall contributed to increased Anopheles and malaria incidence [51,52,53,54].
Our analysis detected P. falciparum + Mixed and P. vivax + Mixed clusters in seven provinces along national borders with Vietnam, Laos, and Thailand. These areas are covered by evergreen broadleaf forests [55] where P. falciparum- and P. vivax-infected Anopheles are prevalent [8]. This geographical distribution of malaria in Cambodia was previously known [50]. However, this analysis provides updated information confirming that these areas, in the elimination phase, are the top priority. The population’s livelihood in the greater incidence regions remains dependent on the forests, including logging precious woods, agriculture activities, and residency on the fringe of the forests [11]. To reduce P. falciparum and P. vivax, the intensification plan, which was effective in 2019, should be maintained, and innovative approaches aimed at blocking transmission from forests to communities should be considered. The main limitation of this cluster analysis is that we can only locate cases at the place of diagnosis, i.e., the location of health facilities, and it is not possible to know where the case contracted malaria. This is a common limitation when using surveillance data.
The overall decrease in malaria has been observed mainly in women, while the annual incidence has fluctuated in men, with the highest rate observed among men aged 15–49 years old. Although the number of cases was similar in men and women in 2006 (8.2 vs. 6.7/1000 population), the gap increased considerably in 2019 (3.5 vs. 0.6/1000 population). The number of cases in women represented only 14% of the total cases. This notable difference in malaria incidence in males and females is likely related to a change in exposure. The evolution of women’s daily activities may have become less related to forests, with adult males being the most exposed to forest activities and, therefore, most at risk of malaria infections [11,12,56]. According to malaria surveys of forest goers, 72.2% in 2004, 70.7% in 2007, 69.9% in 2010, 78.8% in 2013, and 85.3% in 2017 were males aged 15 years old or older while 17.1% in 2004, 20.7% in 2007, 20.6% in 2010, 15.6% in 2013, and 9.7% in 2017 were females aged 15 years old or older [57,58,59,60]. From these figures, the proportion of females who go to forests has decreased by half since 2010. Another possible reason could be the increasing distance between villages and forests due to deforestation. In Cambodia, two primary vectors—Anopheles dirus and Anopheles minimus—play an important role in spreading malaria [61]. Anopheles dirus is more efficient (higher percentage of mosquitoes tested positive for P. falciparum and P. vivax) than Anopheles minimus [61]. They also have different habitats, with Anopheles dirus living in natural forests and forest fringes, and Anopheles minimus living around rice fields and forest fringes [61,62]. With Cambodia undergoing rapid deforestation in recent years [62], the territory of Anopheles dirus could decrease, leading to a lower risk of malaria transmission. However, deforestation activities correspond to an influx of loggers into Anopheles dirus habitat which can lead to increases in malaria cases. This highlights the importance of interventions targeting forest goers to block malaria transmission in forests.
Several key factors may have impacted the malaria trends, although our study design does not provide evidence to prove causality.
First, we observed that all the upward trends in malaria have coincided with periods of high treatment failure rates for P. falciparum, while all of the observed declines in incidence occurred one to two years after a change in first-line treatment. It is plausible that the delayed parasite clearance due to ineffective antimalarial drugs increased malaria transmission among high-risk populations. In 2009, the total malaria cases increased by 29% (n = 16,493) compared to 2008. This rebound occurred after the first-line antimalarial—AS-MQ—showed high treatment failure rates [7,24,25,61,62,63,64,65,66,67,68,69,70,71]. A new first-line antimalarial (DHA-PIP) was introduced nationally in 2010 to replace AS-MQ [7,15,72]. This was followed by an observed reduction of P. falciparum cases in Cambodia until 2013. DHA-PIP treatment failures emerged after three years of usage [27,72,73]. Total malaria cases increased again by 68% (n = 16,679) in 2014 compared to 2013. Malaria notifications continued to increase in Cambodia until 2018. Of note, the ineffective DHA-PIP was still used in several parts of Cambodia due to procurement challenges until 2018 [7,30]. This finding suggests that antimalarial drug efficacy is key to controlling P. falciparum. Timely surveillance of drug efficacy provides critical information to respond to a malaria epidemics promptly.
Another key factor was that the dramatic drop in malaria cases was likely associated with the changing use of the private sector and an important limitation of the HMIS that did not capture data from the private sector. The private sector played a significant role in changing the national malaria trends between 2011 and 2018. However, the most substantial impact would have occurred between 2011 and 2013, when malaria incidence dropped from 4.1 to 1.6/1000 population with a sharper decrease in P. falciparum notifications compared to P. vivax. This drop was possibly due to the Private Public Mix (PPM) project piloted in 2011 and scaled up in 2012 [7,15,16]. It may have favored larger access to early detection and treatment of infections since more than 50% of antimalarial drugs were delivered at private points of care between 2011 and 2013 [16,73]. This is consistent with the findings from national malaria surveys that showed 7.3% of malaria patients used the private sector in 2007, 41.1% in 2010, and 56.4% in 2013, with a reduction to 30.3% in 2017 [57,58,59,60]. Therefore, the decreasing trend of malaria cases recorded in the HMIS between 2011 and 2013 could be related to the increasing use of the private sector, which was not included in the HMIS before it was banned in April 2018. This highlights the importance of comprehensive surveillance with mandatory reporting requirements for public and private sectors to the HMIS to avoid undercounting or misleading trends.
The decreasing trends may have been overestimated due to the inclusion of unconfirmed cases between 2006 and 2013, when the malaria test was not available throughout the country. Unconfirmed cases may include non-malaria cases with similar clinical symptoms. The malaria cases were overestimated if the unconfirmed cases were included and underestimated if the unconfirmed cases were excluded. However, there are no data on the proportion of actual cases among the unconfirmed cases. Changing case definitions and diagnostics has a significant impact on surveillance data and therefore it is important to recognize such events when interpreting surveillance data.
The role of VMWs also played a key role in malaria trends in Cambodia. The malaria incidence increased from 2.7/1000 population in 2014 to 3.3/1000 population in 2015. This increase happened after the scaling up of VMWs in 2014 [7,15,74]. VMWs provided free services in villages and reported to the HMIS through health centers, leading to greater malaria notifications. VMWs administered 41% of all antimalarial drugs in 2015 [73]. The scaling up of VMWs helped change health-seeking behaviors with a noted decrease in private sector services and improved HMIS sensitivity.
In contrast, malaria cases dropped to the lowest point in 2016 but increased again in 2017 and 2018. This occurred after many VMWs were laid off, and all other activities supported by Global Fund funding were suspended between July 2015 and December 2016 [30,75]. The immediate drop in malaria cases in 2016 may be due to the absence of notifications from VMWs. This finding implied that early access to diagnosis and treatment through VMWs or care points is another critical feature to control and eliminate malaria. Funding to maintain VMWs, or other equivalent interventions, should be prioritized until elimination is confirmed.
One more key factor was the intensification plan initiated in 2018 to respond to increasing case notifications in 2017 and 2018. A significant drop in malaria cases from 3.9/1000 in 2018 to 1.9/1000 in 2019 was observed after the plan started. For the first time, both species—P. falciparum and P. vivax—had a dramatic decline. The success of the new intensification plan targeting forest goers provides sufficient evidence that this plan should be maintained.
The final key factor, nationwide crackdowns by police and forest rangers on illegal logging, may also have had an impact on malaria trends. The massive crackdown started in 2019. It prevented the highest risk groups (forest goers) from entering forests where most malaria transmission was occurring [74]. It is of interest that the effectiveness of malaria control programs may be positively affected by external and independent factors such as anti-deforestation activities.
Concerning the Plasmodium species, an increase in the proportion of P. vivax was observed over time. Two reasons may have contributed to these changes. Firstly, the index of species was likely affected by the shift from P. falciparum-only RDT to dual RDT, which could have detected both P. falciparum and P. vivax. As previously shown, P. falciparum notifications were predominant before 2010. This was due to Cambodia using P. falciparum-only RDT before 2010 [76]. In that period, P. vivax cases could only be detected using microscopy, for which availability was limited [77]. According to an outlet survey in 2015, microscopy was available in only 28% of government health facilities included in the study. The dual RDT was only supplied in late 2009 [76]. One year after using the dual RDT, the proportion of P. vivax exceeded the proportion of P. falciparum. However, P. vivax possibly remains under-detected due to the RDT’s low sensitivity (ability to detect true positives) to P. vivax [78]. RDTs do not perform well in individuals with low parasitemia (<100 parasites/μL) [78]. According to a meta-analysis study, the sensitivity of RDT in detecting P. vivax ranges from 57% to 77% compared to a polymerase chain reaction (PCR) as a gold standard [78]. Secondly, the proportion of malaria patients with P. vivax has increased. Relapses may explain these trends. Relapses among P. vivax patients were described as early as 1893 [79]. According to Taylor et al., along the Thailand–Myanmar border, neighboring countries of Cambodia, 75% of P. vivax patients receiving standard treatment (without primaquine radical curative treatment) relapsed within a 12-month follow-up period [80]. One P. vivax patient is thought to have multiple episodes of relapse, but the exact number of relapses per lifetime per person has not yet been confirmed [81]. The accumulation of P. vivax relapses constitutes a large proportion of all P. vivax cases and maintains the source of P. vivax infections. A radical cure for P. vivax, piloted in four provinces of Cambodia (Kampong Speu, Kampong Chhnang, Battambang, and Pailin), should also be prioritized in all P. vivax clusters.
In contrast to P. vivax, the proportion of P. falciparum notifications declined faster over time. The decline in P. falciparum notifications may have resulted from multiple factors. Firstly, the enhanced P. falciparum elimination efforts by the Ministry of Health have likely had a key role to play, such as interventions including expanding testing and treating sites to facilitate the early diagnosis and treatment of cases [6,7,13]. In addition, a focus on the interruption of P. falciparum transmission, or the 1-3-7 strategy (reporting a confirmed case within one day, investigating within three days, and taking measures to prevent further transmission within seven days), in elimination areas may have also contributed to this decline [7,82,83]. Secondly, P. falciparum is curable. The impact from effective interventions in reducing P. falciparum cases should be rapidly seen. Thirdly, biologically, there is evidence that prior exposure to P. vivax suppresses the course of P. falciparum infection [84,85,86,87]. It is possible that this factor contributed to the decline in P. falciparum. However, the proportion by which P. falciparum is reduced by prior exposure to P. vivax is not known.
The overall P. falciparum prevalence decrease was concurrent with the emergence of drug resistance in P. falciparum. The scaling up of improved access to early testing and treatment probably had a greater impact than the emergence of drug resistance. Ongoing surveillance of antimalarial drug resistance also ensured that first-line treatments were promptly changed to efficient regimens.

5. Conclusions

In this study, we used the nationwide surveillance data collected in the last 14 years between 2006 and 2019. There was a noted decrease in notifications of P. falciparum in 2019, suggesting that the intensification plan should be maintained. P. vivax showed a slower but promising declining trend. Interventions aimed at achieving P. falciparum elimination and preventing new P. vivax infections and relapses should be prioritized. In the context of malaria elimination, all cases detected outside the national system should be reported to the national system to avoid misleading trends.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1660-4601/18/4/1960/s1, Figure S1: Proportions of malaria confirmed cases by Plasmodium species in Cambodia between 2006 and 2019, Figure S2: Malaria cases (confirmed + unconfirmed) by month (2006–2019), from all government healthcare facilities in Cambodia.

Author Contributions

Conceptualization, S.C., P.P., T.H., V.H. and B.T.; Data curation, B.T.; Formal analysis, S.C.; Funding acquisition, P.P.; Methodology, S.C., P.P., T.H., V.H. and B.T.; Supervision, P.P., T.H. and V.H.; Validation, S.C.; Visualization, S.C.; Writing—original draft, S.C.; Writing—review & editing, S.C., P.P., T.H., V.H. and B.T. All authors have read and agreed to the published version of the manuscript.

Funding

This article was developed to fulfill partial competency of Master of Philosophy in Applied Epidemiology at the Australian National University. The lead author was funded by an Australian scholarship through the Australian Government’s Department of Foreign Affairs and Trade.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Restrictions apply to the availability of these data. Data was obtained from the National Center for Parasitology Entomology and Malaria Control (CNM) and were made available by Bun Bunkea Tol (co-author, [email protected]) with the permission of CNM.

Acknowledgments

We would like to thank the staff at the Epidemiology Unit of the National Centre for Parasitology, Entomology and Malaria Control for providing the 14-year malaria surveillance data. We are also grateful to the staff at the health facilities who are involved in data collection.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AS-MQ: artesunate-mefloquine, AS-PYR: artesunate-pyronaridine, DHA-PIP: dihydroartemisinin-piperaquine, HMIS: Health Management Information System, MMW: mobile malaria worker, P. falciparum: Plasmodium falciparum, P. vivax: Plasmodium vivax, VMW: village malaria worker, PPM: Private Public Mix.

References

  1. Centers for Disease Control and Prevention (CDC) Malaria. Available online: https://www.cdc.gov/malaria/about/faqs.html (accessed on 15 August 2020).
  2. WHO. World Malaria Report; World Health Organization: Geneva, Switzerland, 2019. [Google Scholar]
  3. World Health Organization (WHO). WHO African Region. Available online: https://www.who.int/immunization/monitoring_surveillance/data/AFR/en/ (accessed on 28 August 2020).
  4. World Health Organization (WHO). WHO Western Pacific Region. Available online: https://www.who.int/westernpacific (accessed on 28 August 2020).
  5. World Health Organization (WHO). Region of the Americas. Available online: https://www.who.int/choice/demography/american_region/en/ (accessed on 29 August 2020).
  6. National Centre for Parasitology Entomology and Malaria Control (CNM). Cambodia Malaria Elimination Action Framework (2016–2020); CNM: Phnom Penh, Cambodia, 2016.
  7. National Centre for Parasitology Entomology and Malaria Control (CNM). National Malaria Porgram Review; CNM: Phnom Penh, Cambodia, 2019.
  8. St. Laurent, B.; Oy, K.; Miller, B.; Gasteiger, E.B.; Lee, E.; Sovannaroth, S.; Gwadz, R.W.; Anderson, J.M.; Fairhurst, R.M. Cow-baited tents are highly effective in sampling diverse Anopheles malaria vectors in Cambodia. Malar. J. 2016, 15, 440. [Google Scholar] [CrossRef] [Green Version]
  9. Durnez, L.; Mao, S.; Denis, L.; Roelants, P.; Sochantha, T.; Coosemans, M. Outdoor malaria transmission in forested villages of Cambodia. Malar. J. 2013, 12, 329. [Google Scholar] [CrossRef] [Green Version]
  10. Hii, J.; Rueda, L.M. Malaria vectors in the Greater Mekong Subregion: Overview of malaria vectors and remaining challenges. Southeast Asian J. Trop. Med. Public Health 2013, 44, 73–165, discussion 306–307. [Google Scholar]
  11. Bannister-Tyrrell, M.; Gryseels, C.; Sokha, S.; Dara, L.; Sereiboth, N.; James, N.; Thavrin, B.; Ly, P.; Soy Ty, K.; Peeters Grietens, K.; et al. Forest Goers and Multidrug-Resistant Malaria in Cambodia: An Ethnographic Study. Am. J. Trop. Med. Hyg. 2019, 100, 1170–1178. [Google Scholar] [CrossRef] [Green Version]
  12. Incardona, S.; Vong, S.; Chiv, L.; Lim, P.; Nhem, S.; Sem, R.; Khim, N.; Doung, S.; Mercereau-Puijalon, O.; Fandeur, T. Large-scale malaria survey in Cambodia: Novel insights on species distribution and risk factors. Malar. J. 2007, 6, 37. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. National Centre for Parasitology Entomology and Malaria Control (CNM). The National Strategic Plan for Elimination of Malaria in the Kingdom of Cambodia 2011–2025; CNM: Phnom Penh, Cambodia, 2011.
  14. WHO. World Malaria Report; World Health Organization: Geneva, Switzerland, 2018. [Google Scholar]
  15. ACTwatch Cambodia 2009–2015 Reference Document. Available online: http://www.actwatch.info/projects/actwatch/cambodia (accessed on 21 November 2019).
  16. Group, A.C.; Phok, S.; Phanalasy, S.; Thein, S.T.; Likhitsup, A. Private sector opportunities and threats to achieving malaria elimination in the Greater Mekong Subregion: Results from malaria outlet surveys in Cambodia, the Lao PDR, Myanmar, and Thailand. Malar. J. 2017, 16, 180. [Google Scholar]
  17. Dondorp, A.M.; Fairhurst, R.M.; Slutsker, L.; Macarthur, J.R.; Breman, J.G.; Guerin, P.J.; Wellems, T.E.; Ringwald, P.; Newman, R.D.; Plowe, C.V. The threat of artemisinin-resistant malaria. N. Engl. J. Med. 2011, 365, 1073–1075. [Google Scholar] [CrossRef] [Green Version]
  18. Mahase, E. Malaria drugs left ineffective by spread of multidrug resistant parasites in southeast Asia. BMJ 2019, 366, l4807. [Google Scholar] [CrossRef]
  19. van der Pluijm, R.W.; Imwong, M.; Chau, N.H.; Hoa, N.T.; Thuy-Nhien, N.T.; Thanh, N.V.; Jittamala, P.; Hanboonkunupakarn, B.; Chutasmit, K.; Saelow, C.; et al. Determinants of dihydroartemisinin-piperaquine treatment failure in Plasmodium falciparum malaria in Cambodia, Thailand, and Vietnam: A prospective clinical, pharmacological, and genetic study. Lancet Infect. Dis. 2019, 19, 952–961. [Google Scholar] [CrossRef] [Green Version]
  20. Hamilton, W.L.; Amato, R.; van der Pluijm, R.W.; Jacob, C.G.; Quang, H.H.; Thuy-Nhien, N.T.; Hien, T.T.; Hongvanthong, B.; Chindavongsa, K.; Mayxay, M.; et al. Evolution and expansion of multidrug-resistant malaria in southeast Asia: A genomic epidemiology study. Lancet Infect. Dis. 2019, 19, 943–951. [Google Scholar] [CrossRef]
  21. Imwong, M.; Hien, T.T.; Thuy-Nhien, N.T.; Dondorp, A.M.; White, N.J. Spread of a single multidrug resistant malaria parasite lineage (PfPailin) to Vietnam. Lancet Infect. Dis. 2017, 17, 1022–1023. [Google Scholar] [CrossRef]
  22. Imwong, M.; Suwannasin, K.; Kunasol, C.; Sutawong, K.; Mayxay, M.; Rekol, H.; Smithuis, F.M.; Hlaing, T.M.; Tun, K.M.; van der Pluijm, R.W.; et al. The spread of artemisinin-resistant Plasmodium falciparum in the Greater Mekong subregion: A molecular epidemiology observational study. Lancet Infect. Dis. 2017, 17, 491–497. [Google Scholar] [CrossRef] [Green Version]
  23. Denis, M.B.; Tsuyuoka, R.; Poravuth, Y.; Narann, T.S.; Seila, S.; Lim, C.; Incardona, S.; Lim, P.; Sem, R.; Socheat, D.; et al. Surveillance of the efficacy of artesunate and mefloquine combination for the treatment of uncomplicated falciparum malaria in Cambodia. Trop. Med. Int. Health TM IH 2006, 11, 1360–1366. [Google Scholar] [CrossRef] [PubMed]
  24. Noedl, H.; Se, Y.; Schaecher, K.; Smith, B.L.; Socheat, D.; Fukuda, M.M. Evidence of artemisinin-resistant malaria in western Cambodia. N. Engl. J. Med. 2008, 359, 2619–2620. [Google Scholar] [CrossRef]
  25. Lim, P.; Chim, P.; Sem, R.; Nemh, S.; Poravuth, Y.; Lim, C.; Seila, S.; Tsuyuoka, R.; Denis, M.B.; Socheat, D.; et al. In vitro monitoring of Plasmodium falciparum susceptibility to artesunate, mefloquine, quinine and chloroquine in Cambodia: 2001–2002. Acta Trop. 2005, 93, 31–40. [Google Scholar] [CrossRef]
  26. Bustos, M.D.; Wongsrichanalai, C.; Delacollette, C.; Burkholder, B. Monitoring antimalarial drug efficacy in the Greater Mekong Subregion: An overview of in vivo results from 2008 to 2010. Southeast Asian J. Trop. Med. Public Health 2013, 44, 201–230, discussion 306–307. [Google Scholar]
  27. Amaratunga, C.; Lim, P.; Suon, S.; Sreng, S.; Mao, S.; Sopha, C.; Sam, B.; Dek, D.; Try, V.; Amato, R.; et al. Dihydroartemisinin—piperaquine resistance in Plasmodium falciparum malaria in Cambodia: A multisite prospective cohort study. Lancet Infect. Dis. 2016, 16, 357–365. [Google Scholar] [CrossRef] [Green Version]
  28. Spring, M.D.; Lin, J.T.; Manning, J.E.; Vanachayangkul, P.; Somethy, S.; Bun, R.; Se, Y.; Chann, S.; Ittiverakul, M.; Sia-ngam, P.; et al. Dihydroartemisinin-piperaquine failure associated with a triple mutant including kelch13 C580Y in Cambodia: An observational cohort study. Lancet Infect. Dis. 2015, 15, 683–691. [Google Scholar] [CrossRef]
  29. Lon, C.; Manning, J.E.; Vanachayangkul, P.; So, M.; Sea, D.; Se, Y.; Gosi, P.; Lanteri, C.; Chaorattanakawee, S.; Sriwichai, S.; et al. Efficacy of two versus three-day regimens of dihydroartemisinin-piperaquine for uncomplicated malaria in military personnel in northern Cambodia: An open-label randomized trial. PLoS ONE 2014, 9, e93138. [Google Scholar] [CrossRef] [PubMed]
  30. Boyle, D. Dispute delays release of malaria money in Cambodia. Lancet 2015, 386, 1811. [Google Scholar] [CrossRef]
  31. World Health Organization (WHO). Global Technical Strategy for Malaria 2016–2030; WHO: Geneva, Switzerland, 2015. [Google Scholar]
  32. Thacker, S.B.; Choi, K.; Brachman, P.S. The Surveillance of Infectious Diseases. JAMA 1983, 249, 1181–1185. [Google Scholar] [CrossRef]
  33. Aregawi, M.; Lynch, M.; Bekele, W.; Kebede, H.; Jima, D.; Taffese, H.S.; Yenehun, M.A.; Lilay, A.; Williams, R.; Thomson, M.; et al. Time series analysis of trends in malaria cases and deaths at hospitals and the effect of antimalarial interventions, 2001–2011, Ethiopia. PLoS ONE 2014, 9, e106359. [Google Scholar] [CrossRef] [PubMed]
  34. Comfort, A.B.; van Dijk, J.H.; Mharakurwa, S.; Stillman, K.; Gabert, R.; Korde, S.; Nachbar, N.; Derriennic, Y.; Musau, S.; Hamazakaza, P.; et al. Hospitalizations and costs incurred at the facility level after scale-up of malaria control: Pre-post comparisons from two hospitals in Zambia. Am. J. Trop. Med. Hyg. 2014, 90, 20–32. [Google Scholar] [CrossRef]
  35. Cissé, B.; Ba, E.H.; Sokhna, C.; Ndiaye, J.L.; Gomis, J.F.; Dial, Y.; Pitt, C.; Ndiaye, M.; Cairns, M.; Faye, E.; et al. Effectiveness of Seasonal Malaria Chemoprevention in Children under Ten Years of Age in Senegal: A Stepped-Wedge Cluster-Randomised Trial. PLoS Med. 2016, 13, e1002175. [Google Scholar] [CrossRef] [Green Version]
  36. Katureebe, A.; Zinszer, K.; Arinaitwe, E.; Rek, J.; Kakande, E.; Charland, K.; Kigozi, R.; Kilama, M.; Nankabirwa, J.; Yeka, A.; et al. Measures of Malaria Burden after Long-Lasting Insecticidal Net Distribution and Indoor Residual Spraying at Three Sites in Uganda: A Prospective Observational Study. PLoS Med. 2016, 13, e1002167. [Google Scholar] [CrossRef] [Green Version]
  37. Aregawi, M.; Malm, K.L.; Wahjib, M.; Kofi, O.; Allotey, N.K.; Yaw, P.N.; Abba-Baffoe, W.; Segbaya, S.; Owusu-Antwi, F.; Kharchi, A.T.; et al. Effect of anti-malarial interventions on trends of malaria cases, hospital admissions and deaths, 2005–2015, Ghana. Malar. J. 2017, 16, 177. [Google Scholar] [CrossRef]
  38. Ssempiira, J.; Kissa, J.; Nambuusi, B.; Kyozira, C.; Rutazaana, D.; Mukooyo, E.; Opigo, J.; Makumbi, F.; Kasasa, S.; Vounatsou, P. The effect of case management and vector-control interventions on space-time patterns of malaria incidence in Uganda. Malar. J. 2018, 17, 162. [Google Scholar] [CrossRef]
  39. Kenangalem, E.; Poespoprodjo, J.R.; Douglas, N.M.; Burdam, F.H.; Gdeumana, K.; Chalfein, F.; Prayoga; Thio, F.; Devine, A.; Marfurt, J.; et al. Malaria morbidity and mortality following introduction of a universal policy of artemisinin-based treatment for malaria in Papua, Indonesia: A longitudinal surveillance study. PLoS Med. 2019, 16, e1002815. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Lechthaler, F.; Matthys, B.; Lechthaler-Felber, G.; Likwela, J.L.; Mavoko, H.M.; Rika, J.M.; Mutombo, M.M.; Ruckstuhl, L.; Barczyk, J.; Shargie, E.; et al. Trends in reported malaria cases and the effects of malaria control in the Democratic Republic of the Congo. PLoS ONE 2019, 14, e0219853. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Tugume, A.; Muneza, F.; Oporia, F.; Kiconco, A.; Kihembo, C.; Kisakye, A.N.; Nsubuga, P.; Deogratias, S.; Yeka, A. Effects and factors associated with indoor residual spraying with Actellic 300 CS on malaria morbidity in Lira District, Northern Uganda. Malar. J. 2019, 18, 44. [Google Scholar] [CrossRef]
  42. Department of Planning and Health Information (DPHI). Health Information System Master Plan 2016–2020; DPHI: Phnom Penh, Cambodia, 2017. [Google Scholar]
  43. National Centre for Parasitology Entomology and Malaria Control (CNM). Malaria Information System; CNM: Phnom Penh, Cambodia, 2020.
  44. World Health Organization (WHO). WHO Malaria Terminology; WHO: Geneva, Switzerland, 2019. [Google Scholar]
  45. OCHA. Cambodia Administrative Level 0–3 Population Statistics; OCHA: New York, NY, USA, 2019. [Google Scholar]
  46. National Institute of Statistics (NIS). Population Projections for Cambodia, 2008–2030; NIS: Phnom Penh, Cambodia, 2008.
  47. National Center for Parasitology Entomology and Malaria Control (CNM). National Treatment Guidelines for Malaria in Cambodia; CNM: Phnom Penh, Cambodia, 2014.
  48. Kulldorff, M. A spatial scan statistic. Commun. Stat. Theory Methods 1997, 26, 1481–1496. [Google Scholar] [CrossRef]
  49. Martin, K. SaTScan™, a Free Software That Analyzes Spatial, Temporal and Space-Time Data Using the Spatial, Temporal, or Space-Time Scan Statistics. Available online: https://www.satscan.org/ (accessed on 20 March 2019).
  50. Maude, R.J.; Nguon, C.; Ly, P.; Bunkea, T.; Ngor, P.; Canavati de la Torre, S.E.; White, N.J.; Dondorp, A.M.; Day, N.P.J.; White, L.J.; et al. Spatial and temporal epidemiology of clinical malaria in Cambodia 2004–2013. Malar. J. 2014, 13, 385. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Gao, H.-W.; Wang, L.-P.; Liang, S.; Liu, Y.-X.; Tong, S.-L.; Wang, J.-J.; Li, Y.-P.; Wang, X.-F.; Yang, H.; Ma, J.-Q.; et al. Change in Rainfall Drives Malaria Re-Emergence in Anhui Province, China. PLoS ONE 2012, 7, e43686. [Google Scholar] [CrossRef]
  52. Wardrop, N.A.; Barnett, A.G.; Atkinson, J.-A.; Clements, A.C.A. Plasmodium vivax malaria incidence over time and its association with temperature and rainfall in four counties of Yunnan Province, China. Malar. J. 2013, 12, 452. [Google Scholar] [CrossRef] [Green Version]
  53. Huang, F.; Zhou, S.; Zhang, S.; Wang, H.; Tang, L. Temporal correlation analysis between malaria and meteorological factors in Motuo County, Tibet. Malar. J. 2011, 10, 54. [Google Scholar] [CrossRef] [Green Version]
  54. Bomblies, A. Modeling the role of rainfall patterns in seasonal malaria transmission. Clim. Chang. 2012, 112, 673–685. [Google Scholar] [CrossRef]
  55. US Department of Interior. United States Geological Survey (USGS); U.S. Department of Interior: Washington, DC, USA, 2020.
  56. Sluydts, V.; Somony, H.; Coosemans, M.; Van Roey, K.; Gryseels, C.; Canier, L.; Kim, S.; Khim, N.; Sovannaroth, S.; Mean, V.; et al. Spatial clustering and risk factors of malaria infections in Ratanakiri Province, Cambodia. Malar. J. 2014, 13, 387. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. An, U.S.; Mao, B.; Saphonn, V.; Bruce, J.; Meek, S.; Lines, J.; Cox, J. Cambodia Malaria Survey 2007; National Centre for Parasitology, Entomology and Malaria Control (CNM): Phnom Penh, Cambodia, 2007.
  58. Dysoley, L.; Rithea, L.; Bunkea, T.; Babu, S.; Sim, K.; Nguon, C.; Socheat, D.; Thompson, M.; Bruce, J.; de Beyl, C.Z.; et al. Cambodia Malaria Survey 2010; National Centre for Parasitology, Entomology and Malaria Control (CNM): Phnom Penh, Cambodia, 2010.
  59. Rekol, H.; Sovannaroth, S.; Dysoley, L.; Vanna, C.; Sokun, C.; Yeang, O.; Thavrin, B.K.; Montha, M.; Yom, Y.; Theasy, H.C.; et al. Cambodia Malaria Survey 2013; National Centre for Parasitology, Entomology and Malaria Control (CNM): Phnom Penh, Cambodia, 2013.
  60. National Centre for Parasitology, Entomology and Malaria Control (CNM). Cambodia Malaria Survey 2017; CNM: Phnom Penh, Cambodia, 2017.
  61. Jambou, R.; Legrand, E.; Niang, M.; Khim, N.; Lim, P.; Volney, B.; Ekala, M.T.; Bouchier, C.; Esterre, P.; Fandeur, T.; et al. Resistance of Plasmodium falciparum field isolates to in-vitro artemether and point mutations of the SERCA-type PfATPase6. Lancet 2005, 366, 1960–1963. [Google Scholar] [CrossRef]
  62. Khim, N.; Bouchier, C.; Ekala, M.T.; Incardona, S.; Lim, P.; Legrand, E.; Jambou, R.; Doung, S.; Puijalon, O.M.; Fandeur, T. Countrywide survey shows very high prevalence of Plasmodium falciparum multilocus resistance genotypes in Cambodia. Antimicrob. Agents Chemother. 2005, 49, 3147–3152. [Google Scholar] [CrossRef] [Green Version]
  63. Denis, M.B.; Tsuyuoka, R.; Lim, P.; Lindegardh, N.; Yi, P.; Top, S.N.; Socheat, D.; Fandeur, T.; Annerberg, A.; Christophel, E.M.; et al. Efficacy of artemether-lumefantrine for the treatment of uncomplicated falciparum malaria in northwest Cambodia. Trop. Med. Int. Health TM IH 2006, 11, 1800–1807. [Google Scholar] [CrossRef] [PubMed]
  64. Resistance to artemisinin derivatives along the Thai-Cambodian border. Relev. Epidemiol. Hebd. 2007, 82, 360.
  65. Shah, N.K.; Alker, A.P.; Sem, R.; Susanti, A.I.; Muth, S.; Maguire, J.D.; Duong, S.; Ariey, F.; Meshnick, S.R.; Wongsrichanalai, C. Molecular surveillance for multidrug-resistant Plasmodium falciparum, Cambodia. Emerg. Infect. Dis. 2008, 14, 1637–1640. [Google Scholar] [CrossRef]
  66. Wongsrichanalai, C.; Meshnick, S.R. Declining artesunate-mefloquine efficacy against falciparum malaria on the Cambodia-Thailand border. Emerg. Infect. Dis. 2008, 14, 716–719. [Google Scholar] [CrossRef] [PubMed]
  67. Antimalarial drug resistance, Thai-Cambodian border. Relev. Epidemiol. Hebd. 2009, 84, 94–95.
  68. Dondorp, A.M.; Nosten, F.; Yi, P.; Das, D.; Phyo, A.P.; Tarning, J.; Lwin, K.M.; Ariey, F.; Hanpithakpong, W.; Lee, S.J.; et al. Artemisinin resistance in Plasmodium falciparum malaria. N. Engl. J. Med. 2009, 361, 455–467. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Müller, O.; Sié, A.; Meissner, P.; Schirmer, R.H.; Kouyaté, B. Artemisinin resistance on the Thai-Cambodian border. Lancet 2009, 374, 1419. [Google Scholar] [CrossRef]
  70. Rogers, W.O.; Sem, R.; Tero, T.; Chim, P.; Lim, P.; Muth, S.; Socheat, D.; Ariey, F.; Wongsrichanalai, C. Failure of artesunate-mefloquine combination therapy for uncomplicated Plasmodium falciparum malaria in southern Cambodia. Malar. J. 2009, 8, 10. [Google Scholar] [CrossRef] [Green Version]
  71. Yeung, S.; Socheat, D.; Moorthy, V.S.; Mills, A.J. Artemisinin resistance on the Thai-Cambodian border. Lancet 2009, 374, 1418–1419. [Google Scholar] [CrossRef]
  72. Duru, V.; Witkowski, B.; Ménard, D. Plasmodium falciparum Resistance to Artemisinin Derivatives and Piperaquine: A Major Challenge for Malaria Elimination in Cambodia. Am. J. Trop. Med. Hyg. 2016, 95, 1228–1238. [Google Scholar] [CrossRef]
  73. Population Services International and ACTwatch. Actwatch Outlet Survey Results: Cambodia 2009–2015; Population Services International and ACTwatch: Phnom Penh, Cambodia, 2016. [Google Scholar]
  74. Siv, S.; Roca-Feltrer, A.; Vinjamuri, S.B.; Bouth, D.M.; Lek, D.; Rashid, M.A.; By, N.P.; Popovici, J.; Huy, R.; Menard, D. Plasmodium vivax Malaria in Cambodia. Am. J. Trop. Med. Hyg. 2016, 95, 97–107. [Google Scholar] [CrossRef] [Green Version]
  75. Office of the Inspector General (OIG). Investigation Report: Global Fund Grants to CambodiaNational Centre for Parasitology Entomology and Malaria Control (CNM); OIG: Geneva, Switzerland, 2017. Available online: https://www.theglobalfund.org/media/5433/oig_gf-oig-17-004_report_en.pdf?u=637278311890000000 (accessed on 5 May 2020).
  76. Titthara, M. Mondulkiri logging crackdowns stun illegal industry. Khmer Times, 31 July 2019. [Google Scholar]
  77. Akulayi, L.; Alum, A.; Andrada, A.; Archer, J.; Arogundade, E.D.; Auko, E.; Badru, A.R.; Bates, K.; Bouanchaud, P.; Bruce, M.; et al. Evidence on anti-malarial and diagnostic markets in Cambodia to guide malaria elimination strategies and policies. Malar. J. 2017, 16, 171. [Google Scholar]
  78. Abba, K.; Kirkham, A.J.; Olliaro, P.L.; Deeks, J.J.; Donegan, S.; Garner, P.; Takwoingi, Y. Rapid diagnostic tests for diagnosing uncomplicated non-falciparum or Plasmodium vivax malaria in endemic countries. Cochrane Database System. Rev. 2014, 2014, Cd011431. [Google Scholar] [CrossRef] [Green Version]
  79. Chu, C.S.; White, N.J. Management of relapsing Plasmodium vivax malaria. Expert Rev. Anti Infect. Ther. 2016, 14, 885–900. [Google Scholar] [CrossRef] [Green Version]
  80. Taylor, A.R.; Watson, J.A.; Chu, C.S.; Puaprasert, K.; Duanguppama, J.; Day, N.P.J.; Nosten, F.; Neafsey, D.E.; Buckee, C.O.; Imwong, M.; et al. Resolving the cause of recurrent Plasmodium vivax malaria probabilistically. Nat. Commun. 2019, 10, 5595. [Google Scholar] [CrossRef]
  81. World Health Organization (WHO). Control and Elimination of Plasmodium Vivax Malaria: Technical Brief; WHO: Geneva, Switzerland, 2015. [Google Scholar]
  82. Lu, G.; Liu, Y.; Beiersmann, C.; Feng, Y.; Cao, J.; Müller, O. Challenges in and lessons learned during the implementation of the 1-3-7 malaria surveillance and response strategy in China: A qualitative study. Infect. Dis. Poverty 2016, 5, 94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Kheang, S.T.; Sovannaroth, S.; Barat, L.M.; Dysoley, L.; Kapella, B.K.; Po, L.; Nguon, S.; Gimnig, J.; Slot, R.; Samphornarann, T.; et al. Malaria elimination using the 1-3-7 approach: Lessons from Sampov Loun, Cambodia. BMC Public Health 2020, 20, 544. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Bruce, M.C.; Donnelly, C.A.; Alpers, M.P.; Galinski, M.R.; Barnwell, J.W.; Walliker, D.; Day, K.P. Cross-Species Interactions Between Malaria Parasites in Humans. Science 2000, 287, 845. [Google Scholar] [CrossRef]
  85. Nagao, Y.; Kimura-Sato, M.; Chavalitshewinkoon-Petmitr, P.; Thongrungkiat, S.; Wilairatana, P.; Ishida, T.; Tan-Ariya, P.; de Souza, J.B.; Krudsood, S.; Looareesuwan, S. Suppression of Plasmodium falciparum by serum collected from a case of Plasmodium vivax infection. Malar. J. 2008, 7, 113. Available online: http://europepmc.org/abstract/MED/18582375; https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pmc/articles/pmid/18582375/pdf/?tool=EBI (accessed on 30 August 2020). [CrossRef] [PubMed] [Green Version]
  86. Whitehorn, J.; Coltart, C.; Manser, D.; Doherty, T. A mixed malaria infection: Is Plasmodium vivax good for you? Trans. R. Soc. Trop. Med. Hyg. 2010, 104, 240–241. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Haghdoost, A.A.; Alexander, N. Systematic review and meta-analysis of the interaction between Plasmodium falciparum and Plasmodium vivax in humans. J. Vector Borne Dis. 2007, 44, 33–43. [Google Scholar]
Figure 1. Number of confirmed (by species) and unconfirmed malaria cases, incidence (confirmed and unconfirmed) per 1000 population, and untoward events, in 2006–2019, Cambodia. Abbreviations: AS-MQ, artesunate–mefloquine; DHA-PIP, dihydroartemisinin–piperaquine; AS-PYR, artesunate–pyronaridine; VMW, village malaria worker.
Figure 1. Number of confirmed (by species) and unconfirmed malaria cases, incidence (confirmed and unconfirmed) per 1000 population, and untoward events, in 2006–2019, Cambodia. Abbreviations: AS-MQ, artesunate–mefloquine; DHA-PIP, dihydroartemisinin–piperaquine; AS-PYR, artesunate–pyronaridine; VMW, village malaria worker.
Ijerph 18 01960 g001
Figure 2. Malaria incidence (confirmed + unconfirmed cases) per 1000 males (a) and females (b), in 2006–2019, Cambodia.
Figure 2. Malaria incidence (confirmed + unconfirmed cases) per 1000 males (a) and females (b), in 2006–2019, Cambodia.
Ijerph 18 01960 g002
Figure 3. Here: Malaria incidence (confirmed + unconfirmed cases) per 1000 by age groups (a) 0–4, (b) 5–14, (c) 15–49, (d) ≥50, in 2006–2019, Cambodia.
Figure 3. Here: Malaria incidence (confirmed + unconfirmed cases) per 1000 by age groups (a) 0–4, (b) 5–14, (c) 15–49, (d) ≥50, in 2006–2019, Cambodia.
Ijerph 18 01960 g003
Figure 4. Malaria incidence (confirmed + unconfirmed) per 100,000 populations per month (2006–2019), in 2006–2019, Cambodia.
Figure 4. Malaria incidence (confirmed + unconfirmed) per 100,000 populations per month (2006–2019), in 2006–2019, Cambodia.
Ijerph 18 01960 g004
Figure 5. Confirmed P. falciparum + Mixed case distribution by district, 2006–2019, Cambodia. Abbreviation: NA, not applicable (no data).
Figure 5. Confirmed P. falciparum + Mixed case distribution by district, 2006–2019, Cambodia. Abbreviation: NA, not applicable (no data).
Ijerph 18 01960 g005
Figure 6. Clusters of confirmed P. falciparum + Mixed cases, 2006–2019, Cambodia. Abbreviation: RR, relative risk.
Figure 6. Clusters of confirmed P. falciparum + Mixed cases, 2006–2019, Cambodia. Abbreviation: RR, relative risk.
Ijerph 18 01960 g006
Figure 7. Confirmed P. vivax + Mixed cases by district, 2006–2019, Cambodia. Abbreviation: NA, not applicable.
Figure 7. Confirmed P. vivax + Mixed cases by district, 2006–2019, Cambodia. Abbreviation: NA, not applicable.
Ijerph 18 01960 g007
Figure 8. Clusters of confirmed P. vivax + Mixed cases, 2006–2019, Cambodia. Abbreviation: RR, relative risk.
Figure 8. Clusters of confirmed P. vivax + Mixed cases, 2006–2019, Cambodia. Abbreviation: RR, relative risk.
Ijerph 18 01960 g008
Table 1. Number of health centers and hospitals by year, reporting to the Health Management Information System (HMIS), Cambodia, 2016 and 2019.
Table 1. Number of health centers and hospitals by year, reporting to the Health Management Information System (HMIS), Cambodia, 2016 and 2019.
YearNumber of Health CentersNumber of HospitalsTotal Malaria Case NotificationsPopulation (a)
2006108796118313,474,489
2007108796118313,676,693
2008108796118313,880,509
2009108796118314,090,208
2010108796118314,308,740
2011108796118314,537,886
2012108796118314,776,866
2013108796118315,022,692
20141138104124215,270,790
20151148107125515,517,635
20161168111127915,762,370
20171195116131116,005,373
20181213121133416,245,454
20191225125135016,489,135
Population data were only available for 2016 [45]. To obtain population data for other years, we projected backward and forward using the National Institute of Statistics’ national population growth rate [46].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chhim, S.; Piola, P.; Housen, T.; Herbreteau, V.; Tol, B. Malaria in Cambodia: A Retrospective Analysis of a Changing Epidemiology 2006–2019. Int. J. Environ. Res. Public Health 2021, 18, 1960. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18041960

AMA Style

Chhim S, Piola P, Housen T, Herbreteau V, Tol B. Malaria in Cambodia: A Retrospective Analysis of a Changing Epidemiology 2006–2019. International Journal of Environmental Research and Public Health. 2021; 18(4):1960. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18041960

Chicago/Turabian Style

Chhim, Srean, Patrice Piola, Tambri Housen, Vincent Herbreteau, and Bunkea Tol. 2021. "Malaria in Cambodia: A Retrospective Analysis of a Changing Epidemiology 2006–2019" International Journal of Environmental Research and Public Health 18, no. 4: 1960. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18041960

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop